We present a recent Chandra observation of the quiescent low-mass X-ray binary containing a neutron star, located in the globular cluster M30. We fit the thermal emission from the neutron star to extract its mass and radius. We find no evidence of flux variability between the two observations taken in 2001 and 2017, nor between individual 2017 observations, so we analyse them together to increase the signal to noise. We perform simultaneous spectral fits using standard light-element composition atmosphere models (hydrogen or helium), including absorption by the interstellar medium, correction for pile-up of X-ray photons on the detector, and a power-law for count excesses at high photon energy. Using a Markov-chain Monte Carlo approach, we extract mass and radius credible intervals for both chemical compositions of the atmosphere: $R_{\textrm{NS}}=7.94^{+0.76}_{-1.21}$ km and $M_{\textrm{NS}}<1.19$ M$_{\odot}$ assuming pure hydrogen, and $R_{\textrm{NS}}=10.50^{+2.88}_{-2.03}$ km and $M_{\textrm{NS}}<1.78$ M$_{\odot}$ for helium, where the uncertainties represent the 90% credible regions. For H, the small radius is difficult to reconcile with most current nuclear physics models (especially for nucleonic equations of state) and with other measurements of neutron star radii, with recent preferred values generally in the 11-14 km range. Whereas for He, the measured radius is consistent with this range. We discuss possible sources of systematic uncertainty that may result in an underestimation of the radius, identifying the presence of surface temperature inhomogeneities as the most relevant bias. According to this, we conclude that either the atmosphere is composed of He, or it is a H atmosphere with a significant contribution of hot spots to the observed radiation.

The Galactic Center black hole Sgr A* shows significant variability and flares in the submillimeter, infrared, and X-ray wavelengths. Owing to its exquisite resolution in the IR bands, the GRAVITY experiment for the first time spatially resolved the locations of three flares and showed that a bright region moves in ellipse-like trajectories close to but offset from the black hole over the course of each event. We present a model for plasmoids that form during reconnection events and orbit in the coronal region around a black hole to explain these observations. We utilize general-relativistic radiative transfer calculations that include effects from finite light travel time, plasmoid motion, particle acceleration, and synchrotron cooling and obtain a rich structure in the flare lightcurves. This model can naturally account for the observed motion of the bright regions observed by the GRAVITY experiment and the offset between the center of the centroid motion and the position of the black hole. It also explains why some flares may be double-peaked while others have only a single peak and uncovers a correlation between the structure in the lightcurve and the location of the flare. Finally, we make predictions for future observations of flares from the inner accretion flow of Sgr A* that will provide a test of this model.

The Event Horizon Telescope (EHT) is taking the first images of black holes resolved at horizon scales to measure their shadows and probe accretion physics. A promising avenue for testing the hypothesis that astrophysical black holes are described by the Kerr solution to Einstein's equations is to compare the size and shape of the shadow a black hole casts on the surrounding emission to the predictions of the Kerr metric. We develop here an efficient parametric framework to perform this test. We carry out ray-tracing simulations for several parametrized non-Kerr metrics to create a large data set of non-Kerr shadows that probe the allowed parameter space for the free parameters of each metric. We then perform principal components analysis (PCA) on this set of shadows and show that only a small number of components are needed to accurately reconstruct all shadows within the set. We further show that the amplitude of the PCA components are smoothly related to the free parameters in the metrics, and therefore, that these PCA components can be fit to EHT observations in order to place constraints on the free parameters of these metrics that will help quantify any potential deviations from the Kerr solution.

Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

Thermal X-ray emission from rotation-powered pulsars is believed to originate from localized "hotspots" on the stellar surface occurring where large-scale currents from the magnetosphere return to heat the atmosphere. Lightcurve modeling has primarily been limited to simple models, such as circular antipodal emitting regions with constant temperature. We calculate more realistic temperature distributions within the polar caps, taking advantage of recent advances in magnetospheric theory, and we consider their effect on the predicted lightcurves. The emitting regions are non-circular even for a pure dipole magnetic field, and the inclusion of an aligned magnetic quadrupole moment introduces a north-south asymmetry. As the aligned quadrupole moment is increased, one hotspot grows in size before becoming a thin ring surrounding the star. For the pure dipole case, moving to the more realistic model changes the lightcurves by $5-10\%$ for millisecond pulsars, helping to quantify the systematic uncertainty present in current dipolar models. Including the quadrupole gives considerable freedom in generating more complex lightcurves. We explore whether these simple dipole+quadrupole models can account for the qualitative features of the lightcurve of PSR J0437$-$4715.

Electron acceleration during magnetic reconnection is thought to play a key role in time-variable high-energy emission from astrophysical systems. By means of particle-in-cell simulations of trans-relativistic reconnection, we investigate electron injection and acceleration mechanisms in low-$\beta$ electron-proton plasmas. We set up a diversity of density and field structures (e.g., X-points and plasmoids) by varying the guide field strength and choosing whether to trigger reconnection or let it spontaneously evolve. We show that the number of X-points and plasmoids controls the efficiency of electron acceleration, with more X-points leading to a higher efficiency. Using on-the-fly acceleration diagnostics, we also show that the non-ideal electric fields associated with X-points play a critical role in the first stages of electron acceleration. As a further diagnostic, we include two populations of test particles that selectively experience only certain components of electric fields. We find that the out-of-plane component of the parallel electric field determines the hardness of the high-energy tail of the electron energy distribution. These results further our understanding of electron acceleration in this regime of magnetic reconnection and have implications for realistic models of black hole accretion flows.

Massive 10^6-10^10 Msun black holes (BHs) are ubiquitous in local galactic nuclei. They were common by the time the Universe is several Gyr old, and many of them were in place within the first 1~Gyr after the Big Bang. Their quick assembly has been attributed to mechanisms such as the rapid collapse of gas into the nuclei of early protogalaxies, accretion and mergers of stellar-mass BHs accompanying structure formation at early times, and the runaway collapse of early, ultra-dense stellar clusters. The origin of the early massive BHs remains an intriguing and long-standing unsolved puzzle in astrophysics. Here we discuss strategies for discerning between BH seeding models using electromagnetic observations. We argue that the most direct answers will be obtained through detection of BHs with masses M<10^5 Msun at redshifts z>10, where we expect them to first form. Reaching out to these redshifts and down to these masses is crucial, because BHs are expected to lose the memory of their initial assembly by the time they grow well above 10^5 Msun and are incorporated into higher-mass galaxies. The best way to detect 10^4-10^5 Msun BHs at high redshifts is by a sensitive X-ray survey. Critical constraining power is augmented by establishing the properties and the environments of their host galaxies in deep optical/IR imaging surveys. Required OIR data can be obtained with the JWST and WFIRST missions. The required X-ray flux limits (down to 10^{-19} erg/s/cm^2) are accessible only with a next-generation X-ray observatory which has both high (sub-1") angular resolution and high throughput. A combination of deep X-ray and OIR surveys will be capable of probing several generic markers of the BH seed scenarios, and resolving the long-stanging puzzle of their origin. These electromagnetic observations are also highly synergistic with the information from LISA on high-z BH mergers.

We explore the use of principal component analysis (PCA) to characterize high-fidelity simulations and interferometric observations of the millimeter emission that originates near the horizons of accreting black holes. We show mathematically that the Fourier transforms of eigenimages derived from PCA applied to an ensemble of images in the spatial-domain are identical to the eigenvectors of PCA applied to the ensemble of the Fourier transforms of the images, which suggests that this approach may be applied to modeling the sparse interferometric Fourier-visibilities produced by an array such as the Event Horizon Telescope (EHT). We also show that the simulations in the spatial domain themselves can be compactly represented with a PCA-derived basis of eigenimages allowing for detailed comparisons between variable observations and time-dependent models, as well as for detection of outliers or rare events within a time series of images. Furthermore, we demonstrate that the spectrum of PCA eigenvalues is a diagnostic of the power spectrum of the structure and, hence, of the underlying physical processes in the simulated and observed images.

Gravitational waves from the coalescence of two neutron stars were recently detected for the first time by the LIGO-Virgo collaboration, in event GW170817. This detection placed an upper limit on the effective tidal deformability of the two neutron stars and tightly constrained the chirp mass of the system. We report here on a new simplification that arises in the effective tidal deformability of the binary, when the chirp mass is specified. We find that, in this case, the effective tidal deformability of the binary is surprisingly independent of the component masses of the individual neutron stars, and instead depends primarily on the ratio of the chirp mass to the neutron star radius. Thus, a measurement of the effective tidal deformability can be used to directly measure the neutron star radius. We find that the upper limit on the effective tidal deformability from GW170817 implies that the radius cannot be larger than ~13km, at the 90% level, independent of the assumed masses for the component stars. The result can be applied generally, to probe the stellar radii in any neutron star-neutron star merger with a measured chirp mass. The approximate mass-independence disappears for neutron star-black hole mergers. Finally, we discuss a Bayesian inference of the equation of state that uses the measured chirp mass and tidal deformability from GW170817 combined with nuclear and astrophysical priors and discuss possible statistical biases in this inference.

Non-thermal electron acceleration via magnetic reconnection is thought to play an important role in powering the variable X-ray emission from radiatively inefficient accretion flows around black holes. The trans-relativistic regime of magnetic reconnection, where the magnetization $\sigma$, defined as the ratio of magnetic energy density to enthalpy density, is $\sim 1$, is frequently encountered in such flows. By means of a large suite of two-dimensional particle-in-cell simulations, we investigate electron and proton acceleration in the trans-relativistic regime. We focus on the dependence of the electron energy spectrum on $\sigma$ and the proton $\beta$ (i.e., the ratio of proton thermal pressure to magnetic pressure). We find that the electron spectrum in the reconnection region is non-thermal and can be generally modeled as a power law. At $\beta \lesssim 3 \times 10^{-3}$, the slope, $p$, is independent of $\beta$ and it hardens with increasing $\sigma$ as $p\simeq 1.8 +0.7/\sqrt{\sigma}$. Electrons are primarily accelerated by the non-ideal electric field at X-points, either in the initial current layer or in current sheets generated in between merging magnetic islands. At higher values of $\beta$, the electron power law steepens for all values of $\sigma$. At values of $\beta$ near $\beta_{\rm max}\approx1/4\sigma$, when both electrons and protons are relativistically hot prior to reconnection, the spectra of both species display an additional component at high energies, containing a few percent of particles. These particles are accelerated via a Fermi-like process by bouncing in between the reconnection outflow and a stationary magnetic island. We provide an empirical prescription for the dependence of the power-law slope and the acceleration efficiency on $\beta$ and $\sigma$, which can be used in global simulations of collisionless accretion disks.

The mass distribution of compact objects provides a fossil record that can be studied to uncover information on the late stages of massive star evolution, the supernova explosion mechanism, and the dense matter equation of state. Observations of neutron star masses indicate a bimodal Gaussian distribution, while the observed black hole mass distribution decays exponentially for stellar-mass black holes. We use these observed distributions to directly confront the predictions of stellar evolution models and the neutrino-driven supernova simulations of Sukhbold et al. (2016). We find excellent agreement between the black hole and low-mass neutron star distributions created by these simulations and the observations. We show that a large fraction of the stellar envelope must be ejected, either during the formation of stellar-mass black holes or prior to the implosion through tidal stripping due to a binary companion, in order to reproduce the observed black hole mass distribution. We also determine the origins of the bimodal peaks of the neutron star mass distribution, finding that the low-mass peak (centered at ~1.4 M_sun) originates from progenitors with M_zams ~ 9-18 M_sun. The simulations fail to reproduce the observed peak of high-mass neutron stars (centered at ~1.8 M_sun) and we explore several possible explanations. We argue that the close agreement between the observed and predicted black hole and low-mass neutron star mass distributions provides new promising evidence that these stellar evolution and explosion models are accurately capturing the relevant stellar, nuclear, and explosion physics involved in the formation of compact objects.

Non-ideal MHD effects may play a significant role in determining the dynamics, thermal properties, and observational signatures of radiatively inefficient accretion flows onto black holes. In particular, particle acceleration during magnetic reconnection events may influence black hole spectra and flaring properties. We use representative GRMHD simulations of black hole accretion flows to identify and explore the structures and properties of current sheets as potential sites of magnetic reconnection. In the case of standard and normal (SANE) disks, we find that, in the reconnection sites, the plasma beta ranges from $0.1$ to $1000$, the magnetization ranges from $10^{-4}$ to $1$, and the guide fields are weak compared to the reconnecting fields. In magnetically arrested (MAD) disks, we find typical values for plasma beta from $10^{-2}$ to $10^3$, magnetizations from $10^{-3}$ to $10$, and typically stronger guide fields, with strengths comparable to or greater than the reconnecting fields. These are critical parameters that govern the electron energy distribution resulting from magnetic reconnection and can be used in the context of plasma simulations to provide microphysics inputs to global simulations. We also find that ample magnetic energy is available in the reconnection regions to power the fluence of bright X-ray flares observed from the black hole in the center of the Milky Way.

One of the key goals of observing neutron stars is to infer the equation of state (EoS) of the cold, ultradense matter in their interiors. We present here a Bayesian statistical method of inferring the pressures at five fixed densities, from a sample of mock neutron star masses and radii. We show that while five polytropic segments are needed for maximum flexibility in the absence of any prior knowledge of the EoS, regularizers are also necessary to ensure that simple underlying EoS are not over-parametrized. For ideal data with small measurement uncertainties, we show that the pressure at roughly twice the nuclear saturation density, rho_sat, can be inferred to within 0.3 dex for many realizations of potential sources of uncertainties. The pressures of more complicated EoS with significant phase transitions can also be inferred to within ~30%. We also find that marginalizing the multi-dimensional parameter space of pressure to infer a mass-radius relation can lead to biases of nearly 1 km in radius, towards larger radii. Using the full, five-dimensional posterior likelihoods avoids this bias.

The Event Horizon Telescope (EHT) is a millimeter-wavelength, very-long-baseline interferometry (VLBI) experiment that is capable of observing black holes with horizon-scale resolution. Early observations have revealed variable horizon-scale emission in the Galactic Center black hole, Sagittarius A* (Sgr A*). Comparing such observations to time-dependent general relativistic magnetohydrodynamic (GRMHD) simulations requires statistical tools that explicitly consider the variability in both the data and the models. We develop here a Bayesian method to compare time-resolved simulation images to variable VLBI data, in order to infer model parameters and perform model comparisons. We use mock EHT data based on GRMHD simulations to explore the robustness of this Bayesian method and contrast it to approaches that do not consider the effects of variability. We find that time-independent models lead to offset values of the inferred parameters with artificially reduced uncertainties. Moreover, neglecting the variability in the data and the models often leads to erroneous model selections. We finally apply our method to the early EHT data on Sgr A*.

The increasing number and precision of measurements of neutron star masses, radii, and, in the near future, moments of inertia offer the possibility of precisely determining the neutron star equation of state. One way to facilitate the mapping of observables to the equation of state is through a parametrization of the latter. We present here a generic method for optimizing the parametrization of any physically allowed EoS. We use mock equations of state that incorporate physically diverse and extreme behavior to test how well our parametrization reproduces the global properties of the stars, by minimizing the errors in the observables mass, radius, and the moment of inertia. We find that using piecewise polytropes and sampling the EoS with five fiducial densities between ~1-8 times the nuclear saturation density results in optimal errors for the smallest number of parameters. Specifically, it recreates the radii of the assumed EoS to within less than 0.5 km for the extreme mock equations of state and to within less than 0.12 km for 95% of a sample of 42 proposed, physically-motivated equations of state. Such a parametrization is also able to reproduce the maximum mass to within 0.04 M_sun and the moment of inertia of a 1.338 M_sun neutron star to within less than 10% for 95% of the proposed sample of equations of state.

The observable quantities that carry the most information regarding the structures of the images of black holes in the interferometric observations with the Event Horizon Telescope are the closure phases along different baseline triangles. We use long time span, high cadence, GRMHD+radiative transfer models of Sgr A$^*$ to investigate the expected variability of closure phases in such observations. We find that, in general, closure phases along small baseline triangles show little variability, except in the cases when one of the triangle vertices crosses one of a small regions of low visibility amplitude. The closure phase variability increases with the size of the baseline triangle, as larger baselines probe the small-scale structures of the images, which are highly variable. On average, the jet-dominated MAD models show less closure phase variability than the disk-dominated SANE models, even in the large baseline triangles, because the images from the latter are more sensitive to the turbulence in the accretion flow. Our results suggest that image reconstruction techniques need to explicitly take into account the closure phase variability, especially if the quality and quantity of data allow for a detailed characterization of the nature of variability. This also implies that, if image reconstruction techniques that rely on the assumption of a static image are utilized, regions of the $u-v$ space that show a high level of variability will need to be identified and excised.

We investigate black hole-host galaxy scaling relations in cosmological simulations with a self-consistent black hole growth and feedback model. The sub-grid accretion model captures the key scalings governing angular momentum transport from galactic scales down to parsec scales, while our kinetic feedback implementation enables the injection of outflows with properties chosen to match observed nuclear outflows. We show that "quasar mode" feedback can have a large impact on the thermal properties of the intergalactic medium and the growth of galaxies and massive black holes for kinetic feedback efficiencies as low as 0.1% relative to the bolometric luminosity. Nonetheless, our simulations suggest that the black hole-host scaling relations are only weakly dependent on the effects of black hole feedback on galactic scales, owing to feedback suppressing the growth of galaxies and massive black holes by a similar amount. In contrast, the rate at which gravitational torques feed the central black hole relative to the host galaxy star formation rate governs the slope and normalization of the black hole-host correlations. Our results suggest that a common gas supply regulated by gravitational torques is the primary driver of the observed co-evolution of black holes and galaxies.

We present Chandra ACIS-S sub-array observations of the quiescent neutron star low-mass X-ray binaries X7 and X5 in the globular cluster 47 Tuc. The large reduction in photon pile-up compared to previous deep exposures enables a substantial improvement in the spectroscopic determination of the neutron star radius and mass of these neutron stars. Modeling the thermal emission from the neutron star surface with a non-magnetized hydrogen atmosphere and accounting for numerous sources of uncertainties, we obtain for the neutron star in X7 a radius of $R=11.1^{+0.8}_{-0.7}$ km for an assumed stellar mass of $M=1.4$ M$_{\odot}$ (68% C.L.). We argue, based on astrophysical grounds, that the presence of a He atmosphere is unlikely for this source. Due to eclipses and variable absorption, the quiescent low-mass X-ray binary X5 provides less stringent constraints, leading to a radius of $R=9.6^{+0.9}_{-1.1}$ km, assuming a hydrogen atmosphere and a mass of $M=1.4$ M$_{\odot}$. When combined with all other existing spectroscopic radius measurements, these measurements strongly favor radii in the 9.9-11.2 km range for a ~1.5 M$_{\odot}$ neutron star and point to a dense matter equation of state that is somewhat softer than the nucleonic ones that are consistent with laboratory experiments at low densities.

The Galactic Center black hole Sagittarius A* (Sgr A*) is a prime observing target for the Event Horizon Telescope (EHT), which can resolve the 1.3 mm emission from this source on angular scales comparable to that of the general relativistic shadow. Previous EHT observations have used visibility amplitudes to infer the morphology of the millimeter-wavelength emission. Potentially much richer source information is contained in the phases. We report on 1.3 mm phase information on Sgr A* obtained with the EHT on a total of 13 observing nights over 4 years. Closure phases, the sum of visibility phases along a closed triangle of interferometer baselines, are used because they are robust against phase corruptions introduced by instrumentation and the rapidly variable atmosphere. The median closure phase on a triangle including telescopes in California, Hawaii, and Arizona is nonzero. This result conclusively demonstrates that the millimeter emission is asymmetric on scales of a few Schwarzschild radii and can be used to break 180-degree rotational ambiguities inherent from amplitude data alone. The stability of the sign of the closure phase over most observing nights indicates persistent asymmetry in the image of Sgr A* that is not obscured by refraction due to interstellar electrons along the line of sight.

One of the primary science goals of the next generation of hard X-ray timing instruments is to determine the equation of state of the matter at supranuclear densities inside neutron stars, by measuring the radius of neutron stars with different masses to accuracies of a few percent. Three main techniques can be used to achieve this goal. The first involves waveform modelling. The flux we observe from a hotspot on the neutron star surface offset from the rotational pole will be modulated by the star's rotation, giving rise to a pulsation. Information about mass and radius is encoded into the pulse profile via relativistic effects, and tight constraints on mass and radius can be obtained. The second technique involves characterising the spin distribution of accreting neutron stars. The most rapidly rotating stars provide a very clean constraint, since the mass-shedding limit is a function of mass and radius. However the overall spin distribution also provides a guide to the torque mechanisms in operation and the moment of inertia, both of which can depend sensitively on dense matter physics. The third technique is to search for quasi-periodic oscillations in X-ray flux associated with global seismic vibrations of magnetars (the most highly magnetized neutron stars), triggered by magnetic explosions. The vibrational frequencies depend on stellar parameters including the dense matter equation of state. We illustrate how these complementary X-ray timing techniques can be used to constrain the dense matter equation of state, and discuss the results that might be expected from a 10m$^2$ instrument. We also discuss how the results from such a facility would compare to other astronomical investigations of neutron star properties. [Modified for arXiv]

We combine cosmological hydrodynamic simulations with analytic models to evaluate the role of galaxy-scale gravitational torques on the evolution of massive black holes at the centers of star-forming galaxies. We confirm and extend our earlier results to show that torque-limited growth yields black holes and host galaxies evolving on average along the Mbh-Mbulge relation from early times down to z = 0 and that convergence onto the scaling relation occurs independent of the initial conditions and with no need for mass averaging through mergers or additional self-regulation processes. Smooth accretion dominates the long-term evolution, with black hole mergers with mass ratios >1:5 representing typically a small fraction of the total growth. Winds from the accretion disk are required to eject significant mass to suppress black hole growth, but there is no need for coupling this wind to galactic-scale gas to regulate black holes in a non-linear feedback loop. Torque-limited growth yields a close-to-linear relation for the star formation rate and the black hole accretion rate averaged over galaxy evolution time scales. However, the SFR-AGN connection has significant scatter owing to strong variability of black hole accretion at all resolved time scales. Eddington ratios can be described by a broad lognormal distribution with median value evolving roughly as (1 + z)^1.9, suggesting a main sequence for black hole growth similar to the cosmic evolution of specific SFRs. Our results offer an attractive scenario consistent with available observations in which cosmological gas infall and transport of angular momentum in the galaxy by gravitational instabilities regulate the long-term co-evolution of black holes and star-forming galaxies.

Recent models of black hole growth in a cosmological context have forwarded a paradigm in which the growth is self-regulated by feedback from the black hole itself. Here we use cosmological zoom simulations of galaxy formation down to z = 2 to show that such strong self-regulation is required in the popular spherical Bondi accretion model, but that a plausible alternative model in which black hole growth is limited by galaxy-scale torques does not require self-regulation. Instead, this torque-limited accretion model yields black holes and galaxies evolving on average along the observed scaling relations by relying only on a fixed, 5% mass retention rate onto the black hole from the radius at which the accretion flow is fed. Feedback from the black hole may (and likely does) occur, but does not need to couple to galaxy-scale gas in order to regulate black hole growth. We show that this result is insensitive to variations in the initial black hole mass, stellar feedback, or other implementation details. The torque-limited model allows for high accretion rates at very early epochs (unlike the Bondi case), which if viable can help explain the rapid early growth of black holes, while by z = 2 it yields Eddington factors of 1%-10%. This model also yields a less direct correspondence between major merger events and rapid phases of black hole growth. Instead, growth is more closely tied to cosmological disk feeding, which may help explain observational studies showing that, at least at z > 1, active galaxies do not preferentially show merger signatures.

We use high-resolution cosmological zoom simulations with ~200 pc resolution at z = 2 and various prescriptions for galactic outflows in order to explore the impact of winds on the morphological, dynamical, and structural properties of eight individual galaxies with halo masses ~ 10^11--2x10^12 Msun at z = 2. We present a detailed comparison to spatially and spectrally resolved H{\alpha} and other observations of z ~ 2 galaxies. We find that simulations without winds produce massive, compact galaxies with low gas fractions, super-solar metallicities, high bulge fractions, and much of the star formation concentrated within the inner kpc. Strong winds are required to maintain high gas fractions, redistribute star-forming gas over larger scales, and increase the velocity dispersion of simulated galaxies, more in agreement with the large, extended, turbulent disks typical of high-redshift star-forming galaxies. Winds also suppress early star formation to produce high-redshift cosmic star formation efficiencies in better agreement with observations. Sizes, rotation velocities, and velocity dispersions all scale with stellar mass in accord with observations. Our simulations produce a diversity of morphological characteristics - among our three most massive galaxies, we find a quiescent grand-design spiral, a very compact star-forming galaxy, and a clumpy disk undergoing a minor merger; the clumps are evident in H{\alpha} but not in the stars. Rotation curves are generally slowly rising, particularly when calculated using azimuthal velocities rather than enclosed mass. Our results are broadly resolution-converged. These results show that cosmological simulations including outflows can produce disk galaxies similar to those observed during the peak epoch of cosmic galaxy growth.

We use a radiation hydrodynamic simulation of the hydrogen reionization epoch to study OI absorbers at z~6. The intergalactic medium (IGM) is reionized before it is enriched, hence OI absorption originates within dark matter halos. The predicted abundance of OI absorbers is in reasonable agreement with observations. At z=10, roughly 70% of sightlines through atomically-cooled halos encounter a visible (N_OI > 10^14 cm^-2) column. Reionization ionizes and removes gas from halos less massive than 10^8.4 M_0, but 20% of sightlines through more massive halos encounter visible columns even at z=5. The mass scale of absorber host halos is 10-100 times smaller than the halos of Lyman break galaxies and Lyman-alpha emitters, hence absorption probes the dominant ionizing sources more directly. OI absorbers have neutral hydrogen columns of 10^19-10^21 cm^-2, suggesting a close resemblance between objects selected in OI and HI absorption. Finally, the absorption in the foreground of the z=7.085 quasar ULASJ1120+0641 cannot originate in a dark matter halo because halo gas at the observed HI column density is enriched enough to violate the upper limits on the OI column. By contrast, gas at less than one third the cosmic mean density satisfies the constraints. Hence the foreground absorption likely originates in the IGM.

We use a suite of cosmological hydrodynamic simulations including a self-consistent treatment for inhomogeneous reionisation to study the impact of galactic outflows and photoionisation heating on the volume-averaged recombination rate of the intergalactic medium (IGM). By incorporating an evolving ionising escape fraction and a treatment for self-shielding within Lyman limit systems, we have run the first simulations of "photon-starved" reionisation scenarios that simultaneously reproduce observations of the abundance of galaxies, the optical depth to electron scattering of cosmic microwave background photons \tau, and the effective optical depth to Lyman\alpha absorption at z=5. We confirm that an ionising background reduces the clumping factor C by more than 50% by smoothing moderately-overdense (\Delta=1--100) regions. Meanwhile, outflows increase clumping only modestly. The clumping factor of ionised gas is much lower than the overall baryonic clumping factor because the most overdense gas is self-shielded. Photoionisation heating further suppresses recombinations if reionisation heats gas above the canonical 10,000 K. Accounting for both effects within our most realistic simulation, C rises from <1 at z>10 to 3.3 at z=6. We show that incorporating temperature- and ionisation-corrected clumping factors into an analytical reionisation model reproduces the numerical simulation's \tau to within 10%. Finally, we explore how many ionising photons are absorbed during the process of heating filaments by considering the overall photon cost of reionisation in analytical models that assume that the IGM is heated at different redshifts. For reionisation redshifts of 9--10, cold filaments boost the reionisation photon budget by ~1 photon per hydrogen atom.

We explore the evolution of the mass distribution of dust in collision-dominated debris disks, using the collisional code introduced in our previous paper. We analyze the equilibrium distribution and its dependence on model parameters by evolving over 100 models to 10 Gyr. With our numerical models, we confirm that systems reach collisional equilibrium with a mass distribution that is steeper than the traditional solution by Dohnanyi (1969). Our model yields a quasi steady-state slope of n(m) ~ m^{-1.88} [n(a) ~ a^{-3.65}] as a robust solution for a wide range of possible model parameters. We also show that a simple power-law function can be an appropriate approximation for the mass distribution of particles in certain regimes. The steeper solution has observable effects in the submillimeter and millimeter wavelength regimes of the electromagnetic spectrum. We assemble data for nine debris disks that have been observed at these wavelengths and, using a simplified absorption efficiency model, show that the predicted slope of the particle mass distribution generates SEDs that are in agreement with the observed ones.

We carry out a new suite of cosmological radiation hydrodynamic simulations and explore the relative impacts on reionization-epoch star formation of galactic outflows and photoionization heating. By itself, an extragalactic ultraviolet background (EUVB) suppresses the luminosity function by less than 50% at z=6, overproducing the observed galaxy abundance by a factor of 3-5. Galactic outflows restore agreement with observations without preventing Population II star formation from reionizing the Universe by z=6. The resulting EUVB suppresses star formation in halos with virial temperatures below 10^5K but has a weaker impact in more massive halos. Nonetheless, the low-mass halos contribute up to 50% of all ionizing photons owing to the EUVB's inhomogeneity. Overall, star formation rate scales as halo mass M_h to the 1.3-1.4 in halos with $M_h=10^{8.2--10.2}\msun$. This is a steeper dependence than is often assumed in reionization models, boosting the expected power spectrum of 21 centimeter fluctuations on large scales. The luminosity function rises steeply to at least M_1600=-13, indicating that reionization was driven by faint galaxies (M_1600 >= -15) that have not yet been observed. Our models cannot simultaneously explain observations of galaxies, the cosmic microwave background, and the intergalactic medium. Increased dynamic range will alleviate the existing discrepancies, but observations may still require additional physics such as a variable ionizing escape fraction (abridged).