Showing 1–50 of 488 results for author: Weinberg, D

Modeling the Ages and Chemical Abundances of Elliptical Galaxies

Authors: Nicole Marcelina Gountanis, David H. Weinberg, Aliza G. Beverage, Nathan R. Sandford, Charlie Conroy, Mariska Kriek

Abstract: Spectroscopic studies of elliptical galaxies show that their stellar population ages, mean metallicity, and $α$-enhancement traced by [Mg/Fe] all increase with galaxy stellar mass or velocity dispersion. We use one-zone galactic chemical evolution (GCE) models with a flexible star formation history (SFH) to model the age, [Mg/H], and [Mg/Fe] inferred from simple stellar population (SSP) fits to ob… ▽ More Spectroscopic studies of elliptical galaxies show that their stellar population ages, mean metallicity, and $α$-enhancement traced by [Mg/Fe] all increase with galaxy stellar mass or velocity dispersion. We use one-zone galactic chemical evolution (GCE) models with a flexible star formation history (SFH) to model the age, [Mg/H], and [Mg/Fe] inferred from simple stellar population (SSP) fits to observed ellipticals at $z \sim 0$ and $z \sim 0.7$. We show that an SSP fit to the spectrum computed from a full GCE model gives ages and abundances close to the light-weighted, logarithmically averaged values of the composite stellar population, <age>, <[Mg/H]>, and <[Mg/Fe]>. With supernova Mg and Fe yields fixed to values motivated by Milky Way stellar populations, we find that predicted <[Mg/H]>-<age> and <[Mg/Fe]>-<age> relations are surprisingly insensitive to SFH parameters: older galaxies have higher <[Mg/Fe]>, but the detailed form of the SFH has limited impact. The star formation efficiency and outflow efficiency affect the early and late evolution of <[Mg/H]>, respectively; explaining observed trends requires higher star formation efficiency and lower outflows in more massive galaxies. With core collapse supernova yields calibrated to the plateau [Mg/Fe]$_{\rm cc} \approx0.45$ observed in many Milky Way studies, our models underpredict the observed <[Mg/Fe]> ratios of ellipticals by 0.05-0.1 dex. Increasing the core collapse yield ratio to [Mg/Fe]$_{\rm cc} = 0.55$ improves the agreement, though the models still lie below the data. We discuss potential resolutions of this discrepancy, including the possibility that many ellipticals terminate their star formation with a self-enriching, terminating burst that reduces the light-weighted age and boosts <[Mg/Fe]>. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Carbon and Iron Deficiencies in Quiescent Galaxies at z=1-3 from JWST-SUSPENSE: Implications for the Formation Histories of Massive Galaxies

Authors: Aliza G. Beverage, Martje Slob, Mariska Kriek, Charlie Conroy, Guillermo Barro, Rachel Bezanson, Gabriel Brammer, Chloe M. Cheng, Anna de Graaff, Natascha M. Förster Schreiber, Marijn Franx, Brian Lorenz, Pavel E. Mancera Piña, Danilo Marchesini, Adam Muzzin, Andrew B. Newman, Sedona H. Price, Alice E. Shapley, Mauro Stefanon, Katherine A. Suess, Pieter van Dokkum, David Weinberg, Daniel R. Weisz

Abstract: We present the stellar metallicities and multi-element abundances (C, Mg, Si, Ca, Ti, Cr, and Fe) of 15 massive (log M/M$_\odot$=10.2-11.2) quiescent galaxies at z=1-3, derived from ultradeep JWST-SUSPENSE spectra. Compared to quiescent galaxies at z~0, these galaxies exhibit a deficiency of 0.25 dex in [C/H], 0.16 dex in [Fe/H], and 0.07 dex in [Mg/H], implying rapid formation and quenching befor… ▽ More We present the stellar metallicities and multi-element abundances (C, Mg, Si, Ca, Ti, Cr, and Fe) of 15 massive (log M/M$_\odot$=10.2-11.2) quiescent galaxies at z=1-3, derived from ultradeep JWST-SUSPENSE spectra. Compared to quiescent galaxies at z~0, these galaxies exhibit a deficiency of 0.25 dex in [C/H], 0.16 dex in [Fe/H], and 0.07 dex in [Mg/H], implying rapid formation and quenching before significant enrichment from asymptotic giant branch stars and Type Ia supernovae. Additionally, we find that galaxies that form at higher redshift have higher [Mg/Fe] and lower [Fe/H] and [Mg/H], irrespective of their observed redshift. The evolution in [Fe/H] and [C/H] is therefore primarily explained by lower redshift samples naturally including galaxies with longer star-formation timescales. On the other hand, the lower [Mg/H] can be explained by galaxies forming at earlier epochs expelling larger gas reservoirs during their quenching phase. Consequently, the mass-metallicity relation, primarily reflecting [Mg/H], is also lower at z=1-3 compared to the lower redshift relation, though the slopes are similar. Finally, we compare our results to standard stellar population modeling approaches employing solar abundance patterns and non-parametric star-formation histories (using Prospector). Our SSP-equivalent ages agree with the mass-weighted ages from Prospector, while the metallicities disagree significantly. Nonetheless, the metallicities better reflect [Fe/H] than total [Z/H]. We also find that star-formation timescales inferred from elemental abundances are significantly shorter than those from Prospector, and we discuss the resulting implications for the early formation of massive galaxies. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: Submitted to ApJ; 18 pages, 6 figures, 1 table

Static Subspace Approximation for Random Phase Approximation Correlation Energies: Implementation and Performance

Authors: Daniel Weinberg, Olivia A. Hull, Jacob M. Clary, Ravishankar Sundararaman, Derek Vigil-Fowler, Mauro Del Ben

Abstract: Developing theoretical understanding of complex reactions and processes at interfaces requires using methods that go beyond semilocal density functional theory to accurately describe the interactions between solvent, reactants and substrates. Methods based on many-body perturbation theory, such as the random phase approximation (RPA), have previously been limited due to their computational complex… ▽ More Developing theoretical understanding of complex reactions and processes at interfaces requires using methods that go beyond semilocal density functional theory to accurately describe the interactions between solvent, reactants and substrates. Methods based on many-body perturbation theory, such as the random phase approximation (RPA), have previously been limited due to their computational complexity. However, this is now a surmountable barrier due to the advances in computational power available, in particular through modern GPU-based supercomputers. In this work, we describe the implementation of RPA calculations within BerkeleyGW and show its favorable computational performance on large complex systems relevant for catalysis and electrochemistry applications. Our implementation builds off of the static subspace approximation which, by employing a compressed representation of the frequency dependent polarizability, enables the evaluation of the RPA correlation energy with significant acceleration and systematically controllable accuracy. We find that the computational cost of calculating the RPA correlation energy scales only linearly with system size for systems containing up to 50 thousand bands, and is expected to scale quadratically thereafter. We also show excellent strong scaling results across several supercomputers, demonstrating the performance and portability of this implementation. △ Less

Submitted 30 May, 2024; originally announced May 2024.

Comments: 10 pages, 5 figures

CMB lensing and Lyα forest cross bispectrum from DESI's first-year quasar sample

Authors: N. G. Karaçaylı, P. Martini, D. H. Weinberg, S. Ferraro, R. de Belsunce, J. Aguilar, S. Ahlen, E. Armengaud, D. Brooks, T. Claybaugh, A. de la Macorra, B. Dey, P. Doel, K. Fanning, J. E. Forero-Romero, S. Gontcho A Gontcho, A. X. Gonzalez-Morales, G. Gutierrez, J. Guy, K. Honscheid, D. Kirkby, T. Kisner, A. Kremin, A. Lambert, M. Landriau , et al. (28 additional authors not shown)

Abstract: The squeezed cross-bispectrum \bispeconed\ between the gravitational lensing in the Cosmic Microwave Background and the 1D \lya\ forest power spectrum can constrain bias parameters and break degeneracies between $σ_8$ and other cosmological parameters. We detect \bispeconed\ with $4.8σ$ significance at an effective redshift $z_\mathrm{eff}=2.4$ using Planck PR3 lensing map and over 280,000 quasar… ▽ More The squeezed cross-bispectrum \bispeconed\ between the gravitational lensing in the Cosmic Microwave Background and the 1D \lya\ forest power spectrum can constrain bias parameters and break degeneracies between $σ_8$ and other cosmological parameters. We detect \bispeconed\ with $4.8σ$ significance at an effective redshift $z_\mathrm{eff}=2.4$ using Planck PR3 lensing map and over 280,000 quasar spectra from the Dark Energy Spectroscopic Instrument's first-year data. We test our measurement against metal contamination and foregrounds such as Galactic extinction and clusters of galaxies by deprojecting the thermal Sunyaev-Zeldovich effect. We compare our results to a tree-level perturbation theory calculation and find reasonable agreement between the model and measurement. △ Less

Submitted 23 May, 2024; originally announced May 2024.

Comments: 13 pages excluding references, 8 figures

Modeling the Galactic Chemical Evolution of Helium

Authors: Miqaela K. Weller, David H. Weinberg, James W. Johnson

Abstract: We examine the galactic chemical evolution (GCE) of $^4$He in one-zone and multi-zone models, with particular attention to theoretical predictions and empirical constraints on IMF-averaged yields. Published models of massive star winds and core collapse supernovae span a factor of 2 -- 3 in the IMF-averaged $^4$He yield, $y\mathrm{_{He}^{CC}}$. Published models of intermediate mass, asymptotic gia… ▽ More We examine the galactic chemical evolution (GCE) of $^4$He in one-zone and multi-zone models, with particular attention to theoretical predictions and empirical constraints on IMF-averaged yields. Published models of massive star winds and core collapse supernovae span a factor of 2 -- 3 in the IMF-averaged $^4$He yield, $y\mathrm{_{He}^{CC}}$. Published models of intermediate mass, asymptotic giant branch (AGB) stars show better agreement on the IMF-averaged yield, $y\mathrm{_{He}^{AGB}}$, and they predict that more than half of this yield comes from stars with $M=4-8 M_\odot$, making AGB $^4$He enrichment rapid compared to Fe enrichment from Type Ia supernovae. Although our GCE models include many potentially complicating effects, the short enrichment time delay and mild metallicity dependence of the predicted yields makes the results quite simple: across a wide range of metallicity and age, the non-primordial $^4$He mass fraction $ΔY = Y-Y_{\mathrm{P}}$ is proportional to the abundance of promptly produced $α$-elements, like oxygen, with $ΔY/Z_{\mathrm{O}} \approx (y\mathrm{_{He}^{CC}}+y\mathrm{_{He}^{AGB}})/y\mathrm{_{O}^{CC}}$. Reproducing solar abundances with our fiducial choice of the oxygen yield $y\mathrm{_{O}^{CC}}=0.0071$ implies $y\mathrm{_{He}^{CC}}+y\mathrm{_{He}^{AGB}} \approx 0.022$, i.e., $0.022M_\odot$ of net $^4$He production per solar mass of star formation. Our GCE models with this yield normalization are consistent with most available observations, though the implied $y\mathrm{_{He}^{CC}}$ is low compared to most of the published massive star models. More precise measurements of $ΔY$ in stars and gas across a wide range of metallicity and [$α$/Fe] ratio could test our models more stringently, either confirming the simple picture suggested by our calculations or revealing surprises in the evolution of the second most abundant element. △ Less

Submitted 12 April, 2024; originally announced April 2024.

Comments: 19 Pages, 18 Figures, submitted to MNRAS

DESI 2024 VI: Cosmological Constraints from the Measurements of Baryon Acoustic Oscillations

Authors: DESI Collaboration, A. G. Adame, J. Aguilar, S. Ahlen, S. Alam, D. M. Alexander, M. Alvarez, O. Alves, A. Anand, U. Andrade, E. Armengaud, S. Avila, A. Aviles, H. Awan, B. Bahr-Kalus, S. Bailey, C. Baltay, A. Bault, J. Behera, S. BenZvi, A. Bera, F. Beutler, D. Bianchi, C. Blake, R. Blum , et al. (178 additional authors not shown)

Abstract: We present cosmological results from the measurement of baryon acoustic oscillations (BAO) in galaxy, quasar and Lyman-$α$ forest tracers from the first year of observations from the Dark Energy Spectroscopic Instrument (DESI), to be released in the DESI Data Release 1. DESI BAO provide robust measurements of the transverse comoving distance and Hubble rate, or their combination, relative to the s… ▽ More We present cosmological results from the measurement of baryon acoustic oscillations (BAO) in galaxy, quasar and Lyman-$α$ forest tracers from the first year of observations from the Dark Energy Spectroscopic Instrument (DESI), to be released in the DESI Data Release 1. DESI BAO provide robust measurements of the transverse comoving distance and Hubble rate, or their combination, relative to the sound horizon, in seven redshift bins from over 6 million extragalactic objects in the redshift range $0.1<z<4.2$. DESI BAO data alone are consistent with the standard flat $Λ$CDM cosmological model with a matter density $Ω_\mathrm{m}=0.295\pm 0.015$. Paired with a BBN prior and the robustly measured acoustic angular scale from the CMB, DESI requires $H_0=(68.52\pm0.62)$ km/s/Mpc. In conjunction with CMB anisotropies from Planck and CMB lensing data from Planck and ACT, we find $Ω_\mathrm{m}=0.307\pm 0.005$ and $H_0=(67.97\pm0.38)$ km/s/Mpc. Extending the baseline model with a constant dark energy equation of state parameter $w$, DESI BAO alone require $w=-0.99^{+0.15}_{-0.13}$. In models with a time-varying dark energy equation of state parametrized by $w_0$ and $w_a$, combinations of DESI with CMB or with SN~Ia individually prefer $w_0>-1$ and $w_a<0$. This preference is 2.6$σ$ for the DESI+CMB combination, and persists or grows when SN~Ia are added in, giving results discrepant with the $Λ$CDM model at the $2.5σ$, $3.5σ$ or $3.9σ$ levels for the addition of Pantheon+, Union3, or DES-SN5YR datasets respectively. For the flat $Λ$CDM model with the sum of neutrino mass $\sum m_ν$ free, combining the DESI and CMB data yields an upper limit $\sum m_ν< 0.072$ $(0.113)$ eV at 95% confidence for a $\sum m_ν>0$ $(\sum m_ν>0.059)$ eV prior. These neutrino-mass constraints are substantially relaxed in models beyond $Λ$CDM. [Abridged.] △ Less

Submitted 24 April, 2024; v1 submitted 3 April, 2024; originally announced April 2024.

Comments: This DESI Collaboration Key Publication is part of the 2024 publication series using the first year of observations (see https://data.desi.lbl.gov/doc/papers). Typos corrected and a new figure and discussion added to Appendix A

DESI 2024 IV: Baryon Acoustic Oscillations from the Lyman Alpha Forest

Authors: DESI Collaboration, A. G. Adame, J. Aguilar, S. Ahlen, S. Alam, D. M. Alexander, M. Alvarez, O. Alves, A. Anand, U. Andrade, E. Armengaud, S. Avila, A. Aviles, H. Awan, S. Bailey, C. Baltay, A. Bault, J. Bautista, J. Behera, S. BenZvi, F. Beutler, D. Bianchi, C. Blake, R. Blum, S. Brieden , et al. (174 additional authors not shown)

Abstract: We present the measurement of Baryon Acoustic Oscillations (BAO) from the Lyman-$α$ (Ly$α$) forest of high-redshift quasars with the first-year dataset of the Dark Energy Spectroscopic Instrument (DESI). Our analysis uses over $420\,000$ Ly$α$ forest spectra and their correlation with the spatial distribution of more than $700\,000$ quasars. An essential facet of this work is the development of a… ▽ More We present the measurement of Baryon Acoustic Oscillations (BAO) from the Lyman-$α$ (Ly$α$) forest of high-redshift quasars with the first-year dataset of the Dark Energy Spectroscopic Instrument (DESI). Our analysis uses over $420\,000$ Ly$α$ forest spectra and their correlation with the spatial distribution of more than $700\,000$ quasars. An essential facet of this work is the development of a new analysis methodology on a blinded dataset. We conducted rigorous tests using synthetic data to ensure the reliability of our methodology and findings before unblinding. Additionally, we conducted multiple data splits to assess the consistency of the results and scrutinized various analysis approaches to confirm their robustness. For a given value of the sound horizon ($r_d$), we measure the expansion at $z_{\rm eff}=2.33$ with 2\% precision, $H(z_{\rm eff}) = (239.2 \pm 4.8) (147.09~{\rm Mpc} /r_d)$ km/s/Mpc. Similarly, we present a 2.4\% measurement of the transverse comoving distance to the same redshift, $D_M(z_{\rm eff}) = (5.84 \pm 0.14) (r_d/147.09~{\rm Mpc})$ Gpc. Together with other DESI BAO measurements at lower redshifts, these results are used in a companion paper to constrain cosmological parameters. △ Less

Submitted 12 April, 2024; v1 submitted 3 April, 2024; originally announced April 2024.

Comments: This DESI Collaboration Key Publication is part of the 2024 publication series using the first year of observations (see https://data.desi.lbl.gov/doc/papers)

DESI 2024 III: Baryon Acoustic Oscillations from Galaxies and Quasars

Authors: DESI Collaboration, A. G. Adame, J. Aguilar, S. Ahlen, S. Alam, D. M. Alexander, M. Alvarez, O. Alves, A. Anand, U. Andrade, E. Armengaud, S. Avila, A. Aviles, H. Awan, S. Bailey, C. Baltay, A. Bault, J. Behera, S. BenZvi, F. Beutler, D. Bianchi, C. Blake, R. Blum, S. Brieden, A. Brodzeller , et al. (171 additional authors not shown)

Abstract: We present the DESI 2024 galaxy and quasar baryon acoustic oscillations (BAO) measurements using over 5.7 million unique galaxy and quasar redshifts in the range 0.1<z<2.1. Divided by tracer type, we utilize 300,017 galaxies from the magnitude-limited Bright Galaxy Survey with 0.1<z<0.4, 2,138,600 Luminous Red Galaxies with 0.4<z<1.1, 2,432,022 Emission Line Galaxies with 0.8<z<1.6, and 856,652 qu… ▽ More We present the DESI 2024 galaxy and quasar baryon acoustic oscillations (BAO) measurements using over 5.7 million unique galaxy and quasar redshifts in the range 0.1<z<2.1. Divided by tracer type, we utilize 300,017 galaxies from the magnitude-limited Bright Galaxy Survey with 0.1<z<0.4, 2,138,600 Luminous Red Galaxies with 0.4<z<1.1, 2,432,022 Emission Line Galaxies with 0.8<z<1.6, and 856,652 quasars with 0.8<z<2.1, over a ~7,500 square degree footprint. The analysis was blinded at the catalog-level to avoid confirmation bias. All fiducial choices of the BAO fitting and reconstruction methodology, as well as the size of the systematic errors, were determined on the basis of the tests with mock catalogs and the blinded data catalogs. We present several improvements to the BAO analysis pipeline, including enhancing the BAO fitting and reconstruction methods in a more physically-motivated direction, and also present results using combinations of tracers. We present a re-analysis of SDSS BOSS and eBOSS results applying the improved DESI methodology and find scatter consistent with the level of the quoted SDSS theoretical systematic uncertainties. With the total effective survey volume of ~ 18 Gpc$^3$, the combined precision of the BAO measurements across the six different redshift bins is ~0.52%, marking a 1.2-fold improvement over the previous state-of-the-art results using only first-year data. We detect the BAO in all of these six redshift bins. The highest significance of BAO detection is $9.1σ$ at the effective redshift of 0.93, with a constraint of 0.86% placed on the BAO scale. We find our measurements are systematically larger than the prediction of Planck-2018 LCDM model at z<0.8. We translate the results into transverse comoving distance and radial Hubble distance measurements, which are used to constrain cosmological models in our companion paper [abridged]. △ Less

Submitted 3 April, 2024; originally announced April 2024.

Comments: This DESI Collaboration Key Publication is part of the 2024 publication series using the first year of observations (see https://data.desi.lbl.gov/doc/papers)

Vision-language models for decoding provider attention during neonatal resuscitation

Authors: Felipe Parodi, Jordan Matelsky, Alejandra Regla-Vargas, Elizabeth Foglia, Charis Lim, Danielle Weinberg, Konrad Kording, Heidi Herrick, Michael Platt

Abstract: Neonatal resuscitations demand an exceptional level of attentiveness from providers, who must process multiple streams of information simultaneously. Gaze strongly influences decision making; thus, understanding where a provider is looking during neonatal resuscitations could inform provider training, enhance real-time decision support, and improve the design of delivery rooms and neonatal intensi… ▽ More Neonatal resuscitations demand an exceptional level of attentiveness from providers, who must process multiple streams of information simultaneously. Gaze strongly influences decision making; thus, understanding where a provider is looking during neonatal resuscitations could inform provider training, enhance real-time decision support, and improve the design of delivery rooms and neonatal intensive care units (NICUs). Current approaches to quantifying neonatal providers' gaze rely on manual coding or simulations, which limit scalability and utility. Here, we introduce an automated, real-time, deep learning approach capable of decoding provider gaze into semantic classes directly from first-person point-of-view videos recorded during live resuscitations. Combining state-of-the-art, real-time segmentation with vision-language models (CLIP), our low-shot pipeline attains 91\% classification accuracy in identifying gaze targets without training. Upon fine-tuning, the performance of our gaze-guided vision transformer exceeds 98\% accuracy in gaze classification, approaching human-level precision. This system, capable of real-time inference, enables objective quantification of provider attention dynamics during live neonatal resuscitation. Our approach offers a scalable solution that seamlessly integrates with existing infrastructure for data-scarce gaze analysis, thereby offering new opportunities for understanding and refining clinical decision making. △ Less

Submitted 1 April, 2024; originally announced April 2024.

Comments: 9 pages, 4 figures

The Impact of Feedback-driven Outflows on Bar Formation

Authors: Martin D. Weinberg

Abstract: We investigate the coupling between the temporal variation from galaxy-formation feedback and the bar instability. We show that fluctuations from mass outflow on star-formation time scales affect the radial motion of disk orbits. The resulting incoherence in orbital phase leads to the disruption of the bar-forming dynamics. Bar formation is suppressed in starburst galaxies that have fluctuation ti… ▽ More We investigate the coupling between the temporal variation from galaxy-formation feedback and the bar instability. We show that fluctuations from mass outflow on star-formation time scales affect the radial motion of disk orbits. The resulting incoherence in orbital phase leads to the disruption of the bar-forming dynamics. Bar formation is suppressed in starburst galaxies that have fluctuation time scales within the range 10 Myr to 200 Myr with repeated events with wind mass 15% of the disk within 0.5 scale lengths or 1.4% of the total disk mass. The work done by feedback is capable of reducing the amplitude or, with enough amplitude, destroying an existing bar. AGN feedback with similar amplitude and timescales would have similar behavior. To model the dynamics of the coupling and interpret the results of the full N-body simulations, we introduce a generalization of the Hamiltonian mean-field (HMF) model, drawing inspiration from the Lynden-Bell (1979) mechanism for bar growth. Our non-linear 'BarHMF' model is designed to reproduce linear perturbation theory in the low-amplitude limit. Notably, without star-formation feedback, this model exhibits exponential growth whose rate depends on disk mass and reproduces the expected saturation of bar growth observed in N-body simulations. We describe several promising applications of the BarHMF model beyond this study. △ Less

Submitted 18 March, 2024; originally announced March 2024.

Comments: 23 pages, 17 figures, uses AMS and openjournal styles

Chemical Cartography with APOGEE: Two-process Parameters and Residual Abundances for 288,789 Stars from Data Release 17

Authors: Tawny Sit, David H. Weinberg, Adam Wheeler, Christian R. Hayes, Sten Hasselquist, Thomas Masseron, Jennifer Sobeck

Abstract: Stellar abundance measurements are subject to systematic errors that induce extra scatter and artificial correlations in elemental abundance patterns. We derive empirical calibration offsets to remove systematic trends with surface gravity $\log(g)$ in 17 elemental abundances of 288,789 evolved stars from the SDSS APOGEE survey. We fit these corrected abundances as the sum of a prompt process trac… ▽ More Stellar abundance measurements are subject to systematic errors that induce extra scatter and artificial correlations in elemental abundance patterns. We derive empirical calibration offsets to remove systematic trends with surface gravity $\log(g)$ in 17 elemental abundances of 288,789 evolved stars from the SDSS APOGEE survey. We fit these corrected abundances as the sum of a prompt process tracing core-collapse supernovae and a delayed process tracing Type Ia supernovae, thus recasting each star's measurements into the amplitudes $A_{\text{cc}}$ and $A_{\text{Ia}}$ and the element-by-element residuals from this two-parameter fit. As a first application of this catalog, which is $8\times$ larger than that of previous analyses that used a restricted $\log(g)$ range, we examine the median residual abundances of 14 open clusters, nine globular clusters, and four dwarf satellite galaxies. Relative to field Milky Way disk stars, the open clusters younger than 2 Gyr show $\approx 0.1-0.2$ dex enhancements of the neutron-capture element Ce, and the two clusters younger than 0.5 Gyr also show elevated levels of C+N, Na, S, and Cu. Globular clusters show elevated median abundances of C+N, Na, Al, and Ce, and correlated abundance residuals that follow previously known trends. The four dwarf satellites show similar residual abundance patterns despite their different star formation histories, with $\approx 0.2-0.3$ dex depletions in C+N, Na, and Al and $\approx 0.1$ dex depletions in Ni, V, Mn, and Co. We provide our catalog of corrected APOGEE abundances, two-process amplitudes, and residual abundances, which will be valuable for future studies of abundance patterns in different stellar populations and of additional enrichment processes that affect galactic chemical evolution. △ Less

Submitted 24 May, 2024; v1 submitted 12 March, 2024; originally announced March 2024.

Comments: Updated to ApJ accepted version. 35 pages, 18 figures, 5 tables. Catalog available at https://doi.org/10.5281/zenodo.10659205

Nature vs. Nurture: Distinguishing Effects from Stellar Processing and Chemical Evolution on Carbon and Nitrogen in Red Giant Stars

Authors: John D. Roberts, Marc H. Pinsonneault, Jennifer A. Johnson, Joel C. Zinn, David H. Weinberg, Mathieu Vrard, Jamie Tayar, Dennis Stello, Benoît Mosser, James W. Johnson, Kaili Cao, Keivan G. Stassun, Guy S. Stringfellow, Aldo Serenelli, Savita Mathur, Saskia Hekker, Rafael A. García, Yvonne P. Elsworth, Enrico Corsaro

Abstract: 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-$α$ and high-$α$ populations. We demonstrate that these… ▽ More 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-$α$ and high-$α$ 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-$α$ giants, but there is no clear evidence of mixing in the low-$α$ 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. △ Less

Submitted 5 March, 2024; originally announced March 2024.

Comments: 19 pages, 19 figures

Deciphering baryonic feedback with galaxy clusters

Authors: Chun-Hao To, Shivam Pandey, Elisabeth Krause, Nihar Dalal, Dhayaa Anbajagane, David H. Weinberg

Abstract: Upcoming cosmic shear analyses will precisely measure the cosmic matter distribution at low redshifts. At these redshifts, the matter distribution is affected by galaxy formation physics, primarily baryonic feedback from star formation and active galactic nuclei. Employing measurements from the Magneticum and IllustrisTNG simulations and a dark matter + baryon (DMB) halo model, this paper demonstr… ▽ More Upcoming cosmic shear analyses will precisely measure the cosmic matter distribution at low redshifts. At these redshifts, the matter distribution is affected by galaxy formation physics, primarily baryonic feedback from star formation and active galactic nuclei. Employing measurements from the Magneticum and IllustrisTNG simulations and a dark matter + baryon (DMB) halo model, this paper demonstrates that Sunyaev-Zel'dovich (SZ) effect observations of galaxy clusters, whose masses have been calibrated using weak gravitational lensing, can constrain the baryonic impact on cosmic shear with statistical and systematic errors subdominant to the measurement errors of DES-Y3 and LSST-Y1. We further dissect the contributions from different scales and halos with different masses to cosmic shear, highlighting the dominant role of SZ clusters at scales critical for cosmic shear analyses. These findings suggest a promising avenue for future joint analyses of Cosmic Microwave Background (CMB) and lensing surveys. △ Less

Submitted 31 January, 2024; originally announced February 2024.

Comments: To be submitted to JCAP. Comments are welcome!

Forecasting the constraints on optical selection bias and projection effects of galaxy cluster lensing with multiwavelength data

Authors: Conghao Zhou, Hao-Yi Wu, Andrés N. Salcedo, Sebastian Grandis, Tesla Jeltema, Alexie Leauthaud, Matteo Costanzi, Tomomi Sunayama, David H. Weinberg, Tianyu Zhang, Eduardo Rozo, Chun-Hao To, Sebastian Bocquet, Tamas Varga, Matthew Kwiecien

Abstract: Galaxy clusters identified with optical imaging tend to suffer from projection effects, which impact richness (the number of member galaxies in a cluster) and lensing coherently. Physically unassociated galaxies can be mistaken as cluster members due to the significant uncertainties in their line-of-sight distances, thereby changing the observed cluster richness; at the same time, projection effec… ▽ More Galaxy clusters identified with optical imaging tend to suffer from projection effects, which impact richness (the number of member galaxies in a cluster) and lensing coherently. Physically unassociated galaxies can be mistaken as cluster members due to the significant uncertainties in their line-of-sight distances, thereby changing the observed cluster richness; at the same time, projection effects alter the weak gravitational lensing signals of clusters, leading to a correlated scatter between richness and lensing at a given halo mass. As a result, the lensing signals for optically selected clusters tend to be biased high. This optical selection bias problem of cluster lensing is one of the key challenges in cluster cosmology. Fortunately, recently available multiwavelength observations of clusters provide a solution. We analyze a simulated data set mimicking the observed lensing of clusters identified by both optical photometry and gas properties, aiming to constrain this selection bias. Assuming a redMaPPer sample from the Dark Energy Survey with South Pole Telescope Sunyaev-Zeldovich effect observations, we find that an overlapping survey of 1300 square deg, 0.2 < z < 0.65, can constrain the average lensing bias to an accuracy of 5 percent. This provides an exciting opportunity for directly constraining optical selection bias from observations. We further show that our approach can remove the optical selection bias from the lensing signal, paving the way for future optical cluster cosmology analyses. △ Less

Submitted 18 December, 2023; originally announced December 2023.

Comments: 16 pages, 5 figures. Submitted to PRD

Ultrafast 3-D Super Resolution Ultrasound using Row-Column Array specific Coherence-based Beamforming and Rolling Acoustic Sub-aperture Processing: In Vitro, In Vivo and Clinical Study

Authors: Joseph Hansen-Shearer, Jipeng Yan, Marcelo Lerendegui, Biao Huang, Matthieu Toulemonde, Kai Riemer, Qingyuan Tan, Johanna Tonko, Peter D. Weinberg, Chris Dunsby, Meng-Xing Tang

Abstract: The row-column addressed array is an emerging probe for ultrafast 3-D ultrasound imaging. It achieves this with far fewer independent electronic channels and a wider field of view than traditional 2-D matrix arrays, of the same channel count, making it a good candidate for clinical translation. However, the image quality of row-column arrays is generally poor, particularly when investigating tissu… ▽ More The row-column addressed array is an emerging probe for ultrafast 3-D ultrasound imaging. It achieves this with far fewer independent electronic channels and a wider field of view than traditional 2-D matrix arrays, of the same channel count, making it a good candidate for clinical translation. However, the image quality of row-column arrays is generally poor, particularly when investigating tissue. Ultrasound localisation microscopy allows for the production of super-resolution images even when the initial image resolution is not high. Unfortunately, the row-column probe can suffer from imaging artefacts that can degrade the quality of super-resolution images as `secondary' lobes from bright microbubbles can be mistaken as microbubble events, particularly when operated using plane wave imaging. These false events move through the image in a physiologically realistic way so can be challenging to remove via tracking, leading to the production of 'false vessels'. Here, a new type of rolling window image reconstruction procedure was developed, which integrated a row-column array-specific coherence-based beamforming technique with acoustic sub-aperture processing for the purposes of reducing `secondary' lobe artefacts, noise and increasing the effective frame rate. Using an {\it{in vitro}} cross tube, it was found that the procedure reduced the percentage of `false' locations from $\sim$26\% to $\sim$15\% compared to traditional orthogonal plane wave compounding. Additionally, it was found that the noise could be reduced by $\sim$7 dB and that the effective frame rate could be increased to over 4000 fps. Subsequently, {\it{in vivo}} ultrasound localisation microscopy was used to produce images non-invasively of a rabbit kidney and a human thyroid. △ Less

Submitted 15 November, 2023; originally announced November 2023.

Realistic simulated galaxies form [$α$/Fe]-[Fe/H] knees due to a sustained decline in their star formation rates

Authors: Andrew C. Mason, Robert A. Crain, Ricardo P. Schiavon, David H. Weinberg, Joel Pfeffer, Joop Schaye, Matthieu Schaller, Tom Theuns

Abstract: We examine the stellar [$α$/Fe]-[Fe/H] distribution of $\simeq1000$ present-day galaxies in a high-resolution EAGLE simulation. A slight majority of galaxies exhibit the canonical distribution, characterised by a sequence of low-metallicity stars with high [$α$/Fe] that transitions at a "knee" to a sequence of declining [$α$/Fe] with increasing metallicity. This population yields a knee metallicit… ▽ More We examine the stellar [$α$/Fe]-[Fe/H] distribution of $\simeq1000$ present-day galaxies in a high-resolution EAGLE simulation. A slight majority of galaxies exhibit the canonical distribution, characterised by a sequence of low-metallicity stars with high [$α$/Fe] that transitions at a "knee" to a sequence of declining [$α$/Fe] with increasing metallicity. This population yields a knee metallicity - galaxy-mass relation similar to that observed in Local Group galaxies, both in slope and scatter. However, many simulated galaxies lack a knee or exhibit more complicated distributions. Knees are found only in galaxies with star formation histories (SFHs) featuring a sustained decline from an early peak ($t\simeq7~{\rm Gyr}$), which enables enrichment by Type Ia supernovae to dominate that due to Type II supernovae (SN II), reducing [$α$/Fe] in the interstellar gas. The simulation thus indicates that, contrary to the common interpretation implied by analytic galactic chemical evolution (GCE) models, knee formation is not a consequence of the onset of enrichment by SN Ia. We use the SFH of a simulated galaxy exhibiting a knee as input to the VICE GCE model, finding it yields an $α$-rich plateau enriched only by SN II, but the plateau comprises little stellar mass and the galaxy forms few metal-poor ([Fe/H]$\lesssim - 1$) stars. This follows from the short, constant gas consumption timescale typically assumed by GCEs, which implies the presence of a readily-enriched, low-mass gas reservoir. When an initially longer, evolving consumption timescale is adopted, VICE reproduces the simulated galaxy's track through the [$α$/Fe]-[Fe/H] plane and its metallicity distribution function. △ Less

Submitted 31 October, 2023; originally announced November 2023.

Comments: 17 pages, 12 figures, 1 table, submitted to MNRAS

Dark Energy Survey Year 1 Clusters are Consistent with Planck

Authors: Andrés N. Salcedo, Hao-Yi Wu, Eduardo Rozo, David H. Weinberg, Chun-Hao To, Tomomi Sunayama, Andy Lee

Abstract: The recent Dark Energy Survey Year 1 (DES-Y1) analysis of galaxy cluster abundances and weak lensing produced $Ω_{\rm m}$ and $σ_8$ constraints in 5.6$σ$ tension with Planck. It is suggested in that work that this tension is driven by unmodelled systematics in optical cluster selection. We present a novel simulation-based forward modeling framework that explicitly incorporates cluster selection in… ▽ More The recent Dark Energy Survey Year 1 (DES-Y1) analysis of galaxy cluster abundances and weak lensing produced $Ω_{\rm m}$ and $σ_8$ constraints in 5.6$σ$ tension with Planck. It is suggested in that work that this tension is driven by unmodelled systematics in optical cluster selection. We present a novel simulation-based forward modeling framework that explicitly incorporates cluster selection into its model predictions. Applying this framework to the DES-Y1 data we find consistency with Planck, resolving the tension found in the DES-Y1 analysis. An extension of this approach to the final DES data set will produce robust constraints on $Λ$CDM parameters and correspondingly strong tests of cosmological models. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: 6 pages, 2 figures, 1 table, Supplemental material with 2 figures. Submitted to Physical Review Letters

The Scale of Stellar Yields: Implications of the Measured Mean Iron Yield of Core Collapse Supernovae

Authors: David H. Weinberg, Emily J. Griffith, James W. Johnson, Todd A. Thompson

Abstract: The scale of alpha-element yields is difficult to predict from theory because of uncertainties in massive star evolution, supernova physics, and black hole formation, and it is difficult to constrain empirically because the impact of higher yields can be compensated by greater metal loss in galactic winds. We use a recent measurement of the mean iron yield of core collapse supernovae (CCSN) by Rod… ▽ More The scale of alpha-element yields is difficult to predict from theory because of uncertainties in massive star evolution, supernova physics, and black hole formation, and it is difficult to constrain empirically because the impact of higher yields can be compensated by greater metal loss in galactic winds. We use a recent measurement of the mean iron yield of core collapse supernovae (CCSN) by Rodriguez et al. (RMN23), $\bar{y}_{\rm Fe}^{\rm cc} =0.058 \pm 0.007 M_\odot$, to infer the scale of alpha-element yields by assuming that the plateau of [alpha/Fe] abundance ratios observed in low metallicity stars represents the yield ratio of CCSN. For a Kroupa IMF and a plateau at [alpha/Fe]=0.45, we find that the population-averaged yields of O and Mg per unit mass of star formation are about equal to the mass fractions of these elements in the sun. The inferred O and Fe yields agree with predictions of the Sukhbold et al. (2016) CCSN models assuming their Z9.6+N20 neutrino-driven engine, a scenario in which many progenitors with $M<40M_\odot$ implode to black holes rather than exploding. The yields are lower than assumed in some models of galactic chemical evolution (GCE) and the galaxy mass-metallicity relation, reducing the level of outflows needed to match observed abundances. For straightforward assumptions, we find that one-zone GCE models with mass-loading factor $η\approx 0.6$ evolve to solar metallicity at late times. By requiring that models reach [alpha/Fe]=0 at late times, and assuming a mean Fe yield of $0.7M_\odot$ per Type Ia supernova, we infer a Hubble-time integrated SNIa rate of $1.1\times 10^{-3} M_\odot^{-1}$, compatible with estimates from supernova surveys. The RMN23 measurement provides one of the few empirical anchors for the absolute scale of nucleosynthetic yields, with wide-ranging implications for stellar and galactic astrophysics. △ Less

Submitted 11 September, 2023; originally announced September 2023.

Comments: 19 pp., 9 figs, submitted to AAS journals

A Tale of Two Disks: Mapping the Milky Way with the Final Data Release of APOGEE

Authors: Julie Imig, Cathryn Price, Jon A. Holtzman, Alexander Stone-Martinez, Steven R. Majewski, David H. Weinberg, Jennifer A. Johnson, Carlos Allende Prieto, Rachael L. Beaton, Timothy C. Beers, Dmitry Bizyaev, Michael R. Blanton, Joel R. Brownstein, Katia Cunha, José G. Fernández-Trincado, Diane K. Feuillet, Sten Hasselquist, Christian R. Hayes, Henrik Jönsson, Richard R. Lane, Jianhui Lian, Szabolcs Mészáros, David L. Nidever, Annie C. Robin, Matthew Shetrone , et al. (2 additional authors not shown)

Abstract: We present new maps of the Milky Way disk showing the distribution of metallicity ([Fe/H]), $α$-element abundances ([Mg/Fe]), and stellar age, using a sample of 66,496 red giant stars from the final data release (DR17) of the Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey. We measure radial and vertical gradients, quantify the distribution functions for age and metallicity,… ▽ More We present new maps of the Milky Way disk showing the distribution of metallicity ([Fe/H]), $α$-element abundances ([Mg/Fe]), and stellar age, using a sample of 66,496 red giant stars from the final data release (DR17) of the Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey. We measure radial and vertical gradients, quantify the distribution functions for age and metallicity, and explore chemical clock relations across the Milky Way for the low-$α$ disk, high-$α$ disk, and total population independently. The low-$α$ disk exhibits a negative radial metallicity gradient of $-0.06 \pm 0.001$ dex kpc$^{-1}$, which flattens with distance from the midplane. The high-$α$ disk shows a flat radial gradient in metallicity and age across nearly all locations of the disk. The age and metallicity distribution functions shift from negatively skewed in the inner Galaxy to positively skewed at large radius. Significant bimodality in the [Mg/Fe]-[Fe/H] plane and in the [Mg/Fe]-age relation persist across the entire disk. The age estimates have typical uncertainties of $\sim0.15$ in $\log$(age) and may be subject to additional systematic errors, which impose limitations on conclusions drawn from this sample. Nevertheless, these results act as critical constraints on galactic evolution models, constraining which physical processes played a dominant role in the formation of the Milky Way disk. We discuss how radial migration predicts many of the observed trends near the solar neighborhood and in the outer disk, but an additional more dramatic evolution history, such as the multi-infall model or a merger event, is needed to explain the chemical and age bimodality elsewhere in the Galaxy. △ Less

Submitted 25 July, 2023; originally announced July 2023.

Comments: 41 pages, 32 figures, accepted to ApJ

Journal ref: ApJ 954 124 (2023)

KPM: A Flexible and Data-Driven K-Process Model for Nucleosynthesis

Authors: Emily J. Griffith, David W. Hogg, Julianne J. Dalcanton, Sten Hasselquist, Bridget Ratcliffe, Melissa Ness, David H. Weinberg

Abstract: The element abundance pattern found in Milky Way disk stars is close to two-dimensional, dominated by production from one prompt process and one delayed process. This simplicity is remarkable, since the elements are produced by a multitude of nucleosynthesis mechanisms operating in stars with a wide range of progenitor masses. We fit the abundances of 14 elements for 48,659 red-giant stars from AP… ▽ More The element abundance pattern found in Milky Way disk stars is close to two-dimensional, dominated by production from one prompt process and one delayed process. This simplicity is remarkable, since the elements are produced by a multitude of nucleosynthesis mechanisms operating in stars with a wide range of progenitor masses. We fit the abundances of 14 elements for 48,659 red-giant stars from APOGEE DR17 using a flexible, data-driven K-process model -- dubbed KPM. In our fiducial model, with $K=2$, each abundance in each star is described as the sum of a prompt and a delayed process contribution. We find that KPM with $K=2$ is able to explain the abundances well, recover the observed abundance bimodality, and detect the bimodality over a greater range in metallicity than previously has been possible. We compare to prior work by Weinberg et al. (2022), finding that KPM produces similar results, but that KPM better predicts stellar abundances, especially for elements C+N and Mn and for stars at super-solar metallicities. The model fixes the relative contribution of the prompt and delayed process to two elements to break degeneracies and improve interpretability; we find that some of the nucleosynthetic implications are dependent upon these detailed choices. We find that moving to four processes adds flexibility and improves the model's ability to predict the stellar abundances, but doesn't qualitatively change the story. The results of KPM will help us to interpret and constrain the formation of the Galaxy disk, the relationship between abundances and ages, and the physics of nucleosynthesis. △ Less

Submitted 8 December, 2023; v1 submitted 11 July, 2023; originally announced July 2023.

Comments: 33 pages, 13 figures, 3 tables

New Dimensions of Galactic Chemical Evolution

Authors: David H. Weinberg

Abstract: Dramatic recent progress in understanding galactic chemical evolution (GCE) has been driven partly by direct observations of the distant past with HST and JWST of stellar abundances from giant high-resolution spectroscopic surveys (APOGEE, GALAH) and the complementary power of Gaia astrometry and photometry. Focusing on archaeology, I give a rapid-fire, and I hope synthesizing, review of work my c… ▽ More Dramatic recent progress in understanding galactic chemical evolution (GCE) has been driven partly by direct observations of the distant past with HST and JWST of stellar abundances from giant high-resolution spectroscopic surveys (APOGEE, GALAH) and the complementary power of Gaia astrometry and photometry. Focusing on archaeology, I give a rapid-fire, and I hope synthesizing, review of work my collaborators and I have done on theoretical modeling and observational interpretation. I discuss (1) the interleaved but distinguishable roles of stellar scale astrophysics and galactic scale astrophysics in governing GCE, (2) the use of abundance ratio trends to empirically infer nucleosynthetic yields, (3) the uncertainty in the overall scale of yields and its degeneracy with the importance of galactic outflows, (4) the emergence of equilibrium in GCE, (5) the dimensionality of the stellar distribution in chemical abundance space, and (6) insights from chemical abundances on the early history of the Milky Way, including measurements of the intrinsic scatter of abundance ratios in metal-poor stars (-2 < [Fe/H] < -1), suggesting that a typical halo star at this metallicity is enriched by the products of $N \sim 50$ supernovae mixed over $\sim 10^5 M_\odot$ of star-forming gas. △ Less

Submitted 15 June, 2023; originally announced June 2023.

Comments: 10 pp. running lightly over a lot of ground. Proceedings of IAU Symposium 377, Early Disk-Galaxy Formation: From JWST to the Milky Way, eds. F. Tabatabaei, B. Barbuy, and Y. Ting. Reactions welcome!

The Brain Tumor Segmentation (BraTS-METS) Challenge 2023: Brain Metastasis Segmentation on Pre-treatment MRI

Authors: Ahmed W. Moawad, Anastasia Janas, Ujjwal Baid, Divya Ramakrishnan, Rachit Saluja, Nader Ashraf, Leon Jekel, Raisa Amiruddin, Maruf Adewole, Jake Albrecht, Udunna Anazodo, Sanjay Aneja, Syed Muhammad Anwar, Timothy Bergquist, Evan Calabrese, Veronica Chiang, Verena Chung, Gian Marco Marco Conte, Farouk Dako, James Eddy, Ivan Ezhov, Ariana Familiar, Keyvan Farahani, Juan Eugenio Iglesias, Zhifan Jiang , et al. (206 additional authors not shown)

Abstract: The translation of AI-generated brain metastases (BM) segmentation into clinical practice relies heavily on diverse, high-quality annotated medical imaging datasets. The BraTS-METS 2023 challenge has gained momentum for testing and benchmarking algorithms using rigorously annotated internationally compiled real-world datasets. This study presents the results of the segmentation challenge and chara… ▽ More The translation of AI-generated brain metastases (BM) segmentation into clinical practice relies heavily on diverse, high-quality annotated medical imaging datasets. The BraTS-METS 2023 challenge has gained momentum for testing and benchmarking algorithms using rigorously annotated internationally compiled real-world datasets. This study presents the results of the segmentation challenge and characterizes the challenging cases that impacted the performance of the winning algorithms. Untreated brain metastases on standard anatomic MRI sequences (T1, T2, FLAIR, T1PG) from eight contributed international datasets were annotated in stepwise method: published UNET algorithms, student, neuroradiologist, final approver neuroradiologist. Segmentations were ranked based on lesion-wise Dice and Hausdorff distance (HD95) scores. False positives (FP) and false negatives (FN) were rigorously penalized, receiving a score of 0 for Dice and a fixed penalty of 374 for HD95. Eight datasets comprising 1303 studies were annotated, with 402 studies (3076 lesions) released on Synapse as publicly available datasets to challenge competitors. Additionally, 31 studies (139 lesions) were held out for validation, and 59 studies (218 lesions) were used for testing. Segmentation accuracy was measured as rank across subjects, with the winning team achieving a LesionWise mean score of 7.9. Common errors among the leading teams included false negatives for small lesions and misregistration of masks in space.The BraTS-METS 2023 challenge successfully curated well-annotated, diverse datasets and identified common errors, facilitating the translation of BM segmentation across varied clinical environments and providing personalized volumetric reports to patients undergoing BM treatment. △ Less

Submitted 17 June, 2024; v1 submitted 1 June, 2023; originally announced June 2023.

Measuring the dynamical length of galactic bars

Authors: Michael S. Petersen, Martin D. Weinberg, Neal Katz

Abstract: We define a physically-motivated measure for galactic bar length, called the dynamical length. The dynamical length of the bar corresponds to the radial extent of the orbits that are the backbone supporting the bar feature. We propose a direct observational technique using integral field unit spectroscopy to measure it. Identifying these orbits and using the dynamical length is a more faithful tra… ▽ More We define a physically-motivated measure for galactic bar length, called the dynamical length. The dynamical length of the bar corresponds to the radial extent of the orbits that are the backbone supporting the bar feature. We propose a direct observational technique using integral field unit spectroscopy to measure it. Identifying these orbits and using the dynamical length is a more faithful tracer of the secular evolution and influence of the bar. We demonstrate the success of the metric for recovering the maximal bar-parenting orbit in a range of simulations, and to show its promise we perform its measurement on a real galaxy. We also study the difference between traditionally used ellipse fit approaches to determine bar length and the dynamical length proposed here in a wide range of bar-forming N-body simulations of a stellar disc and dark matter halo. We find that ellipse fitting may severely overestimate measurements of the bar length by a factor of 1.5-2.5 relative to the extent of the orbits that are trapped and actually comprise the bar. This bias leads to overestimates of both bar mass and the ratio of corotation radius to bar length, i.e. the bar speed, affecting inferences about the evolution of bars in the real universe. △ Less

Submitted 22 May, 2023; originally announced May 2023.

Comments: Submitted to MNRAS, comments welcome. arXiv admin note: text overlap with arXiv:1903.08203

On the Use of Singular Value Decomposition as a Clutter Filter for Ultrasound Flow Imaging

Authors: Kai Riemer, Marcelo Lerendegui, Matthieu Toulemonde, Jiaqi Zhu, Christopher Dunsby, Peter D. Weinberg, Meng-Xing Tang

Abstract: Filtering based on Singular Value Decomposition (SVD) provides substantial separation of clutter, flow and noise in high frame rate ultrasound flow imaging. The use of SVD as a clutter filter has greatly improved techniques such as vector flow imaging, functional ultrasound and super-resolution ultrasound localization microscopy. The removal of clutter and noise relies on the assumption that tissu… ▽ More Filtering based on Singular Value Decomposition (SVD) provides substantial separation of clutter, flow and noise in high frame rate ultrasound flow imaging. The use of SVD as a clutter filter has greatly improved techniques such as vector flow imaging, functional ultrasound and super-resolution ultrasound localization microscopy. The removal of clutter and noise relies on the assumption that tissue, flow and noise are each represented by different subsets of singular values, so that their signals are uncorrelated and lay on orthogonal sub-spaces. This assumption fails in the presence of tissue motion, for near-wall or microvascular flow, and can be influenced by an incorrect choice of singular value thresholds. Consequently, separation of flow, clutter and noise is imperfect, which can lead to image artefacts not present in the original data. Temporal and spatial fluctuation in intensity are the commonest artefacts, which vary in appearance and strengths. Ghosting and splitting artefacts are observed in the microvasculature where the flow signal is sparsely distributed. Singular value threshold selection, tissue motion, frame rate, flow signal amplitude and acquisition length affect the prevalence of these artefacts. Understanding what causes artefacts due to SVD clutter and noise removal is necessary for their interpretation. △ Less

Submitted 25 April, 2023; originally announced April 2023.

Comments: 10 pages, 7 figures

Buzzard to Cardinal: Improved Mock Catalogs for Large Galaxy Surveys

Authors: Chun-Hao To, Joseph DeRose, Risa H. Wechsler, Eli Rykoff, Hao-Yi Wu, Susmita Adhikari, Elisabeth Krause, Eduardo Rozo, David H. Weinberg

Abstract: We present the Cardinal mock galaxy catalogs, a new version of the Buzzard simulation that has been updated to support ongoing and future cosmological surveys, including DES, DESI, and LSST. These catalogs are based on a one-quarter sky simulation populated with galaxies out to a redshift of $z=2.35$ to a depth of $m_{\rm{r}}=27$. Compared to the Buzzard mocks, the Cardinal mocks include an update… ▽ More We present the Cardinal mock galaxy catalogs, a new version of the Buzzard simulation that has been updated to support ongoing and future cosmological surveys, including DES, DESI, and LSST. These catalogs are based on a one-quarter sky simulation populated with galaxies out to a redshift of $z=2.35$ to a depth of $m_{\rm{r}}=27$. Compared to the Buzzard mocks, the Cardinal mocks include an updated subhalo abundance matching (SHAM) model that considers orphan galaxies and includes mass-dependent scatter between galaxy luminosity and halo properties. This model can simultaneously fit galaxy clustering and group--galaxy cross-correlations measured in three different luminosity threshold samples. The Cardinal mocks also feature a new color assignment model that can simultaneously fit color-dependent galaxy clustering in three different luminosity bins. We have developed an algorithm that uses photometric data to improve the color assignment model further and have also developed a novel method to improve small-scale lensing below the ray-tracing resolution. These improvements enable the Cardinal mocks to accurately reproduce the abundance of galaxy clusters and the properties of lens galaxies in the Dark Energy Survey data. As such, these simulations will be a valuable tool for future cosmological analyses based on large sky surveys. The cardinal mock will be released upon publication at https://chunhaoto.com/cardinalsim. △ Less

Submitted 21 March, 2023; originally announced March 2023.

Comments: 37 pages; 20 figures. See figure 1 for summary and appendix K for a list of main differences between Cardinal and Buzzard. To be submitted to APJ; Comments welcome. Highlights can be found at https://chunhaoto.com/cardinalsim

Investigating biomechanical determinants of endothelial permeability in a modified hollow fibre bioreactor

Authors: Stephen G Gray, Peter D Weinberg

Abstract: Effects of mechanical stress on the permeability of vascular endothelium are important to normal physiology and may be critical in the development of atherosclerosis, where they can account for the patchy arterial distribution of the disease. Such properties are frequently investigated in vitro. Here we evaluate and use the hollow fibre bioreactor for this purpose; in this system, endothelial cell… ▽ More Effects of mechanical stress on the permeability of vascular endothelium are important to normal physiology and may be critical in the development of atherosclerosis, where they can account for the patchy arterial distribution of the disease. Such properties are frequently investigated in vitro. Here we evaluate and use the hollow fibre bioreactor for this purpose; in this system, endothelial cells form a confluent monolayer lining numerous plastic capillaries with porous walls, contained in a cartridge. The capillaries were perfused with a near-aortic waveform by an external pump, and permeability was assessed by the movement of rhodamine-labelled albumin from the intracapillary space to the extracapillary space. Confluence and quiescence of the cells was confirmed by electron microscopy and measurements of glucose consumption and permeability. The system was able to detect previously established influences on permeability: tracer transport was increased by acute application of shear stress and decreased by chronic shear stress compared to a static control, and was increased by thrombin or an NO synthase inhibitor under chronic shear. Increasing viscosity by addition of xanthan gum reduced permeability under both acute and chronic shear. Addition of damping chambers to reduce flow pulsatility increased permeability. Modifying the cartridge to allow chronic convection across the monolayer increased effective permeability more than could be explained the addition of convective transport alone, indicating that it caused an increase in permeability. The off-the-shelf hollow fibre bioreactor provides an excellent system for investigating the biomechanics of endothelial permeability and its potential is increased by simple modifications. △ Less

Submitted 2 March, 2023; originally announced March 2023.

Comments: 33 pages, 13 figures

Communication Trade-offs in Federated Learning of Spiking Neural Networks

Authors: Soumi Chaki, David Weinberg, Ayca Özcelikkale

Abstract: Spiking Neural Networks (SNNs) are biologically inspired alternatives to conventional Artificial Neural Networks (ANNs). Despite promising preliminary results, the trade-offs in the training of SNNs in a distributed scheme are not well understood. Here, we consider SNNs in a federated learning setting where a high-quality global model is created by aggregating multiple local models from the client… ▽ More Spiking Neural Networks (SNNs) are biologically inspired alternatives to conventional Artificial Neural Networks (ANNs). Despite promising preliminary results, the trade-offs in the training of SNNs in a distributed scheme are not well understood. Here, we consider SNNs in a federated learning setting where a high-quality global model is created by aggregating multiple local models from the clients without sharing any data. We investigate federated learning for training multiple SNNs at clients when two mechanisms reduce the uplink communication cost: i) random masking of the model updates sent from the clients to the server; and ii) client dropouts where some clients do not send their updates to the server. We evaluated the performance of the SNNs using a subset of the Spiking Heidelberg digits (SHD) dataset. The results show that a trade-off between the random masking and the client drop probabilities is crucial to obtain a satisfactory performance for a fixed number of clients. △ Less

Submitted 27 February, 2023; originally announced March 2023.

Comments: 5 pages, 7 figures, 1 algorithm, 1 table

Size-dependent lattice symmetry breaking determines the exciton fine structure of perovskite nanocrystals

Authors: Daniel Weinberg, Yoonjae Park, David T. Limmer, Eran Rabani

Abstract: The ordering of optically bright and dark excitonic states in lead-halide perovskite nanocrystals has been a matter of some debate. It has been proposed that the unusually short radiative lifetimes in these materials is due to an optically bright excitonic ground state, a unique situation among all nanomaterials. This proposal was based on the influence of the Rashba effect driven by lattice-induc… ▽ More The ordering of optically bright and dark excitonic states in lead-halide perovskite nanocrystals has been a matter of some debate. It has been proposed that the unusually short radiative lifetimes in these materials is due to an optically bright excitonic ground state, a unique situation among all nanomaterials. This proposal was based on the influence of the Rashba effect driven by lattice-induced inversion symmetry breaking. Direct measurement of the excitonic emission under magnetic fields has shown the signature of a dark ground state, bringing the role of the Rashba effect into question. Here, we use a fully atomistic theory to model the exciton fine structure of perovskite nanocrystals accounting for the realistic lattice distortion at the nanoscale. We calculate optical gaps and exciton fine structure that compare favorably with a wide range of experimental works. We find a non-monotonic dependence of the exciton fine structure splittings due to a size dependence structural transition between cubic and orthorhombic phases. In addition, the excitonic ground state is found to be dark with nearly pure spin triplet character resulting from a small Rashba coupling. We additionally explore the intertwined effects of lattice distortion and nanocrystal shape on the fine structure splittings, clarifying observations on poly-disperse nanocrystals. △ Less

Submitted 1 March, 2023; originally announced March 2023.

Comments: 17 pages, 11 figures

A Comprehensive Perturbative Formalism for Phase Mixing in Perturbed Disks. II. Phase Spirals in an Inhomogeneous Disk Galaxy with a Non-responsive Dark Matter Halo

Authors: Uddipan Banik, Frank C. van den Bosch, Martin D. Weinberg

Abstract: We develop a linear perturbative formalism to compute the response of an inhomogeneous stellar disk embedded in a non-responsive dark matter halo to perturbations like bars, spiral arms and satellite galaxy encounters. Without self-gravity to reinforce it, the response of a Fourier mode phase mixes away due to an intrinsic spread in the vertical ($Ω_z$), radial ($Ω_r$) and azimuthal ($Ω_φ$) freque… ▽ More We develop a linear perturbative formalism to compute the response of an inhomogeneous stellar disk embedded in a non-responsive dark matter halo to perturbations like bars, spiral arms and satellite galaxy encounters. Without self-gravity to reinforce it, the response of a Fourier mode phase mixes away due to an intrinsic spread in the vertical ($Ω_z$), radial ($Ω_r$) and azimuthal ($Ω_φ$) frequencies, giving rise to local phase-space spirals. Collisional diffusion due to scattering of stars by structures like giant molecular clouds causes super-exponential damping of the phase-spiral amplitude. The $z-v_z$ phase-spiral is 1-armed (2-armed) for vertically anti-symmetric (symmetric) bending (breathing) modes. Only transient perturbations with timescales ($τ_{\mathrm{P}}$) comparable to the vertical oscillation period ($τ_z \sim 1/Ω_z$) trigger $z-v_z$ phase-spirals. Each $(n,l,m)$ mode of the response to impulsive ($τ_{\mathrm{P}}<τ=1/(nΩ_z+lΩ_r+mΩ_φ)$) perturbations is power law ($\sim τ_{\mathrm{P}}/τ$) suppressed, but that to adiabatic ($τ_{\mathrm{P}}>τ$) perturbations is exponentially weak ($\sim \exp{\left[-\left(τ_{\mathrm{P}}/τ\right)^α\right]}$) except resonant ($τ\to \infty$) modes. Slower ($τ_{\mathrm{P}}>τ_z$) perturbations, e.g., distant encounters with satellite galaxies, induce stronger bending modes. If the Gaia phase-spiral was triggered by a satellite, Sagittarius is the leading contender as it dominates the Solar neighborhood response of the Milky Way disk to satellite encounters. However, survival against collisional damping necessitates that the impact occurred within $\sim 0.6-0.7$ Gyr ago. We discuss how the detailed galactic potential dictates the phase-spiral shape: phase mixing occurs slower and phase-spirals are less wound in the outer disk and in presence of an ambient halo. △ Less

Submitted 28 February, 2023; originally announced March 2023.

Comments: 10 figures, 1 table; Submitted to ApJ

A framework to measure the properties of intergalactic metal systems with two-point flux statistics

Authors: Naim Göksel Karaçaylı, Paul Martini, David H. Weinberg, Vid Iršič, J. Aguilar, S. Ahlen, D. Brooks, A. de la Macorra, A. Font-Ribera, S. Gontcho A Gontcho, J. Guy, T. Kisner, R. Miquel, C. Poppett, C. Ravoux, M. Schubnell, G. Tarlé, B. A. Weaver, Z. Zhou

Abstract: The abundance, temperature, and clustering of metals in the intergalactic medium are important parameters for understanding their cosmic evolution and quantifying their impact on cosmological analysis with the Ly $α$ forest. The properties of these systems are typically measured from individual quasar spectra redward of the quasar's Ly $α$ emission line, yet that approach may provide biased result… ▽ More The abundance, temperature, and clustering of metals in the intergalactic medium are important parameters for understanding their cosmic evolution and quantifying their impact on cosmological analysis with the Ly $α$ forest. The properties of these systems are typically measured from individual quasar spectra redward of the quasar's Ly $α$ emission line, yet that approach may provide biased results due to selection effects. We present an alternative approach to measure these properties in an unbiased manner with the two-point statistics commonly employed to quantify large-scale structure. Our model treats the observed flux of a large sample of quasar spectra as a continuous field and describes the one-dimensional, two-point statistics of this field with three parameters per ion: the abundance (column density distribution), temperature (Doppler parameter) and clustering (cloud-cloud correlation function). We demonstrate this approach on multiple ions (e.g., C IV, Si IV, Mg II) with early data from the Dark Energy Spectroscopic Instrument (DESI) and high-resolution spectra from the literature. Our initial results show some evidence that the C IV abundance is higher than previous measurements and evidence for abundance evolution over time. The first full year of DESI observations will have over an order of magnitude more quasar spectra than this study. In a future paper we will use those data to measure the growth of clustering and its impact on the Ly $α$ forest, as well as test other DESI analysis infrastructure such as the pipeline noise estimates and the resolution matrix. △ Less

Submitted 6 May, 2023; v1 submitted 14 February, 2023; originally announced February 2023.

Comments: 15 pages, 14 figures, accepted to MNRAS

The Non-Axisymmetric Influence: Radius and Angle-Dependent Trends in a Barred Galaxy

Authors: Carrie Filion, Rachel L. McClure, Martin D. Weinberg, Elena D'Onghia, Kathryne J. Daniel

Abstract: Many disc galaxies host galactic bars, which exert time-dependent, non-axisymmetric forces that can alter the orbits of stars. There should be both angle and radius-dependence in the resulting radial rearrangement of stars ('radial mixing') due to a bar; we present here novel results and trends through analysis of the joint impact of these factors. We use an N-body simulation to investigate the ch… ▽ More Many disc galaxies host galactic bars, which exert time-dependent, non-axisymmetric forces that can alter the orbits of stars. There should be both angle and radius-dependence in the resulting radial rearrangement of stars ('radial mixing') due to a bar; we present here novel results and trends through analysis of the joint impact of these factors. We use an N-body simulation to investigate the changes in the radial locations of star particles in a disc after a bar forms by quantifying the change in orbital radii in a series of annuli at different times post bar-formation. We find that the bar induces both azimuth angle- and radius- dependent trends in the median distance that stars have travelled to enter a given annulus. Angle-dependent trends are present at all radii we consider, and the radius-dependent trends roughly divide the disc into three 'zones'. In the inner zone, stars generally originated at larger radii and their orbits evolved inwards. Stars in the outer zone likely originated at smaller radii and their orbits evolved outwards. In the intermediate zone, there is no net inwards or outwards evolution of orbits. We adopt a simple toy model of a radius-dependent initial metallicity gradient and discuss recent observational evidence for angle-dependent stellar metallicity variations in the Milky Way in the context of this model. We briefly comment on the possibility of using observed angle-dependent metallicity trends to learn about the initial metallicity gradient(s) and the radial rearrangement that occurred in the disc. △ Less

Submitted 3 November, 2023; v1 submitted 2 February, 2023; originally announced February 2023.

Comments: Published in MNRAS, updated to match final refereed version in which a discussion on external galaxies was added

Journal ref: MNRAS, Volume 524, 2023, pp.276-287

Dynamical Data Mining Captures Disc-Halo Couplings that Structure Galaxies

Authors: Alexander Johnson, Michael S. Petersen, Kathryn V. Johnston, Martin D. Weinberg

Abstract: Studying coupling between different galactic components is a challenging problem in galactic dynamics. Using basis function expansions (BFEs) and multichannel singular spectrum analysis (mSSA) as a means of dynamical data mining, we discover evidence for two multi-component disc-halo dipole modes in a Milky-Way-like simulated galaxy. One of the modes grows throughout the simulation, while the othe… ▽ More Studying coupling between different galactic components is a challenging problem in galactic dynamics. Using basis function expansions (BFEs) and multichannel singular spectrum analysis (mSSA) as a means of dynamical data mining, we discover evidence for two multi-component disc-halo dipole modes in a Milky-Way-like simulated galaxy. One of the modes grows throughout the simulation, while the other decays throughout the simulation. The multi-component disc-halo modes are driven primarily by the halo, and have implications for the structural evolution of galaxies, including observations of lopsidedness and other non-axisymmetric structure. In our simulation, the modes create surface density features up to 10 per cent relative to the equilibrium model stellar disc. While the simulated galaxy was constructed to be in equilibrium, BFE+mSSA also uncovered evidence of persistent periodic signals incited by aphysical initial conditions disequilibrium, including rings and weak two-armed spirals, both at the 1 per cent level. The method is sensitive to distinct evolutionary features at and even below the 1 per cent level of surface density variation. The use of mSSA produced clean signals for both modes and disequilibrium, efficiently removing variance owing to estimator noise from the input BFE time series. The discovery of multi-component halo-disc modes is strong motivation for application of BFE+mSSA to the rich zoo of dynamics of multi-component interacting galaxies. △ Less

Submitted 5 January, 2023; originally announced January 2023.

Comments: Submitted to MNRAS, comments welcome

Strong Outflows and Inefficient Star Formation in the Reionization-era Ultra-faint Dwarf Galaxy Eridanus II

Authors: Nathan R. Sandford, David H. Weinberg, Daniel R. Weisz, Sal Wanying Fu

Abstract: We present novel constraints on the underlying galaxy formation physics (e.g., mass loading factor, star formation history, metal retention) at $z\gtrsim7$ for the low-mass ($M_*\sim10^5$ M$_\odot$) Local Group ultra-faint dwarf galaxy (UFD) Eridanus {\sc II} (Eri II). Using a hierarchical Bayesian framework, we apply a one-zone chemical evolution model to Eri II's CaHK-based photometric metallici… ▽ More We present novel constraints on the underlying galaxy formation physics (e.g., mass loading factor, star formation history, metal retention) at $z\gtrsim7$ for the low-mass ($M_*\sim10^5$ M$_\odot$) Local Group ultra-faint dwarf galaxy (UFD) Eridanus {\sc II} (Eri II). Using a hierarchical Bayesian framework, we apply a one-zone chemical evolution model to Eri II's CaHK-based photometric metallicity distribution function (MDF; [Fe/H]) and find that the evolution of Eri II is well-characterized by a short, exponentially declining star-formation history ($τ_\text{SFH}=0.39\pm_{0.13}^{0.18}$ Gyr), a low star-formation efficiency ($τ_\text{SFE}=27.56\pm_{12.92}^{25.14}$ Gyr), and a large mass-loading factor ($η=194.53\pm_{42.67}^{33.37}$). Our results are consistent with Eri II forming the majority of its stars before the end of reionization. The large mass-loading factor implies strong outflows in the early history of Eri II and is in good agreement with theoretical predictions for the mass-scaling of galactic winds. It also results in the ejection of $>$90\% of the metals produced in Eri II. We make predictions for the distribution of [Mg/Fe]-[Fe/H] in Eri II as well as the prevalence of ultra metal-poor stars, both of which can be tested by future chemical abundance measurements. Spectroscopic follow-up of the highest metallicity stars in Eri II ($\text{[Fe/H]} > -2$) will greatly improve model constraints. Our new framework can readily be applied to all UFDs throughout the Local Group, providing new insights into the underlying physics governing the evolution of the faintest galaxies in the reionization era. △ Less

Submitted 31 August, 2023; v1 submitted 30 October, 2022; originally announced October 2022.

Comments: 20 pages; 12 figures, submitted to MNRAS

Untangling the Sources of Abundance Dispersion in Low-Metallicity Stars

Authors: Emily J. Griffith, Jennifer A. Johnson, David H. Weinberg, Ilya Ilyin, James W. Johnson, Romy Rodriguez-Martinez, Klaus G. Strassmeier

Abstract: We measure abundances of 12 elements (Na, Mg, Si, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni) in a sample of 86 metal-poor ($-2 \lesssim \text{[Fe/H]} \lesssim -1$) subgiant stars in the solar neighborhood. Abundances are derived from high-resolution spectra taken with the Potsdam Echelle Polarimetric and Spectroscopic Instrument on the Large Binocular Telescope, modeled using iSpec and MOOG. By carefully… ▽ More We measure abundances of 12 elements (Na, Mg, Si, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni) in a sample of 86 metal-poor ($-2 \lesssim \text{[Fe/H]} \lesssim -1$) subgiant stars in the solar neighborhood. Abundances are derived from high-resolution spectra taken with the Potsdam Echelle Polarimetric and Spectroscopic Instrument on the Large Binocular Telescope, modeled using iSpec and MOOG. By carefully quantifying the impact of photon-noise ($< 0.05$ dex for all elements) we robustly measure the intrinsic scatter of abundance ratios. At fixed [Fe/H] the RMS intrinsic scatter in [X/Fe] ranges from 0.04 dex (Cr) to 0.16 dex (Na), with a median of 0.08 dex. Scatter in [X/Mg] is similar, and accounting for [$α$/Fe] only reduces the overall scatter moderately. We consider several possible origins of the intrinsic scatter with particular attention to fluctuations in the relative enrichment by core-collapse supernovae (CCSN) and Type Ia supernovae (SNIa) and stochastic sampling of the CCSN progenitor mass distribution. The stochastic sampling scenario provides a good quantitative explanation of our data if the effective number of CCSN contributing to the enrichment of a typical sample star is $N \sim 50$. At the median metallicity of our sample, this interpretation implies that the CCSN ejecta are mixed over a gas mass $\sim 10^5 M_{\odot}$ before forming stars. The scatter of elemental abundance ratios is a powerful diagnostic test for simulations of star formation, feedback, and gas mixing in the early phases of the Galaxy. △ Less

Submitted 14 October, 2022; v1 submitted 4 October, 2022; originally announced October 2022.

Comments: 28 pages, 15 figures, 5 tables

New dipole instabilities in spherical stellar systems

Authors: Martin D. Weinberg

Abstract: Spherical stellar systems have weakly-damped response modes. The dipole modes are seiche modes. The quadrupole are zero pattern-speed prolate modes, the stable precursors to the radial-orbit instability (ROI). We demonstrate that small wiggles in the distribution function (DF) can destabilise the dipole modes and describe the newly identified instabilities in NFW-like dark-matter (DM) halos and ot… ▽ More Spherical stellar systems have weakly-damped response modes. The dipole modes are seiche modes. The quadrupole are zero pattern-speed prolate modes, the stable precursors to the radial-orbit instability (ROI). We demonstrate that small wiggles in the distribution function (DF) can destabilise the dipole modes and describe the newly identified instabilities in NFW-like dark-matter (DM) halos and other power-law spherical systems. The modes were identified in N-body simulations using multivariate singular spectrum analysis (MSSA) and corroborated using linear-response theory. The new mode peaks inside the half-mass radius but has a pattern speed typical of an outer-halo orbit. As it grows, the radial angle of the eccentric orbits that make up the mode correlate and lose angular momentum by a resonant couple to outer-halo orbits. This leads to an unsteady pattern with a density enhancement that swings from one side of the halo to another along a diameter, like the orbits that comprise the instability. In this way, the dipole mode is similar to the ROI. Since the DF found in Nature is unlikely to be smooth and isotropic with $df(E)/dE<0$ necessary for Antonov stability, these modes may be ubiquitous albeit slowly growing. Halos that are less extended than NFW, such as the Hernquist model, tend to be stable to this dipole instability. We present the critical stability exponents for one- and two-power models. These different critical outer power-law exponents illustrate that the gravitational coupling between the inner and outer DM halo depends on the global shape of density profile. △ Less

Submitted 24 August, 2023; v1 submitted 14 September, 2022; originally announced September 2022.

Comments: 16 pages, 19 figures

The Poor Old Heart of the Milky Way

Authors: Hans-Walter Rix, Vedant Chandra, René Andrae, Adrian M. Price-Whelan, David H. Weinberg, Charlie Conroy, Morgan Fouesneau, David W. Hogg, Francesca De Angeli, Rohan P. Naidu, Maosheng Xiang, Daniela Ruz-Mieres

Abstract: Massive disk galaxies like our Milky Way should host an ancient, metal-poor, and centrally concentrated stellar population. This population reflects the star formation and enrichment in the few most massive progenitor components that coalesced at high redshift to form the proto-Galaxy. While metal-poor stars are known to reside in the inner few kiloparsecs of our Galaxy, current data do not yet pr… ▽ More Massive disk galaxies like our Milky Way should host an ancient, metal-poor, and centrally concentrated stellar population. This population reflects the star formation and enrichment in the few most massive progenitor components that coalesced at high redshift to form the proto-Galaxy. While metal-poor stars are known to reside in the inner few kiloparsecs of our Galaxy, current data do not yet provide a comprehensive picture of such a metal-poor "heart" of the Milky Way. We use information from Gaia DR3, especially the XP spectra, to construct a sample of 2 million bright (BP $<15.5$ mag) giant stars within $30^\circ$ of the Galactic Center with robust [M/H] estimates, $δ$ [M/H] $\lesssim 0.1$. For most sample members we can calculate orbits based on Gaia RVS velocities and astrometry. This sample reveals an extensive, ancient, and metal-poor population that includes $\sim 18,000$ stars with $-2.7<$ [M/H] $<-1.5$, representing a stellar mass of $\gtrsim 5\times 10^7$ M$_\odot$. The spatial distribution of these [M/H] $<-1.5$ stars has a Gaussian extent of only $σ_{\mathrm{R_{GC}}} \sim 2.7$ kpc around the Galactic center, with most of these orbits being confined to the inner Galaxy. At high orbital eccentricities, there is clear evidence for accreted halo stars in their pericentral orbit phase. Stars with [M/H] $< -2$ show no net rotation, whereas those with [M/H] $\sim -1$ are rotation dominated. Most of the tightly bound stars show $[α/\text{Fe}]$-enhancement and [Al/Fe]-[Mn/Fe] abundance patterns expected for an origin in the more massive portions of the proto-Galaxy. These central, metal-poor stars most likely predate the oldest part of the disk ($τ_{\text{age}}\approx 12.5$ Gyrs), which implies that they formed at $z\gtrsim 5$, forging the proto-Milky Way. △ Less

Submitted 6 September, 2022; originally announced September 2022.

Comments: 21 pages, 9 figures, submitted to ApJ

3D Super-Resolution Ultrasound with Adaptive Weight-Based Beamforming

Authors: Jipeng Yan, Bingxue Wang, Kai Riemer, Joseph Hansen-Shearer, Marcelo Lerendegui, Matthieu Toulemonde, Christopher J Rowlands, Peter D. Weinberg, Meng-Xing Tang

Abstract: Super-resolution ultrasound (SRUS) imaging through localising and tracking sparse microbubbles has been shown to reveal microvascular structure and flow beyond the wave diffraction limit. Most SRUS studies use standard delay and sum (DAS) beamforming, where large main lobe and significant side lobes make separation and localisation of densely distributed bubbles challenging, particularly in 3D due… ▽ More Super-resolution ultrasound (SRUS) imaging through localising and tracking sparse microbubbles has been shown to reveal microvascular structure and flow beyond the wave diffraction limit. Most SRUS studies use standard delay and sum (DAS) beamforming, where large main lobe and significant side lobes make separation and localisation of densely distributed bubbles challenging, particularly in 3D due to the typically small aperture of matrix array probes. This study aims to improve 3D SRUS by implementing a low-cost 3D coherence beamformer based on channel signal variance, as well as two other adaptive weight-based coherence beamformers: nonlinear beamforming with p-th root compression and coherence factor. The 3D coherence beamformers, together with DAS, are compared in computer simulation, on a microflow phantom, and in vivo. Simulation results demonstrate that the adaptive weight-based beamformers can significantly narrow the main lobe and suppress the side lobes for modest computational cost. Significantly improved 3D SR images of microflow phantom and a rabbit kidney are obtained through the adaptive weight-based beamformers. The proposed variance-based beamformer performs best in simulations and experiments. △ Less

Submitted 25 August, 2022; originally announced August 2022.

Comments: Ultrasound localisation microscopy (ULM), super-resolution, contrast-enhanced ultrasound, 3D beamforming

Overview of the DESI Milky Way Survey

Authors: Andrew P. Cooper, Sergey E. Koposov, Carlos Allende Prieto, Christopher J. Manser, Namitha Kizhuprakkat, Adam D. Myers, Arjun Dey, Boris T. Gaensicke, Ting S. Li, Constance Rockosi, Monica Valluri, Joan Najita, Alis Deason, Anand Raichoor, Mei-Yu Wang, Yuan-Sen Ting, Bokyoung Kim, Andreia Carrillo, Wenting Wang, Leandro Beraldo e Silva, Jiwon Jesse Han, Jiani Ding, Miguel Sanchez-Conde, Jessica N. Aguilar, Steven Ahlen , et al. (40 additional authors not shown)

Abstract: We describe the Milky Way Survey (MWS) that will be undertaken with the Dark Energy Spectroscopic Instrument (DESI) on the Mayall 4m telescope at the Kitt Peak National Observatory. Over the next 5 yr DESI MWS will observe approximately seven million stars at Galactic latitudes |b|>20 degrees, with an inclusive target selection scheme focused on the thick disk and stellar halo. MWS will also inclu… ▽ More We describe the Milky Way Survey (MWS) that will be undertaken with the Dark Energy Spectroscopic Instrument (DESI) on the Mayall 4m telescope at the Kitt Peak National Observatory. Over the next 5 yr DESI MWS will observe approximately seven million stars at Galactic latitudes |b|>20 degrees, with an inclusive target selection scheme focused on the thick disk and stellar halo. MWS will also include several high-completeness samples of rare stellar types, including white dwarfs, low-mass stars within 100pc of the Sun, and horizontal branch stars. We summarize the potential of DESI to advance understanding of Galactic structure and stellar evolution. We introduce the final definitions of the main MWS target classes and estimate the number of stars in each class that will be observed. We describe our pipelines for deriving radial velocities, atmospheric parameters, and chemical abundances. We use ~500,000 spectra of unique stellar targets from the DESI Survey Validation program (SV) to demonstrate that our pipelines can measure radial velocities to ~1 km/s and [Fe/H] accurate to ~0.2 dex for typical stars in our main sample. We find the stellar parameter distributions from ~100 sq. deg of SV observations with >90% completeness on our main sample are in good agreement with expectations from mock catalogs and previous surveys. △ Less

Submitted 20 February, 2023; v1 submitted 17 August, 2022; originally announced August 2022.

Comments: Accepted for publication in ApJ, 44 pages, 25 figures, 4 tables, one of a suite of 8 papers detailing targeting for DESI; v2 added links to data shown in figures, added citations to other DESI papers, corrected author list and minor typos; v3 fixed minor errors in Fig. 6 and clarified associated text; v4 updated to include minor changes in response to review

DESI Bright Galaxy Survey: Final Target Selection, Design, and Validation

Authors: ChangHoon Hahn, Michael J. Wilson, Omar Ruiz-Macias, Shaun Cole, David H. Weinberg, John Moustakas, Anthony Kremin, Jeremy L. Tinker, Alex Smith, Risa H. Wechsler, Steven Ahlen, Shadab Alam, Stephen Bailey, David Brooks, Andrew P. Cooper, Tamara M. Davis, Kyle Dawson, Arjun Dey, Biprateep Dey, Sarah Eftekharzadeh, Daniel J. Eisenstein, Kevin Fanning, Jaime E. Forero-Romero, Carlos S. Frenk, Enrique Gaztañaga , et al. (35 additional authors not shown)

Abstract: Over the next five years, the Dark Energy Spectroscopic Instrument (DESI) will use 10 spectrographs with 5000 fibers on the 4m Mayall Telescope at Kitt Peak National Observatory to conduct the first Stage-IV dark energy galaxy survey. At $z < 0.6$, the DESI Bright Galaxy Survey (BGS) will produce the most detailed map of the Universe during the dark energy dominated epoch with redshifts of >10 mil… ▽ More Over the next five years, the Dark Energy Spectroscopic Instrument (DESI) will use 10 spectrographs with 5000 fibers on the 4m Mayall Telescope at Kitt Peak National Observatory to conduct the first Stage-IV dark energy galaxy survey. At $z < 0.6$, the DESI Bright Galaxy Survey (BGS) will produce the most detailed map of the Universe during the dark energy dominated epoch with redshifts of >10 million galaxies over 14,000 deg$^2$. In this work, we present and validate the final BGS target selection and survey design. From the Legacy Surveys, BGS will target a $r < 19.5$ magnitude-limited sample (BGS Bright); a fainter $19.5 < r < 20.175$ sample, color-selected to have high redshift efficiency (BGS Faint); and a smaller low-z quasar sample. BGS will observe these targets using exposure times, scaled to achieve uniform completeness, and visit each point on the footprint three times. We use observations from the Survey Validation programs conducted prior to the main survey along with realistic simulations to show that BGS can complete its strategy and make optimal use of `bright' time. We demonstrate that BGS targets have stellar contamination <1% and that their densities do not depend strongly on imaging properties. We also confirm that BGS Bright will achieve >80% fiber assignment efficiency. Finally, we show that BGS Bright and Faint will achieve >95% redshift success rates with no significant dependence on observing conditions. BGS meets the requirements for an extensive range of scientific applications. BGS will yield the most precise Baryon Acoustic Oscillations and Redshift-Space Distortions measurements at $z < 0.4$. It also presents opportunities to exploit new methods that require highly complete and dense galaxy samples (e.g. N-point statistics, multi-tracers). BGS further provides a powerful tool to study galaxy populations and the relations between galaxies and dark matter. △ Less

Submitted 17 August, 2022; originally announced August 2022.

Comments: AJ, submitted, 34 pages, 22 figures, one of a suite of 8 papers detailing targeting for DESI

A Comprehensive Perturbative Formalism for Phase-Mixing in Perturbed Disks. I. Phase spirals in an Infinite, Isothermal Slab

Authors: Uddipan Banik, Martin D. Weinberg, Frank C. van den Bosch

Abstract: Galactic disks are highly responsive systems that often undergo external perturbations and subsequent collisionless equilibration, predominantly via phase-mixing. We use linear perturbation theory to study the response of infinite isothermal slab analogues of disks to perturbations with diverse spatio-temporal characteristics. Without self-gravity of the response, the dominant Fourier modes that g… ▽ More Galactic disks are highly responsive systems that often undergo external perturbations and subsequent collisionless equilibration, predominantly via phase-mixing. We use linear perturbation theory to study the response of infinite isothermal slab analogues of disks to perturbations with diverse spatio-temporal characteristics. Without self-gravity of the response, the dominant Fourier modes that get excited in a disk are the bending and breathing modes, which, due to vertical phase-mixing, trigger local phase-space spirals that are one- and two-armed, respectively. We demonstrate how the lateral streaming motion of slab stars causes phase spirals to damp out over time. The ratio of the perturbation timescale ($τ_{\mathrm{P}}$) to the local, vertical oscillation time ($τ_z$) ultimately decides which of the two modes is excited. Faster, more impulsive ($τ_{\mathrm{P}} < τ_z$) and slower, more adiabatic ($τ_{\mathrm{P}} > τ_z$) perturbations excite stronger breathing and bending modes, respectively, although the response to very slow perturbations is exponentially suppressed. For encounters with satellite galaxies, this translates to more distant and more perpendicular encounters triggering stronger bending modes. We compute the direct response of the Milky Way disk to several of its satellite galaxies, and find that recent encounters with all of them excite bending modes in the Solar neighborhood. The encounter with Sagittarius triggers a response that is at least $1-2$ orders of magnitude larger than that due to any other satellite, including the Large Magellanic Cloud. We briefly discuss how ignoring the presence of a dark matter halo and the self-gravity of the response might impact our conclusions. △ Less

Submitted 9 August, 2022; originally announced August 2022.

Comments: Accepted for publication in ApJ; 7 figures, 1 table

Dynamically constraining the length of the Milky Way bar

Authors: Madeline Lucey, Sarah Pearson, Jason A. S. Hunt, Keith Hawkins, Melissa Ness, Michael S. Petersen, Adrian M. Price-Whelan, Martin D. Weinberg

Abstract: We present a novel method for constraining the length of the Galactic bar using 6D phase space information to directly integrate orbits. We define a pseudo-length for the Galactic bar, named $R_{Freq}$, based on the maximal extent of trapped bar orbits. We find the $R_{Freq}$ measured from orbits is consistent with the $R_{Freq}$ of the assumed potential only when the length of the bar and pattern… ▽ More We present a novel method for constraining the length of the Galactic bar using 6D phase space information to directly integrate orbits. We define a pseudo-length for the Galactic bar, named $R_{Freq}$, based on the maximal extent of trapped bar orbits. We find the $R_{Freq}$ measured from orbits is consistent with the $R_{Freq}$ of the assumed potential only when the length of the bar and pattern speed of said potential is similar to the model from which the initial phase-space coordinates of the orbits are derived. Therefore, one can measure the model's or the Milky Way's bar length from 6D phase-space coordinates by determining which assumed potential leads to a self-consistent measured $R_{Freq}$. When we apply this method to $\approx$210,000 stars in APOGEE DR17 and $Gaia$ eDR3 data, we find a consistent result only for potential models with a dynamical bar length of $\approx$3.5 kpc. We find the Milky Way's trapped bar orbits extend out to only $\approx$3.5 kpc, but there is also an overdensity of stars at the end of the bar out to 4.8 kpc which could be related to an attached spiral arm. We also find that the measured orbital structure of the bar is strongly dependent on the properties of the assumed potential. △ Less

Submitted 2 February, 2023; v1 submitted 3 June, 2022; originally announced June 2022.

Comments: 15 pages, 8 figures, 2 tables, accepted to MNRAS, comments welcome

Simulations of nonradiative processes in semiconductor nanocrystals

Authors: Dipti Jasrasaria, Daniel Weinberg, John P. Philbin, Eran Rabani

Abstract: The description of carrier dynamics in spatially confined semiconductor nanocrystals (NCs), which have enhanced electron-hole and exciton-phonon interactions, is a great challenge for modern computational science. These NCs typically contain thousands of atoms and tens of thousands of valence electrons with a discrete spectrum at low excitation energies, similar to atoms and molecules, that conver… ▽ More The description of carrier dynamics in spatially confined semiconductor nanocrystals (NCs), which have enhanced electron-hole and exciton-phonon interactions, is a great challenge for modern computational science. These NCs typically contain thousands of atoms and tens of thousands of valence electrons with a discrete spectrum at low excitation energies, similar to atoms and molecules, that converges to the continuum bulk limit at higher energies. Computational methods developed for molecules are limited to very small nanoclusters, and methods for bulk systems with periodic boundary conditions are not suitable due to the lack of translational symmetry in NCs. This perspective focuses on our recent efforts in developing a unified atomistic model based on the semiempirical pseudopotential approach, which is parametrized by first-principles calculations and validated against experimental measurements, to describe two of the main nonradiative relaxation processes of quantum confined excitons: exciton cooling and Auger recombination. We focus on the description of both electron-hole and exciton-phonon interactions in our approach and discuss the role of size, shape, and interfacing on the electronic properties and dynamics for II-VI and III-V semiconductor NCs. △ Less

Submitted 30 October, 2022; v1 submitted 13 April, 2022; originally announced April 2022.

Journal ref: J. Chem. Phys. 157, 020901 (2022)

LyMAS reloaded: improving the predictions of the large-scale Lyman-α forest statistics from dark matter density and velocity fields

Authors: S. Peirani, S. Prunet, S. Colombi, C. Pichon, D. H. Weinberg, C. Laigle, G. Lavaux, Y. Dubois, J. Devriendt

Abstract: We present LyMAS2, an improved version of the "Lyman-α Mass Association Scheme" aiming at predicting the large-scale 3d clustering statistics of the Lyman-α forest (Ly-α) from moderate resolution simulations of the dark matter (DM) distribution, with prior calibrations from high resolution hydrodynamical simulations of smaller volumes. In this study, calibrations are derived from the Horizon-AGN s… ▽ More We present LyMAS2, an improved version of the "Lyman-α Mass Association Scheme" aiming at predicting the large-scale 3d clustering statistics of the Lyman-α forest (Ly-α) from moderate resolution simulations of the dark matter (DM) distribution, with prior calibrations from high resolution hydrodynamical simulations of smaller volumes. In this study, calibrations are derived from the Horizon-AGN suite simulations, (100 Mpc/h)^3 comoving volume, using Wiener filtering, combining information from dark matter density and velocity fields (i.e. velocity dispersion, vorticity, line of sight 1d-divergence and 3d-divergence). All new predictions have been done at z=2.5 in redshift-space, while considering the spectral resolution of the SDSS-III BOSS Survey and different dark matter smoothing (0.3, 0.5 and 1.0 Mpc/h comoving). We have tried different combinations of dark matter fields and found that LyMAS2, applied to the Horizon-noAGN dark matter fields, significantly improves the predictions of the Ly-α 3d clustering statistics, especially when the DM overdensity is associated with the velocity dispersion or the vorticity fields. Compared to the hydrodynamical simulation trends, the 2-point correlation functions of pseudo-spectra generated with LyMAS2 can be recovered with relative differences of ~5% even for high angles, the flux 1d power spectrum (along the light of sight) with ~2% and the flux 1d probability distribution function exactly. Finally, we have produced several large mock BOSS spectra (1.0 and 1.5 Gpc/h) expected to lead to much more reliable and accurate theoretical predictions. △ Less

Submitted 13 April, 2022; originally announced April 2022.

Comments: Submitted to MNRAS (one iteration with the referee) - 29 pages

The chemical characterisation of halo substructure in the Milky Way based on APOGEE

Authors: Danny Horta, Ricardo P. Schiavon, J. Ted Mackereth, David H. Weinberg, Sten Hasselquist, Diane Feuillet, Robert W. O'Connell, Borja Anguiano, Carlos Allende-Prieto, Rachael L. Beaton, Dmitry Bizyaev, Katia Cunha, Doug Geisler, D. A. García-Hernández, Jon Holtzman, Henrik Jönsson, Richard R. Lane, Steve R. Majewski, Szabolcs Mészáros, Dante Minniti, Christian Nitschelm, Matthew Shetrone, Verne V. Smith, Gail Zasowski

Abstract: Galactic haloes in a $Λ$-CDM universe are predicted to host today a swarm of debris resulting from cannibalised dwarf galaxies. The chemo-dynamical information recorded in their stellar populations helps elucidate their nature, constraining the assembly history of the Galaxy. Using data from APOGEE and \textit{Gaia}, we examine the chemical properties of various halo substructures, considering ele… ▽ More Galactic haloes in a $Λ$-CDM universe are predicted to host today a swarm of debris resulting from cannibalised dwarf galaxies. The chemo-dynamical information recorded in their stellar populations helps elucidate their nature, constraining the assembly history of the Galaxy. Using data from APOGEE and \textit{Gaia}, we examine the chemical properties of various halo substructures, considering elements that sample various nucleosynthetic pathways. The systems studied are Heracles, \textit{Gaia}-Enceladus/Sausage (GES), the Helmi stream, Sequoia, Thamnos, Aleph, LMS-1, Arjuna, I'itoi, Nyx, Icarus, and Pontus. Abundance patterns of all substructures are cross-compared in a statistically robust fashion. Our main findings include: {\it i)} the chemical properties of most substructures studied match qualitatively those of dwarf Milky Way satellites, such as the Sagittarius dSph. Exceptions are Nyx and Aleph, which are chemically similar to disc stars, implying that these substructures were likely formed \textit{in situ}; {\it ii)} Heracles differs chemically from {\it in situ} populations such as Aurora and its inner halo counterparts in a statistically significant way. The differences suggest that the star formation rate was lower in Heracles than in the early Milky Way; {\it iii)} the chemistry of Arjuna, LMS-1, and I'itoi is indistinguishable from that of GES, suggesting a possible common origin; {\it iv)} all three Sequoia samples studied are qualitatively similar. However, only two of those samples present chemistry that is consistent with GES in a statistically significant fashion; {\it v)} the abundance patterns of the Helmi stream and Thamnos are different from all other halo substructures. △ Less

Submitted 2 November, 2022; v1 submitted 8 April, 2022; originally announced April 2022.

Comments: Accepted for publication in MNRAS. The revised version contains a detailed comparison of the chemical composition of Heracles against Aurora and its inner halo in situ counterparts as well as a revised discussion of the chemistry of various Sequoia samples in contrast with Gaia-Enceladus/Sausage

Birth of the Galactic Disk Revealed by the H3 Survey

Authors: Charlie Conroy, David H. Weinberg, Rohan P. Naidu, Tobias Buck, James W. Johnson, Phillip Cargile, Ana Bonaca, Nelson Caldwell, Vedant Chandra, Jiwon Jesse Han, Benjamin D. Johnson, Joshua S. Speagle, Yuan-Sen Ting, Turner Woody, Dennis Zaritsky

Abstract: We use chemistry ([alpha/Fe] and [Fe/H]), main sequence turnoff ages, and kinematics determined from H3 Survey spectroscopy and Gaia astrometry to identify the birth of the Galactic disk. We separate in-situ and accreted stars on the basis of angular momenta and eccentricities. The sequence of high-alpha in-situ stars persists down to at least [Fe/H]=-2.5 and shows unexpected non-monotonic behavio… ▽ More We use chemistry ([alpha/Fe] and [Fe/H]), main sequence turnoff ages, and kinematics determined from H3 Survey spectroscopy and Gaia astrometry to identify the birth of the Galactic disk. We separate in-situ and accreted stars on the basis of angular momenta and eccentricities. The sequence of high-alpha in-situ stars persists down to at least [Fe/H]=-2.5 and shows unexpected non-monotonic behavior: with increasing metallicity the population first declines in [alpha/Fe], then increases over the range -1.3<[Fe/H]<-0.7, and then declines again at higher metallicities. The number of stars in the in-situ population rapidly increases above [Fe/H]=-1. The average kinematics of these stars are hot and independent of metallicity at [Fe/H]<-1 and then become increasingly cold and disk-like at higher metallicities. The ages of the in-situ, high-alpha stars are uniformly very old (13 Gyr) at [Fe/H]<-1.3, and span a wider range (8-12 Gyr) at higher metallicities. Interpreting the chemistry with a simple chemical evolution model suggests that the non-monotonic behavior is due to a significant increase in star formation efficiency, which began 13 Gyr ago. These results support a picture in which the first 1 Gyr of the Galaxy was characterized by a "simmering phase" in which the star formation efficiency was low and the kinematics had substantial disorder with some net rotation. The disk then underwent a dramatic transformation to a "boiling phase", in which the star formation efficiency increased substantially, the kinematics became disk-like, and the number of stars formed increased tenfold. We interpret this transformation as the birth of the Galactic disk at z~4. The physical origin of this transformation is unclear and does not seem to be reproduced in current galaxy formation models. △ Less

Submitted 6 April, 2022; originally announced April 2022.

Comments: 13 pages, 10 figures, submitted to ApJ

Primordial helium-3 redux: The helium isotope ratio of the Orion nebula

Authors: Ryan J. Cooke, Pasquier Noterdaeme, James W. Johnson, Max Pettini, Louise Welsh, Celine Peroux, Michael T. Murphy, David H. Weinberg

Abstract: We report the first direct measurement of the helium isotope ratio, 3He/4He, outside of the Local Interstellar Cloud, as part of science verification observations with the upgraded CRyogenic InfraRed Echelle Spectrograph (CRIRES). Our determination of 3He/4He is based on metastable HeI* absorption along the line-of-sight towards Tet02 Ori A in the Orion Nebula. We measure a value 3He/4He=(1.77+/-0… ▽ More We report the first direct measurement of the helium isotope ratio, 3He/4He, outside of the Local Interstellar Cloud, as part of science verification observations with the upgraded CRyogenic InfraRed Echelle Spectrograph (CRIRES). Our determination of 3He/4He is based on metastable HeI* absorption along the line-of-sight towards Tet02 Ori A in the Orion Nebula. We measure a value 3He/4He=(1.77+/-0.13)x10^{-4}, which is just ~40 per cent above the primordial relative abundance of these isotopes, assuming the Standard Model of particle physics and cosmology, (3He/4He)_p = (1.257+/-0.017)x10^-4. We calculate a suite of galactic chemical evolution simulations to study the Galactic build up of these isotopes, using the yields from Limongi & Chieffi (2018) for stars in the mass range M=8-100 M_sun and Lagarde (2011,2012) for M=0.8-8 M_sun. We find that these simulations simultaneously reproduce the Orion and protosolar 3He/4He values if the calculations are initialized with a primordial ratio (3He/4He)_p=(1.043+/-0.089)x10^-4. Even though the quoted error does not include the model uncertainty, this determination agrees with the Standard Model value to within ~2sigma. We also use the present-day Galactic abundance of deuterium (D/H), helium (He/H), and 3He/4He to infer an empirical limit on the primordial 3He abundance, (3He/H)_p < (1.09+/-0.18)x10^-5, which also agrees with the Standard Model value. We point out that it is becoming increasingly difficult to explain the discrepant primordial 7Li/H abundance with non-standard physics, without breaking the remarkable simultaneous agreement of three primordial element ratios (D/H, 4He/H, and 3He/4He) with the Standard Model values. △ Less

Submitted 21 March, 2022; originally announced March 2022.

Comments: 23 pages, 9 figures, Resubmitted to the Astrophysical Journal after addressing referee comments

Optical selection bias and projection effects in stacked galaxy cluster weak lensing

Authors: Hao-Yi Wu, Matteo Costanzi, Chun-Hao To, Andrés N. Salcedo, David H. Weinberg, James Annis, Sebastian Bocquet, Maria Elidaiana da Silva Pereira, Joseph DeRose, Johnny Esteves, Arya Farahi, Sebastian Grandis, Eduardo Rozo, Eli S. Rykoff, Tamás N. Varga, Risa H. Wechsler, Chenxiao Zeng, Yuanyuan Zhang, Zhuowen Zhang

Abstract: Cosmological constraints from current and upcoming galaxy cluster surveys are limited by the accuracy of cluster mass calibration. In particular, optically identified galaxy clusters are prone to selection effects that can bias the weak lensing mass calibration. We investigate the selection bias of the stacked cluster lensing signal associated with optically selected clusters, using clusters ident… ▽ More Cosmological constraints from current and upcoming galaxy cluster surveys are limited by the accuracy of cluster mass calibration. In particular, optically identified galaxy clusters are prone to selection effects that can bias the weak lensing mass calibration. We investigate the selection bias of the stacked cluster lensing signal associated with optically selected clusters, using clusters identified by the redMaPPer algorithm in the Buzzard simulations as a case study. We find that at a given cluster halo mass, the residuals of redMaPPer richness and weak lensing signal are positively correlated. As a result, for a given richness selection, the stacked lensing signal is biased high compared with what we would expect from the underlying halo mass probability distribution. The cluster lensing selection bias can thus lead to overestimated mean cluster mass and biased cosmology results. We show that the lensing selection bias exhibits a strong scale-dependence and is approximately 20 to 60 percent for $ΔΣ$ at large scales. This selection bias largely originates from spurious member galaxies within +/- 20 to 60 Mpc/h along the line of sight, highlighting the importance of quantifying projection effects associated with the broad redshift distribution of member galaxies in photometric cluster surveys. While our results qualitatively agree with those in the literature, accurate quantitative modelling of the selection bias is needed to achieve the goals of cluster lensing cosmology and will require synthetic catalogues covering a wide range of galaxy-halo connection models. △ Less

Submitted 9 August, 2022; v1 submitted 10 March, 2022; originally announced March 2022.

Comments: 16 pages, 16 figures; replaced to match published version

Fast and selective super-resolution ultrasound in vivo with sono-switchable nanodroplets

Authors: Kai Riemer, Matthieu Toulemonde, Jipeng Yan, Marcelo Lerendegui, Eleanor Stride, Peter D. Weinberg, Christopher Dunsby, Meng-Xing Tang

Abstract: Perfusion by the microcirculation is key to the development, maintenance and pathology of tissue. Its measurement with high spatiotemporal resolution is consequently valuable but remains a challenge in deep tissue. Ultrasound Localization Microscopy (ULM) provides very high spatiotemporal resolution but the use of microbubbles requires low contrast agent concentrations, a long acquisition time, an… ▽ More Perfusion by the microcirculation is key to the development, maintenance and pathology of tissue. Its measurement with high spatiotemporal resolution is consequently valuable but remains a challenge in deep tissue. Ultrasound Localization Microscopy (ULM) provides very high spatiotemporal resolution but the use of microbubbles requires low contrast agent concentrations, a long acquisition time, and gives little control over the spatial and temporal distribution of the bubbles. The present study is the first to demonstrate Acoustic Wave Sparsely-Activated Localization Microscopy (AWSALM) and fast-AWSALM for in vivo super-resolution ultrasound imaging, offering contrast on demand and vascular selectivity. Three different formulations of sono-switchable contrast agents were tested. We demonstrate their use with ultrasound mechanical indices well within recommended safety limits to enable fast on-demand sparse switching at very high agent concentrations. We produce super-localization maps of the rabbit renal vasculature with acquisition times between 5.5 s and 0.25 s, and an 4-fold improvement in spatial resolution. We present the unique selectivity of AWSALM in visualizing specific vascular branches and downstream microvasculature, and we show super-localized kidney structures in systole and diastole with fast-AWSALM. In conclusion we demonstrate the feasibility of fast and selective measurement of microvascular dynamics in vivo with subwavelength resolution using ultrasound and sono-switchable nanodroplets. △ Less

Submitted 8 March, 2022; originally announced March 2022.

Comments: phase-change contrast agent, low-boiling point nanodroplet, acoustic vaporization, droplet activation, microcirculation, contrast enhanced ultrasound, plane wave

Empirical Constraints on the Nucleosynthesis of Nitrogen

Authors: James W. Johnson, David H. Weinberg, Fiorenzo Vincenzo, Jonathan C. Bird, Emily J. Griffith

Abstract: We derive empirical constraints on the nucleosynthetic yields of nitrogen by incorporating N enrichment into our previously developed and empirically tuned multi-zone galactic chemical evolution model. We adopt a metallicity-independent ("primary") N yield from massive stars and a metallicity-dependent ("secondary") N yield from AGB stars. In our model, galactic radial zones do not evolve along th… ▽ More We derive empirical constraints on the nucleosynthetic yields of nitrogen by incorporating N enrichment into our previously developed and empirically tuned multi-zone galactic chemical evolution model. We adopt a metallicity-independent ("primary") N yield from massive stars and a metallicity-dependent ("secondary") N yield from AGB stars. In our model, galactic radial zones do not evolve along the observed [N/O]-[O/H] relation, but first increase in [O/H] at roughly constant [N/O], then move upward in [N/O] via secondary N production. By $t\approx5$ Gyr, the model approaches an equilibrium [N/O]-[O/H] relation, which traces the radial oxygen gradient. We find good agreement with the [N/O]-[O/H] trend observed in extra-galactic systems if we adopt an IMF-averaged massive star yield $y_\text{N}^\text{CC}=3.6\times10^{-4}$, consistent with predictions for rapidly rotating progenitors, and a fractional AGB yield that is linear in mass and metallicity $y_\text{N}^\text{AGB}=(9\times10^{-4})(M_*/M_\odot)(Z_*/Z_\odot)$. This model reproduces the [N/O]-[O/H] relation found for Milky Way stars in the APOGEE survey, and it reproduces (though imperfectly) the trends of stellar [N/O] with age and [O/Fe]. The metallicity-dependent yield plays the dominant role in shaping the gas-phase [N/O]-[O/H] relation, but the AGB time-delay is required to match the APOGEE stellar age and [O/Fe] trends. If we add $\sim$40\% oscillations to the star formation rate, the model reproduces the scatter in gas-phase [N/O] vs. [O/H] observed in external galaxies by MaNGA. We also construct models using published AGB yields and examine their empirical successes and shortcomings. For all AGB yields we consider, simple stellar populations release half their N after only $\sim$250 Myr. △ Less

Submitted 9 February, 2022; originally announced February 2022.

Comments: 21 pages, 10 figures. See Figs. 5, 6, and 9 for key results. Submitted to MNRAS

Exploring the s-process history in the Galactic disk: Cerium abundances and gradients in Open Clusters from the OCCAM/APOGEE sample

Authors: J. V. Sales-Silva, S. Daflon, K. Cunha, D. Souto, V. V. Smith, C. Chiappini, J. Donor, P. M. Frinchaboy, D. A. García-Hernández, C. Hayes, S. R. Majewski, T. Masseron, R. P. Schiavon, D. H. Weinberg, R. L. Beaton, J. G. Fernández-Trincado, H. Jönsson, R. R. Lane, D. Minniti, A. Manchado, C. Moni Bidin, C. Nitschelm, J. O'Connell, S. Villanova

Abstract: The APOGEE Open Cluster Chemical Abundances and Mapping (OCCAM) survey is used to probe the chemical evolution of the s-process element cerium in the Galactic disk. Cerium abundances were derived from measurements of Ce II lines in the APOGEE spectra using the Brussels Automatic Code for Characterizing High Accuracy Spectra (BACCHUS) in 218 stars belonging to 42 open clusters. Our results indicate… ▽ More The APOGEE Open Cluster Chemical Abundances and Mapping (OCCAM) survey is used to probe the chemical evolution of the s-process element cerium in the Galactic disk. Cerium abundances were derived from measurements of Ce II lines in the APOGEE spectra using the Brussels Automatic Code for Characterizing High Accuracy Spectra (BACCHUS) in 218 stars belonging to 42 open clusters. Our results indicate that, in general, for Ages $<$ 4 Gyr, younger open clusters have higher [Ce/Fe] and [Ce/$α$-element] ratios than older clusters. In addition, metallicity segregates open clusters in the [Ce/X]-Age plane (where X can be H, Fe, and the $α$-elements O, Mg, Si, or Ca). These metallicity-dependant relations result in [Ce/Fe] and [Ce/$α$] ratios with age that are not universal clocks. Radial gradients of [Ce/H] and [Ce/Fe] ratios in open clusters, binned by age, were derived for the first time, with d[Ce/H]dR$_{GC}$ being negative, while d[Ce/Fe]/dR$_{GC}$ is positive. [Ce/H] and [Ce/Fe] gradients are approximately constant over time, with the [Ce/Fe] gradient becoming slightly steeper, changing by $\sim$+0.009 dex-kpc$^{-1}$-Gyr$^{-1}$. Both the [Ce/H] and [Ce/Fe] gradients are shifted to lower values of [Ce/H] and [Ce/Fe] for older open clusters. The chemical pattern of Ce in open clusters across the Galactic disk is discussed within the context of s-process yields from AGB stars, $\sim$Gyr time delays in Ce enrichment of the interstellar medium, and the strong dependence of Ce nucleosynthesis on the metallicity of its AGB stellar sources. △ Less

Submitted 12 December, 2021; v1 submitted 3 December, 2021; originally announced December 2021.

Comments: 20 pages, 10 figures, Accepted to ApJ