Researcher profile

Benjamin D. Wibking

Benjamin D. Wibking contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 19 - UnverifiedVerification L1Unclaimed author
5works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

5 published item(s)

preprint2026arXiv

QED V: Variations in metal loading of galactic winds with element nucleosynthetic origin

Type Ia supernovae, type II supernovae, and asymptotic giant branch (AGB) stars are important sites of stellar nucleosynthesis, but they differ greatly in their rates, their location within a galaxy, and the mean thermal energy and abundance distribution of their ejecta. In earlier papers in this series we have shown that a significant fraction of metals newly synthesized by type II supernovae are promptly lost to galactic winds -- i.e., galactic winds are metal loaded. Here we investigate whether the elements returned by type Ia supernovae and AGB stars are similarly metal loaded, or whether metal loading varies significantly with nucleosynthetic site. We use a series of high-resolution ``tall box'' simulations of the interstellar medium with the \quokka~GPU-accelerated code, within which we systematically vary the galaxy gas surface density, metallicity, and the scale heights and relative rates of the different nucleosynthetic sources. We show that the metal loadings of galactic winds differ substantially between metals produced by different sources, with typical variations at the level of $\approx 0.3$ dex, a phenomenon we term differential metal loading. Which set of metals suffers preferential loss from this phenomenon varies depending on the galactic environment, and is not easily predictable \textit{a priori}. Our findings call into question the the interpretation of diagnostics of galaxy formation, for example star formation timescales and initial mass functions, based on abundance diagnostics, since the abundance variations upon which these techniques rely are often at levels comparable to those we show can be induced by differential metal loading.

preprint2022arXiv

Quokka: A code for two-moment AMR radiation hydrodynamics on GPUs

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.

preprint2022arXiv

VETTAM: A scheme for radiation hydrodynamics with adaptive mesh refinement using the variable Eddington tensor method

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.

preprint2021arXiv

Exploiting Non-linear Scales in Galaxy-Galaxy Lensing and Galaxy Clustering: A Forecast for the Dark Energy Survey

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.

preprint2019arXiv

Cosmology with galaxy-galaxy lensing on non-perturbative scales: Emulation method and application to BOSS LOWZ

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.