Search

We compute encounters of a realistic white dwarf model with a massive black hole in the regime where relativistic effects are important, using a three-dimensional, finite-difference, Eulerian, piecewise parabolic method (PPM) hydrodynamical code. Both disruptive and nondisruptive encounters are considered. We identify and discuss relativistic effects important for the problem: relativistic shift of the pericenter distance, time delay, relativistic precession, and the tensorial structure of the tidal forces. In the nondisruptive case, stripping of matter takes place. In the surface layers of the surviving core, complicated hydrodynamical phenomena are revealed. In both disruptive and nondispruptive encounters, material flows out in the form of two thin, S-shaped, supersonic jets. Our results provide realistic initial conditions for the subsequent investigation of the dynamics of the debris in the field of the black hole. We evaluate the critical conditions for complete disruption of the white dwarf, and compare our results with the corresponding results for nonrelativistic encounters.

An approximation for the moment of inertia of a neutron star in terms of only its mass and radius is presented, and insight into it is obtained by examining the behavior of the relativistic structural equations. The approximation is accurate to approximately 10% for a variety of nuclear equations of state, for all except very low mass stars. It is combined with information about the neutron-star crust to obtain a simple expression (again in terms only of mass and radius) for the fractional moment of inertia of the crust.

We show that, for physically reasonable bulk and surface properties, the lowest energy state of dense matter consists of quark matter coexisting with nuclear matter in the presence of an essentially uniform background of electrons. We estimate the size and nature of spatial structure in this phase, and show that at the lowest densities the quark matter forms droplets embedded in nuclear matter, whereas at higher densities it can exhibit a variety of different topologies. A finite fraction of the interior of neutron stars could consist of matter in this new phase, which would provide new mechanisms for glitches and cooling.

The quasiparticle decay rates for quarks and gluons in quark-gluon plasmas are calculated by solving the kinetic equation. Introducing an infrared cutoff to allow for nonperturbative effects, we evaluate the quasiparticle lifetime at momenta greater than the inverse Debye screening length to leading order in the coupling constant.

We calculate properties of neutron star matter at subnuclear densities using an improved nuclear Hamiltonian. Nuclei disappear and the matter becomes uniform at a density of about 0.6n(s), where n(s) of about 0.16/cu fm is the saturation density of nuclear matter. As a consequence, the mass of matter in the crusts of neutron stars is only about half as large as previously estimated. In about half of that crustal mass, nuclear matter occurs in shapes very different from the roughly spherical nuclei familiar at lower densities. The thinner crust and the unusual nuclear shape have important consequences for theories of the rotational and thermal evolution of neutron stars, especialy theories of glitches.

It is at present impossible to predict the interior constitution of neutron stars based on theory and results from laboratory studies. It has been proposed that it is possible to obtain information on neutron star interiors by studying thermal radiation from their surfaces, because neutrino emission rates, and hence the temperature of the central part of a neutron star, depend on the properties of dense matter. The theory predicts that neutron stars cool relatively slowly if their cores are made up of nucleons, and cool faster if the matter is in an exotic state, such as a pion condensate, a kaon condensate, or quark matter. This view has recently been questioned by the discovery of a number of other processes that could lead to copious neutrino emission and rapid cooling.

Direct Urca processes with hyperons and/or nucleon isobars can occur in dense matter as long as the concentration of Lambda hyperons exceeds a critical value that is less than 3 percent and is typically about 0.1 percent. The neutrino luminosities from the hyperon Urca processes are about 5-100 times less than the typical luminosity from the nucleon direct Urca process, if the latter process is not forbidden, but they are larger than those expected from other sources. These direct Urca processes provide avenues for rapid cooling of neutron stars which invoke neither exotic states nor the large proton fraction (of order 0.11-0.15) required for the nucleon direct Urca process.

The processes leading to tidal capture of stars by a massive black hole and the consequences of these processes in a dense stellar cluster are discussed in detail. When the amplitude of a tide and the subsequent oscillations are sufficiently large, the energy deposited in a star after periastron passage and formation of a bound orbit cannot be estimated directly using the linear theory of oscillations of a spherical star, but rather numerical estimates must be used. The evolution of a star after tidal capture is discussed. The maximum ratio R of the cross-section for tidal capture to that for tidal disruption is about 3 for real systems. For the case of a stellar system with an empty capture loss cone, even in the case when the impact parameter for tidal capture only slightly exceeds the impact parameter for direct tidal disruption, tidal capture would be much more important than tidal disruption.

It is shown that the direct URCA process can occur in neutron stars if the proton concentration exceeds some critical value in the range 11-15 percent. The proton concentration, which is determined by the poorly known symmetry energy of matter above nuclear density, exceeds the critical value in many current calculations. If it occurs, the direct URCA process enhances neutrino emission and neutron star cooling rates by a large factor compared to any process considered previously.

The equation of state of cold dense matter in strong magnetic fields is calculated in the Thomas-Fermi and Thomas-Fermi-Dirac approximations. For use in the latter calculation, a new expression is derived for the exchange energy of the uniform electron gas in a strong magnetic field. Detailed calculations of the density profile in the surface region of a neutron star are described for a variety of equations of state, and these show that the surface density profile is strongly affected by the magnetic field, irrespective of whether or not matter in a magnetic field has a condensed state bound with respect to isolated atoms. It is also shown that, as a consequence of the field dependence of the screening potential, magnetic fields can significantly increase nuclear reaction rates.

The envelopes of nonmagnetic neutron stars are studied using the best available opacities and equation of state. The general relativistic equations of the structure and evolution of nonmagnetic neutron stars are discussed, and it is shown that they can be reduced to a single equation for calculating the thermal structure of neutron star envelopes. The physical input needed to solve the thermal structure equation is investigated and the numerical results of envelope model calculations are presented. It is shown that the thermal structure of neutron star envelopes is determined by the single parameter T(s) to the 4th/g(s), where T(s) is the effective surface temperature and g(s) the surface gravity of the star. This result is used to derive a number of other scaling relations, and the effects of general relativity on the envelope thermal structure are examined in detail. The results of a sensitivity analysis of the regional opacity needed to obtain a reliable relationship between the temperatures of the inner and outer boundaries of the envelope is presented.

A method is given for determining the sensitivity of certain types of calculations to the uncertainties in the input physics or model parameters; this method is applicable to problems that involve solutions to coupled, ordinary differential equations. In particular the sensitivity of calculations of the thermal structure of neutron star envelopes to uncertainties in the opacity and equation of state is examined. It is found that the uncertainties in the relationship between the surface and interior temperatures of a neutron star are due almost entirely to the imprecision in the values of the conductive opacity in the region where the ions form a liquid; here the conductive opacity is, for the most part, due to the scattering of electrons from ions.

It is pointed out that the separate study of neutron star envelopes is a valuable technique for gaining insights into neutron star cooling calculations. It is shown that, to an excellent approximation, the temperature at the inner boundary of the envelope, T sub b, is a function only of T sub s to the 4th power over g sub s, where T sub s is the effective temperature and g sub s is the surface gravity. This result makes it possible to investigate very simply the effects of general relativity and of opacity changes on the T sub b - T sub s relation and to compare calculations by different workers, even though they use different neutron star models. It is also shown that the most important opacity source for obtaining a reliable T sub b - T sub s relation is the conductive opacity of matter when the ions are liquid. Calculations are presented in which particular attention is given to this region and the best available opacities are used.