We present genetIC, a new code for generating initial conditions for cosmological N-body simulations. The code allows precise, user-specified alterations to be made to arbitrary regions of the simulation (while maintaining consistency with the statistical ensemble). These "genetic modifications" allow, for example, the history, mass, or environment of a target halo to be altered in order to study the effect on their evolution. The code natively supports initial conditions with nested zoom regions at progressively increasing resolution. Modifications in the high-resolution region must propagate self-consistently onto the lower resolution grids; to enable this while maintaining a small memory footprint, we introduce a Fourier-space filtering approach to generating fields at variable resolution. Due to a close correspondence with modifications, constrained initial conditions can also be produced by genetIC (for example with the aim of matching structures in the local Universe). We test the accuracy of modifications performed within zoom initial conditions. The code achieves sub-percent precision, which is easily sufficient for current applications in galaxy formation.

Comets spend most of their lives at large distances from any star, during which time their interior compositions remain relatively unaltered. Cometary observations can therefore provide direct insight into the chemistry that occurred during their birth at the time of planet formation. To-date, there have been no confirmed observations of parent volatiles (gases released directly from the nucleus) of a comet from any planetary system other than our own. Here we present high-resolution, interferometric observations of 2I/Borisov, the first confirmed interstellar comet, obtained using the Atacama Large Millimeter/submillimeter Array (ALMA) on 15th-16th December 2019. Our observations reveal emission from hydrogen cyanide (HCN), and carbon monoxide (CO), coincident with the expected position of 2I/Borisov's nucleus, with production rates Q(HCN)=$(7.0\pm1.1)\times10^{23}$ s$^{-1}$ and Q(CO)=$(4.4\pm0.7)\times10^{26}$ s$^{-1}$. While the HCN abundance relative to water (0.06-0.16%) appears similar to that of typical, previously observed comets in our Solar System, the abundance of CO (35-105%) is among the highest observed in any comet within 2 au of the Sun. This shows that 2I/Borisov must have formed in a relatively CO-rich environment - probably beyond the CO ice-line in the very cold, outer regions of a distant protoplanetary accretion disk, as part of a population of small, icy bodies analogous to our Solar System's own proto-Kuiper Belt.

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.

While tidal disruption events (TDEs) have long been heralded as laboratories for the study of quiescent black holes, the small number of known TDEs and uncertainties in their emission mechanism have hindered progress towards this promise. Here present 17 new TDEs that have been detected recently by the Zwicky Transient Facility along with Swift UV and X-ray follow-up observations. Our homogeneous analysis of the optical/UV light curves, including 22 previously known TDEs from the literature, reveals a clean separation of light curve properties with spectroscopic class. The TDEs with Bowen fluorescence features in their optical spectra have smaller blackbody radii, as well as longer rise times and higher disruption rates compared to the rest of the sample. The Bowen fluorescence mechanism requires a high density which can be reached at smaller radii, which in turn yields longer diffusion timescales. Thus, the difference in rise times suggests the pre-peak TDE light curves are governed not by the fallback timescale, but instead by the diffusion of photons through the tidal debris. The small subset of TDEs that show only helium emission lines in their spectra have the longest rise times, the highest luminosities and the lowest rates. We also report, for the first time, the detection of soft X-ray flares from a TDE on day timescales. Based on the fact the flares peak at a luminosity similar to the optical/UV blackbody luminosity, we attribute them to brief glimpses through a reprocessing layer that otherwise obscures the inner accretion flow.

Comet C/2016 R2 (PanSTARRS) has a peculiar volatile composition, with CO being the dominant volatile as opposed to H$_2$O and one of the largest N$_2$/CO ratios ever observed in a comet. Using observations obtained with the \textit{Spitzer Space Telescope}, NASA's Infrared Telescope Facility, the 3.5-meter ARC telescope at Apache Point Observatory, the Discovery Channel Telescope at Lowell Observatory, and the Arizona Radio Observatory 10-m Submillimeter Telescope we quantified the abundances of 12 different species in the coma of R2 PanSTARRS. We confirm the high abundances of CO and N$_2$ and heavy depletions of H$_2$O, HCN, CH$_3$OH, and H$_2$CO compared to CO reported by previous studies. We provide the first measurements (or most sensitive measurements/constraints) on H$_2$O, CO$_2$, CH$_4$, C$_2$H$_6$, OCS, C$_2$H$_2$, and NH$_3$, all of which are depleted relative to CO by at least one to two orders of magnitude compared to values commonly observed in comets. The observed species also show strong enhancements relative to H$_2$O, and even when compared to other species like CH$_4$ or CH$_3$OH most species show deviations from typical comets by at least a factor of two to three. The only mixing ratios found to be close to typical are CH$_3$OH/CO$_2$ and CH$_3$OH/CH$_4$. While R2 PanSTARRS was located at a heliocentric distance of 2.8 AU at the time of our observations in January/February 2018, we argue that this alone cannot account for the peculiar observed composition of this comet and therefore must reflect its intrinsic composition. We discuss possible implications for this clear outlier in compositional studies of comets obtained to date, and encourage future dynamical and chemical modeling in order to better understand what the composition of R2 PanSTARRS tells us about the early Solar System.

We report the discovery of non-stellar hydrogen Balmer and metastable helium absorption lines accompanying a transient, high-velocity (0.05$c$) broad absorption line (BAL) system in the optical spectra of the tidal disruption event (TDE) AT2018zr ($z=0.071$). In the HST UV spectra, absorption of high- and low-ionization lines are also present at this velocity, making AT2018zr resemble a low-ionization broad absorption line (LoBAL) QSO. We conclude that these transient absorption features are more likely to arise in fast outflows produced by the TDE than absorbed by the unbound debris. In accordance with the outflow picture, we are able to reproduce the flat-topped H$\alpha$ emission in a spherically expanding medium, without invoking the typical prescription of an elliptical disk. We also report the appearance of narrow ($\sim$1000~km~s$^{-1}$) NIII$\lambda$4640, HeII$\lambda4686$, H$\alpha$, and H$\beta$, emission in the late-time optical spectra of AT2018zr, which may be a result of UV continuum hardening at late time as observed by Swift. Including AT2018zr, we find a high association rate (3 out of 4) of BALs in the UV spectra of TDEs. This suggests that outflows may be ubiquitous among TDEs and may be less sensitive to viewing angle effects compared to QSO outflows.

The field of time-domain astrophysics has entered the era of Multi-messenger Astronomy (MMA). One key science goal for the next decade (and beyond) will be to characterize gravitational wave (GW) and neutrino sources using the next generation of Extremely Large Telescopes (ELTs). These studies will have a broad impact across astrophysics, informing our knowledge of the production and enrichment history of the heaviest chemical elements, constrain the dense matter equation of state, provide independent constraints on cosmology, increase our understanding of particle acceleration in shocks and jets, and study the lives of black holes in the universe. Future GW detectors will greatly improve their sensitivity during the coming decade, as will near-infrared telescopes capable of independently finding kilonovae from neutron star mergers. However, the electromagnetic counterparts to high-frequency (LIGO/Virgo band) GW sources will be distant and faint and thus demand ELT capabilities for characterization. ELTs will be important and necessary contributors to an advanced and complete multi-messenger network.

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.

We present Zwicky Transient Facility (ZTF) observations of the tidal disruption flare AT2018zr/PS18kh reported by Holoien et al. and detected during ZTF commissioning. The ZTF light curve of the tidal disruption event (TDE) samples the rise-to-peak exceptionally well, with 50 days of g- and r-band detections before the time of maximum light. We also present our multi-wavelength follow-up observations, including the detection of a thermal (kT~100 eV) X-ray source that is two orders of magnitude fainter than the contemporaneous optical/UV blackbody luminosity, and a stringent upper limit to the radio emission. We use observations of 128 known active galactic nuclei (AGN) to assess the quality of the ZTF astrometry, finding a median host-flare distance of 0.2" for genuine nuclear flares. Using ZTF observations of variability from known AGN and supernovae we show how these sources can be separated from TDEs. A combination of light-curve shape, color, and location in the host galaxy can be used to select a clean TDE sample from multi-band optical surveys such as ZTF or LSST.

We present the first extensive radio to gamma-ray observations of a fast-rising blue optical transient (FBOT), AT2018cow, over its first ~100 days. AT2018cow rose over a few days to a peak luminosity $L_{pk}\sim4\times 10^{44}$ erg/s exceeding those of superluminous supernovae (SNe), before declining as $\propto t^{-2}$. Initial spectra at $\lesssim 15$ days were mostly featureless and indicated large expansion velocities v~0.1c and temperatures reaching 30000 K. Later spectra revealed a persistent optically-thick photosphere and the emergence of H and He emission features with v~sim 4000 km/s with no evidence for ejecta cooling. Our broad-band monitoring revealed a hard X-ray spectral component at $E\ge 10$ keV, in addition to luminous and highly variable soft X-rays, with properties unprecedented among astronomical transients. An abrupt change in the X-ray decay rate and variability appears to accompany the change in optical spectral properties. AT2018cow showed bright radio emission consistent with the interaction of a blastwave with $v_{sh}$~0.1c with a dense environment ($\dot M\sim10^{-3}-10^{-4}\,M_{\odot}yr^{-1}$ for $v_w=1000$ km\s). While these properties exclude Ni-powered transients, our multi-wavelength analysis instead indicates that AT2018cow harbored a "central engine", either a compact object (magnetar or black hole) or an embedded internal shock produced by interaction with a compact, dense circumstellar medium. The engine released $\sim10^{50}-10^{51.5}$ erg over $\sim10^3-10^5$ s and resides within low-mass fast-moving material with equatorial-polar density asymmetry ($M_{ej,fast}\lesssim0.3\,\rm{M_{\odot}}$). Successful SNe from low-mass H-rich stars (like electron-capture SNe) or failed explosions from blue supergiants satisfy these constraints. Intermediate-mass black-holes are disfavored by the large environmental density probed by the radio observations.

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.

We investigate line formation in gas that is outflowing and optically thick to electron scattering, as may be expected following the tidal disruption of a star by a super-massive black hole. Using radiative transfer calculations, we show that the optical line profiles produced by expanding TDE outflows are most likely primarily emission features, rather than the P-Cygni profiles seen in most supernova spectra. This is a result of the high line excitation temperatures in the highly irradiated TDE gas. The outflow kinematics cause the emission peak to be blueshifted and to have an asymmetric red wing. Such features have been observed in some TDE spectra, and we propose that these may be signatures of outflows. We also show that non-coherent scattering off of hot electrons can broaden the emission lines by $\sim 10000$~km~s$^{-1}$, such that in some TDEs the line width may be set by the electron scattering optical depth rather than the gas kinematics. The scattering broadened line profiles produce distinct, wing-shaped profiles that are similar to those observed in some TDE spectra. The narrowing of the emission lines over time in these observed events may be related to a drop in density rather than a drop in line-of-sight velocity.

We report the discovery by the intermediate Palomar Transient Factory (iPTF) of a candidate tidal disruption event (TDE) iPTF16axa at $z=0.108$, and present its broadband photometric and spectroscopic evolution from 3 months of follow-up observations with ground-based telescopes and Swift. The light curve is well fitted with a $t^{-5/3}$ decay, and we constrain the rise-time to peak to be $<$49 rest-frame days after disruption, which is roughly consistent with the fallback timescale expected for the $\sim 5\times$10$^{6}$ $M_\odot$ black hole inferred from the stellar velocity dispersion of the host galaxy. The UV and optical spectral energy distribution (SED) is well described by a constant blackbody temperature of T$\sim$ 3$\times$10$^4$ K over the monitoring period, with an observed peak luminosity of 1.1$\times$10$^{44}$ erg s$^{-1}$. The optical spectra are characterized by a strong blue continuum and broad HeII and H$\alpha$ lines characteristic of TDEs. We compare the photometric and spectroscopic signatures of iPTF16axa with 11 TDE candidates in the literature with well-sampled optical light curves. Based on a single-temperature fit to the optical and near-UV photometry, most of these TDE candidates have peak luminosities confined between log(L [erg s$^{-1}$]) = 43.4-44.4, with constant temperatures of a few $\times 10^{4}$ K during their power-law declines, implying blackbody radii on the order of ten times the tidal disruption radius, that decrease monotonically with time. For TDE candidates with hydrogen and helium emission, the high helium-to-hydrogen ratios suggest that the emission arises from high-density gas, where nebular arguments break down. We find no correlation between the peak luminosity and the black hole mass, contrary to the expectations for TDEs to have $\dot{M} \propto M_{\rm BH}^{-1/2}$.

We show how the interplay between active galactic nuclei (AGN) and merger history determines whether a galaxy quenches star formation at high redshift. We first simulate, in a full cosmological context, a galaxy of total dynamical mass $10^{12}\,M_{\odot}$ at $z=2$. Then we systematically alter the accretion history of the galaxy by minimally changing the linear overdensity in the initial conditions. This "genetic modification" approach allows the generation of three sets of $\Lambda$CDM initial conditions leading to maximum merger ratios of 1:10, 1:5 and 2:3 respectively. The changes leave the final halo mass, large scale structure and local environment unchanged, providing a controlled numerical experiment. Interaction between the AGN physics and mergers in the three cases lead respectively to a star-forming, temporarily-quenched and permanently-quenched galaxy. However the differences do not primarily lie in the black hole accretion rates, but in the kinetic effects of the merger: the galaxy is resilient against AGN feedback unless its gaseous disk is first disrupted. Typical accretion rates are comparable in the three cases, falling below $0.1\,M_{\odot}$ yr$^{-1}$, equivalent to around $2\%$ of the Eddington rate or $10^{-3}$ times the pre-quenching star formation rate, in agreement with observations. This low level of black hole accretion can be sustained even when there is insufficient dense cold gas for star formation. Conversely, supernova feedback is too distributed to generate outflows in high-mass systems, and cannot maintain quenching over periods longer than the halo gas cooling time.

Observations of luminous flares resulting from the possible tidal disruption of stars by supermassive black holes have raised a number of puzzles. Outstanding questions include the origin of the optical and ultraviolet (UV) flux, the weakness of hydrogen lines in the spectrum, and the occasional simultaneous observation of x-rays. Here we study the emission from tidal disruption events (TDEs) produced as radiation from black hole accretion propagates through an extended, optically thick envelope formed from stellar debris. We analytically describe key physics controlling spectrum formation, and present detailed radiative transfer calculations that model the spectral energy distribution (SED) and optical line strengths of TDEs near peak brightness. The steady-state transfer is coupled to a non local thermodynamic equilibrium treatment of the excitation and ionization states of hydrogen, helium and oxygen (as a representative metal). Our calculations show how an extended envelope can reprocess a fraction of soft x-rays and produce the observed optical fluxes of order 10^43 ergs per second. Variations in the mass or size of the envelope may help explain how the optical flux changes over time with roughly constant color. For high enough accretion luminosities, x-rays can highly ionize the reprocessing region and escape to be observed simultaneously with the optical flux, producing an SED not described by a single blackbody. Due to optical depth effects, hydrogen Balmer line emission is often strongly suppressed relative to helium line emission (with HeII-to-H line ratios of at least 5:1 in some cases) even in the disruption of a solar-composition star. We discuss the implications of our results to understanding the type of stars destroyed in TDEs and the physical processes responsible for producing the observed flares.

We present a Hubble Space Telescope STIS spectrum of ASASSN-14li, the first rest-frame UV spectrum of a tidal disruption flare (TDF). The underlying continuum is well fit by a blackbody with $T_{\mathrm{UV}} = 3.5 \times 10^{4}$ K, an order of magnitude smaller than the temperature inferred from X-ray spectra (and significantly more precise than previous efforts based on optical and near-UV photometry). Super-imposed on this blue continuum, we detect three classes of features: narrow absorption from the Milky Way (probably a high-velocity cloud), and narrow absorption and broad (FWHM $\approx 2000$-8000 km s$^{-1}$) emission lines at/near the systemic host velocity. The absorption lines are blueshifted with respect to the emission lines by $\Delta v = -(250$-400) km s$^{-1}$. Due both to this velocity offset and the lack of common low-ionization features (Mg II, Fe II), we argue these arise from the same absorbing material responsible for the low-velocity outflow discovered at X-ray wavelengths. The broad nuclear emission lines display a remarkable abundance pattern: N III], N IV], He II are quite prominent, while the common quasar emission lines of C III] and Mg II are weak or entirely absent. Detailed modeling of this spectrum will help elucidate fundamental questions regarding the nature of the emission process(es) at work in TDFs, while future UV spectroscopy of ASASSN-14li would help to confirm (or refute) the previously proposed connection between TDFs and "N-rich" quasars.

We introduce and explore "paired" cosmological simulations. A pair consists of an A and B simulation with initial conditions related by the inversion $\delta_A(x, t_{initial})=-\delta_B(x,t_{initial})$ (underdensities substituted for overdensities and vice versa). We argue that the technique is valuable for improving our understanding of cosmic structure formation. The A and B fields are by definition equally likely draws from {\Lambda}CDM initial conditions, and in the linear regime evolve identically up to the overall sign. As non-linear evolution takes hold, a region that collapses to form a halo in simulation A will tend to expand to create a void in simulation B. Applications include (i) contrasting the growth of A-halos and B-voids to test excursion-set theories of structure formation; (ii) cross-correlating the density field of the A and B universes as a novel test for perturbation theory; and (iii) canceling error terms by averaging power spectra between the two boxes. Generalizations of the method to more elaborate field transformations are suggested.

We propose a method to generate `genetically-modified' (GM) initial conditions for high-resolution simulations of galaxy formation in a cosmological context. Building on the Hoffman-Ribak algorithm, we start from a reference simulation with fully random initial conditions, then make controlled changes to specific properties of a single halo (such as its mass and merger history). The algorithm demonstrably makes minimal changes to other properties of the halo and its environment, allowing us to isolate the impact of a given modification. As a significant improvement over previous work, we are able to calculate the abundance of the resulting objects relative to the reference simulation. Our approach can be applied to a wide range of cosmic structures and epochs; here we study two problems as a proof-of-concept. First, we investigate the change in density profile and concentration as the collapse time of three individual halos are varied at fixed final mass, showing good agreement with previous statistical studies using large simulation suites. Second, we modify the $z=0$ mass of halos to show that our theoretical abundance calculations correctly recover the halo mass function. The results demonstrate that the technique is robust, opening the way to controlled experiments in galaxy formation using hydrodynamic zoom simulations.

Galaxies and the dark matter halos that host them are not spherically symmetric, yet spherical symmetry is a helpful simplifying approximation for idealised calculations and analysis of observational data. The assumption leads to an exact conservation of angular momentum for every particle, making the dynamics unrealistic. But how much does that inaccuracy matter in practice for analyses of stellar distribution functions, collisionless relaxation, or dark matter core-creation? We provide a general answer to this question for a wide class of aspherical systems; specifically, we consider distribution functions that are "maximally stable", i.e. that do not evolve at first order when external potentials (which arise from baryons, large scale tidal fields or infalling substructure) are applied. We show that a spherically-symmetric analysis of such systems gives rise to the false conclusion that the density of particles in phase space is ergodic (a function of energy alone). Using this idea we are able to demonstrate that: (a) observational analyses that falsely assume spherical symmetry are made more accurate by imposing a strong prior preference for near-isotropic velocity dispersions in the centre of spheroids; (b) numerical simulations that use an idealised spherically-symmetric setup can yield misleading results and should be avoided where possible; and (c) triaxial dark matter halos (formed in collisionless cosmological simulations) nearly attain our maximally-stable limit, but their evolution freezes out before reaching it.

We derive robust constraints on primordial non-Gaussianity (PNG) using the clustering of 800,000 photometric quasars from the Sloan Digital Sky Survey in the redshift range $0.5<z<3.5$. These measurements rely on the novel technique of {\it extended mode projection} to control the impact of spatially-varying systematics in a robust fashion, making use of blind analysis techniques. This allows the accurate measurement of quasar halo bias at the largest scales, while discarding as little as possible of the data. The standard local-type PNG parameters $f_\mathrm{NL}$ and $g_\mathrm{NL}$ both imprint a $k^{-2}$ scale-dependent effect in the bias. Constraining these individually, we obtain $-49<f_\mathrm{NL}<31$ and $-2.7\times10^5<g_\mathrm{NL}<1.9\times10^5$, while their joint constraints lead to $-105<f_\mathrm{NL}<72$ and $-4.0\times10^5<g_\mathrm{NL}<4.9\times10^5$ (all at 95% CL) . Introducing a running parameter $n_{f_\mathrm{NL}}$ to constrain $b(k) \propto k^{-2+n_{f_\mathrm{NL}}}$ and a generalised PNG amplitude $\tilde{f}_\mathrm{NL}$, we obtain $-45.5 \exp({3.7\, n_{f_\mathrm{NL}}}) < \tilde{f}_\mathrm{NL} < 34.4 \exp({3.3\, n_{f_\mathrm{NL}}})$ at 95% CL. These results incorporate uncertainties in the cosmological parameters, redshift distributions, shot noise, and the bias prescription used to relate the quasar clustering to the underlying dark matter. These are the strongest constraints obtained to date on PNG using a single population of large-scale structure tracers, and are already at the level of pre-{\it Planck} constraints from the cosmic microwave background. A conservative forecast for a {\it Large Synoptic Survey Telescope}-like survey incorporating mode projection yields $\sigma(f_\mathrm{NL}) \sim 5$ -- competitive with the {\it Planck} result -- highlighting the power of upcoming large scale structure surveys to probe the initial conditions of the universe.

We explore the application of Monte Carlo transport methods to solving coupled radiation-hydrodynamics problems. We use a time-dependent, frequency-dependent, 3-dimensional radiation transport code, that is special relativistic and includes some detailed microphysical interactions such as resonant line scattering. We couple the transport code to two different 1-dimensional (non-relativistic) hydrodynamics solvers: a spherical Lagrangian scheme and a Eulerian Godunov solver. The gas-radiation energy coupling is treated implicitly, allowing us to take hydrodyanimcal time-steps that are much longer than the radiative cooling time. We validate the code and assess its performance using a suite of radiation hydrodynamical test problems, including ones in the radiation energy dominated regime. We also develop techniques that reduce the noise of the Monte Carlo estimated radiation force by using the spatial divergence of the radiation pressure tensor. The results suggest that Monte Carlo techniques hold promise for simulating the multi-dimensional radiation hydrodynamics of astrophysical systems.

Many ultracompact HII regions exhibit a cometary morphology in radio continuum emission. In such regions, a young massive star is probably ablating, through its ultraviolet radiation, the molecular cloud clump that spawned it. On one side of the star, the radiation drives an ionization front that stalls in dense molecular gas. On the other side, ionized gas streams outward into the more rarefied environment. This wind is underpressured with respect to the neutral gas. The difference in pressure draws in more cloud material, feeding the wind until the densest molecular gas is dissipated. Recent, time-dependent simulations of massive stars turning on within molecular gas show the system evolving in a direction similar to that just described. Here, we explore a semi-analytic model in which the wind is axisymmetric and has already achieved a steady state. Adoption of this simplified picture allows us to study the dependence of both the wind and its bounding ionization front on the stellar luminosity, the peak molecular density, and the displacement of the star from the center of the clump. For typical parameter values, the wind accelerates transonically to a speed of about 15 km / s, and transports mass outward at a rate of 10^-4 solar masses per year. Stellar radiation pressure acts to steepen the density gradient of the wind.

Observational and theoretical arguments suggest that the momentum carried in mass outflows from AGN can reach several times L / c, corresponding to outflow rates of hundreds of solar masses per year. Radiation pressure on lines alone may not be sufficient to provide this momentum deposition, and the transfer of reprocessed IR radiation in dusty nuclear gas has been postulated to provide the extra enhancement. The efficacy of this mechanism, however, will be sensitive to multi-dimensional effects such as the tendency for the reprocessed radiation to preferentially escape along sight-lines of lower column density. We use Monte Carlo radiative transfer calculations to determine the radiation force on dusty gas residing within approximately 10 parsecs from an accreting super-massive black hole. We calculate the net rate of momentum deposition in the surrounding gas and estimate the mass-loss rate in the resulting outflow as a function of solid angle for different black hole luminosities, sightline-averaged column densities, clumping parameters, and opening angles of the dusty gas. We find that these dust-driven winds carry momentum fluxes of 1-5 times L / c and correspond to mass-loss rates of 10-100 solar masses per year for a 10^8 solar mass black hole radiating at or near its Eddington limit. These results help to explain the origin of high velocity molecular and atomic outflows in local ULIRGs, and can inform numerical simulations of galaxy evolution including AGN feedback.

The scale-dependent galaxy bias generated by primordial non-Gaussianity (PNG) can be used to detect and constrain deviations from standard single-field inflation. The strongest signal is expected in the local model for PNG, where the amplitude of non-Gaussianity can be expressed by a set of parameters (fnl, gnl, ...). Current observational constraints from galaxy clustering on fnl and gnl assume that the others PNG parameters are vanishing. Using two sets of cosmological N-body simulations where both fnl and gnl are non-zero, we show that this strong assumption generally leads to biased estimates and spurious redshift dependencies of the parameters. Additionally, if the signs of fnl and gnl are opposite, the amplitude of the scale-dependent bias is reduced, possibly leading to a false null detection. Finally we show that model selection techniques like the Bayesian evidence can (and should) be used to determine if more than one PNG parameter is required by the data.

We test third-order standard perturbation theory (SPT) as an approximation to non-linear cosmological structure formation. A novel approach is used to numerically calculate the three-dimensional dark matter density field using SPT from the initial conditions of two high-resolution cosmological simulations. The calculated density field is compared to the non-linear dark matter field of the simulations both point-by-point and statistically. For smoothing scales above 8 Mpc/h it shows a good agreement up to redshift 0. We present a simple fitting formula to relate the linear and non-linear density contrast that accurately recovers the non-linear time evolution for 0 <= z <= 10 at the per cent level. To address the problem of biasing between the matter field and the haloes identified in the simulation, we employ the Eulerian local bias model (ELB), including non-linear bias up to the third order. The bias parameters are obtained by fitting a scatter plot of halo and matter density (both from the simulation and from SPT). Using these bias parameters, we can reconstruct the halo density field. We find that this reconstruction is not able to capture all the details of the halo distribution. We investigate how well the large scale bias can be described by a constant and if it corresponds to the linear bias parameter b_1 of the local bias model. We also discuss how well the halo-halo power spectrum and the halo-mass cross spectrum from the reconstructed halo density field agree with the corresponding statistics from the simulation. The results show that while SPT is an excellent approximation for the matter field for suitably large smoothing scales even at redshift 0, the ELB model can only account for some of the properties of the halo density field.

PSpectRe is a C++ program that uses Fourier-space pseudo-spectral methods to evolve interacting scalar fields in an expanding universe. PSpectRe is optimized for the analysis of parametric resonance in the post-inflationary universe, and provides an alternative to finite differencing codes, such as Defrost and LatticeEasy. PSpectRe has both second- (Velocity-Verlet) and fourth-order (Runge-Kutta) time integrators. Given the same number of spatial points and/or momentum modes, PSpectRe is not significantly slower than finite differencing codes, despite the need for multiple Fourier transforms at each timestep, and exhibits excellent energy conservation. Further, by computing the post-resonance equation of state, we show that in some circumstances PSpectRe obtains reliable results while using substantially fewer points than a finite differencing code. PSpectRe is designed to be easily extended to other problems in early-universe cosmology, including the generation of gravitational waves during phase transitions and pre-inflationary bubble collisions. Specific applications of this code will be pursued in future work.