Fast and accurate computational tools for gravitational waveforms from binary stars with any orbital eccentricity

The relevance of orbital eccentricity in the detection of gravitational radiation from (steady state) binary stars is emphasized. Computationally effective (fast and accurate) tools for constructing gravitational wave templates from binary stars with any orbital eccentricity are introduced including tight estimation criteria of the pertinent truncation and approximation errors.     

On the interaction of jet plasma with dark matter subhaloes in active galaxies

The interaction of fragmented plasma of active galactic nuclei jets with galactic haloes via gravitational scattering and lensing by dark matter subhaloes is studied using analytical calculations and numerical Monte-Carlo method. The lensing of jet radiation by halo masses is found to be negligible and unobservable. Moving through a galactic halo jet plasma fragments are sequentially deflected on hyperbolic orbits by gravitational field of subhaloes and deviates at some angles when leaving halo, causing widening of the jet. Based on this model jet opening angles are calculated numerically for various values of jet and halo characteristics. Though these angles are very small, gravitational scattering by halo masses results in specific radial profile of jet radiation intensity, that does not depend on halo mass distribution and jet properties. The intensity of jet radiation, obeying the derived profile, decreases by reasonable observable factors giving possibility to probe the presence of dark matter subhaloes.     

Laser Interferometer Space Antenna double black holes: dynamics in gaseous nuclear discs

We study the inspiral of double black holes, with masses in the Laser Interferometer Space Antenna ( LISA ) window of detectability, orbiting inside a massive circumnuclear, rotationally supported gaseous disc. Using high-resolution smoothed particle hydrodynamics simulations, we follow the black hole dynamics in the early phase when gas-dynamical friction acts on the black holes individually, and continue our simulation until they form a close binary. We find that in the early sinking the black holes lose memory of their initial orbital eccentricity if they corotate with the gaseous disc. As a consequence, the massive black holes bind forming a binary with a low eccentricity, consistent with zero within our numerical resolution limit. The cause of circularization resides in the rotation present in the gaseous background where dynamical friction operates. Circularization may hinder gravitational waves from taking over and leading the binary to coalescence. In the case of counter-rotating orbits, the initial eccentricity (if present) does not decrease, and the black holes may bind forming an eccentric binary. When dynamical friction has subsided, for equal mass black holes and regardless their initial eccentricity, angular momentum loss, driven by the gravitational torque exerted on the binary by surrounding gas, is nevertheless observable down to the smallest scale probed (≃1 pc). In the case of unequal masses, dynamical friction remains efficient down to our resolution limit, and there is no sign of formation of any ellipsoidal gas distribution that may further harden the binary. During inspiral, gravitational capture of gas by the black holes occurs mainly along circular orbits; eccentric orbits imply high relative velocities and weak gravitational focusing. Thus, the active galactic nucleus activity may be excited during the black hole pairing process and double active nuclei may form when circularization is completed, on distance scales of tens of parsecs.     

The possible role of r-modes in post-glitch relaxation of the Crab pulsar

The loss of angular momentum through gravitational radiation, driven by the excitation of r-modes, is considered for neutron stars that have rotation frequencies lower than the associated critical frequency. We find that for reasonable values of the initial amplitudes of such pulsation modes of the star, being excited at the event of a glitch in a pulsar, the total post-glitch losses correspond to a negligible fraction of the initial rise of the spin frequency in the case of Vela and older pulsars. However, for the Crab pulsar the same effect would result, within a few months, in a decrease in its spin frequency by an amount larger than its glitch-induced frequency increase. This could provide an explanation for the peculiar behaviour observed in the post-glitch relaxations of the Crab pulsar.     

Simulating a stochastic background of gravitational waves from neutron star formation at cosmological distances

We develop a temporal simulation of the potentially detectable gravitational wave background from neutron star formation at cosmological distances. By using a recent model for the evolving star formation rate, we investigate the statistical distribution of gravitational wave amplitudes due to supernovae that result in neutron star formation in the Einstein–de Sitter cosmology. We find that the gravitational wave amplitude distribution in our frame is highly skewed, with skewness related to the distribution of sources, and that the potentially detectable gravitational wave strain is dominated by sources at a redshift of     . Time traces of the simulation, using selected waveforms, are presented graphically and are also made available as web-based audio files. The method developed can readily be extended to different cosmologies, as well as to incorporate other waveforms and source types. This type of simulation will be useful in testing and optimizing detection strategies for gravitational wave backgrounds due to various types of individually undetectable astrophysical sources.     

Newtonian hydrodynamics of the coalescence of black holes with neutron stars – IV. Irrotational binaries with a soft equation of state

We present the results of three-dimensional hydrodynamical simulations of the final stages of in-spiral in a black hole–neutron star binary, when the separation is comparable to the stellar radius. We use a Newtonian smooth particle hydrodynamics (SPH) code to model the evolution of the system, and take the neutron star to be a polytrope with a soft (adiabatic indices     and     equation of state and the black hole to be a Newtonian point mass. The only non-Newtonian effect we include is a gravitational radiation back reaction force, computed in the quadrupole approximation for point masses. We use irrotational binaries as initial conditions for our dynamical simulations, which are begun when the system is on the verge of initiating mass transfer and followed for approximately 23 ms. For all the cases studied we find that the star is disrupted on a dynamical time-scale, and forms a massive     accretion torus around the spinning (Kerr) black hole. The rotation axis is clear of baryons (less than 10⁻⁵ M_⊙ within 10°) to an extent that would not preclude the formation of a relativistic fireball capable of powering a cosmological gamma-ray burst. Some mass (the specific amount is sensitive to the stiffness of the equation of state) may be dynamically ejected from the system during the coalescence and could undergo r-process nucleosynthesis. We calculate the waveforms, luminosities and energy spectra of the gravitational radiation signal, and show how they reflect the global outcome of the coalescence process.     

Dynamical evolution of Saturn's F ring dust particles

Saturn's F ring has been the subject of study due to its peculiar structure and the proximity to two satellites, named Prometheus (interior) and Pandora (exterior to the ring), which cause perturbations to the ring particles. Early results from Voyager data have proposed that the ring is populated with centimetre- and micrometre-sized particles. The Cassini spacecraft also detected a less dense part in the ring with width of 700 km. Small particles suffer the effects of solar radiation. Burns et al. showed that due to effects of one component of the solar radiation, the Poynting–Robertson drag, a ring particle will decay in the direction of the planet in a time much shorter than the age of the Solar system. In this work, we have analysed a sample of dust particles (1, 3, 5 and 10 μm) under the effects of solar radiation, the Poynting–Robertson drag and the radiation pressure components and the gravitational effects of the satellites Prometheus and Pandora. In this case, the high increase of the eccentricity of the particles leads almost all of them to collide with the outer edge of the A ring. The inclusion of the oblateness of Saturn in this system significantly changes the outcome, since the large variation of the eccentricity is reduced by the oblateness effect. As a result, there is an increase in the lifetime of the particle in the envelope region. Our results show that even the small dust particles, which are very sensitive to the effects of solar radiation, have an orbital evolution similar to larger particles located in the F ring. The fate of all particles is a collision with Prometheus or Pandora in less than 30 years. On the other hand, collisions of these particles with moonlets/clumps present in the F ring could change this scenario.     

Small-Amplitude Disturbances in a Radiating and Scattering Grey MediumIII. Gravitational Effects on the Solutions of Given Real Wave Number k

We study the fundamental modes of radiation hydrodynamic linear waves that arise from one-dimensional small-amplitude initial fluctuations with wave number k in a radiating and scattering grey medium by taking into account the gravitational effects. The equation of radiative acoustics is derived from three hydrodynamic equations, Poisson’s equation, and two moment equations of radiation, by assuming a spherical symmetry for the matter and radiation and by using the Eddington approximation. We solve the dispersion relation as a quintic function of angular frequency ω, the wave number k being a real parameter. Numerical results reveal that wave patterns of five solutions are distinguished into three types: the radiation-dominated, type 1, and type 2 matter-dominated cases. In the case of no gravitaional effects (Kaneko et al., 2005), the following wave modes appear: radiation wave, conservative radiation wave, entropy wave, Newtonian-cooling wave, opacity-damped and cooling-damped waves, constant-volume and constant-pressure diffusions, adiabatic sound wave, cooling-damped and drag-force-damped isothermal sound waves, isentropic radiation-acoustic wave, and gap mode. Meanwhile, the gravitaional effects being taken into account, the growing gravo-diffusion mode newly arises from the constant-pressure diffusion at the point that k agrees with Jeans’ wave number specified by the isothermal sound speed. This mode changes to the growing radiation-acoustic gravity mode near the point that k becomes Jeans’ wave number specified by the isentropic radiation-acoustic speed. In step with a transition between them, the isentropic radiation-acoustic wave splits into the damping radiation-acoustic gravity mode and constant-volume diffusion. The constant-volume diffusion emerges twice if the gravitational effects are taken into account. Since analytic solutions are derived for all wave modes, we discuss their physical significance. The critical conditions are given which distinguish between radiation-dominated and type 1 matter-dominated cases, and between type 1 and type 2 matter-dominated cases. Waves in a self-gravitating scattering grey medium are also analyzed, which provides us some hints for the effects of energy and momentum exchange between matter and radiation.     

Supernovae as a Powerful Source of Gravitational Radiation

We consider gravitational radiation in the nonspherical symmetrical evolution of supernovae. Unlike the gravitational radiation analysis of Imshennik & Popov performed using the formalism of Peters & Mathews, our analysis uses the (PN)5/2 approximation by Damour & Deruelle and Lincoln & Will. We show that the eccentricity is larger than 0.1 at the time when the Roche lobe is filled by the low-mass component, demonstrating that the conclusion of Imshennik and Popov that the final eccentricity is less than 0.1 is incorrect. If the supernova lies in the Large Magellanic Cloud (R = 50 kpc), we estimate the amplitude of gravitational waves to be h 8 × 10-20. The frequency of the emitted gravitational waves is about 1 kHz.     

Orbital Decay of a Binary System Consisting of Spin-Down Neutron Stars

{W}e consider the gravitational radiation from two time variable mass stars, orbiting around each other under the influence of gravity. The total rates of the variation of the energy, angular momentum, semimajor axis, eccentricity and orbital period are obtained. The results could be important for the understanding of general relativistic effects in the case of the variation of the gravitational mass due to spinning down of the compact stars, which sensitively depends on the equations of state. The cases of the binary systems PSR 1913+16 and PSR 1534+12 are analyzed in detail, and, for different equations of state of nuclear matter, the corrections to the orbital decay due to gravitational radiation and to the spinning down of the pulsars are calculated. The results show that a future significant improvement in the observational techniques could lead to the observation of the specific general relativistic effect of mass variation of pulsars due to spinning down, via the study of orbital decay, even in slowly rotating binary systems.     

A Research on the Gravitational Wave Radiation of OJ 287

The research on quasar OJ 287 has lasted over 100 years. OJ 287 exhibits the phenomenon of periodic two-peak outbursts with the eruptive period of 12 years. Observations are rather well interpreted with the black hole binary model. In this model, the secondary black hole moves around the primary black hole and crashes against the accretion disk of the primary black hole, causing outbursts. This model reasonably explains the light curves of OJ 287 and correctly predicts the time of future outbursts. These indirectly justify the precessional effect of general relativity and the existence of gravitational waves. The massive black hole in the center of galaxy is an important kind of gravitational wave source. However, the number of the galaxies with precisely determined kinematical equations of inner components is quite small. The precise kinematic orbits of black holes are provided by the black hole binary model, so the radiation of gravitational waves can be studied on the basis of these kinematic orbits. Based on the existing work, the evolutionary relations of the radiation power and waveform of gravitational waves with time are first derived by using the post-Newtonian approximation method. According to the current progress of the detection equipment of gravitational waves, i.e., IPTA (International Pulsar Timing Array), the direct detection of gravitational waves from OJ 287 may be possible within the future decade.     

Gravitational Radiation of Relativistic <Emphasis Type="Italic">n</Emphasis> = 1 Polytropes

The gravitational radiation of n = 1 polytropes undergoing quasiradial pulsations is examined. The intensity of the gravitational radiation and the gravitational wave amplitudes are calculated for polytropic models of white dwarfs and neutron stars when the energy of rotation of the object serves as the source of the radiated energy. Calculations of h₀ show that objects with a polytropic equation of state can describe the expected gravitational radiation from white dwarfs and neutron stars. The gravitational radiation of polytropic models of galactic nuclei and quasars is also examined. These objects can create a high enough background of gravitational radiation at frequencies of 10^-8–10^-11 Hz for gravitational wave detectors operating in this frequency range. __________ Translated from Astrofizika, Vol. 48, No. 4, pp. 603–612 (November 2005).     

Structure and evolution of low-mass W Ursae Majoris type systems – III. The effects of the spins of the stars

In a previous paper, using Eggleton's stellar evolution code, we have discussed the structure and evolution of low-mass W Ursae Majoris (W UMa) type contact binaries with angular momentum loss owing to gravitational radiation or magnetic braking. We find that gravitational radiation is almost insignificant for cyclic evolution of low-mass W UMa type systems, and it is possible for angular momentum to be lost from W UMa systems in a magnetic stellar wind. The weaker magnetic activity shown by observations in W UMa systems is likely caused by the lower mass of the convective envelopes in these systems than in similar but non-contact binaries. The spin angular momentum cannot be neglected at any time for W UMa type systems, especially for those with extreme mass ratios. The spin angular momenta of both components are included in this paper and they are found to have a significant influence on the cyclic evolution of W UMa systems. We investigate the influence of the energy transfer on the common convective envelopes of both components in detail. We find that the mass of the convective envelope of the primary in contact evolution is slightly more than that in poor thermal contact evolution, and that the mass of the convective envelope of the secondary in contact evolution is much less than that in poor thermal contact evolution. Meanwhile, the rate of angular momentum loss of W UMa type systems is much lower than that of poor thermal contact systems. This is indeed caused by the lower masses of the convective envelopes of the components in W UMa type systems. Although the models with angular momentum loss for W UMa systems exhibit cyclic evolution, they seem to show that a W UMa system cannot continue this type of cyclic evolution indefinitely, and it might coalesce into a fast-rotating star after about 1200 cycles of evolution (about  7.0 × 10⁹ yr  ).     

The onset of General Relativity: gravitationally redshifted emission lines

We study and quantify gravitational redshift by means of relativistic ray tracing simulations of emission lines. The emitter model is based on thin, Keplerian rotating rings in the equatorial plane of a rotating black hole. Emission lines are characterised by a generalized fully relativistic Doppler factor or redshift associated with the line core. Two modes of gravitational redshift, shift and distortion, become stronger with the emitting region closer to the Kerr black hole. Shifts of the line cores reveal an effect at levels of 0.0015 to 60% at gravitational radii ranging from 10⁵ to 2. The corresponding Doppler factors range from 0.999985 to 0.4048. Line shape distortion by strong gravity, i.e. very skewed and asymmetric lines occur at radii smaller than roughly ten gravitational radii. Gravitational redshift decreases with distance to the black hole but remains finite due to the asymptotical flatness of black hole space–time. The onset of gravitational redshift can be tested observationally with sufficient spectral resolution. Assuming a resolving power of ∼100000, yielding a resolution of ≈0.1 Å for optical and near‐infrared broad emission lines, the gravitational redshift can be probed out to approximately 75000 gravitational radii. In general, gravitational redshift is an indicator of black hole mass and spin as well as for the inclination angle of the emitter, e.g. an accretion disk. We suggest to do multi‐wavelength observations because all redshifted features should point towards the same central mass. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

银河系球状星团轨道参数分布的Monte Carlo模拟

选取前文中所列出的29个累积光谱型为F型的球状星团中的3个作为样本，深入研究了初始观测资料的不确定性和选用不同的银河系引力势模型，对样本星团轨道参数的影响。首先采用Monte Carlo方法产生3个样本球状星团的模拟初始观测数据，而后，以这些模拟数据为初始条件，在3种不同的银河系引力势模型下进行轨道计算，得到此3个样本的模拟轨道参数。模拟计算的结果表明：根据模拟初始数据生成的样本轨道参数分布形态大致可分为高斯分布、准高斯分布和非高斯分布等3类；初始观测数据的不确定性对样本轨道参数分布的影响，与样本星团的选择和轨道参数的类型有关；选用不同的银河系引力势模型，对3个样本星团的各个轨道参数的分布和形态结构也会产生不同程度的影响。该工作的结果，可供深入研究球状星团的整体运动和动力学性质等问题参考。     

On the rotation of co-orbital bodies in eccentric orbits

We investigate the resonant rotation of co-orbital bodies in eccentric and planar orbits. We develop a simple analytical model to study the impact of the eccentricity and orbital perturbations on the spin dynamics. This model is relevant in the entire domain of horseshoe and tadpole orbit, for moderate eccentricities. We show that there are three different families of spin–orbit resonances, one depending on the eccentricity, one depending on the orbital libration frequency, and another depending on the pericenter’s dynamics. We can estimate the width and the location of the different resonant islands in the phase space, predicting which are the more likely to capture the spin of the rotating body. In some regions of the phase space the resonant islands may overlap, giving rise to chaotic rotation.     

Interplanetary dust dynamics: I. Long-term gravitational effects of the inner planets on zodiacal dust

Whereas the inner planets' perturbations on meteoroids' and larger interplanetary bodies' orbits have been studied extensively, they are usually neglected in studies of the dynamics of smaller particles producing the zodiacal light through scattering of sunlight. Forces acting on these dust particles are fairly well known and include radiation forces and interaction with the solar wind. This article is the first in a series aimed at improving our knowledge of the dynamical evolution of dust in interplanetary space by studying the combined effects of these perturbations including gravitational perturbations by the planets Venus, Earth, Mars, and Jupiter. The necessity of including effects of the inner planets in dust dynamics investigations is established. Sample trajectories are presented to illustrate commonly occurring phenomenae, such as nonmonotonic changes in semimajor axis, eccentricity, inclination, and in the line of nodes. These perturbations are shown to be due to the inner planets as opposted to Jupiter or nongravitational forces.     

Production of massive spin — 1/2 particles in Robertson-Walker universes with external electromagnetic fields

In this paper we investigate the combined influence of both cosmological and electromagnetic particle creation mechanisms upon massive particles with spin 1/2 on the basis of general covariant Dirac theory.Curved space-time, a radiation-dominated Friedmann universe, is treated as an unquantized gravitational field and the low-frequency part of the 2.7 K background radiation is approximated by homogeneous, constant, and parallel external electric and magnetic fields. We calculate the number density of spin 1/2 particles with massm which are created under the influence of both these external fields.We find that the electric field and the magnetic field both amplify the genuine, purely gravitational particle production. This influence of the magnetic field, which is in contrast to its reducing effect as far as the creation of spin-zero particles is concerned, can clearly be traced back to its coupling to the spin of the particles.Under certain conditions the electromagnetic fields in the early universe can influence the particle creation process even more than the gravitational field.     

Chaotic zone boundary for low free eccentricity particles near an eccentric planet

We consider particles with low free or proper eccentricity that are orbiting near planets on eccentric orbits. Through collisionless particle integration, we numerically find the location of the boundary of the chaotic zone in the planet's corotation region. We find that the distance in semimajor axis between the planet and boundary depends on the planet mass to the 2/7 power and is independent of the planet eccentricity, at least for planet eccentricities below 0.3. Our integrations reveal a similarity between the dynamics of particles at zero eccentricity near a planet in a circular orbit and with zero free eccentricity particles near an eccentric planet. The 2/7th law has been previously explained by estimating the semimajor at which the first-order mean motion resonances are large enough to overlap. Orbital dynamics near an eccentric planet could differ due to first-order corotation resonances that have strength proportional to the planet's eccentricity. However, we find that the corotation resonance width at low free eccentricity is small; also the first-order resonance width at zero free eccentricity is the same as that for a zero-eccentricity particle near a planet in a circular orbit. This accounts for insensitivity of the chaotic zone width to planet eccentricity. Particles at zero free eccentricity near an eccentric planet have similar dynamics to those at zero eccentricity near a planet in a circular orbit.     

The equilibrium configuration of rapid rotating compact stars and some gravitational effects

In this paper, the equilibrium configurations of rapid rotating compact stars and some gravitational effects are studied within the general relativity by use of the Harrison-Wheeler equation of state and by the self-consistent field method. Numerical calculations show that the equilibrium configuration of a rotating star is a spheroid. For large spin velocities, say, ω > 3.0 × 10² sec⁻¹ the eccentricity and mass increase very rapidly as the angular velocity increases, for the critical angular velocity of the rotating star, the eccentricity is about 0.7, the increase in mass is about 10–35%. The difference of the gravitational redshifts at the surface of the star caused by rotation, and the difference of the light bending when the beam moves in the direction of rotation or in the opposite direction are obvious.