Source author record

Mark L. A. Richardson

Mark L. A. Richardson appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

10works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

10 published item(s)

preprint2022arXiv

Towards convergence of turbulent dynamo amplification in cosmological simulations of galaxies

Our understanding of the process through which magnetic fields reached their observed strengths in present-day galaxies remains incomplete. One of the advocated solutions is a turbulent dynamo mechanism that rapidly amplifies weak magnetic field seeds to the order of ${\sim}μ$G. However, simulating the turbulent dynamo is a very challenging computational task due to the demanding span of spatial scales and the complexity of the required numerical methods. In particular, turbulent velocity and magnetic fields are extremely sensitive to the spatial discretisation of simulated domains. To explore how refinement schemes affect galactic turbulence and amplification of magnetic fields in cosmological simulations, we compare two refinement strategies. A traditional quasi-Lagrangian adaptive mesh refinement approach focusing spatial resolution on dense regions, and a new refinement method that resolves the entire galaxy with a high resolution quasi-uniform grid. Our new refinement strategy yields much faster magnetic energy amplification than the quasi-Lagrangian method, which is also significantly greater than the adiabatic compressional estimate indicating that the extra amplification is produced through stretching of magnetic field lines. Furthermore, with our new refinement the magnetic energy growth factor scales with resolution following $\propto \Dres^{-1/2}$, in much better agreement with small-scale turbulent box simulations. Finally, we find evidence suggesting most magnetic amplification in our simulated galaxies occurs in the warm phase of their interstellar medium, which has a better developed turbulent field with our new refinement strategy.

preprint2020arXiv

Simulating gas kinematic studies of high-redshift galaxies with the HARMONI Integral Field Spectrograph

We present simulated observations of gas kinematics in galaxies formed in 10 pc resolution cosmological simulations with the hydrodynamical + N-body code RAMSES, using the new RAMSES2HSIM pipeline with the simulated observing pipeline (HSIM) for the ELT HARMONI IFU spectrograph. We post-process the galaxy's gas kinematics and Hα line emission for each simulation cell, and integrate the emission to produce an extinction-corrected input cube. We then simulate observations of the input cube with HARMONI, for a range of exposure times, spatial sampling, and spectral resolution. We analyze the mock observations to recover galaxy properties such as its kinematics and compare with the known simulation values. We investigate the cause of biases between the 'real' and 'observed' kinematic values, demonstrating the sensitivity of the inferred rotation curve to knowledge of the instrument's point spread function.

preprint2016arXiv

Comparing Simulations of AGN Feedback

We perform adaptive mesh refinement (AMR) and smoothed particle hydrodynamics (SPH) cosmological zoom simulations of a region around a forming galaxy cluster, comparing the ability of the methods to handle successively more complex baryonic physics. In the simplest, non-radiative case, the two methods are in good agreement with each other, but the SPH simulations generate central cores with slightly lower entropies and virial shocks at slightly larger radii, consistent with what has been seen in previous studies. The inclusion of radiative cooling, star formation, and stellar feedback leads to much larger differences between the two methods. Most dramatically, at z=5, rapid cooling in the AMR case moves the accretion shock well within the virial radius, while this shock remains near the virial radius in the SPH case, due to excess heating, coupled with poorer capturing of the shock width. On the other hand, the addition of feedback from active galactic nuclei (AGN) to the simulations results in much better agreement between the methods. In this case both simulations display halo gas entropies of 100 keV cm^2, similar decrements in the star-formation rate, and a drop in the halo baryon content of roughly 30%. This is consistent with AGN growth being self-regulated, regardless of the numerical method. However, the simulations with AGN feedback continue to differ in aspects that are not self-regulated, such that in SPH a larger volume of gas is impacted by feedback, and the cluster still has a lower entropy central core.

preprint2014arXiv

Galactic winds and stellar populations in Lyman $α$ emitting galaxies at z ~ 3.1

We present a sample of 33 spectroscopically confirmed z ~ 3.1 Ly$α$-emitting galaxies (LAEs) in the Cosmological Evolution Survey (COSMOS) field. This paper details the narrow-band survey we conducted to detect the LAE sample, the optical spectroscopy we performed to confirm the nature of these LAEs, and a new near-infrared spectroscopic detection of the [O III] 5007 Å line in one of these LAEs. This detection is in addition to two [O III] detections in two z ~ 3.1 LAEs we have reported on previously (McLinden et al 2011). The bulk of the paper then presents detailed constraints on the physical characteristics of the entire LAE sample from spectral energy distribution (SED) fitting. These characteristics include mass, age, star-formation history, dust content, and metallicity. We also detail an approach to account for nebular emission lines in the SED fitting process - wherein our models predict the strength of the [O III] line in an LAE spectrum. We are able to study the success of this prediction because we can compare the model predictions to our actual near-infrared observations both in galaxies that have [O III] detections and those that yielded non-detections. We find a median stellar mass of 6.9 $\times$ 10$^8$ M$_{\odot}$ and a median star formation rate weighted stellar population age of 4.5 $\times$ 10$^6$ yr. In addition to SED fitting, we quantify the velocity offset between the [O III] and Ly$α$ lines in the galaxy with the new [O III] detection, finding that the Ly$α$ line is shifted 52 km s$^{-1}$ redwards of the [O III] line, which defines the systemic velocity of the galaxy.

preprint2013arXiv

[NeIII]/[OII] as an Ionization Parameter Diagnostic in Star-Forming Galaxies

We present our parameterizations of the log([NeIII]3869/[OII]3727) (Ne3O2) and log([OIII]5007/[OII]3727) ratios as comparable and effective diagnostics of ionization parameter in star-forming galaxies. Our calibrations are based on the most recent generations of the Starburst99/Mappings III photoionization models, which extend up to the extremely high values of ionization parameter found in high-redshift galaxies. While similar calibrations have been presented previously for O3O2, this is the first such calibration of Ne3O2. We illustrate the tight correlation between these two ratios for star-forming galaxies and discuss the underlying physics that dictates their very similar evolution. Based on this work, we propose the Ne3O2 ratio as a new and useful diagnostic of ionization parameter for star-forming galaxies. Given the Ne3O2 ratio's relative insensitivity to reddening, this ratio is particularly valuable for use with galaxies that have uncertain amounts of extinction. The short wavelengths of the Ne3O2 ratio can also be applied out to very high redshifts, extending studies of galaxies' ionization parameters out to z ~ 1.6 with optical spectroscopy and z ~ 5.2 with ground-based near-infrared spectra.

preprint2013arXiv

Formation of Compact Clusters from High Resolution Hybrid Cosmological Simulations

The early Universe hosted a large population of small dark matter `minihalos' that were too small to cool and form stars on their own. These existed as static objects around larger galaxies until acted upon by some outside influence. Outflows, which have been observed around a variety of galaxies, can provide this influence in such a way as to collapse, rather than disperse the minihalo gas. Gray & Scannapieco performed an investigation in which idealized spherically-symmetric minihalos were struck by enriched outflows. Here we perform high-resolution cosmological simulations that form realistic minihalos, which we then extract to perform a large suite of simulations of outflow-minihalo interactions including non-equilibrium chemical reactions. In all models, the shocked minihalo forms molecules through non-equilibrium reactions, and then cools to form dense chemically homogenous clumps of star-forming gas. The formation of these high-redshift clusters will be observable with the next generation of telescopes, and the largest of them should survive to the present day, having properties similar to halo globular clusters.

preprint2013arXiv

High-Redshift Galaxies with Large Ionization Parameters

Motivated by recent observations of galaxies dominated by emission lines, which show evidence of being metal poor with young stellar populations, we present calculations of multiple model grids with a range of abundances, ionization parameters, and stellar ages, finding that the predicted spectral line diagnostics are heavily dependent on all three parameters. These new model grids extend the ionization parameter to larger values than typically explored. We compare these model predictions with previous observations of such objects, including two new Lyman-$α$ emitting galaxies (LAE) that we have observed. Our models give improved constraints on the metallicity and ionization parameter of these previously studied objects, as we are now able to consider high ionization parameter models. However, similar to previous work, these models have difficulty predicting large line diagnostics for high ionization potential species, requiring future work refining the modelling of FUV photons. Our model grids are also able to constrain the metallicity and ionization parameter of our LAEs, and give constraints on their Ly$α$ escape fractions, all of which are consistent with recent lower redshift studies of LAEs.

preprint2013arXiv

Hybrid Cosmological Simulations with Stream Velocities

In the early universe, substantial relative "stream" velocities between the gas and dark matter arise due to radiation pressure and persist after recombination. To asses the impact of these velocities on high-redshift structure formation, we carry out a suite of high-resolution Adaptive Mesh Refinement (AMR) cosmological simulations, which use Smoothed Particle Hydrodynamic datasets as initial conditions, converted using a new tool developed for this work. These simulations resolve structures with masses as small as a few 100 M$_\odot$, and we focus on the $10^6$ M$_\odot$ "mini-halos" in which the first stars formed. At $z \approx 17,$ the presence of stream velocities has only a minor effect on the number density of halos below $10^6$ M$_\odot$, but it greatly suppresses gas accretion onto all halos and the dark matter structures around them. Stream velocities lead to significantly lower halo gas fractions, especially for $\approx 10^5$ M$_\odot$ objects, an effect that is likely to depend on the orientation of a halo's accretion lanes. This reduction in gas density leads to colder, more compact radial profiles, and it substantially delays the redshift of collapse of the largest halos, leading to delayed star formation and possibly delayed reionization. These many differences suggest that future simulations of early cosmological structure formation should include stream velocities to properly predict gas evolution, star-formation, and the epoch of reionization.

preprint2013arXiv

The Dynamical Masses, Densities, and Star Formation Scaling Relations of Lyman Alpha Galaxies

We present the first dynamical mass measurements for Lyman alpha galaxies at high redshift, based on velocity dispersion measurements from rest-frame optical emission lines and size measurements from HST imaging, for a sample of nine galaxies drawn from four surveys. These measurements enable us to study the nature of Lyman alpha galaxies in the context of galaxy scaling relations. The resulting dynamical masses range from 1e9 to 1e10 solar masses. We also fit stellar population models to our sample, and use them to plot the Lyman alpha sample on a stellar mass vs. line width relation. Overall, the Lyman alpha galaxies follow well the scaling relation established by observing star forming galaxies at lower redshift (and without regard for Lyman alpha emission), though in 1/3 of the Lyman alpha galaxies, lower-mass fits are also acceptable. In all cases, the dynamical masses agree with established stellarmass-linewidth relation. Using the dynamical masses as an upper limit on gas mass, we show that Lyman alpha galaxies resemble starbursts (rather than "normal" galaxies) in the relation between gas mass surface density and star formation activity, in spite of relatively modest star formation rates. Finally, we examine the mass densities of these galaxies, and show that their future evolution likely requires dissipational ("wet") merging. In short, we find that Lyman alpha galaxies are low mass cousins of larger starbursts.

preprint2011arXiv

First spectroscopic measurements of [OIII] emission from Lyman-alpha selected field galaxies at z ~ 3.1

We present the first spectroscopic measurements of the [OIII] 5007 A line in two z ~ 3.1 Lyman-alpha emitting galaxies (LAEs) using the new near-infrared instrument LUCIFER1 on the 8.4m Large Binocular Telescope (LBT). We also describe the optical imaging and spectroscopic observations used to identify these Lya emitting galaxies. Using the [OIII] line we have measured accurate systemic redshifts for these two galaxies, and discovered a velocity offset between the [OIII] and Ly-alpha lines in both, with the Lya line peaking 342 and 125 km/s redward of the systemic velocity. These velocity offsets imply that there are powerful outflows in high-redshift LAEs. They also ease the transmission of Lya photons through the interstellar medium and intergalactic medium around the galaxies. By measuring these offsets directly, we can refine both Lya-based tests for reionization, and Lya luminosity function measurements where the Lya forest affects the blue wing of the line. Our work also provides the first direct constraints on the strength of the [OIII] line in high-redshift LAEs. We find [OIII] fluxes of 7 and 36 x 10^-17 erg s^-1 cm^-2 in two z ~ 3.1 LAEs. These lines are strong enough to dominate broad-band flux measurements that include the line (in thiscase, K_s band photometry). Spectral energy distribution fits that do not account for the lines would therefore overestimate the 4000 A (and/or Balmer) break strength in such galaxies, and hence also the ages and stellar masses of such high-z galaxies.