We present two prescriptions for broadband (~77 - 252 GHz), millimeter-wave antireflection coatings for cryogenic, sintered polycrystalline aluminum oxide optics: one for large-format (700 mm diameter) planar and plano-convex elements, the other for densely packed arrays of quasi-optical elements, in our case 5 mm diameter half-spheres (called "lenslets"). The coatings comprise three layers of commercially-available, polytetrafluoroethylene-based, dielectric sheet material. The lenslet coating is molded to fit the 150 mm diameter arrays directly while the large-diameter lenses are coated using a tiled approach. We review the fabrication processes for both prescriptions then discuss laboratory measurements of their transmittance and reflectance. In addition, we present the inferred refractive indices and loss tangents for the coating materials and the aluminum oxide substrate. We find that at 150 GHz and 300 K the large-format coating sample achieves (97 +/- 2)% transmittance and the lenslet coating sample achieves (94 +/- 3)% transmittance.

The cosmic microwave background (CMB) encodes information about the content and evolution of the universe. The presence of light, weakly interacting particles impacts the expansion history of the early universe, which alters the temperature and polarization anisotropies of the CMB. In this way, current measurements of the CMB place interesting constraints on the neutrino energy density and mass, as well as on the abundance of other possible light relativistic particle species. We present the status of an on-going 1500 sq. deg. survey with the SPT-3G receiver, a new mm-wavelength camera on the 10-m diameter South Pole Telescope (SPT). The SPT-3G camera consists of 16,000 superconducting transition edge sensors, a 10x increase over the previous generation camera, which allows it to map the CMB with an unprecedented combination of sensitivity and angular resolution. We highlight projected constraints on the abundance of sterile neutrinos and the sum of the neutrino masses for the SPT-3G survey, which could help determine the neutrino mass hierarchy.

SPT-3G is a polarization-sensitive receiver, installed on the South Pole Telescope, that measures the anisotropy of the cosmic microwave background (CMB) from degree to arcminute scales. The receiver consists of ten 150~mm-diameter detector wafers, containing a total of 16,000 transition-edge sensor (TES) bolometers observing at 95, 150, and 220 GHz. During the 2018-2019 austral summer, one of these detector wafers was replaced by a new wafer fabricated with Al-Mn TESs instead of the Ti/Au design originally deployed for SPT-3G. We present the results of in-lab characterization and on-sky performance of this Al-Mn wafer, including electrical and thermal properties, optical efficiency measurements, and noise-equivalent temperature. In addition, we discuss and account for several calibration-related systematic errors that affect measurements made using frequency-domain multiplexing readout electronics.

Frequency-domain multiplexing (fMux) is an established technique for the readout of large arrays of transition edge sensor (TES) bolometers. Each TES in a multiplexing module has a unique AC voltage bias that is selected by a resonant filter. This scheme enables the operation and readout of multiple bolometers on a single pair of wires, reducing thermal loading onto sub-Kelvin stages. The current receiver on the South Pole Telescope, SPT-3G, uses a 68x fMux system to operate its large-format camera of $\sim$16,000 TES bolometers. We present here the successful implementation and performance of the SPT-3G readout as measured on-sky. Characterization of the noise reveals a median pair-differenced 1/f knee frequency of 33 mHz, indicating that low-frequency noise in the readout will not limit SPT-3G's measurements of sky power on large angular scales. Measurements also show that the median readout white noise level in each of the SPT-3G observing bands is below the expectation for photon noise, demonstrating that SPT-3G is operating in the photon-noise-dominated regime.

We use hydrodynamic cosmological simulations in a (600 Mpc)^3 volume to study the observability of baryon acoustic oscillations (BAO) in the intergalactic medium as probed by Lyman alpha forest (LAF) absorption. The large scale separation between the wavelength of the BAO mode (~150 Mpc) and the size of LAF absorbers (~100 kpc) makes this a numerically challenging problem. We report on several 2048^3 simulations of the LAF using the ENZO code. We adopt WMAP5 concordance cosmological parameters and power spectrum including BAO perturbations. 5000 synthetic HI absorption line spectra are generated randomly piercing the box face. We calculate the cross-correlation function between widely separated pairs. We detect the BAO signal at z=3 where theory predicts to moderate statistical significance.

We consider a PDE system comprising compressible hydrodynamics, flux-limited diffusion radiation transport and chemical ionization kinetics in a cosmologically-expanding universe. Under an operator-split framework, the cosmological hydrodynamics equations are solved through the Piecewise Parabolic Method, as implemented in the Enzo community hydrodynamics code. The remainder of the model, including radiation transport, chemical ionization kinetics, and gas energy feedback, form a stiff coupled PDE system, which we solve using a fully-implicit inexact Newton approach, and which forms the crux of this paper. The inner linear Newton systems are solved using a Schur complement formulation, and employ a multigrid-preconditioned conjugate gradient solver for the inner Schur systems. We describe this approach and provide results on a suite of test problems, demonstrating its accuracy, robustness, and scalability to very large problems.

We give a comprehensive statistical description of the Lyman-alpha absorption from the intergalactic medium in a hydrodynamic simulation at redshifts 0.1-1.6, the range of redshifts covered by HST spectra of QSOs. We use the ENZO code to make a 76 comoving Mpc cube simulation using 75 kpc cells, for a Hubble constant of 71 km/s/Mpc. The best prior work, by \citet{dave99},used an SPH simulation in a 15.6 Mpc box with an effective resolution of 245 kpc and slightly different cosmological parameters. At redshifts z=2 this simulation is different from data. \citet{tytler07b} found that the simulated spectra at z=2 have too little power on large scales, Lyman-alpha lines are too wide, there is a lack high column density lines, and there is a lack of pixels with low flux. Here we present statistics at z<1.6, including the flux distribution, the mean flux, the effective opacity, and the power and correlation of the flux. We also give statistics of the lyman alpha lines including the line width distribution, the column density distribution, the number of lines per unit equivalent width and redshift, and the correlation between the line width and column density. We find that the mean amount of absorption in the simulated spectra changes smoothly with redshift with DA(z)=0.01(1+z)^{2.25}. Both the trend and absolute values are close to measurements of HST spectra by \citet{kirkman07a}. The column density and line width distributions are also close to those measured from HST spectra by \citet{janknecht06a}, except for the mode of the line width distribution which is smaller in the HST spectra. Although some differences that we saw at z=2 are too subtle to be seen in existing HST spectra, overall, the simulation gives an good description of HST spectra at 0.1<z<1.6.

We study the effects of box size on ENZO simulations of the intergalactic medium (IGM) at z = 2. We follow statistics of the cold dark matter (CDM) and the Lya absorption. We find that the larger boxes have fewer pixels with significant absorption (flux < 0.96) and more pixels in longer stretches with little or no absorption, and they have wider Lya lines. We trace these effect back to the additional power in larger boxes from longer wavelength modes. The IGM in our larger boxes is hotter, from increased pressure heating due to faster hydrodynamical infall. When we increase the photoheating in smaller boxes to compensate, their Lya statistics change to mimic those of a box of twice the size. Statistics converge towards their value in the largest (76.8 Mpc) box, except for the most common value of the CDM density which continues to rise. When we compare to errors with data, we find that our 76.8 Mpc box is larger than we need for the mean flux, barely large enough for the column density distribution and the power spectrum of the flux, and too small for the line widths. This box with 75 kpc cells has approximately the same mean flux as QSO spectra, but the Lya lines are too wide by 2.6 km/s, there are too few lines with log H I column densities > 10^17 cm^-2, and the power of the flux is too low by 20 - 50%, from small to large scales. Four times smaller cell size does not resolve these differences, nor do simple changes to the ultraviolet background that drives the H and He II ionization. It is hard to see how simulations using popular cosmological and astrophysical parameters can match Lyman-alpha forest data at z=2.

We present the first large sample of absorption systems in paired QSOs consisting of 691 absorption systems in the spectra of 310 QSOs including 170 pairings. All these absorption systems have metal lines, usually C IV or Mg II. We see 17 cases of absorption in one line-of-sight within 200 km/s (1 Mpc) of absorption in the paired line-of-sight with the probability at least approx 50% at 100kpc, declining rapidly to 23% at 100 - 200 kpc. We detect clustering on 0.5Mpc scales and see a hint of the "fingers of God" redshift-space distortion. The distribution matches absorbers arising in galaxies at z=2 with a normal correlation function and systematic infall velocities but unusually low random pair-wise velocity differences. Absorption in gas flowing out from galaxies at a mean velocity of 250 km/s would produce vastly more elongation than we see. The UV absorption from fast winds that Adelberger et al. 2005 see in spectra of LBGs is not representative of the absorption that we see. Either the winds are confined to LBGs, or they can not extend to 40 kpc with large velocities, while continuing to make UV absorption we see, implying most metals were in place in the IGM long before z=2. Separately, when we examine the absorption seen when a sight line passes a second QSO, we see 19 absorbers within 400 km/s of the partner QSO. The probability of seeing absorption is approximately constant for impact parameters 0.1 - 1.5 Mpc. Perhaps we do not see a rapid rise in the probability at small impact parameters because the UV from QSOs destroys some absorbers near to the QSOs. The 3D distribution of 64 absorbers around 313 QSOs is to first order isotropic, with just a hint of the anisotropy expected if the QSO UV emission is beamed, or alternatively QSOs might emit UV isotropically but for a surprisingly short time of only 0.3Myr.

We study the ionization and thermal evolution of the intergalactic medium during the epoch of \heii reionization by means of radiation hydrodynamical cosmological simulations. We post-process baryonic density fields from a standard optically-thin IGM simulation with a homogeneous galaxy-dominated UV background (UVB) which reionizes \hi and \hei at z=6.5 but does not have any contribution to the ionization of \heii. Quasars with luminosities proportional to the mass of the host halos are then introduced as point sources throughout the 100 Mpc simulation volume consistent with the Pei luminosity function. We evolve the spatial distribution of the \heii ionizing radiation field using a time-implicit variable tensor Eddington factor radiative transfer scheme. Simultaneously, we also solve for the local ionization of \heii to \heii and the associated photoheating of the gas. We find that the percolation of the \heiii regions is essentially complete by z=2.5. When comparing to a self-consistent optically thin simulation we find that in optically thick calculation the gas temperature is higher by a factor of approximately 1.7 at the mean gas density level. We use 300 random lines of sight to compute at $\bar{z} = 2.5 \pm 0.1$ a mean \heii \lya line transmission of $\bar{F} = 0.304 \pm 0.002$. We compare the broadening width of the \hi and \heii \lya lines to the results from the self-consistent optically thin simulation and find a shift by approximately 1.25 km/s of the b-parameter distribution. Estimating the relative broadening width between the two forests shows that the \heii median b-parameter is about 0.8 times the median \hi broadening width. This implies that the \heii absorbers are physically extended consistent with conclusions from observed lines of sight.

We use high quality echelle spectra of 24 QSOs to provide a calibrated measurement of the total amount of Lyman-alpha forest absorption (DA) over the redshift range 2.2 < z < 3.2. Our measurement of DA excludes absorption from metal lines or the Lyman-alpha lines of Lyman limit systems and damped Lyman-alpha systems. We use artificial spectra with realistic flux calibration errors to show that we are able to place continuum levels that are accurate to better than 1%. When we combine our results with our previous results between 1.6 < z < 2.2, we find that the redshift evolution of DA is well described over 1.6 < z < 3.2 as A*(1+z)^gamma, where A = 0.0062 and gamma = 2.75. We detect no significant deviations from a smooth power law evolution over the redshift range studied. We find less H I absorption than expected at z=3, implying that the UV background is about 40% higher than expected. Our data appears to be consistent with an H I ionization rate of Gamma ~ 1.4E-12 1/s.

We present 40 fully hydrodynamical numerical simulations of the intergalactic gas that gives rise to the Ly-alpha forest. We make artificial spectra from each and measure three output parameters: the mean flux, a measure of the most common Ly-alpha line width (b-value), and the 1D power spectrum of the flux on scales from 0.01 - 0.1 s/km. We investigate how these output parameters are correlated with the input parameters that describe the simulations. These input parameters include the main cosmological parameters; and two astrophysical parameters gammah and gammahe. The parameter gammah controls the rate of ionization of H I, He I and He II while gammahe controls the rate of heating. We compare the output from the simulations to data from HIRES and UVES spectra and we give a new measurement of the flux power from these spectra at z=1.95. We find that simulations with a wide variety of sigma_8 values, from at least 0.8 - 1.1, can fit the small scale flux power and b-values when we adjust gammahe to compensate for the sigma_8 change. We can also use gammah to adjust the H I ionization rate to simultaneously match the mean flux. We derive scaling relations that give the output parameter values expected for a variety of input parameters with errors of 2% for optical depth and 1.4% for b-balues. We can break the degeneracy between gammahe and sigma_8 using large scale power or other data to fix sigma_8.

We use hydrodynamic cosmological simulations in a 9.6 Mpc box to explore the evolution of the intergalactic medium (IGM) transmissivity from z=2 through the epoch of reionization. Reionization is achieved through an ultraviolet background (UVB) that includes evolving stellar and QSO source populations. We construct and analyze noiseless synthetic HI Lya absorption spectra along lines of sight through our continuously evolving box and find a smooth evolution of the effective optical depth under a power law up to the epoch of reionization. Crossing into the epoch of reionization, both the mean transmitted flux (MTF) and variance to the mean transmitted flux sharply deviate from a smooth evolution and hence such observables have large margins of error. Despite the statistical uncertainty in inferring the reionization profile from spectra, the end of an opacity phase transition of the IGM correlates well with the redshift when both the mean and variance of the transmitted flux rapidly deviate from the post-reionization profile. However, an unobtainable number of lines of sight is needed to allow an estimate of the MTF with less than 10% relative margin of error. In addition to optical transmission, we compare the predicted dark gap length distribution with observations. We show that this statistic is sensitive to spectral resolution at reionization redshifts, but overall in agreement with results by Songaila and Cowie (2002). Finally, we derive a positive correlation between the mean optical depth within a gap and the size of the gap, in attempt to relate "transmission statistics" to "dark gap statistics" in high redshift studies of the IGM.

Using H I absorption alone, we attempt to separate H I absorption lines in quasar spectra into two categories; HDLs (Higher Density Lines) and LDLs (Lower Density Lines), and we discuss the difference in their physical properties. We deblend and fit all H I lines with Voigt profiles, and make an unbiased sample of H I lines covering a wide column density range (12 < log N_HI < 19 cm^-2). To reduce the influence of line blending, we simultaneously fit several Lyman series lines. As a result of a two-point correlation analysis, we found that higher column density H I lines are clustering at $\Delta$ v < 200 km/s, while lower ones at $\Delta$ v < 100 km/s. We define HDLs as H I lines with 15 < log N_HI < 19 cm^-2 and all H I lines within +/- 200 km/s of a line with log N_HI > 15 cm^-2, and LDLs as others with 12 < log N_HI < 15 cm^-2. We found that the HDLs have smaller minimum b-values for a given column density than the LDLs. This difference is successfully reproduced by our Hydrodynamic simulation. The LDLs seem to be cool or shock-heated diffuse IGM gas, while the HDLs are likely to be cooler dense gas near to galaxies.

We identify a concordant model for the intergalactic medium (IGM) at redshift z=1.9 that uses popular values for cosmological and astrophysical parameters and accounts for all baryons with an uncertainty of 6%. We have measured the amount of absorption, DA, in the Ly-alpha forest at redshift 1.9 in spectra of 77 QSO from the Kast spectrograph. We calibrated the continuum fits with realistic artificial spectra, and we found that averaged over all 77 QSOs the mean continuum level is within 1-2% of the correct value. Absorption from all lines in the Ly-alpha forest at z=1.9 removes DA=15.1 +/- 0.7% of the flux between 1070 and 1170 (rest) Angstroms. This is the first measurement using many QSOs at this z, and the first calibrated measurement at any redshift. Metal lines absorb 2.3 +/- 0.5%, and LLS absorb 1.0 +/- 0.4% leaving 11.8 +/- 1.0% from the lower density bulk of the IGM. Averaging over Delta z=0.1 or 154 Mpc, the dispersion is 6.1 +/- 0.3% including LLS and metal lines, or 3.9 (+0.5, -0.7)% for the lower density IGM alone, consistent with the usual description of large scale structure. LLS and metal lines are major contributors to the variation in the mean flux, and they make the flux field significantly non-Gaussian. We find that a hydrodynamic simulation on a 1024 cubed grid in a 75.7 Mpc box reproduces the observed DA from the low density IGM with parameters values H_o=71 km/s/Mpc, Omega_Lambda=0.73, Omega_m=0.27, Omega_b=0.044, sigma_8=0.9 and a UV background that has an ionization rate that is 1.08 +/- 0.08 times the prediction by Madau, Haardt & Rees (1999).

We examine the dynamical evolution and statistical properties of the supernova ejecta of massive primordial stars in a cosmological framework to determine whether this first population of stars could have enriched the universe to the levels and dispersions seen by the most recent observations of the Lyman-Alpha forest. We evolve a lambda CDM model in a 1 Mpc^3 volume to a redshift of z = 15 and add ``bubbles'' of metal corresponding to the supernova ejecta of the first generation of massive stars in all dark matter halos with masses greater than 5 times 10^5 solar masses. These initial conditions are then evolved to z = 3 and the distribution and levels of metals are compared to observations. In the absence of further star formation the primordial metal is initially contained in halos and filaments. Photoevaporation of metal-enriched gas due to the metagalactic ultraviolet background radiation at the epoch of reionization (z ~ 6) causes a sharp increase of the metal volume filling factor. At z = 3, ~ 2.5% of the simulation volume (approx. 20% of the total gas mass) is filled with gas enriched above a metallicity of 10^-4 Z_solar, and less than 0.6% of the volume is enriched above a metallicity of 10^-3 Z_solar. This suggests that, even with the most optimistic prescription for placement of primordial supernova and the amount of metals produced by each supernova, this population of stars cannot entirely be responsible for the enrichment of the Lyman-$\alpha$ forest to the levels and dispersions seen by current observations unless we have severely underestimated the duration of the Pop III epoch. However, comparison to observations show that Pop III supernovae can be significant contributors to the very low overdensity Lyman-Alpha forest.

We describe an approach for incorporating radiative transfer into 3D hydrodynamic cosmological simulations. The method, while approximate, allows for a self-consistent treatment of self-shielding and shadowing, diffuse and point sources of radiation, and frequency dependent transfer. Applications include photodissociation, photoheating, and photoionization of the IGM.