Simulating Hydrodynamics in Cosmology with CRK-HACC by Nicholas Frontiere et al. on Wednesday 23 November

We introduce CRK-HACC, an extension of the Hardware/Hybrid Accelerated

Cosmology Code (HACC), to resolve gas hydrodynamics in large-scale structure

formation simulations of the universe. The new framework couples the HACC

gravitational N-body solver with a modern smoothed particle hydrodynamics (SPH)

approach called CRKSPH. $\underline{\text{C}}$onservative

$\underline{\text{R}}$eproducing $\underline{\text{K}}$ernel

$\underline{\text{SPH}}$ utilizes smoothing functions that exactly interpolate

linear fields while manifestly preserving conservation laws (momentum, mass,

and energy). The CRKSPH method has been incorporated to accurately model

baryonic effects in cosmology simulations - an important addition targeting the

generation of precise synthetic sky predictions for upcoming observational

surveys. CRK-HACC inherits the codesign strategies of the HACC solver and is

built to run on modern GPU-accelerated supercomputers. In this work, we

summarize the primary solver components and present a number of standard

validation tests to demonstrate code accuracy, including idealized hydrodynamic

and cosmological setups, as well as self-similarity measurements.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2202.02840v2

Redshift requirements for cosmic shear with intrinsic alignment by Silvan Fischbacher et al. on Wednesday 23 November

Intrinsic alignment (IA) modelling and photometric redshift estimation are

two of the main sources of systematic uncertainty in weak lensing surveys. We

investigate the impact of redshift errors and their interplay with different IA

models. Generally, errors on the mean $\delta_z$ and on the width $\sigma_z$ of

the redshift bins can both lead to biases in cosmological constraints. We find

that such biases can, however, only be partially resolved by marginalizing over

$\delta_z$ and $\sigma_z$. For Stage-III surveys, $\delta_z$ and $\sigma_z$

cannot be well constrained due to limited statistics. The resulting biases are

thus sensitive to prior volume effects. For Stage-IV surveys, we observe that

marginalizing over the redshift parameters has an impact and reduces the bias.

We derive requirements on the uncertainty of $\sigma_z$ and $\delta_z$ for both

Stage-III and Stage-IV surveys. We assume that the redshift systematic errors

on $S_8$ should be less than half of the statistical errors, and the median

bias should be smaller than $0.25\sigma$. We find that the uncertainty on

$\delta_z$ has to be $\lesssim0.025$ for the NLA IA model with a Stage-III

survey. For $\sigma_z$, the requirement is met even for large uncertainties

$\leq0.3$. For the TATT IA model, the uncertainty on $\delta_z$ has to be

$\lesssim0.02$ and the uncertainty on $\sigma_z$ has to be $\lesssim0.2$. For

Stage-IV surveys, the uncertainty on $\delta_z$ has to be $\lesssim0.005$ and

the uncertainty on $\sigma_z$ should be $\lesssim0.1$, with no significant

dependence on the IA model. This required high precision will be a challenge

for the redshift calibration of these future surveys. Finally, we investigate

whether the interplay between redshift systematics and IA modelling can explain

the $S_8$-tension between cosmic shear results and CMB measurements. We find

that this is unlikely to explain the current $S_8$-tension.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2207.01627v2

Dark matter halos and scaling relations of extremely massive spiral galaxies from extended HI rotation curves by Enrico M. Di Teodoro et al. on Wednesday 23 November

We present new and archival atomic hydrogen (\hi) observations of \galnum\ of

the most massive spiral galaxies in the local Universe ($M_\star>10^{11} \,

\mathrm{M}_\odot$). From 3D kinematic modeling of the datacubes, we derive

extended \hi\ rotation curves, and from these, we estimate masses of the dark

matter halos and specific angular momenta of the discs. We confirm that massive

spiral galaxies lie at the upper ends of the Tully-Fisher relation (mass vs

velocity, $M \propto V^{4}$) and Fall relation (specific angular momentum vs

mass, $j \propto M^{0.6}$), in both stellar and baryonic forms, with no

significant deviations from single power laws. We study the connections between

baryons and dark matter through the stellar (and baryon)-to-halo ratios of mass

$f_\mathrm{M} \equiv M_\star/M_\mathrm{h}$ and specific angular momentum

$f_\mathrm{j} \equiv j_\star/j_\mathrm{h}$ and $f_\mathrm{j,bar} \equiv

j_\mathrm{bar}/j_\mathrm{h}$. Combining our sample with others from the

literature for less massive disc-dominated galaxies, we find that

$f_\mathrm{M}$ rises monotonically with $M_\star$ and $M_\mathrm{h}$ (instead

of the inverted-U shaped $f_\mathrm{M}$ for spheroid-dominated galaxies), while

$f_\mathrm{j}$ and $f_\mathrm{j,bar}$ are essentially constant near unity over

four decades in mass. Our results indicate that disc galaxies constitute a

self-similar population of objects closely linked to the self-similarity of

their dark halos. This picture is reminiscent of early analytical models of

galaxy formation wherein discs grow by relatively smooth and gradual inflow,

isolated from disruptive events such as major mergers and strong AGN feedback,

in contrast to the more chaotic growth of spheroids.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2207.02906v2

3D Simulations of Strongly Magnetised Non-Rotating Supernovae: Explosion Dynamics and Remnant Properties by Vishnu Varma et al. on Wednesday 23 November

We investigate the impact of strong initial magnetic fields in core-collapse

supernovae of non-rotating progenitors by simulating the collapse and explosion

of a 16.9 Msun star for a strong- and weak-field case assuming a twisted-torus

field with initial central field strengths of ~1012 G and ~106 G. The

strong-field model has been set up with a view to the fossil-field scenario for

magnetar formation and emulates a pre-collapse field configuration that may

occur in massive stars formed by a merger. This model undergoes shock revival

already 100 s after bounce and reaches an explosion energy of 9.3 x 10^50 erg

at 310 ms, in contrast to a more delayed and less energetic explosion in the

weak-field model. The strong magnetic fields help trigger a neutrino-driven

explosion early on, which results in a rapid rise and saturation of the

explosion energy. Dynamically, the strong initial field leads to a fast

build-up of magnetic fields in the gain region to 40% of kinetic equipartition

and also creates sizable pre-shock ram pressure perturbations that are known to

be conducive to asymmetric shock expansion. For the strong-field model, we find

an extrapolated neutron star kick of ~350 km/s, a spin period of ~70 ms, and no

spin-kick alignment. The dipole field strength of the proto-neutron star is 2 x

10^14 G by the end of the simulation with a declining trend. Surprisingly, the

surface dipole field in the weak-field model is stronger, which argues against

a straightforward connection between pre-collapse fields and the birth magnetic

fields of neutron stars.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2204.11009v3

The Science Performance of JWST as Characterized in Commissioning by Jane Rigby et al. on Wednesday 23 November

This paper characterizes the actual science performance of the James Webb

Space Telescope (JWST), as determined from the six month commissioning period.

We summarize the performance of the spacecraft, telescope, science instruments,

and ground system, with an emphasis on differences from pre-launch

expectations. Commissioning has made clear that JWST is fully capable of

achieving the discoveries for which it was built. Moreover, almost across the

board, the science performance of JWST is better than expected; in most cases,

JWST will go deeper faster than expected. The telescope and instrument suite

have demonstrated the sensitivity, stability, image quality, and spectral range

that are necessary to transform our understanding of the cosmos through

observations spanning from near-earth asteroids to the most distant galaxies.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2207.05632v3

Using wavelets to capture deviations from smoothness in galaxy-scale strong lenses by Aymeric Galan et al. on Wednesday 23 November

Modeling the mass distribution of galaxy-scale strong gravitational lenses is

a task of increasing difficulty. The high-resolution and depth of imaging data

now available render simple analytical forms ineffective at capturing lens

structures spanning a large range in spatial scale, mass scale, and morphology.

In this work, we address the problem with a novel multiscale method based on

wavelets. We tested our method on simulated Hubble Space Telescope (HST)

imaging data of strong lenses containing the following different types of mass

substructures making them deviate from smooth models: (1) a localized small

dark matter subhalo, (2) a Gaussian random field (GRF) that mimics a

nonlocalized population of subhalos along the line of sight, and (3)

galaxy-scale multipoles that break elliptical symmetry. We show that wavelets

are able to recover all of these structures accurately. This is made

technically possible by using gradient-informed optimization based on automatic

differentiation over thousands of parameters, which also allow us to sample the

posterior distributions of all model parameters simultaneously. By

construction, our method merges the two main modeling paradigms - analytical

and pixelated - with machine-learning optimization techniques into a single

modular framework. It is also well-suited for the fast modeling of large

samples of lenses. All methods presented here are publicly available in our new

Herculens package.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2207.05763v2

Short GRB Host Galaxies II A Legacy Sample of Redshifts, Stellar Population Properties, and Implications for their Neutron Star Merger Origins by Anya E. Nugent et al. on Wednesday 23 November

We present the stellar population properties of 69 short gamma-ray burst

(GRB) host galaxies, representing the largest uniformly-modeled sample to-date.

Using the Prospector stellar population inference code, we jointly fit

photometry and/or spectroscopy of each host galaxy. We find a population median

redshift of $z=0.64^{+0.83}_{-0.32}$ ($68\%$ confidence), including 10 new or

revised photometric redshifts at $z\gtrsim1$. We further find a median

mass-weighted age of $t_m=0.8^{+2.71}_{-0.53}$Gyr, stellar mass of

$\log(M_*/M_\odot)=9.69^{+0.75}_{-0.65}$, star formation rate of

SFR=$1.44^{+9.37}_{-1.35}M_\odot$yr$^{-1}$, stellar metallicity of

$\log(Z_*/Z_\odot)=-0.38^{+0.44}_{-0.42}$, and dust attenuation of

$A_V=0.43^{+0.85}_{-0.36}$~mag (68\% confidence). Overall, the majority of

short GRB hosts are star-forming ($\approx84\%$), with small fractions that are

either transitioning ($\approx6\%$) or quiescent ($\approx10\%$); however, we

observe a much larger fraction ($\approx40\%$) of quiescent and transitioning

hosts at $z\lesssim0.25$, commensurate with galaxy evolution. We find that

short GRB hosts populate the star-forming main sequence of normal field

galaxies, but do not include as many high-mass galaxies, implying that their

binary neutron star (BNS) merger progenitors are dependent on a combination of

host star formation and stellar mass. The distribution of ages and redshifts

implies a broad delay-time distribution, with a fast-merging channel at $z>1$

and a decreased BNS formation efficiency at lower redshifts. If short GRB hosts

are representative of BNS merger hosts within the horizon of current

gravitational wave detectors, these results can inform future searches for

electromagnetic counterparts. All of the data and modeling products are

available on the BRIGHT website.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2206.01764v3

Constraints on the Optical Depth to Reionization from Balloon-Borne CMB Measurements by Josquin Errard et al. on Wednesday 23 November

We assess the uncertainty with which a balloon-borne experiment, nominally

called Tau Surveyor ($\tau S$), can measure the optical depth to reionization

$\sigma(\tau)$ with given realistic constraints of instrument noise and

foreground emissions. Using a $\tau S$ fiducial design with six frequency bands

between 150 and 380 GHz with white and uniform map noise of 7 $\mu$K arcmin,

achievable with a single mid-latitude flight, and including Planck's 30 and 44

GHz data we assess the error $\sigma(\tau)$ obtained with three foreground

models and as a function of sky fraction $f_{\rm sky}$ between 40% and 54%. We

carry out the analysis using both parametric and blind foreground separation

techniques. We compare $\sigma(\tau)$ values to those obtained with low

frequency and high frequency versions of the experiment called $\tau S$-lf and

$\tau S$-hf that have only four and up to eight frequency bands with narrower

and wider frequency coverage, respectively. We find that with $\tau S$ the

lowest constraint is $\sigma(\tau)=0.0034$, obtained for one of the foreground

models with $f_{\rm sky}$=54%. $\sigma(\tau)$ is larger, in some cases by more

than a factor of 2, for smaller sky fractions, with $\tau S$-lf, or as a

function of foreground model. The $\tau S$-hf configuration does not lead to

significantly tighter constraints. Exclusion of the 30 and 44 GHz data, which

give information about synchrotron emission, leads to significant $\tau$

mis-estimates. Decreasing noise by an ambitious factor of 10 while keeping

$f_{\rm sky}$=40% gives $\sigma(\tau) =0.0031$. The combination of

$\sigma(\tau) =0.0034$, BAO data from DESI, and future CMB B-mode lensing data

from CMB-S3/S4 experiments could give $\sigma(\sum m_{\nu}) = 17$ meV.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2206.03389v2

The Rotation Period Distribution in the Young Open Cluster NGC 6709 by E. M. Cole-Kodikara et al. on Wednesday 23 November

Open clusters serve as a useful tool for calibrating models of the

relationship between mass, rotation, and age for stars with an outer convection

zone due to the homogeneity of the stars within the cluster. Cluster to cluster

comparisons are essential to determine whether the universality of spin down

relations holds. NGC 6709 is selected as a young open cluster for which no

rotation periods of members have previously been obtained. This cluster is at a

distance of over 1 kpc and has two red giant members. Isochrones place the age

of the cluster at around 150 Myr, or approximately the same age as the

Pleiades. Photometry is obtained over a multi-month observing season at the

robotic observatory STELLA. After basic processing, PSF photometry was derived

using Daophot II, and a suite of related software allowed us to create time

series of relative magnitude changes for each star. Four time series analysis

methods are then applied to these light curves to obtain rotation periods for

members stars. We obtain for the first time rotation periods for 45 FGK cluster

members of NGC 6709. We compare our rotation periods to Gaia EDR3 colors and

find a slow-rotating sequence with increasing rotation periods towards redder

stars and a smaller clump of rapid rotators that have not yet joined this

sequence. NGC 6709 has rotation periods very similar to that of another

Pleiades-age open cluster, NGC 2516.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2207.11063v2

EMPRESS IX Extremely Metal-Poor Galaxies are Very Gas-Rich Dispersion-Dominated Systems: Will JWST Witness Gaseous Turbulent High-z Primordial Galaxies? by Yuki Isobe et al. on Wednesday 23 November

We present kinematics of 6 local extremely metal-poor galaxies (EMPGs) with

low metallicities ($0.016-0.098\ Z_{\odot}$) and low stellar masses

($10^{4.7}-10^{7.6} M_{\odot}$). Taking deep medium-high resolution

($R\sim7500$) integral-field spectra with 8.2-m Subaru, we resolve the small

inner velocity gradients and dispersions of the EMPGs with H$\alpha$ emission.

Carefully masking out sub-structures originated by inflow and/or outflow, we

fit 3-dimensional disk models to the observed H$\alpha$ flux, velocity, and

velocity-dispersion maps. All the EMPGs show rotational velocities ($v_{\rm

rot}$) of 5--23 km s$^{-1}$ smaller than the velocity dispersions

($\sigma_{0}$) of 17--31 km s$^{-1}$, indicating dispersion-dominated ($v_{\rm

rot}/\sigma_{0}=0.29-0.80<1$) systems affected by inflow and/or outflow. Except

for two EMPGs with large uncertainties, we find that the EMPGs have very large

gas-mass fractions of $f_{\rm gas}\simeq 0.9-1.0$. Comparing our results with

other H$\alpha$ kinematics studies, we find that $v_{\rm rot}/\sigma_{0}$

decreases and $f_{\rm gas}$ increases with decreasing metallicity, decreasing

stellar mass, and increasing specific star-formation rate. We also find that

simulated high-$z$ ($z\sim 7$) forming galaxies have gas fractions and dynamics

similar to the observed EMPGs. Our EMPG observations and the simulations

suggest that primordial galaxies are gas-rich dispersion-dominated systems,

which would be identified by the forthcoming James Webb Space Telescope (JWST)

observations at $z\sim 7$.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2206.04709v3