-
Notes on Gravitational Physics
Authors:
John L. Friedman
Abstract:
These notes are self-contained, with the first 7 chapters used in a one-semester course with recommended texts by Wald, by Misner, Thorne and Wheeler, and by Schutz. In its treatment of topics covered in these standard texts, the presentation here typically includes steps skipped in Wald or MTW. Treatments of gravitational waves, particle orbits in black-hole backgrounds, the Teukolsky equation, a…
▽ More
These notes are self-contained, with the first 7 chapters used in a one-semester course with recommended texts by Wald, by Misner, Thorne and Wheeler, and by Schutz. In its treatment of topics covered in these standard texts, the presentation here typically includes steps skipped in Wald or MTW. Treatments of gravitational waves, particle orbits in black-hole backgrounds, the Teukolsky equation, and the initial value equations are motivated in part by the discoveries of gravitational waves from the inspiral and coalescence of binary black holes and neutron stars, advances in numerical relativity, and the expected LISA space-based observatory. The notes begin with a detailed presentation of special relativity with a geometrical orientation, starting with with time dilation and length contraction and including relativistic particles, fluids, electromagnetism, and curvilinear coordinates. Chaps. 2-5 cover curvature, the Einstein equation, relativistic stars, and black holes. Chap. 6, on gravitational waves, includes a discussion of detection and noise. Chap. 7 is a brief introduction to cosmology, deriving the metrics of homogeneous isotropic space, the equations governing a universe with matter, radiation and vacuum energy, and their solutions, and discussions of the cosmological redshift and on using gravitational waves to measure the Hubble constant.
Chap. 8, on the initial value problem, has a section on the form of the equations used in numerical relativity. The Newman-Penrose formalism and the Teukolsky equation are covered in Chap. 9. Following that is a chapter on black-hole thermodynamics and a final chapter on the gravitational action and on conserved quantities for asymptotically flat spacetimes, using Noether's theorem. An appendix covers forms, densities, integration, and Cartan calculus.
△ Less
Submitted 21 May, 2024; v1 submitted 18 August, 2023;
originally announced August 2023.
-
Science opportunities with solar sailing smallsats
Authors:
Slava G. Turyshev,
Darren Garber,
Louis D. Friedman,
Andreas M. Hein,
Nathan Barnes,
Konstantin Batygin,
G. David Brin,
Michael E. Brown,
Leroy Cronin,
Artur Davoyan,
Amber Dubill,
T. Marshall Eubanks,
Sarah Gibson,
Donald M. Hassler,
Noam R. Izenberg,
Pierre Kervella,
Philip D. Mauskopf,
Neil Murphy,
Andrew Nutter,
Carolyn Porco,
Dario Riccobono,
James Schalkwyk,
Kevin B. Stevenson,
Mark V. Sykes,
Mahmooda Sultana
, et al. (3 additional authors not shown)
Abstract:
Recently, we witnessed how the synergy of small satellite technology and solar sailing propulsion enables new missions. Together, small satellites with lightweight instruments and solar sails offer affordable access to deep regions of the solar system, also making it possible to realize hard-to-reach trajectories that are not constrained to the ecliptic plane. Combining these two technologies can…
▽ More
Recently, we witnessed how the synergy of small satellite technology and solar sailing propulsion enables new missions. Together, small satellites with lightweight instruments and solar sails offer affordable access to deep regions of the solar system, also making it possible to realize hard-to-reach trajectories that are not constrained to the ecliptic plane. Combining these two technologies can drastically reduce travel times within the solar system, while delivering robust science. With solar sailing propulsion capable of reaching the velocities of ~5-10 AU/yr, missions using a rideshare launch may reach the Jovian system in two years, Saturn in three. The same technologies could allow reaching solar polar orbits in less than two years. Fast, cost-effective, and maneuverable sailcraft that may travel outside the ecliptic plane open new opportunities for affordable solar system exploration, with great promise for heliophysics, planetary science, and astrophysics. Such missions could be modularized to reach different destinations with different sets of instruments. Benefiting from this progress, we present the "Sundiver" concept, offering novel possibilities for the science community. We discuss some of the key technologies, the current design of the Sundiver sailcraft vehicle and innovative instruments, along with unique science opportunities that these technologies enable, especially as this exploration paradigm evolves. We formulate policy recommendations to allow national space agencies, industry, and other stakeholders to establish a strong scientific, programmatic, and commercial focus, enrich and deepen the space enterprise and broaden its advocacy base by including the Sundiver paradigm as a part of broader space exploration efforts.
△ Less
Submitted 19 July, 2023; v1 submitted 27 March, 2023;
originally announced March 2023.
-
Solar Sail Propulsion by 2050: An Enabling Capability for Heliophysics Missions
Authors:
Les Johnson,
Nathan Barnes,
Matteo Ceriotti,
Thomas Y. Chen,
Artur Davoyan,
Louis Friedman,
Darren Garber,
Roman Kezerashvili,
Ken Kobayashi,
Greg Matloff,
Colin McInnes,
Pat Mulligan,
Grover Swartzlander,
Slava G. Turyshev
Abstract:
Solar sails enable missions to observe the solar environment from unique vantage points, such as sustained observations away from the Sun-Earth line; sub-L1 station keeping; high inclination solar orbits; Earth polar-sitting and polar-viewing observatories; fast transit missions to study heliosphere to interstellar medium transition, as well as missions of interest across a broad user community. R…
▽ More
Solar sails enable missions to observe the solar environment from unique vantage points, such as sustained observations away from the Sun-Earth line; sub-L1 station keeping; high inclination solar orbits; Earth polar-sitting and polar-viewing observatories; fast transit missions to study heliosphere to interstellar medium transition, as well as missions of interest across a broad user community. Recent and planned demonstration missions make this technology ready for use on near-term science missions.
△ Less
Submitted 2 January, 2023;
originally announced January 2023.
-
A Mission to Nature's Telescope for High-Resolution Imaging of an Exoplanet
Authors:
Louis D. Friedman,
Darren Garber,
Slava G. Turyshev,
Henry Helvajian,
Thomas Heinshiemer,
John McVey,
Artur R. Davoyan
Abstract:
The solar gravitational lens (SGL) provides a factor of $10^{11}$ amplification for viewing distant point sources beyond our solar system. As such, it may be used for resolved imaging of extended sources, such as exoplanets, not possible otherwise. To use the SGL, a spacecraft carrying a modest telescope and a coronagraph must reach the SGLs focal region, that begins at $\sim$550 astronomical unit…
▽ More
The solar gravitational lens (SGL) provides a factor of $10^{11}$ amplification for viewing distant point sources beyond our solar system. As such, it may be used for resolved imaging of extended sources, such as exoplanets, not possible otherwise. To use the SGL, a spacecraft carrying a modest telescope and a coronagraph must reach the SGLs focal region, that begins at $\sim$550 astronomical units (AU) from the Sun and is oriented outward along the line connecting the distant object and the Sun. No spacecraft has ever reached even a half of that distance; and to do so within a reasonable mission lifetime (e.g., less than 25 years) and affordable cost requires a new type of mission design, using solar sails and microsats ($<100$~kg). The payoff is high -- using the SGL is the only practical way we can ever get a high-resolution, multi-pixel image of an Earth-like exoplanet, one that we identify as potentially habitable. This paper describes a novel mission design starting with a rideshare launch from the Earth, spiraling in toward the Sun, and then flying around it to achieve solar system exit speeds of over $20$ AU/year. A new sailcraft design is used to make possible high area to mass ratio for the sailcraft. The mission design enables other fast solar system missions, starting with a proposed very low cost technology demonstration mission (TDM) to prove the functionality and operation of the microsat-solar sail design and then, building on the TDM, missions to explore distant regions of the solar system, and those to study Kuiper Belt objects (KBOs) and the recently discovered interstellar objects (ISOs) are also possible.
△ Less
Submitted 23 July, 2021;
originally announced July 2021.
-
A Fast Response Mission to Rendezvous with an Interstellar Object
Authors:
Darren Garber,
Louis D. Friedman,
Artur R. Davoyan,
Slava G. Turyshev,
Nahum Melamed,
John McVey,
Todd F. Sheerin
Abstract:
A solar sail propelled small satellite mission concept to intercept and potentially rendezvous with newly discovered transient interstellar objects (ISOs) is described. The mission concept derives from the proposal for a technology demonstration mission for exiting the solar system at high velocity, eventually to reach the focal region of the solar gravitational lens. The ISO mission concept is to…
▽ More
A solar sail propelled small satellite mission concept to intercept and potentially rendezvous with newly discovered transient interstellar objects (ISOs) is described. The mission concept derives from the proposal for a technology demonstration mission for exiting the solar system at high velocity, eventually to reach the focal region of the solar gravitational lens. The ISO mission concept is to fly a solar sail toward a holding orbit around the Sun and when the ISO orbit is confirmed, target the sailcraft to reach an escape velocity of over 6 AU/year. This would permit rapid response to a new ISO discovery and an intercept within 10 AU from the Sun. Two new proven interplanetary technologies are utilized to enable such a mission: i) interplanetary smallsats, such as those demonstrated by the MarCO mission, and ii) solar sails, such as demonstrated by LightSail and IKAROS missions and developed for NEA Scout and Solar Cruiser missions. Current technology work suggests that already within this decade such a mission could fly and reach an ISO moving through the solar system. It might enable the first encounter with an ISO to allow for imaging and spectroscopy, measurements of size and mass, potentially giving a unique information about the object's origin and composition. A similar approach may be used to allow for a sample return.
△ Less
Submitted 8 February, 2022; v1 submitted 27 June, 2021;
originally announced June 2021.
-
Exploring the Outer Solar System with Solar Sailing Smallsats on Fast-Transit Trajectories and In-Flight Autonomous Assembly of Advanced Science Payloads
Authors:
Slava G. Turyshev,
Henry Helvajian,
Louis D. Friedman,
Tom Heinsheimer,
Darren Garber,
Artur Davoyan,
Viktor T. Toth
Abstract:
We discuss the in-flight autonomous assembly as the means to build advanced planetary science payloads to explore the outer regions of the solar system. These payloads are robotically constructed from modular parts delivered by a group of smallsats (< 20 kg) which are placed on fast solar system transfer trajectories while being accelerated by solar sail propulsion to velocities of ~10 AU/yr. This…
▽ More
We discuss the in-flight autonomous assembly as the means to build advanced planetary science payloads to explore the outer regions of the solar system. These payloads are robotically constructed from modular parts delivered by a group of smallsats (< 20 kg) which are placed on fast solar system transfer trajectories while being accelerated by solar sail propulsion to velocities of ~10 AU/yr. This concept provides the planetary science community with inexpensive, frequent access to distant regions of the solar system with flexible, reconfigurable instruments and systems that are assembled in flight. It permits faster revisit times, rapid replenishment and technology insertions, longer mission capability with lower costs. It also increases the science capabilities of smallsats via the use of modular, redundant architectures and allows for proliferation of sensing instrumentation throughout the solar system.
△ Less
Submitted 10 July, 2020;
originally announced July 2020.
-
Astrophysical Implications of Neutron Star Inspiral and Coalescence
Authors:
John L. Friedman,
Nikolaos Stergioulas
Abstract:
The first inspiral of two neutron stars observed in gravitational waves was remarkably close, allowing the kind of simultaneous gravitational wave and electromagnetic observation that had not been expected for several years. Their merger, followed by a gamma-ray burst and a kilonova, was observed across the spectral bands of electromagnetic telescopes. These GW and electromagnetic observations hav…
▽ More
The first inspiral of two neutron stars observed in gravitational waves was remarkably close, allowing the kind of simultaneous gravitational wave and electromagnetic observation that had not been expected for several years. Their merger, followed by a gamma-ray burst and a kilonova, was observed across the spectral bands of electromagnetic telescopes. These GW and electromagnetic observations have led to dramatic advances in understanding short gamma-ray bursts; determining the origin of the heaviest elements; and determining the maximum mass of neutron stars. From the imprint of tides on the gravitational waveforms and from observations of X-ray binaries, one can extract the radius and deformability of inspiraling neutron stars. Together, the radius, maximum mass, and causality constrain the neutron-star equation of state, and future constraints can come from observations of post-merger oscillations. We selectively review these results, filling in some of the physics with derivations and estimates.
△ Less
Submitted 9 January, 2022; v1 submitted 28 May, 2020;
originally announced May 2020.
-
Direct Multipixel Imaging and Spectroscopy of an Exoplanet with a Solar Gravity Lens Mission
Authors:
Slava G. Turyshev,
Michael Shao,
Viktor T. Toth,
Louis D. Friedman,
Leon Alkalai,
Dmitri Mawet,
Janice Shen,
Mark R. Swain,
Hanying Zhou,
Henry Helvajian,
Tom Heinsheimer,
Siegfried Janson,
Zigmond Leszczynski,
John McVey,
Darren Garber,
Artur Davoyan,
Seth Redfield,
Jared R. Males
Abstract:
We examined the solar gravitational lens (SGL) as the means to produce direct high-resolution, multipixel images of exoplanets. The properties of the SGL are remarkable: it offers maximum light amplification of ~1e11 and angular resolution of ~1e-10 arcsec. A probe with a 1-m telescope in the SGL focal region can image an exoplanet at 30 pc with 10-kilometer resolution on its surface, sufficient t…
▽ More
We examined the solar gravitational lens (SGL) as the means to produce direct high-resolution, multipixel images of exoplanets. The properties of the SGL are remarkable: it offers maximum light amplification of ~1e11 and angular resolution of ~1e-10 arcsec. A probe with a 1-m telescope in the SGL focal region can image an exoplanet at 30 pc with 10-kilometer resolution on its surface, sufficient to observe seasonal changes, oceans, continents, surface topography. We reached and exceeded all objectives set for our study: We developed a new wave-optical approach to study the imaging of exoplanets while treating them as extended, resolved, faint sources at large but finite distances. We properly accounted for the solar corona brightness. We developed deconvolution algorithms and demonstrated the feasibility of high-quality image reconstruction under realistic conditions. We have proven that multipixel imaging and spectroscopy of exoplanets with the SGL are feasible. We have developed a new mission concept that delivers an array of optical telescopes to the SGL focal region relying on three innovations: i) a new way to enable direct exoplanet imaging, ii) use of smallsats solar sails fast transit through the solar system and beyond, iii) an open architecture to take advantage of swarm technology. This approach enables entirely new missions, providing a great leap in capabilities for NASA and the greater aerospace community. Our results are encouraging as they lead to a realistic design for a mission that will be able to make direct resolved images of exoplanets in our stellar neighborhood. It could allow exploration of exoplanets relying on the SGL capabilities decades, if not centuries, earlier than possible with other extant technologies. The architecture and mission concepts for a mission to the strong interference region of the SGL are promising and should be explored further.
△ Less
Submitted 19 March, 2020; v1 submitted 26 February, 2020;
originally announced February 2020.
-
Long-baseline horizontal radio-frequency transmission through polar ice
Authors:
P. Allison,
S. Archambault,
J. J. Beatty,
D. Z. Besson,
C. C. Chen,
C. H. Chen,
P. Chen,
A. Christenson,
B. A. Clark,
W. Clay,
A. Connolly,
L. Cremonesi,
C. Deaconu,
M. Duvernois,
L. Friedman,
R. Gaior,
J. Hanson,
K. Hanson,
J. Haugen,
K. D. Hoffman,
E. Hong,
S. Y. Hsu,
L. Hu,
J. J. Huang,
A. M. -H. Huang
, et al. (41 additional authors not shown)
Abstract:
We report on analysis of englacial radio-frequency (RF) pulser data received over horizontal baselines of 1--5 km, based on broadcasts from two sets of transmitters deployed to depths of up to 1500 meters at the South Pole. First, we analyze data collected usingtwo RF bicone transmitters 1400 meters below the ice surface, and frozen into boreholes drilled for the IceCube experiment in 2011. Additi…
▽ More
We report on analysis of englacial radio-frequency (RF) pulser data received over horizontal baselines of 1--5 km, based on broadcasts from two sets of transmitters deployed to depths of up to 1500 meters at the South Pole. First, we analyze data collected usingtwo RF bicone transmitters 1400 meters below the ice surface, and frozen into boreholes drilled for the IceCube experiment in 2011. Additionally, in Dec., 2018, a fat-dipole antenna, fed by one of three high-voltage (~1 kV), fast (~(1-5 ns)) signal generators was lowered into the 1700-m deep icehole drilled for the South Pole Ice Core Experiment (SPICE), approximately 3 km from the geographic South Pole. Signals from transmitters were recorded on the five englacial multi-receiver ARA stations, with receiver depths between 60--200 m. We confirm the long, >1 km RF electric field attenuation length, test our observed signal arrival timing distributions against models, and measure birefringent asymmetries at the 0.15% level.
△ Less
Submitted 19 January, 2021; v1 submitted 28 August, 2019;
originally announced August 2019.
-
Recognizing the Value of the Solar Gravitational Lens for Direct Multipixel Imaging and Spectroscopy of an Exoplanet
Authors:
Slava G. Turyshev,
Michael Shao,
Janice Shen,
Hanying Zhou,
Viktor T. Toth,
Louis Friedman,
Leon Alkalai,
Nitin Arora,
Darren D. Garber,
Henry Helvajian,
Thomas Heinsheimer,
Siegfried W. Janson,
Les Johnson,
Jared R. Males,
Roy Nakagawa,
Seth Redfield,
Nathan Strange,
Mark R. Swain,
David Van Buren,
John L. West,
Stacy Weinstein-Weiss
Abstract:
The Solar Gravitational Lens (SGL) allows for major brightness amplification ($\sim 10^{11}$ at wavelength of $1~μ$m) and extreme angular resolution ($\sim10^{-10}$ arcsec) within a narrow field of view. A meter-class telescope, with a modest coronagraph to block solar light with 1e-6 suppression placed in the focal area of the SGL, can image an exoplanet at a distance of 30 parsec with few kilome…
▽ More
The Solar Gravitational Lens (SGL) allows for major brightness amplification ($\sim 10^{11}$ at wavelength of $1~μ$m) and extreme angular resolution ($\sim10^{-10}$ arcsec) within a narrow field of view. A meter-class telescope, with a modest coronagraph to block solar light with 1e-6 suppression placed in the focal area of the SGL, can image an exoplanet at a distance of 30 parsec with few kilometer-scale resolution on its surface. Notably, spectroscopic broadband SNR is $\sim 10^{-6}$ in two weeks of integration time, providing this instrument with incredible remote sensing capabilities. A mission capable of exploiting the remarkable optical properties of the SGL allows for direct high-resolution imaging/spectroscopy of a potentially habitable exoplanet. Such missions could allow exploration of exoplanets relying on the SGL capabilities decades, if not centuries, earlier than possible with other extant technologies.
△ Less
Submitted 12 March, 2018;
originally announced March 2018.
-
Magnetic field amplification by the r-mode instability
Authors:
Andrey I. Chugunov,
John L. Friedman,
Lee Lindblom,
Luciano Rezzolla
Abstract:
We discuss magnetic field enhancement by unstable r-modes (driven by the gravitational radiation reaction force) in rotating stars. In the absence of a magnetic field, gravitational radiation exponentially increases the r-mode amplitude $α$, and accelerates differential rotation (secular motion of fluid elements). For a magnetized star, differential rotation enhances the magnetic field energy. Rez…
▽ More
We discuss magnetic field enhancement by unstable r-modes (driven by the gravitational radiation reaction force) in rotating stars. In the absence of a magnetic field, gravitational radiation exponentially increases the r-mode amplitude $α$, and accelerates differential rotation (secular motion of fluid elements). For a magnetized star, differential rotation enhances the magnetic field energy. Rezzolla et al. (2000--2001) argued that if the magnetic energy grows faster than the gravitational radiation reaction force pumps energy into the r-modes, then the r-mode instability is suppressed. Chugunov (2015) demonstrated that without gravitational radiation, differential rotation can be treated as a degree of freedom decoupled from the r-modes and controlled by the back reaction of the magnetic field. In particular, the magnetic field windup does not damp r-modes. Here we discuss the effect of the back reaction of the magnetic field on differential rotation of unstable r-modes, and show that it limits the generated magnetic field and the magnetic energy growth rate preventing suppression of the r-mode instability by magnetic windup at low saturation amplitudes, $α\ll 1$, predicted by current models.
△ Less
Submitted 26 December, 2017;
originally announced December 2017.
-
Measurement of the real dielectric permittivity epsilon_r of glacial ice
Authors:
P. Allison,
S. Archambault,
J. Auffenberg,
R. Bard,
J. J. Beatty,
M. Beheler-Amass,
D. Z. Besson,
M. Beydler,
C. Brabec,
C. -C. Chen,
C. -H. Chen,
P. Chen,
A. Christenson,
B. A. Clark,
A. Connolly,
L. Cremonesi,
C. Deaconu,
M. Duvernois,
L. Friedman,
R. Gaior,
P. W. Gorham,
J. Hanson,
K. Hanson,
J. Haugen,
K. D. Hoffman
, et al. (44 additional authors not shown)
Abstract:
Using data collected by the Askaryan Radio Array (ARA) experiment at the South Pole, we have used long-baseline propagation of radio-frequency signals to extract information on the radio-frequency index-of-refraction in South Polar ice. Owing to the increasing ice density over the upper 150--200 meters, rays are observed along two, nearly parallel paths, one of which is direct and a second which r…
▽ More
Using data collected by the Askaryan Radio Array (ARA) experiment at the South Pole, we have used long-baseline propagation of radio-frequency signals to extract information on the radio-frequency index-of-refraction in South Polar ice. Owing to the increasing ice density over the upper 150--200 meters, rays are observed along two, nearly parallel paths, one of which is direct and a second which refracts through an inflection point, with differences in both arrival time and arrival angle that can be used to constrain the neutrino properties. We also observe indications, for the first time, of radio-frequency ice birefringence for signals propagating along predominantly horizontal trajectories, corresponding to an asymmetry of order 0.1% between the ordinary and extra-ordinary paths, numerically compatible with previous measurements of birefringent asymmetries for vertically-propagating radio-frequency signals at South Pole. Taken together, these effects offer the possibility of redundantly measuring the range from receiver to a neutrino interaction in Antarctic ice, if receiver antennas are deployed at shallow (25 m<z<100 m) depths. Such range information is essential in determining both the neutrino energy, as well as the incident neutrino direction.
△ Less
Submitted 14 January, 2019; v1 submitted 8 December, 2017;
originally announced December 2017.
-
Limits on Magnetic Field Amplification from the r-Mode Instability
Authors:
John L. Friedman,
Lee Lindblom,
Luciano Rezzolla,
Andrey I. Chugunov
Abstract:
At second order in perturbation theory, the unstable r-mode of a rotating star includes growing differential rotation whose form and growth rate are determined by gravitational-radiation reaction. With no magnetic field, the angular velocity of a fluid element grows exponentially until the mode reaches its nonlinear saturation amplitude and remains nonzero after saturation. With a background magne…
▽ More
At second order in perturbation theory, the unstable r-mode of a rotating star includes growing differential rotation whose form and growth rate are determined by gravitational-radiation reaction. With no magnetic field, the angular velocity of a fluid element grows exponentially until the mode reaches its nonlinear saturation amplitude and remains nonzero after saturation. With a background magnetic field, the differential rotation winds up and amplifies the field, and previous work where large mode amplitudes were considered suggests that the amplification may damp out the instability. A background magnetic field, however, turns the saturated time-independent perturbations corresponding to adding differential rotation into perturbations whose characteristic frequencies are of order the Alfvén frequency. As found in previous studies, we argue that magnetic- field growth is sharply limited by the saturation amplitude of an unstable mode. In contrast to previous work, however, we show that if the amplitude is small, i.e., of order 10^(-4), then the limit on the magnetic-field growth is stringent enough to prevent the loss of energy to the magnetic field from damping or significantly altering an unstable r-mode in nascent neutron stars with normal interiors and in cold stars whose interiors are type II superconductors. We show this result first for a toy model, and we then obtain an analogous upper limit on magnetic field growth using a more realistic model of a rotating neutron star. Our analysis depends on the assumption that there are no marginally unstable perturbations, and this may not hold when differential rotation leads to a magnetorotational instability.
△ Less
Submitted 2 March, 2018; v1 submitted 28 July, 2017;
originally announced July 2017.
-
Upper limit set by causality on the tidal deformability of a neutron star
Authors:
Eric D. Van Oeveren,
John L. Friedman
Abstract:
A principal goal of gravitational-wave astronomy is to constrain the neutron star equation of state (EOS) by measuring the tidal deformability of neutron stars. The tidally induced departure of the waveform from that of point-particle (or spinless binary black hole (BBH)) increases with the stiffness of the EOS. We show that causality (the requirement that the speed of sound is less than the speed…
▽ More
A principal goal of gravitational-wave astronomy is to constrain the neutron star equation of state (EOS) by measuring the tidal deformability of neutron stars. The tidally induced departure of the waveform from that of point-particle (or spinless binary black hole (BBH)) increases with the stiffness of the EOS. We show that causality (the requirement that the speed of sound is less than the speed of light for a perfect fluid satisfying a one-parameter equation of state) places an upper bound on tidal deformability as a function of mass. Like the upper mass limit, the limit on deformabity is obtained by using an EOS with $v_{sound} = c$ for high densities and matching to a low density (candidate) EOS at a matching density of order nuclear saturation density. We use these results and those of [B.D. Lackey et al., Phys. Rev. D 89, 043009 (2014)] to estimate the resulting upper limit on the gravitational-wave phase shift of a black hole-neutron star (BHNS) binary relative to a BBH. Even for assumptions weak enough to allow an maximum mass of $4\ M_\odot$ (a match at nuclear saturation density to an unusually stiff low-density candidate EOS), the upper limit on dimensionless tidal deformability is stringent. It leads to a still more stringent estimated upper limit on the maximum tidally induced phase shift prior to merger. We comment in an appendix on the relation between causality, the condition $v_{sound} < c$, and the condition $dp/dε< 1$ for the effective EOS governing the equilibrium star.
△ Less
Submitted 3 May, 2017; v1 submitted 13 January, 2017;
originally announced January 2017.
-
Differential rotation of the unstable nonlinear r-modes
Authors:
John L. Friedman,
Lee Lindblom,
Keith H. Lockitch
Abstract:
At second order in perturbation theory, the $r$-modes of uniformly rotating stars include an axisymmetric part that can be identified with differential rotation of the background star. If one does not include radiation-reaction, the differential rotation is constant in time and has been computed by Sá. It has a gauge dependence associated with the family of time-independent perturbations that add…
▽ More
At second order in perturbation theory, the $r$-modes of uniformly rotating stars include an axisymmetric part that can be identified with differential rotation of the background star. If one does not include radiation-reaction, the differential rotation is constant in time and has been computed by Sá. It has a gauge dependence associated with the family of time-independent perturbations that add differential rotation to the unperturbed equilibrium star: For stars with a barotropic equation of state, one can add to the time-independent second-order solution arbitrary differential rotation that is stratified on cylinders (that is a function of distance $\varpi$ to the axis of rotation). We show here that the gravitational radiation-reaction force that drives the $r$-mode instability removes this gauge freedom: The expontially growing differential rotation of the unstable second-order $r$-mode is unique. We derive a general expression for this rotation law for Newtonian models and evaluate it explicitly for slowly rotating models with polytropic equations of state.
△ Less
Submitted 5 March, 2018; v1 submitted 30 March, 2015;
originally announced March 2015.
-
Matter effects on binary neutron star waveforms
Authors:
Jocelyn S. Read,
Luca Baiotti,
Jolien D. E. Creighton,
John L. Friedman,
Bruno Giacomazzo,
Koutarou Kyutoku,
Charalampos Markakis,
Luciano Rezzolla,
Masaru Shibata,
Keisuke Taniguchi
Abstract:
Using an extended set of equations of state and a multiple-group multiple-code collaborative effort to generate waveforms, we improve numerical-relativity-based data-analysis estimates of the measurability of matter effects in neutron-star binaries. We vary two parameters of a parameterized piecewise-polytropic equation of state (EOS) to analyze the measurability of EOS properties, via a parameter…
▽ More
Using an extended set of equations of state and a multiple-group multiple-code collaborative effort to generate waveforms, we improve numerical-relativity-based data-analysis estimates of the measurability of matter effects in neutron-star binaries. We vary two parameters of a parameterized piecewise-polytropic equation of state (EOS) to analyze the measurability of EOS properties, via a parameter Λ that characterizes the quadrupole deformability of an isolated neutron star. We find that, to within the accuracy of the simulations, the departure of the waveform from point-particle (or spinless double black-hole binary) inspiral increases monotonically with Λ, and changes in the EOS that did not change Λ are not measurable. We estimate with two methods the minimal and expected measurability of Λ in second- and third- generation gravitational-wave detectors. The first estimate, using numerical waveforms alone, shows two EOS which vary in radius by 1.3km are distinguishable in mergers at 100Mpc. The second estimate relies on the construction of hybrid waveforms by matching to post-Newtonian inspiral, and estimates that the same EOS are distinguishable in mergers at 300Mpc. We calculate systematic errors arising from numerical uncertainties and hybrid construction, and we estimate the frequency at which such effects would interfere with template-based searches.
△ Less
Submitted 17 June, 2013;
originally announced June 2013.
-
Neutron star equation of state via gravitational wave observations
Authors:
Charalampos Markakis,
Jocelyn S. Read,
Masaru Shibata,
Koji Uryu,
Jolien D. E. Creighton,
John L. Friedman,
Benjamin D. Lackey
Abstract:
Gravitational wave observations can potentially measure properties of neutron star equations of state by measuring departures from the point-particle limit of the gravitational waveform produced in the late inspiral of a neutron star binary. Numerical simulations of inspiraling neutron star binaries computed for equations of state with varying stiffness are compared. As the stars approach their fi…
▽ More
Gravitational wave observations can potentially measure properties of neutron star equations of state by measuring departures from the point-particle limit of the gravitational waveform produced in the late inspiral of a neutron star binary. Numerical simulations of inspiraling neutron star binaries computed for equations of state with varying stiffness are compared. As the stars approach their final plunge and merger, the gravitational wave phase accumulates more rapidly if the neutron stars are more compact. This suggests that gravitational wave observations at frequencies around 1 kHz will be able to measure a compactness parameter and place stringent bounds on possible neutron star equations of state. Advanced laser interferometric gravitational wave observatories will be able to tune their frequency band to optimize sensitivity in the required frequency range to make sensitive measures of the late-inspiral phase of the coalescence.
△ Less
Submitted 18 October, 2011; v1 submitted 17 October, 2011;
originally announced October 2011.
-
Extracting equation of state parameters from black hole-neutron star mergers. I. Nonspinning black holes
Authors:
Benjamin D. Lackey,
Koutarou Kyutoku,
Masaru Shibata,
Patrick R. Brady,
John L. Friedman
Abstract:
The late inspiral, merger, and ringdown of a black hole-neutron star (BHNS) system can provide information about the neutron-star equation of state (EOS). Candidate EOSs can be approximated by a parametrized piecewise-polytropic EOS above nuclear density, matched to a fixed low-density EOS; and we report results from a large set of BHNS inspiral simulations that systematically vary two parameters.…
▽ More
The late inspiral, merger, and ringdown of a black hole-neutron star (BHNS) system can provide information about the neutron-star equation of state (EOS). Candidate EOSs can be approximated by a parametrized piecewise-polytropic EOS above nuclear density, matched to a fixed low-density EOS; and we report results from a large set of BHNS inspiral simulations that systematically vary two parameters. To within the accuracy of the simulations, we find that, apart from the neutron-star mass, a single physical parameter Lambda, describing its deformability, can be extracted from the late inspiral, merger, and ringdown waveform. This parameter is related to the radius, mass, and l=2 Love number, k_2, of the neutron star by Lambda = 2k_2 R^5/3M_{NS}^5, and it is the same parameter that determines the departure from point-particle dynamics during the early inspiral. Observations of gravitational waves from BHNS inspiral thus restrict the EOS to a surface of constant Lambda in the parameter space, thickened by the measurement error. Using various configurations of a single Advanced LIGO detector, we find that neutron stars are distinguishable from black holes of the same mass and that Lambda^{1/5} or equivalently R can be extracted to 10-40% accuracy from single events for mass ratios of Q=2 and 3 at a distance of 100 Mpc, while with the proposed Einstein Telescope, EOS parameters can be extracted to accuracy an order of magnitude better.
△ Less
Submitted 15 September, 2011;
originally announced September 2011.
-
Inferring the neutron star equation of state from binary inspiral waveforms
Authors:
Charalampos Markakis,
Jocelyn S. Read,
Masaru Shibata,
Koji Uryu,
Jolien D. E. Creighton,
John L. Friedman
Abstract:
The properties of neutron star matter above nuclear density are not precisely known. Gravitational waves emitted from binary neutron stars during their late stages of inspiral and merger contain imprints of the neutron-star equation of state. Measuring departures from the point-particle limit of the late inspiral waveform allows one to measure properties of the equation of state via gravitational…
▽ More
The properties of neutron star matter above nuclear density are not precisely known. Gravitational waves emitted from binary neutron stars during their late stages of inspiral and merger contain imprints of the neutron-star equation of state. Measuring departures from the point-particle limit of the late inspiral waveform allows one to measure properties of the equation of state via gravitational wave observations. This and a companion talk by J. S. Read reports a comparison of numerical waveforms from simulations of inspiraling neutron-star binaries, computed for equations of state with varying stiffness. We calculate the signal strength of the difference between waveforms for various commissioned and proposed interferometric gravitational wave detectors and show that observations at frequencies around 1 kHz will be able to measure a compactness parameter and constrain the possible neutron-star equations of state.
△ Less
Submitted 10 August, 2010;
originally announced August 2010.
-
Non-conformally flat initial data for binary compact objects
Authors:
Koji Uryu,
Francois Limousin,
John L. Friedman,
Eric Gourgoulhon,
Masaru Shibata
Abstract:
A new method is described for constructing initial data for a binary neutron-star (BNS) system in quasi-equilibrium circular orbit. Two formulations for non-conformally flat data, waveless (WL) and near-zone helically symmetric (NHS), are introduced; in each formulation, the Einstein-Euler system, written in 3+1 form on an asymptotically flat spacelike hypersurface, is exactly solved for all met…
▽ More
A new method is described for constructing initial data for a binary neutron-star (BNS) system in quasi-equilibrium circular orbit. Two formulations for non-conformally flat data, waveless (WL) and near-zone helically symmetric (NHS), are introduced; in each formulation, the Einstein-Euler system, written in 3+1 form on an asymptotically flat spacelike hypersurface, is exactly solved for all metric components, including the spatially non-conformally flat potentials, and for irrotational flow. A numerical method applicable to both formulations is explained with an emphasis on the imposition of a spatial gauge condition. Results are shown for solution sequences of irrotational BNS with matter approximated by parametrized equations of state that use a few segments of polytropic equations of state. The binding energy and total angular momentum of solution sequences computed within the conformally flat -- Isenberg-Wilson-Mathews (IWM) -- formulation are closer to those of the third post-Newtonian (3PN) two point particles up to the closest orbits, for the more compact stars, whereas sequences resulting from the WL/NHS formulations deviate from the 3PN curve even more for the sequences with larger compactness. We think it likely that this correction reflects an overestimation in the IWM formulation as well as in the 3PN formula, by $\sim 1$ cycle in the gravitational wave phase during the last several orbits. The work suggests that imposing spatial conformal flatness results in an underestimate of the quadrupole deformation of the components of binary neutron-star systems in the last few orbits prior to merger.
△ Less
Submitted 8 October, 2009; v1 submitted 5 August, 2009;
originally announced August 2009.
-
Iteration Stability for Simple Newtonian Stellar Systems
Authors:
Richard H. Price,
Charalampos Markakis,
John L. Friedman
Abstract:
For an equation of state in which pressure is a function only of density, the analysis of Newtonian stellar structure is simple in principle if the system is axisymmetric, or consists of a corotating binary. It is then required only to solve two equations: one stating that the "injection energy", $κ$, a potential, is constant throughout the stellar fluid, and the other being the integral over th…
▽ More
For an equation of state in which pressure is a function only of density, the analysis of Newtonian stellar structure is simple in principle if the system is axisymmetric, or consists of a corotating binary. It is then required only to solve two equations: one stating that the "injection energy", $κ$, a potential, is constant throughout the stellar fluid, and the other being the integral over the stellar fluid to give the gravitational potential. An iterative solution of these equations generally diverges if $κ$ is held fixed, but converges with other choices. We investigate the mathematical reason for this convergence/divergence by starting the iteration from an approximation that is perturbatively different from the actual solution. A cycle of iteration is then treated as a linear "updating" operator, and the properties of the linear operator, especially its spectrum, determine the convergence properties. For simplicity, we confine ourselves to spherically symmetric models in which we analyze updating operators both in the finite dimensional space corresponding to a finite difference representation of the problem, and in the continuum, and we find that the fixed-$κ$ operator is self-adjoint and generally has an eigenvalue greater than unity; in the particularly important case of a polytropic equation of state with index greater than unity, we prove that there must be such an eigenvalue. For fixed central density, on the other hand, we find that the updating operator has only a single eigenvector, with zero eigenvalue, and is nilpotent in finite dimension, thereby giving a convergent solution.
△ Less
Submitted 17 March, 2009;
originally announced March 2009.
-
Measuring the neutron star equation of state with gravitational wave observations
Authors:
Jocelyn S. Read,
Charalampos Markakis,
Masaru Shibata,
Koji Uryu,
Jolien D. E. Creighton,
John L. Friedman
Abstract:
We report the results of a first study that uses numerical simulations to estimate the accuracy with which one can use gravitational wave observations of double neutron star inspiral to measure parameters of the neutron-star equation of state. The simulations use the evolution and initial-data codes of Shibata and Uryu to compute the last several orbits and the merger of neutron stars, with matt…
▽ More
We report the results of a first study that uses numerical simulations to estimate the accuracy with which one can use gravitational wave observations of double neutron star inspiral to measure parameters of the neutron-star equation of state. The simulations use the evolution and initial-data codes of Shibata and Uryu to compute the last several orbits and the merger of neutron stars, with matter described by a parametrized equation of state. Previous work suggested the use of an effective cutoff frequency to place constraints on the equation of state. We find, however, that greater accuracy is obtained by measuring departures from the point-particle limit of the gravitational waveform produced during the late inspiral.
As the stars approach their final plunge and merger, the gravitational wave phase accumulates more rapidly for smaller values of the neutron star compactness (the ratio of the mass of the neutron star to its radius). We estimate that realistic equations of state will lead to gravitational waveforms that are distinguishable from point particle inspirals at an effective distance (the distance to an optimally oriented and located system that would produce an equivalent waveform amplitude) of 100 Mpc or less. As Lattimer and Prakash observed, neutron-star radius is closely tied to the pressure at density not far above nuclear. Our results suggest that broadband gravitational wave observations at frequencies between 500 and 1000 Hz will constrain this pressure, and we estimate the accuracy with which it can be measured. Related first estimates of radius measurability show that the radius can be determined to an accuracy of ~1 km at 100 Mpc.
△ Less
Submitted 21 January, 2009;
originally announced January 2009.
-
Constraints on a phenomenologically parameterized neutron-star equation of state
Authors:
Jocelyn S. Read,
Benjamin D. Lackey,
Benjamin J. Owen,
John L. Friedman
Abstract:
We introduce a parameterized high-density equation of state (EOS) in order to systematize the study of constraints placed by astrophysical observations on the nature of neutron-star matter. To obtain useful constraints, the number of parameters should be smaller than the number of neutron-star properties that have been measured or will have been measured in the next several years. And the set mu…
▽ More
We introduce a parameterized high-density equation of state (EOS) in order to systematize the study of constraints placed by astrophysical observations on the nature of neutron-star matter. To obtain useful constraints, the number of parameters should be smaller than the number of neutron-star properties that have been measured or will have been measured in the next several years. And the set must be large enough to accurately approximate the large set of candidate EOSs. We find that a parameterized EOS based on piecewise polytropes with 3 free parameters matches to about 4% rms error an extensive set of candidate EOSs at densities below the central density of 1.4 solar mass stars. Adding observations of more massive stars constrains the higher density part of the EOS and requires an additional parameter. We obtain constraints on the allowed parameter space set by causality and by present and near-future astronomical observations. In particular, we emphasize potentially stringent constraints on the EOS parameter space associated with two measured properties of a single star; and we find that a measurement of the moment of inertia of PSR J0737-3039A can strongly constrain the maximum neutron-star mass. We also present in an appendix a more efficient algorithm than has previously been used for finding points of marginal stability and the maximum angular velocity of stable stars.
△ Less
Submitted 11 December, 2008;
originally announced December 2008.
-
Binary neutron stars: Equilibrium models beyond spatial conformal flatness
Authors:
Koji Uryu,
Francois Limousin,
John L. Friedman,
Eric Gourgoulhon,
Masaru Shibata
Abstract:
Equilibria of binary neutron stars in close circular orbits are computed numerically in a waveless formulation: The full Einstein-relativistic-Euler system is solved on an initial hypersurface to obtain an asymptotically flat form of the 4-metric and an extrinsic curvature whose time derivative vanishes in a comoving frame. Two independent numerical codes are developed, and solution sequences th…
▽ More
Equilibria of binary neutron stars in close circular orbits are computed numerically in a waveless formulation: The full Einstein-relativistic-Euler system is solved on an initial hypersurface to obtain an asymptotically flat form of the 4-metric and an extrinsic curvature whose time derivative vanishes in a comoving frame. Two independent numerical codes are developed, and solution sequences that model inspiraling binary neutron stars during the final several orbits are successfully computed. The binding energy of the system near its final orbit deviates from earlier results of third post-Newtonian and of spatially conformally flat calculations. The new solutions may serve as initial data for merger simulations and as members of quasiequilibrium sequences to generate gravitational wave templates, and may improve estimates of the gravitational-wave cutoff frequency set by the last inspiral orbit.
△ Less
Submitted 21 October, 2006; v1 submitted 25 November, 2005;
originally announced November 2005.
-
Post-Minkowski action for point particles and a helically symmetric binary solution
Authors:
John L. Friedman,
Koji Uryu
Abstract:
Two Fokker actions and corresponding equations of motion are obtained for two point particles in a post-Minkowski framework, in which the field of each particle is given by the half-retarded + half-advanced solution to the linearized Einstein equations. The first action is parametrization invariant, the second a generalization of the affinely parametrized quadratic action for a relativistic part…
▽ More
Two Fokker actions and corresponding equations of motion are obtained for two point particles in a post-Minkowski framework, in which the field of each particle is given by the half-retarded + half-advanced solution to the linearized Einstein equations. The first action is parametrization invariant, the second a generalization of the affinely parametrized quadratic action for a relativistic particle. Expressions for a conserved 4-momentum and angular momentum tensor are obtained in terms of the particles' trajectories in this post-Minkowski approximation. A formal solution to the equations of motion is found for a binary system with circular orbits. For a bound system of this kind, the post-Minkowski solution is a toy model that omits nonlinear terms of relevant post-Newtonian order; and we also obtain a Fokker action that is accurate to first post-Newtonian order, by adding to the post-Minkowski action a term cubic in the particle masses. Curiously, the conserved energy and angular momentum associated with the Fokker action are each finite and well-defined for this bound 2-particle system despite the fact that the total energy and angular momentum of the radiation field diverge. Corresponding solutions and conserved quantities are found for two scalar charges (for electromagnetic charges we exhibit the solution found by Schild). For a broad class of parametrization-invariant Fokker actions and for the affinely parametrized action, binary systems with circular orbits satisfy the relation $dE = ΩdL$ (a form of the first law of thermodynamics), relating the energy, angular velocity and angular momentum of nearby equilibrium configurations.
△ Less
Submitted 5 May, 2006; v1 submitted 1 October, 2005;
originally announced October 2005.
-
The rotational modes of relativistic stars: Numerical results
Authors:
K. H. Lockitch,
J. L. Friedman,
N. Andersson
Abstract:
We study the inertial modes of slowly rotating, fully relativistic compact stars. The equations that govern perturbations of both barotropic and non-barotropic models are discussed, but we present numerical results only for the barotropic case. For barotropic stars all inertial modes are a hybrid mixture of axial and polar perturbations. We use a spectral method to solve for such modes of variou…
▽ More
We study the inertial modes of slowly rotating, fully relativistic compact stars. The equations that govern perturbations of both barotropic and non-barotropic models are discussed, but we present numerical results only for the barotropic case. For barotropic stars all inertial modes are a hybrid mixture of axial and polar perturbations. We use a spectral method to solve for such modes of various polytropic models. Our main attention is on modes that can be driven unstable by the emission of gravitational waves. Hence, we calculate the gravitational-wave growth timescale for these unstable modes and compare the results to previous estimates obtained in Newtonian gravity (i.e. using post-Newtonian radiation formulas). We find that the inertial modes are slightly stabilized by relativistic effects, but that previous conclusions concerning eg. the unstable r-modes remain essentially unaltered when the problem is studied in full general relativity.
△ Less
Submitted 29 October, 2002;
originally announced October 2002.
-
Thermodynamics of binary black holes and neutron stars
Authors:
John L. Friedman,
Koji Uryu,
Masaru Shibata
Abstract:
We consider compact binary systems, modeled in general relativity as vacuum or perfect-fluid spacetimes with a helical Killing vector k^α, heuristically, the generator of time-translations in a corotating frame. Systems that are stationary in this sense are not asymptotically flat, but have asymptotic behavior corresponding to equal amounts of ingoing and outgoing radiation. For black-hole binar…
▽ More
We consider compact binary systems, modeled in general relativity as vacuum or perfect-fluid spacetimes with a helical Killing vector k^α, heuristically, the generator of time-translations in a corotating frame. Systems that are stationary in this sense are not asymptotically flat, but have asymptotic behavior corresponding to equal amounts of ingoing and outgoing radiation. For black-hole binaries, a rigidity theorem implies that the Killing vector lies along the horizon's generators, and from this one can deduce the zeroth law (constant surface gravity of the horizon). Remarkably, although the mass and angular momentum of such a system are not defined, there is an exact first law, relating the change in the asymptotic Noether charge to the changes in the vorticity, baryon mass, and entropy of the fluid, and in the area of black holes.
Binary systems with MΩsmall have an approximate asymptopia in which one can write the first law in terms of the asymptotic mass and angular momentum. Asymptotic flatness is precise in two classes of solutions used to model binary systems: spacetimes satisfying the post-Newtonian equations, and solutions to a modified set of field equations that have a spatially conformally flat metric. (The spatial conformal flatness formalism with helical symmetry, however, is consistent with maximal slicing only if replaces the extrinsic curvature in the field equations by an artificially tracefree expression in terms of the shift vector.) For these spacetimes, nearby equilibria whose stars have the same vorticity obey the relation δM = ΩδJ, from which one can obtain a turning point criterion that governs the stability of orbits.
△ Less
Submitted 12 December, 2001; v1 submitted 28 August, 2001;
originally announced August 2001.
-
Hydrostatic Expansion and Spin Changes During Type I X-Ray Bursts
Authors:
Andrew Cumming,
Sharon M. Morsink,
Lars Bildsten,
John L. Friedman,
Daniel E. Holz
Abstract:
We present calculations of the spin-down of a neutron star atmosphere due to hydrostatic expansion during a Type I X-ray burst. We show that (i) Cumming and Bildsten overestimated the spin-down of rigidly-rotating atmospheres by a factor of two, and (ii) general relativity has a small (5-10%) effect on the angular momentum conservation law. We rescale our results to different neutron star masses…
▽ More
We present calculations of the spin-down of a neutron star atmosphere due to hydrostatic expansion during a Type I X-ray burst. We show that (i) Cumming and Bildsten overestimated the spin-down of rigidly-rotating atmospheres by a factor of two, and (ii) general relativity has a small (5-10%) effect on the angular momentum conservation law. We rescale our results to different neutron star masses, rotation rates and equations of state, and present some detailed rotational profiles. Comparing with recent observations of large frequency shifts in MXB 1658-298 and 4U 1916-053, we find that the spin-down expected if the atmosphere rotates rigidly is a factor of two to three less than the observed values. If differential rotation is allowed to persist, we find that the upper layers of the atmosphere spin down by an amount comparable to the observed values; however, there is no compelling reason to expect the observed spin frequency to be that of only the outermost layers. We conclude that hydrostatic expansion and angular momentum conservation alone cannot account for the largest frequency shifts observed during Type I bursts.
△ Less
Submitted 1 August, 2001;
originally announced August 2001.
-
Implications of the r-mode instability of rotating relativistic stars
Authors:
John L. Friedman,
Keith H. Lockitch
Abstract:
Several recent surprises appear dramatically to have improved the likelihood that the spin of rapidly rotating, newly formed neutron stars (and, possibly, of old stars spun up by accretion) is limited by a nonaxisymmetric instability driven by gravitational waves. Except for the earliest part of the spin-down, the axial l=m=2 mode (an r-mode) dominates the instability, and the emitted waves may…
▽ More
Several recent surprises appear dramatically to have improved the likelihood that the spin of rapidly rotating, newly formed neutron stars (and, possibly, of old stars spun up by accretion) is limited by a nonaxisymmetric instability driven by gravitational waves. Except for the earliest part of the spin-down, the axial l=m=2 mode (an r-mode) dominates the instability, and the emitted waves may be observable by detectors with the sensitivity of LIGO II. A review of these hopeful results is followed by a discussion of constraints on the instability set by dissipative mechanisms, including viscosity, nonlinear saturation, and energy loss to a magnetic field driven by differential rotation.
△ Less
Submitted 1 March, 2001; v1 submitted 28 February, 2001;
originally announced February 2001.
-
The rotational modes of relativistic stars: Analytic results
Authors:
Keith H. Lockitch,
Nils Andersson,
John L. Friedman
Abstract:
We study the r-modes and rotational ``hybrid'' modes of relativistic stars. As in Newtonian gravity, the spectrum of low-frequency rotational modes is highly sensitive to the stellar equation of state. If the star and its perturbations obey the same one-parameter equation of state (as with isentropic stars), there exist {\it no pure r-modes at all} - no modes whose limit, for a star with zero an…
▽ More
We study the r-modes and rotational ``hybrid'' modes of relativistic stars. As in Newtonian gravity, the spectrum of low-frequency rotational modes is highly sensitive to the stellar equation of state. If the star and its perturbations obey the same one-parameter equation of state (as with isentropic stars), there exist {\it no pure r-modes at all} - no modes whose limit, for a star with zero angular velocity, is an axial-parity perturbation. Rotating stars of this kind similarly have no pure g-modes, no modes whose spherical limit is a perturbation with polar parity and vanishing perturbed pressure and density. We compute the post-Newtonian corrections to the $l=m$ r-modes of isentropic and non-isentropic uniform density stars.
△ Less
Submitted 8 August, 2000;
originally announced August 2000.
-
Gravitational-wave driven instability of rotating relativistic stars
Authors:
John L. Friedman,
Keith H. Lockitch
Abstract:
A brief review of the stability of rotating relativistic stars is followed by a more detailed discussion of recent work on an instability of r-modes, modes of rotating stars that have axial parity in the slow-rotation limit. These modes may dominate the spin-down of neutron stars that are rapidly rotating at birth, and the gravitational waves they emit may be detectable.
A brief review of the stability of rotating relativistic stars is followed by a more detailed discussion of recent work on an instability of r-modes, modes of rotating stars that have axial parity in the slow-rotation limit. These modes may dominate the spin-down of neutron stars that are rapidly rotating at birth, and the gravitational waves they emit may be detectable.
△ Less
Submitted 31 August, 1999;
originally announced August 1999.
-
Where are the r-modes of isentropic stars?
Authors:
Keith H. Lockitch,
John L. Friedman
Abstract:
Almost none of the r-modes ordinarily found in rotating stars exist, if the star and its perturbations obey the same one-parameter equation of state; and rotating relativistic stars with one-parameter equations of state have no pure r-modes at all, no modes whose limit, for a star with zero angular velocity, is a perturbation with axial parity. Similarly (as we show here) rotating stars of this…
▽ More
Almost none of the r-modes ordinarily found in rotating stars exist, if the star and its perturbations obey the same one-parameter equation of state; and rotating relativistic stars with one-parameter equations of state have no pure r-modes at all, no modes whose limit, for a star with zero angular velocity, is a perturbation with axial parity. Similarly (as we show here) rotating stars of this kind have no pure g-modes, no modes whose spherical limit is a perturbation with polar parity and vanishing perturbed pressure and density. Where have these modes gone?
In spherical stars of this kind, r-modes and g-modes form a degenerate zero-frequency subspace. We find that rotation splits the degeneracy to zeroth order in the star's angular velocity $Ω$, and the resulting modes are generically hybrids, whose limit as $Ω\to 0$ is a stationary current with axial and polar parts. Because each mode has definite parity, its axial and polar parts have alternating values of $l$. We show that each mode belongs to one of two classes, axial-led or polar-led, depending on whether the spherical harmonic with lowest value of $l$ that contributes to its velocity field is axial or polar. We numerically compute these modes for slowly rotating polytropes and for Maclaurin spheroids, using a straightforward method that appears to be novel and robust. Timescales for the gravitational-wave driven instability and for viscous damping are computed using assumptions appropriate to neutron stars.
△ Less
Submitted 19 March, 1999; v1 submitted 4 December, 1998;
originally announced December 1998.
-
Axial instability of rotating relativistic stars
Authors:
John L. Friedman,
Sharon M. Morsink
Abstract:
Perturbations of rotating relativistic stars can be classified by their behavior under parity. For axial perturbations (r-modes), initial data with negative canonical energy is found with angular dependence $e^{imφ}$ for all values of $m\geq 2$ and for arbitrarily slow rotation. This implies instability (or marginal stability) of such perturbations for rotating perfect fluids. This low $m$-insta…
▽ More
Perturbations of rotating relativistic stars can be classified by their behavior under parity. For axial perturbations (r-modes), initial data with negative canonical energy is found with angular dependence $e^{imφ}$ for all values of $m\geq 2$ and for arbitrarily slow rotation. This implies instability (or marginal stability) of such perturbations for rotating perfect fluids. This low $m$-instability is strikingly different from the instability to polar perturbations, which sets in first for large values of $m$. The timescale for the axial instability appears, for small angular velocity $Ω$, to be proportional to a high power of $Ω$. As in the case of polar modes, viscosity will again presumably enforce stability except for hot, rapidly rotating neutron stars. This work complements Andersson's numerical investigation of axial modes in slowly rotating stars.
△ Less
Submitted 11 June, 1998; v1 submitted 23 June, 1997;
originally announced June 1997.
-
Nonaxisymmetric Neutral Modes in Rotating Relativistic Stars
Authors:
Nikolaos Stergioulas,
John L. Friedman
Abstract:
We study nonaxisymmetric perturbations of rotating relativistic stars. modeled as perfect-fluid equilibria. Instability to a mode with angular dependence $\exp(imφ)$ sets in when the frequency of the mode vanishes. The locations of these zero-frequency modes along sequences of rotating stars are computed in the framework of general relativity. We consider models of uniformly rotating stars with…
▽ More
We study nonaxisymmetric perturbations of rotating relativistic stars. modeled as perfect-fluid equilibria. Instability to a mode with angular dependence $\exp(imφ)$ sets in when the frequency of the mode vanishes. The locations of these zero-frequency modes along sequences of rotating stars are computed in the framework of general relativity. We consider models of uniformly rotating stars with polytropic equations of state, finding that the relativistic models are unstable to nonaxisymmetric modes at significantly smaller values of rotation than in the Newtonian limit. Most strikingly, the m=2 bar mode can become unstable even for soft polytropes of index $N \leq 1.3$, while in Newtonian theory it becomes unstable only for stiff polytropes of index $N \leq 0.808$. If rapidly rotating neutron stars are formed by the accretion-induced collapse of white dwarfs, instability associated with these nonaxisymmetric, gravitational-wave driven modes may set an upper limit on neutron-star rotation. Consideration is restricted to perturbations that correspond to polar perturbations of a spherical star. A study of axial perturbations is in progress.
△ Less
Submitted 21 May, 1997; v1 submitted 21 May, 1997;
originally announced May 1997.
-
Upper Limit Set by Causality on the Rotation and Mass of Uniformly Rotating Relativistic Stars
Authors:
Scott Koranda,
Nikolaos Stergioulas,
John L. Friedman
Abstract:
Causality alone suffices to set a lower bound on the period of rotation of relativistic stars as a function of their maximum observed mass. That is, by assuming a one-parameter equation of state (EOS) that satisfies v_sound < c and that allows stars with masses as large as the largest observed neutron-star mass, M_sph^max, we find P[ms] > 0.282 + 0.196 ( M_sph^max/M_odot-1.442). The limit does n…
▽ More
Causality alone suffices to set a lower bound on the period of rotation of relativistic stars as a function of their maximum observed mass. That is, by assuming a one-parameter equation of state (EOS) that satisfies v_sound < c and that allows stars with masses as large as the largest observed neutron-star mass, M_sph^max, we find P[ms] > 0.282 + 0.196 ( M_sph^max/M_odot-1.442). The limit does not assume that the EOS agrees with a known low-density form for ordinary matter, but if one adds that assumption, the minimum period is raised by a few percent. Thus the current minimum period of uniformly rotating stars, set by causality, is 0.28ms (0.29ms for stars with normal crust). The minimizing EOS yields models with a maximally soft exterior supported by a maximally stiff core. An analogous upper limit set by causality on the maximum mass of rotating neutron stars requires a low-density match and the limit depends on the matching density, epsilon_m. We recompute it, obtaining a slightly revised value, M_rot^max \sim 6.1( 2 * 10^14 g/cm^3 epsilon_m )^1/2 M_odot.
△ Less
Submitted 27 August, 1996;
originally announced August 1996.
-
Comparing Models of Rapidly Rotating Relativistic Stars Constructed by Two Numerical Methods
Authors:
Nikolaos Stergioulas,
John L. Friedman
Abstract:
We present the first direct comparison of codes based on two different numerical methods for constructing rapidly rotating relativistic stars. A code based on the Komatsu-Eriguchi-Hachisu (KEH) method (Komatsu et al. 1989), written by Stergioulas, is compared to the Butterworth-Ipser code (BI), as modified by Friedman, Ipser and Parker. We compare models obtained by each method and evaluate the…
▽ More
We present the first direct comparison of codes based on two different numerical methods for constructing rapidly rotating relativistic stars. A code based on the Komatsu-Eriguchi-Hachisu (KEH) method (Komatsu et al. 1989), written by Stergioulas, is compared to the Butterworth-Ipser code (BI), as modified by Friedman, Ipser and Parker. We compare models obtained by each method and evaluate the accuracy and efficiency of the two codes. The agreement is surprisingly good. A relatively large discrepancy recently reported (Eriguchi et al. 1994) is found to arise from the use of two different versions of the equation of state. We find, for a given equation of state, that equilibrium models with maximum values of mass, baryon mass, and angular momentum are (generically) all distinct and either all unstable to collapse or are all stable. Our implementation of the KEH method will be available as a public domain program for interested users.
△ Less
Submitted 8 November, 1994;
originally announced November 1994.