Cataclysmic variables as multimessenger sources: the gravitational wave emission

The direct detection of gravitational waves (GW) from merging binary black holes and neutron stars mark the beginning of a new era in gravitational physics, and it brings forth new opportunities to test theories of gravity. To this end, it is crucial to search for anomalous deviations from general relativity in a model-independent way, irrespective of gravity theories, GW sources, and background spacetimes. In this paper, we propose a new universal framework for testing gravity with GW, based on the generalized propagation of a GW in an effective field theory that describes modification of gravity at cosmological scales. Then we perform a parameter estimation study, showing how well the future observation of GW can constrain the model parameters in the generalized models of GW propagation.

The direct detection of gravitational waves (GWs) by pulsar timing arrays (PTAs) is very likely within the next decade. While the stochastic GW background is a promising candidate for detection it is also possible that single resolvable sources may be detectable as well. In this work we will focus on the detection and characterization of single GW sources from supermassive black hole binaries. We introduce a fully Bayesian data analysis pipeline that is meant to carry out a search, characterization, and evaluation phase. This will allow us to rapidly locate the global maxima in parameter space, map out the posterior, and finally weigh the evidence of a GW detection through a Bayes factor. Here we will make use of an adaptive metropolis algorithm and parallel tempering. We test this algorithm on realistic simulated data that are representative of modern PTAs.

Current bounds from the polarization of the CMB predict the scale-invariant gravitational wave (GW) background of inflation to be out of reach for upcoming GW interferometers. This prospect dramatically changes if the inflaton is a pseudoscalar, in which case its generic coupling to any abelian gauge field provides a new source of GWs, directly related to the dynamics of inflation. This opens up new ways of probing the scalar potential responsible for cosmic inflation. Dividing inflation models into universality classes, we analyze the possible observational signatures. One of the most promising scenarios is Starobinsky inflation, which may lead to observational signatures both in direct GW detection as well as in upcoming CMB detectors. In this case, the complementarity between the CMB and direct GW detection, as well as the possibility of a multi-frequency analysis with upcoming ground and space based GW interferometers, may provide a first clue to the microphysics of inflation.

We report on PTF1 J191905.19+481506.2, a newly discovered, partially eclipsing, outbursting AM CVn system found in the Palomar Transient Factory synoptic survey. This is only the second known eclipsing AM CVn system. We use high-speed photometric observations and phase-resolved spectroscopy to establish an orbital period of 22.4559(3) min. We also present a long-term light curve and report on the normal and super-outbursts regularly seen in this system, including a super-outburst recurrence time of 36.8(4) d. We use the presence of the eclipse to place upper and lower limits on the inclination of the system and discuss the number of known eclipsing AM CVn systems versus what would be expected.

We extend recent work that included the effect of pressure forces to derive the precession rate of eccentric accretion discs in cataclysmic variables to the case of double degenerate systems. We find that the logical scaling of the pressure force in such systems results in predictions of unrealistically high primary masses. Using the prototype AM CVn as a calibrator for the magnitude of the effect, we find that there is no scaling that applies consistently to all the systems in the class. We discuss the reasons for the lack of a superhump period to mass ratio relationship analogous to that known for SU UMa systems and suggest that this is because these secondaries do not have a single valued mass–radius relationship. We highlight the unreliability of mass ratios derived by applying the SU UMa expression to the AM CVn binaries.

We present radial velocity observations of four extremely low-mass (0.2 Msol) white dwarfs. All four stars show peak-to-peak radial velocity variations of 540 - 710 km/s with 1.0 - 5.9 hr periods. The optical photometry rules out main-sequence companions. In addition, no milli-second pulsar companions are detected in radio observations. Thus the invisible companions are most likely white dwarfs. Two of the systems are the shortest period binary white dwarfs yet discovered. Due to the loss of angular momentum through gravitational radiation, three of the systems will merge within 500 Myr. The remaining system will merge within a Hubble time. The mass functions for three of the systems imply companions more massive than 0.46 Msol; thus those are carbon/oxygen core white dwarfs. The unknown inclination angles prohibit a definitive conclusion about the future of these systems. However, the chance of a supernova Ia event is only 1% to 5%. These systems are likely to form single R Coronae Borealis stars, providing evidence for a white dwarf + white dwarf merger mechanism for these unusual objects. One of the systems, SDSS J105353.89+520031.0 has a 70% chance of having a low-mass white dwarf companion. This system will probably form a single helium-enriched subdwarf O star. All four white dwarf systems have unusal mass ratios of < 0.2-0.8 that may also lead to the formation of AM CVn systems.

The direct detection of gravitational waves (GWs) opened a new chapter in the modern cosmology to probe possible deviations from the general relativity (GR) theory. In the present work, we investigate for the first time the modified GW form propagation from the inspiraling of compact binary systems within the context of $f(T)$ gravity in order to obtain new forecasts/constraints on the free parameter of the theory. First, we show that the modified waveform differs from the GR waveform essentially due to induced corrections on the GWs amplitude. Then, we discuss the forecasts on the $f(T)$ gravity assuming simulated sources of GWs as black hole binaries, neutron star binaries and black hole--neutron star binary systems, which emit GWs in the frequency band of the Advanced LIGO (aLIGO) interferometer and of the third generation Einstein Telescope (ET). We show that GW sources detected within the aLIGO sensitivity can return estimates of the same order of magnitude of the current cosmological observations. On the other hand, detection within the ET sensitivity can improve by up to 2 orders of magnitude the current bound on the $f(T)$ gravity. Therefore, the statistical accuracy that can be achieved by future ground based GW observations, mainly with the ET detector (and planed detectors with a similar sensitivity), can allow strong bounds on the free parameter of the theory, and can be decisive to test the theory of gravitation.

The direct detection of gravitational waves will provide valuable astrophysical information about many celestial objects. The most promising sources of gravitational waves are neutron stars and black holes. These objects emit waves in a very wide spectrum of frequencies determined by their quasi-normal modes oscillations. In this work we are concerned with the information we can extract from f and pI-modes when a candidate leaves its signature in the resonant mass detectors ALLEGRO, EXPLORER, NAUTILUS, MiniGrail and SCHENBERG. Using the empirical equations, that relate the gravitational wave frequency and damping time with the mass and radii of the source, we have calculated the radii of the stars for a given interval of masses M in the range of frequencies that include the bandwidth of all resonant mass detectors. With these values we obtain diagrams of mass-radii for different frequencies that have allowed to determine the better candidates to future detection taking in account the compactness of the source. Finally, to determine which are the models of compact stars that emit gravitational waves in the frequency band of the mass resonant detectors, we compare the mass-radii diagrams obtained by different neutron stars sequences from several relativistic hadronic equations of state (GM1, GM3, TM1, NL3) and quark matter equations of state (NJL, MTI bag model). We verify that quark stars obtained from MIT bag model with bag constant equal to 170 MeV and quark matter in color-superconductivity phase are the best candidates for mass resonant detectors.

Large-area sky surveys show that massive galaxies undergo at least one major merger in a Hubble time. Ongoing pulsar timing array (PTA) experiments are aimed at measuring the gravitational-wave (GW) emission from binary supermassive black holes (SMBHs) at the centres of galaxy merger remnants. In this paper, using the latest observational estimates for a range of galaxy properties and scaling relations, we predict the amplitude of the GW background generated by the binary SMBH population. We also predict the numbers of individual binary SMBH GW sources. We predict the characteristic strain amplitude of the GW background to lie in the range 5.1 × 10−16 < Ayr < 2.4 × 10−15 at a frequency of (1 yr)−1, with 95 per cent confidence. Higher values within this range, which correspond to the more commonly preferred choice of galaxy merger time-scale, will fall within the expected sensitivity ranges of existing PTA projects in the next few years. In contrast, we find that a PTA consisting of at least 100 pulsars observed with next-generation radio telescopes will be required to detect continuous-wave GWs from binary SMBHs. We further suggest that GW memory bursts from coalescing SMBH pairs are not viable sources for PTAs. Both the GW background and individual GW source counts are dominated by binaries formed in mergers between early-type galaxies of masses ≳5 × 1010 M⊙ at redshifts ≲1.5. Uncertainties in the galaxy merger time-scale and the SMBH mass–galaxy bulge mass relation dominate the uncertainty in our predictions.

It is believed that cores of neutron stars provide a natural laboratory where exotic high baryon density phases of quantum chromo dynamics (QCD) may exist. In fact, the theoretically well-established neutron superfluid phase is also believed to be found only inside neutron stars. Focus on neutrons stars has tremendously intensified in recent years with the direct detection of gravitational waves (GWs) by LIGO/Virgo from binary neutron star (BNS) merger events which has allowed the possibility of directly probing the properties of the interior of a neutron star. A truly remarkable phenomenon manifested by rapidly rotating neutron stars is in their avatar as Pulsars. The accuracy of pulsar timing can reach the level of one part in 10[Formula: see text], comparable to that of atomic clocks. Indeed, it was such a great accuracy which had allowed the first indirect detection of GWs from a BNS system. Such an incredible accuracy of pulse timings points to a very interesting possibility. Any deformation of the pulsar, even if it is extremely tiny, has the potential of leaving its imprints on the pulses through introduction of tiny perturbations in the entire moment of inertia (MI) tensor. While, the diagonal components of perturbed MI tensor affect the pulse timings, the off-diagonal components lead to wobbling of pulsar, directly affecting the pulse profile. This opens up a new window of opportunity for exploring various phase transitions occurring inside a pulsar core, through induced density fluctuations, which may be observable as perturbations in the pulse timing as well as its profile. Such perturbations also naturally induce a rapidly changing quadrupole moment of the star, thereby providing a new source of GW emission. Another remarkable possibility arises when we consider the effect of an external GW on neutron star. With the possibility of detecting any minute changes in its configuration through pulse observations, the neutron star has the potential of performing as a Weber detector of GW. This brief review will focus on these specific aspects of a pulsar. Specifically, the focus will be on the type of physics which can be probed by utilizing the effect of changes in the MI tensor of the pulsar on pulse properties.

Gravitational wave signal templates, pattern recognition, and reciprocal Eulerian gamma functions

The evolution and the stability of mass transfer of CO+He white dwarf (WD) binaries are not well understood. Observationally they may emerge as AM CVn binaries and are important gravitational wave (GW) emitters. In this work, we have modeled the evolution of double WD binaries with accretor masses of 0.50–1.30 M ⊙ and donor masses of 0.17–0.45 M ⊙ using the detailed stellar evolution code mesa. We find that the evolution of binaries with same donor masses but different accretor masses is very similar and binaries with same accretor masses but larger He donor masses have larger maximum mass transfer rates and smaller minimum orbital periods. We also demonstrate that the GW signal from AM CVn binaries can be detected by spaceborne GW observatories, such as LISA and TianQin. There is a linear relation between the donor mass and gravitational wave frequency during the mass transfer phase. In our calculation, all binaries can have dynamically stable mass transfer, which is very different from previous studies. The threshold donor mass of Eddington-limited mass transfer for a given accretor WD mass is lower than previous studies. Assuming that a binary may enter a common envelope if the mass transfer rate exceeds the maximum stable burning rate of He, we provide a new criterion for double WDs surviving mass transfer, which is below the threshold of the Eddington limit. Finally, we find that some systems with oxygen–neon (ONe) WDs in our calculation may evolve into detached binaries consisting of neutron stars and extremely low-mass He WDs, and further ultracompact X-ray binaries.

We estimate the long-lasting gravitational wave (GW) emission of compact dark objects following a binary neutron-star (NS) merger. We consider compact dark objects, which initially reside in the centers of NSs and which may consist of self-interacting dark matter (DM). By approximating the compact dark objects as test particles, we model the merging of NS binaries hosting DM components with three-dimensional relativistic simulations. Our simulation results suggest that the DM components remain gravitationally bound and orbit inside the merger remnant with orbital separations of typically a few km. The subsequent orbital motion of the DM components generates a GW signal with frequencies in the range of a few kHz. When considering a range of different binary masses and high-density equations of state (EOS), we find that the GW frequency of the orbiting DM components scales with the compactness of NSs. Similarly, we find relations between the DM GW frequency and the dominant postmerger GW frequency of the stellar fluid or the tidal deformability, which quantifies EOS effects during the binary inspiral. Hence, a measurement of these quantities can be used to specify the frequency range of the GW emission by DM. Under the assumption that GW back reaction is the only relevant dissipative process, the GW signal may last between seconds and years depending on the mass of the DM component. We estimate the detectibility of the GW signals and find that DM components in NS mergers may only be detectable with existing and projected GW instruments if the dark objects are as massive as about 0.01 to $0.1{M}_{\ensuremath{\bigodot}}$. We emphasize that the GW emission is limited by the lifetime of the remnant. A forming black hole will immediately swallow the DM objects because their orbits are smaller than the innermost stable circular orbit of the black hole.

Three periodically variable stars have recently been discovered (V407 Vul, P=9.5 min; ES Cet, P=10.3 min; RX J0806.3+1527, P=5.3 min) with properties that suggest that their photometric periods are also their orbital periods, making them the most compact binary stars known. If true, this might indicate that close, detached, double white dwarfs are able to survive the onset of mass transfer caused by gravitational wave radiation and emerge as the semi-detached, hydrogen-deficient stars known as the AM CVn stars. The accreting white dwarfs in such systems are large compared to the orbital separations. This has two effects. First, it makes it likely that the mass-transfer stream can hit the accretor directly. Secondly, it causes a loss of angular momentum from the orbit which can destabilize the mass transfer unless the angular momentum lost to the accretor can be transferred back to the orbit. The effect of the destabilization is to reduce the number of systems which survive mass transfer by as much as one hundredfold. In this paper we analyse this destabilization and the stabilizing effect of a dissipative torque between the accretor and the binary orbit. We obtain analytical criteria for the stability of both disc-fed and direct impact accretion, and we carry out numerical integrations to assess the importance of secondary effects, the chief one being that otherwise stable systems can exceed the Eddington accretion rate. We show that to have any effect upon survival rates, the synchronizing torque must act on a time-scale of the order of 1000 yr or less. If synchronization torques are this strong, then they will play a significant role in the spin rates of white dwarfs in cataclysmic variable stars as well.

Quickly localizing the identified white dwarf (WD) binaries is the basic requirement for the space-based gravitational wave (GW) detection. In fact, the amplitude of GW signals are modulated by the periodic motion of GW detectors on the solar orbit. The intensity of the observed signals is enhanced according to the observation time beyond a year to enhance a high signal to noise ratio (SNR). As data gap exists, the completeness of the data observed for a long time depends on filling gaps in the data. Actually, in a year period, the GW sources have a best observation orbit position of GW detectors, where the intensity of GW is maximum. Thus, the best positions can be searched for the verified GW sources of the sky map to enhance SNR too, which avoids filling data gaps. For the three arms response intensity of the GW signals changing more clearly with the location of the GW sources relative to the detector, the noises and the suppression of noise by time delay interferometer are ignored. As a verification case, the four WD binaries are chosen, whose best observation orbit positions of the GW detectors are related to the direction of WD binaries perpendicular to the detection arms. The two verification binaries: J0806 and V407 Vul are observed at the best orbit positions by TAIJI for the short time of 2 and 3 days respectively. The corresponding intensities of GWs are above the values of the TAIJI sensitivity curve, significantly. Location parameter estimation of the verification WD binaries are performed using the Metropolis–Hastings MCMC method. The confidence level of the parameters obtained in the best position is significantly several times higher than that in the worst position where the direction of WD binaries are almost parallel to the detection arms. Compared with a single detector, the network of two detectors improve slightly the accuracy of location of the verification binaries. These results imply that the searching of GW signals and parameter estimation of GW sources from the experimental data of the space-based mission do not ignore the orbit positions relevant to GW sources.