Showing 1–50 of 129 results for author: Thompson, R

Two methods to analyse radial diffusion ensembles: the peril of space- and time- dependent diffusion

Authors: Sarah N. Bentley, J. Stout, Rhys Thompson, Daniel J. Ratliff, Clare E. J. Watt

Abstract: Particle dynamics in Earth's outer radiation belt can be modelled using a diffusion framework, where large-scale electron movements are captured by a diffusion equation across a single adiabatic invariant, $L^{*}$ $``(L)"$. While ensemble models are promoted to represent physical uncertainty, as yet there is no validated method to analyse radiation belt ensembles. Comparisons are complicated by th… ▽ More Particle dynamics in Earth's outer radiation belt can be modelled using a diffusion framework, where large-scale electron movements are captured by a diffusion equation across a single adiabatic invariant, $L^{*}$ $``(L)"$. While ensemble models are promoted to represent physical uncertainty, as yet there is no validated method to analyse radiation belt ensembles. Comparisons are complicated by the domain dependent diffusion, since diffusion coefficient $D_{LL}$ is dependent on $L$. We derive two tools to analyse ensemble members: time to monotonicity $t_m$ and mass/energy moment quantities $\mathcal{N}, \mathcal{E}$. We find that the Jacobian ($1/L^2$) is necessary for radiation belt error metrics. Components of $\partial\mathcal{E}/\partial t$ are explicitly calculated to compare the effects of outer and inner boundary conditions, and loss, on the ongoing diffusion. Using $t_m$, $\mathcal{N}$ and $\mathcal{E}$, we find that: (a) different physically motivated choices of outer boundary condition and location result in different final states and different rates of evolution; (b) the gradients of the particle distribution affect evolution more significantly than $D_{LL}$; (c) the enhancement location, and the amount of initial background particles, are both significant factors determining system evolution; (d) loss from pitch-angle scattering is generally dominant; it mitigates but does not remove the influence of both initial conditions and outer boundary settings, which are due to the $L$-dependence of $D_{LL}$. We anticipate this study will promote renewed focus on the distribution gradients, on the location and nature of the outer boundary in radiation belt modelling, and provide a foundation for systematic ensemble modelling. △ Less

Submitted 5 July, 2024; originally announced July 2024.

Adaptive sampling with PIXL on the Mars Perseverance rover

Authors: Peter R. Lawson, Tanya V. Kizovski, Michael M. Tice, Benton C. Clark III, Scott J. VanBommel, David R. Thompson, Lawrence A. Wade, Robert W. Denise, Christopher M. Heirwegh, W. Timothy Elam, Mariek E. Schmidt, Yang Liu, Abigail C. Allwood, Martin S. Gilbert, Benjamin J. Bornstein

Abstract: Planetary rovers can use onboard data analysis to adapt their measurement plan on the fly, improving the science value of data collected between commands from Earth. This paper describes the implementation of an adaptive sampling algorithm used by PIXL, the X-ray fluorescence spectrometer of the Mars 2020 Perseverance rover. PIXL is deployed using the rover arm to measure X-ray spectra of rocks wi… ▽ More Planetary rovers can use onboard data analysis to adapt their measurement plan on the fly, improving the science value of data collected between commands from Earth. This paper describes the implementation of an adaptive sampling algorithm used by PIXL, the X-ray fluorescence spectrometer of the Mars 2020 Perseverance rover. PIXL is deployed using the rover arm to measure X-ray spectra of rocks with a scan density of several thousand points over an area of typically 5 x 7 mm. The adaptive sampling algorithm is programmed to recognize points of interest and to increase the signal-to-noise ratio at those locations by performing longer integrations. Two approaches are used to formulate the sampling rules based on past quantification data: 1) Expressions that isolate particular regions within a ternary compositional diagram, and 2) Machine learning rules that threshold for a high weight percent of particular compounds. The design of the rulesets are outlined and the performance of the algorithm is quantified using measurements from the surface of Mars. To our knowledge, PIXL's adaptive sampling represents the first autonomous decision-making based on real-time compositional analysis by a spacecraft on the surface of another planet. △ Less

Submitted 23 May, 2024; originally announced May 2024.

Comments: 24 pages including 11 figures and 7 tables. Submitted for publication to the journal Icarus

Evolution of Cosmological Parameters and Fundamental Constants in a Flat Quintessence Cosmology: A Dynamical Alternative to ΛCDM

Authors: Rodger I. Thompson

Abstract: The primary purpose of this work is the provision of accurate, analytic, evolutionary templates for cosmological parameters and fundamental constants in a dynamical cosmology. A flat quintessence cosmology with a dark energy potential that has the mathematical form of the Higgs potential is the specific cosmology and potential addressed in this work. These templates, based on the physics of the co… ▽ More The primary purpose of this work is the provision of accurate, analytic, evolutionary templates for cosmological parameters and fundamental constants in a dynamical cosmology. A flat quintessence cosmology with a dark energy potential that has the mathematical form of the Higgs potential is the specific cosmology and potential addressed in this work. These templates, based on the physics of the cosmology and potential are intended to replace the parameterizations currently used to determine the likelihoods of dynamical cosmologies. Acknowledging that, unlike ΛCDM, the evolutions are dependent on both the specific cosmology and the dark energy potential the templates are referred to as Specific Cosmology and Potential, SCP, templates. The requirements set for the SCP templates are that they must be accurate, analytic functions of an observable such as the scale factor or redshift. This is achieved through the utilization of a modified beta function formalism that is based on a physically motivated dark energy potential to calculate the beta function. The methodology developed here is designed to be adaptable to other cosmologies and dark energy potentials. The SCP templates are essential tools in determining the relative likelihoods of a range of dynamical cosmologies and potentials. An ultimate purpose is the determination whether dark energy is dynamical or static in a quantitative manner. It is suggested that the SCP templates calculated in this work can serve as fiducial dynamical templates in the same manner as ΛCDM serves for static dark energy. △ Less

Submitted 15 May, 2023; originally announced May 2023.

Comments: Published in the Journal Universe

Journal ref: Universe, 2023, 9, 172-199

3D hydrodynamics simulations of internal gravity waves in red giant branch stars

Authors: Simon Blouin, Huaqing Mao, Falk Herwig, Pavel Denissenkov, Paul R. Woodward, William R. Thompson

Abstract: We present the first 3D hydrodynamics simulations of the excitation and propagation of internal gravity waves (IGWs) in the radiative interiors of low-mass stars on the red giant branch (RGB). We use the PPMstar explicit gas dynamics code to simulate a portion of the convective envelope and all the radiative zone down to the hydrogen-burning shell of a 1.2$M_{\odot}$ upper RGB star. We perform sim… ▽ More We present the first 3D hydrodynamics simulations of the excitation and propagation of internal gravity waves (IGWs) in the radiative interiors of low-mass stars on the red giant branch (RGB). We use the PPMstar explicit gas dynamics code to simulate a portion of the convective envelope and all the radiative zone down to the hydrogen-burning shell of a 1.2$M_{\odot}$ upper RGB star. We perform simulations for different grid resolutions (768$^3$, 1536$^3$ and 2880$^3$), a range of driving luminosities, and two different stratifications (corresponding to the bump luminosity and the tip of the RGB). Our RGB tip simulations can be directly performed at the nominal luminosity, circumventing the need for extrapolations to lower luminosities. A rich, continuous spectrum of IGWs is observed, with a significant amount of total power contained at high wavenumbers. By following the time evolution of a passive dye in the stable layers, we find that IGW mixing in our simulations is weaker than predicted by a simple analytical prescription based on shear mixing and not efficient enough to explain the missing RGB extra mixing. However, we may be underestimating the efficiency of IGW mixing given that our simulations include a limited portion of the convective envelope. Quadrupling its radial extent compared to our fiducial setup increases convective velocities by up to a factor 2 and IGW velocities by up to a factor 4. We also report the formation of a $\sim 0.2\,H_P$ penetration zone and evidence that IGWs are excited by plumes that overshoot into the stable layers. △ Less

Submitted 13 March, 2023; originally announced March 2023.

Comments: 20 pages, 29 figures, movies at https://www.ppmstar.org, resubmitted to MNRAS

3D hydrodynamic simulations of massive main-sequence stars. I. Dynamics and mixing of convection and internal gravity waves

Authors: Falk Herwig, Paul R. Woodward, Huaqing Mao, William R. Thompson, Pavel Denissenkov, Josh Lau, Simon Blouin, Robert Andrassy, Adam Paul

Abstract: We performed 3D hydrodynamic simulations of the inner $\approx 50\%$ radial extent of a $25\ \mathrm{M_\odot}$ star in the early phase of the main sequence and investigate core convection and internal gravity waves in the core-envelope boundary region. Simulations for different grid resolutions and driving luminosities establish scaling relations to constrain models of mixing for 1D applications.… ▽ More We performed 3D hydrodynamic simulations of the inner $\approx 50\%$ radial extent of a $25\ \mathrm{M_\odot}$ star in the early phase of the main sequence and investigate core convection and internal gravity waves in the core-envelope boundary region. Simulations for different grid resolutions and driving luminosities establish scaling relations to constrain models of mixing for 1D applications. As in previous works, the turbulent mass entrainment rate extrapolated to nominal heating is unrealistically high ($1.58\times 10^{-4}\ \mathrm{M_\odot/yr}$), which is discussed in terms of the non-equilibrium response of the simulations to the initial stratification. We measure quantitatively the effect of mixing due to internal gravity waves excited by core convection interacting with the boundary in our simulations. The wave power spectral density as a function of frequency and wavelength agrees well with the GYRE eigenmode predictions based on the 1D spherically averaged radial profile. A diffusion coefficient profile that reproduces the spherically averaged abundance distribution evolution is determined for each simulation. Through a combination of eigenmode analysis and scaling relations it is shown that in the $N^2$-peak region, mixing is due to internal gravity waves and follows the scaling relation $D_\mathrm{IGW-hydro} \propto L^{4/3}$ over a $\gtrapprox 2\ \mathrm{dex}$ range of heating factors. Different extrapolations of the mixing efficiency down to nominal heating are discussed. If internal gravity wave mixing is due to thermally-enhanced shear mixing, an upper limit is $D_\mathrm{IGW} \lessapprox 2$ to $3\times 10^4\ \mathrm{cm^2/s}$ at nominal heating in the $N^2$-peak region above the convective core. △ Less

Submitted 19 July, 2023; v1 submitted 9 March, 2023; originally announced March 2023.

Comments: MNRAS. Clarifications and improvements following referee recommendations. Movies and data access https://www.ppmstar.org. 38 Figs. In v2 the bottom panel of Fig 28 is showing the wrong case (M116 instead of M114). Corrected in v3

Dynamical templates for comparison to Lambda CDN: Static or dynamical dark energy?

Authors: Rodger I. Thompson

Abstract: This study forges new tools to discriminate between dynamical and static dark energy. It provides accurate evolutionary templates of dynamical cosmological parameters and fundamental constants as analytic functions of the scale factor. They are designed to replace the commonly used parameterizations in likelihood calculations with evolutionary templates based on the physics of specific dynamical c… ▽ More This study forges new tools to discriminate between dynamical and static dark energy. It provides accurate evolutionary templates of dynamical cosmological parameters and fundamental constants as analytic functions of the scale factor. They are designed to replace the commonly used parameterizations in likelihood calculations with evolutionary templates based on the physics of specific dynamical cosmologies and dark energy potentials. Thus they are termed Cosmology and Potential Specific, CPS, templates. A suite of CPS templates are calculated for a flat quintessence cosmology with a dark energy potential of the same mathematical form as the Higgs potential. This Higgs inspired, HI, polynomial potential produces a rich set of evolutions unlike most monomial potentials. The study produces CPS templates that are analytic functions of the scale factor. It uses a recently developed beta function formalism that provides a differential function for the scalar in terms of the scale factor. This establishes a methodology for easily producing CPS templates for other dark energy potentials and cosmologies to determine the likelihoods of dynamical cosmologies relative to Lambda CDM. The study also examines the evolution of fundamental constants such as the proton to electron mass ratio and the fine structure constant involving an intersection between particle physics and cosmology. Appendix A displays an abridged suite of CPS templates for flat quintessence and the HI dark energy potential. △ Less

Submitted 4 April, 2022; originally announced April 2022.

Comments: For submission to the Journal of Cosmology and Astroparticle Physics

Measuring the stability of fundamental constants with a network of clocks

Authors: G. Barontini, L. Blackburn, V. Boyer, F. Butuc-Mayer, X. Calmet, J. R. Crespo Lopez-Urrutia, E. A. Curtis, B. Darquie, J. Dunningham, N. J. Fitch, E. M. Forgan, K. Georgiou, P. Gill, R. M. Godun, J. Goldwin, V. Guarrera, A. C. Harwood, I. R. Hill, R. J. Hendricks, M. Jeong, M. Y. H. Johnson, M. Keller, L. P. Kozhiparambil Sajith, F. Kuipers, H. S. Margolis , et al. (19 additional authors not shown)

Abstract: The detection of variations of fundamental constants of the Standard Model would provide us with compelling evidence of new physics, and could lift the veil on the nature of dark matter and dark energy. In this work, we discuss how a network of atomic and molecular clocks can be used to look for such variations with unprecedented sensitivity over a wide range of time scales. This is precisely the… ▽ More The detection of variations of fundamental constants of the Standard Model would provide us with compelling evidence of new physics, and could lift the veil on the nature of dark matter and dark energy. In this work, we discuss how a network of atomic and molecular clocks can be used to look for such variations with unprecedented sensitivity over a wide range of time scales. This is precisely the goal of the recently launched QSNET project: A network of clocks for measuring the stability of fundamental constants. QSNET will include state-of-the-art atomic clocks, but will also develop next-generation molecular and highly charged ion clocks with enhanced sensitivity to variations of fundamental constants. We describe the technological and scientific aims of QSNET and evaluate its expected performance. We show that in the range of parameters probed by QSNET, either we will discover new physics, or we will impose new constraints on violations of fundamental symmetries and a range of theories beyond the Standard Model, including dark matter and dark energy models. △ Less

Submitted 11 May, 2022; v1 submitted 20 December, 2021; originally announced December 2021.

Journal ref: EPJ Quantum Technology volume 9, Article number: 12 (2022)

Impact of Scene-Specific Enhancement Spectra on Matched Filter Greenhouse Gas Retrievals from Imaging Spectroscopy

Authors: Markus D. Foote, Philip E. Dennison, Patrick R. Sullivan, Kelly B. O'Neill, Andrew K. Thorpe, David R. Thompson, Daniel H. Cusworth, Riley Duren, Sarang C. Joshi

Abstract: Matched filter (MF) techniques have been widely used for retrieval of greenhouse gas enhancements (enh.) from imaging spectroscopy datasets. While multiple algorithmic techniques and refinements have been proposed, the greenhouse gas target spectrum used for concentration enh. estimation has remained largely unaltered since the introduction of quantitative MF retrievals. The magnitude of retrieved… ▽ More Matched filter (MF) techniques have been widely used for retrieval of greenhouse gas enhancements (enh.) from imaging spectroscopy datasets. While multiple algorithmic techniques and refinements have been proposed, the greenhouse gas target spectrum used for concentration enh. estimation has remained largely unaltered since the introduction of quantitative MF retrievals. The magnitude of retrieved methane and carbon dioxide enh., and thereby integrated mass enh. (IME) and estimated flux of point-source emitters, is heavily dependent on this target spectrum. Current standard use of molecular absorption coefficients to create unit enh. target spectra does not account for absorption by background concentrations of greenhouse gases, solar and sensor geometry, or atmospheric water vapor absorption. We introduce geometric and atmospheric parameters into the generation of scene-specific (SS) unit enh. spectra to provide target spectra that are compatible with all greenhouse gas retrieval MF techniques. For methane plumes, IME resulting from use of standard, generic enh. spectra varied from -22 to +28.7% compared to SS enh. spectra. Due to differences in spectral shape between the generic and SS enh. spectra, differences in methane plume IME were linked to surface spectral characteristics in addition to geometric and atmospheric parameters. IME differences for carbon dioxide plumes, with generic enh. spectra producing integrated mass enh. -76.1 to -48.1% compared to SS enh. spectra. Fluxes calculated from these integrated enh. would vary by the same %s, assuming equivalent wind conditions. Methane and carbon dioxide IME were most sensitive to changes in solar zenith angle and ground elevation. SS target spectra can improve confidence in greenhouse gas retrievals and flux estimates across collections of scenes with diverse geometric and atmospheric conditions. △ Less

Submitted 10 August, 2021; v1 submitted 25 June, 2021; originally announced July 2021.

Comments: 13 pages, 5 figures, 3 tables

Journal ref: Remote Sensing of Environment, Volume 264, October 2021, 112574

Powderday: Dust Radiative Transfer for Galaxy Simulations

Authors: Desika Narayanan, Matthew J. Turk, Thomas Robitaille, Ashley J. Kelly, B. Connor McClellan, Ray S. Sharma, Prerak Garg, Matthew Abruzzo, Ena Choi, Charlie Conroy, Benjamin D. Johnson, Benjamin Kimock, Qi Li, Christopher C. Lovell, Sidney Lower, George C. Privon, Jonathan Roberts, Snigdaa Sethuram, Gregory F. Snyder, Robert Thompson, John H. Wise

Abstract: We present Powderday, a flexible, fast, open-source dust radiative transfer package designed to interface with galaxy formation simulations. Powderday builds on FSPS population synthesis models, Hyperion dust radiative transfer, and employs yt to interface between different software packages. We include our stellar population synthesis modeling on the fly, which allows for significant run-time fle… ▽ More We present Powderday, a flexible, fast, open-source dust radiative transfer package designed to interface with galaxy formation simulations. Powderday builds on FSPS population synthesis models, Hyperion dust radiative transfer, and employs yt to interface between different software packages. We include our stellar population synthesis modeling on the fly, which allows for significant run-time flexibility in the assumed stellar physics. We include a model for nebular line emission that can employ either precomputed Cloudy lookup tables (for efficiency), or direct photoionization calculations for all young stars (for flexibility). The dust content follows either observationally-motivated prescriptions, direct modeling from galaxy formation simulations, or a novel approach that includes the dust content via learning-based algorithms from the SIMBA cosmological galaxy formation simulation. AGN can additionally be included via a range of prescriptions. The output of these models are broadband SEDs, as well as filter-convolved images. Powderday is designed to eliminate last-mile efforts by researchers that employ different hydrodynamic galaxy formation models, and seamlessly interfaces with GIZMO, AREPO, GASOLINE, CHANGA, and ENZO. We demonstrate the capabilities of the code via three applications: a model for the star formation rate (SFR) - infrared luminosity relation in galaxies (including the impact of AGN); the impact of circumstellar dust around AGB stars on the mid-infrared emission from galaxy SEDs; and the impact of galaxy inclination angle on dust attenuation laws. △ Less

Submitted 18 June, 2020; originally announced June 2020.

Comments: Submitted to ApJ; comments welcome

Testing bound dark energy with cosmological parameter and fundamental constant evolution

Authors: Rodger I. Thompson

Abstract: A new bound dark energy, BDE, cosmology has been proposed where the dark energy is the binding energy between light meson fields that condense a few tens of years after the big bang. It is reported that the correct dark energy density emerges using particle physics without fine tuning. This alone makes the BDE cosmology worthy of further investigation. This work looks at the late time BDE predicti… ▽ More A new bound dark energy, BDE, cosmology has been proposed where the dark energy is the binding energy between light meson fields that condense a few tens of years after the big bang. It is reported that the correct dark energy density emerges using particle physics without fine tuning. This alone makes the BDE cosmology worthy of further investigation. This work looks at the late time BDE predictions of the evolution of cosmological parameters and the values of fundamental constants to determine whether the cosmology's predictions are consistent with observation. The work considers the time period between a scale factor of 0.1 and 1.0. A model BDE cosmology is considered with current day values of the cosmological parameters well within the observational limits. The calculations use three different values of the current day dark energy equation of state close to minus one. All three cases produce evolutions of the cosmological parameters and fundamental constants consistent with the observational constraints. Analytic relations between the BDE and cosmological parameters are developed to insure a consistent set of parameters. △ Less

Submitted 10 October, 2019; originally announced October 2019.

Comments: Accepted for publication in MNRAS

Learning Radiative Transfer Models for Climate Change Applications in Imaging Spectroscopy

Authors: Shubhankar Deshpande, Brian D. Bue, David R. Thompson, Vijay Natraj, Mario Parente

Abstract: According to a recent investigation, an estimated 33-50% of the world's coral reefs have undergone degradation, believed to be as a result of climate change. A strong driver of climate change and the subsequent environmental impact are greenhouse gases such as methane. However, the exact relation climate change has to the environmental condition cannot be easily established. Remote sensing methods… ▽ More According to a recent investigation, an estimated 33-50% of the world's coral reefs have undergone degradation, believed to be as a result of climate change. A strong driver of climate change and the subsequent environmental impact are greenhouse gases such as methane. However, the exact relation climate change has to the environmental condition cannot be easily established. Remote sensing methods are increasingly being used to quantify and draw connections between rapidly changing climatic conditions and environmental impact. A crucial part of this analysis is processing spectroscopy data using radiative transfer models (RTMs) which is a computationally expensive process and limits their use with high volume imaging spectrometers. This work presents an algorithm that can efficiently emulate RTMs using neural networks leading to a multifold speedup in processing time, and yielding multiple downstream benefits. △ Less

Submitted 8 June, 2019; originally announced June 2019.

Comments: Accepted to International Conference on Machine Learning (ICML) 2019 Workshop: Climate Change: How Can AI Help?

Beta Function Quintessence Cosmological Parameters and Fundamental Constants II: Exponential and Logarithmic Dark Energy Potentials

Authors: Rodger I. Thompson

Abstract: This paper uses the beta function formalism to extend the analysis of quintessence cosmological parameters to the logarithmic and exponential dark energy potentials. The previous paper (Thompson 2018) demonstrated the formalism using power and inverse power potentials. The essentially identical evolution of the Hubble parameter for all of the quintessence cases and LambdaCDM is attributed to the f… ▽ More This paper uses the beta function formalism to extend the analysis of quintessence cosmological parameters to the logarithmic and exponential dark energy potentials. The previous paper (Thompson 2018) demonstrated the formalism using power and inverse power potentials. The essentially identical evolution of the Hubble parameter for all of the quintessence cases and LambdaCDM is attributed to the flatness of the quintessence dark energy potentials in the dark energy dominated era. The Hubble parameter is therefore incapable of discriminating between static and dynamic dark energy. Unlike the other three potentials considered in the two papers the logarithmic dark energy potential requires a numerical integration in the formula for the superpotential rather than being an analytic function. The dark energy equation of state and the fundamental constants continue to be good discriminators between static and dynamical dark energy. A new analysis of quintessence with all four of the potentials relative the swampland conjectures indicates that the conjecture on the change in the scalar field is satisfied but that the conjecture on the change of the potential is not. △ Less

Submitted 7 November, 2018; originally announced November 2018.

Comments: Accepted for publication in the MNRAS,

The Robustness of Cosmological Hydrodynamic Simulation Predictions to Changes in Numerics and Cooling Physics

Authors: Shuiyao Huang, Neal Katz, Romeel Davé, Mark Fardal, Juna Kollmeier, Benjamin D. Oppenheimer, Molly S. Peeples, Shawn Roberts, David H. Weinberg, Philip F. Hopkins, Robert Thompson

Abstract: We test and improve the numerical schemes in our smoothed particle hydrodynamics (SPH) code for cosmological simulations, including the pressure-entropy formulation (PESPH), a time-dependent artificial viscosity, a refined timestep criterion, and metal-line cooling that accounts for photoionisation in the presence of a recently refined Haardt \& Madau (2012) model of the ionising background. The P… ▽ More We test and improve the numerical schemes in our smoothed particle hydrodynamics (SPH) code for cosmological simulations, including the pressure-entropy formulation (PESPH), a time-dependent artificial viscosity, a refined timestep criterion, and metal-line cooling that accounts for photoionisation in the presence of a recently refined Haardt \& Madau (2012) model of the ionising background. The PESPH algorithm effectively removes the artificial surface tension present in the traditional SPH formulation, and in our test simulations it produces better qualitative agreement with mesh-code results for Kelvin-Helmholtz instability and cold cloud disruption. Using a set of cosmological simulations, we examine many of the quantities we have studied in previous work. Results for galaxy stellar and HI mass functions, star formation histories, galaxy scaling relations, and statistics of the Ly$α$ forest are robust to the changes in numerics and microphysics. As in our previous simulations, cold gas accretion dominates the growth of high-redshift galaxies and of low mass galaxies at low redshift, and recycling of winds dominates the growth of massive galaxies at low redshift. However, the PESPH simulation removes spurious cold clumps seen in our earlier simulations, and the accretion rate of hot gas increases by up to an order of magnitude at some redshifts. The new numerical model also influences the distribution of metals among gas phases, leading to considerable differences in the statistics of some metal absorption lines, most notably NeVIII. △ Less

Submitted 9 January, 2020; v1 submitted 30 October, 2018; originally announced October 2018.

Comments: 29 pages, 25 figures, accepted by MNRAS

Mass Measurements of Neutron-Rich Gallium Isotopes Refine Production of Nuclei of the First r-Process Abundance Peak in Neutron Star Merger Calculations

Authors: M. P. Reiter, S. Ayet San Andrés, S. Nikas, J. Lippuner, C. Andreoiu, C. Babcock, B. R. Barquest, J. Bollig, T. Brunner, T. Dickel, J. Dilling, I. Dillmann, E. Dunling, G. Gwinner, L. Graham, C. Hornung, R. Klawitter, B. Kootte, A. A. Kwiatkowski, Y. Lan, D. Lascar, K. G. Leach, E. Leistenschneider, G. Martínez-Pinedo, J. E. McKay , et al. (11 additional authors not shown)

Abstract: We report mass measurements of neutron-rich Ga isotopes $^{80-85}$Ga with TRIUMF's Ion Trap for Atomic and Nuclear science (TITAN). The measurements determine the masses of $^{80-83}$Ga in good agreement with previous measurements. The masses of $^{84}$Ga and $^{85}$Ga were measured for the first time. Uncertainties between $25-48$ keV were reached. The new mass values reduce the nuclear uncertain… ▽ More We report mass measurements of neutron-rich Ga isotopes $^{80-85}$Ga with TRIUMF's Ion Trap for Atomic and Nuclear science (TITAN). The measurements determine the masses of $^{80-83}$Ga in good agreement with previous measurements. The masses of $^{84}$Ga and $^{85}$Ga were measured for the first time. Uncertainties between $25-48$ keV were reached. The new mass values reduce the nuclear uncertainties associated with the production of A $\approx$ 84 isotopes by the \emph{r}-process for astrophysical conditions that might be consistent with a binary neutron star (BNS) merger producing a blue kilonova. Our nucleosynthesis simulations confirm that BNS merger may contribute to the first abundance peak under moderate neutron-rich conditions with electron fractions $Y_e=0.35-0.38$. △ Less

Submitted 30 January, 2020; v1 submitted 26 October, 2018; originally announced October 2018.

Journal ref: Phys. Rev. C 101, 025803 (2020)

Beta Function Quintessence Cosmological Parameters and Fundamental Constants I: Power and Inverse Power Law Dark Energy Potentials

Authors: Rodger I. Thompson

Abstract: This investigation explores using the beta function formalism to calculate analytic solutions for the observable parameters in rolling scalar field cosmologies. The beta function in this case is the derivative of the scalar $φ$ with respect to the natural log of the scale factor $a$, $β(φ)=\frac{d φ}{d \ln(a)}$. Once the beta function is specified, modulo a boundary condition, the evolution of t… ▽ More This investigation explores using the beta function formalism to calculate analytic solutions for the observable parameters in rolling scalar field cosmologies. The beta function in this case is the derivative of the scalar $φ$ with respect to the natural log of the scale factor $a$, $β(φ)=\frac{d φ}{d \ln(a)}$. Once the beta function is specified, modulo a boundary condition, the evolution of the scalar $φ$ as a function of the scale factor is completely determined. A rolling scalar field cosmology is defined by its action which can contain a range of physically motivated dark energy potentials. The beta function is chosen so that the associated "beta potential" is an accurate, but not exact, representation of the appropriate dark energy model potential. The basic concept is that the action with the beta potential is so similar to the action with the model potential that solutions using the beta action are accurate representations of solutions using the model action. The beta function provides an extra equation to calculate analytic functions of the cosmologies parameters as a function of the scale factor that are that are not calculable using only the model action. As an example this investigation uses a quintessence cosmology to demonstrate the method for power and inverse power law dark energy potentials. An interesting result of the investigation is that the Hubble parameter H is almost completely insensitive to the power of the potentials and that $Λ$CDM is part of the family of quintessence cosmology power law potentials with a power of zero. △ Less

Submitted 6 April, 2018; originally announced April 2018.

Comments: Accepted for publication by the Monthly Notices of the Royal Astronomical Society

Identifying Mergers Using Quantitative Morphologies in Zoom Simulations of High-Redshift Galaxies

Authors: Matthew W. Abruzzo, Desika Narayanan, Romeel Davé, Robert Thompson

Abstract: Non-parametric morphology measures are a powerful tool for identifying galaxy mergers at low redshifts. We employ cosmological zoom simulations using Gizmo with the Mufasa feedback scheme, post-processed using 3D dust radiative transfer into mock observations, to study whether common morphological measures Gini G, M20, concentration C, and asymmetry A are effective at identifying major galaxy merg… ▽ More Non-parametric morphology measures are a powerful tool for identifying galaxy mergers at low redshifts. We employ cosmological zoom simulations using Gizmo with the Mufasa feedback scheme, post-processed using 3D dust radiative transfer into mock observations, to study whether common morphological measures Gini G, M20, concentration C, and asymmetry A are effective at identifying major galaxy mergers at z ~ 2 - 4, i.e. "Cosmic Noon". Our zoom suite covers galaxies with 10^8.6 < M_* < 10^11 M_sun at z ~ 2, and broadly reproduces key global galaxy observations. Our primary result is that these morphological measures are unable to robustly pick out galaxies currently undergoing mergers during Cosmic Noon, typically performing no better than a random guess. This improves only marginally if we consider whether galaxies have undergone a merger within the last Gyr. When also considering minor mergers, galaxies display no trend of moving towards the merger regime with increasing merger ratio. From z = 4 -> 2, galaxies move from the non-merger towards the merger regime in all statistics, but this is primarily an effect of mass: Above a given noise level, higher mass galaxies display a more complex outer morphology induced by their clustered environment. We conclude that during Cosmic Noon, these morphological statistics are of limited value in identifying galaxy mergers. △ Less

Submitted 6 March, 2018; originally announced March 2018.

Comments: submitted to MNRAS, comments welcome

Investigating the physical properties of transiting hot Jupiters with the 1.5-m Kuiper Telescope

Authors: Jake D. Turner, Robin M. Leiter, Lauren I. Biddle, Kyle A. Pearson, Kevin K. Hardegree-Ullman, Robert M. Thompson, Johanna K. Teske, Ian T. Cates, Kendall L. Cook, Michael P. Berube, Megan N. Nieberding, Christen K. Jones, Brandon Raphael, Spencer Wallace, Zachary T. Watson, Robert E. Johnson

Abstract: We present new photometric data of 11 hot Jupiter transiting exoplanets (CoRoT-12b, HAT-P-5b, HAT-P-12b, HAT-P-33b, HAT-P-37b, WASP-2b, WASP-24b, WASP-60b, WASP-80b, WASP-103b, XO-3b) in order to update their planetary parameters and to constrain information about their atmospheres. These observations of CoRoT-12b, HAT-P-37b and WASP-60b are the first follow-up data since their discovery. Addition… ▽ More We present new photometric data of 11 hot Jupiter transiting exoplanets (CoRoT-12b, HAT-P-5b, HAT-P-12b, HAT-P-33b, HAT-P-37b, WASP-2b, WASP-24b, WASP-60b, WASP-80b, WASP-103b, XO-3b) in order to update their planetary parameters and to constrain information about their atmospheres. These observations of CoRoT-12b, HAT-P-37b and WASP-60b are the first follow-up data since their discovery. Additionally, the first near-UV transits of WASP-80b and WASP-103b are presented. We compare the results of our analysis with previous work to search for transit timing variations (TTVs) and a wavelength dependence in the transit depth. TTVs may be evidence of a third body in the system and variations in planetary radius with wavelength can help constrain the properties of the exoplanet's atmosphere. For WASP-103b and XO-3b, we find a possible variation in the transit depths that may be evidence of scattering in their atmospheres. The B-band transit depth of HAT-P-37b is found to be smaller than its near-IR transit depth and such a variation may indicate TiO/VO absorption. These variations are detected from 2-4.6$σ$, so follow-up observations are needed to confirm these results. Additionally, a flat spectrum across optical wavelengths is found for 5 of the planets (HAT-P-5b, HAT-P-12b, WASP-2b, WASP-24b, WASP-80b), suggestive that clouds may be present in their atmospheres. We calculate a refined orbital period and ephemeris for all the targets, which will help with future observations. No TTVs are seen in our analysis with the exception of WASP-80b and follow-up observations are needed to confirm this possible detection. △ Less

Submitted 25 August, 2017; originally announced August 2017.

Comments: 18 pages, 7 figures, 9 Tables. Light Curves available online. Accepted to MNRAS (2017 August 25)

SIGAME simulations of the [CII], [OI] and [OIII] line emission from star forming galaxies at z ~ 6

Authors: K. P. Olsen, T. R. Greve, D. Narayanan, R. Thompson, R. Davé, L. N. Rios, S. Stawinski

Abstract: Of the almost 40 star forming galaxies at z>~5 (not counting QSOs) observed in [CII] to date, nearly half are either very faint in [CII], or not detected at all, and fall well below expectations based on locally derived relations between star formation rate (SFR) and [CII] luminosity. Combining cosmological zoom simulations of galaxies with SIGAME (SImulator of GAlaxy Millimeter/submillimeter Emis… ▽ More Of the almost 40 star forming galaxies at z>~5 (not counting QSOs) observed in [CII] to date, nearly half are either very faint in [CII], or not detected at all, and fall well below expectations based on locally derived relations between star formation rate (SFR) and [CII] luminosity. Combining cosmological zoom simulations of galaxies with SIGAME (SImulator of GAlaxy Millimeter/submillimeter Emission) we have modeled the multi-phased interstellar medium (ISM) and its emission in [CII], [OI] and [OIII], from 30 main sequence galaxies at z~6 with star formation rates ~3-23Msun/yr, stellar masses ~(0.7-8)x10^9Msun, and metallicities ~(0.1-0.4)xZsun. The simulations are able to reproduce the aforementioned [CII]-faintness at z>5, match two of the three existing z>~5 detections of [OIII], and are furthermore roughly consistent with the [OI] and [OIII] luminosity relations with SFR observed for local starburst galaxies. We find that the [CII] emission is dominated by the diffuse ionized gas phase and molecular clouds, which on average contribute ~66% and ~27%, respectively. The molecular gas, which constitutes only ~10% of the total gas mass is thus a more efficient emitter of [CII] than the ionized gas making up ~85% of the total gas mass. A principal component analysis shows that the [CII] luminosity correlates with the star formation activity as well as average metallicity. The low metallicities of our simulations together with their low molecular gas mass fractions can account for their [CII]-faintness, and we suggest these factors may also be responsible for the [CII]-faint normal galaxies observed at these early epochs. △ Less

Submitted 16 August, 2017; originally announced August 2017.

Comments: 24 pages, 14 figures. Accepted for publication in the Astrophysical Journal

Baseline Metal Enrichment from Population III Star Formation in Cosmological Volume Simulations

Authors: Jason Jaacks, Robert Thompson, Steven L. Finkelstein, Volker Bromm

Abstract: We utilize the hydrodynamic and N-body code {\small GIZMO} coupled with our newly developed sub-grid Population~III (Pop~III) Legacy model, designed specifically for cosmological volume simulations, to study the baseline metal enrichment from Pop~III star formation at $z>7$. In this idealized numerical experiment, we only consider Pop~III star formation. We find that our model Pop~III star formati… ▽ More We utilize the hydrodynamic and N-body code {\small GIZMO} coupled with our newly developed sub-grid Population~III (Pop~III) Legacy model, designed specifically for cosmological volume simulations, to study the baseline metal enrichment from Pop~III star formation at $z>7$. In this idealized numerical experiment, we only consider Pop~III star formation. We find that our model Pop~III star formation rate density (SFRD), which peaks at $\sim 10^{-3}\ {\rm M_\odot yr^{-1} Mpc^{-1}}$ near $z\sim10$, agrees well with previous numerical studies and is consistent with the observed estimates for Pop~II SFRDs. The mean Pop~III metallicity rises smoothly from $z=25-7$, but does not reach the critical metallicity value, $Z_{\rm crit}=10^{-4}\ Z_\odot$, required for the Pop~III to Pop~II transition in star formation mode until $z\simeq7$. This suggests that, while individual halos can suppress in-situ Pop~III star formation, the external enrichment is insufficient to globally terminate Pop~III star formation. The maximum enrichment from Pop~III star formation in star forming dark matter halos is $Z\sim10^{-2}\ Z_\odot$, whereas the minimum found in externally enriched haloes is $Z\gtrsim10^{-7}\ Z_\odot$. Finally, mock observations of our simulated IGM enriched with Pop~III metals produce equivalent widths similar to observations of an extremely metal poor damped Lyman alpha (DLA) system at $z=7.04$, which is thought to be enriched by Pop~III star formation only. △ Less

Submitted 17 January, 2018; v1 submitted 22 May, 2017; originally announced May 2017.

Comments: Published 1/2018 17 pages, 14 figures

The IRX-Beta Dust Attenuation Relation in Cosmological Galaxy Formation Simulations

Authors: Desika Narayanan, Romeel Dave, Benjamin Johnson, Robert Thompson, Charlie Conroy, James E. Geach

Abstract: We utilise a series of high-resolution cosmological zoom simulations of galaxy formation to investigate the relationship between the ultraviolet (UV) slope, beta, and the ratio of the infrared luminosity to UV luminosity (IRX) in the spectral energy distributions (SEDs) of galaxies. We employ dust radiative transfer calculations in which the SEDs of the stars in galaxies propagate through the dust… ▽ More We utilise a series of high-resolution cosmological zoom simulations of galaxy formation to investigate the relationship between the ultraviolet (UV) slope, beta, and the ratio of the infrared luminosity to UV luminosity (IRX) in the spectral energy distributions (SEDs) of galaxies. We employ dust radiative transfer calculations in which the SEDs of the stars in galaxies propagate through the dusty interstellar medium. Our main goals are to understand the origin of, and scatter in the IRX-beta relation; to assess the efficacy of simplified stellar population synthesis screen models in capturing the essential physics in the IRX-beta relation; and to understand systematic deviations from the canonical local IRX-beta relations in particular populations of high-redshift galaxies. Our main results follow. Galaxies that have young stellar populations with relatively cospatial UV and IR emitting regions and a Milky Way-like extinction curve fall on or near the standard Meurer relation. This behaviour is well captured by simplified screen models. Scatter in the IRX-beta relation is dominated by three major effects: (i) older stellar populations drive galaxies below the relations defined for local starbursts due to a reddening of their intrinsic UV SEDs; (ii) complex geometries in high-z heavily star forming galaxies drive galaxies toward blue UV slopes owing to optically thin UV sightlines; (iii) shallow extinction curves drive galaxies downward in the IRX-beta plane due to lowered NUV/FUV extinction ratios. We use these features of the UV slopes of galaxies to derive a fitting relation that reasonably collapses the scatter back toward the canonical local relation. Finally, we use these results to develop an understanding for the location of two particularly enigmatic populations of galaxies in the IRX-beta plane: z~2-4 dusty star forming galaxies, and z>5 star forming galaxies. △ Less

Submitted 16 May, 2017; originally announced May 2017.

Comments: Submitted to MNRAS; comments welcome

MUFASA: The assembly of the red sequence

Authors: Romeel Davé, Mika H. Rafieferantsoa, Robert J. Thompson

Abstract: We examine the growth and evolution of quenched galaxies in the Mufasa cosmological hydrodynamic simulations that include an evolving halo mass-based quenching prescription, with galaxy colours computed accounting for line-of-sight extinction to individual star particles. Mufasa reproduces the observed present-day red sequence reasonably well, including its slope, amplitude, and scatter. In Mufasa… ▽ More We examine the growth and evolution of quenched galaxies in the Mufasa cosmological hydrodynamic simulations that include an evolving halo mass-based quenching prescription, with galaxy colours computed accounting for line-of-sight extinction to individual star particles. Mufasa reproduces the observed present-day red sequence reasonably well, including its slope, amplitude, and scatter. In Mufasa, the red sequence slope is driven entirely by the steep stellar mass-stellar metallicity relation, which independently agrees with observations. High-mass star-forming galaxies blend smoothly onto the red sequence, indicating the lack of a well-defined green valley at M*>10^10.5 Mo. The most massive galaxies quench the earliest and then grow very little in mass via dry merging; they attain their high masses at earlier epochs when cold inflows more effectively penetrate hot halos. To higher redshifts, the red sequence becomes increasingly contaminated with massive dusty star-forming galaxies; UVJ selection subtly but effectively separates these populations. We then examine the evolution of the mass functions of central and satellite galaxies split into passive and star-forming via UVJ. Massive quenched systems show good agreement with observations out to z~2, despite not including a rapid early quenching mode associated with mergers. However, low-mass quenched galaxies are far too numerous at z<1 in Mufasa, indicating that Mufasa strongly over-quenches satellites. A challenge for hydrodynamic simulations is to devise a quenching model that produces enough early massive quenched galaxies and keeps them quenched to z=0, while not being so strong as to over-quench satellites; Mufasa's current scheme fails at the latter. △ Less

Submitted 4 July, 2017; v1 submitted 4 April, 2017; originally announced April 2017.

Comments: 19 pages, MNRAS accepted. Several interesting figures added based on referee's comments

The Relation Between Fundamental Constants and Particle Physics Parameters

Authors: Rodger I. Thompson

Abstract: The observed constraints on the variability of the proton to electron mass ratio $μ$ and the fine structure constant $α$ are used to establish constraints on the variability of the Quantum Chromodynamic Scale and a combination of the Higgs Vacuum Expectation Value and the Yukawa couplings. Further model dependent assumptions provide constraints on the Higgs VEV and the Yukawa couplings separately.… ▽ More The observed constraints on the variability of the proton to electron mass ratio $μ$ and the fine structure constant $α$ are used to establish constraints on the variability of the Quantum Chromodynamic Scale and a combination of the Higgs Vacuum Expectation Value and the Yukawa couplings. Further model dependent assumptions provide constraints on the Higgs VEV and the Yukawa couplings separately. A primary conclusion is that limits on the variability of dimensionless fundamental constants such as $μ$ and $α$ provide important constraints on the parameter space of new physics and cosmologies. △ Less

Submitted 23 February, 2017; originally announced February 2017.

Comments: Published in the proceedings of the Conference on Varying Constants and Fundamental Cosmology VARCOSMOFUN16. Modified from the Universe style to process properly in arXiv

Journal ref: Universe (2017) 3, 1, 6, 1-8

Fundamental Constant Observational Bounds on the Variability of the QCD Scale

Authors: Rodger I. Thompson

Abstract: Many physical theories beyond the Standard Model predict time variations of basic physics parameters. Direct measurement of the time variations of these parameters is very difficult or impossible to achieve. By contrast, measurements of fundamental constants are relatively easy to achieve, both in the laboratory and by astronomical spectra of atoms and molecules in the early universe. In this work… ▽ More Many physical theories beyond the Standard Model predict time variations of basic physics parameters. Direct measurement of the time variations of these parameters is very difficult or impossible to achieve. By contrast, measurements of fundamental constants are relatively easy to achieve, both in the laboratory and by astronomical spectra of atoms and molecules in the early universe. In this work measurements of the proton to electron mass ratio $μ$ and the fine structure constant $α$ are combined to place mildly model dependent limits on the fractional variation of the Quantum Chromodynamic Scale and the sum of the fractional variations of the Higgs Vacuum Expectation Value and the Yukawa couplings on time scales of more than half the age of the universe. The addition of another model parameter allows the fractional variation of the Higgs VEV and the Yukawa couplings to be computed separately. Limits on their variation are found at the level of less than $5 \times 10^{-5}$ over the past seven gigayears. A model dependent relation between the expected fractional variation of $α$ relative to $μ$ tightens the limits to $10^{-7}$ over the same time span. Limits on the present day rate of change of the constants and parameters are then calculated using slow roll quintessence. A primary result of this work is that studies of the dimensionless fundamental constants such as $α$ and $μ$, whose values depend on the values of the physics parameters, are excellent monitors of the limits on the time variation of these parameters. △ Less

Submitted 20 February, 2017; originally announced February 2017.

Comments: Accepted for publication in the Monthly Notices of the Royal Astronomical Society, 8 pages, 5 figures

The AGORA High-Resolution Galaxy Simulations Comparison Project. II: Isolated Disk Test

Authors: Ji-hoon Kim, Oscar Agertz, Romain Teyssier, Michael J. Butler, Daniel Ceverino, Jun-Hwan Choi, Robert Feldmann, Ben W. Keller, Alessandro Lupi, Thomas Quinn, Yves Revaz, Spencer Wallace, Nickolay Y. Gnedin, Samuel N. Leitner, Sijing Shen, Britton D. Smith, Robert Thompson, Matthew J. Turk, Tom Abel, Kenza S. Arraki, Samantha M. Benincasa, Sukanya Chakrabarti, Colin DeGraf, Avishai Dekel, Nathan J. Goldbaum , et al. (18 additional authors not shown)

Abstract: Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the st… ▽ More Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt-Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly-formed stellar clump mass functions show more significant variation (difference by up to a factor of ~3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low density region, and between more diffusive and less diffusive schemes in the high density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes. △ Less

Submitted 8 March, 2018; v1 submitted 10 October, 2016; originally announced October 2016.

Comments: 28 pages, 35 figures, Accepted for publication in the Astrophysical Journal, Image resolution greatly reduced, High-resolution version of this article is available at http://www.jihoonkim.org/agora/AGORA_Paper4_draft.pdf, The first paper of the AGORA Initiative is at http://adsabs.harvard.edu/abs/2014ApJS..210...14K, More information on AGORA is at http://www.AGORAsimulations.org/

Journal ref: Astrophysical Journal, 833, 202, 2016

MUFASA: Galaxy star formation, gas, and metal properties across cosmic time

Authors: Romeel Davé, Mika H. Rafieferantsoa, Robert J. Thompson, Philip F. Hopkins

Abstract: We examine galaxy star formation rates (SFRs), metallicities, and gas contents predicted by the MUFASA cosmological hydrodynamic simulations, which employ meshless hydrodynamics and novel feedback prescriptions that yield a good match to observed galaxy stellar mass assembly. We combine 50, 25, and 12.5 Mpc/h boxes with a quarter billion particles each to show that MUFASA broadly reproduces a wide… ▽ More We examine galaxy star formation rates (SFRs), metallicities, and gas contents predicted by the MUFASA cosmological hydrodynamic simulations, which employ meshless hydrodynamics and novel feedback prescriptions that yield a good match to observed galaxy stellar mass assembly. We combine 50, 25, and 12.5 Mpc/h boxes with a quarter billion particles each to show that MUFASA broadly reproduces a wide range of relevant observations, including SFR and specific SFR functions, the mass-metallicity relation, HI and H2 fractions, HI (21 cm) and CO luminosity functions, and cosmic gas density evolution. There are mild but significant discrepancies, such as too many high-SFR galaxies, overly metal-rich and HI-poor galaxies at M*>10^{10} Mo, and sSFRs that are too low at z~1-2. The HI mass function increases by x2 out to z~1 then steepens to higher redshifts, while the CO luminosity function computed using the Narayanan et al. conversion factor shows a rapid increase of CO-bright galaxies out to z~2 in accord with data. Omega_HI and Omega_H2 both scale roughly as (1+z)^0.7 out to z~3, comparable to the rise in HI and H2 fractions. MUFASA galaxies with high SFR at a given M* have lower metallicities and higher HI and H2 fractions, following observed trends; we make quantitative predictions for how fluctuations in the baryon cycle drive correlated scatter around galaxy scaling relations. Most of these trends are well converged with numerical resolution. These successes highlight MUFASA as a viable platform to study many facets of cosmological galaxy evolution. △ Less

Submitted 5 October, 2016; originally announced October 2016.

Comments: 19 pages, MNRAS submitted

gPhoton: The GALEX Photon Data Archive

Authors: Chase Million, Scott W. Fleming, Bernie Shiao, Mark Seibert, Parke Loyd, Michael Tucker, Myron Smith, Randy Thompson, Richard L. White

Abstract: gPhoton is a new database product and software package that enables analysis of GALEX ultraviolet data at the photon level. The project's stand-alone, pure-Python calibration pipeline reproduces the functionality of the original mission pipeline to reduce raw spacecraft data to lists of time-tagged, sky-projected photons, which are then hosted in a publicly available database by the Mikulski Archi… ▽ More gPhoton is a new database product and software package that enables analysis of GALEX ultraviolet data at the photon level. The project's stand-alone, pure-Python calibration pipeline reproduces the functionality of the original mission pipeline to reduce raw spacecraft data to lists of time-tagged, sky-projected photons, which are then hosted in a publicly available database by the Mikulski Archive at Space Telescope (MAST). This database contains approximately 130 terabytes of data describing approximately 1.1 trillion sky-projected events with a timestamp resolution of five milliseconds. A handful of Python and command line modules serve as a front-end to interact with the database and to generate calibrated light curves and images from the photon-level data at user-defined temporal and spatial scales. The gPhoton software and source code are in active development and publicly available under a permissive license. We describe the motivation, design, and implementation of the calibration pipeline, database, and tools, with emphasis on divergence from prior work, as well as challenges created by the large data volume. We summarize the astrometric and photometric performance of gPhoton relative to the original mission pipeline. For a brief example of short time domain science capabilities enabled by gPhoton, we show new flares from the known M dwarf flare star CR Draconis. The gPhoton software has permanent object identifiers with the ASCL (ascl:1603.004) and DOI (doi:10.17909/T9CC7G). This paper describes the software as of version v1.27.2. △ Less

Submitted 29 September, 2016; originally announced September 2016.

Comments: 17 pages, 13 figures, software repo at https://github.com/cmillion/gPhoton

The Minimum Halo Mass for Star Formation at z = 6 - 8

Authors: K. Finlator, M. K. M. Prescott, B. D. Oppenheimer, R. Davé, E. Zackrisson, R. C. Livermore, S. L. Finkelstein, R. Thompson, S. Huang

Abstract: Recent analysis of strongly-lensed sources in the Hubble Frontier Fields indicates that the rest-frame UV luminosity function of galaxies at $z=$6--8 rises as a power law down to $M_\mathrm{UV}=-15$, and possibly as faint as -12.5. We use predictions from a cosmological radiation hydrodynamic simulation to map these luminosities onto physical space, constraining the minimum dark matter halo mass a… ▽ More Recent analysis of strongly-lensed sources in the Hubble Frontier Fields indicates that the rest-frame UV luminosity function of galaxies at $z=$6--8 rises as a power law down to $M_\mathrm{UV}=-15$, and possibly as faint as -12.5. We use predictions from a cosmological radiation hydrodynamic simulation to map these luminosities onto physical space, constraining the minimum dark matter halo mass and stellar mass that the Frontier Fields probe. While previously-published theoretical studies have suggested or assumed that early star formation was suppressed in halos less massive than $10^9$--$10^{11} M_\odot$, we find that recent observations demand vigorous star formation in halos at least as massive as (3.1, 5.6, 10.5)$\times10^9 M_\odot$ at $z=(6,7,8)$. Likewise, we find that Frontier Fields observations probe down to stellar masses of (8.1, 18, 32)$\times10^6 M_\odot$; that is, they are observing the likely progenitors of analogues to Local Group dwarfs such as Pegasus and M32. Our simulations yield somewhat different constraints than two complementary models that have been invoked in similar analyses, emphasizing the need for further observational constraints on the galaxy-halo connection. △ Less

Submitted 20 September, 2016; originally announced September 2016.

Comments: 8 pages, 5 figures, re-submitted to MNRAS after incorporating referee's comments

A Machine Learning Classifier for Fast Radio Burst Detection at the VLBA

Authors: Kiri L. Wagstaff, Benyang Tang, David R. Thompson, Shakeh Khudikyan, Jane Wyngaard, Adam T. Deller, Divya Palaniswamy, Steven J. Tingay, Randall B. Wayth

Abstract: Time domain radio astronomy observing campaigns frequently generate large volumes of data. Our goal is to develop automated methods that can identify events of interest buried within the larger data stream. The V-FASTR fast transient system was designed to detect rare fast radio bursts (FRBs) within data collected by the Very Long Baseline Array. The resulting event candidates constitute a signifi… ▽ More Time domain radio astronomy observing campaigns frequently generate large volumes of data. Our goal is to develop automated methods that can identify events of interest buried within the larger data stream. The V-FASTR fast transient system was designed to detect rare fast radio bursts (FRBs) within data collected by the Very Long Baseline Array. The resulting event candidates constitute a significant burden in terms of subsequent human reviewing time. We have trained and deployed a machine learning classifier that marks each candidate detection as a pulse from a known pulsar, an artifact due to radio frequency interference, or a potential new discovery. The classifier maintains high reliability by restricting its predictions to those with at least 90% confidence. We have also implemented several efficiency and usability improvements to the V-FASTR web-based candidate review system. Overall, we found that time spent reviewing decreased and the fraction of interesting candidates increased. The classifier now classifies (and therefore filters) 80-90% of the candidates, with an accuracy greater than 98%, leaving only the 10-20% most promising candidates to be reviewed by humans. △ Less

Submitted 28 June, 2016; originally announced June 2016.

Comments: Published in PASP

Limits on Fast Radio Bursts from Four Years of the V-FASTR Experiment

Authors: Sarah Burke-Spolaor, Cathryn M. Trott, Walter F. Brisken, Adam T. Deller, Walid A. Majid, Divya Palaniswamy, David R. Thompson, Steven J. Tingay, Kiri L. Wagstaff, Randall B. Wayth

Abstract: The V-FASTR experiment on the Very Long Baseline Array was designed to detect dispersed pulses of milliseconds duration, such as fast radio bursts (FRBs). We use all V-FASTR data through February 2015 to report V-FASTR's upper limits on the rates of FRBs, and compare these with re-derived rates from Parkes FRB detection experiments. V-FASTR's operation at lambda=20 cm allows direct comparison with… ▽ More The V-FASTR experiment on the Very Long Baseline Array was designed to detect dispersed pulses of milliseconds duration, such as fast radio bursts (FRBs). We use all V-FASTR data through February 2015 to report V-FASTR's upper limits on the rates of FRBs, and compare these with re-derived rates from Parkes FRB detection experiments. V-FASTR's operation at lambda=20 cm allows direct comparison with the 20 cm Parkes rate, and we derive a power-law limit of γ<-0.4 (95% confidence limit) on the index of FRB source counts, N(>S)\propto S^γ. Using the previously measured FRB rate and the unprecedented amount of survey time spent searching for FRBs at a large range of wavelengths (0.3 cm > λ> 90 cm), we also place frequency-dependent limits on the spectral distribution of FRBs. The most constraining frequencies place two-point spectral index limits of α_{20cm}^{4cm} < 5.8 and α_{90cm}^{20cm} > -7.6, where fluence F \propto f^αif we assume true the burst rate reported by Champion et al. (2016) of R(F~0.6 Jy ms) = 7 x 10^3 sky^{-1} day^{-1} (for bursts of ~3 ms duration). This upper limit on αsuggests that if FRBs are extragalactic but non-cosmological, that on average they are not experiencing excessive free-free absorption due to a medium with high optical depth (assuming temperature ~8,000 K), which excessively invert their low-frequency spectrum. This in turn implies that the dispersion of FRBs arises in either or both of the intergalactic medium or the host galaxy, rather than from the source itself. △ Less

Submitted 23 May, 2016; originally announced May 2016.

Comments: Accepted for publication in ApJ

MUFASA: Galaxy Formation Simulations With Meshless Hydrodynamics

Authors: Romeel Davé, Robert J. Thompson, Philip F. Hopkins

Abstract: We present the MUFASA suite of cosmological hydrodynamic simulations, which employs the GIZMO meshless finite mass (MFM) code including H2-based star formation, nine-element chemical evolution, two-phase kinetic outflows following scalings from the Feedback in Realistic Environments zoom simulations, and evolving halo mass-based quenching. Our fiducial (50 Mpc/h)^3 volume is evolved to z=0 with a… ▽ More We present the MUFASA suite of cosmological hydrodynamic simulations, which employs the GIZMO meshless finite mass (MFM) code including H2-based star formation, nine-element chemical evolution, two-phase kinetic outflows following scalings from the Feedback in Realistic Environments zoom simulations, and evolving halo mass-based quenching. Our fiducial (50 Mpc/h)^3 volume is evolved to z=0 with a quarter billion particles, The predicted galaxy stellar mass functions (GSMF) reproduce observations from z=4-0 to <1.2sigma in cosmic variance, providing an unprecedented match to this key diagnostic. The cosmic star formation history and stellar mass growth show general agreement with data, with a strong archaeological downsizing trend such that dwarf galaxies form the majority of their stars after z~1. We run 25 Mpc/h and 12.5 Mpc/h volumes to z=2 with identical feedback prescriptions, the latter resolving all hydrogen-cooling halos, and the three runs display fair resolution convergence. The specific star formation rates broadly agree with data at z=0, but are underpredicted at z~2 by a factor of three, re-emphasizing a longstanding puzzle in galaxy evolution models. We compare runs using MFM and two flavours of Smoothed Particle Hydrodynamics, and show that the GSMF is sensitive to hydrodynamics methodology at the ~x2 level, which is sub-dominant to choices for parameterising feedback. △ Less

Submitted 4 August, 2016; v1 submitted 5 April, 2016; originally announced April 2016.

Comments: 22 pages, MNRAS, accepted

The Soft, Fluctuating UVB at $z\sim6$ as Traced by C IV, SiIV, and CII

Authors: K. Finlator, B. D. Oppenheimer, R. Davé, E. Zackrisson, R. Thompson, S. Huang

Abstract: The sources that drove cosmological reionization left clues regarding their identity in the slope and inhomogeneity of the ultraviolet ionizing background (UVB): Bright quasars (QSOs) generate a hard UVB with predominantly large-scale fluctuations while Population II stars generate a softer one with smaller-scale fluctuations. Metal absorbers probe the UVB's slope because different ions are sensit… ▽ More The sources that drove cosmological reionization left clues regarding their identity in the slope and inhomogeneity of the ultraviolet ionizing background (UVB): Bright quasars (QSOs) generate a hard UVB with predominantly large-scale fluctuations while Population II stars generate a softer one with smaller-scale fluctuations. Metal absorbers probe the UVB's slope because different ions are sensitive to different energies. Likewise, they probe spatial fluctuations because they originate in regions where a galaxy-driven UVB is harder and more intense. We take a first step towards studying the reionization-epoch UVB's slope and inhomogeneity by comparing observations of 12 metal absorbers at $z\sim6$ versus predictions from a cosmological hydrodynamic simulation using three different UVBs: a soft, spatially-inhomogeneous "galaxies+QSOs" UVB; a homogeneous "galaxies+QSOs" UVB (Haardt & Madau 2012); and a QSOs-only model. All UVBs reproduce the observed column density distributions of CII, SiIV, and CIV reasonably well although high-column, high-ionization absorbers are underproduced, reflecting numerical limitations. With upper limits treated as detections, only a soft, fluctuating UVB reproduces both the observed SiIV/CIV and CII/CIV distributions. The QSOs-only UVB overpredicts both CIV/CII and CIV/SiIV, indicating that it is too hard. The Haardt & Madau (2012) UVB underpredicts CIV/SiIV, suggesting that it lacks amplifications near galaxies. Hence current observations prefer a soft, fluctuating UVB as expected from a predominantly Population II background although they cannot rule out a harder one. Future observations probing a factor of two deeper in metal column density will distinguish between the soft, fluctuating and QSOs-only UVBs. △ Less

Submitted 21 March, 2016; originally announced March 2016.

Comments: 13 pages, 10 figures, re-submitted to MNRAS after incorporating referee's suggestions

Ground-based near-UV observations of 15 transiting exoplanets: Constraints on their atmospheres and no evidence for asymmetrical transits

Authors: Jake D. Turner, Kyle A. Pearson, Lauren I. Biddle, Brianna M. Smart, Robert T. Zellem, Johanna K. Teske, Kevin K. Hardegree-Ullman, Caitlin C. Griffith, Robin M. Leiter, Ian T. Cates, Megan N. Nieberding, Carter-Thaxton W. Smith, Robert M. Thompson, Ryan Hofmann, Michael P. Berube, Chi H. Nguyen, Lindsay C. Small, Blythe C. Guvenen, Logan Richardson, Allison McGraw, Brandon Raphael, Benjamin E. Crawford, Amy N. Robertson, Ryan Tombleson, Timothy M. Carleton , et al. (15 additional authors not shown)

Abstract: Transits of exoplanets observed in the near-UV have been used to study the scattering properties of their atmospheres and possible star-planet interactions. We observed the primary transits of 15 exoplanets (CoRoT-1b, GJ436b, HAT-P-1b, HAT-P-13b, HAT-P-16b, HAT-P-22b, TrES-2b, TrES-4b, WASP-1b, WASP-12b, WASP-33b, WASP-36b, WASP-44b, WASP-48b, and WASP-77Ab) in the near-UV and several optical phot… ▽ More Transits of exoplanets observed in the near-UV have been used to study the scattering properties of their atmospheres and possible star-planet interactions. We observed the primary transits of 15 exoplanets (CoRoT-1b, GJ436b, HAT-P-1b, HAT-P-13b, HAT-P-16b, HAT-P-22b, TrES-2b, TrES-4b, WASP-1b, WASP-12b, WASP-33b, WASP-36b, WASP-44b, WASP-48b, and WASP-77Ab) in the near-UV and several optical photometric bands to update their planetary parameters, ephemerides, search for a wavelength dependence in their transit depths to constrain their atmospheres, and determine if asymmetries are visible in their light curves. Here we present the first ground-based near-UV light curves for 12 of the targets (CoRoT-1b, GJ436b, HAT-P-1b, HAT-P-13b, HAT-P-22b, TrES-2b, TrES-4b, WASP-1b, WASP-33b, WASP-36b, WASP-48b, and WASP-77Ab). We find that none of the near-UV transits exhibit any non-spherical asymmetries, this result is consistent with recent theoretical predictions by Ben-Jaffel et al. and Turner et al. The multi-wavelength photometry indicates a constant transit depth from near-UV to optical wavelengths in 10 targets (suggestive of clouds), and a varying transit depth with wavelength in 5 targets (hinting at Rayleigh or aerosol scattering in their atmospheres). We also present the first detection of a smaller near-UV transit depth than that measured in the optical in WASP-1b and a possible opacity source that can cause such radius variations is currently unknown. WASP-36b also exhibits a smaller near-UV transit depth at 2.6$σ$. Further observations are encouraged to confirm the transit depth variations seen in this study. △ Less

Submitted 8 March, 2016; originally announced March 2016.

Comments: 35 pages, 16 tables, 10 figures. Accepted for publication in MNRAS (March 8, 2016). Light curves are available online. Updated EXOMOP (transit modeling software) and is also available online

Fundamental Constants as Monitors of the Universe

Authors: Rodger I. Thompson

Abstract: Astronomical observations have a unique ability to determine the laws of physics at distant times in the universe. They, therefore, have particular relevance in answering the basic question as to whether the laws of physics are invariant with time. The dimesionless fundamental constants, such as the proton to electron mass ratio and the fine structure constant are key elements in the investigation… ▽ More Astronomical observations have a unique ability to determine the laws of physics at distant times in the universe. They, therefore, have particular relevance in answering the basic question as to whether the laws of physics are invariant with time. The dimesionless fundamental constants, such as the proton to electron mass ratio and the fine structure constant are key elements in the investigation. If they vary with time then the answer is clearly that the laws of physics are not invariant with time and significant new physics must be developed to describe the universe. Limits on their variance, on the other hand, constrains the parameter space available to new physics that requires a variation with time of basic physical law. There are now observational constraints on the time variation of the proton to electron mass ratio mu at the 1.E-7 level. Constraints on the variation of the fine structure constant alpha are less rigorous, 1E-5, but are imposed at higher redshift. The implications of these limits on new cosmologies that require rolling scalar fields has already had its first investigations. Here we address the implications on basic particle physics. The proton to electron mass ratio is obviously dependent on the particle physics parameters that set the mass of the proton and the electron. To first order the ratio is dependent on a combination of the Quantum Chromodynamic scale, the Yukawa couplings, and the Higgs Vacuum Expectation Value. Here that relationship is quantitative defined for the first time. When coupled with previous determinations of the relation of the fine structure constant to the same parameters two constraints exist on the fractional variation of these parameters with time. A third independent constraint involving only the three parameters could set the stage for constraints on their individual fractional variation. △ Less

Submitted 9 February, 2016; originally announced February 2016.

Comments: To appear in the Proceedings of the 14th Grossmann Conference. Seven pages, two figures

The AIMSS Project III: The Stellar Populations of Compact Stellar Systems

Authors: Joachim Janz, Mark A. Norris, Duncan A. Forbes, Avon Huxor, Aaron J. Romanowsky, Matthias J. Frank, Carlos G. Escudero, Favio R. Faifer, Juan Carlos Forte, Sheila J. Kannappan, Claudia Maraston, Jean P. Brodie, Jay Strader, Bradley R. Thompson

Abstract: In recent years, a growing zoo of compact stellar systems (CSSs) have been found whose physical properties (mass, size, velocity dispersion) place them between classical globular clusters (GCs) and true galaxies, leading to debates about their nature. Here we present results using a so far underutilised discriminant, their stellar population properties. Based on new spectroscopy from 8-10m telesco… ▽ More In recent years, a growing zoo of compact stellar systems (CSSs) have been found whose physical properties (mass, size, velocity dispersion) place them between classical globular clusters (GCs) and true galaxies, leading to debates about their nature. Here we present results using a so far underutilised discriminant, their stellar population properties. Based on new spectroscopy from 8-10m telescopes, we derive ages, metallicities, and [α/Fe] of 29 CSSs. These range from GCs with sizes of merely a few parsec to compact ellipticals larger than M32. Together with a literature compilation, this provides a panoramic view of the stellar population characteristics of early-type systems. We find that the CSSs are predominantly more metal rich than typical galaxies at the same stellar mass. At high mass, the compact ellipticals (cEs) depart from the mass-metallicity relation of massive early-type galaxies, which forms a continuous sequence with dwarf galaxies. At lower mass, the metallicity distribution of ultra-compact dwarfs (UCDs) changes at a few times $10^7$ M$_{\odot}$, which roughly coincides with the mass where luminosity function arguments previously suggested the GC population ends. The highest metallicities in CSSs are paralleled only by those of dwarf galaxy nuclei and the central parts of massive early types. These findings can be interpreted as CSSs previously being more massive and undergoing tidal interactions to obtain their current mass and compact size. Such an interpretation is supported by CSSs with direct evidence for tidal stripping, and by an examination of the CSS internal escape velocities. △ Less

Submitted 3 December, 2015; v1 submitted 10 November, 2015; originally announced November 2015.

Comments: 18 pages, 10 figures, accepted for publication in MNRAS

The Formation of Submillimetre-Bright Galaxies from Gas Infall over a Billion Years

Authors: Desika Narayanan, Matthew Turk, Robert Feldmann, Thomas Robitaille, Philip Hopkins, Robert Thompson, Christopher Hayward, David Ball, Claude-Andre Faucher-Giguere, Dusan Keres

Abstract: Submillimetre-luminous galaxies at high-redshift are the most luminous, heavily star-forming galaxies in the Universe, and are characterised by prodigious emission in the far-infrared at 850 microns (S850 > 5 mJy). They reside in halos ~ 10^13Msun, have low gas fractions compared to main sequence disks at a comparable redshift, trace complex environments, and are not easily observable at optical w… ▽ More Submillimetre-luminous galaxies at high-redshift are the most luminous, heavily star-forming galaxies in the Universe, and are characterised by prodigious emission in the far-infrared at 850 microns (S850 > 5 mJy). They reside in halos ~ 10^13Msun, have low gas fractions compared to main sequence disks at a comparable redshift, trace complex environments, and are not easily observable at optical wavelengths. Their physical origin remains unclear. Simulations have been able to form galaxies with the requisite luminosities, but have otherwise been unable to simultaneously match the stellar masses, star formation rates, gas fractions and environments. Here we report a cosmological hydrodynamic galaxy formation simulation that is able to form a submillimetre galaxy which simultaneously satisfies the broad range of observed physical constraints. We find that groups of galaxies residing in massive dark matter halos have rising star formation histories that peak at collective rates ~ 500-1000 Msun/yr at z=2-3, by which time the interstellar medium is sufficiently enriched with metals that the region may be observed as a submillimetre-selected system. The intense star formation rates are fueled in part by a reservoir gas supply enabled by stellar feedback at earlier times, not through major mergers. With a duty cycle of nearly a gigayear, our simulations show that the submillimetre-luminous phase of high-z galaxies is a drawn out one that is associated with significant mass buildup in early Universe proto-clusters, and that many submillimetre-luminous galaxies are actually composed of numerous unresolved components (for which there is some observational evidence). △ Less

Submitted 21 September, 2015; originally announced September 2015.

Comments: Nature (in press), To be published September 24th, 2015, PDF version includes Methods and Supplementary Figures

Identifying mergers using non-parametric morphological classification at high redshifts

Authors: Robert Thompson, Romeel Davé, Shuiyao Huang, Neal Katz

Abstract: We investigate the time evolution of non-parametric morphological quantities and their relationship to major mergers between $4\geq z \geq 2$ in high-resolution cosmological zoom simulations of disk galaxies that implement kinetic wind feedback, $H_2$-based star formation, and minimal ISM pressurisation. We show that the resulting galaxies broadly match basic observed physical properties of… ▽ More We investigate the time evolution of non-parametric morphological quantities and their relationship to major mergers between $4\geq z \geq 2$ in high-resolution cosmological zoom simulations of disk galaxies that implement kinetic wind feedback, $H_2$-based star formation, and minimal ISM pressurisation. We show that the resulting galaxies broadly match basic observed physical properties of $z\sim 2$ objects. We measure the galaxies' concentrations ($C$), asymmetries ($A$), and $Gini$ ($G$) and $M_{20}$ coefficients, and correlate these with major merger events identified from the mass growth history. We find that high values of asymmetry provide the best indicator for identifying major mergers of $>1:4$ mass ratio within our sample, with $Gini$-$M_{20}\,$ merger classification only as effective for face-on systems and much less effective for edge-on or randomly-oriented galaxies. The canonical asymmetry cut of $A\geq0.35$, however, is only able to correctly identify major mergers $\sim 10\%$ of the time, while a higher cut of $A\geq 0.8$ more efficiently picks out mergers at this epoch. We further examine the temporal correlation between morphological statistics and mergers, and show that for randomly-oriented galaxies, half the galaxies with $A\geq0.8$ undergo a merger within $\pm0.2\,{\rm Gyr}$, whereas $Gini$-$M_{20}\,$ identification only identifies about a third correctly. The fraction improves further using $A\geq 1.5$, but about the half the mergers are missed by this stringent cut. △ Less

Submitted 8 August, 2015; originally announced August 2015.

Comments: 19 pages, 19 figures, submitted

Simulator of Galaxy Millimeter/Submillimeter Emission (SIGAME): The [CII]-SFR Relationship of Massive z=2 Main Sequence Galaxies

Authors: Karen P. Olsen, Thomas R. Greve, Desika Narayanan, Robert Thompson, Sune Toft, Christian Brinch

Abstract: We present SÍGAME simulations of the [CII]157.7$μ$ fine structure line emission from cosmological smoothed particle hydrodynamics (SPH) simulations of seven main sequence galaxies at z=2. Using sub-grid physics prescriptions the gas in our simulations is modeled as a multi-phased interstellar medium (ISM) comprised of molecular gas residing in giant molecular clouds, an atomic gas phase associated… ▽ More We present SÍGAME simulations of the [CII]157.7$μ$ fine structure line emission from cosmological smoothed particle hydrodynamics (SPH) simulations of seven main sequence galaxies at z=2. Using sub-grid physics prescriptions the gas in our simulations is modeled as a multi-phased interstellar medium (ISM) comprised of molecular gas residing in giant molecular clouds, an atomic gas phase associated with photo-dissociation regions (PDRs) at the cloud surfaces, and a diffuse, ionized gas phase. Adopting logotropic cloud density profiles and accounting for heating by the local FUV radiation field and cosmic rays by scaling both with local star formation rate (SFR) volume density, we calculate the [CII] emission using a photon escape probability formalism. The [CII] emission peaks in the central $\lesssim$1 kpc of our galaxies as do the SFR radial profiles, with most [CII] ($\gtrsim$70%) originating in the molecular gas phase, whereas further out ($\gtrsim$2 kpc), the atomic/PDR gas dominates ($\gtrsim$90%) the [CII] emission, no longer tracing on-going star formation. Throughout, the ionized gas contribution is negligible ($\lesssim$3%). The [CII] luminosity vs. SFR ([CII]-SFR) relationship, integrated as well as spatially resolved (on scales of 1 kpc), delineated by our simulated galaxies is in good agreement with the corresponding relations observed locally and at high redshifts. In our simulations, the molecular gas dominates the [CII] budget at SFR$\gtrsim$20 M$_{\odot}$yr$^{-1}$ (SFR density $\gtrsim$0.5 M$_{\odot}$yr$^{-1}$kpc$^{-2}$), while atomic/PDR gas takes over at lower SFRs, suggesting a picture in which [CII] predominantly traces the molecular gas in high-density/pressure regions where star formation is on-going, and otherwise reveals the atomic/PDR gas phase. △ Less

Submitted 23 November, 2015; v1 submitted 1 July, 2015; originally announced July 2015.

Comments: 23 pages, 14 figures, 2 tables; Accepted for publication in the Astrophysical Journal

The Reionization of Carbon

Authors: K. Finlator, R. Thompson, S. Huang, R. Davé, E. Zackrisson, B. D. Oppenheimer

Abstract: Observations suggest that CII was more abundant than CIV in the intergalactic medium towards the end of the hydrogen reionization epoch. This transition provides a unique opportunity to study the enrichment history of intergalactic gas and the growth of the ionizing background (UVB) at early times. We study how carbon absorption evolves from z=10-5 using a cosmological hydrodynamic simulation that… ▽ More Observations suggest that CII was more abundant than CIV in the intergalactic medium towards the end of the hydrogen reionization epoch. This transition provides a unique opportunity to study the enrichment history of intergalactic gas and the growth of the ionizing background (UVB) at early times. We study how carbon absorption evolves from z=10-5 using a cosmological hydrodynamic simulation that includes a self-consistent multifrequency UVB as well as a well-constrained model for galactic outflows to disperse metals. Our predicted UVB is within 2-4 times that of Haardt & Madau (2012), which is fair agreement given the uncertainties. Nonetheless, we use a calibration in post-processing to account for Lyman-alpha forest measurements while preserving the predicted spectral slope and inhomogeneity. The UVB fluctuates spatially in such a way that it always exceeds the volume average in regions where metals are found. This implies both that a spatially-uniform UVB is a poor approximation and that metal absorption is not sensitive to the epoch when HII regions overlap globally even at column densites of 10^{12} cm^{-2}. We find, consistent with observations, that the CII mass fraction drops to low redshift while CIV rises owing the combined effects of a growing UVB and continued addition of carbon in low-density regions. This is mimicked in absorption statistics, which broadly agree with observations at z=6-3 while predicting that the absorber column density distributions rise steeply to the lowest observable columns. Our model reproduces the large observed scatter in the number of low-ionization absorbers per sightline, implying that the scatter does not indicate a partially-neutral Universe at z=6. △ Less

Submitted 8 January, 2015; v1 submitted 15 December, 2014; originally announced December 2014.

Comments: 16 pages, 14 figures, accepted to MNRAS

The rise and fall of a challenger: the Bullet Cluster in $Λ$ Cold Dark Matter simulations

Authors: Robert Thompson, Romeel Davé, Kentaro Nagamine

Abstract: The Bullet Cluster has provided some of the best evidence for the $Λ$ cold dark matter ($Λ\rm{CDM}$) model via direct empirical proof of the existence of collisionless dark matter, while posing a serious challenge owing to the unusually high inferred pairwise velocities of its progenitor clusters. Here we investigate the probability of finding such a high-velocity pair in large-volume N-body simul… ▽ More The Bullet Cluster has provided some of the best evidence for the $Λ$ cold dark matter ($Λ\rm{CDM}$) model via direct empirical proof of the existence of collisionless dark matter, while posing a serious challenge owing to the unusually high inferred pairwise velocities of its progenitor clusters. Here we investigate the probability of finding such a high-velocity pair in large-volume N-body simulations, particularly focusing on differences between halo finding algorithms. We find that algorithms that do not account for the kinematics of infalling groups yield vastly different statistics and probabilities. When employing the ROCKSTAR halo finder that considers particle velocities, we find numerous Bullet-like pair candidates that closely match not only the high pairwise velocity, but also the mass, mass ratio, separation distance, and collision angle of the initial conditions that have been shown to produce the Bullet Cluster in non-cosmological hydrodynamic simulations. The probability of finding a high pairwise velocity pair among haloes with $M_{\rm halo}\geq10^{14} M_{\odot}$ is $4.6\times 10^{-4}$ using ROCKSTAR, while it is $\approx 34\times$ lower using a friends-of-friends (FOF) based approach as in previous studies. This is because the typical spatial extent of Bullet progenitors is such that FOF tends to group them into a single halo despite clearly distinct kinematics. Further requiring an appropriately high average mass among the two progenitors, we find the comoving number density of potential Bullet-like candidates to be on the order of $\approx10^{-10} \,{\rm Mpc}^{-3}$. Our findings suggest that $Λ\rm{CDM}$ straightforwardly produces massive, high relative velocity halo pairs analogous to Bullet Cluster progenitors, and hence the Bullet Cluster does not present a challenge to the $Λ\rm{CDM}$ model. △ Less

Submitted 29 June, 2015; v1 submitted 27 October, 2014; originally announced October 2014.

Comments: 10 pages, 6 figures, 3 tables, accepted to MNRAS

Iris: an Extensible Application for Building and Analyzing Spectral Energy Distributions

Authors: Omar Laurino, Jamie Budynkiewicz, Raffaele D'Abrusco, Nina Bonaventura, Ivo Busko, Mark Cresitello-Dittmar, Stephen M. Doe, Rick Ebert, Janet D. Evans, Patrick Norris, Olga Pevunova, Brian Refsdal, Brian Thomas, Randy Thompson

Abstract: Iris is an extensible application that provides astronomers with a user-friendly interface capable of ingesting broad-band data from many different sources in order to build, explore, and model spectral energy distributions (SEDs). Iris takes advantage of the standards defined by the International Virtual Observatory Alliance, but hides the technicalities of such standards by implementing differen… ▽ More Iris is an extensible application that provides astronomers with a user-friendly interface capable of ingesting broad-band data from many different sources in order to build, explore, and model spectral energy distributions (SEDs). Iris takes advantage of the standards defined by the International Virtual Observatory Alliance, but hides the technicalities of such standards by implementing different layers of abstraction on top of them. Such intermediate layers provide hooks that users and developers can exploit in order to extend the capabilities provided by Iris. For instance, custom Python models can be combined in arbitrary ways with the Iris built-in models or with other custom functions. As such, Iris offers a platform for the development and integration of SED data, services, and applications, either from the user's system or from the web. In this paper we describe the built-in features provided by Iris for building and analyzing SEDs. We also explore in some detail the Iris framework and software development kit, showing how astronomers and software developers can plug their code into an integrated SED analysis environment. △ Less

Submitted 25 July, 2014; originally announced July 2014.

Comments: 18 pages, 8 figures, accepted for publication in Astronomy & Computing

The Magellan Uniform Survey of Damped Lyman alpha Systems II: Paucity of Strong Molecular Hydrogen Absorption

Authors: Regina A. Jorgenson, Michael T. Murphy, Rodger Thompson, Robert F. Carswell

Abstract: We present the first large, blind and uniformly selected survey for molecular hydrogen (H_2) in damped Lyman alpha systems (DLAs) with moderate-to-high resolution spectra. 86 DLAs were searched for absorption in the many Lyman and Werner H_2 transitions, with approximately 79% completeness for H_2 column densities above N(H_2)=10^17.5 cm^-2 for an assumed Doppler broadening parameter b=2 km s^-1.… ▽ More We present the first large, blind and uniformly selected survey for molecular hydrogen (H_2) in damped Lyman alpha systems (DLAs) with moderate-to-high resolution spectra. 86 DLAs were searched for absorption in the many Lyman and Werner H_2 transitions, with approximately 79% completeness for H_2 column densities above N(H_2)=10^17.5 cm^-2 for an assumed Doppler broadening parameter b=2 km s^-1. Only a single strong H_2 absorber was found -- a system detected previously in VLT/UVES spectra. Given our distribution of N(H_2) upper limits, this ~1% detection rate is smaller than expected from previous surveys at 99.8% confidence. Assuming the N(H_2) distribution shape from previous surveys, our detection rate implies a covering factor of ~1% for N(H_2) >= 10^17.5 cm^-2 gas in DLAs (<6% at 95% confidence). We obtained new Magellan/MagE spectra for 53 DLAs; 8 km s^-1 -resolution spectra were available for 27 DLAs. MagE's moderate resolution (~71 km s^-1) yields weaker N(H_2) upper limits and makes them dependent on the assumed Doppler parameter. For example, half the (relevant) previous H_2 detections have N(H_2) >= 10^18.1 cm^-2, a factor of just 3 higher than our median upper limit. Nevertheless, several tests suggest our upper limits are accurate, and they would need to be increased by 1.8 dex to bring our detection rate within 95% confidence of previous surveys. △ Less

Submitted 3 July, 2014; originally announced July 2014.

Comments: 21 pages, 20 figures, Accepted by MNRAS

Warm Ice Giant GJ 3470b. II Revised Planetary and Stellar Parameters from Optical to Near-infrared Transit Photometry

Authors: Lauren I. Biddle, Kyle A. Pearson, Ian J. M. Crossfield, Benjamin J. Fulton, Simona Ciceri, Jason Eastman, Travis Barman, Andrew W. Mann, Gregory W. Henry, Andrew W. Howard, Michael H. Williamson, Evan Sinukoff, Diana Dragomir, Laura Vican, Luigi Mancini, John Southworth, Adam Greenberg, Jake D. Turner, Robert Thompson, Brian W. Taylor, Stephen E. Levine, Matthew W. Webber

Abstract: It is important to explore the diversity of characteristics of low-mass, low-density planets to understand the nature and evolution of this class of planets. We present a homogeneous analysis of 12 new and 9 previously published broadband photometric observations of the Uranus-sized extrasolar planet GJ 3470b, which belongs to the growing sample of sub-Jovian bodies orbiting M dwarfs. The consiste… ▽ More It is important to explore the diversity of characteristics of low-mass, low-density planets to understand the nature and evolution of this class of planets. We present a homogeneous analysis of 12 new and 9 previously published broadband photometric observations of the Uranus-sized extrasolar planet GJ 3470b, which belongs to the growing sample of sub-Jovian bodies orbiting M dwarfs. The consistency of our analysis explains some of the discrepancies between previously published results and provides updated constraints on the planetary parameters. Our data are also consistent with previous transit observations of this system. We also provide new spectroscopic measurements of GJ 3470 from 0.33 to 2.42 $μ$$m$ to aid our analysis. We find $R_{\star}$ = 0.48$\pm$0.04 $R_{\odot}$, $M_{\star}$ = 0.51$\pm$0.06 $M_{\odot}$, and $T_{\rm eff}$ = 3652$\pm$50 K for GJ 3470, along with a rotation period of $20.70\pm{0.15}$ d and an R-band amplitude of 0.01 mag, which is small enough that current transit measurements should not be strongly affected by stellar variability. We also present the most precise orbital ephemeris for this system: T$_{o}$ = 2455983.70472$\pm$0.00021 BJD$_{TDB}$, P = 3.3366487$^{+0.0000043}_{-0.0000033}$ d, and we see no evidence for transit timing variations greater than 1 minute. Our reported planet to star radius ratio is 0.07642$\pm$0.00037. The physical parameters of this planet are $R_{p}$ = 3.88$\pm$0.32 $R_{\oplus}$, and $M_{p}$ = 13.73$\pm$1.61 $M_{\oplus}$. Because of our revised stellar parameters, the planetary radius we present is smaller than previously reported values. We also perform a second analysis of the transmission spectrum of the entire ensemble of transit observations to date, supporting the existence of a H$_{2}$ dominated atmosphere exhibiting a strong Rayleigh scattering slope. △ Less

Submitted 14 July, 2014; v1 submitted 24 June, 2014; originally announced June 2014.

Comments: Submitted to MNRAS Dec 2013, Accepted June 2014, 27 Pages, 6 Figures, 5 Tables, Corrected typo

Confronting Cosmology and New Physics with Fundamental Constants

Authors: Rodger I. Thompson

Abstract: The values of the fundamental constants such as $μ= m_P/m_e$, the proton to electron mass ratio and $α$, the fine structure constant, are sensitive to the product $\sqrt{ζ_x^2(w+1)}$ where $ζ_x$ is a coupling constant between a rolling scalar field responsible for the acceleration of the expansion of the universe and the electromagnetic field with x standing for either $μ$ or $α$. The dark energy… ▽ More The values of the fundamental constants such as $μ= m_P/m_e$, the proton to electron mass ratio and $α$, the fine structure constant, are sensitive to the product $\sqrt{ζ_x^2(w+1)}$ where $ζ_x$ is a coupling constant between a rolling scalar field responsible for the acceleration of the expansion of the universe and the electromagnetic field with x standing for either $μ$ or $α$. The dark energy equation of state $w$ can assume values different than $-1$ in cosmologies where the acceleration of the expansion is due to a scalar field. In this case the value of both $μ$ and $α$ changes with time. The values of the fundamental constants, therefore, monitor the equation of state and are a valuable tool for determining $w$ as a function of redshift. In fact the rolling of the fundamental constants is one of the few definitive discriminators between acceleration due to a cosmological constant and acceleration due to a quintessence rolling scalar field. $w$ is often given in parameterized form for comparison with observations. In this manuscript the predicted evolution of $μ$, is calculated for a range of parameterized equation of state models and compared to the observational constraints on $Δμ/ μ$. We find that the current limits on $Δμ/ μ$ place significant constraints on linear equation of state models and on thawing models where $w$ deviates from $-1$ at late times. They also constrain non-dynamical models that have a constant $w$ not equal to $-1$. These constraints are an important compliment to geometric tests of $w$ in that geometric tests are sensitive to the evolution of the universe before the epoch of observation while fundamental constants are sensitive to the evolution of the universe after the observational epoch. Abstract truncated. △ Less

Submitted 17 December, 2013; originally announced December 2013.

Comments: To appear in the conference proceedings of the Sesto Conference on Fundamental Constants and Couplings

WFIRST Science with a Probe Class Mission

Authors: Rodger Thompson, James Green, George Rieke, Dennis Zaritsky. Brant Robertson, Glenn Schneider, Daniel Stark, Buell Jannuzi, Dennis Ebbets, Michael Kaplan, Renee Gracey

Abstract: WFIRST is the highest priority space mission of the Decadal review, however, it is unlikely to begin in this decade primarily due to a anticipated NASA budget that is unlikely to have sufficient resources to fund such a mission. For this reason we present a lower cost mission that accomplishes all of the WFIRST science as described in the Design Reference Mission 1 with a probe class design. This… ▽ More WFIRST is the highest priority space mission of the Decadal review, however, it is unlikely to begin in this decade primarily due to a anticipated NASA budget that is unlikely to have sufficient resources to fund such a mission. For this reason we present a lower cost mission that accomplishes all of the WFIRST science as described in the Design Reference Mission 1 with a probe class design. This is effort is motivated by a desire to begin WFIRST in a timely manner and within a budget that can fit within the assets available to NASA on a realistic basis. The design utilizes dichroics to form four focal planes all having the same field of view to use the majority of available photons from a 1.2 meter telescope. △ Less

Submitted 16 December, 2013; originally announced December 2013.

Comments: 5 pages 3 figures

Dust properties of Lyman break galaxies in cosmological simulations

Authors: Hidenobu Yajima, Kentaro Nagamine, Robert Thompson, Jun-Hwan Choi

Abstract: Recent observations have indicated the existence of dust in high-redshift galaxies, however, the dust properties in them are still unknown. Here we present theoretical constraints on dust properties in Lyman break galaxies (LBGs) at z=3 by post-processing a cosmological smoothed particle hydrodynamics simulation with radiative transfer calculations. We calculate the dust extinction in 2800 dark ma… ▽ More Recent observations have indicated the existence of dust in high-redshift galaxies, however, the dust properties in them are still unknown. Here we present theoretical constraints on dust properties in Lyman break galaxies (LBGs) at z=3 by post-processing a cosmological smoothed particle hydrodynamics simulation with radiative transfer calculations. We calculate the dust extinction in 2800 dark matter halos using the metallicity information of individual gas particles in our simulation. We use only bright galaxies with rest-frame UV magnitude M_1700 < -20 mag, and study the dust size, dust-to-metal mass ratio, and dust composition. From the comparison of calculated color excess between B and V-band (i.e., E(B-V)) and the observations, we constrain the typical dust size, and show that the best-fitting dust grain size is ~ 0.05 micron, which is consistent with the results of theoretical dust models for Type-II supernova. Our simulation with the dust extinction effect can naturally reproduce the observed rest-frame UV luminosity function of LBGs at z=3 without assuming an ad hoc constant extinction value. In addition, in order to reproduce the observed mean E(B-V), we find that the dust-to-metal mass ratio needs to be similar to that of the local galaxies, and that the graphite dust is dominant or at least occupy half of dust mass. △ Less

Submitted 21 January, 2014; v1 submitted 27 September, 2013; originally announced September 2013.

Comments: 13 pages, 13 figures, accepted for publication in MNRAS

Tracing Inflows and Outflows with Absorption Lines in Circumgalactic Gas

Authors: Amanda Brady Ford, Romeel Davé, Benjamin D. Oppenheimer, Neal Katz, Juna A. Kollmeier, Robert Thompson, David H. Weinberg

Abstract: We examine how HI and metal absorption lines within low-redshift galaxy halos trace the dynamical state of circumgalactic gas, using cosmological hydrodynamic simulations that include a well-vetted heuristic model for galactic outflows. We categorize inflowing, outflowing, and ambient gas based on its history and fate as tracked in our simulation. Following our earlier work showing that the ionisa… ▽ More We examine how HI and metal absorption lines within low-redshift galaxy halos trace the dynamical state of circumgalactic gas, using cosmological hydrodynamic simulations that include a well-vetted heuristic model for galactic outflows. We categorize inflowing, outflowing, and ambient gas based on its history and fate as tracked in our simulation. Following our earlier work showing that the ionisation level of absorbers was a primary factor in determining the physical conditions of absorbing gas, we show here that it is also a governing factor for its dynamical state. Low-ionisation metal absorbers (e.g. MgII) tend to arise in gas that will fall onto galaxies within several Gyr, while high-ionisation metal absorbers (e.g. OVI) generally trace material that was deposited by outflows many Gyr ago. Inflowing gas is dominated by enriched material that was previously ejected in an outflow, hence accretion at low redshifts is typically substantially enriched. Recycling wind material is preferentially found closer to galaxies, and is more dominant in lower-mass halos since high-mass halos have more hot gas that is able to support itself against infall. Low-mass halos also tend to re-eject more of their accreted material, owing to our outflow prescription that employs higher mass loading factors for lower-mass galaxies. Typical HI absorbers trace unenriched ambient material that is not participating in the baryon cycle, but stronger HI absorbers arise in cool, enriched inflowing gas. Instantaneous radial velocity measures of absorbers are generally poor at distinguishing between inflowing and outflowing gas, except in the case of very recent outflows. These results suggest that probing halo gas using a range of absorbers can provide detailed information about the amount and physical conditions of material that is participating in the baryon cycle. △ Less

Submitted 11 October, 2013; v1 submitted 23 September, 2013; originally announced September 2013.

Comments: 25 pages, 13 Figures, submitted to MNRAS. Updated citations from original submission

Novel Probes of Gravity and Dark Energy

Authors: Bhuvnesh Jain, Austin Joyce, Rodger Thompson, Amol Upadhye, James Battat, Philippe Brax, Anne-Christine Davis, Claudia de Rham, Scott Dodelson, Adrienne Erickcek, Gregory Gabadadze, Wayne Hu, Lam Hui, Dragan Huterer, Marc Kamionkowski, Justin Khoury, Kazuya Koyama, Baojiu Li, Eric Linder, Fabian Schmidt, Roman Scoccimarro, Glenn Starkman, Chris Stubbs, Masahiro Takada, Andrew Tolley , et al. (7 additional authors not shown)

Abstract: The discovery of cosmic acceleration has stimulated theorists to consider dark energy or modifications to Einstein's General Relativity as possible explanations. The last decade has seen advances in theories that go beyond smooth dark energy -- modified gravity and interactions of dark energy. While the theoretical terrain is being actively explored, the generic presence of fifth forces and dark s… ▽ More The discovery of cosmic acceleration has stimulated theorists to consider dark energy or modifications to Einstein's General Relativity as possible explanations. The last decade has seen advances in theories that go beyond smooth dark energy -- modified gravity and interactions of dark energy. While the theoretical terrain is being actively explored, the generic presence of fifth forces and dark sector couplings suggests a set of distinct observational signatures. This report focuses on observations that differ from the conventional probes that map the expansion history or large-scale structure. Examples of such novel probes are: detection of scalar fields via lab experiments, tests of modified gravity using stars and galaxies in the nearby universe, comparison of lensing and dynamical masses of galaxies and clusters, and the measurements of fundamental constants at high redshift. The observational expertise involved is very broad as it spans laboratory experiments, high resolution astronomical imaging and spectroscopy and radio observations. In the coming decade, searches for these effects have the potential for discovering fundamental new physics. We discuss how the searches can be carried out using experiments that are already under way or with modest adaptations of existing telescopes or planned experiments. The accompanying paper on the Growth of Cosmic Structure describes complementary tests of gravity with observations of large-scale structure. △ Less

Submitted 24 September, 2013; v1 submitted 20 September, 2013; originally announced September 2013.

Comments: Report from the "Dark Energy and CMB" working group for the American Physical Society's Division of Particles and Fields long-term planning exercise ("Snowmass")

The AGORA High-Resolution Galaxy Simulations Comparison Project

Authors: Ji-hoon Kim, Tom Abel, Oscar Agertz, Greg L. Bryan, Daniel Ceverino, Charlotte Christensen, Charlie Conroy, Avishai Dekel, Nickolay Y. Gnedin, Nathan J. Goldbaum, Javiera Guedes, Oliver Hahn, Alexander Hobbs, Philip F. Hopkins, Cameron B. Hummels, Francesca Iannuzzi, Dusan Keres, Anatoly Klypin, Andrey V. Kravtsov, Mark R. Krumholz, Michael Kuhlen, Samuel N. Leitner, Piero Madau, Lucio Mayer, Christopher E. Moody , et al. (21 additional authors not shown)

Abstract: We introduce the AGORA project, a comprehensive numerical study of well-resolved galaxies within the LCDM cosmology. Cosmological hydrodynamic simulations with force resolutions of ~100 proper pc or better will be run with a variety of code platforms to follow the hierarchical growth, star formation history, morphological transformation, and the cycle of baryons in and out of 8 galaxies with halo… ▽ More We introduce the AGORA project, a comprehensive numerical study of well-resolved galaxies within the LCDM cosmology. Cosmological hydrodynamic simulations with force resolutions of ~100 proper pc or better will be run with a variety of code platforms to follow the hierarchical growth, star formation history, morphological transformation, and the cycle of baryons in and out of 8 galaxies with halo masses M_vir ~= 1e10, 1e11, 1e12, and 1e13 Msun at z=0 and two different ("violent" and "quiescent") assembly histories. The numerical techniques and implementations used in this project include the smoothed particle hydrodynamics codes GADGET and GASOLINE, and the adaptive mesh refinement codes ART, ENZO, and RAMSES. The codes will share common initial conditions and common astrophysics packages including UV background, metal-dependent radiative cooling, metal and energy yields of supernovae, and stellar initial mass function. These are described in detail in the present paper. Subgrid star formation and feedback prescriptions will be tuned to provide a realistic interstellar and circumgalactic medium using a non-cosmological disk galaxy simulation. Cosmological runs will be systematically compared with each other using a common analysis toolkit, and validated against observations to verify that the solutions are robust - i.e., that the astrophysical assumptions are responsible for any success, rather than artifacts of particular implementations. The goals of the AGORA project are, broadly speaking, to raise the realism and predictive power of galaxy simulations and the understanding of the feedback processes that regulate galaxy "metabolism." The proof-of-concept dark matter-only test of the formation of a galactic halo with a z=0 mass of M_vir ~= 1.7e11 Msun by 9 different versions of the participating codes is also presented to validate the infrastructure of the project. △ Less

Submitted 24 December, 2013; v1 submitted 12 August, 2013; originally announced August 2013.

Comments: 21 pages, 6 figures, Accepted for publication in the Astrophysical Journal Supplement, Image resolution greatly reduced, High-resolution version of this article and more information about the AGORA Project including the science goals of the 13 Working Groups are available at http://www.AGORAsimulations.org/ and at http://sites.google.com/site/santacruzcomparisonproject/details/

Journal ref: Astrophysical Journal Supplement Series, 210, 14, 2014

The Magellan Uniform Survey of Damped Lyman alpha Systems I: Cosmic Metallicity Evolution

Authors: Regina A. Jorgenson, Michael T. Murphy, Rodger Thompson

Abstract: We present the chemical abundance measurements of the first large, medium-resolution, uniformly selected damped Lyman alpha system (DLA) survey. The sample contains 99 DLAs towards 89 quasars selected from the SDSS DR5 DLA sample in a uniform way. We analyze the metallicities and kinematic diagnostics, including the velocity width of 90% of the optical depth, Δv_{90}, and the equivalent widths of… ▽ More We present the chemical abundance measurements of the first large, medium-resolution, uniformly selected damped Lyman alpha system (DLA) survey. The sample contains 99 DLAs towards 89 quasars selected from the SDSS DR5 DLA sample in a uniform way. We analyze the metallicities and kinematic diagnostics, including the velocity width of 90% of the optical depth, Δv_{90}, and the equivalent widths of the SiII 1526, CIV 1548 and MgII 2796 transitions. To avoid strong line-saturation effects on the metallicities measured in medium-resolution spectra (FWHM~71 km/s), we derived metallicities from metal transitions which absorbed at most 35% of the quasar continuum flux. We find the evolution in cosmic mean metallicity of the sample, <Z> = (-0.04 +- 0.013)z -(1.06 +- 0.36), consistent with no evolution over the redshift range z ~ [2.2, 4.4], but note that the majority of our sample falls at z ~ [2.2, 3.5]. The apparent lack of metallicity evolution with redshift is also seen in a lack of evolution in the median Δv_{90} and SiII 1526 equivalent width values. While this result may seem to conflict with other large surveys that have detected significant metallicity evolution, such as Rafelski et al. 2012 who found <Z> = (-0.22 +- 0.03)z - (0.65 +- 0.09) over z ~ [0, 5], several tests show that these surveys are not inconsistent with our new result. However, over the smaller redshift range covered by our uniformly-selected sample, the true evolution of the cosmic mean metallicity in DLAs may be somewhat flatter than the Rafelski et al. 2012 estimate. △ Less

Submitted 16 July, 2013; originally announced July 2013.

Comments: 20 pages, 21 figures, accepted for publication in MNRAS

A new substantive proton to electron mass ratio constraint on rolling scalar field cosmologies

Authors: Rodger I. Thompson

Abstract: New PKS1830-211 radio frequency observations of methanol at a redshift of 0.88582 have established the most stringent limits on changes in the proton to electron mass ratio mu to date. The observations place the limit of (delta mu)/mu </= (0.0 +/- 1.0) x 10^{-7} which is approximately a factor of four lower than the previous lowest limit at a redshift of 0.6742. This stringent limit at a look back… ▽ More New PKS1830-211 radio frequency observations of methanol at a redshift of 0.88582 have established the most stringent limits on changes in the proton to electron mass ratio mu to date. The observations place the limit of (delta mu)/mu </= (0.0 +/- 1.0) x 10^{-7} which is approximately a factor of four lower than the previous lowest limit at a redshift of 0.6742. This stringent limit at a look back time of roughly half the age of the universe has profound implications for rolling scalar field cosmologies and the new physics that they require. Many of these cosmologies invoke a scalar field phi that is also coupled to the electromagnetic field causing the values of the fundamental constants, mu and the fine structure constant alpha to roll with time. If the lowest expected value of the coupling to mu, zeta_{mu}$ is invoked the new limit requires a limit on the dark energy equation of state parameter w such that w+1 </= 0.001 at a redshift of 0.88582. This eliminates almost all of the expected parameter space for such cosmologies and new physics that have a coupling to the electromagnetic field. In these cases the limit requires that w must be extremely close to -1 for the last half of the age of the universe or that the coupling of the rolling scalar field to mu and the electromagnetic field be significantly below or at the limit of its expected range. The new observations solidify the role of fundamental constants in providing probes of the possible cosmologies and new physics to explain the acceleration of the expansion of the universe. △ Less

Submitted 22 February, 2013; originally announced February 2013.

Comments: 5 pages, 5 figures, accepted for publication in MNRAS