Showing 1–50 of 410 results for author: Bressan, A

Two Nitsche-based mixed finite element discretizations for the seepage problem in Richards' equation

Authors: Federico Gatti, Andrea Bressan, Alessio Fumagalli, Domenico Gallipoli, Leonardo Maria Lalicata, Simone Pittaluga, Lorenzo Tamellini

Abstract: This paper proposes two algorithms to impose seepage boundary conditions in the context of Richards' equation for groundwater flows in unsaturated media. Seepage conditions are non-linear boundary conditions, that can be formulated as a set of unilateral constraints on both the pressure head and the water flux at the ground surface, together with a complementarity condition: these conditions in pr… ▽ More This paper proposes two algorithms to impose seepage boundary conditions in the context of Richards' equation for groundwater flows in unsaturated media. Seepage conditions are non-linear boundary conditions, that can be formulated as a set of unilateral constraints on both the pressure head and the water flux at the ground surface, together with a complementarity condition: these conditions in practice require switching between Neumann and Dirichlet boundary conditions on unknown portions on the boundary. Upon realizing the similarities of these conditions with unilateral contact problems in mechanics, we take inspiration from that literature to propose two approaches: the first method relies on a strongly consistent penalization term, whereas the second one is obtained by an hybridization approach, in which the value of the pressure on the surface is treated as a separate set of unknowns. The flow problem is discretized in mixed form with div-conforming elements so that the water mass is preserved. Numerical experiments show the validity of the proposed strategy in handling the seepage boundary conditions on geometries with increasing complexity. △ Less

Submitted 10 July, 2024; originally announced July 2024.

The PLATO Mission

Authors: Heike Rauer, Conny Aerts, Juan Cabrera, Magali Deleuil, Anders Erikson, Laurent Gizon, Mariejo Goupil, Ana Heras, Jose Lorenzo-Alvarez, Filippo Marliani, Cesar Martin-Garcia, J. Miguel Mas-Hesse, Laurence O'Rourke, Hugh Osborn, Isabella Pagano, Giampaolo Piotto, Don Pollacco, Roberto Ragazzoni, Gavin Ramsay, Stéphane Udry, Thierry Appourchaux, Willy Benz, Alexis Brandeker, Manuel Güdel, Eduardo Janot-Pacheco , et al. (801 additional authors not shown)

Abstract: PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observati… ▽ More PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observations from the ground, planets will be characterised for their radius, mass, and age with high accuracy (5 %, 10 %, 10 % for an Earth-Sun combination respectively). PLATO will provide us with a large-scale catalogue of well-characterised small planets up to intermediate orbital periods, relevant for a meaningful comparison to planet formation theories and to better understand planet evolution. It will make possible comparative exoplanetology to place our Solar System planets in a broader context. In parallel, PLATO will study (host) stars using asteroseismology, allowing us to determine the stellar properties with high accuracy, substantially enhancing our knowledge of stellar structure and evolution. The payload instrument consists of 26 cameras with 12cm aperture each. For at least four years, the mission will perform high-precision photometric measurements. Here we review the science objectives, present PLATO's target samples and fields, provide an overview of expected core science performance as well as a description of the instrument and the mission profile at the beginning of the serial production of the flight cameras. PLATO is scheduled for a launch date end 2026. This overview therefore provides a summary of the mission to the community in preparation of the upcoming operational phases. △ Less

Submitted 8 June, 2024; originally announced June 2024.

Intermediate Domains for Scalar Conservation Laws

Authors: Fabio Ancona, Alberto Bressan, Elio Marconi, Luca Talamini

Abstract: For a scalar conservation law with strictly convex flux, by Oleinik's estimates the total variation of a solution with initial data $\overline{u}\in \bf{L}^\infty(\mathbb R)$ decays like $t^{-1}$. This paper introduces a class of intermediate domains $\mathcal P_α$, $0<α<1$, such that for $\overline u\in \mathcal P_α$ a faster decay rate is achieved:… ▽ More For a scalar conservation law with strictly convex flux, by Oleinik's estimates the total variation of a solution with initial data $\overline{u}\in \bf{L}^\infty(\mathbb R)$ decays like $t^{-1}$. This paper introduces a class of intermediate domains $\mathcal P_α$, $0<α<1$, such that for $\overline u\in \mathcal P_α$ a faster decay rate is achieved: $\mathrm{Tot.Var.}\bigl\{ u(t,\cdot)\bigr\}\sim t^{α-1}$. A key ingredient of the analysis is a ``Fourier-type" decomposition of $\overline u$ into components which oscillate more and more rapidly. The results aim at extending the theory of fractional domains for analytic semigroups to an entirely nonlinear setting. △ Less

Submitted 16 April, 2024; originally announced April 2024.

MSC Class: 35L65

Discovery of a dormant 33 solar-mass black hole in pre-release Gaia astrometry

Authors: Gaia Collaboration, P. Panuzzo, T. Mazeh, F. Arenou, B. Holl, E. Caffau, A. Jorissen, C. Babusiaux, P. Gavras, J. Sahlmann, U. Bastian, Ł. Wyrzykowski, L. Eyer, N. Leclerc, N. Bauchet, A. Bombrun, N. Mowlavi, G. M. Seabroke, D. Teyssier, E. Balbinot, A. Helmi, A. G. A. Brown, A. Vallenari, T. Prusti, J. H. J. de Bruijne , et al. (390 additional authors not shown)

Abstract: Gravitational waves from black-hole merging events have revealed a population of extra-galactic BHs residing in short-period binaries with masses that are higher than expected based on most stellar evolution models - and also higher than known stellar-origin black holes in our Galaxy. It has been proposed that those high-mass BHs are the remnants of massive metal-poor stars. Gaia astrometry is exp… ▽ More Gravitational waves from black-hole merging events have revealed a population of extra-galactic BHs residing in short-period binaries with masses that are higher than expected based on most stellar evolution models - and also higher than known stellar-origin black holes in our Galaxy. It has been proposed that those high-mass BHs are the remnants of massive metal-poor stars. Gaia astrometry is expected to uncover many Galactic wide-binary systems containing dormant BHs, which may not have been detected before. The study of this population will provide new information on the BH-mass distribution in binaries and shed light on their formation mechanisms and progenitors. As part of the validation efforts in preparation for the fourth Gaia data release (DR4), we analysed the preliminary astrometric binary solutions, obtained by the Gaia Non-Single Star pipeline, to verify their significance and to minimise false-detection rates in high-mass-function orbital solutions. The astrometric binary solution of one source, Gaia BH3, implies the presence of a 32.70 \pm 0.82 M\odot BH in a binary system with a period of 11.6 yr. Gaia radial velocities independently validate the astrometric orbit. Broad-band photometric and spectroscopic data show that the visible component is an old, very metal-poor giant of the Galactic halo, at a distance of 590 pc. The BH in the Gaia BH3 system is more massive than any other Galactic stellar-origin BH known thus far. The low metallicity of the star companion supports the scenario that metal-poor massive stars are progenitors of the high-mass BHs detected by gravitational-wave telescopes. The Galactic orbit of the system and its metallicity indicate that it might belong to the Sequoia halo substructure. Alternatively, and more plausibly, it could belong to the ED-2 stream, which likely originated from a globular cluster that had been disrupted by the Milky Way. △ Less

Submitted 19 April, 2024; v1 submitted 16 April, 2024; originally announced April 2024.

Comments: 23 pages, accepted fro publication in A&A Letters. New version with small fixes

Generic Properties of Conjugate Points in Optimal Control Problems

Authors: Alberto Bressan, Marco Mazzola, Khai T. Nguyen

Abstract: The first part of the paper studies a class of optimal control problems in Bolza form, where the dynamics is linear w.r.t.~the control function. A necessary condition is derived, for the optimality of a trajectory which starts at a conjugate point. The second part is concerned with a classical problem in the Calculus of Variations, with free terminal point. For a generic terminal cost… ▽ More The first part of the paper studies a class of optimal control problems in Bolza form, where the dynamics is linear w.r.t.~the control function. A necessary condition is derived, for the optimality of a trajectory which starts at a conjugate point. The second part is concerned with a classical problem in the Calculus of Variations, with free terminal point. For a generic terminal cost $ψ\in \C^4(\mathbb{R}^n)$, applying the previous necessary condition we show that the set of conjugate points is contained in the image of an $(n-2)$-dimensional manifold, and has locally bounded $(n-2)$-dimensional Hausdorff measure. △ Less

Submitted 2 April, 2024; originally announced April 2024.

Comments: 13 pages

MSC Class: 49K05; 49L12

Parallelization in time by diagonalization

Authors: Andrea Bressan, Alen Kushova, Gabriele Loli, Monica Montardini, Giancarlo Sangalli, Mattia Tani

Abstract: This is a review of preconditioning techniques based on fast-diagonalization methods for space-time isogeometric discretization of the heat equation. Three formulation are considered: the Galerkin approach, a discrete least-square and a continuous least square. For each formulation the heat differential operator is written as a sum of terms that are kronecker products of uni-variate operators. The… ▽ More This is a review of preconditioning techniques based on fast-diagonalization methods for space-time isogeometric discretization of the heat equation. Three formulation are considered: the Galerkin approach, a discrete least-square and a continuous least square. For each formulation the heat differential operator is written as a sum of terms that are kronecker products of uni-variate operators. These are used to speed-up the application of the operator in iterative solvers and to construct a suitable preconditioner. Contrary to the fast-diagonalization technique for the Laplace equation where all uni-variate operators acting on the same direction can be simultaneously diagonalized in the case of the heat equation this is not possible. Luckily this can be done up to an additional term that has low rank allowing for the utilization of arrow-head like factorization or inversion by Sherman-Morrison formula. The proposed preconditioners work extremely well on the parametric domain and, when the domain is parametrized or when the equation coefficients are not constant, they can be adapted and retain good performance characteristics. △ Less

Submitted 30 November, 2023; originally announced March 2024.

Comments: arXiv admin note: substantial text overlap with arXiv:1909.07309, arXiv:2311.18461

MSC Class: 65

Constraining the Initial Mass function in the Epoch of Reionization from Astrophysical and Cosmological data

Authors: A. Lapi, G. Gandolfi, L. Boco, F. Gabrielli, M. Massardi, B. S. Haridasu, C. Baccigalupi, A. Bressan, L. Danese

Abstract: [abridged] We aim to constrain the stellar initial mass function (IMF) during the epoch of reionization. To this purpose, we build up a semi-empirical model for the reionization history of the Universe, based on various ingredients: the latest determination of the UV galaxy luminosity function from JWST out to redshift $z\lesssim 12$; data-inferred and simulation-driven assumptions on the redshift… ▽ More [abridged] We aim to constrain the stellar initial mass function (IMF) during the epoch of reionization. To this purpose, we build up a semi-empirical model for the reionization history of the Universe, based on various ingredients: the latest determination of the UV galaxy luminosity function from JWST out to redshift $z\lesssim 12$; data-inferred and simulation-driven assumptions on the redshift-dependent escape fraction of ionizing photons from primordial galaxies; a simple yet flexible parameterization of the IMF $φ(m_\star)\sim m_\star^ξ\, e^{-m_{\star,\rm c}/m_\star}$ in terms of a high-mass end slope $ξ<0$ and of a characteristic mass $m_{\star,\rm c}$ below which a flattening or a bending sets in; the PARSEC stellar evolution code to compute the UV and ionizing emission from different star's masses as a function of age and metallicity; a few physical constraints related to stellar and galaxy formation in faint galaxies at the reionization redshifts. We compare our model outcomes with the reionization observables from different astrophysical and cosmological probes, and perform Bayesian inference on the IMF parameters. We find that the IMF slope $ξ$ is within the range from $-2.8$ to $-2.3$, while appreciably flatter slopes are excluded at great significance. However, the bestfit value of the IMF characteristic mass $m_{\star,\rm c}\sim$ a few $M_\odot$ implies a suppression in the formation of small stellar masses, at variance with the IMF in the local Universe; this may be induced by the thermal background $\sim 20-30$ K provided by CMB photons at the reionization redshifts. Finally, we investigate the implications of our reconstructed IMF on the recent JWST detections of massive galaxies at and beyond the reionization epoch, showing that any putative tension with the standard cosmological framework is substantially alleviated. △ Less

Submitted 13 March, 2024; v1 submitted 12 March, 2024; originally announced March 2024.

Comments: 21 pages, 11 figures, typos corrected, in press on Universe

On the Contribution of Very Massive Stars to the Sulfur Abundance in Star-Forming Galaxies: the Role of PISN

Authors: S. Goswami, J. M. Vilchez, B. Perez-Diaz, L. Silva, A. Bressan, E. Perez-Montero

Abstract: Recent work presented increasing evidence of high, non-constant S/O abundance ratios observed in star-forming metal-poor galaxies, showing deviations from the constant canonical S/O across a large range of O/H abundance. Similar peculiar high Fe/O ratios have been also recently detected. We investigate whether these high S/O ratios at low metallicities could be explained taking into consideration… ▽ More Recent work presented increasing evidence of high, non-constant S/O abundance ratios observed in star-forming metal-poor galaxies, showing deviations from the constant canonical S/O across a large range of O/H abundance. Similar peculiar high Fe/O ratios have been also recently detected. We investigate whether these high S/O ratios at low metallicities could be explained taking into consideration the process of Pair Instability Supernovae (PISN) in chemical modelling through which similar behaviour observed for Fe/O ratios was successfully reproduced. We use chemical evolution models which take into account the stages of PISN in the yields published by Goswami et al. 2022, and adopt a suitable initial mass function (IMF) to characterize this evolutionary stage .appropriately. The peculiar high values and the behaviour of the observed S/O versus O/H relation can be reproduced when the ejecta of very massive stars that go through the process of PISN are taken into account. Additionally, a bi-modal top-heavy IMF and an initial strong burst of star formation are required to attain the reported high S/O values. We show that the role of very massive stars going through the process of PISN should be taken into account when explaining the chemical enrichment of sulfur and oxygen in metal-poor star-forming regions. △ Less

Submitted 20 February, 2024; originally announced February 2024.

Comments: 8 pages, 3 figures, accepted for publication in A&A

Journal ref: A&A 685, A81 (2024)

GalaPy, the highly optimised C++/Python spectral modelling tool for galaxies -- I. Library presentation and photometric fitting

Authors: Tommaso Ronconi, Andrea Lapi, Martina Torsello, Alessandro Bressan, Darko Donevski, Lara Pantoni, Meriem Behiri, Lumen Boco, Andrea Cimatti, Quirino D'Amato, Luigi Danese, Marika Giulietti, Francesca Perrotta, Laura Silva, Margherita Talia, Marcella Massardi

Abstract: Fostered by upcoming data from new generation observational campaigns, we are about to enter a new era for the study of how galaxies form and evolve. The unprecedented quantity of data that will be collected, from distances only marginally grasped up to now, will require analysis tools designed to target the specific physical peculiarities of the observed sources and handle extremely large dataset… ▽ More Fostered by upcoming data from new generation observational campaigns, we are about to enter a new era for the study of how galaxies form and evolve. The unprecedented quantity of data that will be collected, from distances only marginally grasped up to now, will require analysis tools designed to target the specific physical peculiarities of the observed sources and handle extremely large datasets. One powerful method to investigate the complex astrophysical processes that govern the properties of galaxies is to model their observed spectral energy distribution (SED) at different stages of evolution and times throughout the history of the Universe. To address these challenges, we have developed GalaPy, a new library for modelling and fitting SEDs of galaxies from the X-ray to the radio band, as well as the evolution of their components and dust attenuation/reradiation. GalaPy incorporates both empirical and physically-motivated star formation histories, state-of-the-art single stellar population synthesis libraries, a two-component dust model for attenuation, an age-dependent energy conservation algorithm to compute dust reradiation, and additional sources of stellar continuum such as synchrotron, nebular/free-free emission and X-ray radiation from low and high mass binary stars. GalaPy has a hybrid implementation that combines the high performance of compiled C++ with the flexibility of Python, and exploits an object-oriented design. It generates models on the fly without relying on templates, and exploits fully Bayesian parameter space sampling. In this first work, we introduce the project and showcase the photometric SED fitting tools already available to users. The library is available on the Python Package Index (PyPI) and comes with extensive online documentation and tutorials. △ Less

Submitted 19 February, 2024; originally announced February 2024.

Comments: 41 pages, 31 figures, 7 tables, to be published on A&A, links to documentation and PyPI available in the PDF, comments are very welcome

Journal ref: A&A 685, A161 (2024)

The JWST Resolved Stellar Populations Early Release Science Program V. DOLPHOT Stellar Photometry for NIRCam and NIRISS

Authors: Daniel R. Weisz, Andrew E. Dolphin, Alessandro Savino, Kristen B. W. McQuinn, Max J. B. Newman, Benjamin F. Williams, Nitya Kallivayalil, Jay Anderson, Martha L. Boyer, Matteo Correnti, Marla C. Geha, Karin M. Sandstrom, Andrew A. Cole, Jack T. Warfield, Evan D. Skillman, Roger E. Cohen, Rachael Beaton, Alessandro Bressan, Alberto Bolatto, Michael Boylan-Kolchin, Alyson M. Brooks, James S. Bullock, Charlie Conroy, Michael C. Cooper, Julianne J. Dalcanton , et al. (16 additional authors not shown)

Abstract: We present NIRCam and NIRISS modules for DOLPHOT, a widely-used crowded field stellar photometry package. We describe details of the modules including pixel masking, astrometric alignment, star finding, photometry, catalog creation, and artificial star tests (ASTs). We tested these modules using NIRCam and NIRISS images of M92 (a Milky Way globular cluster), Draco II (an ultra-faint dwarf galaxy),… ▽ More We present NIRCam and NIRISS modules for DOLPHOT, a widely-used crowded field stellar photometry package. We describe details of the modules including pixel masking, astrometric alignment, star finding, photometry, catalog creation, and artificial star tests (ASTs). We tested these modules using NIRCam and NIRISS images of M92 (a Milky Way globular cluster), Draco II (an ultra-faint dwarf galaxy), and WLM (a star-forming dwarf galaxy). DOLPHOT's photometry is highly precise and the color-magnitude diagrams are deeper and have better definition than anticipated during original program design in 2017. The primary systematic uncertainties in DOLPHOT's photometry arise from mismatches in the model and observed point spread functions (PSFs) and aperture corrections, each contributing $\lesssim0.01$ mag to the photometric error budget. Version 1.2 of WebbPSF models, which include charge diffusion and interpixel capacitance effects, significantly reduced PSF-related uncertainties. We also observed minor ($\lesssim0.05$ mag) chip-to-chip variations in NIRCam's zero points, which will be addressed by the JWST flux calibration program. Globular cluster observations are crucial for photometric calibration. Temporal variations in the photometry are generally $\lesssim0.01$ mag, although rare large misalignment events can introduce errors up to 0.08 mag. We provide recommended DOLPHOT parameters, guidelines for photometric reduction, and advice for improved observing strategies. Our ERS DOLPHOT data products are available on MAST, complemented by comprehensive online documentation and tutorials for using DOLPHOT with JWST imaging data. △ Less

Submitted 5 February, 2024; originally announced February 2024.

Comments: 30 pages, 17 figures. Accepted to ApJS. Data products to be hosted on MAST. For DOLPHOT/JWST tutorials, see https://dolphot-jwst.readthedocs.io/en/latest/ . For more program and DOLPHOT info, see https://ers-stars.github.io

The Impact of the Third Dredge-up and Mass Loss in Shaping the Initial-Final Mass Relation of White Dwarfs

Authors: Francesco Addari, Paola Marigo, Alessandro Bressan, Guglielmo Costa, Kendall Shepherd, Guglielmo Volpato

Abstract: The initial-final mass relation (IFMR) plays a crucial role in understanding stellar structure and evolution by linking a star's initial mass to the mass of the resulting white dwarf. This study explores the IFMR in the initial mass range $0.8 \leq M_\mathrm{ini} / M_\odot \leq 4$ using full PARSEC evolutionary calculations supplemented with COLIBRI computations to complete the ejection of the env… ▽ More The initial-final mass relation (IFMR) plays a crucial role in understanding stellar structure and evolution by linking a star's initial mass to the mass of the resulting white dwarf. This study explores the IFMR in the initial mass range $0.8 \leq M_\mathrm{ini} / M_\odot \leq 4$ using full PARSEC evolutionary calculations supplemented with COLIBRI computations to complete the ejection of the envelope and obtain the final core mass. Recent works have shown that the supposed monotonicity of the IFMR is interrupted by a kink in the initial mass range $M_\mathrm{ini} \approx 1.65-2.10 M_\odot$, due to the interaction between recurrent dredge-up episodes and stellar winds in carbon stars evolving on the thermally-pulsing asymptotic giant branch phase. To reproduce the IFMR non-monotonic behavior we investigate the role of convective overshooting efficiency applied to the base of the convective envelope ($f_\mathrm{env}$) and to the borders of the pulse-driven convective zone ($f_\mathrm{pdcz}$), as well as its interplay with mass loss. We compare our models to observational data and find that $f_\mathrm{env}$ must vary with initial mass in order to accurately reproduce the IFMR's observed kink and slopes. We find some degeneracy between the overshooting parameters when only the IFMR information is used. Nonetheless, this analysis provides valuable insights into the internal mixing processes during the TP-AGB phase. △ Less

Submitted 18 January, 2024; originally announced January 2024.

Comments: Accepted for publication in ApJ. 17 pages, 12 figures

High-statistics measurement of Collins and Sivers asymmetries for transversely polarised deuterons

Authors: G. D. Alexeev, M. G. Alexeev, C. Alice, A. Amoroso, V. Andrieux, V. Anosov, S. Asatryan, K. Augsten, W. Augustyniak, C. D. R. Azevedo, B. Badelek, J. Barth, R. Beck, J. Beckers, Y. Bedfer, J. Bernhard, M. Bodlak, F. Bradamante, A. Bressan, W. -C. Chang, C. Chatterjee, M. Chiosso, A. G. Chumakov, S. -U. Chung, A. Cicuttin , et al. (162 additional authors not shown)

Abstract: New results are presented on a high-statistics measurement of Collins and Sivers asymmetries of charged hadrons produced in deep inelastic scattering of muons on a transversely polarised $^6$LiD target. The data were taken in 2022 with the COMPASS spectrometer using the 160 \gevv\ muon beam at CERN, balancing the existing data on transversely polarised proton targets. The first results from about… ▽ More New results are presented on a high-statistics measurement of Collins and Sivers asymmetries of charged hadrons produced in deep inelastic scattering of muons on a transversely polarised $^6$LiD target. The data were taken in 2022 with the COMPASS spectrometer using the 160 \gevv\ muon beam at CERN, balancing the existing data on transversely polarised proton targets. The first results from about two-thirds of the new data have total uncertainties smaller by up to a factor of three compared to the previous deuteron measurements. Using all the COMPASS proton and deuteron results, both the transversity and the Sivers distribution functions of the $u$ and $d$ quark, as well as the tensor charge in the measured $x$-range are extracted. In particular, the accuracy of the $d$ quark results is significantly improved. △ Less

Submitted 30 December, 2023; originally announced January 2024.

Report number: CERN-EP-2023-308

Final COMPASS results on the transverse-spin-dependent azimuthal asymmetries in the pion-induced Drell-Yan process

Authors: G. D. Alexeev, M. G. Alexeev, C. Alice, A. Amoroso, V. Andrieux, V. Anosov, K. Augsten, W. Augustyniak, C. D. R. Azevedo, B. Badelek, J. Barth, R. Beck, J. Beckers, Y. Bedfer, J. Bernhard, M. Bodlak, F. Bradamante, A. Bressan, W. -C. Chang, C. Chatterjee, M. Chiosso, A. G. Chumakov, S. -U. Chung, A. Cicuttin, P. M. M. Correia , et al. (159 additional authors not shown)

Abstract: The COMPASS Collaboration performed measurements of the Drell-Yan process in 2015 and 2018 using a 190 GeV/c $π^{-}$ beam impinging on a transversely polarised ammonia target. Combining the data of both years, we present final results on the amplitudes of the five azimuthal modulations in the dimuon production cross section. Three of these transverse-spin-dependent azimuthal asymmetries (TSAs) pro… ▽ More The COMPASS Collaboration performed measurements of the Drell-Yan process in 2015 and 2018 using a 190 GeV/c $π^{-}$ beam impinging on a transversely polarised ammonia target. Combining the data of both years, we present final results on the amplitudes of the five azimuthal modulations in the dimuon production cross section. Three of these transverse-spin-dependent azimuthal asymmetries (TSAs) probe the nucleon leading-twist Sivers, transversity, and pretzelosity transverse-momentum dependent (TMD) parton distribution functions (PDFs). The other two are induced by subleading effects. These TSAs provide unique new inputs for the study of the nucleon TMD PDFs and their universality properties. In particular, the Sivers TSA observed in this measurement is consistent with the fundamental QCD prediction of a sign change of naive time-reversal-odd TMD PDFs when comparing the Drell-Yan process with semi-inclusive measurements of deep inelastic scattering. Also, within the context of model predictions, the observed transversity TSA is consistent with the expectation of a sign change for the Boer-Mulders function. △ Less

Submitted 28 December, 2023; originally announced December 2023.

Report number: CERN-EP-2023-307

Nebular emission from young stellar populations including binary stars

Authors: M. Lecroq, S. Charlot, A. Bressan, G. Bruzual, G. Costa, G. Iorio, M. Spera, M. Mapelli, Y. Chen, J. Chevallard, M. Dall'Amico

Abstract: We investigate the nebular emission produced by young stellar populations using the new GALSEVN model based on the combination of the SEVN population-synthesis code including binary-star processes and the GALAXEV code for the spectral evolution of stellar populations. Photoionization calculations performed with the CLOUDY code confirm that accounting for binary-star processes strongly influences t… ▽ More We investigate the nebular emission produced by young stellar populations using the new GALSEVN model based on the combination of the SEVN population-synthesis code including binary-star processes and the GALAXEV code for the spectral evolution of stellar populations. Photoionization calculations performed with the CLOUDY code confirm that accounting for binary-star processes strongly influences the predicted emission-line properties of young galaxies. In particular, we find that our model naturally reproduces the strong HeII/Hb ratios commonly observed at high Hb equivalent widths in metal-poor, actively star-forming galaxies, which have proven challenging to reproduce using previous models. Including bursty star formation histories broadens the agreement with observations, while the most extreme HeII equivalent widths can be reproduced by models dominated by massive stars. GALSEVN also enables us to compute, for the first time in a way physically consistent with stellar emission, the emission from accretion discs of X-ray binaries (XRBs) and radiative shocks driven by stellar winds and supernova explosions. We find that these contributions are unlikely to prominently affect the predicted HeII/Hb ratio, and that previous claims of a significant contribution by XRBs to the luminosities of high-ionization lines are based on models predicting improbably high ratios of X-ray luminosity to star formation rate, inconsistent with the observed average luminosity function of XRBs in nearby galaxies. The results presented here provide a solid basis for a more comprehensive investigation of the physical properties of observed galaxies with GALSEVN using Bayesian inference. △ Less

Submitted 13 December, 2023; originally announced December 2023.

Comments: 26 pages, 18 figures. Accepted for publication in MNRAS

Energy conservative isogeometric techniques for the wave equation

Authors: Andrea Bressan, Annalisa Buffa, Alen Kushova, Rafael Vázquez

Abstract: We analyze the wave equation in mixed form, with periodic and/or Dirichlet homogeneous boundary conditions, and nonconstant coefficients that depend on the spatial variable. For the discretization, the weak form of the second equation is replaced by a strong form, written in terms of a projection operator. The system of equations is discretized with B-splines forming a De Rham complex along with s… ▽ More We analyze the wave equation in mixed form, with periodic and/or Dirichlet homogeneous boundary conditions, and nonconstant coefficients that depend on the spatial variable. For the discretization, the weak form of the second equation is replaced by a strong form, written in terms of a projection operator. The system of equations is discretized with B-splines forming a De Rham complex along with suitable commutative projectors for the approximation of the second equation. The discrete scheme is energy conservative when discretized in time with a conservative method such as Crank-Nicolson. We propose a convergence analysis of the method to study the dependence with respect to the mesh size $h$, with focus on the consistency error. Numerical results show optimal convergence of the error in energy norm, and a relative error in energy conservation for long-time simulations of the order of machine precision. △ Less

Submitted 30 November, 2023; originally announced November 2023.

MSC Class: 65

Space-time least squares approximation for Schrödinger equation and efficient solver

Authors: Andrea Bressan, Alen Kushova, Giancarlo Sangalli, Mattia Tani

Abstract: In this work we present a space-time least squares isogeometric discretization of the Schrödinger equation and propose a preconditioner for the arising linear system in the parametric domain. Exploiting the tensor product structure of the basis functions, the preconditioner is written as the sum of Kronecker products of matrices. Thanks to an extension to the classical Fast Diagonalization method,… ▽ More In this work we present a space-time least squares isogeometric discretization of the Schrödinger equation and propose a preconditioner for the arising linear system in the parametric domain. Exploiting the tensor product structure of the basis functions, the preconditioner is written as the sum of Kronecker products of matrices. Thanks to an extension to the classical Fast Diagonalization method, the application of the preconditioner is efficient and robust w.r.t. the polynomial degree of the spline space. The time required for the application is almost proportional to the number of degrees-of-freedom, for a serial execution. △ Less

Submitted 30 November, 2023; originally announced November 2023.

Comments: arXiv admin note: text overlap with arXiv:1909.07309

MSC Class: 65

A Study of Primordial Very Massive Star Evolution. II. Stellar Rotation and Gamma-Ray Burst Progenitors

Authors: Guglielmo Volpato, Paola Marigo, Guglielmo Costa, Alessandro Bressan, Michele Trabucchi, Léo Girardi, Francesco Addari

Abstract: We calculate new evolutionary models of rotating primordial very massive stars, with initial mass from $100\,M_{\odot}$ to $200\,M_{\odot}$, for two values of the initial metallicity ${Z=0}$ and ${Z=0.0002}$. For the first time in this mass range, we consider stellar rotation and pulsation-driven mass loss, along with radiative winds. The models evolve from the zero-age main sequence, until the on… ▽ More We calculate new evolutionary models of rotating primordial very massive stars, with initial mass from $100\,M_{\odot}$ to $200\,M_{\odot}$, for two values of the initial metallicity ${Z=0}$ and ${Z=0.0002}$. For the first time in this mass range, we consider stellar rotation and pulsation-driven mass loss, along with radiative winds. The models evolve from the zero-age main sequence, until the onset of pair instability. We discuss the main properties of the models during their evolution and then focus on the final fate and the possible progenitors of jet-driven events. All tracks that undergo pulsational-pair instability produce successful gamma-ray bursts (GRB) in the collapsar framework, while those that collapse directly to black holes (BH) produce jet-driven supernova events. In these latter cases, the expected black hole mass changes due to the jet propagation inside the progenitor, resulting in different models that should produce BH within the pair-instability black-hole mass gap. Successful GRBs predicted here from zero-metallicity and very metal-poor progenitors may be bright enough to be detected even up to redshift ${\sim20}$ using current telescopes such as the Swift-BAT X-ray detector and the JWST. △ Less

Submitted 19 January, 2024; v1 submitted 30 November, 2023; originally announced November 2023.

Comments: Published in ApJ, 23 pages, 12 figures

Journal ref: ApJ 961 89 (2024)

One Dimensional Hyperbolic Conservation Laws: Past and Future

Authors: Alberto Bressan

Abstract: Aim of these notes is provide a brief review of the current well-posedness theory for hyperbolic systems of conservation laws in one space dimension, also pointing out open problems and possible research directions. They supplement the slides of the short course given by the author in Erice, May 2023, available at: sites.google.com/view/erice23/speakers-and-slides. Aim of these notes is provide a brief review of the current well-posedness theory for hyperbolic systems of conservation laws in one space dimension, also pointing out open problems and possible research directions. They supplement the slides of the short course given by the author in Erice, May 2023, available at: sites.google.com/view/erice23/speakers-and-slides. △ Less

Submitted 25 October, 2023; originally announced October 2023.

Comments: 38 pages, 25 figures

MSC Class: 35L45; 35L65; 35L67

AESOPUS 2.0: Low-Temperature Opacities with Solid Grains

Authors: P. Marigo, P. Woitke, E. Tognelli, L. Girardi, B. Aringer, A. Bressan

Abstract: In this study we compute the equation of state and Rosseland mean opacity from temperatures of T~30000 K down to T~400 K, pushing the capabilities of the AESOPUS code (Marigo et al., 2022; Marigo & Aringer, 2009) into the regime where solid grains can form. The GGchem code (Woitke et al. 2018) is used to solve the chemistry for temperatures less than ~3000 K. Atoms, molecules, and dust grains in t… ▽ More In this study we compute the equation of state and Rosseland mean opacity from temperatures of T~30000 K down to T~400 K, pushing the capabilities of the AESOPUS code (Marigo et al., 2022; Marigo & Aringer, 2009) into the regime where solid grains can form. The GGchem code (Woitke et al. 2018) is used to solve the chemistry for temperatures less than ~3000 K. Atoms, molecules, and dust grains in thermodynamic equilibrium are all included in the equation of state. To incorporate monochromatic atomic and molecular cross sections, an optimized opacity sampling technique is used. The Mie theory is employed to calculate the opacity of 43 grain species. Tables of Rosseland mean opacities for scaled-solar compositions are provided. Based on our computing resources, opacities for other chemical patterns, as well as various grain sizes, porosity, and shapes, can be easily computed upon user request to the corresponding author. △ Less

Submitted 23 October, 2023; originally announced October 2023.

Comments: 14 pages, 8 figures, resubmitted to ApJ following moderate revision

Gaia Focused Product Release: Sources from Service Interface Function image analysis -- Half a million new sources in omega Centauri

Authors: Gaia Collaboration, K. Weingrill, A. Mints, J. Castañeda, Z. Kostrzewa-Rutkowska, M. Davidson, F. De Angeli, J. Hernández, F. Torra, M. Ramos-Lerate, C. Babusiaux, M. Biermann, C. Crowley, D. W. Evans, L. Lindegren, J. M. Martín-Fleitas, L. Palaversa, D. Ruz Mieres, K. Tisanić, A. G. A. Brown, A. Vallenari, T. Prusti, J. H. J. de Bruijne, F. Arenou, A. Barbier , et al. (378 additional authors not shown)

Abstract: Gaia's readout window strategy is challenged by very dense fields in the sky. Therefore, in addition to standard Gaia observations, full Sky Mapper (SM) images were recorded for nine selected regions in the sky. A new software pipeline exploits these Service Interface Function (SIF) images of crowded fields (CFs), making use of the availability of the full two-dimensional (2D) information. This ne… ▽ More Gaia's readout window strategy is challenged by very dense fields in the sky. Therefore, in addition to standard Gaia observations, full Sky Mapper (SM) images were recorded for nine selected regions in the sky. A new software pipeline exploits these Service Interface Function (SIF) images of crowded fields (CFs), making use of the availability of the full two-dimensional (2D) information. This new pipeline produced half a million additional Gaia sources in the region of the omega Centauri ($ω$ Cen) cluster, which are published with this Focused Product Release. We discuss the dedicated SIF CF data reduction pipeline, validate its data products, and introduce their Gaia archive table. Our aim is to improve the completeness of the {\it Gaia} source inventory in a very dense region in the sky, $ω$ Cen. An adapted version of {\it Gaia}'s Source Detection and Image Parameter Determination software located sources in the 2D SIF CF images. We validated the results by comparing them to the public {\it Gaia} DR3 catalogue and external Hubble Space Telescope data. With this Focused Product Release, 526\,587 new sources have been added to the {\it Gaia} catalogue in $ω$ Cen. Apart from positions and brightnesses, the additional catalogue contains parallaxes and proper motions, but no meaningful colour information. While SIF CF source parameters generally have a lower precision than nominal {\it Gaia} sources, in the cluster centre they increase the depth of the combined catalogue by three magnitudes and improve the source density by a factor of ten. This first SIF CF data publication already adds great value to the {\it Gaia} catalogue. It demonstrates what to expect for the fourth {\it Gaia} catalogue, which will contain additional sources for all nine SIF CF regions. △ Less

Submitted 8 November, 2023; v1 submitted 10 October, 2023; originally announced October 2023.

Journal ref: A&A 680, A35 (2023)

Gaia Focused Product Release: A catalogue of sources around quasars to search for strongly lensed quasars

Authors: Gaia Collaboration, A. Krone-Martins, C. Ducourant, L. Galluccio, L. Delchambre, I. Oreshina-Slezak, R. Teixeira, J. Braine, J. -F. Le Campion, F. Mignard, W. Roux, A. Blazere, L. Pegoraro, A. G. A. Brown, A. Vallenari, T. Prusti, J. H. J. de Bruijne, F. Arenou, C. Babusiaux, A. Barbier, M. Biermann, O. L. Creevey, D. W. Evans, L. Eyer, R. Guerra , et al. (376 additional authors not shown)

Abstract: Context. Strongly lensed quasars are fundamental sources for cosmology. The Gaia space mission covers the entire sky with the unprecedented resolution of $0.18$" in the optical, making it an ideal instrument to search for gravitational lenses down to the limiting magnitude of 21. Nevertheless, the previous Gaia Data Releases are known to be incomplete for small angular separations such as those ex… ▽ More Context. Strongly lensed quasars are fundamental sources for cosmology. The Gaia space mission covers the entire sky with the unprecedented resolution of $0.18$" in the optical, making it an ideal instrument to search for gravitational lenses down to the limiting magnitude of 21. Nevertheless, the previous Gaia Data Releases are known to be incomplete for small angular separations such as those expected for most lenses. Aims. We present the Data Processing and Analysis Consortium GravLens pipeline, which was built to analyse all Gaia detections around quasars and to cluster them into sources, thus producing a catalogue of secondary sources around each quasar. We analysed the resulting catalogue to produce scores that indicate source configurations that are compatible with strongly lensed quasars. Methods. GravLens uses the DBSCAN unsupervised clustering algorithm to detect sources around quasars. The resulting catalogue of multiplets is then analysed with several methods to identify potential gravitational lenses. We developed and applied an outlier scoring method, a comparison between the average BP and RP spectra of the components, and we also used an extremely randomised tree algorithm. These methods produce scores to identify the most probable configurations and to establish a list of lens candidates. Results. We analysed the environment of 3 760 032 quasars. A total of 4 760 920 sources, including the quasars, were found within 6" of the quasar positions. This list is given in the Gaia archive. In 87\% of cases, the quasar remains a single source, and in 501 385 cases neighbouring sources were detected. We propose a list of 381 lensed candidates, of which we identified 49 as the most promising. Beyond these candidates, the associate tables in this Focused Product Release allow the entire community to explore the unique Gaia data for strong lensing studies further. △ Less

Submitted 10 October, 2023; originally announced October 2023.

Comments: 35 pages, 60 figures, accepted for publication by Astronomy and Astrophysics

Journal ref: A&A 685, A130 (2024)

Gaia Focused Product Release: Radial velocity time series of long-period variables

Authors: Gaia Collaboration, Gaia Collaboration, M. Trabucchi, N. Mowlavi, T. Lebzelter, I. Lecoeur-Taibi, M. Audard, L. Eyer, P. García-Lario, P. Gavras, B. Holl, G. Jevardat de Fombelle, K. Nienartowicz, L. Rimoldini, P. Sartoretti, R. Blomme, Y. Frémat, O. Marchal, Y. Damerdji, A. G. A. Brown, A. Guerrier, P. Panuzzo, D. Katz, G. M. Seabroke, K. Benson , et al. (382 additional authors not shown)

Abstract: The third Gaia Data Release (DR3) provided photometric time series of more than 2 million long-period variable (LPV) candidates. Anticipating the publication of full radial-velocity (RV) in DR4, this Focused Product Release (FPR) provides RV time series for a selection of LPVs with high-quality observations. We describe the production and content of the Gaia catalog of LPV RV time series, and the… ▽ More The third Gaia Data Release (DR3) provided photometric time series of more than 2 million long-period variable (LPV) candidates. Anticipating the publication of full radial-velocity (RV) in DR4, this Focused Product Release (FPR) provides RV time series for a selection of LPVs with high-quality observations. We describe the production and content of the Gaia catalog of LPV RV time series, and the methods used to compute variability parameters published in the Gaia FPR. Starting from the DR3 LPVs catalog, we applied filters to construct a sample of sources with high-quality RV measurements. We modeled their RV and photometric time series to derive their periods and amplitudes, and further refined the sample by requiring compatibility between the RV period and at least one of the $G$, $G_{\rm BP}$, or $G_{\rm RP}$ photometric periods. The catalog includes RV time series and variability parameters for 9\,614 sources in the magnitude range $6\lesssim G/{\rm mag}\lesssim 14$, including a flagged top-quality subsample of 6\,093 stars whose RV periods are fully compatible with the values derived from the $G$, $G_{\rm BP}$, and $G_{\rm RP}$ photometric time series. The RV time series contain a mean of 24 measurements per source taken unevenly over a duration of about three years. We identify the great most sources (88%) as genuine LPVs, with about half of them showing a pulsation period and the other half displaying a long secondary period. The remaining 12% consists of candidate ellipsoidal binaries. Quality checks against RVs available in the literature show excellent agreement. We provide illustrative examples and cautionary remarks. The publication of RV time series for almost 10\,000 LPVs constitutes, by far, the largest such database available to date in the literature. The availability of simultaneous photometric measurements gives a unique added value to the Gaia catalog (abridged) △ Less

Submitted 9 October, 2023; originally announced October 2023.

Comments: 36 pages, 38 figures

Characterization of LAPPD timing at CERN PS testbeam

Authors: Deb Sankar Bhattacharya, Andrea Bressan, Chandradoy Chatterjee, Silvia Dalla Torre, Mauro Gregori, Alexander Kiselev, Stefano Levorato, Anna Martin, Saverio Minutoli, Mikhail Osipenko, Richa Rai, Marco Ripani, Fulvio Tessarotto, Triloki Triloki

Abstract: Large Area Picosecond PhotoDetectors (LAPPDs) are photosensors based on microchannel plate technology with about 400 cm$^2$ sensitive area. The external readout plane of a capacitively coupled LAPPD can be segmented into pads providing a spatial resolution down to 1 mm scale. The LAPPD signals have about 0.5 ns risetime followed by a slightly longer falltime and their amplitude reaches a few dozen… ▽ More Large Area Picosecond PhotoDetectors (LAPPDs) are photosensors based on microchannel plate technology with about 400 cm$^2$ sensitive area. The external readout plane of a capacitively coupled LAPPD can be segmented into pads providing a spatial resolution down to 1 mm scale. The LAPPD signals have about 0.5 ns risetime followed by a slightly longer falltime and their amplitude reaches a few dozens of mV per single photoelectron. In this article, we report on the measurement of the time resolution of an LAPPD prototype in a test beam exercise at CERN PS. Most of the previous measurements of LAPPD time resolution had been performed with laser sources. In this article we report time resolution measurements obtained through the detection of Cherenkov radiation emitted by high energy hadrons. Our approach has been demonstrated capable of measuring time resolutions as fine as 25-30 ps. The available prototype had performance limitations, which prevented us from applying the optimal high voltage setting. The measured time resolution for single photoelectrons is about 80 ps r.m.s. △ Less

Submitted 26 September, 2023; originally announced September 2023.

Comments: 35 pages, 23 figures

Dissecting the Gaia HR diagram II. The vertical structure of the star formation history across the Solar Cylinder

Authors: Alessandro Mazzi, Léo Girardi, Michele Trabucchi, Julianne J. Dalcanton, Rodrigo Luger, Paola Marigo, Andrea Miglio, Guglielmo Costa, Yang Chen, Giada Pastorelli, Morgan Fouesneau, Simone Zaggia, Alessandro Bressan, Piero Dal Tio

Abstract: Starting from the Gaia DR3 HR diagram, we derive the star formation history (SFH) as a function of distance from the Galactic Plane within a cylinder centred on the Sun with a 200~pc radius and spanning 1.3~kpc above and below the Galaxy's midplane. We quantify both the concentration of the more recent star formation in the Galactic Plane, and the age-related increase in the scale height of the Ga… ▽ More Starting from the Gaia DR3 HR diagram, we derive the star formation history (SFH) as a function of distance from the Galactic Plane within a cylinder centred on the Sun with a 200~pc radius and spanning 1.3~kpc above and below the Galaxy's midplane. We quantify both the concentration of the more recent star formation in the Galactic Plane, and the age-related increase in the scale height of the Galactic Disc stellar component, which is well-described by power-laws with indices ranging from $1/2$ to $2/3$. The vertically-integrated star formation rate falls from $(1.147 \pm 0.039)\times10^{-8}\, \text{M}_{\odot} \text{yr}^{-1} \text{pc}^{-2}$ at earlier times down to $(6.2 \pm 3.0) \times10^{-9}\, \text{M}_{\odot} \text{yr}^{-1} \text{pc}^{-2}$ at present times, but we find a significant peak of star formation in the 2 to 3 Gyr age bin. The total mass of stars formed per unit area over time is $118.7 \pm 6.2\, \text{M}_{\odot} \text{pc}^{-2}$, which is nearly twice the present stellar mass derived from kinematics within 1~kpc from the Galactic Plane, implying a high degree of matter recycling in successive generations of stars. The method is then modified by adopting an age-dependent correlation between the SFH across the different slices, which results in less noisy and more symmetrical results without significantly changing the previously mentioned quantities. This appears to be a promising way to improve SFH recovery in external galaxies. △ Less

Submitted 23 September, 2023; originally announced September 2023.

Comments: Accepted for publication in MNRAS

Long term experience with perfluorobutane in COMPASS RICH

Authors: F. Bradamante, A. Bressan, A. Cicuttin, M. L. Crespo, C. Chatterjee, P. Ciliberti, S. Dalla Torre, W. Florian, L. Garcia Ordonez, M. Gregori, A. Kerbizi, S. Levorato, A. Martin, G. Menon, R. S. Molina, A. Moretti, F. Tessarotto, Triloki, B. Valinoti

Abstract: COMPASS RICH-1 has used high-purity perfluorobutane as radiator gas since 2001. The operation and control of the radiator gas has evolved over years with continuous improvements. We report on the experience gained in the 20 year-long operation of perfluorobutane as COMPASS RICH radiator. Very accurate values for the radiator gas refractive index are needed for high-performance particle identificat… ▽ More COMPASS RICH-1 has used high-purity perfluorobutane as radiator gas since 2001. The operation and control of the radiator gas has evolved over years with continuous improvements. We report on the experience gained in the 20 year-long operation of perfluorobutane as COMPASS RICH radiator. Very accurate values for the radiator gas refractive index are needed for high-performance particle identification. The procedure has evolved over years and the one presently in use, which provides refractive index estimate at the 1 ppm level, is discussed. Perfluorobutane procurement is becoming challenging, and the minimization of material waste is now a priority for the protection of the environment. Commercially available perfluorobutane needs dedicated filtering before usage and typical material losses in the filtering procedure were around 30%. Recent efforts allowed us to reduce them to about 5%. A potential alternative to fluorocarbon radiators in gaseous RICHes is also presented. △ Less

Submitted 3 August, 2023; originally announced August 2023.

Comments: 4 pages, 3 figures 2 tables

Generic Singularities for 2D Pressureless Flow

Authors: Alberto Bressan, Geng Chen, Shoujun Huang

Abstract: In this paper, we consider the Cauchy problem for pressureless gases in two space dimensions with generic smooth initial data (density and velocity). These equations give rise to singular curves, where the mass has positive density w.r.t.~1-dimensional Hausdorff measure. We observe that the system of equations describing these singular curves is not hyperbolic. For analytic data, local solutions a… ▽ More In this paper, we consider the Cauchy problem for pressureless gases in two space dimensions with generic smooth initial data (density and velocity). These equations give rise to singular curves, where the mass has positive density w.r.t.~1-dimensional Hausdorff measure. We observe that the system of equations describing these singular curves is not hyperbolic. For analytic data, local solutions are constructed using a version of the Cauchy-Kovalevskaya theorem. We then study the interaction of two singular curves, in generic position. Finally, for a generic initial velocity field, we investigate the asymptotic structure of the smooth solution up to the first time when a singularity is formed. △ Less

Submitted 21 July, 2023; originally announced July 2023.

Comments: The manuscript contains 4 figures

The first comprehensive Milky Way stellar mock catalogue for the Chinese Space Station Telescope Survey Camera

Authors: Yang Chen, Xiaoting Fu, Chao Liu, Piero Dal Tio, Léo Girardi, Giada Pastorelli, Alessandro Mazzi, Michele Trabucchi, Hao Tian, Dongwei Fan, Paola Marigo, Alessandro Bressan

Abstract: The Chinese Space Station Telescope (CSST) is a cutting-edge two-meter astronomical space telescope currently under construction. Its primary Survey Camera (SC) is designed to conduct large-area imaging sky surveys using a sophisticated seven-band photometric system. The resulting data will provide unprecedented data for studying the structure and stellar populations of the Milky Way. To support t… ▽ More The Chinese Space Station Telescope (CSST) is a cutting-edge two-meter astronomical space telescope currently under construction. Its primary Survey Camera (SC) is designed to conduct large-area imaging sky surveys using a sophisticated seven-band photometric system. The resulting data will provide unprecedented data for studying the structure and stellar populations of the Milky Way. To support the CSST development and scientific projects related to its survey data, we generate the first comprehensive Milky Way stellar mock catalogue for the CSST SC photometric system using the TRILEGAL stellar population synthesis tool. The catalogue includes approximately 12.6 billion stars, covering a wide range of stellar parameters, photometry, astrometry, and kinematics, with magnitude reaching down to $g\,=\,27.5$ mag in the AB magnitude system. The catalogue represents our benchmark understanding of the stellar populations in the Milky Way, enabling a direct comparison with the future CSST survey data. Particularly, it sheds light on faint stars that are hidden from current sky surveys. Our crowding limit analysis based on this catalogue provides compelling evidence for the extension of the CSST Optical Survey (OS) to cover low Galactic latitude regions. The strategic extension of the CSST-OS coverage, combined with this comprehensive mock catalogue, will enable transformative science with the CSST. △ Less

Submitted 19 July, 2023; originally announced July 2023.

Comments: Accepted by Sci. China Phys. Mech. Astron, catalogue at https://nadc.china-vo.org/data/data/csst-trilegal/f, 16 pages, 9 figures

Optimal Solutions for a Class of Set-Valued Evolution Problems

Authors: Stefano Bianchini, Alberto Bressan, Maria Teresa Chiri

Abstract: The paper is concerned with a class of optimization problems for moving sets $t\mapstoΩ(t)\subset\mathbb{R}^2$, motivated by the control of invasive biological populations. Assuming that the initial contaminated set $Ω_0$ is convex, we prove that a strategy is optimal if an only if at each given time $t\in [0,T]$ the control is active along the portion of the boundary $\partial Ω(t)$ where the cur… ▽ More The paper is concerned with a class of optimization problems for moving sets $t\mapstoΩ(t)\subset\mathbb{R}^2$, motivated by the control of invasive biological populations. Assuming that the initial contaminated set $Ω_0$ is convex, we prove that a strategy is optimal if an only if at each given time $t\in [0,T]$ the control is active along the portion of the boundary $\partial Ω(t)$ where the curvature is maximal. In particular, this implies that $Ω(t)$ is convex for all $t\geq 0$. The proof relies on the analysis of a one-step constrained optimization problem, obtained by a time discretization. △ Less

Submitted 7 July, 2023; originally announced July 2023.

Comments: 41 pages, 18 figures

MSC Class: 49J53; 49Q20; 28A75; 49K05

A remark on the uniqueness of solutions to hyperbolic conservation laws

Authors: Alberto Bressan, Camillo De Lellis

Abstract: Given a strictly hyperbolic $n\times n$ system of conservation laws, it is well known that there exists a unique Lipschitz semigroup of weak solutions, defined on a domain of functions with small total variation, which are limits of vanishing viscosity approximations. Aim of this note is to prove that every weak solution taking values in the domain of the semigroup, and whose shocks satisfy the Li… ▽ More Given a strictly hyperbolic $n\times n$ system of conservation laws, it is well known that there exists a unique Lipschitz semigroup of weak solutions, defined on a domain of functions with small total variation, which are limits of vanishing viscosity approximations. Aim of this note is to prove that every weak solution taking values in the domain of the semigroup, and whose shocks satisfy the Liu admissibility conditions, actually coincides with a semigroup trajectory. △ Less

Submitted 26 May, 2023; originally announced May 2023.

Comments: 11 pages, 1 figure. arXiv admin note: text overlap with arXiv:2305.10737

MSC Class: 35L65

Unique Solutions to Hyperbolic Conservation Laws with a Strictly Convex Entropy

Authors: Alberto Bressan, Graziano Guerra

Abstract: Consider a strictly hyperbolic $n\times n$ system of conservation laws, where each characteristic field is either genuinely nonlinear or linearly degenerate. In this standard setting, it is well known that there exists a Lipschitz semigroup of weak solutions, defined on a domain of functions with small total variation. If the system admits a strictly convex entropy, we give a short proof that ever… ▽ More Consider a strictly hyperbolic $n\times n$ system of conservation laws, where each characteristic field is either genuinely nonlinear or linearly degenerate. In this standard setting, it is well known that there exists a Lipschitz semigroup of weak solutions, defined on a domain of functions with small total variation. If the system admits a strictly convex entropy, we give a short proof that every entropy weak solution taking values within the domain of the semigroup coincides with a semigroup trajectory. The result shows that the assumptions of ``Tame Variation" or ``Tame Oscillation", previously used to achieve uniqueness, can be removed in the presence of a strictly convex entropy. △ Less

Submitted 18 May, 2023; originally announced May 2023.

Comments: 15 pages, 3 figures

MSC Class: 35L65

Massive binary black holes from Population II and III stars

Authors: Guglielmo Costa, Michela Mapelli, Giuliano Iorio, Filippo Santoliquido, Gastón J. Escobar, Ralf S. Klessen, Alessandro Bressan

Abstract: Population III stars, born from the primordial gas in the Universe, lose a negligible fraction of their mass via stellar winds and possibly follow a top-heavy mass function. Hence, they have often been regarded as the ideal progenitors of massive black holes (BHs), even above the pair instability mass gap. Here, we evolve a large set of Population III binary stars (metallicity $Z=10^{-11}$) with o… ▽ More Population III stars, born from the primordial gas in the Universe, lose a negligible fraction of their mass via stellar winds and possibly follow a top-heavy mass function. Hence, they have often been regarded as the ideal progenitors of massive black holes (BHs), even above the pair instability mass gap. Here, we evolve a large set of Population III binary stars (metallicity $Z=10^{-11}$) with our population-synthesis code SEVN, and compare them with Population II binary stars ($Z=10^{-4}$). In our models, the lower edge of the pair-instability mass gap corresponds to a BH mass of $\approx{86}$ ($\approx{91}$) M$_\odot$ for single Population III (II) stars. Overall, we find only mild differences between the properties of binary BHs (BBHs) born from Population III and II stars, especially if we adopt the same initial mass function and initial orbital properties. Most BBH mergers born from Population III and II stars have primary BH mass below the pair-instability gap, and the maximum secondary BH mass is $ < 50$ M$_\odot$. Only up to $\approx{3.3}$% ($\approx{0.09}$%) BBH mergers from Population III (II) progenitors have primary mass above the gap. Unlike metal-rich binary stars, the main formation channel of BBH mergers from Population III and II stars involves only stable mass transfer episodes in our fiducial model. △ Less

Submitted 24 August, 2023; v1 submitted 27 March, 2023; originally announced March 2023.

Comments: 15 pages, 17 figures, Comments are welcome. Published in MNRAS

Controlled Traveling Profiles for Models of Invasive Biological Species

Authors: Alberto Bressan, Minyan Zhang

Abstract: We consider a family of controlled reaction-diffusion equations, describing the spatial spreading of an invasive biological species. For a given propagation speed $c\in{I\!\!R}$, we seek a control with minimum cost, which achieves a traveling profile with speed $c$. For various nonlinear models, the existence of a (possibly measure valued) optimal control is proved, together with necessary conditi… ▽ More We consider a family of controlled reaction-diffusion equations, describing the spatial spreading of an invasive biological species. For a given propagation speed $c\in{I\!\!R}$, we seek a control with minimum cost, which achieves a traveling profile with speed $c$. For various nonlinear models, the existence of a (possibly measure valued) optimal control is proved, together with necessary conditions for optimality. In the last section we study a case where the wave speed cannot be modified by any control with finite cost. The present analysis is motivated by the recent results in arXiv:2201.01723 and arXiv:2108.09321, showing how a control problem for a reaction-diffusion equation can be approximated by a simpler problem of optimal control of a moving set. △ Less

Submitted 19 February, 2023; originally announced February 2023.

MSC Class: 35C07; 49J20; 49J15; 49K15

The JWST Resolved Stellar Populations Early Release Science Program II. Survey Overview

Authors: Daniel R. Weisz, Kristen B. W. McQuinn, Alessandro Savino, Nitya Kallivayalil, Jay Anderson, Martha L. Boyer, Matteo Correnti, Marla C. Geha, Andrew E. Dolphin, Karin M. Sandstrom, Andrew A. Cole, Benjamin F. Williams, Evan D. Skillman, Roger E. Cohen, Max J. B. Newman, Rachael Beaton, Alessandro Bressan, Alberto Bolatto, Michael Boylan-Kolchin, Alyson M. Brooks, James S. Bullock, Charlie Conroy, M. C. Cooper, Julianne J. Dalcanton, Aaron L. Dotter , et al. (17 additional authors not shown)

Abstract: We present the JWST Resolved Stellar Populations Early Release Science (ERS) science program. We obtained 27.5 hours of NIRCam and NIRISS imaging of three targets in the Local Group (Milky Way globular cluster M92, ultra-faint dwarf galaxy Draco II, star-forming dwarf galaxy WLM), which span factors of $\sim10^5$ in luminosity, $\sim10^4$ in distance, and $\sim10^5$ in surface brightness. We descr… ▽ More We present the JWST Resolved Stellar Populations Early Release Science (ERS) science program. We obtained 27.5 hours of NIRCam and NIRISS imaging of three targets in the Local Group (Milky Way globular cluster M92, ultra-faint dwarf galaxy Draco II, star-forming dwarf galaxy WLM), which span factors of $\sim10^5$ in luminosity, $\sim10^4$ in distance, and $\sim10^5$ in surface brightness. We describe the survey strategy, scientific and technical goals, implementation details, present select NIRCam color-magnitude diagrams (CMDs), and validate the NIRCam exposure time calculator (ETC). Our CMDs are among the deepest in existence for each class of target. They touch the theoretical hydrogen burning limit in M92 ($<0.08$ $M_{\odot}$; SNR $\sim5$ at $m_{F090W}\sim28.2$; $M_{F090W}\sim+13.6$), include the lowest-mass stars observed outside the Milky Way in Draco II (0.09 $M_{\odot}$; SNR $=10$ at $m_{F090W}\sim29$; $M_{F090W}\sim+12.1$), and reach $\sim1.5$ magnitudes below the oldest main sequence turnoff in WLM (SNR $=10$ at $m_{F090W}\sim29.5$; $M_{F090W}\sim+4.6$). The PARSEC stellar models provide a good qualitative match to the NIRCam CMDs, though are $\sim0.05$ mag too blue compared to M92 F090W$-$F150W data. The NIRCam ETC (v2.0) matches the SNRs based on photon noise from DOLPHOT stellar photometry in uncrowded fields, but the ETC may not be accurate in more crowded fields, similar to what is known for HST. We release beta versions of DOLPHOT NIRCam and NIRISS modules to the community. Results from this ERS program will establish JWST as the premier instrument for resolved stellar populations studies for decades to come. △ Less

Submitted 11 January, 2023; originally announced January 2023.

Comments: 23 pages, 7 Figures, 3 Tables. Submitted to AAS Journals. Comments welcome

A Study of Primordial Very Massive Star Evolution

Authors: Guglielmo Volpato, Paola Marigo, Guglielmo Costa, Alessandro Bressan, Michele Trabucchi, Léo Girardi

Abstract: We present new evolutionary models of primordial very massive stars, with initial masses ranging from $100\,\mathrm{{M}_{\odot}}$ to $1000\,\mathrm{{M}_{\odot}}$, that extend from the main sequence until the onset of dynamical instability caused by the creation of electron-positron pairs during core C, Ne, or O burning, depending on the star's mass and metallicity. Mass loss accounts for radiation… ▽ More We present new evolutionary models of primordial very massive stars, with initial masses ranging from $100\,\mathrm{{M}_{\odot}}$ to $1000\,\mathrm{{M}_{\odot}}$, that extend from the main sequence until the onset of dynamical instability caused by the creation of electron-positron pairs during core C, Ne, or O burning, depending on the star's mass and metallicity. Mass loss accounts for radiation-driven winds as well as pulsation-driven mass-loss on the main sequence and during the red supergiant phase. After examining the evolutionary properties, we focus on the final outcome of the models and associated compact remnants. Stars that avoid the pair-instability supernova channel, should produce black holes with masses ranging from $\approx 40\, \mathrm{{M}_{\odot}}$ to $\approx 1000\,\mathrm{{M}_{\odot}}$. In particular, stars with initial masses of about $100\,\mathrm{{M}_{\odot}}$ could leave black holes of $\simeq 85-90\, \mathrm{{M}_{\odot}}$, values consistent with the estimated primary black hole mass of the GW190521 merger event. Overall, these results may contribute to explain future data from next-generation gravitational-wave detectors, such as the Einstein Telescope and Cosmic Explorer, which will have access to as-yet unexplored BH mass range of $\approx 10^2-10^4\,\mathrm{{M}_{\odot}}$ in the early universe. △ Less

Submitted 13 February, 2023; v1 submitted 19 December, 2022; originally announced December 2022.

Comments: 14 pages, 9 figures; zenodo link for wind ejecta tables

Journal ref: The Astrophysical Journal, 2023, Volume 944, Number 1, 14 pages

Compact object mergers: exploring uncertainties from stellar and binary evolution with SEVN

Authors: Giuliano Iorio, Michela Mapelli, Guglielmo Costa, Mario Spera, Gastón J. Escobar, Cecilia Sgalletta, Alessandro A. Trani, Erika Korb, Filippo Santoliquido, Marco Dall'Amico, Nicola Gaspari, Alessandro Bressan

Abstract: Population-synthesis codes are an unique tool to explore the parameter space of massive binary star evolution and binary compact object (BCO) formation. Most population-synthesis codes are based on the same stellar evolution model, limiting our ability to explore the main uncertainties. Here, we present the new version of the code SEVN, which overcomes this issue by interpolating the main stellar… ▽ More Population-synthesis codes are an unique tool to explore the parameter space of massive binary star evolution and binary compact object (BCO) formation. Most population-synthesis codes are based on the same stellar evolution model, limiting our ability to explore the main uncertainties. Here, we present the new version of the code SEVN, which overcomes this issue by interpolating the main stellar properties from a set of pre-computed evolutionary tracks. We describe the new interpolation and adaptive time-step algorithms of SEVN, and the main upgrades on single and binary evolution. With SEVN, we evolved $1.2\times10^9$ binaries in the metallicity range $0.0001\leq Z \leq 0.03$, exploring a number of models for electron-capture, core-collapse and pair-instability supernovae, different assumptions for common envelope, stability of mass transfer, quasi-homogeneous evolution and stellar tides. We find that stellar evolution has a dramatic impact on the formation of single and binary compact objects. Just by slightly changing the overshooting parameter ($λ_{\rm ov}=0.4,0.5$) and the pair-instability model, the maximum mass of a black hole can vary from $\approx{60}$ to $\approx{100}\ \mathrm{M}_\odot$. Furthermore, the formation channels of BCOs and the merger efficiency we obtain with SEVN show significant differences with respect to the results of other population-synthesis codes, even when the same binary-evolution parameters are used. For example, the main traditional formation channel of BCOs is strongly suppressed in our models: at high metallicity ($Z\gtrsim{0.01}$) only $<20$% of the merging binary black holes and binary neutron stars form via this channel, while other authors found fractions $>70$%. The local BCO merger rate density of our fiducial models is consistent with the most recent estimates by the LIGO--Virgo--KAGRA collaboration. △ Less

Submitted 29 June, 2023; v1 submitted 21 November, 2022; originally announced November 2022.

Comments: Accepted for publication in MNRAS. The SEVN code is available at https://gitlab.com/sevncodes/sevn.git. All the data underlying this article are available in Zenodo at the link https://doi.org/10.5281/zenodo.7794546. All the Jupyter notebooks used to produce the plots in the paper are available in the gitlab repository https://gitlab.com/iogiul/iorio22_plot.git

Collins and Sivers transverse-spin asymmetries in inclusive muoproduction of $ρ^0$ mesons

Authors: G. D. Alexeev, M. G. Alexeev, C. Alice, A. Amoroso, V. Andrieux, V. Anosov, K. Augsten, W. Augustyniak, C. D. R. Azevedo, B. Badelek, J. Barth, R. Beck, Y. Bedfer, J. Bernhard, M. Bodlak, F. Bradamante, A. Bressan, V. E. Burtsev, W. -C. Chang, C. Chatterjee, M. Chiosso, A. G. Chumakov, S. -U. Chung, A. Cicuttin, P. M. M. Correia , et al. (167 additional authors not shown)

Abstract: The production of vector mesons in deep inelastic scattering is an interesting yet scarsely explored channel to study the transverse spin structure of the nucleon and the related phenomena. The COMPASS collaboration has performed the first measurement of the Collins and Sivers asymmetries for inclusively produced $ρ^0$ mesons. The analysis is based on the data set collected in deep inelastic scatt… ▽ More The production of vector mesons in deep inelastic scattering is an interesting yet scarsely explored channel to study the transverse spin structure of the nucleon and the related phenomena. The COMPASS collaboration has performed the first measurement of the Collins and Sivers asymmetries for inclusively produced $ρ^0$ mesons. The analysis is based on the data set collected in deep inelastic scattering in $2010$ using a $160\,\,\rm{GeV}/c$ $μ^+$ beam impinging on a transversely polarized $\rm{NH}_3$ target. The $ρ^{0}$ mesons are selected from oppositely charged hadron pairs, and the asymmetries are extracted as a function of the Bjorken-$x$ variable, the transverse momentum of the pair and the fraction of the energy $z$ carried by the pair. Indications for positive Collins and Sivers asymmetries are observed. △ Less

Submitted 29 July, 2023; v1 submitted 31 October, 2022; originally announced November 2022.

Report number: CERN-EP-2022--234

Spin Density Matrix Elements in Exclusive $ρ^0$ Meson Muoproduction

Authors: G. D. Alexeev, M. G. Alexeev, C. Alice, A. Amoroso, V. Andrieux, V. Anosov, K. Augsten, W. Augustyniak, C. D. R. Azevedo, B. Badelek, J. Barth, R. Beck, Y. Bedfer, J. Bernhard, M. Bodlak, F. Bradamante, A. Bressan, V. E. Burtsev, W. -C. Chang, C. Chatterjee, M. Chiosso, A. G. Chumakov, S. -U. Chung, A. Cicuttin, P. M. M. Correia , et al. (165 additional authors not shown)

Abstract: We report on a measurement of Spin Density Matrix Elements (SDMEs) in hard exclusive $ρ^0$ meson muoproduction at COMPASS using 160~GeV/$c$ polarised $ μ^{+}$ and $ μ^{-}$ beams impinging on a liquid hydrogen target. The measurement covers the kinematic range 5.0~GeV/$c^2$ $< W <$ 17.0~GeV/$c^2$, 1.0 (GeV/$c$)$^2$ $< Q^2 <$ 10.0 (GeV/$c$)$^2$ and 0.01 (GeV/$c$)$^2$ $< p_{\rm{T}}^2 <$ 0.5 (GeV/$c$)… ▽ More We report on a measurement of Spin Density Matrix Elements (SDMEs) in hard exclusive $ρ^0$ meson muoproduction at COMPASS using 160~GeV/$c$ polarised $ μ^{+}$ and $ μ^{-}$ beams impinging on a liquid hydrogen target. The measurement covers the kinematic range 5.0~GeV/$c^2$ $< W <$ 17.0~GeV/$c^2$, 1.0 (GeV/$c$)$^2$ $< Q^2 <$ 10.0 (GeV/$c$)$^2$ and 0.01 (GeV/$c$)$^2$ $< p_{\rm{T}}^2 <$ 0.5 (GeV/$c$)$^2$. Here, $W$ denotes the mass of the final hadronic system, $Q^2$ the virtuality of the exchanged photon, and $p_{\rm{T}}$ the transverse momentum of the $ρ^0$ meson with respect to the virtual-photon direction. The measured non-zero SDMEs for the transitions of transversely polarised virtual photons to longitudinally polarised vector mesons ($γ^*_T \to V^{ }_L$) indicate a violation of $s$-channel helicity conservation. Additionally, we observe a dominant contribution of natural-parity-exchange transitions and a very small contribution of unnatural-parity-exchange transitions, which is compatible with zero within experimental uncertainties. The results provide important input for modelling Generalised Parton Distributions (GPDs). In particular, they may allow one to evaluate in a model-dependent way the role of parton helicity-flip GPDs in exclusive $ρ^0$ production. △ Less

Submitted 29 July, 2023; v1 submitted 30 October, 2022; originally announced October 2022.

Report number: CERN-EP-2022-231

Generic Properties of First Order Mean Field Games

Authors: Alberto Bressan, Khai T. Nguyen

Abstract: We consider a class of deterministic mean field games, where the state associated with each player evolves according to an ODE which is linear w.r.t. the control. Existence, uniqueness, and stability of solutions are studied from the point of view of generic theory. Within a suitable topological space of dynamics and cost functionals, we prove that, for nearly all mean field games(in the Baire cat… ▽ More We consider a class of deterministic mean field games, where the state associated with each player evolves according to an ODE which is linear w.r.t. the control. Existence, uniqueness, and stability of solutions are studied from the point of view of generic theory. Within a suitable topological space of dynamics and cost functionals, we prove that, for nearly all mean field games(in the Baire category sense) the best reply map is single valued for a.e. player. As a consequence, the mean field game admits a strong (not randomized) solution. Examples are given of open sets of games admitting a single solution, and other open sets admitting multiple solutions. Further examples show the existence of an open set of MFG having a unique solution which is asymptotically stable w.r.t. the best reply map, and another open set of MFG having a unique solution which is unstable. We conclude with an example of a MFG with terminal constraints which does not have any solution, not even in the mild sense with randomized strategies. △ Less

Submitted 26 October, 2022; originally announced October 2022.

Comments: 34 pages, 2 figures

ATHENA Detector Proposal -- A Totally Hermetic Electron Nucleus Apparatus proposed for IP6 at the Electron-Ion Collider

Authors: ATHENA Collaboration, J. Adam, L. Adamczyk, N. Agrawal, C. Aidala, W. Akers, M. Alekseev, M. M. Allen, F. Ameli, A. Angerami, P. Antonioli, N. J. Apadula, A. Aprahamian, W. Armstrong, M. Arratia, J. R. Arrington, A. Asaturyan, E. C. Aschenauer, K. Augsten, S. Aune, K. Bailey, C. Baldanza, M. Bansal, F. Barbosa, L. Barion , et al. (415 additional authors not shown)

Abstract: ATHENA has been designed as a general purpose detector capable of delivering the full scientific scope of the Electron-Ion Collider. Careful technology choices provide fine tracking and momentum resolution, high performance electromagnetic and hadronic calorimetry, hadron identification over a wide kinematic range, and near-complete hermeticity. This article describes the detector design and its e… ▽ More ATHENA has been designed as a general purpose detector capable of delivering the full scientific scope of the Electron-Ion Collider. Careful technology choices provide fine tracking and momentum resolution, high performance electromagnetic and hadronic calorimetry, hadron identification over a wide kinematic range, and near-complete hermeticity. This article describes the detector design and its expected performance in the most relevant physics channels. It includes an evaluation of detector technology choices, the technical challenges to realizing the detector and the R&D required to meet those challenges. △ Less

Submitted 13 October, 2022; originally announced October 2022.

Journal ref: JINST 17 (2022) 10, P10019

Low-Temperature Gas Opacities with AESOPUS 2.0

Authors: Paola Marigo, Bernhard Aringer, Leo Girardi, Alessandro Bressan

Abstract: This work introduces new low-temperature gas opacities, in the range 3.2 <= log(T/K) <= 4.5, computed with the AESOPUS code under the assumption of thermodynamic equilibrium (Marigo &_Aringer_2009). In comparison to the previous version AESOPUS 1.0, we updated and expanded molecular absorption to include 80 species, mostly using the recommended line lists currently available from the ExoMol and HI… ▽ More This work introduces new low-temperature gas opacities, in the range 3.2 <= log(T/K) <= 4.5, computed with the AESOPUS code under the assumption of thermodynamic equilibrium (Marigo &_Aringer_2009). In comparison to the previous version AESOPUS 1.0, we updated and expanded molecular absorption to include 80 species, mostly using the recommended line lists currently available from the ExoMol and HITRAN databases. Furthermore, in light of a recent study, we revised the H- photodetachment cross section, added the free-free absorption of other negative ions of atoms and molecules, and updated the collision-induced absorption due to H2/H2, H2/H, H2/He, and H/He pairs. Using the new input physics, we computed tables of Rosseland mean opacities for several scaled-solar chemical compositions, including Magg et al. (2022)'s most recent one, as well as alpha-enhanced mixtures. The differences in opacity between the new AESOPUS 2.0 and the original AESOPUS 1.0 versions, as well as other sets of calculations, are discussed. The new opacities are released to the community via a dedicated web-page that includes both pre-computed tables for widely used chemical compositions, and a web-interface for calculating opacities on-the-fly for any abundance distribution. △ Less

Submitted 16 October, 2022; originally announced October 2022.

Comments: Accepted for publication in ApJ. The paper includes 22 pages, 11 figures, 2 tables

The Gaia-ESO Public Spectroscopic Survey: Motivation, implementation, GIRAFFE data processing, analysis, and final data products

Authors: G. Gilmore, S. Randich, C. C. Worley, A. Hourihane, A. Gonneau, G. G. Sacco, J. R. Lewis, L. Magrini, P. Francois, R. D. Jeffries, S. E. Koposov, A. Bragaglia, E. J. Alfaro, C. Allende Prieto, R. Blomme, A. J. Korn, A. C. Lanzafame, E. Pancino, A. Recio-Blanco, R. Smiljanic, S. Van Eck, T. Zwitter, T. Bensby, E. Flaccomio, M. J. Irwin , et al. (143 additional authors not shown)

Abstract: The Gaia-ESO Public Spectroscopic Survey is an ambitious project designed to obtain astrophysical parameters and elemental abundances for 100,000 stars, including large representative samples of the stellar populations in the Galaxy, and a well-defined sample of 60 (plus 20 archive) open clusters. We provide internally consistent results calibrated on benchmark stars and star clusters, extending a… ▽ More The Gaia-ESO Public Spectroscopic Survey is an ambitious project designed to obtain astrophysical parameters and elemental abundances for 100,000 stars, including large representative samples of the stellar populations in the Galaxy, and a well-defined sample of 60 (plus 20 archive) open clusters. We provide internally consistent results calibrated on benchmark stars and star clusters, extending across a very wide range of abundances and ages. This provides a legacy data set of intrinsic value, and equally a large wide-ranging dataset that is of value for homogenisation of other and future stellar surveys and Gaia's astrophysical parameters. This article provides an overview of the survey methodology, the scientific aims, and the implementation, including a description of the data processing for the GIRAFFE spectra. A companion paper (arXiv:2206.02901) introduces the survey results. Gaia-ESO aspires to quantify both random and systematic contributions to measurement uncertainties. Thus all available spectroscopic analysis techniques are utilised, each spectrum being analysed by up to several different analysis pipelines, with considerable effort being made to homogenise and calibrate the resulting parameters. We describe here the sequence of activities up to delivery of processed data products to the ESO Science Archive Facility for open use. The Gaia-ESO Survey obtained 202,000 spectra of 115,000 stars using 340 allocated VLT nights between December 2011 and January 2018 from GIRAFFE and UVES. The full consistently reduced final data set of spectra was released through the ESO Science Archive Facility in late 2020, with the full astrophysical parameters sets following in 2022. △ Less

Submitted 10 August, 2022; originally announced August 2022.

Comments: 38 pages. A&A in press

Journal ref: A&A 666, A120 (2022)

Simulating the Legacy Survey of Space and Time stellar content with TRILEGAL

Authors: Piero Dal Tio, Giada Pastorelli, Alessandro Mazzi, Michele Trabucchi, Guglielmo Costa, Alice Jacques, Adriano Pieres, Léo Girardi, Yang Chen, Knut A. G. Olsen, Mario Juric, Željko Ivezić, Peter Yoachim, William I. Clarkson, Paola Marigo, Thaise S. Rodrigues, Simone Zaggia, Mauro Barbieri, Yazan Momany, Alessandro Bressan, Robert Nikutta, Luiz Nicolaci da Costa

Abstract: We describe a large simulation of the stars to be observed by the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). The simulation is based on the TRILEGAL code, which resorts to large databases of stellar evolutionary tracks, synthetic spectra, and pulsation models, added to simple prescriptions for the stellar density and star formation histories of the main structures of the Gal… ▽ More We describe a large simulation of the stars to be observed by the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). The simulation is based on the TRILEGAL code, which resorts to large databases of stellar evolutionary tracks, synthetic spectra, and pulsation models, added to simple prescriptions for the stellar density and star formation histories of the main structures of the Galaxy, to generate mock stellar samples through a population synthesis approach. The main bodies of the Magellanic Clouds are also included. A complete simulation is provided for single stars, down to the $r=27.5$ mag depth of the co-added wide-fast-deep survey images. A second simulation is provided for a fraction of the binaries, including the interacting ones, as derived with the BinaPSE module of TRILEGAL. We illustrate the main properties and numbers derived from these simulations, including: comparisons with real star counts; the expected numbers of Cepheids, long-period variables and eclipsing binaries; the crowding limits as a function of seeing and filter; the star-to-galaxy ratios, etc. Complete catalogs are accessible through the NOIRLab Astro Data Lab, while the stellar density maps are incorporated in the LSST metrics analysis framework (MAF). △ Less

Submitted 1 August, 2022; originally announced August 2022.

Comments: Accepted for publication in the LSST focused ApJS issue

Gaia Data Release 3: Summary of the content and survey properties

Authors: Gaia Collaboration, A. Vallenari, A. G. A. Brown, T. Prusti, J. H. J. de Bruijne, F. Arenou, C. Babusiaux, M. Biermann, O. L. Creevey, C. Ducourant, D. W. Evans, L. Eyer, R. Guerra, A. Hutton, C. Jordi, S. A. Klioner, U. L. Lammers, L. Lindegren, X. Luri, F. Mignard, C. Panem, D. Pourbaix, S. Randich, P. Sartoretti, C. Soubiran , et al. (431 additional authors not shown)

Abstract: We present the third data release of the European Space Agency's Gaia mission, GDR3. The GDR3 catalogue is the outcome of the processing of raw data collected with the Gaia instruments during the first 34 months of the mission by the Gaia Data Processing and Analysis Consortium. The GDR3 catalogue contains the same source list, celestial positions, proper motions, parallaxes, and broad band photom… ▽ More We present the third data release of the European Space Agency's Gaia mission, GDR3. The GDR3 catalogue is the outcome of the processing of raw data collected with the Gaia instruments during the first 34 months of the mission by the Gaia Data Processing and Analysis Consortium. The GDR3 catalogue contains the same source list, celestial positions, proper motions, parallaxes, and broad band photometry in the G, G$_{BP}$, and G$_{RP}$ pass-bands already present in the Early Third Data Release. GDR3 introduces an impressive wealth of new data products. More than 33 million objects in the ranges $G_{rvs} < 14$ and $3100 <T_{eff} <14500 $, have new determinations of their mean radial velocities based on data collected by Gaia. We provide G$_{rvs}$ magnitudes for most sources with radial velocities, and a line broadening parameter is listed for a subset of these. Mean Gaia spectra are made available to the community. The GDR3 catalogue includes about 1 million mean spectra from the radial velocity spectrometer, and about 220 million low-resolution blue and red prism photometer BPRP mean spectra. The results of the analysis of epoch photometry are provided for some 10 million sources across 24 variability types. GDR3 includes astrophysical parameters and source class probabilities for about 470 million and 1500 million sources, respectively, including stars, galaxies, and quasars. Orbital elements and trend parameters are provided for some $800\,000$ astrometric, spectroscopic and eclipsing binaries. More than $150\,000$ Solar System objects, including new discoveries, with preliminary orbital solutions and individual epoch observations are part of this release. Reflectance spectra derived from the epoch BPRP spectral data are published for about 60\,000 asteroids. Finally, an additional data set is provided, namely the Gaia Andromeda Photometric Survey (abridged) △ Less

Submitted 30 July, 2022; originally announced August 2022.

Comments: 23 pages, 2 figures

PARSEC V2.0: Stellar tracks and isochrones of low and intermediate mass stars with rotation

Authors: C. T. Nguyen, G. Costa, L. Girardi, G. Volpato, A. Bressan, Y. Chen, P. Marigo, X. Fu, P. Goudfrooij

Abstract: We present a new comprehensive collection of stellar evolutionary tracks and isochrones for rotating low- and intermediate-mass stars assembled with the updated version of PARSEC V2.0. This version includes our recent calibration of the extra mixing from overshooting and rotation, as well as several improvements in nuclear reaction network, treatment of convective zones, mass loss and other physic… ▽ More We present a new comprehensive collection of stellar evolutionary tracks and isochrones for rotating low- and intermediate-mass stars assembled with the updated version of PARSEC V2.0. This version includes our recent calibration of the extra mixing from overshooting and rotation, as well as several improvements in nuclear reaction network, treatment of convective zones, mass loss and other physical input parameters. The initial mass of the stellar models covers the range from 0.09 $M_\odot$ to 14 $M_\odot$, for six sets of initial metallicity, from Z=0.004 to Z=0.017. Rotation is considered for stars above $\sim 1 M_\odot$ with a smooth transition between non rotating and extremely fast-rotating models, based on the initial mass. For stars more massive than $\sim 1.3~M_\odot$ the full rotation range, from low to the critical one, is considered. We adopt the solar-scaled chemical mixture by Caffau et al. with Z$_\odot$ = 0.01524. All the evolutionary phases from the pre-main-sequence to the first few thermal-pulses on the asymptotic-giant-branch or central C exhaustion, are considered. The corresponding theoretical isochrones are further derived with TRILEGAL code and are converted in several photometric systems, taking into account different inclination angles. Besides magnitudes, they also offer many other stellar observables in line with the data that are being provided by current large surveys. The new collection is fully integrated in a user friendly WEB interface for the benefit of easily performing stellar population studies. △ Less

Submitted 22 July, 2022; v1 submitted 18 July, 2022; originally announced July 2022.

Comments: Accepted for publication in A&A. 17 pages. 20 Figures

Journal ref: A&A 665, A126 (2022)

ALMA resolves the first strongly-lensed Optical/NIR-dark galaxy

Authors: M. Giulietti, A. Lapi, M. Massardi, M. Behiri, M. Torsello, Q. D'Amato, T. Ronconi, F. Perrotta, A. Bressan

Abstract: We present high-resolution ($\lesssim0.1$arcsec) ALMA observations of the strongly-lensed galaxy HATLASJ113526.2-01460 at redshift $z\sim3.1$ discovered in the Gama 12$^{\rm th}$ field of the Herschel-ATLAS survey. The gravitationally lensed system is remarkably peculiar in that neither the background source nor the foreground lens show a clearly detected optical/NIR emission. We perform accurate… ▽ More We present high-resolution ($\lesssim0.1$arcsec) ALMA observations of the strongly-lensed galaxy HATLASJ113526.2-01460 at redshift $z\sim3.1$ discovered in the Gama 12$^{\rm th}$ field of the Herschel-ATLAS survey. The gravitationally lensed system is remarkably peculiar in that neither the background source nor the foreground lens show a clearly detected optical/NIR emission. We perform accurate lens modeling and source morphology reconstruction in three different (sub-)mm continuum bands, and in the C[II] and CO(8-7) spectral lines. The modeling indicates a foreground lensing (likely elliptical) galaxy with mass $\gtrsim10^{11}\, M_\odot$ at $z\gtrsim1.5$, while the source (sub-)mm continuum and line emissions are amplified by factors $μ\sim6-13$. We estimate extremely compact sizes $\lesssim0.5$ kpc for the star-forming region and $\lesssim 1$ kpc for the gas component, with no clear evidence of rotation or of ongoing merging events. We perform broadband SED-fitting and retrieve the intrinsic de-magnified physical properties of the source, which is found to feature a very high star-formation rate $\gtrsim10^3\, M_\odot$ yr$^{-1}$, that given the compact sizes is on the verge of the Eddington limit for starbursts; the radio luminosity at 6 cm from available EVLA observations is consistent with the star-formation activity. The galaxy is found to be extremely rich in gas $\sim10^{11}\, M_\odot$ and dust $\gtrsim10^9\, M_\odot$. The stellar content $\lesssim10^{11}\, M_\odot$ places the source well above the main sequence of starforming galaxies, indicating that the starburst is rather young with estimated age $\sim10^8$ yr. Our results indicate that the overall properties of HATLASJ113526.2-01460 are consistently explained by in-situ galaxy formation and evolution scenarios. △ Less

Submitted 1 July, 2022; originally announced July 2022.

Comments: 17 Pages, 12 Figures, 7 Tables. Submitted to ApJ

Gaia Data Release 3: Reflectance spectra of Solar System small bodies

Authors: Gaia Collaboration, L. Galluccio, M. Delbo, F. De Angeli, T. Pauwels, P. Tanga, F. Mignard, A. Cellino, A. G. A. Brown, K. Muinonen, A. Penttila, S. Jordan, A. Vallenari, T. Prusti, J. H. J. de Bruijne, F. Arenou, C. Babusiaux, M. Biermann, O. L. Creevey, C. Ducourant, D. W. Evans, L. Eyer, R. Guerra, A. Hutton, C. Jordi , et al. (422 additional authors not shown)

Abstract: The Gaia mission of the European Space Agency (ESA) has been routinely observing Solar System objects (SSOs) since the beginning of its operations in August 2014. The Gaia data release three (DR3) includes, for the first time, the mean reflectance spectra of a selected sample of 60 518 SSOs, primarily asteroids, observed between August 5, 2014, and May 28, 2017. Each reflectance spectrum was deriv… ▽ More The Gaia mission of the European Space Agency (ESA) has been routinely observing Solar System objects (SSOs) since the beginning of its operations in August 2014. The Gaia data release three (DR3) includes, for the first time, the mean reflectance spectra of a selected sample of 60 518 SSOs, primarily asteroids, observed between August 5, 2014, and May 28, 2017. Each reflectance spectrum was derived from measurements obtained by means of the Blue and Red photometers (BP/RP), which were binned in 16 discrete wavelength bands. We describe the processing of the Gaia spectral data of SSOs, explaining both the criteria used to select the subset of asteroid spectra published in Gaia DR3, and the different steps of our internal validation procedures. In order to further assess the quality of Gaia SSO reflectance spectra, we carried out external validation against SSO reflectance spectra obtained from ground-based and space-borne telescopes and available in the literature. For each selected SSO, an epoch reflectance was computed by dividing the calibrated spectrum observed by the BP/RP at each transit on the focal plane by the mean spectrum of a solar analogue. The latter was obtained by averaging the Gaia spectral measurements of a selected sample of stars known to have very similar spectra to that of the Sun. Finally, a mean of the epoch reflectance spectra was calculated in 16 spectral bands for each SSO. The agreement between Gaia mean reflectance spectra and those available in the literature is good for bright SSOs, regardless of their taxonomic spectral class. We identify an increase in the spectral slope of S-type SSOs with increasing phase angle. Moreover, we show that the spectral slope increases and the depth of the 1 um absorption band decreases for increasing ages of S-type asteroid families. △ Less

Submitted 24 June, 2022; originally announced June 2022.

Comments: 30 pages, 26 figures

The Black Hole Mass Function Across Cosmic Times II. Heavy Seeds and (Super)Massive Black Holes

Authors: A. Sicilia, A. Lapi, L. Boco, F. Shankar, D. M. Alexander, V. Allevato, C. Villforth, M. Massardi, M. Spera, A. Bressan, L. Danese

Abstract: This is the second paper in a series aimed at modeling the black hole (BH) mass function, from the stellar to the (super)massive regime. In the present work we focus on (super)massive BHs and provide an ab-initio computation of their mass function across cosmic times. We consider two main mechanisms to grow the central BH, that are expected to cooperate in the high-redshift star-forming progenitor… ▽ More This is the second paper in a series aimed at modeling the black hole (BH) mass function, from the stellar to the (super)massive regime. In the present work we focus on (super)massive BHs and provide an ab-initio computation of their mass function across cosmic times. We consider two main mechanisms to grow the central BH, that are expected to cooperate in the high-redshift star-forming progenitors of local massive galaxies. The first is the gaseous dynamical friction process, that can cause the migration toward the nuclear regions of stellar-mass BHs originated during the intense bursts of star formation in the gas-rich host progenitor galaxy, and the buildup of a central heavy BH seed $M_\bullet\sim 10^{3-5}\, M_\odot$ within short timescales $\lesssim$ some $10^7$ yr. The second mechanism is the standard Eddington-type gas disk accretion onto the heavy BH seed, through which the central BH can become (super)massive $M_\bullet\sim 10^{6-10}\, M_\odot$ within the typical star-formation duration $\lesssim 1$ Gyr of the host. We validate our semi-empirical approach by reproducing the observed redshift-dependent bolometric AGN luminosity functions and Eddington ratio distributions, and the relationship between the star-formation and the bolometric luminosity of the accreting central BH. We then derive the relic (super)massive BH mass function at different redshifts via a generalized continuity equation approach, and compare it with present observational estimates. Finally, we reconstruct the overall BH mass function from the stellar to the (super)massive regime, over more than ten orders of magnitudes in BH mass. △ Less

Submitted 15 June, 2022; originally announced June 2022.

Comments: 29 pages, 11 figures. Accepted by ApJ

Gaia Data Release 3: Mapping the asymmetric disc of the Milky Way

Authors: Gaia Collaboration, R. Drimmel, M. Romero-Gomez, L. Chemin, P. Ramos, E. Poggio, V. Ripepi, R. Andrae, R. Blomme, T. Cantat-Gaudin, A. Castro-Ginard, G. Clementini, F. Figueras, M. Fouesneau, Y. Fremat, K. Jardine, S. Khanna, A. Lobel, D. J. Marshall, T. Muraveva, A. G. A. Brown, A. Vallenari, T. Prusti, J. H. J. de Bruijne, F. Arenou , et al. (431 additional authors not shown)

Abstract: With the most recent Gaia data release the number of sources with complete 6D phase space information (position and velocity) has increased to well over 33 million stars, while stellar astrophysical parameters are provided for more than 470 million sources, in addition to the identification of over 11 million variable stars. Using the astrophysical parameters and variability classifications provid… ▽ More With the most recent Gaia data release the number of sources with complete 6D phase space information (position and velocity) has increased to well over 33 million stars, while stellar astrophysical parameters are provided for more than 470 million sources, in addition to the identification of over 11 million variable stars. Using the astrophysical parameters and variability classifications provided in Gaia DR3, we select various stellar populations to explore and identify non-axisymmetric features in the disc of the Milky Way in both configuration and velocity space. Using more about 580 thousand sources identified as hot OB stars, together with 988 known open clusters younger than 100 million years, we map the spiral structure associated with star formation 4-5 kpc from the Sun. We select over 2800 Classical Cepheids younger than 200 million years, which show spiral features extending as far as 10 kpc from the Sun in the outer disc. We also identify more than 8.7 million sources on the red giant branch (RGB), of which 5.7 million have line-of-sight velocities, allowing the velocity field of the Milky Way to be mapped as far as 8 kpc from the Sun, including the inner disc. The spiral structure revealed by the young populations is consistent with recent results using Gaia EDR3 astrometry and source lists based on near infrared photometry, showing the Local (Orion) arm to be at least 8 kpc long, and an outer arm consistent with what is seen in HI surveys, which seems to be a continuation of the Perseus arm into the third quadrant. Meanwhile, the subset of RGB stars with velocities clearly reveals the large scale kinematic signature of the bar in the inner disc, as well as evidence of streaming motions in the outer disc that might be associated with spiral arms or bar resonances. (abridged) △ Less

Submitted 5 August, 2022; v1 submitted 13 June, 2022; originally announced June 2022.

Comments: 35 pages, 27 figures, accepted for publication in A&A special Gaia DR3 issue. V2: abstract completed. V3: complete author list and link to data: https://drive.google.com/drive/u/1/folders/1yOJPjYmM7QK5XVsqaiSOTuwDQNti2LlZ

Journal ref: A&A 674, A37 (2023)

Gaia Data Release 3: Pulsations in main sequence OBAF-type stars

Authors: Gaia Collaboration, J. De Ridder, V. Ripepi, C. Aerts, L. Palaversa, L. Eyer, B. Holl, M. Audard, L. Rimoldini, A. G. A. Brown, A. Vallenari, T. Prusti, J. H. J. de Bruijne, F. Arenou, C. Babusiaux, M. Biermann, O. L. Creevey, C. Ducourant, D. W. Evans, R. Guerra, A. Hutton, C. Jordi, S. A. Klioner, U. L. Lammers, L. Lindegren , et al. (423 additional authors not shown)

Abstract: The third Gaia data release provides photometric time series covering 34 months for about 10 million stars. For many of those stars, a characterisation in Fourier space and their variability classification are also provided. This paper focuses on intermediate- to high-mass (IHM) main sequence pulsators M >= 1.3 Msun) of spectral types O, B, A, or F, known as beta Cep, slowly pulsating B (SPB), del… ▽ More The third Gaia data release provides photometric time series covering 34 months for about 10 million stars. For many of those stars, a characterisation in Fourier space and their variability classification are also provided. This paper focuses on intermediate- to high-mass (IHM) main sequence pulsators M >= 1.3 Msun) of spectral types O, B, A, or F, known as beta Cep, slowly pulsating B (SPB), delta Sct, and gamma Dor stars. These stars are often multi-periodic and display low amplitudes, making them challenging targets to analyse with sparse time series. All datasets used in this analysis are part of the Gaia DR3 data release. The photometric time series were used to perform a Fourier analysis, while the global astrophysical parameters necessary for the empirical instability strips were taken from the Gaia DR3 gspphot tables, and the vsini data were taken from the Gaia DR3 esphs tables. We show that for nearby OBAF-type pulsators, the Gaia DR3 data are precise and accurate enough to pinpoint them in the Hertzsprung-Russell diagram. We find empirical instability strips covering broader regions than theoretically predicted. In particular, our study reveals the presence of fast rotating gravity-mode pulsators outside the strips, as well as the co-existence of rotationally modulated variables inside the strips as reported before in the literature. We derive an extensive period-luminosity relation for delta Sct stars and provide evidence that the relation features different regimes depending on the oscillation period. Finally, we demonstrate how stellar rotation attenuates the amplitude of the dominant oscillation mode of delta Sct stars. △ Less

Submitted 16 August, 2022; v1 submitted 13 June, 2022; originally announced June 2022.

Journal ref: A&A 674, A36 (2023)

Gaia Data Release 3: A Golden Sample of Astrophysical Parameters

Authors: Gaia Collaboration, O. L. Creevey, L. M. Sarro, A. Lobel, E. Pancino, R. Andrae, R. L. Smart, G. Clementini, U. Heiter, A. J. Korn, M. Fouesneau, Y. Frémat, F. De Angeli, A. Vallenari, D. L. Harrison, F. Thévenin, C. Reylé, R. Sordo, A. Garofalo, A. G. A. Brown, L. Eyer, T. Prusti, J. H. J. de Bruijne, F. Arenou, C. Babusiaux , et al. (423 additional authors not shown)

Abstract: Gaia Data Release 3 (DR3) provides a wealth of new data products for the astronomical community to exploit, including astrophysical parameters for a half billion stars. In this work we demonstrate the high quality of these data products and illustrate their use in different astrophysical contexts. We query the astrophysical parameter tables along with other tables in Gaia DR3 to derive the samples… ▽ More Gaia Data Release 3 (DR3) provides a wealth of new data products for the astronomical community to exploit, including astrophysical parameters for a half billion stars. In this work we demonstrate the high quality of these data products and illustrate their use in different astrophysical contexts. We query the astrophysical parameter tables along with other tables in Gaia DR3 to derive the samples of the stars of interest. We validate our results by using the Gaia catalogue itself and by comparison with external data. We have produced six homogeneous samples of stars with high quality astrophysical parameters across the HR diagram for the community to exploit. We first focus on three samples that span a large parameter space: young massive disk stars (~3M), FGKM spectral type stars (~3M), and UCDs (~20K). We provide these sources along with additional information (either a flag or complementary parameters) as tables that are made available in the Gaia archive. We furthermore identify 15740 bone fide carbon stars, 5863 solar-analogues, and provide the first homogeneous set of stellar parameters of the Spectro Photometric Standard Stars. We use a subset of the OBA sample to illustrate its usefulness to analyse the Milky Way rotation curve. We then use the properties of the FGKM stars to analyse known exoplanet systems. We also analyse the ages of some unseen UCD-companions to the FGKM stars. We additionally predict the colours of the Sun in various passbands (Gaia, 2MASS, WISE) using the solar-analogue sample. △ Less

Submitted 12 June, 2022; originally announced June 2022.

Comments: 35 pages, (incl 6 pages references, acknowledgements, affiliations), 37 figures, A&A accepted

Journal ref: A&A 674, A39 (2023)