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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Nonlinear evolution of the gravitational radiation (GR) driven instability in the r-modes of neutron stars is studied by full numerical 3D hydrodynamical simulations. The growth of the r-mode instability is found to be limited by the formation of shocks and breaking waves when the dimensionless amplitude of the mode grows to about three in value. This maximum mode amplitude is shown by numerical tests to be rather insensitive to the strength of the GR driving force. Upper limits on the strengths of possible nonlinear mode--mode coupling are inferred. Previously unpublished details of the numerical techniques used are presented, and the results of numerous calibration runs are discussed.

Neutron stars are expected to contain a significant number of hyperons in addition to protons and neutrons in the highest density portions of their cores. Following the work of Jones, we calculate the coefficient of bulk viscosity due to nonleptonic weak interactions involving hyperons in neutron-star cores, including new relativistic and superfluid effects. We evaluate the influence of this new bulk viscosity on the gravitational radiation driven instability in the r-modes. We find that the instability is completely suppressed in stars with cores cooler than a few times 10^9 K, but that stars rotating more rapidly than 10-30% of maximum are unstable for temperatures around 10^10 K. Since neutron-star cores are expected to cool to a few times 10^9 K within seconds (much shorter than the r-mode instability growth time) due to direct Urca processes, we conclude that the gravitational radiation instability will be suppressed in young neutron stars before it can significantly change the angular momentum of the star.

The instability in the r-modes of rotating neutron stars can (in principle) emit substantial amounts of gravitational radiation (GR) which might be detectable by LIGO and similar detectors. Estimates are given here of the detectability of this GR based the non-linear simulations of the r-mode instability by Lindblom, Tohline and Vallisneri. The burst of GR produced by the instability in the rapidly rotating 1.4 solar mass neutron star in this simulation is fairly monochromatic with frequency near 960 Hz and duration about 100 s. A simple analytical expression is derived here for the optimal S/N for detecting the GR from this type of source. For an object located at a distance of 20 Mpc we estimate the optimal S/N to be in the range 1.2 to about 12.0 depending on the LIGO II configuration.

Equilibrium models of differentially rotating nascent neutron stars are constructed, which represent the result of the accretion induced collapse of rapidly rotating white dwarfs. The models are built in a two-step procedure: (1) a rapidly rotating pre-collapse white dwarf model is constructed; (2) a stationary axisymmetric neutron star having the same total mass and angular momentum distribution as the white dwarf is constructed. The resulting collapsed objects consist of a high density central core of size roughly 20 km, surrounded by a massive accretion torus extending over 1000 km from the rotation axis. The ratio of the rotational kinetic energy to the gravitational potential energy of these neutron stars ranges from 0.13 to 0.26, suggesting that some of these objects may have a non-axisymmetric dynamical instability that could emit a significant amount of gravitational radiation.

The evolution of a neutron-star r-mode driven unstable by gravitational radiation (GR) is studied here using numerical solutions of the full non-linear fluid equations. The amplitude of the mode grows to order unity before strong shocks develop which quickly damp the mode. In this simulation the star loses about 40% of its initial angular momentum and 50% of its rotational kinetic energy before the mode is damped. The non-linear evolution causes the fluid to develop strong differential rotation which is concentrated near the surface and especially near the poles of the star.

Gravitational radiation (GR) drives an instability in certain modes of rotating stars. This instability is strong enough in the case of the r-modes to cause their amplitudes to grow on a timescale of tens of seconds in rapidly rotating neutron stars. GR emitted by these modes removes angular momentum from the star at a rate which would spin it down to a relatively small angular velocity within about one year, if the dimensionless amplitude of the mode grows to order unity. A pedagogical level discussion is given here on the mechanism of GR instability in rotating stars, on the relevant properties of the r-modes, and on our present understanding of the dissipation mechanisms that tend to suppress this instability in neutron stars. The astrophysical implications of this GR driven instability are discussed for young neutron stars, and for older systems such as low mass x-ray binaries. Recent work on the non-linear evolution of the r-modes is also presented.

The presence of a viscous boundary layer under the solid crust of a neutron star dramatically increases the viscous damping rate of the fluid r-modes. We improve previous estimates of this damping rate by including the effect of the Coriolis force on the boundary-layer eigenfunction and by using more realistic neutron-star models. If the crust is assumed to be perfectly rigid, the gravitational radiation driven instability in the r-modes is completely suppressed in neutron stars colder than about 1.5 x 10^8 K. Energy generation in the boundary layer will heat the star, and will even melt the crust if the amplitude of the r-mode is large enough. We solve the heat equation explicitly (including the effects of thermal conduction and neutrino emission) and find that the r-mode amplitude needed to melt the crust is approximately a_c = 5 x 10^{-3} for maximally rotating neutron stars. If the r-mode saturates at an amplitude larger than a_c, the heat generated is sufficient to maintain the outer layers of the star in a mixed fluid-solid state analogous to the pack ice on the fringes of the Arctic Ocean. We argue that in young, rapidly rotating neutron stars this effect considerably delays the formation of the crust. By considering the dissipation in the ice flow, we show that the final spin frequency of stars with r-mode amplitude of order unity is close to the value estimated for fluid stars without a crust.

The analogs of r-modes in superfluid neutron stars are studied here. These modes, which are governed primarily by the Coriolis force, are identical to their ordinary-fluid counterparts at the lowest order in the small angular-velocity expansion used here. The equations that determine the next order terms are derived and solved numerically for fairly realistic superfluid neutron-star models. The damping of these modes by superfluid ``mutual friction'' (which vanishes at the lowest order in this expansion) is found to have a characteristic time-scale of about 10^4 s for the m=2 r-mode in a ``typical'' superfluid neutron-star model. This time-scale is far too long to allow mutual friction to suppress the recently discovered gravitational radiation driven instability in the r-modes. However, the strength of the mutual friction damping depends very sensitively on the details of the neutron-star core superfluid. A small fraction of the presently acceptable range of superfluid models have characteristic mutual friction damping times that are short enough (i.e. shorter than about 5 s) to suppress the gravitational radiation driven instability completely.

Stability of the r-modes in rapidly rotating white dwarf stars is investigated. Improved estimates of the growth times of the gravitational-radiation driven instability in the r-modes of the observed DQ Her objects are found to be longer (probably considerably longer) than 6x10^9y. This rules out the possibility that the r-modes in these objects are emitting gravitational radiation at levels that could be detectable by LISA. More generally it is shown that the r-mode instability can only be excited in a very small subset of very hot (T>10^6K), rather massive (M>0.9M_sun) and very rapidly rotating (P_min<P<1.2P_min) white dwarf stars. Further, the growth times of this instability are so long that these conditions must persist for a very long time (t>10^9y) to allow the amplitude to grow to a dynamically significant level. This makes it extremely unlikely that the r-mode instability plays a significant role in any real white dwarf stars.

Techniques are developed here for evaluating the r-modes of rotating neutron stars through second order in the angular velocity of the star. Second-order corrections to the frequencies and eigenfunctions for these modes are evaluated for neutron star models. The second-order eigenfunctions for these modes are determined here by solving an unusual inhomogeneous hyperbolic boundary-value problem. The numerical techniques developed to solve this unusual problem are somewhat non-standard and may well be of interest beyond the particular application here. The bulk-viscosity coupling to the r-modes, which appears first at second order, is evaluated. The bulk-viscosity timescales are found here to be longer than previous estimates for normal neutron stars, but shorter than previous estimates for strange stars. These new timescales do not substantially affect the current picture of the gravitational radiation driven instability of the r-modes either for neutron stars or for strange stars.

Analytical solutions are presented for a class of generalized r-modes of rigidly rotating uniform density stars---the Maclaurin spheroids---with arbitrary values of the angular velocity. Our analysis is based on the work of Bryan; however, we derive the solutions using slightly different coordinates that give purely real representations of the r-modes. The class of generalized r-modes is much larger than the previously studied `classical' r-modes. In particular, for each l and m we find l-m (or l-1 for the m=0 case) distinct r-modes. Many of these previously unstudied r-modes (about 30% of those examined) are subject to a secular instability driven by gravitational radiation. The eigenfunctions of the `classical' r-modes, the l=m+1 case here, are found to have particularly simple analytical representations. These r-modes provide an interesting mathematical example of solutions to a hyperbolic eigenvalue problem.

A new method is presented here for evaluating approximately the pulsation modes of relativistic stellar models. This approximation relies on the fact that gravitational radiation influences these modes only on timescales that are much longer than the basic hydrodynamic timescale of the system. This makes it possible to impose the boundary conditions on the gravitational potentials at the surface of the star rather than in the asymptotic wave zone of the gravitational field. This approximation is tested here by predicting the frequencies of the outgoing non-radial hydrodynamic modes of non-rotating stars. The real parts of the frequencies are determined with an accuracy that is better than our knowledge of the exact frequencies (about 0.01%) except in the most relativistic models where it decreases to about 0.1%. The imaginary parts of the frequencies are determined with an accuracy of approximately M/R, where M is the mass and R is the radius of the star in question.

The magnitudes of the external gravitational perturbations associated with the normal modes of the Sun are evaluated to determine whether these solar oscillations could be observed with the proposed Laser Interferometer Space Antenna (LISA), a network of satellites designed to detect gravitational radiation. The modes of relevance to LISA---the $l=2$, low-order $p$, $f$ and $g$-modes---have not been conclusively observed to date. We find that the energy in these modes must be greater than about $10^{30} \rm{ergs}$ in order to be observable above the LISA detector noise. These mode energies are larger than generally expected, but are much smaller than the current observational upper limits. LISA may be confusion-limited at the relevant frequencies due to the galactic background from short-period white dwarf binaries. Present estimates of the number of these binaries would require the solar modes to have energies above about $10^{33} \rm{ergs}$ to be observable by LISA.