We study light variability of gravitationally magnified high-redshift star clusters induced by a foreground population of microlenses. This arises as the incoherent superposition of light variations from a large number of source stars traversing the random magnification pattern on the source plane. The light curve resembles a scale-invariant, Gaussian process on timescales of years to decades, while exhibits rapid and frequent micro-caustic crossing flares of larger amplitudes on timescales of days to months. For a concrete example, we study a young Lyman-continuum-leaking star cluster recently discovered in the lensed Sunburst Arc at $z=2.37$. We show that one magnified image happens to be intervened by a faint foreground galaxy, and hence should exhibit a variable flux at the $1$--$2\%$ level, which is measurable in space with $\sim 1$--$3\,$ks exposures on the Hubble Space Telescope and more easily with the James Webb Space Telescope, or from the ground using a $\sim$4-meter telescope without adaptive optics. Detailed measurement of this variability will enable us to determine the absolute macro magnification and hence the intrinsic mass and length scales of the star cluster, test synthetic models of stellar population, and probe multiplicity of massive stars. We furthermore suggest that monitoring the other lensed images of the star cluster, which are free from significant intervention by foreground stellar microlenses, will allow us to probe planetary to stellar mass compact objects constituting as little as just a few percent of the dark matter. Given the typical surface density of intracluster stars, we expect this phenomenon to be relevant for many other gravitationally magnified star clusters at Cosmic Noon behind galaxy cluster lenses.

Numerical simulations have historically played a major role in understanding the hydrodynamics of the tidal disruption process. Given the complexity of the geometry of the system, the challenges posed by the problem have indeed stimulated much work on the numerical side. Smoothed Particles Hydrodynamics methods, for example, have seen their very first applications in the context of tidal disruption and still play a major role to this day. Likewise, initial attempts at simulating the evolution of the disrupted star with the so-called affine method have been historically very useful. In this Chapter, we provide an overview of the numerical techniques used in the field and of their limitations, and summarize the work that has been done to simulate numerically the tidal disruption process.

We study the highly magnified arc SGAS J122651.3+215220 caused by a star-forming galaxy at $z_s=2.93$ crossing the lensing caustic cast by the galaxy cluster SDSS J1226+2152 ($z_l=0.43$), using Hubble Space Telescope observations. We report in the arc several asymmetric surface brightness features whose angular separations are a fraction of an arcsecond from the lensing critical curve and appear to be highly but unequally magnified image pairs of underlying compact sources, with one brightest pair having clear asymmetry consistently across four filters. One explanation of unequal magnification is microlensing by intracluster stars, which induces independent flux variations in the images of individual or groups of source stars in the lensed galaxy. For a second possibility, intracluster dark matter subhalos invisible to telescopes effectively perturb lensing magnifications near the critical curve and give rise to persistently unequal image pairs. Our modeling suggests, at least for the most prominent identified image pair, that the microlensing hypothesis is in tension with the absence of notable asymmetry variation over a six-year baseline, while subhalos of $\sim 10^6$--$10^8\,M_\odot$ anticipated from structure formation with Cold Dark Matter typically produce stationary and sizable asymmetries. We judge that observations at additional times and more precise lens models are necessary to stringently constrain temporal variability and robustly distinguish between the two explanations. The arc under this study is a scheduled target of a Director's Discretionary Early Release Science program of the James Webb Space Telescope, which will provide deep images and a high-resolution view with integral field spectroscopy.

We report the detection of new binary black hole merger events in the publicly available data from the second observing run of advanced LIGO and advanced Virgo (O2). The mergers were discovered using the new search pipeline described in Venumadhav et al. [Phys. Rev. D 100, 023011 (2019)], and are above the detection thresholds as defined in Abbott et al. (LIGO Scientific and Virgo Collaborations) [Phys. Rev. X 9, 031040 (2019).] Three of the mergers (GW170121, GW170304, GW170727) have inferred probabilities of being of astrophysical origin $p_{\rm astro} > 0.98$. The remaining three (GW170425, GW170202, GW170403) are less certain, with $p_{\rm astro}$ ranging from 0.5 to 0.8. The newly found mergers largely share the statistical properties of previously reported events, with the exception of GW170403, the least secure event, which has a highly negative effective spin parameter $\chi_{\rm eff}$ . The most secure new event, GW170121 ($p_{\rm astro} > 0.99$), is also notable due to its inferred negative value of $\chi_{\rm eff}$, which is inconsistent with being positive at the ~95.8% confidence level. The new mergers nearly double the sample of gravitational wave events reported from O2, and present a substantial opportunity to explore the statistics of the binary black hole population in the Universe. The number of detected events is not surprising since we estimate that the detection volume of our pipeline may be larger than that of other pipelines by as much as a factor of two (with significant uncertainties in the estimate). The increase in volume is larger when the constituent detectors of the network have very different sensitivities, as is likely to be the case in current and future runs.

We perform a detailed parameter estimation study of binary black hole merger events reported in Zackay et al. and Venumadhav et al.. These are some of the faintest signals reported so far, and hence, relative to the loud events in the GWTC-1 catalog, the data should have lesser constraining power on their intrinsic parameters. Hence we examine the robustness of parameter inference to choices made in the analysis, as well as any potential systematics. We check the impact of different methods of estimating the noise power spectral density, different waveform models, and different priors for the compact object spins. For most of the events, the resulting differences in the inferred values of the parameters are much smaller than their statistical uncertainties. The estimation of the effective spin parameter $\chi_{\mathrm{eff}}$, i.e. the projection of the mass-weighted total spin along the angular momentum, can be sensitive to analysis choices for two of the sources with the largest effective spin magnitudes, GW151216 and GW170403. The primary differences arise from using a 3D isotropic spin prior: the tails of the posterior distributions should be interpreted with care and due consideration of the other data analysis choices.

We present light curves and spectra of the tidal disruption event (TDE) ASASSN-18pg / AT 2018dyb spanning a period of one year. The event shows a plethora of strong emission lines, including the Balmer series, He II, He I and metal lines of O III $\lambda$3760 and N III $\lambda\lambda$ 4100, 4640 (blended with He II). The latter lines are consistent with originating from the Bowen fluorescence mechanism. By analyzing literature spectra of past events, we conclude that these lines are common in TDEs. The spectral diversity of optical TDEs is thus larger than previously thought and includes N-rich events besides H- and He-rich events. We study how the spectral lines evolve with time, by means of their width, relative strength, and velocity offsets. The velocity width of the lines starts at $\sim$ 13000 km s$^{-1}$ and decreases with time. The ratio of He II to N III increases with time. The same is true for ASASSN-14li, which has a very similar spectrum to AT 2018dyb but its lines are narrower by a factor of $>$2. We estimate a black hole mass of $M_{\rm BH}$ = $3.3^{+5.0}_{-2.0}\times 10^6$ $M_{\odot}$ by using the $M$-$\sigma$ relation. This is consistent with the black hole mass derived using the MOSFiT transient fitting code. The detection of strong Bowen lines in the optical spectrum is an indirect proof for extreme ultraviolet and (reprocessed) X-ray radiation and favors an accretion origin for the TDE optical luminosity. A model where photons escape after multiple scatterings through a super-Eddington thick disk and its optically thick wind, viewed at an angle close to the disk plane, is consistent with the observations.

Axions are a viable candidate for Cold Dark Matter (CDM) which should generically form minihalos of sub-planetary masses from white-noise isocurvature density fluctuations if the Peccei-Quinn phase transition occurs after inflation. Despite being denser than the larger halos formed out of adiabatic fluctuations from inflation, axion minihalos have surface densities much smaller than the critical value required for gravitational lensing to produce multiple images or high magnification, and hence are practically undetectable as lenses in isolation. However, their lensing effect can be enhanced when superposed near critical curves of other lenses. We propose a method to detect them through photometric monitoring of recently discovered caustic transiting stars behind cluster lenses, under extreme magnification factors $\mu \gtrsim 10^3$--$10^4$ as the lensed stars cross microlensing caustics induced by intracluster stars. For masses of the first gravitationally collapsed minihalos in the range $\sim 10^{-15}$--$10^{-8}\,h^{-1}\,M_\odot$, we show that axion minihalos in galaxy clusters should collectively produce subtle surface density fluctuations of amplitude $\sim 10^{-4}$--$10^{-3}$ on projected length scales of $\sim 10$--$10^4\,$AU, which imprint irregularities in the microlensing light curves of caustic transiting stars. We estimate that, inside a cluster halo and over the age of the Universe, most of these minihalos are likely to avoid dynamic disruption by encounters with stars or other minihalos.

We introduce a new technique to search for gravitational wave events from compact binary mergers that produce a clear signal only in a single gravitational wave detector, and marginal signals in other detectors. Such a situation can arise when the detectors in a network have different sensitivities, or when sources have unfavorable sky locations or orientations. We start with a short list of loud single-detector triggers from regions of parameter space that are empirically unaffected by glitches (after applying signal-quality vetoes). For each of these triggers, we compute evidence for astrophysical origin from the rest of the detector network by coherently combining the likelihoods from all detectors and marginalizing over extrinsic geometric parameters. We report the discovery of two new binary black hole (BBH) mergers in the second observing run of Advanced LIGO and Virgo (O2), in addition to the ones that were reported in Abbott et al. (2018) and Venumadhav et al. (2019). We estimate that the two events have false alarm rates of one in 19 years (60 O2) and one in 11 years (36 O2). One of the events, GW170817A, has primary and secondary masses $m_1^{\rm src} = 56_{-10}^{+16} \, M_\odot$ and $m_2^{\rm src} = 40_{-11}^{+10} \, M_\odot$ in the source frame. The existence of GW170817A should be very informative about the theoretically predicted upper mass gap for stellar mass black holes. Its effective spin parameter is measured to be $\chi_{\rm eff} = 0.5 \pm 0.2$, which is consistent with the tendency of the heavier detected BBH systems to have large and positive effective spin parameters. The other event, GWC170402, will be discussed thoroughly in future work.

This White Paper outlines the importance of addressing the fundamental science theme <<How are charged particles energized in space plasmas>> through a future ESA mission. The White Paper presents five compelling science questions related to particle energization by shocks, reconnection,waves and turbulence, jets and their combinations. Answering these questions requires resolving scale coupling, nonlinearity and nonstationarity, which cannot be done with existing multi-point observations. In situ measurements from a multi-point, multi-scale L-class plasma observatory consisting of at least 7 spacecraft covering fluid, ion and electron scales are needed. The plasma observatory will enable a paradigm shift in our comprehension of particle energization and space plasma physics in general, with very important impact on solar and astrophysical plasmas. It will be the next logical step following Cluster, THEMIS and MMS for the very large and active European space plasmas community. Being one of the cornerstone missions of the future ESA Voyage 2035-2050 science program, it would further strengthen the European scientific and technical leadership in this important field.

Searches for gravitational waves crucially depend on exact signal processing of noisy strain data from gravitational wave detectors, which are known to exhibit significant non-Gaussian behavior. In this paper, we study two distinct non-Gaussian effects in the LIGO/Virgo data which reduce the sensitivity of searches: first, variations in the noise power spectral density (PSD) on timescales of more than a few seconds; and second, loud and abrupt transient `glitches' of terrestrial or instrumental origin. We derive a simple procedure to correct, at first order, the effect of the variation in the PSD on the search background. Given the knowledge of the existence of localized glitches in particular segments of data, we also develop a method to insulate statistical inference from these glitches, so as to cleanly excise them without affecting the search background in neighboring seconds. We show the importance of applying these methods on the publicly available LIGO data, and measure an increase in the detection volume of at least $15\%$ from the PSD-drift correction alone, due to the improved background distribution.

We introduce an algorithm for placing template waveforms for the search of compact binary mergers in gravitational wave interferometer data. We exploit the smooth dependence of the amplitude and unwrapped phase of the frequency-domain waveform on the parameters of the binary. We group waveforms with similar amplitude profiles and perform a singular value decomposition of the phase profiles to obtain an orthonormal basis for the phase functions. The leading basis functions span a lower-dimensional linear space in which the unwrapped phase of any physical waveform is well approximated. The optimal template placement is given by a regular grid in the space of linear coefficients. The algorithm is applicable to any frequency-domain waveform model and detector sensitivity curve. It is computationally efficient and requires little tuning. Applying this method, we construct a set of template banks suitable for the search of aligned-spin binary neutron star, neutron-star--black-hole and binary black hole mergers in LIGO--Virgo data.

Following the pioneering observations with COBE in the early 1990s, studies of the cosmic microwave background (CMB) have focused on temperature and polarization anisotropies. CMB spectral distortions - tiny departures of the CMB energy spectrum from that of a perfect blackbody - provide a second, independent probe of fundamental physics, with a reach deep into the primordial Universe. The theoretical foundation of spectral distortions has seen major advances in recent years, which highlight the immense potential of this emerging field. Spectral distortions probe a fundamental property of the Universe - its thermal history - thereby providing additional insight into processes within the cosmological standard model (CSM) as well as new physics beyond. Spectral distortions are an important tool for understanding inflation and the nature of dark matter. They shed new light on the physics of recombination and reionization, both prominent stages in the evolution of our Universe, and furnish critical information on baryonic feedback processes, in addition to probing primordial correlation functions at scales inaccessible to other tracers. In principle the range of signals is vast: many orders of magnitude of discovery space could be explored by detailed observations of the CMB energy spectrum. Several CSM signals are predicted and provide clear experimental targets, some of which are already observable with present-day technology. Confirmation of these signals would extend the reach of the CSM by orders of magnitude in physical scale as the Universe evolves from the initial stages to its present form. The absence of these signals would pose a huge theoretical challenge, immediately pointing to new physics.

We combine gravitational wave (GW) and electromagnetic (EM) data to perform a Bayesian parameter estimation of the binary neutron star (NS) merger GW170817. The EM likelihood is constructed from a fit to a large number of numerical relativity simulations which we combine with a lower bound on the mass of the remnant's accretion disk inferred from the modeling of the EM light curve. In comparison with previous works, our analysis yields a more precise determination of the tidal deformability of the binary, for which the EM data provide a lower bound, and of the mass ratio of the binary, with the EM data favoring a smaller mass asymmetry. The 90\% credible interval for the areal radius of a $1.4\ M_\odot$ NS is found to be $12.2^{+1.0}_{-0.8} \pm 0.2\ {\rm km}$ (statistical and systematic uncertainties).

Processes of dissipation of energy released during capture and disruption stars by Galactic central Black hole. This process should periodically happen in the Galactic center and affect Galactic center environment as well as Galactic halo. Routine star disruptions by the central supermassive black hole BH provide enough cumulative energy to form the Fermi Bubbles. A single star of super-Eddington eruption provide an subsonic outflow gas whose total energy is about $10^{52}$ erg. The average rate of disruption events is expected to be $10^{-4}-10^{-5}$yr$^{-1}$ that provides the average power of energy release in the GC $\dot{W}\sim 3\times 10^{41}$ erg s$^{-1}$, just needed to support Fermi Bubbles in the halo. In the exponential atmosphere the energy from the GC propagates perpendicular to the Galactic plane and provides symmetric features in the halo.

We discuss the theory of pulsar-timing and astrometry probes of a stochastic gravitational-wave background with a recently developed "total-angular-momentum" (TAM) formalism for cosmological perturbations. We review the formalism, emphasizing in particular the features relevant for this work and describe the observables we consider (i.e. the pulsar redshift and stellar angular displacement). Using the TAM approach, we calculate the angular power spectra for the observables and from them derive angular auto- and cross-correlation functions. We provide the full set of power spectra and correlation functions not only for the standard transverse-traceless propagating degrees of freedom in general relativity, but also for the four additional non-Einsteinian polarizations that may arise in alternative-gravity theories. We discuss how pulsar-timing and astrometry surveys can complement and serve as cross checks to one another and comment on the importance of testing the chirality of the gravitational-wave background as a tool to understand the nature of its sources. A simple rederivation of the power spectra from the plane-wave formalism is provided in an Appendix.

We report a new binary black hole merger in the publicly available LIGO First Observing Run (O1) data release. The event has an inverse false alarm rate of one per six years in the detector-frame chirp-mass range $\mathcal{M}^{\rm det} \in [20,40]M_\odot$ in a new independent analysis pipeline that we developed. Our best estimate of the probability that the event is of astrophysical origin is $P_{\rm astro} \sim 0.71\, .$ The estimated physical parameters of the event indicate that it is the merger of two massive black holes, $\mathcal{M}^{\rm det} = 31^{+2}_{-3}\,M_\odot$ with an effective spin parameter, $\chi_{\rm eff} = 0.81^{+0.15}_{-0.21}$, making this the most highly spinning merger reported to date. It is also among the two highest redshift mergers observed so far. The high aligned spin of the merger supports the hypothesis that merging binary black holes can be created by binary stellar evolution.

In this paper, we report on the construction of a new and independent pipeline for analyzing the public data from the first observing run of advanced LIGO for mergers of compact binary systems. The pipeline incorporates different techniques and makes independent implementation choices in all its stages including the search design, the method to construct template banks, the automatic routines to detect bad data segments ("glitches") and to insulate good data from them, the procedure to account for the non-stationary nature of the detector noise, the signal-quality vetoes at the single-detector level and the methods to combine results from multiple detectors. Our pipeline enabled us to identify a new binary black-hole merger GW151216 in the public LIGO data. This paper serves as a bird's eye view of the pipeline's important stages. Full details and derivations underlying the various stages will appear in accompanying papers.

We examine a caustic-straddling arc at $z=0.9397$ in the field of the galaxy cluster MACS J0416.1-2403 ($z=0.397$) using archival multiband HST images and show that its surface brightness exhibits anomalies that can be explained by a single highly magnified star undergoing microlensing. First, we show that the surface brightness pattern is not perfectly symmetric across the cluster critical curve, which is inconsistent with a locally smooth lens model; the location of the candidate star exhibits the most significant asymmetry. Second, our analysis indicates that the asymmetric feature has $\sim 30\%$ higher flux in the 2012 visits compared to the Frontier Fields program visits in 2014. Moreover, the variable asymmetric feature shows an anomalous color between the F814W and F105W filters in 2014. These anomalies are naturally explained by microlensing induced variability of a caustic-transiting blue supergiant in a star-forming region, with a mean magnification factor around $\mu \sim 200$. We extend this study to a statistical analysis of the whole arc image and find tentative evidence of the increased mismatch of the two images in the proximity of the critical line. Robust detection of one or multiple caustic-transiting stars in this arc will enable detailed follow-up studies that can shed light on the small-scale structure of the dark matter inside the cluster halo.

The tidal disruption of stars by (super-)massive black holes in galactic nuclei has been discussed in theoretical terms for about 30 years but only in the past decade have we been able to detect such events in substantial numbers. Thus, we are now starting to carry out observational tests of models for the disruption. We are also formulating expectations for the inspiral and disruption of white dwarfs by intermediate-mass black holes with masses $< 10^5\;{\rm M}_\odot$. Such events are very rich with information and open a new window to intermediate-mass black holes, thought to live in dwarf galaxies and star clusters. They can inform us of the demographics of intermediate-mass black holes, stellar populations and dynamics in their immediate vicinity, and the physics of accretion of hydrogen-deficient material. The combination of upcoming transient surveys using ground-based, electromagnetic observatories and low-frequency gravitational wave observations is ideal for exploiting tidal disruptions of white dwarfs. The detection rate of gravitational wave signals, optimistically, may reach a few dozen per year in a volume up to $z\approx 0.1$. Gravitational wave observations are particularly useful because they yield the masses of the objects involved and allow determination of the spin of the black hole, affording tests of physical models for black hole formation and growth. They also give us advance warning of the electromagnetic flares by weeks or more. The right computing infrastructure for modern models for the disruption process and event rates will allow us to make the most of the upcoming observing facilities.

The dynamical assembly of binary black holes (BBHs) in dense star clusters (SCs) is one of the most promising pathways for producing observable gravitational wave (GW) sources, however several other formation scenarios likely operate as well. One of the current outstanding questions is how these different pathways may be distinguished apart. In this paper we suggest a new multi-messenger observable that can be used to constrain the formation of BBH mergers originating from SCs: the electromagnetic signal from tidal disruptions (TDs) of stars by BBHs. Such TDs will show variability in their light curve from the orbital motion of the disruptive BBHs, and can therefore be used to map the BBH orbital period distribution, and thereby also the dynamical mechanisms that eventually drive the BBHs to merger. Using an analytical approach including General Relativistic effects, we find that the orbital period distribution of BBHs within globular clusters peaks on timescales of days, which we argue is unique to this assembly pathway. We propose that the search for variable TDs in current and future EM transient surveys might be used to constrain the merger history of BBHs in SCs.

The Wide-Field Infrared Survey Telescope (WFIRST) is expected to launch in the mid-2020s. With its wide-field near-infrared (NIR) camera, it will survey the sky to unprecedented detail. As part of normal operations and as the result of multiple expected dedicated surveys, WFIRST will produce several relatively wide-field (tens of square degrees) deep (limiting magnitude of 28 or fainter) fields. In particular, a planned supernova survey is expected to image 3 deep fields in the LSST footprint roughly every 5 days over 2 years. Stacking all data, this survey will produce, over all WFIRST supernova fields in the LSST footprint, ~12-25 deg^2 and ~5-15 deg^2 regions to depths of ~28 mag and ~29 mag, respectively. We suggest LSST undertake mini-surveys that will match the WFIRST cadence and simultaneously observe the supernova survey fields during the 2-year WFIRST supernova survey, achieving a stacked depth similar to that of the WFIRST data. We also suggest additional observations of these same regions throughout the LSST survey to get deep images earlier, have long-term monitoring in the fields, and produce deeper images overall. These fields will provide a legacy for cosmology, extragalactic, and transient/variable science.

Luminous stars in background galaxies straddling the lensing caustic of a foreground galaxy cluster can be individually detected due to extreme magnification factors of $\sim 10^2$--$10^3$, as recently observed in deep HST images. We propose a direct method to probe the presence of dark matter subhalos in galaxy clusters by measuring the astrometric perturbation they induce on the image positions of magnified stars or bright clumps: lensing by subhalos breaks the symmetry of a smooth critical curve, traced by the midpoints of close image pairs. For the giant arc at $z = 0.725$ behind the lensing cluster Abell 370 at $z = 0.375$, a promising target for detecting image pairs of stars, we find that subhalos of masses in the range $10^6$--$10^8\,M_\odot$ with the abundance predicted in the cold dark matter theory should typically imprint astrometric distortions at the level of $20$--$80\,{\rm mas}$. We estimate that $\sim 10\,$hr integrations with JWST at $\sim 1$--$3\,\mu{\rm m}$ may uncover several magnified stars whose image doublets will reveal the subhalo-induced structures of the critical curve. This method can probe a dynamic range in the subhalo to cluster halo mass ratio $m/M \sim 10^{-7}$--$10^{-9}$, thereby placing new constraints on the nature of dark matter.

Gravitational waves emitted from compact binary coalescence can be subject to wave diffraction if they are gravitationally lensed by an intervening mass clump whose Schwarzschild timescale matches the wave period. Waves in the ground-based frequency band $f\sim 10$--$10^3\,$Hz are sensitive to clumps with masses $M_E \sim 10^2$--$10^3\,M_\odot$ enclosed within the impact parameter. These can be the central parts of low mass $M_L \sim 10^3$--$10^6\,M_\odot$ dark matter halos, which are predicted in Cold Dark Matter scenarios but are challenging to observe. Neglecting finely-tuned impact parameters, we focus on lenses aligned generally on the Einstein scale for which multiple lensed images may not form in the case of an extended lens. In this case, diffraction induces amplitude and phase modulations whose sizes $\sim 10\%$--$20\%$ are small enough so that standard matched filtering with unlensed waveforms do not degrade, but are still detectable for events with high signal-to-noise ratio. We develop and test an agnostic detection method based on dynamic programming, which does not require a detailed model of the lensed waveforms. For pseudo-Jaffe lenses aligned up to the Einstein radius, we demonstrate that a pair of fully upgraded aLIGO/Virgo detectors can extract diffraction imprints from binary black hole mergers out to $z_s \sim 0.2$--$0.3$. The prospect will improve dramatically for a third-generation detector for which binary black hole mergers out to $z_s \sim 2$--$4$ will all become valuable sources.

Relative binning is a new method for fast and accurate evaluation of the likelihood of gravitational wave strain data. This technique can be used to produce reliable posterior distributions for compact object mergers with very moderate computational resources. We use a fast likelihood evaluation code based on this technique to estimate the parameters of the double neutron-star merger event GW170817 using publicly available LIGO data. We obtain statistically similar posteriors using either Markov-chain Monte-Carlo or nested sampling. The results do not favor non-zero aligned spins at a statistically significant level. There is no significant sign of non-zero tidal deformability (as quantified by the Bayesian evidence), whether or not high-spin or low-spin priors are adopted. Our posterior samples are publicly available, and we also provide a tutorial Python code to implement fast likelihood evaluation using the relative binning method.

We present a method to accelerate the evaluation of the likelihood in gravitational wave parameter estimation. Parameter estimation codes compute likelihoods of similar waveforms, whose phases and amplitudes differ smoothly with frequency. We exploit this by precomputing frequency-binned overlaps of the best-fit waveform with the data. We show how these summary data can be used to approximate the likelihood of any waveform that is sufficiently probable within the required accuracy. We demonstrate that $\simeq 60$ bins suffice to accurately compute likelihoods for strain data at a sampling rate of $4096\,$Hz and duration of $T=2048\,$s around the binary neutron star merger GW170817. Relative binning speeds up parameter estimation for frequency domain waveform models by a factor of $\sim 10^4$ compared to naive matched filtering and $\sim 10$ compared to reduced order quadrature.

We revisit the 21 cm power spectrum from the epoch of cosmic dawn in light of the recent EDGES detection of the 21 cm global signal at frequencies corresponding to $z\sim20$. The shape of the signal suggests that the spin temperature of neutral hydrogen was coupled to the kinetic temperature of the gas relatively rapidly ($19\lesssim z \lesssim 21$). We therefore consider models in which the UV photons were dominantly produced in the rarest and most massive halos ($M\gtrsim 10^9M_\odot$), since their abundance grows fast enough at those redshifts to account for this feature of the signal. We show that these models predict large power spectrum amplitudes during the inhomogeneous coupling, and then inhomogeneous heating by CMB and Lyman-$\alpha$ photons due to the large shot noise associated with the rare sources. The power spectrum is enhanced by more than an order of magnitude compared to previous models which did not include the shot noise contribution, making it a promising target for upcoming radio interferometers that aim to detect high-redshift 21 cm fluctuations.

The intergalactic medium is expected to be at its coldest point before the formation of the first stars in the universe. Motivated by recent results from the EDGES experiment, we revisit the standard calculation of the kinetic temperature of the neutral gas through this period. When the first ultraviolet (UV) sources turn on, photons redshift into the Lyman lines of neutral hydrogen and repeatedly scatter within the Lyman-$\alpha$ line. They heat the gas via atomic recoils, and, through the Wouthuysen-Field effect, set the spin temperature of the 21-cm hyperfine (spin-flip) line of atomic hydrogen in competition with the resonant cosmic microwave background (CMB) photons. We show that the Lyman-$\alpha$ photons also mediate energy transfer between the CMB photons and the thermal motions of the hydrogen atoms. In the absence of X-ray heating, this new mechanism is the major correction to the temperature of the adiabatically cooling gas ($\sim 10 \%$ at $z=17$), and is several times the size of the heating rate found in previous calculations. We also find that the effect is more dramatic in non-standard scenarios that either enhance the radio background above the CMB or invoke new physics to cool the gas in order to explain the EDGES results. The coupling with the radio background can reduce the depth of the 21-cm absorption feature by almost a factor of two relative to the case with no sources of heating, and prevent the feature from developing a flattened bottom. As an inevitable consequence of the UV background that generates the absorption feature, this heating should be accounted for in any theoretical prediction.

In the past few years wide-field optical and UV transient surveys as well as X-ray telescopes have allowed us to identify a few dozen candidate tidal disruption events (TDEs). While in theory the physical processes in TDEs are expected to be ubiquitous, a few distinct classes of TDEs have been observed. Some TDEs radiate mainly in NUV/optical while others produce prominent X-rays. Moreover, relativistic jets have been observed in only a handful of TDEs. This diversity might be related to the details of the super-Eddington accretion and emission physics relevant to TDE disks. In this Letter, we utilize novel three-dimensional general relativistic radiation magnetohydrodynamics simulations to study the super-Eddington compact disk phase expected in TDEs. Consistent with previous studies, geometrically thick disks, wide-angle optically-thick fast outflows and relativistic jets are produced. The outflow density and velocity depend sensitively on the inclination angle, and hence so does the reprocessing of emission produced from the inner disk. We then use Monte-Carlo radiative transfer to calculate the reprocessed spectra and find that that the observed ratio of optical to X-ray fluxes increases with increasing inclination angle. This naturally leads to a unified model for different classes of TDEs in which the spectral properties of the TDE depend mainly on the viewing-angle of the observer with respect to the orientation of the disk.

Recent observations of lensed galaxies at cosmological distances have detected individual stars that are extremely magnified when crossing the caustics of lensing clusters. In idealized cluster lenses with smooth mass distributions, two images of a star of radius $R$ approaching a caustic brighten as $t^{-1/2}$ and reach a peak magnification $\sim 10^{6}\, (10\, R_{\odot}/R)^{1/2}$ before merging on the critical curve. We show that a mass fraction ($\kappa_\star \gtrsim \, 10^{-4.5}$) in microlenses inevitably disrupts the smooth caustic into a network of corrugated microcaustics, and produces light curves with numerous peaks. Using analytical calculations and numerical simulations, we derive the characteristic width of the network, caustic-crossing frequencies, and peak magnifications. For the lens parameters of a recent detection and a population of intracluster stars with $\kappa_\star \sim 0.01$, we find a source-plane width of $\sim 20 \, {\rm pc}$ for the caustic network, which spans $0.2 \, {\rm arcsec}$ on the image plane. A source star takes $\sim 2\times 10^4$ years to cross this width, with a total of $\sim 6 \times 10^4$ crossings, each one lasting for $\sim 5\,{\rm hr}\,(R/10\,R_\odot)$ with typical peak magnifications of $\sim 10^{4} \left( R/ 10\,R_\odot \right)^{-1/2}$. The exquisite sensitivity of caustic-crossing events to the granularity of the lens-mass distribution makes them ideal probes of dark matter components, such as compact halo objects and ultralight axion dark matter.

At only 90 Mpc, ASASSN-14li is one of the nearest tidal disruption event (TDE) ever discovered, and because of this, it has been observed by several observatories at many wavelengths. In this paper, we present new results on archival XMM-Newton observations, three of which were taken at early times (within 40 days of the discovery), and three of which were taken at late times, about one year after the peak. We find that, at early times, in addition to the ~1e5 K blackbody component that dominates the X-ray band, there is evidence for a broad, P Cygni-like absorption feature at around 0.7 keV in all XMM-Newton instruments (CCD detectors and grating spectrometers), and that this feature disappears (or at least diminishes) in the late-time observations. We perform photoionization modelling with XSTAR and interpret this absorption feature as blueshifted OVIII, from an ionized outflow with a velocity of 0.2c. As the TDE transitions from high to low accretion rate, the outflow turns off, thus explaining the absence of the absorption feature during the late-time observations.

Given the possible repetitive nature of fast radio bursts (FRBs), their cosmological origin, and their high occurrence, detection of strongly lensed sources due to intervening galaxy lenses is possible with forthcoming radio surveys. We show that if multiple images of a repeating source are resolved with VLBI, using a method independent of lens modeling, accurate timing could reveal non-uniform motion, either physical or apparent, of the emission spot. This can probe the physical nature of FRBs and their surrounding environments, constraining scenarios including orbital motion around a stellar companion if FRBs require a compact star in a special system, and jet-medium interactions for which the location of the emission spot may randomly vary. The high timing precision possible for FRBs ($\sim {\rm ms}$) compared to the typical time delays between images in galaxy lensing ($\gtrsim 10\, {\rm days}$) enables the measurement of tiny fractional changes in the delays ($\sim 10^{-9}$), and hence the detection of time-delay variations induced by relative motions between the source, the lens, and the Earth. We show that uniform cosmic peculiar velocities only cause the delay time to drift linearly, and that the effect from the Earth's orbital motion can be accurately subtracted, thus enabling a search for non-trivial source motion. For a timing accuracy of $\sim 1\,$ms and a repetition rate (of detected bursts) $\sim 0.05$ per day of a single FRB source, non-uniform displacement $\gtrsim 0.1 - 1\,$AU of the emission spot perpendicular to the line of sight is detectable if repetitions are seen over a period of hundreds of days.

Powerful jets and outflows generated in tidal disruption events (TDEs) around supermassive black holes have been suggested to be possible sites producing high-energy neutrinos, but it is unclear whether such environment can provide the bulk of the neutrinos detected by the IceCube Observatory. In this work, by considering realistic limits on the non-thermal emission power of a TDE jet and the birth rate of the TDEs with jets pointing towards us, we show that it is hard to use the jetted TDE population to explain the flux and the isotropic arrival directions of the observed TeV-PeV neutrinos. Therefore, TDEs cannot be the dominant sources, unless those without aligned jets can produce wide-angle emission of neutrino particles. Supposing that is the case, we list a few recent jetted and non-jetted TDEs that have the best chance to be detected by IceCube, based on their energetics, distances, and directions. A spatial and temporal association of these predicted events with the IceCube data should provide a decisive test on the full TDE population as origins of the IceCube neutrinos.

Strong lensing by intervening galaxies can produce multiple images of gravitational waves from sources at cosmological distances. These images acquire additional phase-shifts as the over-focused wavefront passes through itself along the line of sight. Time domain waveforms of Type-II images (associated with saddle points of the time delay) exhibit a non-trivial distortion from the unlensed waveforms. This phenomenon is in addition to the usual frequency-independent magnification, and happens even in the geometric limit where the wavelength is much shorter than the deflector's gravitational length scale. Similarly, Type-III images preserve the original waveform's shape but exhibit a sign flip. We show that for non-precessing binaries undergoing circular inspiral and merger, these distortions are equivalent to rotating the line of sight about the normal to the orbital plane by $45^\circ$ (Type II) and $90^\circ$ (Type III). This effect will enable us to distinguish between the different topological types among a set of multiple images, and give us valuable insight into the lens model. Furthermore, we show that for eccentric binaries, the waveform of a Type-II image is distorted in a manner that is inequivalent to a change of the source's orbital parameters.

The recent detection of gravitational waves indicates that stellar-mass black hole binaries are likely to be a key population of sources for forthcoming observations. With future upgrades, ground-based detectors could detect merging black hole binaries out to cosmological distances. Gravitational wave bursts from high redshifts ($z \gtrsim 1$) can be strongly magnified by gravitational lensing due to intervening galaxies along the line of sight. In the absence of electromagnetic counterparts, the mergers' intrinsic mass scale and redshift are degenerate with the unknown magnification factor $\mu$. Hence, strongly magnified low-mass mergers from high redshifts appear as higher-mass mergers from lower redshifts. We assess the impact of this degeneracy on the mass-redshift distribution of observable events for generic models of binary black hole formation from normal stellar evolution, Pop III star remnants, or a primordial black hole population. We find that strong magnification ($\mu \gtrsim 3$) generally creates a heavy tail of apparently massive mergers in the event distribution from a given detector. For LIGO and its future upgrades, this tail may dominate the population of intrinsically massive, but unlensed mergers in binary black hole formation models involving normal stellar evolution or primordial black holes. Modeling the statistics of lensing magnification can help account for this magnification bias when testing astrophysical scenarios of black hole binary formation and evolution.

The possible Fermi detection of an electromagnetic counterpart to the double black hole merger GW150914 has inspired many theoretical models, some of which propose that the holes spiraled together inside a massive star. However, we show that the heat produced by the dynamical friction on such black hole orbits can exceed the stellar binding energy by a large factor, which means that this heat could destroy the star and thus make it difficult for enough gas to be near the holes at merger to produce detectable photons. These considerations must be taken into account when models are proposed for electromagnetic counterparts to the coalescence of two stellar-mass black holes. We find that only when the two black holes form very close to the center can the star avoid destruction. In that case, dynamical friction can make the black holes coalesce faster than they would in vacuum, which leads to a modification of the gravitational waveform that is potentially observable by advanced LIGO.

In the next few years, intensity-mapping surveys that target lines such as CO, Ly$\alpha$, and CII stand to provide powerful probes of high-redshift astrophysics. However, these line emissions are highly non-Gaussian, and so the typical power-spectrum methods used to study these maps will leave out a significant amount of information. We propose a new statistic, the probability distribution of voxel intensities, which can access this extra information. Using a model of a CO intensity map at $z\sim3$ as an example, we demonstrate that this voxel intensity distribution (VID) provides substantial constraining power beyond what is obtainable from the power spectrum alone. We find that a future survey similar to the planned COMAP Full experiment could constrain the CO luminosity function to order $\sim10\%$. We also explore the effects of contamination from continuum emission, interloper lines, and gravitational lensing on our constraints and find that the VID statistic retains significant constraining power even in pessimistic scenarios.

The possibility that part of the dark matter is made of massive compact halo objects (MACHOs) remains poorly constrained over a wide range of masses, and especially in the $20-100\, M_\odot$ window. We show that strong gravitational lensing of extragalactic fast radio bursts (FRBs) by MACHOs of masses larger than $\sim20\,M_\odot$ would result in repeated FRBs with an observable time delay. Strong lensing of an FRB by a lens of mass $M_L$ induces two images, separated by a typical time delay $\sim$ few $\times(M_L/30\, M_\odot)$ milliseconds. Considering the expected FRB detection rate by upcoming experiments, such as CHIME, of $10^4$ FRBs per year, we should observe from tens to hundreds of repeated bursts yearly, if MACHOs in this window make up all the dark matter. A null search for echoes with just $10^4$ FRBs, would constrain the fraction $f_{\rm DM}$ of dark matter in MACHOs to $f_{\rm DM}\lesssim 0.08$ for $M_L\gtrsim 20\,M_\odot$.

Our current understanding of the curved space-time around supermassive black holes is based on actively accreting black holes, which make up only ten per cent or less of the overall population. X-ray observations of that small fraction reveal strong gravitational redshifts that indicate that many of these black holes are rapidly rotating; however, selection biases suggest that these results are not necessarily reflective of the majority of black holes in the Universe. Tidal disruption events, where a star orbiting an otherwise dormant black hole gets tidally shredded and accreted onto the black hole, can provide a short, unbiased glimpse at the space-time around the other ninety per cent of black holes. Observations of tidal disruptions have hitherto revealed the formation of an accretion disk and the onset of an accretion-powered jet, but have failed to reveal emission from the inner accretion flow, which enables the measurement of black hole spin. Here we report observations of reverberation arising from gravitationally redshifted iron K photons reflected off the inner accretion flow in the tidal disruption event Swift J1644+57. From the reverberation timescale, we estimate the mass of the black hole to be a few million solar masses, suggesting an accretion rate of 100 times the Eddington limit or more. The detection of reverberation from the relativistic depths of this rare super-Eddington event demonstrates that the X-rays do not arise from the relativistically moving regions of a jet, as previously thought.

The auto-correlation between two members of a galaxy population is symmetric under the interchange of the two galaxies being correlated. The cross-correlation between two different types of galaxies, separated by a vector $\bf{r}$, is not necessarily the same as that for a pair separated by $-\bf{r}$. Local anisotropies in the two-point cross-correlation function may thus indicate a specific direction which when mapped as a function of position trace out a vector field. This vector field can then be decomposed into longitudinal and transverse components, and those transverse components written as positive- and negative-helicity components. A locally asymmetric cross-correlation of the longitudinal type arises naturally in halo clustering, even with Gaussian initial conditions, and could be enhanced with local-type non-Gaussianity. Early-Universe scenarios that introduce a vector field may also give rise to such effects. These antisymmetric cross-correlations also provide a new possibility to seek a preferred cosmic direction correlated with the hemispherical power asymmetry in the cosmic microwave background and to seek a preferred location associated with the CMB cold spot. New ways to seek cosmic parity breaking are also possible.

In the standard inflationary scenario, primordial perturbations are adiabatic. The amplitudes of most types of isocurvature perturbations are generally constrained by current data to be small. If, however, there is a baryon-density perturbation that is compensated by a dark-matter perturbation in such a way that the total matter density is unperturbed, then this compensated isocurvature perturbation (CIP) has no observable consequence in the cosmic microwave background (CMB) at linear order in the CIP amplitude. Here we search for the effects of CIPs on CMB power spectra to quadratic order in the CIP amplitude. An analysis of the Planck temperature data leads to an upper bound $\Delta_{\rm rms}^2 \leq 7.1\times 10^{-3}$, at the 68\% confidence level, to the variance $\Delta_{\rm rms}^2$ of the CIP amplitude. This is then strengthened to $\Delta_{\rm rms}^2\leq 5.0\times 10^{-3}$ if Planck small-angle polarization data are included. A cosmic-variance-limited CMB experiment could improve the $1\sigma$ sensitivity to CIPs to $\Delta^2_{\rm rms} \lesssim 9\times 10^{-4}$. It is also found that adding CIPs to the standard $\Lambda$CDM model can improve the fit of the observed smoothing of CMB acoustic peaks just as much as adding a non-standard lensing amplitude.

Massive black hole binaries (MBHBs) are a natural byproduct of galaxy mergers. Previous studies have shown that flares from stellar tidal disruption events (TDEs) are modified by the presence of a secondary perturber, causing interruptions in the light curve. We study the dynamics of TDE debris in the presence of a milliparsec-separated MBHB by integrating ballistic particle orbits in the time-varying potential of the binary. We find that gaps in the light curve appear when material misses the accretion radius on its first return to pericentre. Subsequent recurrences can be decomposed into "continuous" and "delayed" components, which exhibit different behaviour. We find that this potential can substantially alter the locations of stream self-intersections. When debris is confined to the plane, we find that close encounters with the secondary BH leave noticeable signatures on the fallback rate and can result in significant accretion onto the secondary BH. Tight, equal-mass MBHBs accrete equally, periodically trading the infalling stream.

One of the puzzles associated with tidal disruption event candidates (TDEs) is that there is a dichotomy between the color temperatures of ${\rm few}\times 10^4$~K for TDEs discovered with optical and UV telescopes, and the color temperatures of ${\rm few}\times 10^5 - 10^6$~K for TDEs discovered with X-ray satellites. Here we propose that high-temperature TDEs are produced when the tidal debris of a disrupted star self-intersects relatively close to the supermassive black hole, in contrast to the more distant self-intersection that leads to lower color temperatures. In particular, we note from simple ballistic considerations that greater apsidal precession in an orbit is the key to closer self-intersection. Thus larger values of $\beta$, the ratio of the tidal radius to the pericenter distance of the initial orbit, are more likely to lead to higher temperatures of more compact disks which are super-Eddington and geometrically and optically thick. For a given star and $\beta$, apsidal precession also increases for larger black hole masses, but larger black hole masses imply a lower temperature at the Eddington luminosity. Thus the expected dependence of the temperature on the mass of the black hole is non-monotonic. We find that in order to produce a soft X-ray temperature TDE, a deep plunging stellar orbit with $\beta> 3$ is needed and a black hole mass of $\lesssim 5\times 10^6 M_\odot$ is favored. Although observations of TDEs are comparatively scarce and are likely dominated by selection effects, it is encouraging that both expectations are consistent with current data.

The radiative efficiency of super-Eddington accreting black holes (BHs) is explored for magnetically-arrested disks (MADs), where magnetic flux builds-up to saturation near the BH. Our three-dimensional general relativistic radiation magnetohydrodynamic (GRRMHD) simulation of a spinning BH (spin $a/M=0.8$) accreting at $\sim 50$ times Eddington shows a total efficiency $\sim 50\%$ when time-averaged and total efficiency $\gtrsim 100\%$ in moments. Magnetic compression by the magnetic flux near the rotating BH leads to a thin disk, whose radiation escapes via advection by a magnetized wind and via transport through a low-density channel created by a Blandford-Znajek (BZ) jet. The BZ efficiency is sub-optimal due to inertial loading of field lines by optically thick radiation, leading to BZ efficiency $\sim 40\%$ on the horizon and BZ efficiency $\sim 5\%$ by $r\sim 400r_g$ (gravitational radii) via absorption by the wind. Importantly, radiation escapes at $r\sim 400r_g$ with efficiency $\eta\approx 15\%$ (luminosity $L\sim 50L_{\rm Edd}$), similar to $\eta\approx 12\%$ for a Novikov-Thorne thin disk and beyond $\eta\lesssim 1\%$ seen in prior GRRMHD simulations or slim disk theory. Our simulations show how BH spin, magnetic field, and jet mass-loading affect the radiative and jet efficiencies of super-Eddington accretion.

(abridged version) The separate universe conjecture states that in General Relativity a density perturbation behaves locally (i.e. on scales much smaller than the wavelength of the mode) as a separate universe with different background density and curvature. We prove this conjecture for a spherical compensated tophat density perturbation of arbitrary amplitude and radius in $\Lambda$CDM. We then use Conformal Fermi Coordinates to generalize this result to scalar perturbations of arbitrary configuration and scale. In this case, the separate universe conjecture holds for the isotropic part of the perturbations. The anisotropic part on the other hand is exactly captured by a tidal field in the Newtonian form. We show that the separate universe picture is restricted to scales larger than the sound horizons of all fluid components. We then derive an expression for the locally measured matter bispectrum induced by a long-wavelength mode of arbitrary wavelength. We show that nonlinear gravitational dynamics does not generate observable contributions that scale like local-type non-Gaussianity $f_{\rm NL}^{\rm loc}$, and hence does not contribute to a scale-dependent galaxy bias $\Delta b \propto k^{-2}$ on large scales; rather, the locally measurable long-short mode coupling assumes a form essentially identical to subhorizon perturbation theory results, once the long-mode density perturbation is replaced by the synchronous-comoving gauge density perturbation. Apparent $f_{\rm NL}^{\rm loc}$-type contributions arise through projection effects on photon propagation, which depend on the specific large-scale structure tracer and observable considered, and are in principle distinguishable from the local mode coupling induced by gravity. We conclude that any observation of $f_{\rm NL}^{\rm loc}$ beyond these projection effects signals a departure from standard single-clock inflation.

Fermi Normal Coordinates (FNC) are a useful frame for isolating the locally observable, physical effects of a long-wavelength spacetime perturbation. Their cosmological application, however, is hampered by the fact that they are only valid on scales much smaller than the horizon. We introduce a generalization that we call Conformal Fermi Coordinates (CFC). CFC preserve all the advantages of FNC, but in addition are valid outside the horizon. They allow us to calculate the coupling of long- and short-wavelength modes on all scales larger than the sound horizon of the cosmological fluid, starting from the epoch of inflation until today, by removing the complications of the second order Einstein equations to a large extent, and eliminating all gauge ambiguities. As an application, we present a calculation of the effect of long-wavelength tensor modes on small scale density fluctuations. We recover previous results, but clarify the physical content of the individual contributions in terms of locally measurable effects and "projection" terms.

We investigate the possibility to detect the integrated Sachs-Wolfe (ISW) effect by cross-correlating 21-cm surveys at high redshifts with galaxies, in a way similar to the usual CMB-galaxy cross-correlation. The high-redshift 21-cm signal is dominated by CMB photons that travel freely without interacting with the intervening matter, and hence its late-time ISW signature should correlate extremely well with that of the CMB at its peak frequencies. Using the 21-cm temperature brightness instead of the CMB would thus be a further check of the detection of the ISW effect, measured by different instruments at different frequencies and suffering from different systematics. We also study the ISW effect on the photons that are scattered by HI clouds. We show that a detection of the unscattered photons is achievable with planned radio arrays, while one using scattered photons will require advanced radio interferometers, either an extended version of the planned Square Kilometre Array or futuristic experiments such as a lunar radio array.

Spectral distortions of the cosmic microwave background (CMB) may become a powerful probe of primordial perturbations at small scales. Existing studies of spectral distortions focus almost exclusively on primordial scalar metric perturbations. Similarly, vector and tensor perturbations should source CMB spectral distortions. In this paper, we give general expressions for the effective heating rate caused by these types of perturbations, including previously neglected contributions from polarization states and higher multipoles. We then focus our discussion on the dissipation of tensors, showing that for nearly scale invariant tensor power spectra, the overall distortion is some six orders of magnitudes smaller than from the damping of adiabatic scalar modes. We find simple analytic expressions describing the effective heating rate from tensors using a quasi-tight coupling approximation. In contrast to adiabatic modes, tensors cause heating without additional photon diffusion and thus over a wider range of scales, as recently pointed out by Ota et. al 2014. Our results are in broad agreement with their conclusions, but we find that small-scale modes beyond k< 2x10^4 Mpc^{-1} cannot be neglected, leading to a larger distortion, especially for very blue tensor power spectra. At small scales, also the effect of neutrino damping on the tensor amplitude needs to be included.

Future galaxy surveys will measure baryon acoustic oscillations (BAOs) with high significance, and a complete understanding of the anisotropies of BAOs in redshift space will be important to exploit the cosmological information in BAOs. Here we describe the anisotropies that arise in the redshift-space galaxy two-point correlation function (2PCF) and elucidate the origin of features that arise in the dependence of the BAOs on the angle between the orientation of the galaxy pair and the line of sight. We do so with a derivation of the configuration-space 2PCF using streaming model. We find that, contrary to common belief, the locations of BAO peaks in the redshift-space 2PCF are anisotropic even in the linear theory. Anisotropies in BAO depend strongly on the method of extracting the peak, showing maximum 3 % angular variation. We also find that extracting the BAO peak of $r^2\xi(r,\mu)$ significantly reduces the anisotropy to sub-percent level angular variation. When subtracting the tilt due to the broadband behavior of the 2PCF, the BAO bump is enhanced along the line of sight because of local infall velocities toward the BAO bump. Precise measurement of the angular dependence of the redshift-space 2PCF will allow new geometrical tests of dark energy beyond the BAO.

Evidence from the BICEP2 experiment for a significant gravitational-wave background has focussed attention on inflaton potentials $V(\phi) \propto \phi^\alpha$ with $\alpha=2$ ("chaotic" or "$m^2\phi^2$" inflation) or with smaller values of $\alpha$, as may arise in axion-monodromy models. Here we show that reheating considerations may provide additional constraints to these models. The reheating phase preceding the radiation era is modeled by an effective equation-of-state parameter $w_{\rm re}$. The canonical reheating scenario is then described by $w_{\rm re}=0$. The simplest $\alpha=2$ models are consistent with $w_{\rm re} = 0$ for values of $n_s$ well within the current $1\sigma$ range. Models with $\alpha=1$ or $\alpha=2/3$ require a more exotic reheating phase, with $-1/3<w_{\rm re}<0$, unless $n_s$ falls above the current $1\sigma$ range. Likewise, models with $\alpha=4$ require a physically implausible $w_{\rm re}>1/3$, unless $n_s$ is close to the lower limit of the $2\sigma$ range. For $m^2\phi^2$ inflation and canonical reheating as a benchmark, we derive a relation $\log_{10}\left(T_{\rm re}/10^6\,{\rm GeV} \right) \simeq 2000\,(n_s-0.96)$ between the reheat temperature $T_{\rm re}$ and the scalar spectral index $n_s$. Thus, if $n_s$ is close to its central value, then $T_{\rm re}\lesssim 10^6$~GeV, just above the electroweak scale. If the reheat temperature is higher, as many theorists may prefer, then the scalar spectral index should be closer to $n_s\simeq0.965$ (at the pivot scale $k=0.05\,{\rm Mpc}^{-1}$), near the upper limit of the $1\sigma$ error range. Improved precision in the measurement of $n_s$ should allow $m^2\phi^2$, axion-monodromy, and $\phi^4$ models to be distinguished, even without precise measurement of $r$, and to test the $m^2\phi^2$ expectation of $n_s\simeq0.965$.

Aberration kernels describe how harmonic-space multipole coefficients of cosmic microwave background (CMB) observables transform under Lorentz boosts of the reference frame. For spin-weighted CMB observables, transforming like the CMB temperature (i.e. Doppler weight $d = 1$), we show that the aberration kernels are the matrix elements of a unitary boost operator in harmonic space. Algebraic properties of the rotation and boost generators then give simple, exact recursion relations that allow us to raise or lower the multipole quantum numbers $\ell$ and $m$, and the spin weight $s$. Further recursion relations express kernels of other Doppler weights $d \neq 1$ in terms of the $d = 1$ kernels. From those we show that on the full sky, to all orders in $\beta=v/c$, $E$- and $B$-mode polarization observables do not mix under aberration if and only if $d = 1$. The new relations, fully non-linear in the boost velocity $\beta$, form the basis of a practical recursive algorithm to accurately compute any aberration kernel. In addition, we develop a second, fast algorithm in which aberration kernels are obtained using a set of ordinary differential equations. This system can also be approximately solved at small scales, providing simple asymptotic formulae for the aberration kernels. The results of this work will be useful for further studying the effect of aberration on future CMB temperature and polarization analysis, and might provide a basis for relativistic radiative transfer schemes.

The recent indications for a tensor-to-scalar ratio $ r \simeq 0.2$ from BICEP2 measurements of the cosmic microwave background (CMB) B-mode polarization present some tension with upper limits $r \lesssim 0.1$ from measurements of CMB temperature fluctuations. Here we point out that tensor perturbations can induce B modes in the CMB polarization without inducing any temperature fluctuations nor E-mode polarization whatsoever, but only, at the expense of violating the Copernican principle. We present this mathematical possibility as a new ingredient for the model-builder's toolkit in case the tension between B modes and temperature fluctuations cannot be resolved with more conventional ideas.

For the simplest inflation models, the recent detection of a large primordial B-mode polarization signal by the BICEP2 experiment indicates a slight tension with the upper limit on the tensor-to-scalar ratio, $r$, from the Planck satellite. Here, we discuss spatially varying $r$ as a possible explanation for this discrepancy. This idea seems attractive since it may also explain part of the hemispherical temperature power asymmetry seen by WMAP and Planck at large angular scales. If these two aspects are indeed connected, the model suggests that in the Northern hemisphere $r$ should be much smaller, a hypothesis that could be confirmed with future B-mode experiments, providing a test for the stationarity of primordial tensor contributions across the sky. The BICEP2 measurement furthermore rules out that a simple dipolar modulation of $r$ alone can be responsible for the full hemispherical power asymmetry.

When a cosmic microwave background (CMB) photon travels from the surface of last scatter through spacetime metric perturbations, the polarization vector may rotate about its direction of propagation. This gravitational rotation is distinct from, and occurs in addition to, the lensing deflection of the photon trajectory. This rotation can be sourced by linear vector or tensor metric perturbations and is fully coherent with the curl deflection field. Therefore, lensing corrections to the CMB polarization power spectra as well as the temperature-polarization cross-correlations due to non-scalar perturbations are modified. The rotation does not affect lensing by linear scalar perturbations, but needs to be included when calculations go to higher orders. We present complete results for weak lensing of the full-sky CMB power spectra by general linear metric perturbations, taking into account both deflection of the photon trajectory and rotation of the polarization. For the case of lensing by gravitational waves, we show that the B modes induced by the rotation largely cancel those induced by the curl component of deflection.

Future high-resolution, high-sensitivity Sunyaev-Zeldovich (SZ) observations of individual clusters will allow us to study the dynamical state of the intra cluster medium (ICM). To exploit the full potential of this new observational window, it is crucial to understand the different contributions to the total SZ signal. Here, we investigate the signature caused by multiple scatterings at lowest order of the electron temperature. Previous analytic discussions of this problem used the isotropic scattering approximation (ISA), which even for the simplest cluster geometries is rather rough. We take a step forward and consistently treat the anisotropy of the ambient radiation field caused by the first scattering. We show that the multiple scattering SZ signal directly probes line of sight anisotropies of the ICM, thereby delivering a new set of SZ observables which could be used for 3D cluster-profile reconstruction. The multiple scattering signal should furthermore correlate spatially with the cluster's X-ray and SZ polarization signals, an effect that could allow enhancing the detectability of this contribution.

We study the multiple scattering Sunyaev-Zeldovich (SZ) signature, extending our previous analysis to high-temperature clusters. We consistently treat the anisotropy of the ambient radiation field caused by the first scattering and also consider lowest order kinematic terms. We show that due to temperature corrections monopole through octupole anisotropy of the singly scattered SZ signal attain different spectra in the second scattering. The difference becomes more pronounced at high temperature, and thus could be used to constrain individual line of sight moments of the electron density and temperature profiles. While very challenging from the observational point of view, this further extends the list of possible SZ observables that will be important for 3D cluster-profile reconstruction, possibly helping to break geometric degeneracies caused by projection effects. We also briefly discuss the scattering of primordial CMB anisotropies by SZ clusters.

Our motion relative to the cosmic-microwave-background (CMB) rest frame deflects light rays giving rise to shifts as large as L -> L(1+beta), where beta=0.00123 is our velocity (in units of the speed of light) on measurements CMB fluctuations. Here we present a novel harmonic-space approach to this CMB aberration that improves upon prior work by allowing us to (i) go to higher orders in beta, thus extending the validity of the analysis to measurements at L > 1/beta ~ 800; and (ii) treat the effects of window functions and pixelization in a more accurate and computationally efficient manner. We calculate precisely the magnitude of the systematic bias in the power spectrum inferred from the partial sky, and show that aberration shifts the multipole moment by Delta L/L ~ beta<cos(theta)>, with <cos(theta)> averaged over the survey footprint. Such a shift, if ignored, would bias the measurement of the sound-horizon size theta_* at the 0.01%-level, which is comparable to the measurement uncertainties of Planck. The bias can then propagate into cosmological parameters such as the angular-diameter distance, Hubble parameter and dark-energy equation of state. We study the effect of aberration for current Planck, South Pole Telescope (SPT) and Atacama Cosmology Telescope (ACT) data and show that the bias cannot be neglected. We suggest that the small tension between Planck, ACT, and SPT may be due partially to aberration. An Appendix shows how the near constancy of the full-sky power spectrum under aberration follows from unitarity of the aberration kernel.

We have carried out general relativistic particle simulations of stars tidally disrupted by massive black holes. When a star is disrupted in a bound orbit with moderate eccentricity instead of a parabolic orbit, the temporal behavior of the resulting stellar debris changes qualitatively. The debris is initially all bound, returning to pericenter in a short time ~ the original stellar orbital timescale. The resulting fallback rate can thus be much higher than the Eddington rate. Furthermore if the star is disrupted close to the hole, in a regime where general relativity is important, the stellar and debris orbits display general relativistic precession. Apsidal precession can make the debris stream cross itself after several orbits, likely leading to fast debris energy dissipation. If the star is disrupted in an inclined orbit around a spinning hole, nodal precession reduces the probability of self-intersection, and circularization may take many dynamical timescales, delaying the onset of flare activity. An examination of the particle dynamics suggests that quasi-periodic flares with short durations, produced when the center of the tidal stream passes pericenter, may occur in the early-time light curve. The late-time light curve may still show power-law behavior which is generic to disk accretion processes. The detection triggers for future surveys should be extended to capture such "non-standard" short-term flaring activity before the event enters the asymptotic decay phase, as this activity is likely to be more sensitive to physical parameters such as the black hole spin.

Inflationary models predict a correlation between primordial density perturbations (scalar metric perturbations) and gravitational waves (tensor metric perturbations) in the form of a scalar-scalar-tensor three-point correlation, or bispectrum in Fourier space. The squeezed limit of this bispectrum implies a quadrupolar asymmetry in the observed local power spectrum for matter and galaxies. Here we show (like others before) that an infrared divergence in the amplitude of this power asymmetry predicted in single-field slow-roll models is canceled by projection effects when considering the observed power spectrum. We then further evaluate the nonzero, but finite, residual quadrupolar power asymmetry that remains after the divergences are canceled. While the quadrupolar power asymmetry is small, it is conceptually important. Our calculation moreover clarifies how the predictions for this power asymmetry may change in models with different scalar-scalar-tensor bispectra, and shows that convincing detection of the quadrupolar power asymmetry would rule out the single-field slow-roll models of inflation.

If during inflation the inflaton couples to a "fossil" field, some new scalar, vector, or tensor field, it typically induces a scalar-scalar-fossil bispectrum. Even if the fossil field leaves no direct physical trace after inflation, it gives rise to correlations between different Fourier modes of the curvature or, equivalently, a nonzero curvature trispectrum, but without a curvature bispectrum. Here we quantify the effects of a fossil field on the cosmic microwave background (CMB) temperature fluctuations in terms of bipolar spherical harmonics (BiPoSHs). The effects of vector and tensor fossils can be distinguished geometrically from those of scalars through the parity of the BiPoSHs they induce. However, the two-dimensional nature of the CMB sky does not allow vectors to be distinguished geometrically from tensors. We estimate the detectability of a signal in terms of the scalar-scalar-fossil coupling for scalar, vector, and tensor fossils, assuming a local-type coupling. We comment on a divergence that arises in the quadrupolar BiPoSH from the scalar-scalar-tensor correlation in single-field slow-roll inflation.

Physical models for the hemispherical power asymmetry in the cosmic microwave background (CMB) reported by the Planck Collaboration must satisfy CMB constraints to the homogeneity of the Universe and quasar constraints to power asymmetries. We survey a variety of models for the power asymmetry and show that consistent models include a modulated scale-dependent isocurvature contribution to the matter power spectrum or a modulation of the reionization optical depth, gravitational-wave amplitude, or scalar spectral index. We propose further tests to distinguish between the different scenarios.

Total-angular-momentum (TAM) waves provide a set of basis functions for scalar, vector, and tensor fields that can be used in place of plane waves and that reflect the rotational symmetry of the spherical sky. Here we discuss three-point correlation functions, or bispectra in harmonic space, for scalar, vector, and tensor fields in terms of TAM waves. The Wigner-Eckart theorem dictates that the expectation value, assuming statistical isotropy, of the product of three TAM waves is the product of a Clebsch-Gordan coefficient (or Wigner-3j symbol) times a function only of the total-angular-momentum quantum numbers. Here we show how this works, and we provide explicit expressions relating the bispectra for TAM waves in terms of the more commonly used Fourier-space bispectra. This formalism will be useful to simplify calculations of projections of three-dimensional bispectra onto the spherical sky.

Most calculations in cosmological perturbation theorydecompose those perturbations into plane waves (Fourier modes). However, for some calculations, particularly those involving observations performed on a spherical sky, a decomposition into waves of fixed total angular momentum (TAM) may be more appropriate. Here we introduce TAM waves, solutions of fixed total angular momentum to the Helmholtz equation, for three-dimensional scalar, vector, and tensor fields. The vector TAM waves of given total angular momentum can be decomposed further into a set of three basis functions of fixed orbital angular momentum (OAM), a set of fixed helicity, or a basis consisting of a longitudinal (L) and two transverse (E and B) TAM waves. The symmetric traceless rank-2 tensor TAM waves can be similarly decomposed into a basis of fixed OAM or fixed helicity, or a basis that consists of a longitudinal (L), two vector (VE and VB, of opposite parity), and two tensor (TE and TB, of opposite parity) waves. We show how all of the vector and tensor TAM waves can be obtained by applying derivative operators to scalar TAM waves. This operator approach then allows one to decompose a vector field into three covariant scalar fields for the L, E, and B components and symmetric-traceless-tensor fields into five covariant scalar fields for the L, VE, VB, TE, and TB components. We provide projections of the vector and tensor TAM waves onto vector and tensor spherical harmonics. We provide calculational detail to facilitate the assimilation of this formalism into cosmological calculations. As an example, we calculate the power spectra of the deflection angle for gravitational lensing by density perturbations and by gravitational waves. We comment on an alternative approach to CMB fluctuations based on TAM waves. Our work may have applications elsewhere in field theory and in general relativity.

We have calculated the equilibrium properties of a star in a circular, equatorial orbit about a Super-Massive Black Hole (SMBH), when the star fills and overflows its Roche lobe. The mass transfer time scale is anticipated to be long compared with the dynamical time and short compared with the thermal time of the star, so that the entropy as a function of the interior mass is conserved. We have studied how the stellar entropy, pressure, radius, mean density, and orbital angular momentum vary when the star is evolved adiabatically, for a representative set of stars. We have shown that the stellar orbits change with the stellar mean density. Therefore, sun-like stars, upper main sequence stars and red giants will spiral inward and then outward with respect to the hole in this stable mass transfer process, while lower main sequence stars, brown dwarfs and white dwarfs will always spiral outward.

In this paper we consider Roche accretion in an Extreme Mass-Ratio Inspiral (EMRI) binary system formed by a star orbiting a massive black hole. The ultimate goal is to detect the mass and spin of the black hole and provide a test of general relativity in the strong-field regime from the resultant quasi-periodic signals. Before accretion starts, the stellar orbit is presumed to be circular and equatorial, and shrinks due to gravitational radiation. New fitting formulae are presented for the inspiral time and the radiation-reaction torque in the relativistic regime. If the inspiralling star fills its Roche lobe outside the Innermost Stable Circular Orbit (ISCO) of the hole, gas will flow through the inner Lagrange point (L1) to the hole. We give new relativistic interpolation formulae for the volume enclosed by the Roche lobe. If this mass-transfer happens on a time scale faster than the thermal time scale but slower than the dynamical time scale, the star will evolve adiabatically, and, in most cases, will recede from the hole while filling its Roche lobe. We calculate how the stellar orbital period and mass-transfer rate will change through the "Roche evolution" for various types of stars in the relativistic regime. We envisage that the mass stream eventually hits the accretion disc, where it forms a hot spot orbiting the hole and may ultimately modulate the luminosity with the stellar orbital frequency. The observability of such a modulation is discussed along with a possible interpretation of an intermittent 1 hour period in the X-ray emission of RE J1034+396.

We present simulated results of quasi-periodic flares generated by the inelastic collisions of a star bound to a super-massive black hole (SMBH) and its attendant accretion disc. We show that the behavior of the quasi-periodicity is affected by the mass and spin of the black hole and the orbital elements of the stellar orbit. We also evaluate the possibility of extracting useful information on these parameters and verifying the character of the Kerr metric from such quasi-periodic signals. Comparisons are made with the observed optical outbursts of OJ287, infrared flares from the Galactic center and X-ray variability in RE J1034+396.