Search

Neutron stars provide a unique opportunity to study strongly interacting matter under extreme density conditions. The intricacies of matter inside neutron stars and their equation of state are not directly visible, but determine bulk properties, such as mass and radius, which affect the star's thermal X-ray emissions. However, the telescope spectra of these emissions are also affected by the stellar distance, hydrogen column, and effective surface temperature, which are not always well-constrained. Uncertainties on these nuisance parameters must be accounted for when making a robust estimation of the equation of state. In this study, we develop a novel methodology that, for the first time, can infer the full posterior distribution of both the equation of state and nuisance parameters directly from telescope observations. This method relies on the use of neural likelihood estimation, in which normalizing flows use samples of simulated telescope data to learn the likelihood of the neutron star spectra as a function of these parameters, coupled with Hamiltonian Monte Carlo methods to efficiently sample from the corresponding posterior distribution. Our approach surpasses the accuracy of previous methods, improves the interpretability of the results by providing access to the full posterior distribution, and naturally scales to a growing number of neutron star observations expected in the coming years.

Neutron stars provide a unique opportunity to study strongly interacting matter under extreme density conditions. The intricacies of matter inside neutron stars and their equation of state are not directly visible, but determine bulk properties, such as mass and radius, which affect the star's thermal X-ray emissions. However, the telescope spectra of these emissions are also affected by the stellar distance, hydrogen column, and effective surface temperature, which are not always well-constrained. Uncertainties on these nuisance parameters must be accounted for when making a robust estimation of the equation of state. In this study, we develop a novel methodology that, for the first time, can infer the full posterior distribution of both the equation of state and nuisance parameters directly from telescope observations. This method relies on the use of neural likelihood estimation, in which normalizing flows use samples of simulated telescope data to learn the likelihood of the neutron star spectra as a function of these parameters, coupled with Hamiltonian Monte Carlo methods to efficiently sample from the corresponding posterior distribution. Our approach surpasses the accuracy of previous methods, improves the interpretability of the results by providing access to the full posterior distribution, and naturally scales to a growing number of neutron star observations expected in the coming years.

Neutron stars provide a unique opportunity to study strongly interacting matter under extreme density conditions. The intricacies of matter inside neutron stars and their equation of state are not directly visible, but determine bulk properties, such as mass and radius, which affect the star's thermal X-ray emissions. However, the telescope spectra of these emissions are also affected by the stellar distance, hydrogen column, and effective surface temperature, which are not always well-constrained. Uncertainties on these nuisance parameters must be accounted for when making a robust estimation of the equation of state. In this study, we develop a novel methodology that, for the first time, can infer the full posterior distribution of both the equation of state and nuisance parameters directly from telescope observations. This method relies on the use of neural likelihood estimation, in which normalizing flows use samples of simulated telescope data to learn the likelihood of the neutron star spectra as a function of these parameters, coupled with Hamiltonian Monte Carlo methods to efficiently sample from the corresponding posterior distribution. Our approach surpasses the accuracy of previous methods, improves the interpretability of the results by providing access to the full posterior distribution, and naturally scales to a growing number of neutron star observations expected in the coming years.

The interiors of neutron stars reach densities and temperatures beyond the limits of terrestrial experiments, providing vital laboratories for probing nuclear physics. While the star's interior is not directly observable, its pressure and density determine the star's macroscopic structure which affects the spectra observed in telescopes. The relationship between the observations and the internal state is complex and partially intractable, presenting difficulties for inference. Previous work has focused on the regression from stellar spectra of parameters describing the internal state. We demonstrate a calculation of the full likelihood of the internal state parameters given observations, accomplished by replacing intractable elements with machine learning models trained on samples of simulated stars. Our machine-learning-derived likelihood allows us to perform maximum a posteriori estimation of the parameters of interest, as well as full scans. We demonstrate the technique by inferring stellar mass and radius from an individual stellar spectrum, as well as equation of state parameters from a set of spectra. Our results are more precise than pure regression models, reducing the width of the parameter residuals by 11.8% in the most realistic scenario. The neural networks will be released as a tool for fast simulation of neutron star properties and observed spectra.

The interiors of neutron stars reach densities and temperatures beyond the limits of terrestrial experiments, providing vital laboratories for probing nuclear physics. While the star's interior is not directly observable, its pressure and density determine the star's macroscopic structure which affects the spectra observed in telescopes. The relationship between the observations and the internal state is complex and partially intractable, presenting difficulties for inference. Previous work has focused on the regression from stellar spectra of parameters describing the internal state. We demonstrate a calculation of the full likelihood of the internal state parameters given observations, accomplished by replacing intractable elements with machine learning models trained on samples of simulated stars. Our machine-learning-derived likelihood allows us to perform maximum a posteriori estimation of the parameters of interest, as well as full scans. We demonstrate the technique by inferring stellar mass and radius from an individual stellar spectrum, as well as equation of state parameters from a set of spectra. Our results are more precise than pure regression models, reducing the width of the parameter residuals by 11.8% in the most realistic scenario. The neural networks will be released as a tool for fast simulation of neutron star properties and observed spectra., Comment: 18 pages, 8 figures

We report on a new search for continuous gravitational waves from NS 1987A, the neutron star born in SN 1987A, using open data from Advanced LIGO and Virgo's third observing run (O3). The search covered frequencies from 35-1050 Hz, more than five times the band of the only previous gravitational wave search to constrain NS 1987A [B. J. Owen et al., ApJL 935, L7 (2022)]. It used an improved code and coherently integrated from 5.10 days to 14.85 days depending on frequency. No astrophysical signals were detected. By expanding the frequency range and using O3 data, this search improved on strain upper limits from the previous search and was sensitive at the highest frequencies to ellipticities of 1.6e-5 and r-mode amplitudes of 4.4e-4, both an order of magnitude improvement over the previous search and both well within the range of theoretical predictions., Comment: 7 pages, 1 figure

Spectral representations have been shown to provide an efficient way to represent the poorly understood high-density portion of the neutron-star equation of state. This paper shows how the efficiency and accuracy of those representations can be improved by a very simple change., Comment: 5 pages, 6 figures

Neutron stars provide a unique laboratory for studying matter at extreme pressures and densities. While there is no direct way to explore their interior structure, X-rays emitted from these stars can indirectly provide clues to the equation of state (EOS) of superdense nuclear matter through the inference of the star's mass and radius. However, inference of EOS directly from a star's X-ray spectra is extremely challenging and is complicated by systematic uncertainties. The current state of the art is to use simulation-based likelihoods in a piece-wise method, which first infer the star's mass and radius to reduce the dimensionality of the problem, and from those quantities infer the EOS. We demonstrate a series of enhancements to the state of the art, in terms of realistic uncertainty quantification and improved regression of physical properties with machine learning. We also demonstrate novel inference of the EOS directly from the high-dimensional spectra of observed stars, avoiding the intermediate mass-radius step. Our network is conditioned on the sources of uncertainty of each star, allowing for natural and complete propagation of uncertainties to the EOS., Comment: 24 pages, 20 figures

We report on a search for continuous gravitational waves (GWs) from NS 1987A, the neutron star born in SN 1987A. The search covered a frequency band of 75-275 Hz, included a wide range of spin-down parameters for the first time, and coherently integrated 12.8 days of data below 125 Hz and 8.7 days of data above 125 Hz from the second Advanced LIGO observing run. We found no astrophysical signal. We set upper limits on GW emission as tight as an intrinsic strain of $2\times10^{-25}$ at 90\% confidence. The large spin-down parameter space makes this search the first astrophysically consistent one for continuous GWs from NS 1987A. Our upper limits are the first consistent ones to beat an analog of the spin-down limit based on the age of the neutron star, and hence are the first GW observations to put new constraints on NS 1987A.

Spectral representations have been shown to provide an efficient way to represent the poorly understood high-density portion of the neutron-star equation of state. This paper shows how the efficiency and accuracy of those representations can be improved by a very simple change., Comment: 5 pages, 6 figures

A method is developed here for building differentiable three-dimensional manifolds on multicube structures. This method constructs a sequence of reference metrics that determine differentiable structures on the cubic regions that serve as non-overlapping coordinate charts on these manifolds. It uses solutions to the two- and three-dimensional biharmonic equations in a sequence of steps that increase the differentiability of the reference metrics across the interfaces between cubic regions. This method is algorithmic and has been implemented in a computer code that automatically generates these reference metrics. Examples of three-manifolds constructed in this way are presented here, including representatives from five of the eight Thurston geometrization classes, plus the well-known Hantzsche-Wendt, the Poincare dodecahedral space, and the Seifert-Weber space., Comment: 37 pages, 10 figures, 12 Tables; v2 includes minor revisions and additions to agree with the final published version

We describe directed searches for continuous gravitational waves from twelve well localized non-pulsing candidate neutron stars in young supernova remnants using data from Advanced LIGO's second observing run. We assumed that each neutron star is isolated and searched a band of frequencies from 15 to 150\,Hz, consistent with frequencies expected from known young pulsars. After coherently integrating spans of data ranging from 12.0 to 55.9 days using the F-statistic and applying data-based vetoes, we found no evidence of astrophysical signals. We set upper limits on intrinsic gravitational wave amplitude in some cases stronger than 10^{-25}, generally about a factor of two better than upper limits on the same objects from Advanced LIGO's first observing run., Comment: v2 updated to version accepted for publication in Physical Review D

We describe directed searches for continuous gravitational waves from twelve well localized non-pulsing candidate neutron stars in young supernova remnants using data from Advanced LIGO's second observing run. We assumed that each neutron star is isolated and searched a band of frequencies from 15 to 150\,Hz, consistent with frequencies expected from known young pulsars. After coherently integrating spans of data ranging from 12.0 to 55.9 days using the F-statistic and applying data-based vetoes, we found no evidence of astrophysical signals. We set upper limits on intrinsic gravitational wave amplitude in some cases stronger than 10^{-25}, generally about a factor of two better than upper limits on the same objects from Advanced LIGO's first observing run., Comment: v2 updated to version accepted for publication in Physical Review D

Parameterized representations of the equation of state play an important role in efforts to measure the properties of the matter in the cores of neutron stars using astronomical observations. New representations are presented here that are capable of representing any equation of state to any desired level of accuracy, while automatically imposing causality and thermodynamic stability constraints. Numerical tests are presented that measure how accurately and efficiently these new parameterizations represent a collection of causal nuclear-theory model equations of state., Comment: 10 pages, 8 figures; v2 updated with minor modifications to form accepted for publication in Phys. Rev. D

Gravitational wave detectors in the LIGO/Virgo frequency band are able to measure the individual masses and the composite tidal deformabilities of neutron-star binary systems. This paper demonstrates that high accuracy measurements of these quantities from an ensemble of binary systems can in principle be used to determine the high density neutron-star equation of state exactly. This analysis assumes that all neutron stars have the same thermodynamically stable equation of state, but does not use simplifying approximations for the composite tidal deformability or make additional assumptions about the high density equation of state., Comment: 6 pages, 4 figures; v2 updated to version accepted for publication in Phys. Rev. D

Gravitational wave detectors in the LIGO/Virgo frequency band are able to measure the individual masses and the composite tidal deformabilities of neutron-star binary systems. This paper demonstrates that high accuracy measurements of these quantities from an ensemble of binary systems can in principle be used to determine the high density neutron-star equation of state exactly. This analysis assumes that all neutron stars have the same thermodynamically stable equation of state, but does not use simplifying approximations for the composite tidal deformability or make additional assumptions about the high density equation of state., Comment: 6 pages, 4 figures; v2 updated to version accepted for publication in Phys. Rev. D

Parameterized representations of the equation of state play an important role in efforts to measure the properties of the matter in the cores of neutron stars using astronomical observations. New representations are presented here that are capable of representing any equation of state to any desired level of accuracy, while automatically imposing causality and thermodynamic stability constraints. Numerical tests are presented that measure how accurately and efficiently these new parameterizations represent a collection of causal nuclear-theory model equations of state., Comment: 10 pages, 8 figures; v2 updated with minor modifications to form accepted for publication in Phys. Rev. D

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 $\alpha$, 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, $\alpha \ll 1$, predicted by current models., Comment: 5 pages; Contribution to the proceedings of the Physics of Neutron Stars - 2017 (Saint Petersburg, July 10-14, 2017), associated with a talk by A.I. Chugunov

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\'en 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., Comment: This is essentially the published version

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 $\alpha$, 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, $\alpha \ll 1$, predicted by current models., Comment: 5 pages; Contribution to the proceedings of the Physics of Neutron Stars - 2017 (Saint Petersburg, July 10-14, 2017), associated with a talk by A.I. Chugunov

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\'en 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., Comment: This is essentially the published version

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\'a. 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., Comment: Published (2016) version, but with two updated references

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\'a. 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., Comment: Published (2016) version, but with two updated references

The observable macroscopic properties of relativistic stars (whose equations of state are known) can be predicted by solving the stellar structure equations that follow from Einstein's equation. For neutron stars, however, our knowledge of the equation of state is poor, so the direct stellar structure problem can not be solved without modeling the highest density part of the equation of state in some way. This talk will describe recent work on developing a model independent approach to determining the high-density neutron-star equation of state by solving an inverse stellar structure problem. This method uses the fact that Einstein's equation provides a deterministic relationship between the equation of state and the macroscopic observables of the stars which are composed of that material. This talk illustrates how this method will be able to determine the high-density part of the neutron-star equation of state with few percent accuracy when high quality measurements of the masses and radii of just two or three neutron stars become available. This talk will also show that this method can be used with measurements of other macroscopic observables, like the masses and tidal deformabilities, which can (in principle) be measured by gravitational wave observations of binary neutron-star mergers., Comment: Talk given at the Fifth Leopoldo Garcia-Colin Mexican Meeting on Mathematical and Experimental Physics

The observable macroscopic properties of relativistic stars (whose equations of state are known) can be predicted by solving the stellar structure equations that follow from Einstein’s equation. For neutron stars, however, our knowledge of the equation of state is poor, so the direct stellar structure problem can not be solved without modeling the highest density part of the equation of state in some way. This talk will describe recent work on developing a model independent approach to determining the high-density neutron-star equation of state by solving an inverse stellar structure problem. This method uses the fact that Einstein’s equation provides a deterministic relationship between the equation of state and the macroscopic observables of the stars which are composed of that material. This talk illustrates how this method will be able to determine the high-density part of the neutron-star equation of state with few percent accuracy when high quality measurements of the masses and radii of just two or three neutron stars become available. This talk will also show that this method can be used with measurements of other macroscopic observables, like the masses and tidal deformabilities, which can (in principle) be measured by gravitational wave observations of binary neutron-star mergers.

The observable macroscopic properties of relativistic stars (whose equations of state are known) can be predicted by solving the stellar structure equations that follow from Einstein's equation. For neutron stars, however, our knowledge of the equation of state is poor, so the direct stellar structure problem can not be solved without modeling the highest density part of the equation of state in some way. This talk will describe recent work on developing a model independent approach to determining the high-density neutron-star equation of state by solving an inverse stellar structure problem. This method uses the fact that Einstein's equation provides a deterministic relationship between the equation of state and the macroscopic observables of the stars which are composed of that material. This talk illustrates how this method will be able to determine the high-density part of the neutron-star equation of state with few percent accuracy when high quality measurements of the masses and radii of just two or three neutron stars become available. This talk will also show that this method can be used with measurements of other macroscopic observables, like the masses and tidal deformabilities, which can (in principle) be measured by gravitational wave observations of binary neutron-star mergers., Comment: Talk given at the Fifth Leopoldo Garcia-Colin Mexican Meeting on Mathematical and Experimental Physics

The inverse stellar structure problem determines the equation of state of the matter in stars from a knowledge of their macroscopic observables (e.g. their masses and radii). This problem was solved in a previous paper by constructing a spectral representation of the equation of state whose stellar models match a prescribed set of macroscopic observables. This paper improves and extends that work in two significant ways: i) The method is made more robust by accounting for an unexpected feature of the enthalpy based representations of the equations of state used in this work. After making the appropriate modifications, accurate initial guesses for the spectral parameters are no longer needed so Monte-Carlo techniques can now be used to ensure the best fit to the observables. ii) The method is extended here to use masses and tidal deformabilities (which will be measured by gravitational wave observations of neutron-star mergers) as the macroscopic observables instead of masses and radii. The accuracy and reliability of this extended and more robust spectral method is evaluated in this paper using mock data for observables from stars based on 34 different theoretical models of the high density neutron-star equation of state. In qualitative agreement with earlier work, these tests suggest the high density part of the neutron-star equation of state could be determined at the few-percent accuracy level using high quality measurements of the masses and radii (or masses and tidal deformabilities) of just two or three neutron stars., Comment: 16 pages, 5 figures, 2 tables; v2 updated to published version, includes expanded discussion section, v3 corrected typo in Eq. (C7)

The inverse stellar structure problem determines the equation of state of the matter in stars from a knowledge of their macroscopic observables (e.g. their masses and radii). This problem was solved in a previous paper by constructing a spectral representation of the equation of state whose stellar models match a prescribed set of macroscopic observables. This paper improves and extends that work in two significant ways: i) The method is made more robust by accounting for an unexpected feature of the enthalpy based representations of the equations of state used in this work. After making the appropriate modifications, accurate initial guesses for the spectral parameters are no longer needed so Monte-Carlo techniques can now be used to ensure the best fit to the observables. ii) The method is extended here to use masses and tidal deformabilities (which will be measured by gravitational wave observations of neutron-star mergers) as the macroscopic observables instead of masses and radii. The accuracy and reliability of this extended and more robust spectral method is evaluated in this paper using mock data for observables from stars based on 34 different theoretical models of the high density neutron-star equation of state. In qualitative agreement with earlier work, these tests suggest the high density part of the neutron-star equation of state could be determined at the few-percent accuracy level using high quality measurements of the masses and radii (or masses and tidal deformabilities) of just two or three neutron stars., Comment: 16 pages, 5 figures, 2 tables; v2 updated to published version, includes expanded discussion section, v3 corrected typo in Eq. (C7)

A new method for solving the relativistic inverse stellar structure problem is presented. This method determines a spectral representation of the unknown high density portion of the stellar equation of state from a knowledge of the total masses M and radii R of the stars. Spectral representations of the equation of state are very efficient, generally requiring only a few spectral parameters to achieve good accuracy. This new method is able, therefore, to determine the high density equation of state quite accurately from only a few accurately measured [M,R] data points. This method is tested here by determining the equations of state from mock [M,R] data computed from tabulated "realistic" neutron-star equations of state. The spectral equations of state obtained from these mock data are shown to agree on average with the originals to within a few percent (over the entire high density range of the neutron-star interior) using only two [M,R] data points. Higher accuracies are achieved when more data are used. The accuracies of the equations of state determined in these examples are shown to be nearly optimal, in the sense that their errors are comparable to the errors of the best-fit spectral representations of these realistic equations of state., Comment: 12 pages; 1 table; v2 minor changes to version accepted in Phys. Rev. D

A new method for solving the relativistic inverse stellar structure problem is presented. This method determines a spectral representation of the unknown high density portion of the stellar equation of state from a knowledge of the total masses M and radii R of the stars. Spectral representations of the equation of state are very efficient, generally requiring only a few spectral parameters to achieve good accuracy. This new method is able, therefore, to determine the high density equation of state quite accurately from only a few accurately measured [M,R] data points. This method is tested here by determining the equations of state from mock [M,R] data computed from tabulated "realistic" neutron-star equations of state. The spectral equations of state obtained from these mock data are shown to agree on average with the originals to within a few percent (over the entire high density range of the neutron-star interior) using only two [M,R] data points. Higher accuracies are achieved when more data are used. The accuracies of the equations of state determined in these examples are shown to be nearly optimal, in the sense that their errors are comparable to the errors of the best-fit spectral representations of these realistic equations of state., Comment: 12 pages; 1 table; v2 minor changes to version accepted in Phys. Rev. D

Methods are developed for constructing spectral representations of cold (barotropic) neutron-star equations of state. These representations are faithful in the sense that every physical equation of state has a representation of this type, and conversely every such representation satisfies the minimal thermodynamic stability criteria required of any physical equation of state. These spectral representations are also efficient, in the sense that only a few spectral coefficients are generally required to represent neutron-star equations of state quiet accurately. This accuracy and efficiency is illustrated by constructing spectral fits to a large collection of "realistic" neutron-star equations of state., Comment: 10 pages, 8 figures, 2 tables; updated to accepted Phys. Rev. D version

Methods are developed for constructing spectral representations of cold (barotropic) neutron-star equations of state. These representations are faithful in the sense that every physical equation of state has a representation of this type, and conversely every such representation satisfies the minimal thermodynamic stability criteria required of any physical equation of state. These spectral representations are also efficient, in the sense that only a few spectral coefficients are generally required to represent neutron-star equations of state quiet accurately. This accuracy and efficiency is illustrated by constructing spectral fits to a large collection of "realistic" neutron-star equations of state., Comment: 10 pages, 8 figures, 2 tables; updated to accepted Phys. Rev. D version

The generalized harmonic representation of Einstein's equations is manifestly hyperbolic for a large class of gauge conditions. Unfortunately most of the useful gauges developed over the past several decades by the numerical relativity community are incompatible with the hyperbolicity of the equations in this form. This paper presents a new method of imposing gauge conditions that preserves hyperbolicity for a much wider class of conditions, including as special cases many of the standard ones used in numerical relativity: e.g., K freezing, Gamma freezing, Bona-Massó slicing, conformal Gamma drivers, etc. Analytical and numerical results are presented which test the stability and the effectiveness of this new gauge-driver evolution system.

The generalized harmonic representation of Einstein's equations is manifestly hyperbolic for a large class of gauge conditions. Unfortunately most of the useful gauges developed over the past several decades by the numerical relativity community are incompatible with the hyperbolicity of the equations in this form. This paper presents a new method of imposing gauge conditions that preserves hyperbolicity for a much wider class of conditions, including as special cases many of the standard ones used in numerical relativity: e.g., K freezing, Gamma freezing, Bona-Massó slicing, conformal Gamma drivers, etc. Analytical and numerical results are presented which test the stability and the effectiveness of this new gauge-driver evolution system.

The accuracy and stability of the Caltech-Cornell pseudospectral code is evaluated using the Kidder, Scheel, and Teukolsky (KST) representation of the Einstein evolution equations. The basic “Mexico City tests” widely adopted by the numerical relativity community are adapted here for codes based on spectral methods. Exponential convergence of the spectral code is established, apparently limited only by numerical roundoff error or by truncation error in the time integration. A general expression for the growth of errors due to finite machine precision is derived, and it is shown that this limit is achieved here for the linear plane-wave test.

Binary black hole simulations starting from quasi-circular (i.e., zero radial velocity) initial data have orbits with small but nonzero orbital eccentricities. In this paper, the quasi-equilibrium initial-data method is extended to allow nonzero radial velocities to be specified in binary black hole initial data. New low-eccentricity initial data are obtained by adjusting the orbital frequency and radial velocities to minimize the orbital eccentricity, and the resulting (~5 orbit) evolutions are compared with those of quasi-circular initial data. Evolutions of the quasi-circular data clearly show eccentric orbits, with eccentricity that decays over time. The precise decay rate depends on the definition of eccentricity; if defined in terms of variations in the orbital frequency, the decay rate agrees well with the prediction of Peters (1964 Phys. Rev. 136 1224–32). The gravitational waveforms, which contain ~8 cycles in the dominant l = m = 2 mode, are largely unaffected by the eccentricity of the quasi-circular initial data. The overlap between the dominant mode in the quasi-circular evolution and the same mode in the low-eccentricity evolution is about 0.99.

The accuracy and stability of the Caltech-Cornell pseudospectral code is evaluated using the Kidder, Scheel, and Teukolsky (KST) representation of the Einstein evolution equations. The basic “Mexico City tests” widely adopted by the numerical relativity community are adapted here for codes based on spectral methods. Exponential convergence of the spectral code is established, apparently limited only by numerical roundoff error or by truncation error in the time integration. A general expression for the growth of errors due to finite machine precision is derived, and it is shown that this limit is achieved here for the linear plane-wave test.

Binary black hole simulations starting from quasi-circular (i.e., zero radial velocity) initial data have orbits with small but nonzero orbital eccentricities. In this paper, the quasi-equilibrium initial-data method is extended to allow nonzero radial velocities to be specified in binary black hole initial data. New low-eccentricity initial data are obtained by adjusting the orbital frequency and radial velocities to minimize the orbital eccentricity, and the resulting (~5 orbit) evolutions are compared with those of quasi-circular initial data. Evolutions of the quasi-circular data clearly show eccentric orbits, with eccentricity that decays over time. The precise decay rate depends on the definition of eccentricity; if defined in terms of variations in the orbital frequency, the decay rate agrees well with the prediction of Peters (1964 Phys. Rev. 136 1224–32). The gravitational waveforms, which contain ~8 cycles in the dominant l = m = 2 mode, are largely unaffected by the eccentricity of the quasi-circular initial data. The overlap between the dominant mode in the quasi-circular evolution and the same mode in the low-eccentricity evolution is about 0.99.

A new representation of the Einstein evolution equations is presented that is first order, linearly degenerate and symmetric hyperbolic. This new system uses the generalized harmonic method to specify the coordinates, and exponentially suppresses all small short-wavelength constraint violations. Physical and constraint-preserving boundary conditions are derived for this system, and numerical tests that demonstrate the effectiveness of the constraint suppression properties and the constraint-preserving boundary conditions are presented.

A method is introduced for solving Einstein's equations using two distinct coordinate systems. The coordinate basis vectors associated with one system are used to project out components of the metric and other fields, in analogy with the way fields are projected onto an orthonormal tetrad basis. These field components are then determined as functions of a second independent coordinate system. The transformation to the second coordinate system can be thought of as a mapping from the original inertial coordinate system to the computational domain. This dual-coordinate method is used to perform stable numerical evolutions of a black-hole spacetime using the generalized harmonic form of Einstein's equations in coordinates that rotate with respect to the inertial frame at infinity; such evolutions are found to be generically unstable using a single rotating-coordinate frame. The dual-coordinate method is also used here to evolve binary black-hole spacetimes for several orbits. The great flexibility of this method allows comoving coordinates to be adjusted with a feedback control system that keeps the excision boundaries of the holes within their respective apparent horizons.

A new representation of the Einstein evolution equations is presented that is first order, linearly degenerate and symmetric hyperbolic. This new system uses the generalized harmonic method to specify the coordinates, and exponentially suppresses all small short-wavelength constraint violations. Physical and constraint-preserving boundary conditions are derived for this system, and numerical tests that demonstrate the effectiveness of the constraint suppression properties and the constraint-preserving boundary conditions are presented.

One of the motivations behind numerical relativity is to provide gravitational wave signals of compact objects to observers using the new gravitational wave detectors. Yet, because of the complexities involved, no dependable signals of binary-black hole coalescences have been established. The work in this proceedings is motivated by how numerical relativity can be used today to predict robust features in gravitational wave signals of binary black-hole coalescence by making approximations to the full problem. To illustrate this, we present results from evolving a Klein-Gordon equation on a frozen background. The background is set by a sequence of initial data in which the binary is in quasi-equilibrium. We probe the data resulting from the evolution for the transition between the linear and non-linear regimes using oscillations of the black holes as our guide. This information is used to motivate a qualitative picture of the gravitational signal of a black-hole coalescence.

A method is introduced for solving Einstein's equations using two distinct coordinate systems. The coordinate basis vectors associated with one system are used to project out components of the metric and other fields, in analogy with the way fields are projected onto an orthonormal tetrad basis. These field components are then determined as functions of a second independent coordinate system. The transformation to the second coordinate system can be thought of as a mapping from the original inertial coordinate system to the computational domain. This dual-coordinate method is used to perform stable numerical evolutions of a black-hole spacetime using the generalized harmonic form of Einstein's equations in coordinates that rotate with respect to the inertial frame at infinity; such evolutions are found to be generically unstable using a single rotating-coordinate frame. The dual-coordinate method is also used here to evolve binary black-hole spacetimes for several orbits. The great flexibility of this method allows comoving coordinates to be adjusted with a feedback control system that keeps the excision boundaries of the holes within their respective apparent horizons.

New boundary conditions are constructed and tested numerically for a general first-order form of the Einstein evolution system. These conditions prevent constraint violations from entering the computational domain through timelike boundaries, allow the simulation of isolated systems by preventing physical gravitational waves from entering the computational domain, and are designed to be compatible with the fixed-gauge evolutions used here. These new boundary conditions are shown to be effective in limiting the growth of constraints in 3D nonlinear numerical evolutions of dynamical black-hole spacetimes.

New boundary conditions are constructed and tested numerically for a general first-order form of the Einstein evolution system. These conditions prevent constraint violations from entering the computational domain through timelike boundaries, allow the simulation of isolated systems by preventing physical gravitational waves from entering the computational domain, and are designed to be compatible with the fixed-gauge evolutions used here. These new boundary conditions are shown to be effective in limiting the growth of constraints in 3D nonlinear numerical evolutions of dynamical black-hole spacetimes.

We have modelled the nonlinear development of the secular bar-mode instability that is driven by gravitational radiation-reaction (GRR) forces in rotating neutron stars. In the absence of any competing viscous effects, an initially uniformly rotating, axisymmetric $n=1/2$ polytropic star with a ratio of rotational to gravitational potential energy $T/|W| = 0.181$ is driven by GRR forces to a bar-like structure, as predicted by linear theory. The pattern frequency of the bar slows to nearly zero, that is, the bar becomes almost stationary as viewed from an inertial frame of reference as GRR removes energy and angular momentum from the star. In this ``Dedekind-like'' state, rotational energy is stored as motion of the fluid in highly noncircular orbits inside the bar. However, in less than 10 dynamical times after its formation, the bar loses its initially coherent structure as the ordered flow inside the bar is disrupted by what appears to be a purely hydrodynamical, short-wavelength, ``shearing'' type instability. The gravitational waveforms generated by such an event are determined, and an estimate of the detectability of these waves is presented., Comment: 25 pages, 9 figures, accepted for publication in ApJ, refereed version, updated, for quicktime movie, see http://www.phys.lsu.edu/~ou/movie/fmode/new/fmode.b181.om4.2e5.mov

Techniques are developed for projecting the solutions of symmetric-hyperbolic evolution systems onto the constraint submanifold (the constraint-satisfying subset of the dynamical field space). These optimal projections map a field configuration to the nearest configuration in the constraint submanifold, where distances between configurations are measured with the natural metric on the space of dynamical fields. The construction and use of these projections are illustrated for a new representation of the scalar field equation that exhibits both bulk and boundary generated constraint violations. Numerical simulations on a black hole background show that bulk constraint violations cannot be controlled by constraint-preserving boundary conditions alone, but are effectively controlled by constraint projection. Simulations also show that constraint violations entering through boundaries cannot be controlled by constraint projection alone, but are controlled by constraint-preserving boundary conditions. Numerical solutions to the pathological scalar field system are shown to converge to solutions of a standard representation of the scalar field equation when constraint projection and constraint-preserving boundary conditions are used together.

Motivated by the need to control the exponential growth of constraint violations in numerical solutions of the Einstein evolution equations, two methods are studied here for controlling this growth in general hyperbolic evolution systems. The first method adjusts the evolution equations dynamically, by adding multiples of the constraints, in a way designed to minimize this growth. The second method imposes special constraint preserving boundary conditions on the incoming components of the dynamical fields. The efficacy of these methods is tested by using them to control the growth of constraints in fully dynamical 3D numerical solutions of a particular representation of the Maxwell equations that is subject to constraint violations. The constraint preserving boundary conditions are found to be much more effective than active constraint control in the case of this Maxwell system.

Techniques are developed for projecting the solutions of symmetric-hyperbolic evolution systems onto the constraint submanifold (the constraint-satisfying subset of the dynamical field space). These optimal projections map a field configuration to the nearest configuration in the constraint submanifold, where distances between configurations are measured with the natural metric on the space of dynamical fields. The construction and use of these projections are illustrated for a new representation of the scalar field equation that exhibits both bulk and boundary generated constraint violations. Numerical simulations on a black hole background show that bulk constraint violations cannot be controlled by constraint-preserving boundary conditions alone, but are effectively controlled by constraint projection. Simulations also show that constraint violations entering through boundaries cannot be controlled by constraint projection alone, but are controlled by constraint-preserving boundary conditions. Numerical solutions to the pathological scalar field system are shown to converge to solutions of a standard representation of the scalar field equation when constraint projection and constraint-preserving boundary conditions are used together.

Motivated by the need to control the exponential growth of constraint violations in numerical solutions of the Einstein evolution equations, two methods are studied here for controlling this growth in general hyperbolic evolution systems. The first method adjusts the evolution equations dynamically, by adding multiples of the constraints, in a way designed to minimize this growth. The second method imposes special constraint preserving boundary conditions on the incoming components of the dynamical fields. The efficacy of these methods is tested by using them to control the growth of constraints in fully dynamical 3D numerical solutions of a particular representation of the Maxwell equations that is subject to constraint violations. The constraint preserving boundary conditions are found to be much more effective than active constraint control in the case of this Maxwell system.