Showing 1–50 of 193 results for author: Schmidt, F

Straightening the Ruler: Field-Level Inference of the BAO Scale with LEFTfield

Authors: Ivana Babić, Fabian Schmidt, Beatriz Tucci

Abstract: Current inferences of the BAO scale from galaxy clustering employ a reconstruction technique at fixed cosmology and bias parameters. Here, we present the first consistent joint Bayesian inference of the isotropic BAO scale, jointly varying the initial conditions as well as all bias coefficients, based on the EFT-based field-level forward model $\texttt{LEFTfield}$. We apply this analysis to mock d… ▽ More Current inferences of the BAO scale from galaxy clustering employ a reconstruction technique at fixed cosmology and bias parameters. Here, we present the first consistent joint Bayesian inference of the isotropic BAO scale, jointly varying the initial conditions as well as all bias coefficients, based on the EFT-based field-level forward model $\texttt{LEFTfield}$. We apply this analysis to mock data generated with a much higher cutoff (resolution), resulting in a significant model mismatch between mock data and the model used in the inference. We demonstrate that the remaining systematic bias in the BAO scale is below 2% for all data considered and below 1% when Eulerian bias is used for inference. Furthermore, we find that the inferred error on the BAO scale is up to 1.6 times smaller compared to that from a replication of the standard power-spectrum reconstruction approach, using the same scales as in the field-level inference, with the improvement growing towards smaller scales (higher $k$). Thus, a field-level approach to BAO not only allows for a consistent inference of the BAO scale, but promises to achieve more precise measurements on the same data as well. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Comments: 25 pages, 15 figures

Shockingly Bright Warm Carbon Monoxide Molecular Features in the Supernova Remnant Cassiopeia A Revealed by JWST

Authors: J. Rho, S. -H. Park, R. Arendt, M. Matsuura, D. Milisavljevic, T. Temim, I. De Looze, W. P. Blair, A. Rest, O. Fox, A. P. Ravi, B. -C. Koo, M. Barlow, A. Burrows, R. Chevalier, G. Clayton, R. Fesen, C. Fransson, C. Fryer, H. L. Gomez, H. -T. Janka, F. Kirchschlarger, J. M. Laming, S. Orlando, D. Patnaude , et al. (14 additional authors not shown)

Abstract: We present JWST NIRCam (F356W and F444W filters) and MIRI (F770W) images and NIRSpec- IFU spectroscopy of the young supernova remnant Cassiopeia A (Cas A). We obtained the data as part of a JWST survey of Cas A. The NIRCam and MIRI images map the spatial distributions of synchrotron radiation, Ar-rich ejecta, and CO on both large and small scales, revealing remarkably complex structures. The CO em… ▽ More We present JWST NIRCam (F356W and F444W filters) and MIRI (F770W) images and NIRSpec- IFU spectroscopy of the young supernova remnant Cassiopeia A (Cas A). We obtained the data as part of a JWST survey of Cas A. The NIRCam and MIRI images map the spatial distributions of synchrotron radiation, Ar-rich ejecta, and CO on both large and small scales, revealing remarkably complex structures. The CO emission is stronger at the outer layers than the Ar ejecta, which indicates the reformation of CO molecules behind the reverse shock. NIRSpec-IFU spectra (3 - 5.5 microns) were obtained toward two representative knots in the NE and S fields. Both regions are dominated by the bright fundamental rovibrational band of CO in the two R and P branches, with strong [Ar VI] and relatively weaker, variable strength ejecta lines of [Si IX], [Ca IV], [Ca V] and [Mg IV]. The NIRSpec-IFU data resolve individual ejecta knots and filaments spatially and in velocity space. The fundamental CO band in the JWST spectra reveals unique shapes of CO, showing a few tens of sinusoidal patterns of rovibrational lines with pseudo-continuum underneath, which is attributed to the high-velocity widths of CO lines. The CO also shows high J lines at different vibrational transitions. Our results with LTE modeling of CO emission indicate a temperature of 1080 K and provide unique insight into the correlations between dust, molecules, and highly ionized ejecta in supernovae, and have strong ramifications for modeling dust formation that is led by CO cooling in the early Universe. △ Less

Submitted 5 June, 2024; originally announced June 2024.

Comments: accepted for the ApJ letter (17 pages and 10 figures)

The Renormalization Group for Large-Scale Structure: Primordial non-Gaussianities

Authors: Charalampos Nikolis, Henrique Rubira, Fabian Schmidt

Abstract: The renormalization group for large-scale structure (RG-LSS) describes the evolution of galaxy bias and stochastic parameters as a function of the cutoff $Λ$. In this work, we introduce interaction vertices that describe primordial non-Gaussianity into the Wilson-Polchinski framework, thereby extending the free theory to the interacting case. The presence of these interactions forces us to include… ▽ More The renormalization group for large-scale structure (RG-LSS) describes the evolution of galaxy bias and stochastic parameters as a function of the cutoff $Λ$. In this work, we introduce interaction vertices that describe primordial non-Gaussianity into the Wilson-Polchinski framework, thereby extending the free theory to the interacting case. The presence of these interactions forces us to include new operators and bias coefficients to the bias expansion to ensure closure under renormalization. We recover the previously-derived ``scale-dependent bias'' contributions, as well as a new (subdominant) stochastic contribution. We derive the renormalization group equations governing the RG-LSS for a large class of interactions which account for vertices at linear order in $f_{\rm NL}$ that parametrize interacting scalar and massive spinning fields during inflation. Solving the RG equations, we show the evolution of the non-Gaussian contributions to galaxy clustering as a function of scale. △ Less

Submitted 31 May, 2024; originally announced May 2024.

Comments: 51 pages, 9 figures

Euclid. I. Overview of the Euclid mission

Authors: Euclid Collaboration, Y. Mellier, Abdurro'uf, J. A. Acevedo Barroso, A. Achúcarro, J. Adamek, R. Adam, G. E. Addison, N. Aghanim, M. Aguena, V. Ajani, Y. Akrami, A. Al-Bahlawan, A. Alavi, I. S. Albuquerque, G. Alestas, G. Alguero, A. Allaoui, S. W. Allen, V. Allevato, A. V. Alonso-Tetilla, B. Altieri, A. Alvarez-Candal, A. Amara, L. Amendola , et al. (1086 additional authors not shown)

Abstract: The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14… ▽ More The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue`Euclid on Sky'

A Parameter-Masked Mock Data Challenge for Beyond-Two-Point Galaxy Clustering Statistics

Authors: Beyond-2pt Collaboration, :, Elisabeth Krause, Yosuke Kobayashi, Andrés N. Salcedo, Mikhail M. Ivanov, Tom Abel, Kazuyuki Akitsu, Raul E. Angulo, Giovanni Cabass, Sofia Contarini, Carolina Cuesta-Lazaro, ChangHoon Hahn, Nico Hamaus, Donghui Jeong, Chirag Modi, Nhat-Minh Nguyen, Takahiro Nishimichi, Enrique Paillas, Marcos Pellejero Ibañez, Oliver H. E. Philcox, Alice Pisani, Fabian Schmidt, Satoshi Tanaka, Giovanni Verza , et al. (2 additional authors not shown)

Abstract: The last few years have seen the emergence of a wide array of novel techniques for analyzing high-precision data from upcoming galaxy surveys, which aim to extend the statistical analysis of galaxy clustering data beyond the linear regime and the canonical two-point (2pt) statistics. We test and benchmark some of these new techniques in a community data challenge "Beyond-2pt", initiated during the… ▽ More The last few years have seen the emergence of a wide array of novel techniques for analyzing high-precision data from upcoming galaxy surveys, which aim to extend the statistical analysis of galaxy clustering data beyond the linear regime and the canonical two-point (2pt) statistics. We test and benchmark some of these new techniques in a community data challenge "Beyond-2pt", initiated during the Aspen 2022 Summer Program "Large-Scale Structure Cosmology beyond 2-Point Statistics," whose first round of results we present here. The challenge dataset consists of high-precision mock galaxy catalogs for clustering in real space, redshift space, and on a light cone. Participants in the challenge have developed end-to-end pipelines to analyze mock catalogs and extract unknown ("masked") cosmological parameters of the underlying $Λ$CDM models with their methods. The methods represented are density-split clustering, nearest neighbor statistics, BACCO power spectrum emulator, void statistics, LEFTfield field-level inference using effective field theory (EFT), and joint power spectrum and bispectrum analyses using both EFT and simulation-based inference. In this work, we review the results of the challenge, focusing on problems solved, lessons learned, and future research needed to perfect the emerging beyond-2pt approaches. The unbiased parameter recovery demonstrated in this challenge by multiple statistics and the associated modeling and inference frameworks supports the credibility of cosmology constraints from these methods. The challenge data set is publicly available and we welcome future submissions from methods that are not yet represented. △ Less

Submitted 3 May, 2024; originally announced May 2024.

Comments: New submissions welcome! Challenge data available at https://github.com/ANSalcedo/Beyond2ptMock

The Renormalization Group for Large-Scale Structure: Origin of Galaxy Stochasticity

Authors: Henrique Rubira, Fabian Schmidt

Abstract: The renormalization group equations for large-scale structure (RG-LSS) describe how the bias and stochastic (noise) parameters -- both of matter and biased tracers such as galaxies -- evolve as a function of the cutoff $Λ$ of the effective field theory. In previous work, we derived the RG-LSS equations for the bias parameters using the Wilson-Polchinski framework. Here, we extend these results to… ▽ More The renormalization group equations for large-scale structure (RG-LSS) describe how the bias and stochastic (noise) parameters -- both of matter and biased tracers such as galaxies -- evolve as a function of the cutoff $Λ$ of the effective field theory. In previous work, we derived the RG-LSS equations for the bias parameters using the Wilson-Polchinski framework. Here, we extend these results to include stochastic contributions, corresponding to terms in the effective action that are higher order in the current $J$. We show that the RG equations exhibit an interesting, previously unnoticed structure at all orders in $J$, which implies that a single nonlinear bias term immediately generates all stochastic moments through RG evolution. We then derive the nonlinear RG evolution of the (leading-derivative) stochastic parameters for all $n$-point functions, and show that this evolution is controlled by a different, lower scale than the nonlinear scale. This has implications for the optimal choice of the renormalization scale when comparing the theory with data to obtain cosmological constraints. △ Less

Submitted 25 April, 2024; originally announced April 2024.

Comments: 37 pages, 7 figures

Deep JWST/NIRCam imaging of Supernova 1987A

Authors: Mikako Matsuura, M. Boyer, Richard G. Arendt, J. Larsson, C. Fransson, A. Rest, A. P. Ravi, S. Park, P. Cigan, T. Temim, E. Dwek, M. J. Barlow, P. Bouchet, G. Clayton, R. Chevalier, J. Danziger, J. De Buizer, I. De Looze, G. De Marchi, O. Fox, C. Gall, R. D. Gehrz, H. L. Gomez, R. Indebetouw, T. Kangas , et al. (24 additional authors not shown)

Abstract: JWST/NIRCam obtained high angular-resolution (0.05-0.1''), deep near-infrared 1--5 micron imaging of Supernova (SN) 1987A taken 35 years after the explosion. In the NIRCam images, we identify: 1) faint H2 crescents, which are emissions located between the ejecta and the equatorial ring, 2) a bar, which is a substructure of the ejecta, and 3) the bright 3-5 micron continuum emission exterior to the… ▽ More JWST/NIRCam obtained high angular-resolution (0.05-0.1''), deep near-infrared 1--5 micron imaging of Supernova (SN) 1987A taken 35 years after the explosion. In the NIRCam images, we identify: 1) faint H2 crescents, which are emissions located between the ejecta and the equatorial ring, 2) a bar, which is a substructure of the ejecta, and 3) the bright 3-5 micron continuum emission exterior to the equatorial ring. The emission of the remnant in the NIRCam 1-2.3 micron images is mostly due to line emission, which is mostly emitted in the ejecta and in the hot spots within the equatorial ring. In contrast, the NIRCam 3-5 micron images are dominated by continuum emission. In the ejecta, the continuum is due to dust, obscuring the centre of the ejecta. In contrast, in the ring and exterior to the ring, synchrotron emission contributes a substantial fraction to the continuum. Dust emission contributes to the continuum at outer spots and diffuse emission exterior to the ring, but little within the ring. This shows that dust cooling and destruction time scales are shorter than the synchrotron cooling time scale, and the time scale of hydrogen recombination in the ring is even longer than the synchrotron cooling time scale. With the advent of high sensitivity and high angular resolution images provided by JWST/NIRCam, our observations of SN 1987A demonstrate that NIRCam opens up a window to study particle-acceleration and shock physics in unprecedented details, probed by near-infrared synchrotron emission, building a precise picture of how a SN evolves. △ Less

Submitted 15 April, 2024; originally announced April 2024.

Comments: Accepted for publication in MNRAS; 18 pages

Data availability and requirements relevant for the Ariel space mission and other exoplanet atmosphere applications

Authors: Katy L. Chubb, Séverine Robert, Clara Sousa-Silva, Sergei N. Yurchenko, Nicole F. Allard, Vincent Boudon, Jeanna Buldyreva, Benjamin Bultel, Athena Coustenis, Aleksandra Foltynowicz, Iouli E. Gordon, Robert J. Hargreaves, Christiane Helling, Christian Hill, Helgi Rafn Hrodmarsson, Tijs Karman, Helena Lecoq-Molinos, Alessandra Migliorini, Michaël Rey, Cyril Richard, Ibrahim Sadiek, Frédéric Schmidt, Andrei Sokolov, Stefania Stefani, Jonathan Tennyson , et al. (30 additional authors not shown)

Abstract: The goal of this white paper is to provide a snapshot of the data availability and data needs primarily for the Ariel space mission, but also for related atmospheric studies of exoplanets and brown dwarfs. It covers the following data-related topics: molecular and atomic line lists, line profiles, computed cross-sections and opacities, collision-induced absorption and other continuum data, optical… ▽ More The goal of this white paper is to provide a snapshot of the data availability and data needs primarily for the Ariel space mission, but also for related atmospheric studies of exoplanets and brown dwarfs. It covers the following data-related topics: molecular and atomic line lists, line profiles, computed cross-sections and opacities, collision-induced absorption and other continuum data, optical properties of aerosols and surfaces, atmospheric chemistry, UV photodissociation and photoabsorption cross-sections, and standards in the description and format of such data. These data aspects are discussed by addressing the following questions for each topic, based on the experience of the "data-provider" and "data-user" communities: (1) what are the types and sources of currently available data, (2) what work is currently in progress, and (3) what are the current and anticipated data needs. We present a GitHub platform for Ariel-related data, with the goal to provide a go-to place for both data-users and data-providers, for the users to make requests for their data needs and for the data-providers to link to their available data. Our aim throughout the paper is to provide practical information on existing sources of data whether in databases, theoretical, or literature sources. △ Less

Submitted 2 April, 2024; originally announced April 2024.

Comments: 58 pages, submitted to RAS Techniques and Instruments (RASTI). The authors welcome feedback: corresponding author emails can be found as footnotes on page 2

How much information can be extracted from galaxy clustering at the field level?

Authors: Nhat-Minh Nguyen, Fabian Schmidt, Beatriz Tucci, Martin Reinecke, Andrija Kostić

Abstract: We present optimal Bayesian field-level cosmological constraints from nonlinear tracers of the large-scale structure, specifically the amplitude $σ_8$ of linear matter fluctuations inferred from rest-frame simulated dark matter halos in a comoving volume of $8\,(h^{-1}\mathrm{Gpc})^3$. Our constraint on $σ_8$ is entirely due to nonlinear information, and obtained by explicitly sampling the initial… ▽ More We present optimal Bayesian field-level cosmological constraints from nonlinear tracers of the large-scale structure, specifically the amplitude $σ_8$ of linear matter fluctuations inferred from rest-frame simulated dark matter halos in a comoving volume of $8\,(h^{-1}\mathrm{Gpc})^3$. Our constraint on $σ_8$ is entirely due to nonlinear information, and obtained by explicitly sampling the initial conditions along with bias and noise parameters via a Lagrangian EFT-based forward model, LEFTfield. The comparison with a simulation-based inference analysis employing the power spectrum and bispectrum -- likewise using the LEFTfield forward model -- shows that, when including precisely the same modes of the same data up to $k_{\mathrm{max}}= 0.10\,h\,\mathrm{Mpc}^{-1}$ ($0.12\,h\,\mathrm{Mpc}^{-1}$), the field-level approach yields a factor of 3.5 (5.2) improvement on the $σ_8$ constraint, from 20.0% to 5.7% (17.0% to 3.3%). This study provides direct insights into cosmological information encoded in galaxy clustering beyond low-order $n$-point functions. △ Less

Submitted 2 April, 2024; v1 submitted 5 March, 2024; originally announced March 2024.

Comments: 3750 words + 3 appendices. 17 pages (including references), 12 figures. Main results in Figs. 2 & 3. All comments welcome! v2: Less words and more results. SBI P+B results, comparisons and discussion thereof updated. SBI P+B and FBI now use exactly the same forward model, yielding completely consistent posterior means

Report number: LCTP-24-05

From total destruction to complete survival: Dust processing at different evolutionary stages in the supernova remnant Cassiopeia A

Authors: Florian Kirchschlager, Nina Sartorio, Ilse De Looze, M. J. Barlow, Franziska Schmidt, Felix Priestley

Abstract: The expanding ejecta of supernova remnants (SNRs) are believed to form dust in dense clumps of gas. Before the dust can be expelled into the interstellar medium and contribute to the interstellar dust budget, it has to survive the reverse shock that is generated through the interaction of the preceding supernova blast wave with the surrounding medium. The conditions under which the reverse shock h… ▽ More The expanding ejecta of supernova remnants (SNRs) are believed to form dust in dense clumps of gas. Before the dust can be expelled into the interstellar medium and contribute to the interstellar dust budget, it has to survive the reverse shock that is generated through the interaction of the preceding supernova blast wave with the surrounding medium. The conditions under which the reverse shock hits the clumps change with remnant age and define the dust survival rate. To study the dust destruction in the SNR Cassiopeia A, we conduct magnetohydrodynamical simulations of the evolution of a supernova blast wave and of the reverse shock. In a second step we use these evolving conditions to model clumps that are disrupted by the reverse shock at different remnant ages. Finally, we compute the amount of dust that is destroyed by the impact of the reverse shock. We find that most of the dust in the SNR is hit by the reverse shock within the first 350 yr after the SN explosion. While the dust destruction in the first 200 yr is almost complete, we expect greater dust survival rates at later times and almost total survival for clumps that are first impacted at ages beyond 1000 yr. Integrated over the entire evolution of the SNR, the dust mass shows the lowest survival fraction (17 per cent) for the smallest grains (1 nm) and the highest survival fraction (28 per cent) for the largest grains (1000 nm). △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: Accepted by MNRAS. Author accepted manuscript. Accepted on 31/01/2024. Deposited on 01/02/2024. 12 pages

Constraining $f(R)$ gravity using future galaxy cluster abundance and weak-lensing mass calibration datasets

Authors: Sophie M. L. Vogt, Sebastian Bocquet, Christopher T. Davies, Joseph J. Mohr, Fabian Schmidt

Abstract: We present forecasts for constraints on the Hu \& Sawicki $f(R)$ modified gravity model using realistic mock data representative of future cluster and weak lensing surveys. We create mock thermal Sunyaev-Zel'dovich effect selected cluster samples for SPT-3G and CMB-S4 and the corresponding weak gravitational lensing data from next-generation weak-lensing (ngWL) surveys like Euclid and Rubin. We em… ▽ More We present forecasts for constraints on the Hu \& Sawicki $f(R)$ modified gravity model using realistic mock data representative of future cluster and weak lensing surveys. We create mock thermal Sunyaev-Zel'dovich effect selected cluster samples for SPT-3G and CMB-S4 and the corresponding weak gravitational lensing data from next-generation weak-lensing (ngWL) surveys like Euclid and Rubin. We employ a state-of-the-art Bayesian likelihood approach that includes all observational effects and systematic uncertainties to obtain constraints on the $f(R)$ gravity parameter $\log_{10}|f_{R0}|$. In this analysis we vary the cosmological parameters $[Ω_{\rm m}, Ω_νh^2, h^2, A_s, n_s, \log_{10}|f_{R0}|]$, which allows us to account for possible degeneracies between cosmological parameters and $f(R)$ modified gravity. The analysis accounts for $f(R)$ gravity via its effect on the halo mass function which is enhanced on cluster mass scales compared to the expectations within general relativity (GR). Assuming a fiducial GR model, the upcoming cluster dataset SPT-3G$\times$ngWL is expected to obtain an upper limit of $\log_{10}|f_{R0}| < -5.95$ at $95\,\%$ credibility, which significantly improves upon the current best bounds. The CMB-S4$\times$ngWL dataset is expected to improve this even further to $\log_{10}|f_{R0}| < -6.23$. Furthermore, $f(R)$ gravity models with $\log_{10}|f_{R0}| \geq -6$, which have larger numbers of clusters, would be distinguishable from GR with both datasets. We also report degeneracies between $\log_{10}|f_{R0}|$ and $Ω_{\mathrm{m}}$ as well as $σ_8$ for $\log_{10}|f_{R0}| > -6$ and $\log_{10}|f_{R0}| > -5$ respectively. Our forecasts indicate that future cluster abundance studies of $f(R)$ gravity will enable substantially improved constraints that are competitive with other cosmological probes. △ Less

Submitted 10 June, 2024; v1 submitted 18 January, 2024; originally announced January 2024.

Comments: 24 pages, 12 figures, minor changes in version 2 to match version in PRD, published in PRD: https://journals.aps.org/prd/abstract/10.1103/PhysRevD.109.123503

A JWST Survey of the Supernova Remnant Cassiopeia A

Authors: Dan Milisavljevic, Tea Temim, Ilse De Looze, Danielle Dickinson, J. Martin Laming, Robert Fesen, John C. Raymond, Richard G. Arendt, Jacco Vink, Bettina Posselt, George G. Pavlov, Ori D. Fox, Ethan Pinarski, Bhagya Subrayan, Judy Schmidt, William P. Blair, Armin Rest, Daniel Patnaude, Bon-Chul Koo, Jeonghee Rho, Salvatore Orlando, Hans-Thomas Janka, Moira Andrews, Michael J. Barlow, Adam Burrows , et al. (21 additional authors not shown)

Abstract: We present initial results from a JWST survey of the youngest Galactic core-collapse supernova remnant Cassiopeia A (Cas A), made up of NIRCam and MIRI imaging mosaics that map emission from the main shell, interior, and surrounding circumstellar/interstellar material (CSM/ISM). We also present four exploratory positions of MIRI/MRS IFU spectroscopy that sample ejecta, CSM, and associated dust fro… ▽ More We present initial results from a JWST survey of the youngest Galactic core-collapse supernova remnant Cassiopeia A (Cas A), made up of NIRCam and MIRI imaging mosaics that map emission from the main shell, interior, and surrounding circumstellar/interstellar material (CSM/ISM). We also present four exploratory positions of MIRI/MRS IFU spectroscopy that sample ejecta, CSM, and associated dust from representative shocked and unshocked regions. Surprising discoveries include: 1) a web-like network of unshocked ejecta filaments resolved to 0.01 pc scales exhibiting an overall morphology consistent with turbulent mixing of cool, low-entropy matter from the progenitor's oxygen layer with hot, high-entropy matter heated by neutrino interactions and radioactivity, 2) a thick sheet of dust-dominated emission from shocked CSM seen in projection toward the remnant's interior pockmarked with small (approximately one arcsecond) round holes formed by knots of high-velocity ejecta that have pierced through the CSM and driven expanding tangential shocks, 3) dozens of light echoes with angular sizes between 0.1 arcsecond to 1 arcminute reflecting previously unseen fine-scale structure in the ISM. NIRCam observations place new upper limits on infrared emission from the neutron star in Cas A's center and tightly constrain scenarios involving a possible fallback disk. These JWST survey data and initial findings help address unresolved questions about massive star explosions that have broad implications for the formation and evolution of stellar populations, the metal and dust enrichment of galaxies, and the origin of compact remnant objects. △ Less

Submitted 10 June, 2024; v1 submitted 4 January, 2024; originally announced January 2024.

Comments: 27 pages, 10 figures, now published in ApJL

Journal ref: ApJL 965 (2024) L27 (21pp)

Gravity-Induced Ice Compaction and Subsurface Porosity on Icy Moons

Authors: Cyril Mergny, Frédéric Schmidt

Abstract: Our understanding of the surface porosity of icy moons and its evolution with depth remains limited, including the precise scale at which ice compaction occurs under self-weight pressure. This parameter is of crucial interest for the correct interpretation of current remote sensing data (spectroscopy in the visible, infrared to passive microwave) but also for planetary exploration when designing a… ▽ More Our understanding of the surface porosity of icy moons and its evolution with depth remains limited, including the precise scale at which ice compaction occurs under self-weight pressure. This parameter is of crucial interest for the correct interpretation of current remote sensing data (spectroscopy in the visible, infrared to passive microwave) but also for planetary exploration when designing a lander, a rover or a cryobot. In situ exploration of the ice crust would require knowledge about subsurface porosity. This study employs a compaction model solely driven by overburden pressure based on prior research. The formulation for density as a function of depth, incorporates an essential parameter: the ice compaction coefficient. To determine this coefficient, we fit our depth-dependent density model to existing data obtained from Earth-based measurements of ice cores in Antarctica and North Greenland. Our results yield a typical lengthscale for ice compaction on Earth of approximately 20.1 $\pm$ 0.6 m , consistent with the existing literature. We apply the model to Europa, which due to its lower gravity, has a typical ice compaction scale of 150 $\pm$ 4 m. We compare it with the depths scanned by current spaceborne data and find that porosity can be considered constant when accounting only for gravity-induced compaction. △ Less

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

Journal ref: Icarus, Volume 413, 2024

EFTofLSS meets simulation-based inference: $σ_8$ from biased tracers

Authors: Beatriz Tucci, Fabian Schmidt

Abstract: Cosmological inferences typically rely on explicit expressions for the likelihood and covariance of the data vector, which normally consists of a set of summary statistics. However, in the case of nonlinear large-scale structure, exact expressions for either likelihood or covariance are unknown, and even approximate expressions can become very cumbersome, depending on the scales and summary statis… ▽ More Cosmological inferences typically rely on explicit expressions for the likelihood and covariance of the data vector, which normally consists of a set of summary statistics. However, in the case of nonlinear large-scale structure, exact expressions for either likelihood or covariance are unknown, and even approximate expressions can become very cumbersome, depending on the scales and summary statistics considered. Simulation-based inference (SBI), in contrast, does not require an explicit form for the likelihood but only a prior and a simulator, thereby naturally circumventing these issues. In this paper, we explore how this technique can be used to infer $σ_8$ from a Lagrangian effective field theory (EFT) based forward model for biased tracers. The power spectrum and bispectrum are used as summary statistics to obtain the posterior of the cosmological, bias and noise parameters via neural density estimation. We compare full simulation-based inference with cases where the data vector is drawn from a Gaussian likelihood with sample and analytical covariances. We conclude that, for $k_{\text{max}}=0.1h\text{Mpc}^{-1}$ and $0.2h\text{Mpc}^{-1}$, the form of the covariance is more important than the non-Gaussianity of the likelihood, although this conclusion is expected to depend on the cosmological parameter inferred, the summary statistics considered and range of scales probed. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: 41 pages, 17 figures

Galaxy bias renormalization group

Authors: Henrique Rubira, Fabian Schmidt

Abstract: The effective field theory of large-scale structure allows for a consistent perturbative bias expansion of the rest-frame galaxy density field. In this work, we present a systematic approach to renormalize galaxy bias and stochastic parameters using a finite cutoff scale $Λ$. We derive the differential equations of the Wilson-Polchinski renormalization group that describe the evolution of the fini… ▽ More The effective field theory of large-scale structure allows for a consistent perturbative bias expansion of the rest-frame galaxy density field. In this work, we present a systematic approach to renormalize galaxy bias and stochastic parameters using a finite cutoff scale $Λ$. We derive the differential equations of the Wilson-Polchinski renormalization group that describe the evolution of the finite-scale bias parameters with $Λ$, analogous to the $β$-function running in QFT. We further provide the connection between the finite-cutoff scheme and the renormalization procedure for $n$-point functions that has been used as standard in the literature so far; some inconsistencies in the treatment of renormalized bias in current EFT analyses are pointed out as well. The fixed-cutoff scheme allows us to predict, in a principled way, the finite part of loop contributions which is due to perturbative modes and which, in the standard renormalization approach, is absorbed into counterterms. We expect that this will allow for the robust extraction of (a yet-to-be-determined amount of) additional cosmological information from galaxy clustering, both when using field-level techniques and $n$-point functions. △ Less

Submitted 27 July, 2023; originally announced July 2023.

Comments: 40 page, 6 figure

Report number: TUM-HEP- 1464/23

Halo assembly bias from a deep learning model of halo formation

Authors: Luisa Lucie-Smith, Alexandre Barreira, Fabian Schmidt

Abstract: We build a deep learning framework that connects the local formation process of dark matter halos to the halo bias. We train a convolutional neural network (CNN) to predict the final mass and concentration of dark matter halos from the initial conditions. The CNN is then used as a surrogate model to derive the response of the halos' mass and concentration to long-wavelength perturbations in the in… ▽ More We build a deep learning framework that connects the local formation process of dark matter halos to the halo bias. We train a convolutional neural network (CNN) to predict the final mass and concentration of dark matter halos from the initial conditions. The CNN is then used as a surrogate model to derive the response of the halos' mass and concentration to long-wavelength perturbations in the initial conditions, and consequently the halo bias parameters following the "response bias" definition. The CNN correctly predicts how the local properties of dark matter halos respond to changes in the large-scale environment, despite no explicit knowledge of halo bias being provided during training. We show that the CNN recovers the known trends for the linear and second-order density bias parameters $b_1$ and $b_2$, as well as for the local primordial non-Gaussianity linear bias parameter $b_φ$. The expected secondary assembly bias dependence on halo concentration is also recovered by the CNN: at fixed mass, halo concentration has only a mild impact on $b_1$, but a strong impact on $b_φ$. Our framework opens a new window for discovering which physical aspects of the halo's Lagrangian patch determine assembly bias, which in turn can inform physical models of halo formation and bias. △ Less

Submitted 19 April, 2023; originally announced April 2023.

Comments: 11 pages, 5 figures, to be submitted to MNRAS, comments welcome

Effective Field Theory of Intrinsic Alignments at One Loop Order: a Comparison to Dark Matter Simulations

Authors: Thomas Bakx, Toshiki Kurita, Nora Elisa Chisari, Zvonimir Vlah, Fabian Schmidt

Abstract: We test the regime of validity of the effective field theory (EFT) of intrinsic alignments (IA) at the one-loop level by comparing with 3D halo shape statistics in N-body simulations. This model is based on the effective field theory of large-scale structure (EFT of LSS) and thus a theoretically well-motivated extension of the familiar non-linear alignment (NLA) model and the tidal-alignment-tidal… ▽ More We test the regime of validity of the effective field theory (EFT) of intrinsic alignments (IA) at the one-loop level by comparing with 3D halo shape statistics in N-body simulations. This model is based on the effective field theory of large-scale structure (EFT of LSS) and thus a theoretically well-motivated extension of the familiar non-linear alignment (NLA) model and the tidal-alignment-tidal-torquing (TATT) model. It contains a total of $8$ free bias parameters. Specifically, we measure the dark matter halo shape-shape multipoles $P_{EE}^{(0)}(k), P_{EE}^{(2)}(k), P_{BB}^{(0)}(k), P_{BB}^{(2)}(k)$ as well as the matter-shape multipoles $P_{δE}^{(0)}(k), P_{δE}^{(2)}(k)$ from the simulations and perform a joint fit to determine the largest wavenumber $k_{\text{max}}$ up to which the theory predictions from the EFT of IA are consistent with the measurements. We find that the EFT of IA is able to describe intrinsic alignments of dark matter halos up to $k_\text{max}=0.30\,h/\text{Mpc}$ at $z=0$. This demonstrates a clear improvement over other existing alignment models like NLA and TATT, which are only accurate up to $k_\text{max}=0.05\,h/\text{Mpc}$ . We examine the posterior distributions of the higher-order bias parameters, and show that their inclusion is necessary to describe intrinsic alignments in the quasi-linear regime. Further, the EFT of IA is able to accurately describe the auto-spectrum of intrinsic alignment B-modes, in contrast to the other alignment models considered. △ Less

Submitted 27 March, 2023; originally announced March 2023.

Comments: 35 pages, 2 tables, 6 figures

Cosmology inference at the field level from biased tracers in redshift-space

Authors: Julia Stadler, Fabian Schmidt, Martin Reinecke

Abstract: Cosmology inference of galaxy clustering at the field level with the EFT likelihood in principle allows for extracting all non-Gaussian information from quasi-linear scales, while robustly marginalizing over any astrophysical uncertainties. A pipeline in this spirit is implemented in the \texttt{LEFTfield} code, which we extend in this work to describe the clustering of galaxies in redshift space.… ▽ More Cosmology inference of galaxy clustering at the field level with the EFT likelihood in principle allows for extracting all non-Gaussian information from quasi-linear scales, while robustly marginalizing over any astrophysical uncertainties. A pipeline in this spirit is implemented in the \texttt{LEFTfield} code, which we extend in this work to describe the clustering of galaxies in redshift space. Our main additions are: the computation of the velocity field in the LPT gravity model, the fully nonlinear displacement of the evolved, biased density field to redshift space, and a systematic expansion of velocity bias. We test the resulting analysis pipeline by applying it to synthetic data sets with a known ground truth at increasing complexity: mock data generated from the perturbative forward model itself, sub-sampled matter particles, and dark matter halos in N-body simulations. By fixing the initial-time density contrast to the ground truth, while varying the growth rate $f$, bias coefficients and noise amplitudes, we perform a stringent set of checks. These show that indeed a systematic higher-order expansion of the velocity bias is required to infer a growth rate consistent with the ground truth within errors. Applied to dark matter halos, our analysis yields unbiased constraints on $f$ at the level of a few percent for a variety of halo masses at redshifts $z=0,\,0.5,\,1$ and for a broad range of cutoff scales $0.08\,h/\mathrm{Mpc} \leq Λ\leq 0.20\,h/\mathrm{Mpc}$. Importantly, deviations between true and inferred growth rate exhibit the scaling with halo mass, redshift and cutoff that one expects based on the EFT of Large Scale Structure. Further, we obtain a robust detection of velocity bias through its effect on the redshift-space density field and are able to disentangle it from higher-derivative bias contributions. △ Less

Submitted 30 October, 2023; v1 submitted 17 March, 2023; originally announced March 2023.

Journal ref: JCAP10(2023)069

Approximation of radiative transfer for surface spectral features

Authors: Frédéric Schmidt

Abstract: Remote sensing hyperspectral and more generally spectral instruments are common tools to decipher surface features in Earth and Planetary science. While linear mixture is the most common approximation for compounds detection (mineral, water, ice, etc...), the transfer of light in surface and atmospheric medium are highly non-linear. The exact simulation of non-linearities can be estimated at very… ▽ More Remote sensing hyperspectral and more generally spectral instruments are common tools to decipher surface features in Earth and Planetary science. While linear mixture is the most common approximation for compounds detection (mineral, water, ice, etc...), the transfer of light in surface and atmospheric medium are highly non-linear. The exact simulation of non-linearities can be estimated at very high numerical cost. Here I propose a very simple non-linear form (that includes the regular linear area mixture) of radiative transfer to approximate surface spectral feature. I demonstrate that this analytical form is able to approximate the grain size and intimate mixture dependence of surface features. In addition, the same analytical form can approximate the effect of Martian mineral aerosols. Unfortunately, Earth aerosols are more complex (water droplet, water ice, soot,...) and are not expected to follow the same trend. △ Less

Submitted 13 April, 2023; v1 submitted 6 February, 2023; originally announced February 2023.

Comments: 4 pages, 3 figures, submitted 21st october 2022 to IEEE Geoscience and Remote Sensing Letters

Journal ref: IEEE Geoscience and Remote Sensing Letters, 2023, 20, 1-3

Consistency tests of field level inference with the EFT likelihood

Authors: Andrija Kostić, Nhat-Minh Nguyen, Fabian Schmidt, Martin Reinecke

Abstract: Analyzing the clustering of galaxies at the field level in principle promises access to all the cosmological information available. Given this incentive, in this paper we investigate the performance of field-based forward modeling approach to galaxy clustering using the effective field theory (EFT) framework of large-scale structure (LSS). We do so by applying this formalism to a set of consistenc… ▽ More Analyzing the clustering of galaxies at the field level in principle promises access to all the cosmological information available. Given this incentive, in this paper we investigate the performance of field-based forward modeling approach to galaxy clustering using the effective field theory (EFT) framework of large-scale structure (LSS). We do so by applying this formalism to a set of consistency and convergence tests on synthetic datasets. We explore the high-dimensional joint posterior of LSS initial conditions by combining Hamiltonian Monte Carlo sampling for the field of initial conditions, and slice sampling for cosmology and model parameters. We adopt the Lagrangian perturbation theory forward model from [1], up to second order, for the forward model of biased tracers. We specifically include model mis-specifications in our synthetic datasets within the EFT framework. We achieve this by generating synthetic data at a higher cutoff scale $Λ_0$, which controls which Fourier modes enter the EFT likelihood evaluation, than the cutoff $Λ$ used in the inference. In the presence of model mis-specifications, we find that the EFT framework still allows for robust, unbiased joint inference of a) cosmological parameters - specifically, the scaling amplitude of the initial conditions - b) the initial conditions themselves, and c) the bias and noise parameters. In addition, we show that in the purely linear case, where the posterior is analytically tractable, our samplers fully explore the posterior surface. We also demonstrate convergence in the cases of nonlinear forward models. Our findings serve as a confirmation of the EFT field-based forward model framework developed in [2-7], and as another step towards field-level cosmological analyses of real galaxy surveys. △ Less

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

Comments: 31 + 13 pages, 15 figures; Added 3 new figures, text cleanup and fix typos; matching the version to be published in JCAP

Journal ref: JCAP 07(2023) 063

Dust survival rates in clumps passing through the Cas A reverse shock -- II. The impact of magnetic fields

Authors: Florian Kirchschlager, Franziska D. Schmidt, M. J. Barlow, Ilse De Looze, Nina S. Sartorio

Abstract: Dust grains form in the clumpy ejecta of core-collapse supernovae where they are subject to the reverse shock, which is able to disrupt the clumps and destroy the grains. Important dust destruction processes include thermal and kinetic sputtering as well as fragmentation and grain vaporization. In the present study, we focus on the effect of magnetic fields on the destruction processes. We have pe… ▽ More Dust grains form in the clumpy ejecta of core-collapse supernovae where they are subject to the reverse shock, which is able to disrupt the clumps and destroy the grains. Important dust destruction processes include thermal and kinetic sputtering as well as fragmentation and grain vaporization. In the present study, we focus on the effect of magnetic fields on the destruction processes. We have performed magneto-hydrodynamical simulations using AstroBEAR to model a shock wave interacting with an ejecta clump. The dust transport and destruction fractions are computed using our post-processing code Paperboats in which the acceleration of grains due to the magnetic field and a procedure that allows partial grain vaporization have been newly implemented. For the oxygen-rich supernova remnant Cassiopeia A we found a significantly lower dust survival rate when magnetic fields are aligned perpendicular to the shock direction compared to the non-magnetic case. For a parallel field alignment, the destruction is also enhanced but at a lower level. The survival fractions depend sensitively on the gas density contrast between the clump and the ambient medium and on the grain sizes. For a low-density contrast of $100$, e.g., $5\,$nm silicate grains are completely destroyed while the survival fraction of $1\,μ$m grains is $86\,$per cent. For a high-density contrast of $1000$, $95\,$per cent of the $5\,$nm grains survive while the survival fraction of $1\,μ$m grains is $26\,$per cent. Alternative clump sizes or dust materials (carbon) have non-negligible effects on the survival rate but have a lower impact compared to density contrast, magnetic field strength, and grain size. △ Less

Submitted 16 February, 2023; v1 submitted 13 October, 2022; originally announced October 2022.

Comments: Accepted by MNRAS. Author accepted manuscript. Accepted on 23/01/2023. 24 pages, 21 Figures

Non-linear Eulerian Hydrodynamics of Dark Energy: Riemann problem and Finite Volume Schemes

Authors: Linda Blot, Pier Stefano Corasaniti, Fabian Schmidt

Abstract: Upcoming large-scale-structure surveys can shed new light on the properties of dark energy. In particular, if dark energy is a dynamical component, it must have spatial perturbations. Their behaviour is regulated by the speed of sound parameter, which is currently unconstrained. In this work we present the numerical methods that will allow to perform cosmological simulations of inhomogeneous dark… ▽ More Upcoming large-scale-structure surveys can shed new light on the properties of dark energy. In particular, if dark energy is a dynamical component, it must have spatial perturbations. Their behaviour is regulated by the speed of sound parameter, which is currently unconstrained. In this work we present the numerical methods that will allow to perform cosmological simulations of inhomogeneous dark energy scenarios where the speed of sound is small and non-vanishing. We treat the dark energy component as an effective fluid and build upon established numerical methods for hydrodynamics to construct a numerical solution of the effective continuity and Euler equations. In particular, we develop conservative finite volume schemes that rely on the solution of the Riemann problem, which we provide here in both exact and approximate forms for the case of a dark energy fluid. △ Less

Submitted 21 April, 2023; v1 submitted 10 October, 2022; originally announced October 2022.

Comments: 38 pages, 22 figures, replaced to match the published version on JCAP. Significant changes in Section 4.2 and new advection test in Section 5.4

Journal ref: JCAP, 2023, 029P, 1022

Assembly bias in the local PNG halo bias and its implication for $f_{\rm NL}$ constraints

Authors: Titouan Lazeyras, Alexandre Barreira, Fabian Schmidt, Vincent Desjacques

Abstract: We use $N$-body simulations to study halo assembly bias (i.e., the dependence of halo clustering on properties beyond total mass) in the density and primordial non-Gaussianity (PNG) linear bias parameters $b_1$ and $b_φ$, respectively. We consider concentration, spin and sphericity as secondary halo properties, for which we find a clear detection of assembly bias for $b_1$ and $b_φ$. At fixed tota… ▽ More We use $N$-body simulations to study halo assembly bias (i.e., the dependence of halo clustering on properties beyond total mass) in the density and primordial non-Gaussianity (PNG) linear bias parameters $b_1$ and $b_φ$, respectively. We consider concentration, spin and sphericity as secondary halo properties, for which we find a clear detection of assembly bias for $b_1$ and $b_φ$. At fixed total mass, halo spin and sphericity impact $b_1$ and $b_φ$ in a similar manner, roughly preserving the shape of the linear $b_φ(b_1)$ relation satisfied by the global halo population. Halo concentration, however, drives $b_1$ and $b_φ$ in opposite directions. This induces significant changes to the $b_φ(b_1)$ relation, with higher concentration halos having higher amplitude of $b_φ(b_1)$. For $z=0.5$ and $b_1 \approx 2$ in particular, the population comprising either all halos, those with the $33\%$ lowest or those with the $33\%$ highest concentrations have a PNG bias of $b_φ\approx 3$, $b_φ\approx -1$ and $b_φ\approx 9$, respectively. Varying the halo concentration can make $b_φ$ very small and even change its sign. These results have important ramifications for galaxy clustering constraints of the local PNG parameter $f_{\rm NL}$ that assume fixed forms for the $b_φ(b_1)$ relation. We illustrate the significant impact of halo assembly bias in actual data using the BOSS DR12 galaxy power spectrum: assuming that BOSS galaxies are representative of all halos, the $33\%$ lowest or the $33\%$ highest concentration halos yields $σ_{f_{\rm NL}} = 44, 165, 19$, respectively. Our results suggest taking host halo concentration into account in galaxy selection strategies to maximize the signal-to-noise on $f_{\rm NL}$. They also motivate more simulation-based efforts to study the $b_φ(b_1)$ relation of halos and galaxies. △ Less

Submitted 19 January, 2023; v1 submitted 15 September, 2022; originally announced September 2022.

Comments: 21 pages, 4 figures, comments welcome; V2 minor modifications throughout, published in JCAP

Journal ref: JCAP 01 (2023) 023

Bayesian field-level inference of primordial non-Gaussianity using next-generation galaxy surveys

Authors: Adam Andrews, Jens Jasche, Guilhem Lavaux, Fabian Schmidt

Abstract: Detecting and measuring a non-Gaussian signature of primordial origin in the density field is a major science goal of next-generation galaxy surveys. The signal will permit us to determine primordial physics processes and constrain models of cosmic inflation. While traditional approaches utilise a limited set of statistical summaries of the galaxy distribution to constrain primordial non-Gaussiani… ▽ More Detecting and measuring a non-Gaussian signature of primordial origin in the density field is a major science goal of next-generation galaxy surveys. The signal will permit us to determine primordial physics processes and constrain models of cosmic inflation. While traditional approaches utilise a limited set of statistical summaries of the galaxy distribution to constrain primordial non-Gaussianity, we present a field-level approach by Bayesian forward-modelling the entire three-dimensional galaxy survey. Our method naturally and fully self-consistently exploits the entirety of the large-scale structure, e.g., higher-order statistics, peculiar velocity fields, and scale-dependent galaxy bias, to extract information on the local non-Gaussianity parameter, $\fnl$. We demonstrate the performance of our approach through various tests with mock galaxy data emulating relevant features of the \sdssiii{}-like survey, and additional tests with a \textit{Stage IV} mock data set. These tests reveal that the method infers unbiased values of $\fnl$ by accurately handling survey geometries, noise, and unknown galaxy biases. We demonstrate that our method can achieve constraints of $σ_{\fnl} \approx 8.78$ for \sdssiii{}-like data, an improvement of a factor $\sim 2.5$ over currently published constraints. Tests with next-generation mock data show that significant further improvements are feasible with sufficiently high resolution. Furthermore, the results demonstrate that our method can consistently marginalise all nuisance parameters of the data model. The method further provides an inference of the three-dimensional primordial density field, providing opportunities to explore additional signatures of primordial physics. △ Less

Submitted 16 March, 2022; originally announced March 2022.

Comments: 17 pages, 16 figures. Submitted to MNRAS. Comments are welcome

BAO scale inference from biased tracers using the EFT likelihood

Authors: Ivana Babić, Fabian Schmidt, Beatriz Tucci

Abstract: The physical scale corresponding to baryon acoustic oscillations (BAO), the size of the sound horizon at recombination, is precisely determined by CMB experiments. Measuring the apparent size of the BAO scale imprinted in the clustering of galaxies gives us a direct estimate of the angular-diameter distance and the Hubble parameter as a function of redshift. The BAO feature is damped by non-linear… ▽ More The physical scale corresponding to baryon acoustic oscillations (BAO), the size of the sound horizon at recombination, is precisely determined by CMB experiments. Measuring the apparent size of the BAO scale imprinted in the clustering of galaxies gives us a direct estimate of the angular-diameter distance and the Hubble parameter as a function of redshift. The BAO feature is damped by non-linear structure formation, which reduces the precision with which we can infer the BAO scale from standard galaxy clustering analysis methods. Many methods to undo this damping via the so-called BAO reconstruction have so far been proposed; however, they all rely on backward modeling. In this paper, we present the first results of BAO inference from rest-frame halo catalogs using forward modeling combined with the EFT likelihood, in the case where the initial phases of the density field are fixed. We show that the remaining systematic bias is less than 2% when we consider cutoff values of $Λ\leq 0.25 \,h\,{\rm Mpc}^{-1}$ for all halo samples considered, and below 1% and consistent with zero for all but the most highly biased samples. We also demonstrate that, when compared to the standard power spectrum likelihood approach under the same assumption of fixed phases, the 1$σ$ errors associated to the field level inference of the BAO scale are 1.1 to 3.3 times smaller, depending on the value of the cutoff and the halo sample. Our analysis therefore unveils another promising feature of using field-level inference for high-precision cosmology. △ Less

Submitted 26 August, 2022; v1 submitted 11 March, 2022; originally announced March 2022.

Comments: 22 pages, 12 figures

Journal ref: Journal of Cosmology and Astroparticle Physics, Volume 2022, Issue 08, id.007, 25 pp

Field-level inference of galaxy intrinsic alignment from the SDSS-III BOSS survey

Authors: Eleni Tsaprazi, Nhat-Minh Nguyen, Jens Jasche, Fabian Schmidt, Guilhem Lavaux

Abstract: As a large-scale overdensity collapses, it affects the orientation and shape of galaxies that form, by exerting tidal shear along their axes. Therefore, the shapes of elliptical galaxies align with the tidal field of cosmic structures. This intrinsic alignment provides insights into galaxy formation and the primordial universe, complements late-time cosmological probes and constitutes a significan… ▽ More As a large-scale overdensity collapses, it affects the orientation and shape of galaxies that form, by exerting tidal shear along their axes. Therefore, the shapes of elliptical galaxies align with the tidal field of cosmic structures. This intrinsic alignment provides insights into galaxy formation and the primordial universe, complements late-time cosmological probes and constitutes a significant systematic effect for weak gravitational lensing observations. In the present study, we provide constraints on the linear alignment model using a fully Bayesian field-level approach, using galaxy shape measurements from the SDSS-III BOSS LOWZ sample and three-dimensional tidal fields constrained with the LOWZ and CMASS galaxy samples of the SDSS-III BOSS survey. We find 4$σ$ evidence of intrinsic alignment, with an amplitude of $A_I=2.9 \pm 0.7$ at 20$h^{-1}\;\mathrm{Mpc}$. △ Less

Submitted 29 June, 2022; v1 submitted 8 December, 2021; originally announced December 2021.

Comments: 26 pages, 8 figures, accepted by JCAP

Prospects for Mapping Temporal Height Variations of the Seasonal CO2 Snow/Ice Caps at the Martian Poles by Co-registration of MOLA Profiles

Authors: Haifeng Xiao, Alexander Stark, Gregor Steinbrügge, Robin Thor, Frédéric Schmidt, Jürgen Oberst

Abstract: We investigate the feasibility and demonstrate the merits of using Mars Orbiter Laser Altimeter (MOLA) profiles to retrieve seasonal height variations of CO2 snow/ice cap in Mars' polar areas by applying a co-registration strategy. We present a prototype analysis on the research region of [85.75°S, 86.25°S, 300°E, 330°E] that is located on the residual south polar cap. Our method comprises the rec… ▽ More We investigate the feasibility and demonstrate the merits of using Mars Orbiter Laser Altimeter (MOLA) profiles to retrieve seasonal height variations of CO2 snow/ice cap in Mars' polar areas by applying a co-registration strategy. We present a prototype analysis on the research region of [85.75°S, 86.25°S, 300°E, 330°E] that is located on the residual south polar cap. Our method comprises the recomputation of MOLA footprint coordinates with an updated Mars Global Surveyor (MGS) ephemeris and a revised Mars rotation model. The reprocessed MOLA dataset at the South Pole of Mars (poleward of 78°S) is then self-registered to form a coherent reference digital terrain model (DTM). We co-register segments of reprocessed MOLA profiles to the self-registered MOLA reference DTM to obtain the temporal height differences at either footprints or cross-overs. Subsequently, a two-step Regional Pseudo Cross-over Adjustment (RPCA) procedure is proposed and applied to post-correct the aforementioned temporal height differences for a temporal systematic bias and other residual errors. These pseudo cross-overs are formed by profile pairs that do not necessarily intersect, but are connected through the underlaying DTM. Finally, CO2 snow/ice temporal height variation is obtained by median-filtering those post-corrected temporal height differences. The precision of the derived height change time series is ~4.9 cm. The peak-to-peak height variation is estimated to be ~2 m. In addition, a pronounced "pit" (transient height accumulation) of ~0.5 m in magnitude centered at Ls=210° in southern spring is observed. The proposed method opens the possibility to map the seasonal CO2 snow/ice height variations at the entire North and South polar regions of Mars. △ Less

Submitted 10 September, 2021; originally announced September 2021.

Stellar streams and dark substructure: the diffusion regime

Authors: M. Sten Delos, Fabian Schmidt

Abstract: The cold dark matter picture predicts an abundance of substructure within the Galactic halo. However, most substructures host no stars and can only be detected indirectly. Stellar streams present a promising probe of this dark substructure. These streams arise from tidally stripped star clusters or dwarf galaxies, and their low dynamical temperature and negligible self-gravity give them a sharp me… ▽ More The cold dark matter picture predicts an abundance of substructure within the Galactic halo. However, most substructures host no stars and can only be detected indirectly. Stellar streams present a promising probe of this dark substructure. These streams arise from tidally stripped star clusters or dwarf galaxies, and their low dynamical temperature and negligible self-gravity give them a sharp memory of gravitational perturbations caused by passing dark substructures. For this reason, perturbed stellar streams have been the subject of substantial study. While previous studies have been largely numerical, we show here that in the diffusion regime -- where stream stars are subjected to many small velocity kicks -- stream perturbations can be understood on a fully analytic level. In particular, we derive how the (three-dimensional) power spectrum of the substructure density field determines the power spectrum of the (one-dimensional) density of a stellar stream. Our analytic description supplies a clear picture of the behaviour of stream perturbations in response to a perturbing environment, which may include contributions from both dark and luminous substructure. In particular, stream perturbations grow in amplitude initially, settle into a steady state, and ultimately decay. By directly relating stellar stream perturbations to the surrounding matter distribution, this analytic framework represents a versatile new tool for probing the nature of dark matter through astrophysical observations. △ Less

Submitted 16 May, 2022; v1 submitted 30 August, 2021; originally announced August 2021.

Comments: 27 pages, 11 figures; matches published version. Arguments have been significantly clarified (results are unchanged); includes new discussion of the impact of time-varying substructure properties

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 513, Issue 3, July 2022, Pages 3682-3708

Covariant Decomposition of The Non-Linear Galaxy Number Counts and Their Monopole

Authors: Yonadav Barry Ginat, Vincent Desjacques, Donghui Jeong, Fabian Schmidt

Abstract: We present a fully non-linear and relativistically covariant expression for the observed galaxy density contrast. Building on a null tetrad tailored to the cosmological observer's past light cone, we find a decomposition of the non-linear galaxy over-density into manifestly gauge-invariant quantities, each of which has a clear physical interpretation as a cosmological observable. This ensures that… ▽ More We present a fully non-linear and relativistically covariant expression for the observed galaxy density contrast. Building on a null tetrad tailored to the cosmological observer's past light cone, we find a decomposition of the non-linear galaxy over-density into manifestly gauge-invariant quantities, each of which has a clear physical interpretation as a cosmological observable. This ensures that the monopole of the galaxy over-density field (the mean galaxy density as a function of observed redshift) is properly accounted for. We anticipate that this decomposition will be useful for future work on non-linearities in galaxy number counts, for example, deriving the relativistic expression for the galaxy bispectrum. We then specialise our results to conformal Newtonian gauge, with a Hubble parameter either defined globally or measured locally, illustrating the significance of the different contributions to the observed monopole of the galaxy density. △ Less

Submitted 21 February, 2022; v1 submitted 29 June, 2021; originally announced June 2021.

Comments: (v1) Comments welcome. (v2) Updated to agree with version accepted to JCAP; typos fixed

Assembly bias in quadratic bias parameters of dark matter halos from forward modeling

Authors: Titouan Lazeyras, Alexandre Barreira, Fabian Schmidt

Abstract: We use the forward modeling approach to galaxy clustering combined with the likelihood from the effective-field theory of large-scale structure to measure assembly bias, i.e. the dependence of halo bias on properties beyond the total mass, in the linear ($b_1$) and second order bias parameters ($b_2$ and $b_{K^2}$) of dark matter halos in $N$-body simulations. This is the first time that assembly… ▽ More We use the forward modeling approach to galaxy clustering combined with the likelihood from the effective-field theory of large-scale structure to measure assembly bias, i.e. the dependence of halo bias on properties beyond the total mass, in the linear ($b_1$) and second order bias parameters ($b_2$ and $b_{K^2}$) of dark matter halos in $N$-body simulations. This is the first time that assembly bias in the tidal bias parameter $b_{K^2}$ is measured. We focus on three standard halo properties: the concentration $c$, spin $λ$, and sphericity $s$, for which we find an assembly bias signal in $b_{K^2}$ that is opposite to that in $b_1$. Specifically, at fixed mass, halos that get more (less) positively biased in $b_1$, get less (more) negatively biased in $b_{K^2}$. We also investigate the impact of assembly bias on the $b_2(b_1)$ and $b_{K^2}(b_1)$ relations, and find that while the $b_2(b_1)$ relation stays roughly unchanged, assembly bias strongly impacts the $b_{K^2}(b_1)$ relation. This impact likely extends also to the corresponding relation for galaxies, which motivates future studies to design better priors on $b_{K^2}(b_1)$ for use in cosmological constraints from galaxy clustering data. △ Less

Submitted 8 November, 2021; v1 submitted 28 June, 2021; originally announced June 2021.

Comments: 22 pages, 5 figures ; v2: version published in JCAP

Journal ref: JCAP10(2021)063

Galaxy bias from forward models: linear and second-order bias of IllustrisTNG galaxies

Authors: Alexandre Barreira, Titouan Lazeyras, Fabian Schmidt

Abstract: We use field-level forward models of galaxy clustering and the EFT likelihood formalism to study, for the first time for self-consistently simulated galaxies, the relations between the linear $b_1$ and second-order bias parameters $b_2$ and $b_{K^2}$. The forward models utilize all of the information available in the galaxy distribution up to a given order in perturbation theory, which allows us t… ▽ More We use field-level forward models of galaxy clustering and the EFT likelihood formalism to study, for the first time for self-consistently simulated galaxies, the relations between the linear $b_1$ and second-order bias parameters $b_2$ and $b_{K^2}$. The forward models utilize all of the information available in the galaxy distribution up to a given order in perturbation theory, which allows us to infer these bias parameters with high signal-to-noise, even from relatively small volumes ($L_{\rm box} = 205{\rm Mpc}/h$). We consider galaxies from the IllustrisTNG simulations, and our main result is that the $b_2(b_1)$ and $b_{K^2}(b_1)$ relations obtained from gravity-only simulations for total mass selected objects are broadly preserved for simulated galaxies selected by stellar mass, star formation rate, color and black hole accretion rate. We also find good agreement between the bias relations of the simulated galaxies and a number of recent estimates for observed galaxy samples. The consistency under different galaxy selection criteria suggests that theoretical priors on these bias relations may be used to improve cosmological constraints based on observed galaxy samples. We do identify some small differences between the bias relations in the hydrodynamical and gravity-only simulations, which we show can be linked to the environmental dependence of the relation between galaxy properties and mass. We also show that the EFT likelihood recovers the value of $σ_8$ to percent-level from various galaxy samples (including splits by color and star formation rate) and after marginalizing over 8 bias parameters. This demonstration using simulated galaxies adds to previous works based on halos as tracers, and strengthens further the potential of forward models to infer cosmology from galaxy data. △ Less

Submitted 14 July, 2021; v1 submitted 6 May, 2021; originally announced May 2021.

Comments: 19 pages, 9 figures. Version 2 matches published version and contains a comparison to bias estimates from observations. Comments welcomed!

An $n$-th order Lagrangian Forward Model for Large-Scale Structure

Authors: Fabian Schmidt

Abstract: A forward model of matter and biased tracers at arbitrary order in Lagrangian perturbation theory (LPT) is presented. The forward model contains the complete LPT displacement field at any given order in perturbations, as well as all relevant bias operators at that order and leading order in derivatives. The construction is done for any expansion history and does not rely on the Einstein-de Sitter… ▽ More A forward model of matter and biased tracers at arbitrary order in Lagrangian perturbation theory (LPT) is presented. The forward model contains the complete LPT displacement field at any given order in perturbations, as well as all relevant bias operators at that order and leading order in derivatives. The construction is done for any expansion history and does not rely on the Einstein-de Sitter approximation. A large subset of higher-derivative bias operators is also included. As validation test, we compare the $n$LPT-predicted matter density field and that from N-body simulations using the same initial conditions. For simulations using a cutoff in the initial conditions, we find subpercent agreement up to scales of $k\sim 0.2 h\,{\rm Mpc}^{-1}$. We also find subpercent agreement with full simulations without cutoff, both for the power spectrum and nonlinear $σ_8$-inference, when allowing for the effective sound speed. The application to biased tracers (halos) has already been presented in a recent paper (arXiv:2009.14176). △ Less

Submitted 1 October, 2021; v1 submitted 17 December, 2020; originally announced December 2020.

Comments: 24 pages, 5 figures; v2: minor corrections, clarifications, added references, reflects JCAP published version; code available at https://gitlab.mpcdf.mpg.de/leftfield/release/leftfield

Journal ref: JCAP 04 (2021), 033

Machine Learning for automatic identification of new minor species

Authors: Frederic Schmidt, Guillaume Cruz Mermy, Justin Erwin, Severine Robert, Lori Neary, Ian R. Thomas, Frank Daerden, Bojan Ristic, Manish R. Patel, Giancarlo Bellucci, Jose-Juan Lopez-Moreno, Ann-Carine Vandaele

Abstract: One of the main difficulties to analyze modern spectroscopic datasets is due to the large amount of data. For example, in atmospheric transmittance spectroscopy, the solar occultation channel (SO) of the NOMAD instrument onboard the ESA ExoMars2016 satellite called Trace Gas Orbiter (TGO) had produced $\sim$10 millions of spectra in 20000 acquisition sequences since the beginning of the mission in… ▽ More One of the main difficulties to analyze modern spectroscopic datasets is due to the large amount of data. For example, in atmospheric transmittance spectroscopy, the solar occultation channel (SO) of the NOMAD instrument onboard the ESA ExoMars2016 satellite called Trace Gas Orbiter (TGO) had produced $\sim$10 millions of spectra in 20000 acquisition sequences since the beginning of the mission in April 2018 until 15 January 2020. Other datasets are even larger with $\sim$billions of spectra for OMEGA onboard Mars Express or CRISM onboard Mars Reconnaissance Orbiter. Usually, new lines are discovered after a long iterative process of model fitting and manual residual analysis. Here we propose a new method based on unsupervised machine learning, to automatically detect new minor species. Although precise quantification is out of scope, this tool can also be used to quickly summarize the dataset, by giving few endmembers ("source") and their abundances. We approximate the dataset non-linearity by a linear mixture of abundance and source spectra (endmembers). We used unsupervised source separation in form of non-negative matrix factorization to estimate those quantities. Several methods are tested on synthetic and simulation data. Our approach is dedicated to detect minor species spectra rather than precisely quantifying them. On synthetic example, this approach is able to detect chemical compounds present in form of 100 hidden spectra out of $10^4$, at 1.5 times the noise level. Results on simulated spectra of NOMAD-SO targeting CH$_{4}$ show that detection limits goes in the range of 100-500 ppt in favorable conditions. Results on real martian data from NOMAD-SO show that CO$_{2}$ and H$_{2}$O are present, as expected, but CH$_{4}$ is absent. Nevertheless, we confirm a set of new unexpected lines in the database, attributed by ACS instrument Team to the CO$_{2}$ magnetic dipole. △ Less

Submitted 15 December, 2020; originally announced December 2020.

Comments: 26 pages, 10 figures

Journal ref: Quantitative Spectroscopy and Radiative Transfer, 2021, 259, 107361

Galaxy shape statistics in the effective field theory

Authors: Zvonimir Vlah, Nora Elisa Chisari, Fabian Schmidt

Abstract: Intrinsic galaxy alignments yield an important contribution to the observed statistics of galaxy shapes. The general bias expansion for galaxy sizes and shapes in three dimensions has been recently described by Vlah, Chisari \& Schmidt using the general perturbative effective field theory (EFT) framework, in analogy to the clustering of galaxies. In this work, we present a formalism that uses the… ▽ More Intrinsic galaxy alignments yield an important contribution to the observed statistics of galaxy shapes. The general bias expansion for galaxy sizes and shapes in three dimensions has been recently described by Vlah, Chisari \& Schmidt using the general perturbative effective field theory (EFT) framework, in analogy to the clustering of galaxies. In this work, we present a formalism that uses the properties of spherical tensors to project galaxy shapes onto the observed sky in the flat-sky approximation and compute the two-point functions at next-to-leading order as well as the leading-order three-point functions of galaxy shapes and number counts. The resulting expressions are given in forms that are convenient for efficient numerical implementation. For a source redshift distribution typical of Stage IV surveys, we find that that nonlinear intrinsic alignment contributions to galaxy shape correlations become relevant at angular wavenumbers $l \gtrsim 100$. △ Less

Submitted 7 December, 2020; originally announced December 2020.

Comments: 32 pages, 5 figures

Neutrino mass constraints beyond linear order: cosmology dependence and systematic biases

Authors: Aoife Boyle, Fabian Schmidt

Abstract: We demonstrate the impact on forecasted neutrino mass constraints of extending galaxy clustering and CMB lensing predictions from linear to next-to-leading-order power spectra. The redshift-space 1-loop power spectrum model we adopt requires an additional four free bias parameters, a velocity bias parameter and two new stochastic parameters. These additional nuisance parameters appreciably weaken… ▽ More We demonstrate the impact on forecasted neutrino mass constraints of extending galaxy clustering and CMB lensing predictions from linear to next-to-leading-order power spectra. The redshift-space 1-loop power spectrum model we adopt requires an additional four free bias parameters, a velocity bias parameter and two new stochastic parameters. These additional nuisance parameters appreciably weaken the constraints on $M_ν$. CMB lensing plays a significant role in helping to alleviate these degeneracies and tighten the final constraints. The constraint on the optical depth to reionisation $τ$ has a strong effect on the constraint on $M_ν$, but only when CMB lensing is included in the analysis to keep the degeneracies with the nuisance parameters under control. We also extract constraints when 1) using the BAO signature only as a distance probe, and 2) isolating the scale-dependence of the power spectrum, which, as shown in previous work, provides a cosmology-independent probe of $M_ν$. All constraints except the latter remain strongly sensitive to the assumption of a flat $Λ$CDM universe. We perform an analysis of the magnitude of the shift introduced in the inferred $M_ν$ value when neglecting nonlinear corrections, and show that, for a Euclid-like survey, this shift becomes roughly equal to the 1$σ$ constraint itself even with a conservative cut-off scale of $k_{max} = 0.1~h~{\rm Mpc}^{-1}$. We also perform a calculation of the appropriate expected bias in neutrino mass caused by not including the next, 2-loop order and expect a shift of only about 20% of the 1$σ$ error for $k_{max}=0.2~h~{\rm Mpc}^{-1}$ in a Euclid-like survey. △ Less

Submitted 3 June, 2021; v1 submitted 20 November, 2020; originally announced November 2020.

Comments: Prepared for submission to JCAP

Journal ref: JCAP 04 (2021) 022

Impacts of the physical data model on the forward inference of initial conditions from biased tracers

Authors: Nhat-Minh Nguyen, Fabian Schmidt, Guilhem Lavaux, Jens Jasche

Abstract: We investigate the impact of each ingredient in the employed physical data model on the Bayesian forward inference of initial conditions from biased tracers at the field level. Specifically, we use dark matter halos in a given cosmological simulation volume as tracers of the underlying matter density field. We study the effect of tracer density, grid resolution, gravity model, bias model and likel… ▽ More We investigate the impact of each ingredient in the employed physical data model on the Bayesian forward inference of initial conditions from biased tracers at the field level. Specifically, we use dark matter halos in a given cosmological simulation volume as tracers of the underlying matter density field. We study the effect of tracer density, grid resolution, gravity model, bias model and likelihood on the inferred initial conditions. We find that the cross-correlation coefficient between true and inferred phases reacts weakly to all ingredients above, and is well predicted by the theoretical expectation derived from a Gaussian model on a broad range of scales. The bias in the amplitude of the inferred initial conditions, on the other hand, depends strongly on the bias model and the likelihood. We conclude that the bias model and likelihood hold the key to an unbiased cosmological inference. Together they must keep the systematics -- which arise from the sub-grid physics that are marginalized over -- under control in order to obtain an unbiased inference. △ Less

Submitted 16 March, 2021; v1 submitted 12 November, 2020; originally announced November 2020.

Comments: 29 pages, 11 figures, 1 table; main results in Sec 4.2; comments welcome. v2: Match JCAP accepted version (no major change), update references. v3: Fix links, typos and further update reference info. v4: Match JCAP published version, fix an incorrect link and line spacing

Journal ref: JCAP03(2021)058

On the sensitivity of weak gravitational lensing to the cosmic expansion function

Authors: Christian F. Schmidt, Matthias Bartelmann

Abstract: We analyse the functional derivative of the cosmic-shear power spectrum $C_\ell^γ$ with respect to the cosmic expansion function. Our interest in doing so is two-fold: (i) In view of attempts to detect minor changes of the cosmic expansion function which may be due to a possibly time-dependent dark-energy density, we wish to know how sensitive the weak-lensing power spectrum is to changes in the e… ▽ More We analyse the functional derivative of the cosmic-shear power spectrum $C_\ell^γ$ with respect to the cosmic expansion function. Our interest in doing so is two-fold: (i) In view of attempts to detect minor changes of the cosmic expansion function which may be due to a possibly time-dependent dark-energy density, we wish to know how sensitive the weak-lensing power spectrum is to changes in the expansion function. (ii) In view of recent empirical determinations of the cosmic expansion function from distance measurements, independent of specific cosmological models, we wish to find out how uncertainties in the expansion function translate to uncertainties in the cosmic-shear power spectrum. We find the following answers: Relative changes of the expansion function are amplified by the cosmic-shear power spectrum by a factor $\approx 2-6$, weakly depending on the scale factor where the change is applied, and the current uncertainty of one example for an empirically determined expansion function translates to a relative uncertainty of the cosmic-shear power spectrum of $\approx10\,\%$. △ Less

Submitted 21 December, 2023; v1 submitted 6 November, 2020; originally announced November 2020.

Comments: 8 pages, 8 figures, submitted to MNRAS

Sigma-Eight at the Percent Level: The EFT Likelihood in Real Space

Authors: Fabian Schmidt

Abstract: The effective field theory likelihood for the density field of biased tracers allows for cosmology inference from the clustering of galaxies that consistently uses all available information at a given order in perturbation theory. This paper presents results and implementation details on the real-space (as opposed to Fourier-space) formulation of the likelihood, which allows for the incorporation… ▽ More The effective field theory likelihood for the density field of biased tracers allows for cosmology inference from the clustering of galaxies that consistently uses all available information at a given order in perturbation theory. This paper presents results and implementation details on the real-space (as opposed to Fourier-space) formulation of the likelihood, which allows for the incorporation of survey window functions. The implementation further uses a Lagrangian forward model for biased tracers which automatically accounts for all relevant contributions up to any desired order. Unbiased inference of $σ_8$ is demonstrated at the 2% level for cutoff values $Λ\leq 0.14h\,{\rm Mpc}^{-1}$ for halo samples over a range of masses and redshifts. The inferred value shows the expected convergence to the ground truth in the low-cutoff limit. Apart from the possibility of including observational effects, this represents further substantial improvement over previous results based on the EFT likelihood. △ Less

Submitted 22 April, 2021; v1 submitted 29 September, 2020; originally announced September 2020.

Comments: 31 pages, 10 figures; results summarized in Fig. 5; v2: corrected and expanded discussion of cosmic variance; v3: minor edits and added references, reflects JCAP published version

Journal ref: JCAP 04 (2021), 032

Constraints on effusive cryovolcanic eruptions on Europa using topography obtained from Galileo images

Authors: Elodie Lesage, Frédéric Schmidt, François Andrieu, Hélène Massol

Abstract: Images of Europa's surface taken by the Galileo Solid State Imager (SSI) show smooth features measuring a few kilometers, potentially resulting from eruptions of low-viscosity material such as liquid cryomagma. We estimated the volume of four of these smooth features by producing digital elevation models (DEMs) of four Galileo/SSI images. We used the shape-from-shading technique with special care… ▽ More Images of Europa's surface taken by the Galileo Solid State Imager (SSI) show smooth features measuring a few kilometers, potentially resulting from eruptions of low-viscosity material such as liquid cryomagma. We estimated the volume of four of these smooth features by producing digital elevation models (DEMs) of four Galileo/SSI images. We used the shape-from-shading technique with special care to estimate the uncertainties on the produced DEMs and estimated feature volumes to be between ($5.7*10^{7}$ m$^{3}$ and ($2.7*10^{8}$ m$^{3}$. We discussed the implications for putative sub-surface liquid reservoir dimensions in the case of eruptions induced from freezing reservoirs. Our previous cryovolcanic eruption model was improved by considering a cycle of cryomagma freezing and effusion and by estimating the vaporized cryolava fraction once cryolava spreads onto Europa's surface. Our results show that the cryomagma reservoirs would have to be relatively large to generate these smooth features (1 to 100 km$^{3}$ if the flow features result from a single eruption, and 0.4 to 60 km$^{3}$ for the full lifetime of a reservoir generating cyclic eruptions). The two future missions JUICE (ESA) and Europa Clipper (NASA) should reach Europa during the late 2020s. They shall give more information on those putative cryovolcanic regions which appear as interesting targets that could provide a better understanding of the material exchanges between the surface, sub-surface and ocean of Europa. △ Less

Submitted 24 February, 2021; v1 submitted 29 September, 2020; originally announced September 2020.

Journal ref: Icarus, 2021, 114373

Galaxy imaging surveys as spin-sensitive detector for cosmological colliders

Authors: Kazuhiro Kogai, Kazuyuki Akitsu, Fabian Schmidt, Yuko Urakawa

Abstract: Galaxy imaging surveys provide us with information on both the galaxy distribution and their shapes. In this paper, we systematically investigate the sensitivity of galaxy shapes to new physics in the initial conditions. For this purpose, we decompose the galaxy shape function into spin components, and compute the contributions to each spin component from both intrinsic alignment and weak lensing.… ▽ More Galaxy imaging surveys provide us with information on both the galaxy distribution and their shapes. In this paper, we systematically investigate the sensitivity of galaxy shapes to new physics in the initial conditions. For this purpose, we decompose the galaxy shape function into spin components, and compute the contributions to each spin component from both intrinsic alignment and weak lensing. We then consider the angular-dependent primordial non-Gaussianity, which is generated by a non-zero integer spin particle when active during inflation, and show that a galaxy imaging survey essentially functions as a spin-sensitive detector of such particles in the early universe. We also perform a forecast of the PNG generated from a higher spin particle, considering a Rubin Observatory LSST-like galaxy survey. △ Less

Submitted 18 March, 2021; v1 submitted 11 September, 2020; originally announced September 2020.

Comments: 50 pages, 8 figures

Journal ref: JCAP03(2021)060

Image processing for precise geometry determination

Authors: Ines Belgacem, Grégory Jonniaux, Frédéric Schmidt

Abstract: Reliable spatial information can be difficult to obtain in planetary remote sensing applications because of errors present in the metadata of images taken with space probes. We have designed a pipeline to address this problem on disk-resolved images of Jupiter's moon Europa taken with New Horizons' LOng Range Reconnaissance Imager, Galileo's Solid State Imager and Voyager's Imaging Science Subsyst… ▽ More Reliable spatial information can be difficult to obtain in planetary remote sensing applications because of errors present in the metadata of images taken with space probes. We have designed a pipeline to address this problem on disk-resolved images of Jupiter's moon Europa taken with New Horizons' LOng Range Reconnaissance Imager, Galileo's Solid State Imager and Voyager's Imaging Science Subsystem. We correct for errors in the spacecraft position, pointing and the target's attitude by comparing them to the same reference. We also address ways to correct for distortion prior to any metadata consideration. Finally, we propose a vectorized method to efficiently project images pixels onto an elliptic target and compute the coordinates and geometry of observation at each intercept point. △ Less

Submitted 4 September, 2020; originally announced September 2020.

Journal ref: PSS.2020.105081

Taking measurements of the kinematic Sunyaev-Zel'dovich effect forward: including uncertainties from velocity reconstruction with forward modeling

Authors: Nhat-Minh Nguyen, Jens Jasche, Guilhem Lavaux, Fabian Schmidt

Abstract: We measure the kinematic Sunyaev-Zel'dovich (kSZ) effect, imprinted by maxBCG clusters, on the Planck SMICA map of the Cosmic Microwave Background (CMB). Our measurement, for the first time, directly accounts for uncertainties in the velocity reconstruction step through the process of Bayesian forward modeling. We show that this often neglected uncertainty budget typically increases the final unce… ▽ More We measure the kinematic Sunyaev-Zel'dovich (kSZ) effect, imprinted by maxBCG clusters, on the Planck SMICA map of the Cosmic Microwave Background (CMB). Our measurement, for the first time, directly accounts for uncertainties in the velocity reconstruction step through the process of Bayesian forward modeling. We show that this often neglected uncertainty budget typically increases the final uncertainty on the measured kSZ signal amplitude by $\simeq15\%$ at cluster scale. We observe evidence for the kSZ effect, at a significance of $\simeq2σ$. Our analysis, when applied to future higher-resolution CMB data, together with minor improvements in map-filtering and signal-modeling methods, should yield both significant and unbiased measurements of the kSZ signal, which can then be used to probe and constrain baryonic content of galaxy clusters and galaxy groups. △ Less

Submitted 6 October, 2020; v1 submitted 27 July, 2020; originally announced July 2020.

Comments: 29 pages, 15 figures, 2 tables; main results in Figs. 4, 8, 10; comments welcome. v2: Add arguments and a simplified treatment for photo-z clusters (sec. 3.3), update numerical results accordingly, main results and conclusions remain unchanged. v3: Add statements regarding data and code availability, add citations, update acknowledgement, match JCAP version

Journal ref: JCAP12(2020)011

Regional study of Europa's photometry

Authors: Ines Belgacem, Frédéric Schmidt, Grégory Jonniaux

Abstract: The surface of Europa is geologically young and shows signs of current activity. Studying it from a photometric point of view gives us insight on its physical state. We used a collection of 57 images from Voyager's Imaging Science System and New Horizons' LOng Range Reconnaissance Imager for which we corrected the geometric metadata and projected every pixel to compute photometric information (ref… ▽ More The surface of Europa is geologically young and shows signs of current activity. Studying it from a photometric point of view gives us insight on its physical state. We used a collection of 57 images from Voyager's Imaging Science System and New Horizons' LOng Range Reconnaissance Imager for which we corrected the geometric metadata and projected every pixel to compute photometric information (reflectance and geometry of observation). We studied 20 areas scattered across the surface of Europa and estimated their photometric behavior using the Hapke radiative transfer model and a Bayesian framework in order to estimate their microphysical state. We have found that most of them were consistent with the bright backscattering behavior of Europa, already observed at a global scale, indicating the presence of grains maturated by space weathering. However, we have identified very bright areas showing a narrow forward scattering possibly indicating the presence of fresh deposits that could be attributed to recent cryovolcanism or jets. Overall, we showed that the photometry of Europa's surface is more diverse than previously thought and so is its microphysical state. △ Less

Submitted 21 July, 2020; originally announced July 2020.

Comments: 46 pages, 17 figures, 3 tables

Report number: ISSN 0019-1035

Journal ref: Icarus, Volume 338, 2020, 113525, ISSN 0019-1035

Information content in the redshift-space galaxy power spectrum and bispectrum

Authors: Nishant Agarwal, Vincent Desjacques, Donghui Jeong, Fabian Schmidt

Abstract: We present a Fisher information study of the statistical impact of galaxy bias and selection effects on the estimation of key cosmological parameters from galaxy redshift surveys; in particular, the angular diameter distance, Hubble parameter, and linear growth rate at a given redshift, the cold dark matter density, and the tilt and running of the primordial power spectrum. The line-of-sight-depen… ▽ More We present a Fisher information study of the statistical impact of galaxy bias and selection effects on the estimation of key cosmological parameters from galaxy redshift surveys; in particular, the angular diameter distance, Hubble parameter, and linear growth rate at a given redshift, the cold dark matter density, and the tilt and running of the primordial power spectrum. The line-of-sight-dependent selection contributions we include here are known to exist in real galaxy samples. We determine the maximum wavenumber included in the analysis by requiring that the next-order corrections to the galaxy power spectrum or bispectrum, treated here at next-to-leading and leading order, respectively, produce shifts of $\lesssim 0.25σ$ on each of the six cosmological parameters. With the galaxy power spectrum alone, selection effects can deteriorate the constraints severely, especially on the linear growth rate. Adding the galaxy bispectrum helps break parameter degeneracies significantly. We find that a joint power spectrum-bispectrum analysis of a Euclid-like survey can still measure the linear growth rate to 10% precision after complete marginalization over selection bias. We also discuss systematic parameter shifts arising from ignoring selection effects and/or other bias parameters, and emphasize that it is necessary to either control selection effects at the percent level or marginalize over them. We obtain similar results for the Roman Space Telescope and HETDEX. △ Less

Submitted 5 March, 2021; v1 submitted 8 July, 2020; originally announced July 2020.

Comments: 42 pages, including 3 appendices, 4 tables, 12 figures. Fisher code available at: https://bitbucket.org/nishant_agarwal/fisher_lss_p-b/src/master/. v2: Minor revisions, added discussion in sections 4.3.2 and 6, matches published version

Journal ref: JCAP 03 (2021) 021

Galaxy bias and primordial non-Gaussianity: insights from galaxy formation simulations with IllustrisTNG

Authors: Alexandre Barreira, Giovanni Cabass, Fabian Schmidt, Annalisa Pillepich, Dylan Nelson

Abstract: We study the impact that large-scale perturbations of (i) the matter density and (ii) the primordial gravitational potential with local primordial non-Gaussianity (PNG) have on galaxy formation using the IllustrisTNG model. We focus on the linear galaxy bias $b_1$ and the coefficient $b_φ$ of the scale-dependent bias induced by PNG, which describe the response of galaxy number counts to these two… ▽ More We study the impact that large-scale perturbations of (i) the matter density and (ii) the primordial gravitational potential with local primordial non-Gaussianity (PNG) have on galaxy formation using the IllustrisTNG model. We focus on the linear galaxy bias $b_1$ and the coefficient $b_φ$ of the scale-dependent bias induced by PNG, which describe the response of galaxy number counts to these two types of perturbations, respectively. We perform our study using separate universe simulations, in which the effect of the perturbations is mimicked by changes to the cosmological parameters: modified cosmic matter density for $b_1$ and modified amplitude $\mathcal{A}_s$ of the primordial scalar power spectrum for $b_φ$. We find that the widely used universality relation $b_φ= 2δ_c(b_1 - 1)$ is a poor description of the bias of haloes and galaxies selected by stellar mass $M_*$, which is instead described better by $b_φ(M_*) = 2δ_c(b_1(M_*) - p)$ with $p \in [0.4, 0.7]$. This is explained by the different impact that matter overdensities and local PNG have on the median stellar-to-halo-mass relation. A simple model of this impact allows us to describe the stellar mass dependence of $b_1$ and $b_φ$ fairly well. Our results also show a nontrivial relation between $b_1$ and $b_φ$ for galaxies selected by color and black hole mass accretion rate. Our results provide refined priors on $b_φ$ for local PNG constraints and forecasts using galaxy clustering. Given that the widely used universality relation underpredicts $b_φ(M_*)$, existing analyses may underestimate the true constraining power on local PNG. △ Less

Submitted 3 December, 2020; v1 submitted 16 June, 2020; originally announced June 2020.

Comments: 21 pages, 12 figures, 1 table. Comments welcomed! v2: version published in JCAP

Unbiased Cosmology Inference from Biased Tracers using the EFT Likelihood

Authors: Fabian Schmidt, Giovanni Cabass, Jens Jasche, Guilhem Lavaux

Abstract: We present updates on the cosmology inference using the effective field theory (EFT) likelihood presented previously in Schmidt et al., 2018, Elsner et al., 2019 [1,2]. Specifically, we add a cutoff to the initial conditions that serve as starting point for the matter forward model. We show that this cutoff, which was not employed in any previous related work, is important to regularize loop integ… ▽ More We present updates on the cosmology inference using the effective field theory (EFT) likelihood presented previously in Schmidt et al., 2018, Elsner et al., 2019 [1,2]. Specifically, we add a cutoff to the initial conditions that serve as starting point for the matter forward model. We show that this cutoff, which was not employed in any previous related work, is important to regularize loop integrals that otherwise involve small-scale, non-perturbative modes. We then present results on the inferred value of the linear power spectrum normalization $σ_{8}$ from rest-frame halo catalogs using both second- and third-order bias expansions, imposing uniform priors on all bias parameters. Due to the perfect bias-$σ_{8}$ degeneracy at linear order, constraints on $σ_{8}$ rely entirely on nonlinear information. The results show the expected convergence behavior when lowering the cutoff in wavenumber, $Λ$. When including modes up to $k \leq Λ= 0.1\,h\,{\rm Mpc}^{-1}$ in the second-order case, $σ_{8}$ is recovered to within $\lesssim 6\,\%$ for a range of halo masses and redshifts. The systematic bias shrinks to $4\,\%$ or less for the third-order bias expansion on the same range of scales. Together with additional evidence we provide, this shows that the residual mismatch in $σ_{8}$ can be attributed to higher-order bias contributions. We conclude that the EFT likelihood is able to infer unbiased cosmological constraints, within expected theoretical systematic errors, from physical biased tracers on quasilinear scales △ Less

Submitted 5 November, 2020; v1 submitted 14 April, 2020; originally announced April 2020.

Comments: 36 pages, 10 figures, 2 tables; main results in Figs. 2, 4, 6; v2: Sec. 3 clarified, Sec. 5.4 expanded; matches published JCAP version

Journal ref: JCAP 11 (2020), 008

The Likelihood for LSS: Stochasticity of Bias Coefficients at All Orders

Authors: Giovanni Cabass, Fabian Schmidt

Abstract: In the EFT of biased tracers the noise field $\varepsilon_g$ is not exactly uncorrelated with the nonlinear matter field $δ$. Its correlation with $δ$ is effectively captured by adding stochasticities to each bias coefficient. We show that if these stochastic fields are Gaussian (the impact of their non-Gaussianity being subleading on quasi-linear scales anyway) it is possible to resum exactly the… ▽ More In the EFT of biased tracers the noise field $\varepsilon_g$ is not exactly uncorrelated with the nonlinear matter field $δ$. Its correlation with $δ$ is effectively captured by adding stochasticities to each bias coefficient. We show that if these stochastic fields are Gaussian (the impact of their non-Gaussianity being subleading on quasi-linear scales anyway) it is possible to resum exactly their effect on the conditional likelihood ${\cal P}[δ_g|δ]$ to observe a galaxy field $δ_g$ given an underlying $δ$. This resummation allows to take them into account in EFT-based approaches to Bayesian forward modeling. We stress that the resulting corrections to a purely Gaussian conditional likelihood with white-noise covariance are the most relevant on scales where the EFT is under control: they are more important than any non-Gaussianity of the noise $\varepsilon_g$. △ Less

Submitted 23 July, 2020; v1 submitted 1 April, 2020; originally announced April 2020.

Comments: 15 pages. v2: very minor change in Eqs. (2.6), (2.8), (2.17), (2.19), (3.16), (3.18), (3.19) to show factorization in real space of the normalization of the likelihood, updated acknowledgements. v3: added section on renormalization, matches published version

Measuring the Tidal Response of Structure Formation: Anisotropic Separate Universe Simulations using TreePM

Authors: Jens Stücker, Andreas. S. Schmidt, Simon D. M. White, Fabian Schmidt, Oliver Hahn

Abstract: We present anisotropic "separate universe" simulations which modify the N-body code Gadget-4 in order to represent a large-scale tidal field through an anisotropic expansion factor. These simulations are used to measure the linear, quasi-linear and nonlinear response of the matter power spectrum to a spatially uniform trace-free tidal field up to wavenumber $k = 7 h \text{Mpc}^{-1}$. Together with… ▽ More We present anisotropic "separate universe" simulations which modify the N-body code Gadget-4 in order to represent a large-scale tidal field through an anisotropic expansion factor. These simulations are used to measure the linear, quasi-linear and nonlinear response of the matter power spectrum to a spatially uniform trace-free tidal field up to wavenumber $k = 7 h \text{Mpc}^{-1}$. Together with the response to a large-scale overdensity measured in previous work, this completely describes the nonlinear matter bispectrum in the squeezed limit. We find that the response amplitude does not approach zero on small scales in physical coordinates, but rather a constant value at $z=0$, $R_K\approx 0.5$ for $k \geq 3 h\text{Mpc}^{-1}$ up to the scale where we consider our simulations reliable, $k \leq 7 h\text{Mpc}^{-1}$ at $z=0$. This shows that even the inner regions of haloes are affected by the large-scale tidal field. We also measure directly the alignment of halo shapes with the tidal field, finding a clear signal which increases with halo mass. △ Less

Submitted 13 November, 2020; v1 submitted 13 March, 2020; originally announced March 2020.

Comments: 19 pages, 9 figures

Silicate grain growth due to ion trapping in oxygen-rich supernova remnants like Cassiopeia A

Authors: Florian Kirchschlager, M. J. Barlow, Franziska D. Schmidt

Abstract: Core-collapse supernovae can condense large masses of dust post-explosion. However, sputtering and grain-grain collisions during the subsequent passage of the dust through the reverse shock can potentially destroy a significant fraction of the newly formed dust before it can reach the interstellar medium. Here we show that in oxygen-rich supernova remnants like Cassiopeia A the penetration and tra… ▽ More Core-collapse supernovae can condense large masses of dust post-explosion. However, sputtering and grain-grain collisions during the subsequent passage of the dust through the reverse shock can potentially destroy a significant fraction of the newly formed dust before it can reach the interstellar medium. Here we show that in oxygen-rich supernova remnants like Cassiopeia A the penetration and trapping within silicate grains of the same impinging ions of oxygen, silicon and magnesium that are responsible for grain surface sputtering can significantly reduce the net loss of grain material. We model conditions representative of dusty clumps (density contrast $χ=100$) passing through the reverse shock in the oxygen-rich Cassiopeia A remnant and find that, compared to cases where the effect is neglected, as well as facilitating the formation of grains larger than those that had originally condensed, ion trapping increases the surviving masses of silicate dust by factors of up to two to four, depending on initial grain radii. For higher density contrasts ($χ\gtrsim180$), we find that the effect of gas accretion on the surface of dust grains surpasses ion trapping, and the survival rate increases to ${\sim}55 \%$ of the initial dust mass for $χ=256$. △ Less

Submitted 10 March, 2020; v1 submitted 6 March, 2020; originally announced March 2020.

Comments: Accepted by ApJ. Author accepted manuscript. Accepted on 06/03/2020. Deposited on 06/03/2020. 11 pages

Compensated Isocurvature Perturbations in the Galaxy Power Spectrum

Authors: Alexandre Barreira, Giovanni Cabass, Kaloian D. Lozanov, Fabian Schmidt

Abstract: We investigate the potential of the galaxy power spectrum to constrain compensated isocurvature perturbations (CIPs), primordial fluctuations in the baryon density that are compensated by fluctuations in CDM density to ensure an unperturbed total matter density. We show that CIPs contribute to the galaxy overdensity at linear order, and if they are close to scale-invariant, their effects are nearl… ▽ More We investigate the potential of the galaxy power spectrum to constrain compensated isocurvature perturbations (CIPs), primordial fluctuations in the baryon density that are compensated by fluctuations in CDM density to ensure an unperturbed total matter density. We show that CIPs contribute to the galaxy overdensity at linear order, and if they are close to scale-invariant, their effects are nearly perfectly degenerate with the local PNG parameter $f_{\rm nl}$ if they correlate with the adiabatic perturbations. This degeneracy can however be broken by analyzing multiple galaxy samples with different bias parameters, or by taking CMB priors on $f_{\rm nl}$ into account. Parametrizing the amplitude of the CIP power spectrum as $P_{σσ} = A^2P_{\mathcal{R}\mathcal{R}}$ (where $P_{\mathcal{R}\mathcal{R}}$ is the adiabatic power spectrum) we find, for a number of fiducial galaxy samples in a simplified forecast setup, that constraints on $A$, relative to those on $f_{\rm nl}$, of order $σ_{A}/σ_{f_{\rm nl}} \approx 1-2$ are achievable for CIPs correlated with adiabatic perturbations, and $σ_{A}/σ_{f_{\rm nl}} \approx 5$ for the uncorrelated case. These values are independent of survey volume, and suggest that current galaxy data are already able to improve significantly on the tightest existing constraints on CIPs from the CMB. Future galaxy surveys that aim to achieve $σ_{f_{\rm nl}} \sim 1$ have the potential to place even stronger bounds on CIPs. △ Less

Submitted 21 September, 2020; v1 submitted 28 February, 2020; originally announced February 2020.

Comments: 17 pages + references and appendix; 6 figures; 1 table. Comments welcomed! v2 matches version published in JCAP