The population of Milky Way satellite galaxies is of great interest for cosmology, fundamental physics, and astrophysics. They represent the faint end of the galaxy luminosity function, are the most dark-matter dominated objects in the local Universe, and contain the oldest and most metal-poor stellar populations. Recent surveys have revealed around 60 satellites, but this could represent less than half of the total. Characterization of these systems remains a challenge due to their low luminosity. We consider the gravitational wave observatory LISA as a potential tool for studying these satellites through observations of their short-period double white dwarf populations. LISA will observe the entire sky without selection effects due to dust extinction, complementing optical surveys, and could potentially discover massive satellites hidden behind the disk of the galaxy.

Massive black hole binaries (MBHBs) of $10^5 \, \rm M_\odot - 3 \times 10^7 \, \rm M_\odot $ merging in low redshift galaxies ($z\le4$) are sufficiently loud to be detected weeks before coalescence with the Laser Interferometer Space Antenna (LISA). This allows us to perform the parameter estimation $on$ $the$ $fly$, i.e. as a function of the time to coalescence during the inspiral phase, relevant for early warning of the planned LISA protected periods and for searches of electromagnetic signals. In this work, we study the evolution of the sky position, luminosity distance, chirp mass and mass ratio uncertainties as function of time left before merger. Overall, light systems with total intrinsic mass $\rm M_{\rm tot} = 3 \times 10^5 \, \rm M_\odot$ are characterized by smaller uncertainties than heavy ones ($\rm M_{\rm tot} = 10^7 \, \rm M_\odot$) during the inspiral. Luminosity distance, chirp mass and mass ratio are well constrained at the end of the inspiral. Concerning sky position, at $z=1$, MBHBs with $\rm M_{\rm tot} = 3 \times 10^5 \, \rm M_\odot$ can be localized with a median precision of $\simeq 10^2 \, \rm deg^2 (\simeq 1 \, \rm deg^2)$ at 1 month (1 hour) from merger, while the sky position of heavy MBHBs can be determined to $10 \, \rm deg^2$ only 1 hour before merger. However the uncertainty around the median values broadens with time, ranging in between 0.04 -- 20 $\rm deg^2$ (0.3 -- 3 $\times 10^3 \, \rm deg^2$) for light (heavy) systems at 1 hour before merger. At merger the sky localization improves down to $\simeq 10^{-1} \, \rm deg^2$ for all masses. For the benefit of the observer community, we provide the full set of data from our simulations and simple and ready-to-use analytical fits to describe the time evolution of uncertainties in the aforementioned parameters, valid for systems with total mass between $10^5$--$10^7 \, \rm M_\odot$ and redshift $0.3$--$3$.

Milky Way dwarf satellites are unique objects that encode the early structure formation and therefore represent a window into the high redshift Universe. So far, their study was conducted using electromagnetic waves only. The future Laser Interpreter Space Antenna (LISA) has the potential to reveal Milky Way satellites in gravitational waves emitted by double white dwarf (DWD) binaries. We investigate gravitational wave (GW) signals detectable by LISA as a possible tool for the identification and characterisation of the Milky Way satellites. We use the binary population synthesis technique to model the population of DWDs in dwarf satellites and we assess the impact on the number of LISA detections when making changes to the total stellar mass, distance, star formation history and metallicity of satellites. We calibrate predictions for the known Milky Way satellites on their observed properties. We find that DWDs emitting at frequencies $\gtrsim 3\,$mHz can be detected in Milky Way satellites at large galactocentric distances. The number of these high frequency DWDs per satellite primarily depends on its mass, distance, age and star formation history, and only mildly depends on the other assumptions regarding their evolution such as metallicity. We find that dwarf galaxies with $M_\star>10^6\,$M$_{\odot}$ can host detectable LISA sources with a number of detections that scales linearly with the satellite's mass. We forecast that out of the known satellites, Sagittarius, Fornax, Sculptor and the Magellanic Clouds can be detected with LISA. As an all-sky survey that does not suffer from contamination and dust extinction, LISA will provide observations of the Milky Way and dwarf satellites galaxies valuable for Galactic archaeology and near-field cosmology.

We investigate the ability of the Laser Interferometer Space Antenna (LISA) to measure the center of mass acceleration of stellar-origin black hole binaries emitting gravitational waves. Our analysis is based on the idea that the acceleration of the center of mass induces a time variation in the redshift of the gravitational wave, which in turn modifies its waveform. We confirm that while the cosmological acceleration is too small to leave a detectable imprint on the gravitational waveforms observable by LISA, larger peculiar accelerations may be measurable for sufficiently long lived sources. We focus on stellar mass black hole binaries, which will be detectable at low frequencies by LISA and near coalescence by ground based detectors. These sources may have large peculiar accelerations, for instance, if they form in nuclear star clusters or in AGN accretion disks. If that is the case, we find that in an astrophysical population calibrated to the LIGO-Virgo observed merger rate, LISA will be able to measure the peculiar acceleration of a small but significant fraction of the events if the mission lifetime is extended beyond the nominal duration of 4 years. In this scenario LISA will be able to assess whether black hole binaries form close to galactic centers, particularly in AGN disks, and will thus help discriminate between different formation mechanisms. Although for a nominal 4 years LISA mission the peculiar acceleration effect cannot be measured, a consistent fraction of events may be biased by strong peculiar accelerations which, if present, may imprint large systematic errors on some waveform parameters. In particular, estimates of the luminosity distance could be strongly biased and consequently induce large systematic errors on LISA measurements of the Hubble constant with stellar mass black hole binaries.

One of the central goals of LISA is the detection of gravitational waves from the merger of supermassive black holes. Contrary to stellar-mass black hole mergers, such events are expected to be rich X-ray sources due to the accretion of material from the circumbinary disks onto the black holes. The orbital dynamics before merger is also expected to modulate the resulting X-ray emission via Doppler boosting in a quasi-periodic way, and in a simple phase relation with the gravitational wave from the inspiral of the black holes. Detecting the X-ray source would enable a precise and early localization of the binary, thus allowing many telescopes to observe the very moment of the merger. Although identifying the correct X-ray source in the relatively large LISA sky localization will be challenging due to the presence of many confounding point sources, the quasi-periodic modulation may aid in the identification. We explore the practical feasibility of such idea. We simulate populations of merging supermassive black holes, their detection with LISA and their X-ray lightcurves using a simple model. Taking the parameters of the X-ray Telescope on the proposed NASA Transient Astrophysics Probe, we then design and simulate an observation campaign that searches for the modulated X-ray source while LISA is still observing the inspiral of the black holes. Assuming a fiducial LISA detection rate of $10$ mergers per year at redshift closer than $3.5$, we expect a few detections of modulated X-ray counterparts over the nominal duration of the LISA mission.

The space-based detector LISA may observe gravitational waves from the early inspiral of stellar-mass black hole binaries, some of which could have significant eccentricity. Current gravitational waveform templates are only valid for small orbital velocities (i.e., in a post-Newtonian expansion) and small initial eccentricity $e_0$ ("post-circular" expansion). We conventionally define $e_0$ as the eccentricity corresponding to an orbital frequency of $5 \text{ mHz}$, and we study the convergence properties of frequency-domain inspiral templates that are accurate up to 2PN and order $e_0^6$ in eccentricity. We compute the so-called "unfaithfulness" between the full template and "reduced" templates obtained by dropping some terms in the phasing series; we investigate the conditions under which systematic errors are negligible with respect to statistical errors, and we study the convergence properties of statistical errors. In general, eccentric waveforms lead to larger statistical errors than circular waveforms due to correlations between the parameters, but the error estimates do not change significantly as long as we include terms of order $e_0^2$ or higher in the phasing.

The progenitors of the high-mass black-hole mergers observed by LIGO and Virgo are potential LISA sources and promising candidates for multiband GW observations. In this letter, we consider the minimum signal-to-noise ratio these sources must have to be detected by LISA. Our revised threshold of $\rho_{\rm thr}\sim 15$ is higher than previous estimates, which significantly reduces the expected number of events. We also point out the importance of the detector performance at high-frequencies and the duration of the LISA mission, which both influence the event rate substantially.

Modifications of General Relativity leave their imprint both on the cosmic expansion history through a non-trivial dark energy equation of state, and on the evolution of cosmological perturbations in the scalar and in the tensor sectors. In particular, the modification in the tensor sector gives rise to a notion of gravitational-wave (GW) luminosity distance, different from the standard electromagnetic luminosity distance, that can be studied with standard sirens at GW detectors such as LISA or third-generation ground based experiments. We discuss the predictions for modified GW propagation from some of the best studied theories of modified gravity, such as Horndeski or the more general degenerate higher order scalar-tensor (DHOST) theories, non-local infrared modifications of gravity, bigravity theories and the corresponding phenomenon of GW oscillation, as well as theories with extra or varying dimensions. We show that modified GW propagation is a completely generic phenomenon in modified gravity. We then use a simple parametrization of the effect in terms of two parameters $(\Xi_0,n)$, that is shown to fit well the results from a large class of models, to study the prospects of observing modified GW propagation using supermassive black hole binaries as standard sirens with LISA. We construct mock source catalogs and perform detailed Markov Chain Monte Carlo studies of the likelihood obtained from LISA standard sirens alone, as well as by combining them with CMB, BAO and SNe data to reduce the degeneracies between cosmological parameters. We find that the combination of LISA with the other cosmological datasets allows one to measure the parameter $\Xi_0$ that characterizes modified GW propagation to the percent level accuracy, sufficient to test several modified gravity theories. [Abridged]

We build two families of inspiral waveforms for precessing binaries on eccentric orbits in the Fourier domain. To achieve this, we use a small eccentricity expansion of the waveform amplitudes in order to separate the periastron precession timescale from the orbital timescale, and use a shifted uniform asymptotics transformation to compute the Fourier transform in the presence of spin-induced precession. We show that the resulting waveforms can yield a median faithfulness above 0.993 when compared to an equivalent time domain waveform with an initial eccentricity of $e_0 \approx 0.3$. We also show that when the spins are large, using a circular waveform can potentially lead to significant biases in the recovery of the parameters, even when the system has fully circularized, particularly when the accumulated number of cycles is large. This is an effect of the residual eccentricity present when the objects forming the binary have nonvanishing spin components in the orbital plane.

The Laser Interferometer Space Antenna (LISA) will observe black hole binaries of stellar origin during their gravitational wave inspiral, months to years before coalescence. Due to the long duration of the signal in the LISA band, a faithful waveform is necessary in order to keep track of the binary phase. This is crucial to extract the signal from the data and to perform an unbiased estimation of the source parameters. We consider Post-Newtonian (PN) waveforms, and analyze the PN order needed to keep the bias caused by the PN approximation negligible relative to the statistical parameter estimation error, as a function of the source parameters. By considering realistic population models, we conclude that for $\sim 90\%$ of the stellar black hole binaries detectable by LISA, waveforms at low Post-Newtonian (PN) order (PN $\le 2$) are sufficiently accurate for an unbiased recovery of the source parameters. Our results provide a first estimate of the trade-off between waveform accuracy and information recovery for this class of LISA sources.

The space-based Laser Interferometer Space Antenna (LISA) will be able to observe the gravitational-wave signals from systems comprised of a massive black hole and a stellar-mass compact object. These systems are known as extreme-mass-ratio inspirals (EMRIs) and are expected to complete $\sim 10^4$-$10^5$ cycles in band, thus allowing exquisite measurements of their parameters. In this work, we attempt to quantify the astrophysical uncertainties affecting the predictions for the number of EMRIs detectable by LISA, and find that competing astrophysical assumptions produce a variance of about three orders of magnitude in the expected intrinsic EMRI rate. However, we find that irrespective of the astrophysical model, at least a few EMRIs per year should be detectable by the LISA mission, with up to a few thousands per year under the most optimistic astrophysical assumptions. We also investigate the precision with which LISA will be able to extract the parameters of these sources. We find that typical fractional statistical errors with which the intrinsic parameters (redshifted masses, massive black hole spin and orbital eccentricity) can be recovered are $\sim 10^{-6}$-$10^{-4}$. Luminosity distance (which is required to infer true masses) is inferred to about $10\%$ precision and sky position is localized to a few square degrees, while tests of the multipolar structure of the Kerr metric can be performed to percent-level precision or better.

The coalescence of compact objects is one of the most promising sources of gravitational waves for ground-based interferometric detectors, such as advanced LIGO and Virgo. Generically, com- pact objects in binaries are expected to be spinning with spin angular momenta misaligned with the orbital angular momentum, causing the orbital plane to precess. This precession adds rich structure to the gravitational waves, introducing such complexity that an analytic closed-form description has been unavailable until now. We here construct the first closed-form frequency- domain gravitational waveforms that are valid for generic spin-precessing quasicircular compact binary inspirals. We first construct time-domain gravitational waves by solving the post-Newtonian precession equations of motion with radiation reaction through multiple scale analysis. We then Fourier transform these time-domain waveforms with the method of shifted uniform asymptotics to obtain closed-form expressions for frequency-domain waveforms. We study the accuracy of these analytic, frequency-domain waveforms relative to waveforms obtained by numerically evolving the post-Newtonian equations of motion and find that they are suitable for unbiased parameter estimation for 99.2%(94.6%) of the binary configurations we studied at a signal-to-noise ratio of 10(25). These new frequency-domain waveforms could be used for detection and parameter estimation studies due to their accuracy and low computational cost.

Ultralight scalar fields around spinning black holes can trigger superradiant instabilities, forming a long-lived bosonic condensate outside the horizon. We use numerical solutions of the perturbed field equations and astrophysical models of massive and stellar-mass black hole populations to compute, for the first time, the stochastic gravitational-wave background from these sources. In optimistic scenarios the background is observable by Advanced LIGO and LISA for field masses $m_s$ in the range $\sim [2\times 10^{-13}, 10^{-12}]\,{\rm eV}$ and $\sim 5\times[ 10^{-19}, 10^{-16}]\,{\rm eV}$, respectively, and it can affect the detectability of resolvable sources. Our estimates suggest that an analysis of the stochastic background limits from LIGO O1 might already be used to marginally exclude axions with mass $\sim 10^{-12.5}{\rm eV}$. Semicoherent searches with Advanced LIGO (LISA) should detect $\sim 15~(5)$ to $200~(40)$ resolvable sources for scalar field masses $3\times 10^{-13}$ ($10^{-17}$) eV. LISA measurements of massive BH spins could either rule out bosons in the range $\sim [10^{-18}, 2\times 10^{-13}]$ eV, or measure $m_s$ with ten percent accuracy in the range $\sim[10^{-17}, 10^{-13}]$ eV.

Ultralight bosons can induce superradiant instabilities in spinning black holes, tapping their rotational energy to trigger the growth of a bosonic condensate. Possible observational imprints of these boson clouds include (i) direct detection of the nearly monochromatic (resolvable or stochastic) gravitational waves emitted by the condensate, and (ii) statistically significant evidence for the formation of "holes" at large spins in the spin versus mass plane (sometimes also referred to as "Regge plane") of astrophysical black holes. In this work, we focus on the prospects of LISA and LIGO detecting or constraining scalars with mass in the range $m_s\in [10^{-19},\,10^{-15}]$ eV and $m_s\in [10^{-14},\,10^{-11}]$ eV, respectively. Using astrophysical models of black-hole populations calibrated to observations and black-hole perturbation theory calculations of the gravitational emission, we find that, in optimistic scenarios, LIGO could observe a stochastic background of gravitational radiation in the range $m_s\in [2\times 10^{-13}, 10^{-12}]$ eV, and up to $10^4$ resolvable events in a $4$-year search if $m_s\sim 3\times 10^{-13}\,{\rm eV}$. LISA could observe a stochastic background for boson masses in the range $m_s\in [5\times 10^{-19}, 5\times 10^{-16}]$, and up to $\sim 10^3$ resolvable events in a $4$-year search if $m_s\sim 10^{-17}\,{\rm eV}$. LISA could further measure spins for black-hole binaries with component masses in the range $[10^3, 10^7]~M_\odot$, which is not probed by traditional spin-measurement techniques. A statistical analysis of the spin distribution of these binaries could either rule out scalar fields in the mass range $\sim [4 \times 10^{-18}, 10^{-14}]$ eV, or measure $m_s$ with ten percent accuracy if light scalars in the mass range $\sim [10^{-17}, 10^{-13}]$ eV exist.

Knights Landing (KNL) is the code name for the second-generation Intel Xeon Phi product family. KNL has generated significant interest in the data analysis and machine learning communities because its new many-core architecture targets both of these workloads. The KNL many-core vector processor design enables it to exploit much higher levels of parallelism. At the Lincoln Laboratory Supercomputing Center (LLSC), the majority of users are running data analysis applications such as MATLAB and Octave. More recently, machine learning applications, such as the UC Berkeley Caffe deep learning framework, have become increasingly important to LLSC users. Thus, the performance of these applications on KNL systems is of high interest to LLSC users and the broader data analysis and machine learning communities. Our data analysis benchmarks of these application on the Intel KNL processor indicate that single-core double-precision generalized matrix multiply (DGEMM) performance on KNL systems has improved by ~3.5x compared to prior Intel Xeon technologies. Our data analysis applications also achieved ~60% of the theoretical peak performance. Also a performance comparison of a machine learning application, Caffe, between the two different Intel CPUs, Xeon E5 v3 and Xeon Phi 7210, demonstrated a 2.7x improvement on a KNL node.

We investigate the capability of various configurations of the space interferometer eLISA to probe the late-time background expansion of the universe using gravitational wave standard sirens. We simulate catalogues of standard sirens composed by massive black hole binaries whose gravitational radiation is detectable by eLISA, and which are likely to produce an electromagnetic counterpart observable by future surveys. The main issue for the identification of a counterpart resides in the capability of obtaining an accurate enough sky localisation with eLISA. This seriously challenges the capability of four-link (2 arm) configurations to successfully constrain the cosmological parameters. Conversely, six-link (3 arm) configurations have the potential to provide a test of the expansion of the universe up to $z\sim 8$ which is complementary to other cosmological probes based on electromagnetic observations only. In particular, in the most favourable scenarios, they can provide a significant constraint on $H_0$ at the level of 0.5%. Furthermore, $(\Omega_M, \Omega_\Lambda)$ can be constrained to a level competitive with present SNIa results. On the other hand, the lack of massive black hole binary standard sirens at low redshift allows to constrain dark energy only at the level of few percent.

The recent Advanced LIGO detections of coalescing black hole binaries (BHBs) imply a large population of such systems emitting at milli-Hz frequencies, accessible to the Laser Interferometer Space Antenna (LISA). We show that these systems provide a new class of cosmological standard sirens. Direct LISA luminosity distance -$D_l$- measurements, combined with the inhomogeneous redshift -$z$- distribution of possible host galaxies provide an effective way to populate the $D_l-z$ diagram at $z<0.1$, thus allowing a precise local measurement of the Hubble expansion rate. To be effective, the method requires a sufficiently precise LISA distance determination and sky localization of a sizeable number of BHBs, which is best achieved for a 6-link detector configuration. We find that, for a BHB population consistent with current fiducial LIGO rates, the Hubble constant $H_0$ can be determined at the $\sim$5% and $\sim$2% level (68% confidence) assuming two and five million Km arm-length respectively.

Following the selection of The Gravitational Universe by ESA, and the successful flight of LISA Pathfinder, the LISA Consortium now proposes a 4 year mission in response to ESA's call for missions for L3. The observatory will be based on three arms with six active laser links, between three identical spacecraft in a triangular formation separated by 2.5 million km. LISA is an all-sky monitor and will offer a wide view of a dynamic cosmos using Gravitational Waves as new and unique messengers to unveil The Gravitational Universe. It provides the closest ever view of the infant Universe at TeV energy scales, has known sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales near the horizons of black holes, all the way to cosmological scales. The LISA mission will scan the entire sky as it follows behind the Earth in its orbit, obtaining both polarisations of the Gravitational Waves simultaneously, and will measure source parameters with astrophysically relevant sensitivity in a band from below $10^{-4}\,$Hz to above $10^{-1}\,$Hz.

Binary systems of two compact objects circularize and spiral toward each other via the emission of gravitational waves. The coupling of the spins of each object with the orbital angular momentum causes the orbital plane to precess, which leads to modulation of the gravitational wave signal. Until now, generating frequency-domain waveforms for fully precessing systems for use in gravitational wave data analysis meant numerically integrating the equations of motion, then Fourier transforming the result, which is very computationally intensive for systems that complete hundreds or thousands of cycles in the sensitive band of a detector. Previously, analytic solutions were only available for certain special cases or for simplified models. Here we describe the construction of closed-form, frequency-domain waveforms for fully-precessing, quasi-circular binary inspirals.

Up to hundreds of black hole binaries individually resolvable by eLISA will coalesce in the Advanced LIGO/Virgo band within ten years, allowing for multi-band gravitational wave observations. Binaries formed via dynamical interactions in dense star clusters are expected to have eccentricities $e_0\sim 10^{-3}$-$10^{-1}$ at the frequencies $f_0=10^{-2}$ Hz where eLISA is most sensitive, while binaries formed in the field should have negligible eccentricity in both frequency bands. We estimate that eLISA should always be able to detect a nonzero $e_0$ whenever $e_0\gtrsim 10^{-2}$; if $e_0\sim 10^{-3}$, eLISA should detect nonzero eccentricity for a fraction $\sim 90\%$ ($\sim 25\%$) of binaries when the observation time is $T_{\rm obs}=5$ ($2$) years, respectively. Therefore eLISA observations of BH binaries have the potential to distinguish between field and cluster formation scenarios.

A space-based interferometer such as eLISA could observe few to few thousands progenitors of black hole binaries (BHBs) similar to those recently detected by Advanced LIGO. Gravitational radiation circularizes the orbit during inspiral, but some BHBs retain a measurable eccentricity at the low frequencies where eLISA is most sensitive. The eccentricity of a BHB carries precious information about its formation channel: BHBs formed in the field, in globular clusters, or close to a massive black hole (MBH) have distinct eccentricity distributions in the eLISA band. We generate mock eLISA observations, folding in measurement errors, and using Bayesian model selection we study whether eLISA measurements can identify the BHB formation channel. We find that a handful of observations would suffice to tell whether BHBs were formed in the gravitational field of a MBH. Conversely, several tens of observations are needed to tell apart field formation from globular cluster formation. A five-year eLISA mission with the longest possible armlength is desirable to shed light on BHB formation scenarios.

We compare the science capabilities of different eLISA mission designs, including four-link (two-arm) and six-link (three-arm) configurations with different arm lengths, low-frequency noise sensitivities and mission durations. For each of these configurations we consider a few representative massive black hole formation scenarios. These scenarios are chosen to explore two physical mechanisms that greatly affect eLISA rates, namely (i) black hole seeding, and (ii) the delays between the merger of two galaxies and the merger of the black holes hosted by those galaxies. We assess the eLISA parameter estimation accuracy using a Fisher matrix analysis with spin-precessing, inspiral-only waveforms. We quantify the information present in the merger and ringdown by rescaling the inspiral-only Fisher matrix estimates using the signal-to-noise ratio from non-precessing inspiral-merger-ringdown phenomenological waveforms, and from a reduced set of precessing numerical relativity/post-Newtonian hybrid waveforms. We find that all of the eLISA configurations considered in our study should detect some massive black hole binaries. However, configurations with six links and better low-frequency noise will provide much more information on the origin of black holes at high redshifts and on their accretion history, and they may allow the identification of electromagnetic counterparts to massive black hole mergers.

One century after its formulation, Einstein's general relativity has made remarkable predictions and turned out to be compatible with all experimental tests. Most of these tests probe the theory in the weak-field regime, and there are theoretical and experimental reasons to believe that general relativity should be modified when gravitational fields are strong and spacetime curvature is large. The best astrophysical laboratories to probe strong-field gravity are black holes and neutron stars, whether isolated or in binary systems. We review the motivations to consider extensions of general relativity. We present a (necessarily incomplete) catalog of modified theories of gravity for which strong-field predictions have been computed and contrasted to Einstein's theory, and we summarize our current understanding of the structure and dynamics of compact objects in these theories. We discuss current bounds on modified gravity from binary pulsar and cosmological observations, and we highlight the potential of future gravitational wave measurements to inform us on the behavior of gravity in the strong-field regime.

Gravitational waves from neutron star binary inspirals contain information about the equation of state of supranuclear matter. In the absence of definitive experimental evidence that determines the correct equation of state, a number of diverse models that give the pressure in a neutron star as function of its density have been proposed. These models differ not only in the approximations and techniques they use to solve the many-body Schr\"odinger equation, but also in the neutron star composition they assume. We study whether gravitational wave observations of neutron star binaries in quasicircular inspirals will allow us to distinguish between equations of state of differing internal composition, thereby providing important information about the properties of extremely high density matter. We carry out a Bayesian model selection analysis, and find that second generation gravitational wave detectors can heavily constrain equations of state that contain only quark matter, but hybrid stars containing both normal and quark matter are harder to distinguish from normal matter stars. A gravitational wave detection with a signal-to-noise ratio of 30 and masses around $1.4M_{\odot}$ could either detect or rule out strange quark stars with a 20 to 1 confidence. The presence of kaon condensates or hyperons in neutron star inner cores cannot be easily confirmed. For example, for the equations of state studied in this paper, even a gravitational wave signal with a signal-to-noise ratio as high as 60 would not allow us to claim a detection of kaon condensates or hyperons with confidence greater than 5 to 1. On the other hand, if kaon condensates and hyperons do not form in neutron stars, a gravitational wave signal with similar signal-to-noise ratio would be able to constrain their existence with an 11 to 1 confidence for high-mass systems.

Binary black holes on quasicircular orbits with spins aligned with their orbital angular momentum have been testbeds for analytic and numerical relativity for decades, not least because symmetry ensures that such configurations are equilibrium solutions to the spin-precession equations. In this work, we show that these solutions can be unstable when the spin of the higher-mass black hole is aligned with the orbital angular momentum and the spin of the lower-mass black hole is anti-aligned. Spins in these configurations are unstable to precession to large misalignment when the binary separation $r$ is between the values $r_{\rm ud\pm}= (\sqrt{\chi_1} \pm \sqrt{q \chi_2})^4 (1-q)^{-2} M$, where $M$ is the total mass, $q \equiv m_2/m_1$ is the mass ratio, and $\chi_1$ ($\chi_2$) is the dimensionless spin of the more (less) massive black hole. This instability exists for a wide range of spin magnitudes and mass ratios and can occur in the strong-field regime near merger. We describe the origin and nature of the instability using recently developed analytical techniques to characterize fully generic spin precession. This instability provides a channel to circumvent astrophysical spin alignment at large binary separations, allowing significant spin precession prior to merger affecting both gravitational-wave and electromagnetic signatures of stellar-mass and supermassive binary black holes.

Mergers of compact stellar remnant are prime targets for the LIGO/Virgo gravitational wave detectors. One hopes that the gravitational wave signals from these merger events can be used to study the mass and spin distribution of stellar remnants, and provide information about black hole horizons and the material properties of neutron stars. However, it has been suggested that degeneracies in the way that the star's mass and spin are imprinted in the waveforms may make it impossible to distinguish between black holes and neutron stars. Here we show that the precession of the orbital plane due to spin-orbit coupling breaks the mass-spin degeneracy, and allows us to distinguish between standard neutron stars and alternative possibilities, such as black holes or exotic neutron stars with large masses and spins.

The detection and characterization of gravitational wave signals from compact binary coalescence events relies on accurate waveform templates in the frequency domain. The stationary phase approximation (SPA) can be used to compute closed-form frequency-domain waveforms for non-precessing, quasi-circular binary inspirals. However, until now, no fast frequency-domain waveforms have existed for generic, spin-precessing quasi-circular compact binary inspirals. Templates for these systems have had to be computed via a discrete Fourier transform of finely sampled time-domain waveforms, which is far more computationally expensive than those constructed directly in the frequency-domain, especially for those systems that are dominated by the inspiral part. There are two obstacles to deriving frequency-domain waveforms for precessing systems: (i) the spin-precession equations do not admit closed-form solutions for generic systems; (ii) the SPA fails catastrophically. Presently there is no general solution to the first problem, so we must resort to numerical integration of the spin precession equations. This is not a significant obstacle, as numerical integration on the slow precession timescale adds very little to the computational cost of generating the waveforms. Our main result is to solve the second problem, by providing an alternative to the SPA that we call the method of Shifted Uniform Asymptotics, or SUA, that cures the divergences in the SPA. The construction of frequency-domain templates using the SUA can be orders of magnitude more efficient than the time-domain ones obtained through a discrete Fourier transform. Moreover, this method is very faithful to the discrete Fourier transform, with mismatches on the order of $10^{-5}$.

The angular dependence of the gravitational radiation emitted in compact binary mergers and gravitational collapse is usually separated using spin-weighted spherical harmonics ${}_sY_{\ell m}$ of spin weight $s$, that reduce to the ordinary spherical harmonics $Y_{\ell m}$ when $s=0$. Teukolsky first showed that the perturbations of the Kerr black hole that may be produced as a result of these events are separable in terms of a different set of angular functions: the spin-weighted spheroidal harmonics ${}_sS_{\ell m n}$, where $n$ denotes the "overtone index" of the corresponding Kerr quasinormal mode frequency $\omega_{\ell m n}$. In this paper we compute the complex-valued scalar products of the ${}_sS_{\ell m n}$'s with the ${}_sY_{\ell m}$'s ("spherical-spheroidal mixing coefficients") and with themselves ("spheroidal-spheroidal mixing coefficients") as functions of the dimensionless Kerr parameter $j$. Tables of these coefficients and analytical fits of their dependence on $j$ are available online for use in gravitational-wave source modeling and in other applications of black-hole perturbation theory.

We calculate analytical gravitational waveforms in the time- and frequency-domain for precessing quasi-circular binaries with spins of arbitrary magnitude, but nearly aligned with the orbital angular momentum. We first derive an analytical solution to the precession equations by expanding in the misalignment angle and using multiple scale analysis to separate timescales. We then use uniform asymptotic expansions to analytically Fourier transform the time-domain waveform, thus extending the stationary-phase approximation, which fails when precession is present. The resulting frequency-domain waveform family has a high overlap with numerical waveforms obtained by direct integration of the post-Newtonian equations of motion and discrete Fourier transformations. Such a waveform family lays the foundations for the accurate inclusion of spin precession effects in analytical gravitational waveforms, and thus, it can aid in the detection and parameter estimation of gravitational wave signals from the inspiral phase of precessing binary systems.

Unbound interacting compact binaries emit gravitational radiation in a wide frequency range. Since short burst-like signals are expected in future detectors, such as LISA or advanced LIGO, it is interesting to study their energy spectrum and the position of the frequency peak. Here we derive them for a system of massive objects interacting on hyperbolic orbits within the quadrupole approximation, following the work of Capozziello et al. In particular, we focus on the derivation of an analytic formula for the energy spectrum of the emitted waves. Within numerical approximation our formula is in agreement with the two known limiting cases: for the eccentricity {\epsilon} = 1, the parabolic case, whose spectrum was computed by Berry and Gair, and the large {\epsilon} limit with the formula given by Turner.

We review the expected science performance of the New Gravitational-Wave Observatory (NGO, a.k.a. eLISA), a mission under study by the European Space Agency for launch in the early 2020s. eLISA will survey the low-frequency gravitational-wave sky (from 0.1 mHz to 1 Hz), detecting and characterizing a broad variety of systems and events throughout the Universe, including the coalescences of massive black holes brought together by galaxy mergers; the inspirals of stellar-mass black holes and compact stars into central galactic black holes; several millions of ultracompact binaries, both detached and mass transferring, in the Galaxy; and possibly unforeseen sources such as the relic gravitational-wave radiation from the early Universe. eLISA's high signal-to-noise measurements will provide new insight into the structure and history of the Universe, and they will test general relativity in its strong-field dynamical regime.

The emission of gravitational waves is studied for a system of massive objects interacting on hyperbolic orbits within the quadrupole approximation following the work of Capozziello et al. Here we focus on the derivation of an analytic formula for the energy spectrum of the emitted waves. We checked numerically that our formula is in agreement with the two limiting cases for which results were already available: for the eccentricity {\epsilon} = 1, the parabolic case whose spectrum was computed by Berry and Gair, and the large {\epsilon} limit with the formula given by Turner.

This document introduces the exciting and fundamentally new science and astronomy that the European New Gravitational Wave Observatory (NGO) mission (derived from the previous LISA proposal) will deliver. The mission (which we will refer to by its informal name "eLISA") will survey for the first time the low-frequency gravitational wave band (about 0.1 mHz to 1 Hz), with sufficient sensitivity to detect interesting individual astrophysical sources out to z = 15. The eLISA mission will discover and study a variety of cosmic events and systems with high sensitivity: coalescences of massive black holes binaries, brought together by galaxy mergers; mergers of earlier, less-massive black holes during the epoch of hierarchical galaxy and black-hole growth; stellar-mass black holes and compact stars in orbits just skimming the horizons of massive black holes in galactic nuclei of the present era; extremely compact white dwarf binaries in our Galaxy, a rich source of information about binary evolution and about future Type Ia supernovae; and possibly most interesting of all, the uncertain and unpredicted sources, for example relics of inflation and of the symmetry-breaking epoch directly after the Big Bang. eLISA's measurements will allow detailed studies of these signals with high signal-to-noise ratio, addressing most of the key scientific questions raised by ESA's Cosmic Vision programme in the areas of astrophysics and cosmology. They will also provide stringent tests of general relativity in the strong-field dynamical regime, which cannot be probed in any other way. This document not only describes the science but also gives an overview on the mission design and orbits.

We compute here the spin-orbit and spin-spin couplings needed for an accurate computation of the phasing of gravitational waves emitted by comparable-mass binaries on eccentric orbits at the second post-Newtonian (PN) order. We use a quasi-Keplerian parametrization of the orbit free of divergencies in the zero eccentricity limit. We find that spin-spin couplings induce a residual eccentricity for coalescing binaries at 2PN, of the order of $10^{-4}$-$10^{-3}$ for supermassive black hole binaries in the LISA band. Spin-orbit precession also induces a non-trivial pattern in the evolution of the eccentricity, which could help to reduce the errors on the determination of the eccentricity and spins in a gravitational wave measurement.

With one exception, previous analyses of the measurement accuracy of gravitational wave experiments for comparable-mass binary systems have neglected either spin-precession effects or subdominant harmonics and amplitude modulations. Here we give the first explicit description of how these effects combine to improve parameter estimation. We consider supermassive black hole binaries as expected to be observed with the planned space-based interferometer LISA, and study the measurement accuracy for several astrophysically interesting parameters obtainable taking into account the full 2PN waveform for spinning bodies, as well as spin-precession effects. We find that for binaries with a total mass in the range 10^5 M_Sun < M < 10^7 M_Sun at a redshift of 1, a factor ~1.5 is in general gained in accuracy, with the notable exception of the determination of the individual masses in equal-mass systems, for which a factor ~5 can be gained. We also find, as could be expected, that using the full waveform helps increasing the upper mass limit for detection, which can be as high as M = 10^8 M_Sun at a redshift of 1, as well as the redshift limit where some information can be extracted from a system, which is roughly z = 10 for M < 10^7 M_Sun, 1.5-5 times higher than with the restricted waveform. We computed that the full waveform allows to use supermassive black hole binaries as standard sirens up to a redshift of z = 1.6, about 0.4 larger than what previous studies allowed. We found that for lower unequal-mass binary systems, the measurement accuracy is not as drastically improved as for other systems. This suggests that for these systems, adding parameters such as eccentricity or alternative gravity parameters could be achieved without much loss in the accuracy.