Search

The Nancy Grace Roman Space Telescope, set to launch in 2026, will bring unprecedented precision to measurements of weak gravitational lensing. Because weak lensing is an inherently small signal, it is imperative to minimize systematic errors in measurements as completely as possible; this will ensure that the lensing measurements can be used to their full potential when extracting cosmological information. In this paper, we use laboratory tests of the Roman detectors, simulations of the Roman High Latitude Survey observations, and the proposed Roman image combination pipeline to investigate the magnitude of detector read noise biasing on weak lensing measurements from Roman. First, we combine lab-measured detector noise fields with simulated observations and propagate these images through the Roman image combination pipeline, IMCOM. We characterize the specific signatures of the noise fields in the resultant images and find that noise contributes to the combined images most strongly at scales relevant to physical characteristics of the detector including PSF shape, chip boundaries, and roll angles. We then measure shapes of simulated stars and galaxies and determine the magnitude of noise-induced shear bias on these measurements. We find that star shape correlations satisfy the system noise requirements as defined by the Roman Science Requirements Document. However, for galaxies fainter than $m_{\rm AB}\simeq24$, correction for noise correlations will be needed in order to ensure confidence in shape measurements in any observation band., Comment: 22 pages, 15 figures. Updated to version accepted to PASP

The Nancy Grace Roman Space Telescope will implement a devoted weak gravitational lensing program with its High Latitude Wide Area Survey. For cosmological purposes, a critical step in Roman image processing is to combine dithered undersampled images into unified oversampled images and thus enable high-precision shape measurements. IMCOM is an image coaddition algorithm which offers control over point spread functions in output images. This paper presents the refactored IMCOM software, featuring full object-oriented programming, improved data structures, and alternative linear algebra strategies for determining coaddition weights. Combining these improvements and other acceleration measures, to produce almost equivalent coadded images, the consumption of core-hours has been reduced by about an order of magnitude. We then re-coadd a $16 \times 16 \,{\rm arcmin}^2$ region of our previous image simulations with three linear algebra kernels in four bands, and compare the results in terms of IMCOM optimization goals, properties of coadded noise frames, and measurements of simulated stars. The Cholesky kernel is efficient and relatively accurate, yet its irregular windows for input pixels slightly bias coaddition results. The iterative kernel avoids this issue by tailoring input pixel selection for each output pixel; it yields better noise control, but can be limited by random errors due to finite tolerance. The empirical kernel coadds images using an empirical relation based on geometry; it is inaccurate, but being much faster, it provides a valid option for "quick look" purposes. We fine-tune IMCOM hyperparameters in a companion paper., Comment: 38 pages, 22 figures, submitted to ApJSupp

In the third APOKASC catalog, we present data for the complete sample of 15,808 evolved stars with APOGEE spectroscopic parameters and Kepler asteroseismology. We used ten independent asteroseismic analysis techniques and anchor our system on fundamental radii derived from Gaia $L$ and spectroscopic $T_{\rm eff}$. We provide evolutionary state, asteroseismic surface gravity, mass, radius, age, and the spectroscopic and asteroseismic measurements used to derive them for 12,418 stars. This includes 10,036 exceptionally precise measurements, with median fractional uncertainties in \nmax, \dnu, mass, radius and age of 0.6\%, 0.6\%, 3.8\%, 1.8\%, and 11.1\% respectively. We provide more limited data for 1,624 additional stars which either have lower quality data or are outside of our primary calibration domain. Using lower red giant branch (RGB) stars, we find a median age for the chemical thick disk of $9.14 \pm 0.05 ({\rm ran}) \pm 0.9 ({\rm sys})$ Gyr with an age dispersion of 1.1 Gyr, consistent with our error model. We calibrate our red clump (RC) mass loss to derive an age consistent with the lower RGB and provide asymptotic GB and RGB ages for luminous stars. We also find a sharp upper age boundary in the chemical thin disk. We find that scaling relations are precise and accurate on the lower RGB and RC, but they become more model dependent for more luminous giants and break down at the tip of the RGB. We recommend the usage of multiple methods, calibration to a fundamental scale, and the usage of stellar models to interpret frequency spacings., Comment: 43 pages, 25 figures, submitted ApJSupp. Comments welcome. Data tables available on request from pinsonneault.1@osu.edu

We present the ContEvol (continuous evolution) formalism, a family of implicit numerical methods which only need to solve linear equations and are almost symplectic. Combining values and derivatives of functions, ContEvol outputs allow users to recover full history and render full distributions. Using classic harmonic oscillator as a prototype case, we show that ContEvol methods lead to lower-order errors than two commonly used Runge--Kutta methods. Applying first-order ContEvol to simple celestial mechanics problems, we demonstrate that deviation from equation(s) of motion of ContEvol tracks is still $\mathcal{O}(h^5)$ ($h$ is the step length) by our definition. Numerical experiments with an eccentric elliptical orbit indicate that first-order ContEvol is a viable alternative to classic Runge--Kutta or the symplectic leapfrog integrator. Solving stationary Schr\"odinger equation in quantum mechanics, we manifest ability of ContEvol to handle boundary value or eigenvalue problems. Important directions for future work, including mathematical foundation, higher dimensions, and technical improvements, are discussed at the end of this article., Comment: 74 pages, 22 figures. No journal submission plan. Comments are welcome

The surface [C/N] ratios of evolved giants are strongly affected by the first dredge-up (FDU) of nuclear-processed material from stellar cores. C and N also have distinct nucleosynthetic origins and serve as diagnostics of mixing and mass loss. We use subgiants to find strong trends in the birth [C/N] with [Fe/H], which differ between the low-$\alpha$ and high-$\alpha$ populations. We demonstrate that these birth trends have a strong impact on the surface abundances after the FDU. This effect is neglected in current stellar models, which use solar-scaled C and N. We map out the FDU as a function of evolutionary state, mass, and composition using a large and precisely measured asteroseismic dataset in first-ascent red giant branch (RGB) and core He-burning, or red clump (RC), stars. We describe the domains where [C/N] is a useful mass diagnostic and find that the RC complements the RGB and extends the range of validity to higher mass. We find evidence for extra mixing on the RGB below [Fe/H]= -0.4, matching literature results, for high-$\alpha$ giants, but there is no clear evidence of mixing in the low-$\alpha$ giants. The predicted signal of mass loss is weak and difficult to detect in our sample. We discuss implications for stellar physics and stellar population applications., Comment: 19 pages, 19 figures

Dynamically relaxed galaxy clusters have long played a role in galaxy cluster studies because it is thought their properties can be reconstructed more precisely and with less systematics. As relaxed clusters are desirable, there exist a plethora of criteria for classifying a galaxy cluster as relaxed. In this work, we examine $9$ commonly used observational and theoretical morphological metrics extracted from $54,000$ Mock-X synthetic X-ray images of galaxy clusters taken from the IllustrisTNG, BAHAMAS and MACSIS simulation suites. We find that the simulated criteria distributions are in reasonable agreement with the observed distributions. Many criteria distributions evolve as a function of redshift, cluster mass, numerical resolution and subgrid physics, limiting the effectiveness of a single relaxation threshold value. All criteria are positively correlated with each other, however, the strength of the correlation is sensitive to redshift, mass and numerical choices. Driven by the intrinsic scatter inherent to all morphological metrics and the arbitrary nature of relaxation threshold values, we find the consistency of relaxed subsets defined by the different metrics to be relatively poor. Therefore, the use of relaxed cluster subsets introduces significant selection effects that are non-trivial to resolve., Comment: 21 pages, 6 figures, submitted to MNRAS. Comments welcome

Surveys in the next decade will deliver large samples of galaxy clusters that transform our understanding of their formation. Cluster astrophysics and cosmology studies will become systematics limited with samples of this magnitude. With known properties, hydrodynamical simulations of clusters provide a vital resource for investigating potential systematics. However, this is only realized if we compare simulations to observations in the correct way. Here we introduce the \textsc{Mock-X} analysis framework, a multiwavelength tool that generates synthetic images from cosmological simulations and derives halo properties via observational methods. We detail our methods for generating optical, Compton-$y$ and X-ray images. Outlining our synthetic X-ray image analysis method, we demonstrate the capabilities of the framework by exploring hydrostatic mass bias for the IllustrisTNG, BAHAMAS and MACSIS simulations. Using simulation derived profiles we find an approximately constant bias $b\approx0.13$ with cluster mass, independent of hydrodynamical method or subgrid physics. However, the hydrostatic bias derived from synthetic observations is mass-dependent, increasing to $b=0.3$ for the most massive clusters. This result is driven by a single temperature fit to a spectrum produced by gas with a wide temperature distribution in quasi-pressure equilibrium. The spectroscopic temperature and mass estimate are biased low by cooler gas dominating the emission, due to its quadratic density dependence. The bias and the scatter in estimated mass remain independent of the numerical method and subgrid physics. Our results are consistent with current observations and future surveys will contain sufficient samples of massive clusters to confirm the mass dependence of the hydrostatic bias., Comment: 20 pages, 9 figures, submitted to MNRAS