-
On the limitations of H alpha luminosity as a star formation tracer in spatially resolved observations
Authors:
Zipeng Hu,
Benjamin D. Wibking,
Mark R. Krumholz,
Christoph Federrath
Abstract:
This study examines the limitations of H$α$ luminosity as a tracer of star formation rates (SFR) in spatially resolved observations. We carry out high-resolution simulations of a Milky Way-like galaxy including both supernova and photoionization feedback, and from these we generate synthetic H$α$ emission maps that we compare to maps of the true distribution of young stellar objects (YSOs) on scal…
▽ More
This study examines the limitations of H$α$ luminosity as a tracer of star formation rates (SFR) in spatially resolved observations. We carry out high-resolution simulations of a Milky Way-like galaxy including both supernova and photoionization feedback, and from these we generate synthetic H$α$ emission maps that we compare to maps of the true distribution of young stellar objects (YSOs) on scales from whole-galaxy to individual molecular clouds ($\lesssim 100$ pc). Our results reveal significant spatial mismatches between H$α$ and true YSO maps on sub-100 pc scales, primarily due to ionizing photon leakage, with a secondary contribution from young stars drifting away from their parent molecular clouds. On small scales these effect contribute significantly to the observed anti-correlation between gas and star formation, such that there is noticeably less anti-correlation if we replace an H$α$-based star formation map with a YSO-based one; this in turn implies that previous studies have underestimated the time it takes for young stars to disperse their parent molecular clouds. However, these effects are limited in dense regions with hydrogen columns $N_\mathrm{H} > 3 \times 10^{21}$ cm$^{-2}$, where the H$α$- and YSO-based SFR maps show better agreement. Based on this finding we propose a calibration model that can precisely measure the SFR of large molecular clouds (mean radius > 100 pc) with a combination of H$α$ and CO observations, which provides a foundation for future study of star formation processes in extragalactic molecular clouds.
△ Less
Submitted 16 June, 2024;
originally announced June 2024.
-
An Asymptotically-Correct Implicit-Explicit Time Integration Scheme for Finite Volume Radiation-Hydrodynamics
Authors:
Chong-Chong He,
Benjamin D. Wibking,
Mark R. Krumholz
Abstract:
Numerical radiation-hydrodynamics (RHD) for non-relativistic flows is a challenging problem because it encompasses processes acting over a very broad range of timescales, and where the relative importance of these processes often varies by orders of magnitude across the computational domain. Here we present a new implicit-explicit (IMEX) method for numerical RHD that has a number of desirable prop…
▽ More
Numerical radiation-hydrodynamics (RHD) for non-relativistic flows is a challenging problem because it encompasses processes acting over a very broad range of timescales, and where the relative importance of these processes often varies by orders of magnitude across the computational domain. Here we present a new implicit-explicit (IMEX) method for numerical RHD that has a number of desirable properties that have not previously been combined in a single method. Our scheme is based on moments and allows machine-precision conservation of energy and momentum, making it highly suitable for adaptive mesh refinement applications; it requires no more communication than hydrodynamics and includes no non-local iterative steps, making it highly suitable for massively parallel and GPU-based systems where communication is a bottleneck; and we show that it is asymptotically-accurate in the streaming, static diffusion, and dynamic diffusion limits, including in the so-called asymptotic diffusion regime where the computational grid does not resolve the photon mean free path. We implement our method in the GPU-accelerated RHD code QUOKKA and show that it passes a wide range of numerical tests.
△ Less
Submitted 12 April, 2024;
originally announced April 2024.
-
Quokka-based Understanding of Outflows (QED) I. Metal loading, phase structure, and convergence testing for Solar neighbourhood conditions
Authors:
Aditi Vijayan,
Mark R. Krumholz,
Benjamin D. Wibking
Abstract:
Multiphase galactic outflows, generated by supernova feedback, are likely to be more metal-rich than the interstellar media from which they are driven due to incomplete mixing between supernova ejecta and the ambient ISM. This enrichment is important for shaping galactic metallicities and metallicity gradients, but measuring it quantitatively from simulations requires resolution high enough to res…
▽ More
Multiphase galactic outflows, generated by supernova feedback, are likely to be more metal-rich than the interstellar media from which they are driven due to incomplete mixing between supernova ejecta and the ambient ISM. This enrichment is important for shaping galactic metallicities and metallicity gradients, but measuring it quantitatively from simulations requires resolution high enough to resolve mass, momentum and energy exchanges between the different phases of the outflows. In this context, we present simulations of outflows, driven by SN feedback, conducted using \textsc{Quokka}, a new GPU-optimised AMR radiation-hydrodynamics code. This code allows us to reach combinations of resolution, simulation volume, and simulation duration larger than those that have previously been possible, and to resolve all gas phases from cold neutral medium, $T \sim 100$ K, to hot ionised gas, $T \gtrsim 10^7$ K. In this, a first of a series of papers exploring generation and evolution of multiphase outflows from a wide range of galactic environments and star formation rates, we quantify the extent of selective metal loading in Solar neighbourhood-like environments. We explain the selective metal loading we find as a result of the transport of metals within and between phases, a phenomenon we can study owing to the parsec-scale resolution that our simulations achieve. We also quantify the sensitivity of metal loading studies to numerical resolution, and present convergence criteria for future studies.
△ Less
Submitted 14 September, 2023;
originally announced September 2023.
-
The sub-critical illusion: synthetic Zeeman effect observations from galactic zoom-in simulations
Authors:
Zipeng Hu,
Benjamin D. Wibking,
Mark R. Krumholz
Abstract:
Mass-to-flux ratios measured via the Zeeman effect suggest the existence of a transition from a magnetically sub-critical state in HI clouds to a super-critical state in molecular clouds. However, due to projection, chemical, and excitation effects, Zeeman measurements are subject to a number of biases, and may not reflect the true relations between gravitational and magnetic energies. In this pap…
▽ More
Mass-to-flux ratios measured via the Zeeman effect suggest the existence of a transition from a magnetically sub-critical state in HI clouds to a super-critical state in molecular clouds. However, due to projection, chemical, and excitation effects, Zeeman measurements are subject to a number of biases, and may not reflect the true relations between gravitational and magnetic energies. In this paper, we carry out simulations of the formation of magnetised molecular clouds, zooming in from an entire galaxy to sub-pc scales, which we post-process to produce synthetic HI and OH Zeeman measurements. The mass-to-flux ratios we recover from the simulated observations show a transition in magnetic criticality that closely matches observations, but we find that the gravitational-magnetic energy ratios on corresponding scales are mostly super-critical, even in the HI regime. We conclude that HI clouds in the process of assembling to form molecular clouds are already super-critical even before H_2 forms, and that the apparent transition from sub- to super-criticality between HI and H_2 is primarily an illusion created by chemical and excitation biases affecting the Zeeman measurements.
△ Less
Submitted 22 February, 2023;
originally announced February 2023.
-
A radiation hydrodynamics scheme on adaptive meshes using the Variable Eddington Tensor (VET) closure
Authors:
Shyam H. Menon,
Christoph Federrath,
Mark R. Krumholz,
Rolf Kuiper,
Benjamin D. Wibking,
Manuel Jung
Abstract:
We present a new algorithm to solve the equations of radiation hydrodynamics (RHD) in a frequency-integrated, two-moment formulation. Novel features of the algorithm include i) the adoption of a non-local Variable Eddington Tensor (VET) closure for the radiation moment equations, computed with a ray-tracing method, ii) support for adaptive mesh refinement (AMR), iii) use of a time-implicit Godunov…
▽ More
We present a new algorithm to solve the equations of radiation hydrodynamics (RHD) in a frequency-integrated, two-moment formulation. Novel features of the algorithm include i) the adoption of a non-local Variable Eddington Tensor (VET) closure for the radiation moment equations, computed with a ray-tracing method, ii) support for adaptive mesh refinement (AMR), iii) use of a time-implicit Godunov method for the hyperbolic transport of radiation, and iv) a fixed-point Picard iteration scheme to accurately handle the stiff nonlinear gas-radiation energy exchange. Tests demonstrate that our scheme works correctly, yields accurate rates of energy and momentum transfer between gas and radiation, and obtains the correct radiation field distribution even in situations where more commonly used -- but less accurate -- closure relations like the Flux-limited Diffusion and Moment-1 approximations fail. Our scheme presents an important step towards performing RHD simulations with increasing spatial and directional accuracy, effectively improving their predictive capabilities.
△ Less
Submitted 22 March, 2022;
originally announced March 2022.
-
VETTAM: A scheme for radiation hydrodynamics with adaptive mesh refinement using the variable Eddington tensor method
Authors:
Shyam H. Menon,
Christoph Federrath,
Mark R. Krumholz,
Rolf Kuiper,
Benjamin D. Wibking,
Manuel Jung
Abstract:
We present Variable Eddington Tensor-closed Transport on Adaptive Meshes (\texttt{VETTAM}), a new algorithm to solve the equations of radiation hydrodynamics (RHD) with support for adaptive mesh refinement (AMR) in a frequency-integrated, two-moment formulation. The method is based on a non-local Variable Eddington Tensor (VET) closure computed with a hybrid characteristics scheme for ray tracing.…
▽ More
We present Variable Eddington Tensor-closed Transport on Adaptive Meshes (\texttt{VETTAM}), a new algorithm to solve the equations of radiation hydrodynamics (RHD) with support for adaptive mesh refinement (AMR) in a frequency-integrated, two-moment formulation. The method is based on a non-local Variable Eddington Tensor (VET) closure computed with a hybrid characteristics scheme for ray tracing. We use a Godunov method for the hyperbolic transport of radiation with an implicit backwards-Euler temporal update to avoid the explicit timestep constraint imposed by the light-crossing time, and a fixed-point Picard iteration scheme to handle the nonlinear gas-radiation exchange term, with the two implicit update stages jointly iterated to convergence. We also develop a modified wave-speed correction method for AMR, which we find to be crucial for obtaining accurate results in the diffusion regime. We demonstrate the robustness of our scheme with a suite of pure radiation and RHD tests, and show that it successfully captures the streaming, static diffusion, and dynamic diffusion regimes and the spatial transitions between them, casts sharp shadows, and yields accurate results for rates of momentum and energy exchange between radiation and gas. A comparison between different closures for the radiation moment equations, with the Eddington approximation (0th-moment closure) and the $M_1$ approximation (1st-moment closure), demonstrates the advantages of the VET method (2nd-moment closure) over the simpler closure schemes. \texttt{VETTAM} has been coupled to the AMR \texttt{FLASH} (magneto-)hydrodynamics code and we summarize by reporting performance features and bottlenecks of our implementation.
△ Less
Submitted 17 February, 2022;
originally announced February 2022.
-
Quokka: A code for two-moment AMR radiation hydrodynamics on GPUs
Authors:
Benjamin D. Wibking,
Mark R. Krumholz
Abstract:
We present Quokka, a new subcycling-in-time, block-structured adaptive mesh refinement (AMR) radiation hydrodynamics code optimised for graphics processing units (GPUs). Quokka solves the equations of hydrodynamics with the piecewise parabolic method (PPM) in a method-of-lines formulation, and handles radiative transfer via the variable Eddington tensor (VET) radiation moment equations with a loca…
▽ More
We present Quokka, a new subcycling-in-time, block-structured adaptive mesh refinement (AMR) radiation hydrodynamics code optimised for graphics processing units (GPUs). Quokka solves the equations of hydrodynamics with the piecewise parabolic method (PPM) in a method-of-lines formulation, and handles radiative transfer via the variable Eddington tensor (VET) radiation moment equations with a local closure. We use the AMReX library to handle the adaptive mesh management. In order to maximise GPU performance, we combine explicit-in-time evolution of the radiation moment equations with the reduced speed-of-light approximation. We show results for a wide range of test problems for hydrodynamics, radiation, and coupled radiation hydrodynamics. On uniform grids in 3D on a single GPU, our code achieves > 250 million hydrodynamic updates per second and almost 40 million radiation hydrodynamic updates per second. For radiation hydrodynamics problems on uniform grids in 3D, our code scales from 4 GPUs to 256 GPUs with an efficiency of 76 per cent. The code is publicly released under an open-source license on GitHub.
△ Less
Submitted 30 January, 2022; v1 submitted 4 October, 2021;
originally announced October 2021.
-
Exploiting Non-linear Scales in Galaxy-Galaxy Lensing and Galaxy Clustering: A Forecast for the Dark Energy Survey
Authors:
Andrés N. Salcedo,
David H. Weinberg,
Hao-Yi Wu,
Benjamin D. Wibking
Abstract:
The combination of galaxy-galaxy lensing (GGL) and galaxy clustering is a powerful probe of low redshift matter clustering, especially if it is extended to the non-linear regime. To this end, we extend the N-body and halo occupation distribution (HOD) emulator method of arxiv:1907.06293 to model the redMaGiC sample of colour-selected passive galaxies in the Dark Energy Survey (DES), adding paramet…
▽ More
The combination of galaxy-galaxy lensing (GGL) and galaxy clustering is a powerful probe of low redshift matter clustering, especially if it is extended to the non-linear regime. To this end, we extend the N-body and halo occupation distribution (HOD) emulator method of arxiv:1907.06293 to model the redMaGiC sample of colour-selected passive galaxies in the Dark Energy Survey (DES), adding parameters that describe central galaxy incompleteness, galaxy assembly bias, and a scale-independent multiplicative lensing bias $A_{lens}$. We use this emulator to forecast cosmological constraints attainable from the GGL surface density profile $ΔΣ(r_p)$ and the projected galaxy correlation function $w_{p,gg}(r_p)$ in the final (Year 6) DES data set over scales $r_p=0.3-30h^{-1}$ Mpc. For a $3\%$ prior on $A_{lens}$ we forecast precisions of $1.9\%$, $2.0\%$, and $1.9\%$ on $Ω_m$, $σ_8$, and $S_8 \equiv σ_8Ω_m^{0.5}$, marginalized over all halo occupation distribution (HOD) parameters as well as $A_{lens}$ and a point-mass contribution to $ΔΣ$. Adding scales $r_p=0.3-3h^{-1}$ Mpc improves the $S_8$ precision by a factor of $\sim1.6$ relative to a large scale ($3.0-30.0h^{-1}$ Mpc) analysis, equivalent to increasing the survey area by a factor of ${\sim}2.6$. Sharpening the $A_{lens}$ prior to $1\%$ further improves the $S_8$ precision by a factor of $1.7$ (to $1.1\%$), and it amplifies the gain from including non-linear scales. Our emulator achieves percent-level accuracy similar to the projected DES statistical uncertainties, demonstrating the feasibility of a fully non-linear analysis. Obtaining precise parameter constraints from multiple galaxy types and from measurements that span linear and non-linear clustering offers many opportunities for internal cross-checks, which can diagnose systematics and demonstrate the robustness of cosmological results.
△ Less
Submitted 13 July, 2021;
originally announced July 2021.
-
The global structure of magnetic fields and gas in simulated Milky Way-analogue galaxies
Authors:
Benjamin D. Wibking,
Mark R. Krumholz
Abstract:
We simulate an isolated, magnetised Milky Way-like disc galaxy using a self-consistent model of unresolved star formation and feedback, evolving the system until it reaches statistical steady state. We show that the quasi-steady-state structure is distinctly layered in galactocentric height $z$, with an innermost region having comparable gas and magnetic pressures (plasma beta $β\sim 1$), an outer…
▽ More
We simulate an isolated, magnetised Milky Way-like disc galaxy using a self-consistent model of unresolved star formation and feedback, evolving the system until it reaches statistical steady state. We show that the quasi-steady-state structure is distinctly layered in galactocentric height $z$, with an innermost region having comparable gas and magnetic pressures (plasma beta $β\sim 1$), an outermost region having dominant gas pressures ($β\gg 1$), and an intermediate region between $300$ pc $\lesssim |z| \lesssim 3$ kpc that is dynamically dominated by magnetic fields ($β\ll 1$). We find field strengths, gas surface densities, and star formation rates that agree well with those observed both in the Galactic centre and in the Solar neighbourhood. The most significant dynamical effect of magnetic fields on the global properties of the disc is a reduction of the star formation rate by a factor of 1.5-2 with respect to an unmagnetised control simulation. At fixed star formation rate, there is no significant difference in the mass outflow rates or profiles between the magnetised and non-magnetised simulations. Our results for the global structure of the magnetic field have significant implications for models of cosmic ray-driven winds and cosmic-ray propagation in the Galaxy, and can be tested against observations with the forthcoming Square Kilometre Array and other facilities. Finally, we report the discovery of a physical error in the implementation of neutral gas heating and cooling in the popular GIZMO code, which may lead to qualitatively incorrect phase structures if not corrected.
△ Less
Submitted 10 May, 2021;
originally announced May 2021.
-
Cosmology with galaxy cluster weak lensing: statistical limits and experimental design
Authors:
Hao-Yi Wu,
David H. Weinberg,
Andrés N. Salcedo,
Benjamin D. Wibking
Abstract:
We forecast constraints on the amplitude of matter clustering sigma_8(z) achievable with the combination of cluster weak lensing and number counts, in current and next-generation weak lensing surveys. We advocate an approach, analogous to galaxy-galaxy lensing, in which the observables in each redshift bin are the mean number counts and the mean weak lensing profile of clusters above a mass proxy…
▽ More
We forecast constraints on the amplitude of matter clustering sigma_8(z) achievable with the combination of cluster weak lensing and number counts, in current and next-generation weak lensing surveys. We advocate an approach, analogous to galaxy-galaxy lensing, in which the observables in each redshift bin are the mean number counts and the mean weak lensing profile of clusters above a mass proxy threshold. The primary astrophysical nuisance parameter is the logarithmic scatter between the mass proxy and true mass near the threshold. For surveys similar to the Dark Energy Survey (DES), the Roman Space Telescope High Latitude Survey (HLS), and the Rubin Observatory Legacy Survey of Space and Time (LSST), we forecast aggregate precision on sigma_8 of 0.26%, 0.24%, and 0.10%, respectively, if the mass-observable scatter has an external prior better than 0.01. These constraints would be degraded by about 20% for a 0.05 prior on scatter in the case of DES or HLS and for a 0.016 prior for LSST. A one-month observing program with Roman Space Telescope targeting approximately 2500 massive clusters could achieve a 0.5% constraint on sigma_8(z=0.7) on its own, or a ~0.33% constraint in combination with the HLS. Realizing the constraining power of clusters requires accurate knowledge of the mass-observable relation and stringent control of systematics. We provide analytic approximations to our numerical results that allow easy scaling to other survey assumptions or other methods of cluster mass estimation.
△ Less
Submitted 26 March, 2021; v1 submitted 3 December, 2020;
originally announced December 2020.
-
Covariance matrices for galaxy cluster weak lensing: from virial regime to uncorrelated large-scale structure
Authors:
Hao-Yi Wu,
David H. Weinberg,
Andrés N. Salcedo,
Benjamin D. Wibking,
Ying Zu
Abstract:
Next-generation optical imaging surveys will revolutionise the observations of weak gravitational lensing by galaxy clusters and provide stringent constraints on growth of structure and cosmic acceleration. In these experiments, accurate modelling of covariance matrices of cluster weak lensing plays the key role in obtaining robust measurements of the mean mass of clusters and cosmological paramet…
▽ More
Next-generation optical imaging surveys will revolutionise the observations of weak gravitational lensing by galaxy clusters and provide stringent constraints on growth of structure and cosmic acceleration. In these experiments, accurate modelling of covariance matrices of cluster weak lensing plays the key role in obtaining robust measurements of the mean mass of clusters and cosmological parameters. We use a combination of analytical calculations and high-resolution N-body simulations to derive accurate covariance matrices that span from the virial regime to linear scales of the cluster-matter cross-correlation. We validate this calculation using a public ray-tracing lensing simulation and provide a software package for calculating covariance matrices for a wide range of cluster and source sample choices. We discuss the relative importance of shape noise and density fluctuations, the impact of radial bin size, and the impact of off-diagonal elements. For a weak lensing source density 10 per square arcmin, shape noise typically dominates the variance on comoving scales less than 5 Mpc/h. However, for 60 per square arcmin, potentially achievable with future weak lensing experiments, density fluctuations typically dominate the variance at scales greater than 1 Mpc/h and remain comparable to shape noise on smaller scales.
△ Less
Submitted 15 October, 2019; v1 submitted 15 July, 2019;
originally announced July 2019.
-
Cosmology with galaxy-galaxy lensing on non-perturbative scales: Emulation method and application to BOSS LOWZ
Authors:
Benjamin D. Wibking,
David H. Weinberg,
Andrés N. Salcedo,
Hao-Yi Wu,
Sukhdeep Singh,
Sergio Rodríguez-Torres,
Lehman H. Garrison,
Daniel J. Eisenstein
Abstract:
We describe our nonlinear emulation (i.e., interpolation) framework that combines the halo occupation distribution (HOD) galaxy bias model with $N$-body simulations of nonlinear structure formation, designed to accurately predict the projected clustering and galaxy-galaxy lensing signals from luminous red galaxies (LRGs) in the redshift range $0.16 < z < 0.36$ on comoving scales $0.6 < r_p < 30$ \…
▽ More
We describe our nonlinear emulation (i.e., interpolation) framework that combines the halo occupation distribution (HOD) galaxy bias model with $N$-body simulations of nonlinear structure formation, designed to accurately predict the projected clustering and galaxy-galaxy lensing signals from luminous red galaxies (LRGs) in the redshift range $0.16 < z < 0.36$ on comoving scales $0.6 < r_p < 30$ \hMpc. The interpolation accuracy is $\lesssim 1-2$ per cent across the entire physically plausible range of parameters for all scales considered. We correctly recover the true value of the cosmological parameter $S_8 = ({σ_8}/{0.8228}) ({Ω_{\text{m}}}/{0.3107})^{0.6}$ from mock measurements produced via subhalo abundance matching (SHAM)-based lightcones designed to approximately match the properties of the SDSS LOWZ galaxy sample. Applying our model to Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 14 (DR14) LOWZ galaxy clustering and galaxy-shear cross-correlation measurements made with Sloan Digital Sky Survey (SDSS) Data Release 8 (DR8) imaging, we perform a prototype cosmological analysis marginalizing over $w$CDM cosmological parameters and galaxy HOD parameters. We obtain a 4.4 per cent measurement of $S_8 = 0.847 \pm 0.037$, in $3.5σ$ tension with the Planck cosmological results of $1.00 \pm 0.02$. We discuss the possibility of underestimated systematic uncertainties or astrophysical effects that could explain this discrepancy.
△ Less
Submitted 3 December, 2019; v1 submitted 14 July, 2019;
originally announced July 2019.
-
Cosmology with Stacked Cluster Weak Lensing and Cluster-Galaxy Cross-Correlations
Authors:
Andrés N. Salcedo,
Benjamin D. Wibking,
David H. Weinberg,
Hao-Yi Wu,
Douglas Ferrer,
Daniel Eisenstein,
Philip Pinto
Abstract:
Cluster weak lensing is a sensitive probe of cosmology, particularly the amplitude of matter clustering $σ_8$ and matter density parameter $Ω_m$. The main nuisance parameter in a cluster weak lensing cosmological analysis is the scatter between the true halo mass and the relevant cluster observable, denoted $σ_{\ln Mc}$. We show that combining the cluster weak lensing observable $ΔΣ$ with the proj…
▽ More
Cluster weak lensing is a sensitive probe of cosmology, particularly the amplitude of matter clustering $σ_8$ and matter density parameter $Ω_m$. The main nuisance parameter in a cluster weak lensing cosmological analysis is the scatter between the true halo mass and the relevant cluster observable, denoted $σ_{\ln Mc}$. We show that combining the cluster weak lensing observable $ΔΣ$ with the projected cluster-galaxy cross-correlation function $w_{p,cg}$ and galaxy auto-correlation function $w_{p,gg}$ can break the degeneracy between $σ_8$ and $σ_{\ln Mc}$ to achieve tight, percent-level constraints on $σ_8$. Using a grid of cosmological N-body simulations, we compute derivatives of $ΔΣ$, $w_{p,cg}$, and $w_{p,gg}$ with respect to $σ_8$, $Ω_m$, $σ_{\ln Mc}$ and halo occupation distribution (HOD) parameters describing the galaxy population. We also compute covariance matrices motivated by the properties of the Dark Energy Suvery (DES) cluster and weak lensing survey and the BOSS CMASS galaxy redshift survey. For our fiducial scenario combining $ΔΣ$, $w_{p,cg}$, and $w_{p,gg}$ measured over $0.3-30.0 \; h^{-1} \; \mathrm{Mpc}$, for clusters at $z=0.35-0.55$ above a mass threshold $M_c\approx 2\times 10^{14} \; h^{-1} \; \mathrm{M_{\odot}}$, we forecast a $1.4\%$ constraint on $σ_8$ while marginalizing over $σ_{\ln Mc}$ and all HOD parameters. Reducing the mass threshold to $1\times 10^{14} \; h^{-1} \; \mathrm{M_{\odot}}$ and adding a $z=0.15-0.35$ redshift bin sharpens this constraint to $0.8\%$. The small scale $(r_p < 3.0 \; h^{-1} \; \mathrm{Mpc})$ ``mass function'' and large scale $(r_p > 3.0 \; h^{-1} \; \mathrm{Mpc})$ ``halo-mass cross-correlation'' regimes of $ΔΣ$ have comparable constraining power, allowing internal consistency tests from such an analysis.
△ Less
Submitted 15 June, 2019;
originally announced June 2019.
-
Radiation pressure in galactic disks: stability, turbulence, and winds in the single-scattering limit
Authors:
Benjamin D. Wibking,
Todd A. Thompson,
Mark R. Krumholz
Abstract:
The radiation force on dust grains may be dynamically important in driving turbulence and outflows in rapidly star-forming galaxies. Recent studies focus on the highly optically-thick limit relevant to the densest ultra-luminous galaxies and super star clusters, where reprocessed infrared photons provide the dominant source of electromagnetic momentum. However, even among starburst galaxies, the g…
▽ More
The radiation force on dust grains may be dynamically important in driving turbulence and outflows in rapidly star-forming galaxies. Recent studies focus on the highly optically-thick limit relevant to the densest ultra-luminous galaxies and super star clusters, where reprocessed infrared photons provide the dominant source of electromagnetic momentum. However, even among starburst galaxies, the great majority instead lie in the so-called "single-scattering" limit, where the system is optically-thick to the incident starlight, but optically-thin to the re-radiated infrared. In this paper we present a stability analysis and multidimensional radiation-hydrodynamic simulations exploring the stability and dynamics of isothermal dusty gas columns in this regime. We describe our algorithm for full angle-dependent radiation transport based on the discontinuous Galerkin finite element method. For a range of near-Eddington fluxes, we show that the medium is unstable, producing convective-like motions in a turbulent atmosphere with a scale height significantly inflated compared to the gas pressure scale height and mass-weighted turbulent energy densities of $\sim 0.01-0.1$ of the midplane radiation energy density, corresponding to mass-weighted velocity dispersions of Mach number $\sim 0.5-2$. Extrapolation of our results to optical depths of $10^3$ implies maximum turbulent Mach numbers of $\sim20$. Comparing our results to galaxy-averaged observations, and subject to the approximations of our calculations, we find that radiation pressure does not contribute significantly to the effective supersonic pressure support in star-forming disks, which in general are substantially sub-Eddington. We further examine the time-averaged vertical density profiles in dynamical equilibrium and comment on implications for radiation-pressure-driven galactic winds.
△ Less
Submitted 4 January, 2018;
originally announced January 2018.
-
Emulating galaxy clustering and galaxy-galaxy lensing into the deeply nonlinear regime: methodology, information, and forecasts
Authors:
Benjamin D. Wibking,
Andrés N. Salcedo,
David H. Weinberg,
Lehman H. Garrison,
Douglas Ferrer,
Jeremy Tinker,
Daniel Eisenstein,
Marc Metchnik,
Philip Pinto
Abstract:
The combination of galaxy-galaxy lensing (GGL) with galaxy clustering is one of the most promising routes to determining the amplitude of matter clustering at low redshifts. We show that extending clustering+GGL analyses from the linear regime down to $\sim 0.5 \, h^{-1}$ Mpc scales increases their constraining power considerably, even after marginalizing over a flexible model of non-linear galaxy…
▽ More
The combination of galaxy-galaxy lensing (GGL) with galaxy clustering is one of the most promising routes to determining the amplitude of matter clustering at low redshifts. We show that extending clustering+GGL analyses from the linear regime down to $\sim 0.5 \, h^{-1}$ Mpc scales increases their constraining power considerably, even after marginalizing over a flexible model of non-linear galaxy bias. Using a grid of cosmological N-body simulations, we construct a Taylor-expansion emulator that predicts the galaxy autocorrelation $ξ_{\text{gg}}(r)$ and galaxy-matter cross-correlation $ξ_{\text{gm}}(r)$ as a function of $σ_8$, $Ω_m$, and halo occupation distribution (HOD) parameters, which are allowed to vary with large scale environment to represent possible effects of galaxy assembly bias. We present forecasts for a fiducial case that corresponds to BOSS LOWZ galaxy clustering and SDSS-depth weak lensing (effective source density $\sim 0.3$ arcmin$^{-2}$). Using tangential shear and projected correlation function measurements over $0.5 \leq r_p \leq 30 \, h^{-1}$ Mpc yields a 1.8% constraint on the parameter combination $σ_8Ω_m^{0.58}$, a factor of two better than a constraint that excludes non-linear scales ($r_p > 2 \, h^{-1}$ Mpc, $4 \, h^{-1}$ Mpc for $γ_t,w_p$). Much of this improvement comes from the non-linear clustering information, which breaks degeneracies among HOD parameters that would otherwise degrade the inference of matter clustering from GGL. Increasing the effective source density to $3$ arcmin$^{-2}$ sharpens the constraint on $σ_8Ω_m^{0.58}$ by a further factor of two. With robust modeling into the non-linear regime, low-redshift measurements of matter clustering at the 1-percent level with clustering+GGL alone are well within reach of current data sets such as those provided by the Dark Energy Survey.
△ Less
Submitted 20 September, 2017;
originally announced September 2017.
-
Towards Accurate Modelling of Galaxy Clustering on Small Scales: Testing the Standard $Λ\mathrm{CDM}$ + Halo Model
Authors:
Manodeep Sinha,
Andreas A. Berlind,
Cameron K. McBride,
Roman Scoccimarro,
Jennifer A. Piscionere,
Benjamin D. Wibking
Abstract:
Interpreting the small-scale clustering of galaxies with halo models can elucidate the connection between galaxies and dark matter halos. Unfortunately, the modelling is typically not sufficiently accurate for ruling out models statistically. It is thus difficult to use the information encoded in small scales to test cosmological models or probe subtle features of the galaxy-halo connection. In th…
▽ More
Interpreting the small-scale clustering of galaxies with halo models can elucidate the connection between galaxies and dark matter halos. Unfortunately, the modelling is typically not sufficiently accurate for ruling out models statistically. It is thus difficult to use the information encoded in small scales to test cosmological models or probe subtle features of the galaxy-halo connection. In this paper, we attempt to push halo modelling into the "accurate" regime with a fully numerical mock-based methodology and careful treatment of statistical and systematic errors. With our forward-modelling approach, we can incorporate clustering statistics beyond the traditional two-point statistics. We use this modelling methodology to test the standard $Λ\mathrm{CDM}$ + halo model against the clustering of SDSS DR7 galaxies. Specifically, we use the projected correlation function, group multiplicity function and galaxy number density as constraints. We find that while the model fits each statistic separately, it struggles to fit them simultaneously. Adding group statistics leads to a more stringent test of the model and significantly tighter constraints on model parameters. We explore the impact of varying the adopted halo definition and cosmological model and find that changing the cosmology makes a significant difference. The most successful model we tried (Planck cosmology with Mvir halos) matches the clustering of low luminosity galaxies, but exhibits a 2.3$σ$ tension with the clustering of luminous galaxies, thus providing evidence that the "standard" halo model needs to be extended. This work opens the door to adding interesting freedom to the halo model and including additional clustering statistics as constraints.
△ Less
Submitted 6 August, 2018; v1 submitted 16 August, 2017;
originally announced August 2017.
