ABSTRACT
We simulate an isolated, magnetized 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 a broken power-law structure in Alfven Mach number and plasma beta. Magnetic pressure exceeds turbulent and thermal pressures after the gas is depleted to levels below that of the present-day Galaxy, but is subdominant at higher gas fractions and star formation rates. We find field strengths, gas surface densities, and star formation rates that agree well with those observed 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 unmagnetized control simulation. At a fixed star formation rate of approximately $2 \, {\rm M}_{\odot }$ yr−1, there is no significant difference in the mass outflow rates or profiles between the magnetized and non-magnetized 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.
ABSTRACT
Multiphase galactic outflows, generated by supernova (SN) feedback, are likely to be more metal rich than the interstellar media from which they are driven due to incomplete mixing between ...SN ejecta and the ambient interstellar medium. 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 QED, which are simulations of outflows, driven by SN feedback, conducted using Quokka, a new GPU-optimized adaptive mesh refinement 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 ∼ 100 K, to hot ionized gas, T ≳ 107 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.
ABSTRACT
Mass-to-flux ratios measured via the Zeeman effect suggest the existence of a transition from a magnetically sub-critical state in H i clouds to a supercritical 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 magnetized molecular clouds, zooming in from an entire galaxy to sub-pc scales, which we post-process to produce synthetic H i 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 supercritical, even in the H i regime. We conclude that H i clouds in the process of assembling to form molecular clouds are already supercritical even before H2 forms, and that the apparent transition from sub- to supercriticality between H i and H2 is primarily an illusion created by chemical and excitation biases affecting the Zeeman measurements.
ABSTRACT
We present quokka, a new subcycling-in-time, block-structured adaptive mesh refinement (AMR) radiation hydrodynamics (RHD) code optimized for graphics processing units (GPUs). quokka solves ...the equations of HD 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 AM management. In order to maximize 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 HD, radiation, and coupled RHD. 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 RHD problems on uniform grids in 3D, our code scales from 4 to 256 GPUs with an efficiency of 76 per cent. The code is publicly released under an open-source license on GitHub.
Numerical radiation-hydrodynamics (RHD) for non-relativistic flows is a challenging problem because it encompasses processes acting over a very broad range of time-scales, 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 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 Graphics Processing Unit (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.
ABSTRACT
Next-generation optical imaging surveys will revolutionize 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 ns = 10 arcmin−2, shape noise typically dominates the variance on comoving scales $r_{\rm p}\lesssim 5\ h^{-1} \, \rm Mpc$. However, for ns = 60 arcmin−2, potentially achievable with future weak lensing experiments, density fluctuations typically dominate the variance at $r_{\rm p}\gtrsim 1\ h^{-1} \, \rm Mpc$ and remain comparable to shape noise on smaller scales.
ABSTRACT
We present Variable Eddington Tensor (VET)-closed Transport on Adaptive Meshes (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 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 time-step 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 M1 approximation (1st-moment closure), demonstrates the advantages of the VET method (2nd-moment closure) over the simpler closure schemes. VETTAM has been coupled to the AMR FLASH (magneto-)hydrodynamics code and we summarize by reporting performance features and bottlenecks of our implementation.
ABSTRACT
We describe our non-linear emulation (i.e. interpolation) framework that combines the halo occupation distribution (HOD) galaxy bias model with N-body simulations of non-linear structure ...formation, designed to accurately predict the projected clustering and galaxy–galaxy lensing signals from luminous red galaxies in the redshift range 0.16 < z < 0.36 on comoving scales 0.6 < rp < 30 $h^{-1} \, \text{Mpc}$. The interpolation accuracy is ≲ 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 S8 = (σ8/0.8228)(Ωm/0.3107)0.6 from mock measurements produced via subhalo abundance matching (SHAM)-based light-cones 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 wCDM cosmological parameters and galaxy HOD parameters. We obtain a 4.4 per cent measurement of S8 = 0.847 ± 0.037, in 3.5σ tension with the Planck cosmological results of 1.00 ± 0.02. We discuss the possibility of underestimated systematic uncertainties or astrophysical effects that could explain this discrepancy.