Super-Earths and mini-Neptunes exhibit great diversity in their compositional and orbital properties. Their bulk densities span a large range, from those dense enough to be purely rocky to those needing a substantial contribution from volatiles to their volumes. Their orbital configurations range from compact, circular multi-transiting systems like Kepler-11 to systems like our Solar System's terrestrial planets with wider spacings and modest but significant eccentricities and mutual inclinations. Here we investigate whether a continuum of formation conditions resulting from variation in the amount of solids available in the inner disk can account for the diversity of orbital and compositional properties observed for super Earths, including the apparent dichotomy between single transiting and multiple transiting system. We simulate in situ formation of super-Earths via giant impacts and compare to the observed Kepler sample. We find that intrinsic variations among disks in the amount of solids available for in situ formation can account for the orbital and compositional diversity observed among Kepler's transiting planets. Our simulations can account for the planets' distributions of orbital period ratios, transit duration ratios, and transit multiplicity; higher eccentricities for single than multi transiting planets; smaller eccentricities for larger planets; scatter in the mass-radius relation, including lower densities for planets with masses measured with TTVs than RVs; and similarity in planets' sizes and spacings within each system. Our findings support the theory that variation among super-Earth and mini-Neptune properties is primarily locked in by different in situ formation conditions, rather than arising stochastically through subsequent evolution.

Resolved debris disk features (e.g., warps, offsets, edges and gaps, azimuthal asymmetries, radially thickened rings, scale heights) contain valuable information about the underlying planetary systems, such as the posited planet's mass, semi-major axis, and other orbital parameters. Most existing models assume a single planet is sculpting the disk feature, but recent observations of mature planetary systems (e.g., by radial velocity surveys or \textit{Kepler}) have revealed that many planets reside in multi-planet systems. Here we investigate if/how planet properties inferred from single-planet models are compromised when multiple planets reside in the system. For each disk feature, we build a two-planet model that includes a planet b with fixed parameters and a planet c with a full range of possible parameters. We investigate these two-planet systems and summarize the configurations for which assuming a single planet (i.e., planet b) leads to significantly flawed inferences of that planet's properties. We find that although disk features are usually primarily dominated by a single planet, when using single-planet models we are at risk of misinterpreting planet properties by orders of magnitude in extreme cases. Specifically, we are at high risk of misinterpreting planet properties from disk warps; at moderate risk from disk edges and gaps, radially thickened rings, and scale height features; and at low risk from host star-disk center offsets and azimuthal asymmetries. We summarize situations where we can infer the need to use a multi-planet model instead of a single-planet one from disk morphology dissimilarities.

We present SQUEzE, a software package to classify quasar spectra and estimate their redshifts. SQUEzE is a random forest classifier operating on the parameters of candidate emission peaks identified in the spectra. We test the performance of the algorithm using visually inspected data from BOSS as a truth table. Only 4\% of the sample ($\sim$6,800 quasars and $\sim$11,520 contaminants) is needed for converged training in recommended choices of the confidence threshold ($0.2<p_{\rm min}<0.7$). For an operational mode which balances purity and completeness ($p_{\rm min}=0.28$) we recover a purity of $96.81\pm0.39\%$ ($99.30\pm0.14\%$ for quasars with $z \geq 2.1$) and a completeness of $96.83\pm0.30\%$ ($98.42\pm0.15\%$ for quasars with $z \geq 2.1$). SQUEzE can be used to obtain a $\approx$100\% pure sample of $z \geq 2.1$. quasars (with $\approx$96\% completeness) by using a confidence threshold of $p_{\rm min}=0.7$. The estimated redshift error is $1,500{\rm \thinspace km/s}$ and we recommend that SQUEzE be used in conjunction with an additional step of redshift tuning to achieve maximum precision. We find that SQUEzE achieves the necessary performance to replace visual inspection in BOSS-like spectroscopic surveys of quasars with subsequent publications in this series exploring expectations for future surveys and alternative methods. Keywords: cosmology: observations - quasar: emission lines - quasar: absorption lines

The current disparity in computational knowledge is a critical hindrance to the diversity and success of the field. Recommendations are outlined for policies and funding models to enable the growth and retention of a new generation of computational researchers that reflect the demographics of the undergraduate population in Astronomy and Physics.

We present an analysis combining the simultaneous measurement of intergalactic absorption by hydrogen (H I), helium (He II) and oxygen (O VI) in UV and optical quasar spectra. The combination of the H I and He II Lyman-alpha forests through $\eta$ (the ratio of column densities of singly ionized helium to neutral hydrogen) is thought to be sensitive to large-scale inhomogeneities in the extragalactic UV background. We test this assertion by measuring associated five-times-ionized oxygen (O VI) absorption, which is also sensitive to the UV background. We apply the pixel optical depth technique to O VI absorption in high and low $\eta$ samples filtered on various scales. This filtering scale is intended to represent the dominant scale of any coherent oxygen excess/deficit. We find a $2\sigma$ detection of an O VI opacity excess in the low $\eta$ sample on scales of $\sim$10 cMpc for HE 2347-4342 at $\bar{z}\approx 2.6$, consistent with a large-scale excess in hard UV photons. However, for HS 1700+6416 at $\bar{z}\approx 2.5$ we find that the measured O VI absorption is not sensitive to differences in $\eta$. HS 1700+6416 also shows a relative absence of O VI overall, which is $6\sigma$ inconsistent with that of HE 2347-4342. This implies UV background inhomogeneities on $\gtrsim$200 cMpc scales, hard UV regions having internal ionization structure on $\sim$10 cMpc scales and soft UV regions showing no such structure. Furthermore, we perform the pixel optical depth search for oxygen on the He II Gunn-Peterson trough of HE 2347-4342 and find results consistent with post-He II-reionization conditions.

While many observed debris disks are thought to have gaps suggestive of the presence of planets, direct imaging surveys do not find many high mass planets in these systems. We investigate if divergent migration is a viable mechanism for forming gaps in young debris disks with planets of low enough mass to currently elude detection. We perform numerical integrations of planet pairs embedded in planetesimal disks to assess the conditions for which divergent, planetesimal-driven migration occurs and gaps form within the disk. Gap widths and the migration rate of planets within a pair depend on both disk mass and the degree to which the planets share disk material. We find that planet pairs with planets more massive than Neptune can produce gaps with widths similar to their orbit distance within 10 Myr at orbit separations probed by direct imaging campaigns. Pairs of migrating super-Earths likely cannot form observable gaps on the same time and distance scales, however. Inferring the responsible planet masses from these gaps while neglecting migration could overestimate the mass of planets by more than an order of magnitude.

We present abundance measurements of elements O, C, Si and Fe for three gas-rich galaxies at z~5 using observations from the Very Large Telescope (VLT) and the Keck telescope in order to better constrain the early chemical enrichment of gas-rich galaxies. These galaxies show strong Lyman-{\alpha} absorption in the spectra of background quasars, with neutral hydrogen column densities log N${_H{}_I}$(cm$^{-2}$)=20.10$\pm$0.15, 20.10$\pm$0.15, and 20.80$\pm$0.15. Using the undepleted element O, we find the metallicities [O/H] to be in the range of -2.51 to -2.05 dex. Our study has doubled the existing sample of measurements of undepleted elements at z>4.5. Combining our measurements with those from the literature, we find that the N${_H{}_I}$-weighted mean metallicity of z~5 absorbers is consistent with the prediction based on z<4.5 DLAs. Thus, we find no significant evidence of a sudden drop in metallicity at z>4.7 as reported by some prior studies. We also determine the extent of dust depletion using a combination of both the volatile element O and the refractory elements Si and/or Fe. Some of the absorbers show evidence of depletion of elements on dust grains, e.g. low [Si/O] or [Fe/O]. The relative abundances of these absorbers along with other z~5 absorbers from the literature show some peculiarities, e.g. low [C/O] in several absorbers and high [Si/O] in one absorber. Finally, we find that the metallicity vs. velocity dispersion relation of z~5 absorbers may be different from that of lower-redshift absorbers.

Planets and minor bodies such as asteroids, Kuiper-belt objects and comets are integral components of a planetary system. Interactions among them leave clues about the formation process of a planetary system. The signature of such interactions is most prominent through observations of its debris disk at millimeter wavelengths where emission is dominated by the population of large grains that stay close to their parent bodies. Here we present ALMA 1.3 mm observations of HD 95086, a young early-type star that hosts a directly imaged giant planet b and a massive debris disk with both asteroid- and Kuiper-belt analogs. The location of the Kuiper-belt analog is resolved for the first time. The system can be depicted as a broad ($\Delta R/R \sim$0.84), inclined (30\arcdeg$\pm$3\arcdeg) ring with millimeter emission peaked at 200$\pm$6 au from the star. The 1.3 mm disk emission is consistent with a broad disk with sharp boundaries from 106$\pm$6 to 320$\pm$20 au with a surface density distribution described by a power law with an index of --0.5$\pm$0.2. Our deep ALMA map also reveals a bright source located near the edge of the ring, whose brightness at 1.3 mm and potential spectral energy distribution are consistent with it being a luminous star-forming galaxy at high redshift. We set constraints on the orbital properties of planet b assuming co-planarity with the observed disk.

We present a new method for interferometric imaging that is ideal for the large fields of view and compact arrays common in 21 cm cosmology. We first demonstrate the method with simulations for two very different low frequency interferometers, the Murchison Widefield Array (MWA) and the MIT Epoch of Reionization (MITEoR) Experiment. We then apply the method to the MITEoR data set collected in July 2013 to obtain the first northern sky map from 128 MHz to 175 MHz at about 2 degree resolution, and find an overall spectral index of -2.73+/-0.11. The success of this imaging method bodes well for upcoming compact redundant low-frequency arrays such as HERA. Both the MITEoR interferometric data and the 150 MHz sky map are publicly available at http://space.mit.edu/home/tegmark/omniscope.html.

Element abundances in high-redshift quasar absorbers offer excellent probes of the chemical enrichment of distant galaxies, and can constrain models for population III and early population II stars. Recent observations indicate that the sub-damped Lyman-alpha (sub-DLA) absorbers are more metal-rich than DLA absorbers at redshifts 0$<$$z$$<$3. It has also been suggested that the DLA metallicity drops suddenly at $z$$>$4.7. However, only 3 DLAs at $z$$>$4.5 and none at $z$$>$3.5 have "dust-free" metallicity measurements of undepleted elements. We report the first quasar sub-DLA metallicity measurement at $z$$>$3.5, from detections of undepleted elements in high-resolution data for a sub-DLA at $z$=5.0. We obtain fairly robust abundances of C, O, Si, and Fe, using lines outside the Lyman-alpha forest. This absorber is metal-poor, with O/H]=-2.00$\pm$0.12, which is $\gtrsim$4$\sigma$ below the level expected from extrapolation of the trend for $z$$<$3.5 sub-DLAs. The C/O ratio is 1.8$^{+0.4}_{-0.3}$ times lower than in the Sun. More strikingly, Si/O is 3.2$^{+0.6}_{-0.5}$ times lower than in the Sun, while Si/Fe is nearly (1.2$^{+0.4}_{-0.3}$ times) solar. This absorber does not display a clear alpha/Fe enhancement. Dust depletion may have removed more Si from the gas phase than is common in the Milky Way interstellar medium, which may be expected if high-redshift supernovae form more silicate-rich dust. C/O and Si/O vary substantially between different velocity components, indicating spatial variations in dust depletion and/or early stellar nucleosynethesis (e.g., population III star initial mass function). The higher velocity gas may trace an outflow enriched by early stars.

In the coming years, high contrast imaging surveys are expected to reveal the characteristics of the population of wide-orbit, massive, exoplanets. To date, a handful of wide planetary mass companions are known, but only one such multi-planet system has been discovered: HR8799. For low mass planetary systems, multi-planet interactions play an important role in setting system architecture. In this paper, we explore the stability of these high mass, multi-planet systems. While empirical relationships exist that predict how system stability scales with planet spacing at low masses, we show that extrapolating to super-Jupiter masses can lead to up to an order of magnitude overestimate of stability for massive, tightly packed systems. We show that at both low and high planet masses, overlapping mean motion resonances trigger chaotic orbital evolution, which leads to system instability. We attribute some of the difference in behavior as a function of mass to the increasing importance of second order resonances at high planet-star mass ratios. We use our tailored high mass planet results to estimate the maximum number of planets that might reside in double component debris disk systems, whose gaps may indicate the presence of massive bodies.

We report Keck/ESI and VLT/UVES observations of three super-damped Lyman-alpha quasar absorbers with H I column densities log N(HI) >= 21.7 at redshifts z=2-2.5. All three absorbers show similar metallicities (-1.3 to -1.5 dex), and dust depletion of Fe, Ni, and Mn. Two of the absorbers show supersolar [S/Zn] and [Si/Zn]. We combine our results with those for other DLAs to examine trends between N(HI), metallicity, dust depletion. A larger fraction of the super-DLAs lie close to or above the line [X/H]=20.59-log N(HI) in the metallicity vs. N(HI) plot, compared to the less gas-rich DLAs, suggesting that super-DLAs are more likely to be rich in molecules. Unfortunately, our data for Q0230-0334 and Q0743+1421 do not cover H2 absorption lines. For Q1418+0718, some H2 lines are covered, but not detected. CO is not detected in any of our absorbers. For DLAs with log N(HI) < 21.7, we confirm strong correlation between metallicity and Fe depletion, and find a correlation between metallicity and Si depletion. For super-DLAs, these correlations are weaker or absent. The absorbers toward Q0230-0334 and Q1418+0718 show potential detections of weak Ly-alpha emission, implying star formation rates of about 1.6 and 0.7 solar masses per year, respectively (ignoring dust extinction). Upper limits on the electron densities from C II*/C II or Si II*/Si II are low, but are higher than the median values in less gas-rich DLAs. Finally, systems with log N(HI) > 21.7 may have somewhat narrower velocity dispersions delta v_90 than the less gas-rich DLAs, and may arise in cooler and/or less turbulent gas.

HD 95086 is a young early-type star that hosts (1) a 5 MJ planet at the projected distance of 56 AU revealed by direct imaging, and (2) a prominent debris disk. Here we report the detection of 69 um crystalline olivine feature from the disk using the Spitzer/MIPS-SED data covering 55-95 um. Due to the low resolution of MIPS-SED mode, this feature is not spectrally resolved, but is consistent with the emission from crystalline forsterite contributing 5% of the total dust mass. We also present detailed analysis of the disk SED and re-analysis of resolved images obtained by Herschel. Our results suggest that the debris structure around HD 95086 consists of a warm (175 K) belt, a cold (55 K) disk, and an extended disk halo (up to 800 AU), and is very similar to that of HR 8799. We compare the properties of the three debris components, and suggest that HD 95086 is a young analog of HR 8799. We further investigate and constrain single-planet, two-planet, three-planet and four-planet architectures that can account for the observed debris structure and are compatible with dynamical stability constraints. We find that equal-mass four-planet configurations of geometrically spaced orbits, with each planet of mass 5 MJ, could maintain the gap between the warm and cold debris belts, and also be just marginally stable for timescales comparable to the age of the system.

We investigate the orbital evolution of particles in a planet's chaotic zone to determine their final destinations and their timescales of clearing. There are four possible final states of chaotic particles: collision with the planet, collision with the star, escape, or bounded but non-collision orbits. In our investigations, within the framework of the planar circular restricted three body problem for planet-star mass ratio $\mu$ in the range $10^{-9}$ to $10^{-1.5}$, we find no particles hitting the star. The relative frequencies of escape and collision with the planet are not scale-free, as they depend upon the size of the planet. For planet radius $R_p\ge0.001R_H$ where $R_H$ is the planet's Hill radius, we find that most chaotic zone particles collide with the planet for $\mu\lesssim10^{-5}$; particle scattering to large distances is significant only for higher mass planets. For fixed ratio $R_p/R_H$, the particle clearing timescale, $T_{cl}$, has a broken power-law dependence on $\mu$. A shallower power-law, $T_{cl}\sim \mu^{-{1/3}}$, prevails at small $\mu$ where particles are cleared primarily by collisions with the planet; a steeper power law, $T_{cl}\sim\mu^{-{3/2}}$, prevails at larger $\mu$ where scattering dominates the particle loss. In the limit of vanishing planet radius, we find $T_{cl}\approx0.024\mu^{-{3\over2}}$. The interior and exterior boundaries of the annular zone in which chaotic particles are cleared are increasingly asymmetric about the planet's orbit for larger planet masses; the inner boundary coincides well with the classical first order resonance overlap zone, $\Delta a_{cl,int}\simeq1.2\mu^{0.28}a_p$; the outer boundary is better described by $\Delta a_{cl,ext}\simeq1.7\mu^{0.31}a_p$, where $a_p$ is the planet-star separation.

We report on the MIT Epoch of Reionization (MITEoR) experiment, a pathfinder low-frequency radio interferometer whose goal is to test technologies that improve the calibration precision and reduce the cost of the high-sensitivity 3D mapping required for 21 cm cosmology. MITEoR accomplishes this by using massive baseline redundancy, which enables both automated precision calibration and correlator cost reduction. We demonstrate and quantify the power and robustness of redundancy for scalability and precision. We find that the calibration parameters precisely describe the effect of the instrument upon our measurements, allowing us to form a model that is consistent with $\chi^2$ per degree of freedom < 1.2 for as much as 80% of the observations. We use these results to develop an optimal estimator of calibration parameters using Wiener filtering, and explore the question of how often and how finely in frequency visibilities must be reliably measured to solve for calibration coefficients. The success of MITEoR with its 64 dual-polarization elements bodes well for the more ambitious Hydrogen Epoch of Reionization Array (HERA) project and other next-generation instruments, which would incorporate many identical or similar technologies.

Mapping our universe in 3D by imaging the redshifted 21 cm line from neutral hydrogen has the potential to overtake the cosmic microwave background as our most powerful cosmological probe, because it can map a much larger volume of our Universe, shedding new light on the epoch of reionization, inflation, dark matter, dark energy, and neutrino masses. We report on MITEoR, a pathfinder low-frequency radio interferometer whose goal is to test technologies that greatly reduce the cost of such 3D mapping for a given sensitivity. MITEoR accomplishes this by using massive baseline redundancy both to enable automated precision calibration and to cut the correlator cost scaling from N^2 to NlogN, where N is the number of antennas. The success of MITEoR with its 64 dual-polarization elements bodes well for the more ambitious HERA project, which would incorporate many identical or similar technologies using an order of magnitude more antennas, each with dramatically larger collecting area.

The existence of strong absorption bands of singly deuterated methane (CH3D) at wavelengths where normal methane (CH4) absorbs comparatively weakly could enable remote measurement of D/H ratios in methane ice on outer solar system bodies. We performed laboratory transmission spectroscopy experiments, recording spectra at wavelengths from 1 to 6 \mum to study CH3D bands at 2.47, 2.87, and 4.56 \mum, wavelengths where ordinary methane absorption is weak. We report temperature-dependent absorption coefficients of these bands when the CH3D is diluted in CH4 ice and also when it is dissolved in N2 ice, and describe how these absorption coefficients can be combined with data from the literature to simulate arbitrary D/H ratio absorption coefficients for CH4 ice and for CH4 in N2 ice. We anticipate these results motivating new telescopic observations to measure D/H ratios in CH4 ice on Triton, Pluto, Eris, and Makemake.

Growing interest in 21 cm tomography has led to the design and construction of broadband radio interferometers with low noise, moderate angular resolution, high spectral resolution, and wide fields of view. With characteristics somewhat different from traditional radio instruments, these interferometers may require new calibration techniques in order to reach their design sensitivities. Self-calibration or redundant calibration techniques that allow an instrument to be calibrated off complicated sky emission structures are ideal. In particular, the large number of redundant baselines possessed by these new instruments makes redundant calibration an especially attractive option. In this paper, we explore the errors and biases in existing redundant calibration schemes through simulations, and show how statistical biases can be eliminated. We also develop a general calibration formalism that includes both redundant baseline methods and basic point source calibration methods as special cases, and show how slight deviations from perfect redundancy and coplanarity can be taken into account.