Showing 1–50 of 265 results for author: Silva, L

Minimum-entropy constraints on galactic potentials

Authors: Leandro Beraldo e Silva, Monica Valluri, Eugene Vasiliev, Kohei Hattori, Walter de Siqueira Pedra, Kathryne J. Daniel

Abstract: A tracer sample in a gravitational potential, starting from a generic initial condition, phase-mixes towards a stationary state. This evolution is accompanied by an entropy increase, and the final state is characterized by a distribution function (DF) that depends only on integrals of motion (Jeans theorem). In terms of angle-action variables, the final state is uniform in angles (high entropy) an… ▽ More A tracer sample in a gravitational potential, starting from a generic initial condition, phase-mixes towards a stationary state. This evolution is accompanied by an entropy increase, and the final state is characterized by a distribution function (DF) that depends only on integrals of motion (Jeans theorem). In terms of angle-action variables, the final state is uniform in angles (high entropy) and maximally clustered in actions (low entropy). We present a method exploring this fact to constrain a gravitational potential using a sample that is stationary in it. We estimate the entropy in the action space of trial potentials and recover the true potential by minimizing this entropy. This method avoids assuming a known DF, and may be applicable to other sets of integrals. We provide expressions for the entropy of DFs depending on energy, $f(E)$, energy and angular momentum, $f(E,L)$, or three actions, $f(\vec{J})$, and investigate the bias and fluctuations in their estimates. We show that the method correctly recovers the potential parameters for spherical and axisymmetric potentials. We also present a methodology to characterize the posterior probability distribution of the parameters with an Approximate Bayesian Computation, and indicate a pathway for application to observational data. Using $N=10^4$ tracers with $20\%$-uncertainties in the 6D coordinates, we recover the flattening parameter $q$ of an axisymmetric potential with $σ_q/q\sim 10\%$. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Comments: 17 pages + appendix + references, 13 figs. Comments are very welcome!

GD-1 Stellar Stream and Cocoon in the DESI Early Data Release

Authors: Monica Valluri, Parker Fagrelius, Sergey. E. Koposov, Ting S. Li, Oleg Y. Gnedin, Eric F. Bell, Raymond G. Carlberg, Andrew P. Cooper, Jessia N. Aguilar, Carlos Allende Prieto, Vasily Belokurov, Leandro Beraldo e Silva, David Brooks, Amanda Byström, Todd Claybaugh, Kyle Dawson, Arjun Dey, Peter Doel, Jaime E. Forero-Romero, Enrique Gaztañaga, Satya Gontcho A Gontcho, Klaus Honscheid, T . Kisner, Anthony Kremin, A. Lambert , et al. (27 additional authors not shown)

Abstract: We present ~ 126 new spectroscopically identified members of the GD-1 tidal stream obtained with the 5000-fiber Dark Energy Spectroscopic Instrument (DESI). We confirm the existence of a ``cocoon'' which is broad (FWHM~2.932deg~460pc) and kinematically hot (velocity dispersion, sigma~5-8km/s) component that surrounds a narrower (FWHM~0.353deg~55pc) and colder (sigma~ 2.2-2.6km/s) thin stream compo… ▽ More We present ~ 126 new spectroscopically identified members of the GD-1 tidal stream obtained with the 5000-fiber Dark Energy Spectroscopic Instrument (DESI). We confirm the existence of a ``cocoon'' which is broad (FWHM~2.932deg~460pc) and kinematically hot (velocity dispersion, sigma~5-8km/s) component that surrounds a narrower (FWHM~0.353deg~55pc) and colder (sigma~ 2.2-2.6km/s) thin stream component (based on a median per star velocity precision of 2.7km/s). The cocoon extends over at least a ~ 20deg segment of the stream observed by DESI. The thin and cocoon components have similar mean values of [Fe/H]: -2.54+/- 0.04dex and -2.45+/-0.06dex suggestive of a common origin. The data are consistent with the following scenarios for the origin of the cocoon. The progenitor of the GD-1 stream was an accreted globular cluster (GC) and: (a) the cocoon was produced by pre-accretion tidal stripping of the GC while it was still inside its parent dwarf galaxy; (b) the cocoon is debris from the parent dwarf galaxy; (c) an initially thin GC tidal stream was heated by impacts from dark subhalos in the Milky Way; (d) an initially thin GC stream was heated by a massive Sagittarius dwarf galaxy; or a combination of some these. In the first two cases the velocity dispersion and mean metallicity are consistent with the parent dwarf galaxy having a halo mass of ~0^9\msun. Future DESI spectroscopy and detailed modeling may enable us to distinguish between these possible origins. △ Less

Submitted 8 July, 2024; originally announced July 2024.

Comments: Submitted to ApJ, 23 pages, 13 figures 4 tables

DESI Early Data Release Milky Way Survey Value-Added Catalogue

Authors: Sergey E. Koposov, C. Allende-Prieto, A. P. Cooper, T. S. Li, L. Beraldo e Silva, B. Kim, A. Carrillo, A. Dey, C. J. Manser, F. Nikakhtar, A. H. Riley, C. Rockosi, M. Valluri, J. Aguilar, S. Ahlen, S. Bailey, R. Blum, D. Brooks, T. Claybaugh, S. Cole, A. de la Macorra, B. Dey, J. E. Forero-Romero, E. Gaztañaga, J. Guy , et al. (18 additional authors not shown)

Abstract: We present the stellar value-added catalogue based on the Dark Energy Spectroscopic Instrument (DESI) Early Data Release. The catalogue contains radial velocity and stellar parameter measurements for $\simeq$ 400,000 unique stars observed during commissioning and survey validation by DESI. These observations were made under conditions similar to the Milky Way Survey (MWS) currently carried out by… ▽ More We present the stellar value-added catalogue based on the Dark Energy Spectroscopic Instrument (DESI) Early Data Release. The catalogue contains radial velocity and stellar parameter measurements for $\simeq$ 400,000 unique stars observed during commissioning and survey validation by DESI. These observations were made under conditions similar to the Milky Way Survey (MWS) currently carried out by DESI but also include multiple specially targeted fields, such as those containing well-studied dwarf galaxies and stellar streams. The majority of observed stars have $16<r<20$ with a median signal-to-noise ratio in the spectra of $\sim$ 20. In the paper, we describe the structure of the catalogue, give an overview of different target classes observed, as well as provide recipes for selecting clean stellar samples. We validate the catalogue using external high-resolution measurements and show that radial velocities, surface gravities, and iron abundances determined by DESI are accurate to 1 km/s, $0.3$ dex and $\sim$ 0.15 dex respectively. We also demonstrate possible uses of the catalogue for chemo-dynamical studies of the Milky Way stellar halo and Draco dwarf spheroidal. The value-added catalogue described in this paper is the very first DESI MWS catalogue. The next DESI data release, expected in less than a year, will add the data from the first year of DESI survey operations and will contain approximately 4 million stars, along with significant processing improvements. △ Less

Submitted 8 July, 2024; originally announced July 2024.

Comments: Submitted to MNRAS; Value added catalogue is available at https://data.desi.lbl.gov/public/edr/vac/edr/mws/fuji/

The S-PLUS Ultra-Short Survey: first data release

Authors: Hélio D. Perottoni, Vinicius M. Placco, Felipe Almeida-Fernandes, Fábio R. Herpich, Silvia Rossi, Timothy C. Beers, Rodolfo Smiljanic, João A. S. Amarante, Guilherme Limberg, Ariel Werle, Helio J. Rocha-Pinto, Leandro Beraldo e Silva, Simone Daflon, Alvaro Alvarez-Candal, Gustavo B Oliveira Schwarz, William Schoenell, Tiago Ribeiro, Antonio Kanaan

Abstract: This paper presents the first public data release of the S-PLUS Ultra-Short Survey (USS), a photometric survey with short exposure times, covering approximately 9300 deg$^{2}$ of the Southern sky. The USS utilizes the Javalambre 12-band magnitude system, including narrow and medium-band and broad-band filters targeting prominent stellar spectral features. The primary objective of the USS is to ide… ▽ More This paper presents the first public data release of the S-PLUS Ultra-Short Survey (USS), a photometric survey with short exposure times, covering approximately 9300 deg$^{2}$ of the Southern sky. The USS utilizes the Javalambre 12-band magnitude system, including narrow and medium-band and broad-band filters targeting prominent stellar spectral features. The primary objective of the USS is to identify bright, extremely metal-poor (EMP; [Fe/H] $\leq -3$) and ultra metal-poor (UMP; [Fe/H] $\leq -4$) stars for further analysis using medium- and high-resolution spectroscopy.}{This paper provides an overview of the survey observations, calibration method, data quality, and data products. Additionally, it presents the selection of EMP and UMP candidates.}{The data from the USS were reduced and calibrated using the same methods as presented in the S-PLUS DR2. An additional step was introduced, accounting for the offset between the observed magnitudes off the USS and the predicted magnitudes from the very low-resolution Gaia XP spectra.}{This first release contains data for 163 observed fields totaling $\sim$324 deg$^{2}$ along the Celestial Equator. The magnitudes obtained from the USS are well-calibrated, showing a difference of $\sim 15$ mmag compared to the predicted magnitudes by the GaiaXPy toolkit. By combining colors and magnitudes, 140 candidates for EMP or UMP have been identified for follow-up studies.}{The S-PLUS USS DR1 is an important milestone in the search for bright metal-poor stars, with magnitudes in the range 10 $ < r \leq 14$. The USS is an ongoing survey; in the near future, it will provide many more bright metal-poor candidate stars for spectroscopic follow-up. △ Less

Submitted 6 July, 2024; originally announced July 2024.

Comments: Accepted for publication in A&A. 17 pages 6 figures. Long table at the end of the paper

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.

Spontaneous Lorentz symmetry-breaking constraints in Kalb-Ramond gravity

Authors: Ednaldo L. B. Junior, José Tarciso S. S. Junior, Francisco S. N. Lobo, Manuel E. Rodrigues, Diego Rubiera-Garcia, Luís F. Dias da Silva, Henrique A. Vieira

Abstract: In this work, we study timelike and lightlike geodesics in Kalb-Ramond (KR) gravity around a black hole with the goal of constraining the Lorentz symmetry-breaking parameter $l$. The analysis involves studying the precession of the S2 star periastron orbiting Sgr A* and geodesic precession around the Earth. The ratio of precession frequencies for General Relativity (GR) and KR gravity is computed,… ▽ More In this work, we study timelike and lightlike geodesics in Kalb-Ramond (KR) gravity around a black hole with the goal of constraining the Lorentz symmetry-breaking parameter $l$. The analysis involves studying the precession of the S2 star periastron orbiting Sgr A* and geodesic precession around the Earth. The ratio of precession frequencies for General Relativity (GR) and KR gravity is computed, with Event Horizon Telescope (EHT) results providing a parameter range for the spontaneous symmetry-breaking of $-0.185022 \leq l \leq 0.060938$. Utilizing the geodesic precession frequency from the Gravity Probe B (GP-B), the $l$ parameter is further constrained to $-6.30714 \times 10^{-12} \leq l \leq 3.90708 \times 10^{-12}$, which is consistent with the Schwarzschild limits. Moreover, for timelike geodesics, the innermost circular orbit (ICO) and innermost stable circular orbit (ISCO) are determined and analyzed to illustrate the impact of the symmetry breaking term. Zoom-whirl obstructions are compared with the Schwarzschild solution. Lower and upper limits of the photon sphere for lightlike geodesics are established to demonstrate the influence of KR gravity on the photon sphere. Additionally, the shadow radius is determined for two observers, one situated at a finite distance from the KR black hole, and the other located at an infinite distance, to constrain the symmetry-breaking parameter $l$, with comparisons made to EHT results. The bounds for $l$ derived from constraints on the photon sphere radius for lightlike geodesics yield $-0.0700225 \leq l \leq 0.189785$ using EHT data. The findings of this paper align with experimental results in the $l \rightarrow 0$ limit. △ Less

Submitted 6 May, 2024; originally announced May 2024.

Comments: 11 pages, 4 figures

Gravitational lensing of a Schwarzschild-like black hole in Kalb-Ramond gravity

Authors: Ednaldo L. B. Junior, José Tarciso S. S. Junior, Francisco S. N. Lobo, Manuel E. Rodrigues, Diego Rubiera-Garcia, Luís F. Dias da Silva, Henrique A. Vieira

Abstract: In this paper, we investigate the gravitational lensing effect for the Schwarzschild-like black hole spacetime in the background of a Kalb-Ramond (KR) field proposed in [K. Yang et. al., Phys. Rev. D 108 (2023) 124004]. The solution is characterized by a single extra parameter $l$, which is associated to the Lorentz symmetry breaking induced by the KR field. First, we calculate the exact deflectio… ▽ More In this paper, we investigate the gravitational lensing effect for the Schwarzschild-like black hole spacetime in the background of a Kalb-Ramond (KR) field proposed in [K. Yang et. al., Phys. Rev. D 108 (2023) 124004]. The solution is characterized by a single extra parameter $l$, which is associated to the Lorentz symmetry breaking induced by the KR field. First, we calculate the exact deflection angle of massive and massless particles for finite distances using elliptic integrals. Then we study this effect in the weak and strong field regimes, discussing the correction of the KR parameter on the coefficients of the expansions in both limits. We also find that increasing $l$ decreases the deflection angle. Furthermore, we use the available data from the Sagittarius $A^{\star}$ object, which is believed to be a supermassive black hole at the center of our galaxy, to calculate relevant observables, such as, the image position, luminosity, and delay time. The values found could be potentially measured in the weak field regime, though for strong fields one would have to wait for the next generation of interferometers. △ Less

Submitted 15 May, 2024; v1 submitted 6 May, 2024; originally announced May 2024.

Comments: 16 pages, 8 figures. V2: references added. arXiv admin note: text overlap with arXiv:2309.02658

Ring momentum distributions as a general feature of Vlasov dynamics in the synchrotron dominated regime

Authors: Pablo. J. Bilbao, Robert J. Ewart, Francisco Assunçao, Thales Silva, Luis O. Silva

Abstract: We study how radiation reaction leads plasmas initially in kinetic equilibrium to develop features in momentum space, such as anisotropies and population inversion, resulting in a ring-shaped momentum distribution that can drive kinetic instabilities. We employ the Landau-Lifshiftz radiation reaction model for a plasma in a strong magnetic field, and we obtain the necessary condition for the devel… ▽ More We study how radiation reaction leads plasmas initially in kinetic equilibrium to develop features in momentum space, such as anisotropies and population inversion, resulting in a ring-shaped momentum distribution that can drive kinetic instabilities. We employ the Landau-Lifshiftz radiation reaction model for a plasma in a strong magnetic field, and we obtain the necessary condition for the development of population inversion, we show that isotropic Maxwellian and Maxwell-Jüttner plasmas, with thermal temperature $T>m_e c^2/\sqrt{3}$, will develop a ring-like momentum distribution. The timescales and features for forming ring-shaped momentum distributions, the effect of collisions and non-uniform magnetic fields are disscussed, and compared with typical astrophysical and laboratory plasmas parameters. Our results show the pervasiveness of ring-like momentum distribution functions in synchrotron dominated plasma conditions. △ Less

Submitted 17 April, 2024; originally announced April 2024.

Comments: 10 pages, 3 Figures, Under consideration for PoP

Journal ref: Phys. Plasmas 1 May 2024; 31 (5): 052112

The ďAlembert solution in hyperboloidal foliations

Authors: Juan A. Valiente Kroon, Lidia J. Gomes Da Silva

Abstract: We explicitly construct the analogue of the ďAlembert solution to the 1+1 wave equation in an hyperboloidal setting. This hyperboloidal ďAlembert solution is used, in turn, to gain intuition into the behaviour of solutions to the wave equation in a hyperboloidal foliation and to explain some apparently anomalous behaviour observed in numerically constructed solutions discussed in the literature. We explicitly construct the analogue of the ďAlembert solution to the 1+1 wave equation in an hyperboloidal setting. This hyperboloidal ďAlembert solution is used, in turn, to gain intuition into the behaviour of solutions to the wave equation in a hyperboloidal foliation and to explain some apparently anomalous behaviour observed in numerically constructed solutions discussed in the literature. △ Less

Submitted 11 March, 2024; originally announced March 2024.

Comments: 15 pages, 7 figures. Submitted to Springer's topical collection: "Hyperboloidal foliations in the era of gravitational-wave astronomy: from mathematical relativity to astrophysics". Comments welcomed

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)

Azimuthal metallicity variations, spiral structure, and the failure of radial actions based on assuming axisymmetry

Authors: Victor P. Debattista, Tigran Khachaturyants, Joao A. S. Amarante, Christopher Carr, Leandro Beraldo e Silva, Chervin F. P. Laporte

Abstract: We study azimuthal variations in the mean stellar metallicity, <[Fe/H]>, in a self-consistent, isolated simulation in which all stars form out of gas. We find <[Fe/H]> variations comparable to those observed in the Milky Way and which are coincident with the spiral density waves. The azimuthal variations are present in young and old stars and therefore are not a result of recently formed stars. Si… ▽ More We study azimuthal variations in the mean stellar metallicity, <[Fe/H]>, in a self-consistent, isolated simulation in which all stars form out of gas. We find <[Fe/H]> variations comparable to those observed in the Milky Way and which are coincident with the spiral density waves. The azimuthal variations are present in young and old stars and therefore are not a result of recently formed stars. Similar variations are present in the mean age and alpha-abundance. We measure the pattern speeds of the <[Fe/H]>-variations and find that they match those of the spirals, indicating that they are at the origin of the metallicity patterns. Because younger stellar populations are not just more [Fe/H]-rich and alpha-poor but also dynamically cooler, we expect them to more strongly support spirals, which is indeed the case in the simulation. However, if we measure the radial action, J_R, using the Stackel axisymmetric approximation, we find that the spiral ridges are traced by regions of high J_R, contrary to expectations. Assuming that the passage of stars through the spirals leads to unphysical variations in the measured J_R, we obtain an improved estimate of J_R by averaging over a 1 Gyr time interval. This time-averaged J_R is a much better tracer of the spiral structure, with minima at the spiral ridges. We conclude that the errors incurred by the axisymmetric approximation introduce correlated deviations large enough to render the instantaneous radial action inadequate for tracing spirals. △ Less

Submitted 13 February, 2024; originally announced February 2024.

Comments: 21 pages. Submitted to MNRAS. Comments welcomed

Laboratory study of magnetic reconnection in lunar-relevant mini-magnetospheres

Authors: Lucas Rovige, Filipe D. Cruz, Robert S. Dorst, Jessica J. Pilgram, Carmen G. Constantin, Stephen Vincena, Fábio Cruz, Luis O. Silva, Christoph Niemann, Derek B. Schaeffer

Abstract: Mini-magnetospheres are small ion-scale structures that are well-suited to studying kinetic-scale physics of collisionless space plasmas. Such ion-scale magnetospheres can be found on local regions of the Moon, associated with the lunar crustal magnetic field. In this paper, we report on the laboratory experimental study of magnetic reconnection in laser-driven, lunar-like ion-scale magnetospheres… ▽ More Mini-magnetospheres are small ion-scale structures that are well-suited to studying kinetic-scale physics of collisionless space plasmas. Such ion-scale magnetospheres can be found on local regions of the Moon, associated with the lunar crustal magnetic field. In this paper, we report on the laboratory experimental study of magnetic reconnection in laser-driven, lunar-like ion-scale magnetospheres on the Large Plasma Device (LAPD) at the University of California - Los Angeles. In the experiment, a high-repetition rate (1 Hz), nanosecond laser is used to drive a fast moving, collisionless plasma that expands into the field generated by a pulsed magnetic dipole embedded into a background plasma and magnetic field. The high-repetition rate enables the acquisition of time-resolved volumetric data of the magnetic and electric fields to characterize magnetic reconnection and calculate the reconnection rate. We notably observe the formation of Hall fields associated with reconnection. Particle-in-cell simulations reproducing the experimental results were performed to study the micro-physics of the interaction. By analyzing the generalized Ohm's law terms, we find that the electron-only reconnection is driven by kinetic effects, through the electron pressure anisotropy. These results are compared to recent satellite measurements that found evidence of magnetic reconnection near the lunar surface. △ Less

Submitted 7 February, 2024; originally announced February 2024.

DiscoTEX: Discontinuous collocation and implicit-turned-explicit (IMTEX) integration symplectic, symmetric numerical algorithms with higher order jumps for differential equations with numerical black hole perturbation theory applications

Authors: Lidia J. Gomes Da Silva

Abstract: Dirac delta distributionally sourced differential equations emerge in many dynamical physical systems from neuroscience to black hole perturbation theory. Most of these lack exact analytical solutions and are thus best tackled numerically. This work describes a generic numerical algorithm which constructs discontinuous spatial and temporal discretisations by operating on discontinuous Lagrange and… ▽ More Dirac delta distributionally sourced differential equations emerge in many dynamical physical systems from neuroscience to black hole perturbation theory. Most of these lack exact analytical solutions and are thus best tackled numerically. This work describes a generic numerical algorithm which constructs discontinuous spatial and temporal discretisations by operating on discontinuous Lagrange and Hermite interpolation formulae recovering higher order accuracy. It is shown by solving the distributionally sourced wave equation, which has analytical solutions, that numerical weak-form solutions can be recovered to high order accuracy by solving a first-order reduced system of ordinary differential equations. The method-of-lines framework is applied to the DiscoTEX algorithm i.e through discontinuous collocation with implicit-turned-explicit (IMTEX) integration methods which are symmetric and conserve symplectic structure. Furthermore, the main application of the algorithm is proved, for the first-time, by calculating the amplitude at any desired location within the numerical grid, including at the position (and at its right and left limit) where the wave- (or wave-like) equation is discontinuous via interpolation using DiscoTEX. This is shown, firstly by solving the wave- (or wave-like) equation and comparing the numerical weak-form solution to the exact solution. Finally, one shows how to reconstruct the scalar and gravitational metric perturbations from weak-form numerical solutions of a non-rotating black hole, which do not have known exact analytical solutions, and compare against state-of-the-art frequency domain results. One concludes by motivating how DiscoTEX, and related algorithms, open a promising new alternative Extreme-Mass-Ratio-Inspiral (EMRI)s waveform generation route via a self-consistent evolution for the gravitational self-force programme in the time-domain. △ Less

Submitted 16 January, 2024; originally announced January 2024.

Comments: 44 pages, 17 figures, 8 tables. Includes overview of previous numerical methods implemented in the time-domain for the modelling of Extreme-Mass-Ratio-Inspirals with suitability checks on Table 8. Comments are welcome

OSIRIS-GR: General relativistic activation of the polar cap of a compact neutron star

Authors: R. Torres, T. Grismayer, F. Cruz, R. A. Fonseca, L. O. Silva

Abstract: We present ab initio global general-relativistic Particle-in-cell (GR-PIC) simulations of compact millisecond neutron star magnetospheres in the axisymmetric aligned rotator configuration. We investigate the role of GR and plasma supply on the polar cap particle acceleration efficiency - the precursor of coherent radio emission - employing a new module for the PIC code OSIRIS, designed to model pl… ▽ More We present ab initio global general-relativistic Particle-in-cell (GR-PIC) simulations of compact millisecond neutron star magnetospheres in the axisymmetric aligned rotator configuration. We investigate the role of GR and plasma supply on the polar cap particle acceleration efficiency - the precursor of coherent radio emission - employing a new module for the PIC code OSIRIS, designed to model plasma dynamics around compact objects with fully self-consistent GR effects. We provide a detailed description of the main sub-algorithms of the novel PIC algorithm, including a charge-conserving current deposit scheme for curvilinear coordinates. We demonstrate efficient particle acceleration in the polar caps of compact neutron stars with denser magnetospheres, numerically validating the spacelike current extension provided by force-free models. We show that GR relaxes the minimum required poloidal magnetospheric current for the transition of the polar cap to the accelerator regime, thus justifying the observation of weak pulsars beyond the expected death line. We denote that spin-down luminosity intermittency and radio pulse nullings for older pulsars might arise from the interplay between the polar and outer gaps. Also, narrower radio beams are expected for weaker low-obliquity pulsars. △ Less

Submitted 1 June, 2024; v1 submitted 5 January, 2024; originally announced January 2024.

Comments: 23 pages, 8 figures, 2 tables, accepted to New Astronomy

Journal ref: New Astronomy, Volume 112, June 2024, Pages 102261

Observations on the massive particle surface method

Authors: Ednaldo L. B. Junior, José Tarciso S. S. Junior, Francisco S. N. Lobo, Manuel E. Rodrigues, Luís F. Dias da Silva, Henrique A. Vieira

Abstract: The geodesic method has played a crucial role in understanding the circular orbits generated by compact objects, culminating in the definition of the photon sphere, which was later generalized to a photon surface in arbitrary spacetimes. This new formulation extends the concept of the photon sphere in a broader sense, including dynamical spacetimes, as shown by the Vaidya solution. The photon surf… ▽ More The geodesic method has played a crucial role in understanding the circular orbits generated by compact objects, culminating in the definition of the photon sphere, which was later generalized to a photon surface in arbitrary spacetimes. This new formulation extends the concept of the photon sphere in a broader sense, including dynamical spacetimes, as shown by the Vaidya solution. The photon surface essentially defines the null geodesics, which are originally tangent to the temporal surface, and keeps them confined to this surface. However, this formalism does not cover all classes of particles, and to overcome this limitation, a more comprehensive approach, denoted as the "massive particle surface", has been proposed that also accounts for charged massive particles. Indeed, the photon surface concept is recovered when the charge and mass of the particles are zero. In this work, we use these three formalisms to check the consistency of the results for the values of the radius of the photon sphere ($r_{ps}$) and the radius of the "innermost stable circular orbit" (ISCO) ($r_{\rm ISCO}$) for some gravitational models. In our results, the first model is described by conformal gravity, with the peculiarity that $g_{00}\neq-g_{11}^{-1}$. The second model, i.e. the Culetu solution, is developed by coupling General Relativity with nonlinear electrodynamics, which requires the consideration of the effective metric ($g_{\rm eff}^{μν}$) for geodesic approaches. Furthermore, we have also analysed the expressions for $r_{ps}$ and $r_{\rm ISCO}$ in a general static and spherically symmetric metric. Under these circumstances, we have found a discrepancy of $r_{ps}$ and $r_{\rm ISCO}$ obtained by the massive particle surface formalism as compared to the geodesic and photon surface formalisms. △ Less

Submitted 3 January, 2024; originally announced January 2024.

Comments: 15 pages, 1 figure; comments are welcome

Laboratory realization of relativistic pair-plasma beams

Authors: C. D. Arrowsmith, P. Simon, P. Bilbao, A. F. A. Bott, S. Burger, H. Chen, F. D. Cruz, T. Davenne, I. Efthymiopoulos, D. H. Froula, A. M. Goillot, J. T. Gudmundsson, D. Haberberger, J. Halliday, T. Hodge, B. T. Huffman, S. Iaquinta, F. Miniati, B. Reville, S. Sarkar, A. A. Schekochihin, L. O. Silva, R. Simpson, V. Stergiou, R. M. G. M. Trines , et al. (4 additional authors not shown)

Abstract: Relativistic electron-positron plasmas are ubiquitous in extreme astrophysical environments such as black holes and neutron star magnetospheres, where accretion-powered jets and pulsar winds are expected to be enriched with such pair plasmas. Their behaviour is quite different from typical electron-ion plasmas due to the matter-antimatter symmetry of the charged components and their role in the dy… ▽ More Relativistic electron-positron plasmas are ubiquitous in extreme astrophysical environments such as black holes and neutron star magnetospheres, where accretion-powered jets and pulsar winds are expected to be enriched with such pair plasmas. Their behaviour is quite different from typical electron-ion plasmas due to the matter-antimatter symmetry of the charged components and their role in the dynamics of such compact objects is believed to be fundamental. So far, our experimental inability to produce large yields of positrons in quasi-neutral beams has restricted the understanding of electron-positron pair plasmas to simple numerical and analytical studies which are rather limited. We present first experimental results confirming the generation of high-density, quasi-neutral, relativistic electron-positron pair beams using the 440 GeV/c beam at CERN's Super Proton Synchrotron (SPS) accelerator. The produced pair beams have a volume that fills multiple Debye spheres and are thus able to sustain collective plasma oscillations. Our work opens up the possibility of directly probing the microphysics of pair plasmas beyond quasi-linear evolution into regimes that are challenging to simulate or measure via astronomical observations. △ Less

Submitted 8 December, 2023; originally announced December 2023.

Comments: 14 pages, 8 figures

Characterisation of high velocity stars in the S-PLUS internal fourth data release

Authors: F. Quispe-Huaynasi, F. Roig, V. M. Placco, L. Beraldo e Silva, S. Daflon, C. B. Pereira, A. Kanaan, C. Mendes de Oliveira, T. Ribeiro, W. Schoenell

Abstract: In general, the atypical high velocity of some stars in the Galaxy can only be explained by invoking acceleration mechanisms related to extreme astrophysical events in the Milky Way. Using astrometric data from Gaia and the photometric information in 12 filters of the S-PLUS, we performed a kinematic, dynamical, and chemical analysis of 64 stars with galactocentric velocities higher than 400… ▽ More In general, the atypical high velocity of some stars in the Galaxy can only be explained by invoking acceleration mechanisms related to extreme astrophysical events in the Milky Way. Using astrometric data from Gaia and the photometric information in 12 filters of the S-PLUS, we performed a kinematic, dynamical, and chemical analysis of 64 stars with galactocentric velocities higher than 400 $\mathrm{km\,s}^{-1}$. All the stars are gravitationally bound to the Galaxy and exhibit halo kinematics. Some of the stars could be remnants of structures such as the Sequoia and the Gaia-Sausage/Enceladus. Supported by orbital and chemical analysis, we identified Gaia DR3 5401875170994688896 as a star likely to be originated at the centre of the Galaxy. Application of a machine learning technique to the S-PLUS photometric data allows us to obtain very good estimates of magnesium abundances for this sample of high velocity stars. △ Less

Submitted 20 November, 2023; originally announced November 2023.

Particle-in-cell simulations of pulsar magnetospheres: transition between electrosphere and force-free regimes

Authors: Fábio Cruz, Thomas Grismayer, Alexander Y. Chen, Anatoly Spitkovsky, Ricardo A. Fonseca, Luis O. Silva

Abstract: Global particle-in-cell (PIC) simulations of pulsar magnetospheres are performed with a volume, surface and pair production-based plasma injection schemes to systematically investigate the transition between electrosphere and force-free pulsar magnetospheric regimes. A new extension of the PIC code OSIRIS to model pulsar magnetospheres using a two-dimensional axisymmetric spherical grid is present… ▽ More Global particle-in-cell (PIC) simulations of pulsar magnetospheres are performed with a volume, surface and pair production-based plasma injection schemes to systematically investigate the transition between electrosphere and force-free pulsar magnetospheric regimes. A new extension of the PIC code OSIRIS to model pulsar magnetospheres using a two-dimensional axisymmetric spherical grid is presented. The sub-algorithms of the code and thorough benchmarks are presented in detail, including a new first-order current deposition scheme that conserves charge to machine precision. It is shown that all plasma injection schemes produce a range of magnetospheric regimes. Active solutions can be obtained with surface and volume injection schemes when using artificially large plasma injection rates, and with pair production-based plasma injection for sufficiently large separation between kinematic and pair production energy scales. △ Less

Submitted 9 September, 2023; originally announced September 2023.

Comments: 19 pages, 18 figures, submitted to A&A

Seasonal thaws under mid-to-low pressure atmospheres on Early Mars

Authors: Paolo Simonetti, Giovanni Vladilo, Stavro L. Ivanovski, Laura Silva, Lorenzo Biasiotti, Michele Maris, Giuseppe Murante, Erica Bisesi, Sergio Monai

Abstract: Despite decades of scientific research on the subject, the climate of the first 1.5 Gyr of Mars history has not been fully understood yet. Especially challenging is the need to reconcile the presence of liquid water for extended periods of time on the martian surface with the comparatively low insolation received by the planet, a problem which is known as the Faint Young Sun (FYS) Paradox. In this… ▽ More Despite decades of scientific research on the subject, the climate of the first 1.5 Gyr of Mars history has not been fully understood yet. Especially challenging is the need to reconcile the presence of liquid water for extended periods of time on the martian surface with the comparatively low insolation received by the planet, a problem which is known as the Faint Young Sun (FYS) Paradox. In this paper we use ESTM, a latitudinal energy balance model with enhanced prescriptions for meridional heat diffusion, and the radiative transfer code EOS to investigate how seasonal variations of temperature can give rise to local conditions which are conductive to liquid water runoffs. We include the effects of the martian dichotomy, a northern ocean with either 150 or 550 m of Global Equivalent Layer (GEL) and simplified CO$_2$ or H$_2$O clouds. We find that 1.3-to-2.0 bar CO$_2$-dominated atmospheres can produce seasonal thaws due to inefficient heat redistribution, provided that the eccentricity and the obliquity of the planet are sufficiently different from zero. We also studied the impact of different values for the argument of perihelion. When local favorable conditions exist, they nearly always persist for $>15\%$ of the martian year. These results are obtained without the need for additional greenhouse gases (e.g. H$_2$, CH$_4$) or transient heat-injecting phenomena (e.g. asteroid impacts, volcanic eruptions). Moderate amounts (0.1 to 1\%) of CH$_4$ significantly widens the parameter space region in which seasonal thaws are possible. △ Less

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

Comments: Second and final version, 27 pages, 10 figures, accepted for publication in ApJ

Photon rings as tests for alternative spherically symmetric geometries with thin accretion disks

Authors: Luís F. Dias da Silva, Francisco S. N. Lobo, Gonzalo J. Olmo, Diego Rubiera-Garcia

Abstract: The imaging by the Event Horizon Telescope (EHT) of the supermassive central objects at the heart of the M87 and Milky Way (Sgr A$^\star$) galaxies, has marked the first step into peering at the photon rings and central brightness depression that characterize the optical appearance of black holes surrounded by an accretion disk. Recently, Vagnozzi et. al. [S.~Vagnozzi, \textit{et al.} arXiv:2205.0… ▽ More The imaging by the Event Horizon Telescope (EHT) of the supermassive central objects at the heart of the M87 and Milky Way (Sgr A$^\star$) galaxies, has marked the first step into peering at the photon rings and central brightness depression that characterize the optical appearance of black holes surrounded by an accretion disk. Recently, Vagnozzi et. al. [S.~Vagnozzi, \textit{et al.} arXiv:2205.07787 [gr-qc]] used the claim by the EHT that the size of the {\it shadow} of Sgr A$^\star$ can be inferred by calibrated measurements of the bright ring enclosing it, to constrain a large number of spherically symmetric space-time geometries. In this work we use this result to study some features of the first and second photon rings of a restricted pool of such geometries in thin accretion disk settings. The emission profile of the latter is described by calling upon three analytic samples belonging to the family introduced by Gralla, Lupsasca and Marrone, in order to characterize such photon rings using the Lyapunov exponent of nearly bound orbits and discuss its correlation with the luminosity extinction rate between the first and second photon rings. We finally elaborate on the chances of using such photon rings as observational discriminators of alternative black hole geometries using very long baseline interferometry. △ Less

Submitted 9 October, 2023; v1 submitted 13 July, 2023; originally announced July 2023.

Comments: 17 pages, 58 figures/images; v2: many comments and additions, version accepted for publication in Physical Review D

Journal ref: Phys. Rev. D 108, 084055 (2023)

Exploring the Evolution of Massive Clumps in Simulations that Reproduce the Observed Milky Way α-element Abundance Bimodality

Authors: Bethany R. Garver, David L. Nidever, Victor P. Debattista, Leandro Beraldo e Silva, Tigran Khachaturyants

Abstract: The Milky Way stellar disk has both a thin and a thick component. The thin disk is composed mostly of younger stars ($\lesssim$8 Gyr) with a lower abundance of $α$ elements, while the thick disk contains predominantly older stars ($\gtrsim$8--12 Gyr) with a higher $α$ abundance, giving rise to an $α$-bimodality most prominent at intermediate metallicities. A proposed explanation for the bimodality… ▽ More The Milky Way stellar disk has both a thin and a thick component. The thin disk is composed mostly of younger stars ($\lesssim$8 Gyr) with a lower abundance of $α$ elements, while the thick disk contains predominantly older stars ($\gtrsim$8--12 Gyr) with a higher $α$ abundance, giving rise to an $α$-bimodality most prominent at intermediate metallicities. A proposed explanation for the bimodality is an episode of clumpy star formation, where high-$α$ stars form in massive clumps that appear in the first few Gyrs of the Milky Way's evolution, while low-$α$ stars form throughout the disk and over a longer time span. To better understand the evolution of clumps, we track them and their constituent stars in two clumpy Milky Way simulations that reproduce the $α$-abundance bimodality, one with 10% and the other with 20% supernova feedback efficiency. We investigate the paths that these clumps take in the chemical space ([O/Fe]--[Fe/H]) as well as their mass, star formation rate (SFR), formation location, lifetime, and merger history. The clumps in the simulation with lower feedback last longer on average, with several lasting hundreds of Myr. Some of the clumps do not reach high-$α$, but the ones that do on average had a higher SFR, longer lifetime, greater mass, and form closer to the galactic center than the ones that do not. Most clumps that reach high-$α$ merge with others and eventually spiral into the galactic center, but shed stars along the way to form most of the thick disk component. △ Less

Submitted 21 June, 2023; originally announced June 2023.

Comments: 14 pages, 14 figures, accepted for publication in ApJ

RomAndromeda: The Roman Survey of the Andromeda Halo

Authors: Arjun Dey, Joan Najita, Carrie Filion, Jiwon Jesse Han, Sarah Pearson, Rosemary Wyse, Adrien C. R. Thob, Borja Anguiano, Miranda Apfel, Magda Arnaboldi, Eric F. Bell, Leandro Beraldo e Silva, Gurtina Besla, Aparajito Bhattacharya, Souradeep Bhattacharya, Vedant Chandra, Yumi Choi, Michelle L. M. Collins, Emily C. Cunningham, Julianne J. Dalcanton, Ivanna Escala, Hayden R. Foote, Annette M. N. Ferguson, Benjamin J. Gibson, Oleg Y. Gnedin , et al. (28 additional authors not shown)

Abstract: As our nearest large neighbor, the Andromeda Galaxy provides a unique laboratory for investigating galaxy formation and the distribution and substructure properties of dark matter in a Milky Way-like galaxy. Here, we propose an initial 2-epoch ($Δt\approx 5$yr), 2-band Roman survey of the entire halo of Andromeda, covering 500 square degrees, which will detect nearly every red giant star in the ha… ▽ More As our nearest large neighbor, the Andromeda Galaxy provides a unique laboratory for investigating galaxy formation and the distribution and substructure properties of dark matter in a Milky Way-like galaxy. Here, we propose an initial 2-epoch ($Δt\approx 5$yr), 2-band Roman survey of the entire halo of Andromeda, covering 500 square degrees, which will detect nearly every red giant star in the halo (10$σ$ detection in F146, F062 of 26.5, 26.1AB mag respectively) and yield proper motions to $\sim$25 microarcsec/year (i.e., $\sim$90 km/s) for all stars brighter than F146 $\approx 23.6$ AB mag (i.e., reaching the red clump stars in the Andromeda halo). This survey will yield (through averaging) high-fidelity proper motions for all satellites and compact substructures in the Andromeda halo and will enable statistical searches for clusters in chemo-dynamical space. Adding a third epoch during the extended mission will improve these proper motions by $\sim t^{-1.5}$, to $\approx 11$ km/s, but this requires obtaining the first epoch in Year 1 of Roman operations. In combination with ongoing and imminent spectroscopic campaigns with ground-based telescopes, this Roman survey has the potential to yield full 3-d space motions of $>$100,000 stars in the Andromeda halo, including (by combining individual measurements) robust space motions of its entire globular cluster and most of its dwarf galaxy satellite populations. It will also identify high-velocity stars in Andromeda, providing unique information on the processes that create this population. These data offer a unique opportunity to study the immigration history, halo formation, and underlying dark matter scaffolding of a galaxy other than our own. △ Less

Submitted 21 June, 2023; originally announced June 2023.

Comments: Submitted in response to the call for Roman Space Telescope Core Community Survey white papers

NANCY: Next-generation All-sky Near-infrared Community surveY

Authors: Jiwon Jesse Han, Arjun Dey, Adrian M. Price-Whelan, Joan Najita, Edward F. Schlafly, Andrew Saydjari, Risa H. Wechsler, Ana Bonaca, David J Schlegel, Charlie Conroy, Anand Raichoor, Alex Drlica-Wagner, Juna A. Kollmeier, Sergey E. Koposov, Gurtina Besla, Hans-Walter Rix, Alyssa Goodman, Douglas Finkbeiner, Abhijeet Anand, Matthew Ashby, Benedict Bahr-Kalus, Rachel Beaton, Jayashree Behera, Eric F. Bell, Eric C Bellm , et al. (184 additional authors not shown)

Abstract: The Nancy Grace Roman Space Telescope is capable of delivering an unprecedented all-sky, high-spatial resolution, multi-epoch infrared map to the astronomical community. This opportunity arises in the midst of numerous ground- and space-based surveys that will provide extensive spectroscopy and imaging together covering the entire sky (such as Rubin/LSST, Euclid, UNIONS, SPHEREx, DESI, SDSS-V, GAL… ▽ More The Nancy Grace Roman Space Telescope is capable of delivering an unprecedented all-sky, high-spatial resolution, multi-epoch infrared map to the astronomical community. This opportunity arises in the midst of numerous ground- and space-based surveys that will provide extensive spectroscopy and imaging together covering the entire sky (such as Rubin/LSST, Euclid, UNIONS, SPHEREx, DESI, SDSS-V, GALAH, 4MOST, WEAVE, MOONS, PFS, UVEX, NEO Surveyor, etc.). Roman can uniquely provide uniform high-spatial-resolution (~0.1 arcsec) imaging over the entire sky, vastly expanding the science reach and precision of all of these near-term and future surveys. This imaging will not only enhance other surveys, but also facilitate completely new science. By imaging the full sky over two epochs, Roman can measure the proper motions for stars across the entire Milky Way, probing 100 times fainter than Gaia out to the very edge of the Galaxy. Here, we propose NANCY: a completely public, all-sky survey that will create a high-value legacy dataset benefiting innumerable ongoing and forthcoming studies of the universe. NANCY is a pure expression of Roman's potential: it images the entire sky, at high spatial resolution, in a broad infrared bandpass that collects as many photons as possible. The majority of all ongoing astronomical surveys would benefit from incorporating observations of NANCY into their analyses, whether these surveys focus on nearby stars, the Milky Way, near-field cosmology, or the broader universe. △ Less

Submitted 20 June, 2023; originally announced June 2023.

Comments: Submitted to the call for white papers for the Roman Core Community Survey (June 16th, 2023), and to the Bulletin of the AAS

Early Growing Supermassive Black Holes Strengthen Bars and Boxy/Peanut Bulges

Authors: Vance Wheeler, Monica Valluri, Leandro Beraldo e Silva, Shashank Dattathri, Victor P. Debattista

Abstract: Using N-body simulations we explore the effects of growing a supermassive black hole (SMBH) prior to or during the formation of a stellar bar. Keeping the final mass and growth rate of the SMBH fixed, we show that if it is introduced before or while the bar is still growing, the SMBH does not cause a decrease in bar amplitude. Rather, in most cases, it is strengthened. In addition early growing SM… ▽ More Using N-body simulations we explore the effects of growing a supermassive black hole (SMBH) prior to or during the formation of a stellar bar. Keeping the final mass and growth rate of the SMBH fixed, we show that if it is introduced before or while the bar is still growing, the SMBH does not cause a decrease in bar amplitude. Rather, in most cases, it is strengthened. In addition early growing SMBHs always either decreases the buckling amplitude, delay buckling, or both. This weakening of buckling is caused by an increase in the disk vertical velocity dispersion at radii well beyond the nominal black hole sphere-of-influence. While we find considerable stochasticity and sensitivity to initial conditions, the only case where the SMBH causes a decrease in bar amplitude is when it is introduced after the bar has attained a steady state. In this case we confirm previous findings that the decrease in bar strength is a result of scattering of bar-supporting orbits with small pericenter radii. By heating the inner disk both radially and vertically, an early growing SMBH increases the fraction of stars that can be captured by the Inner Lindblad Resonance (ILR) and the vertical ILR, thereby strengthening both the bar and the boxy peanut shaped bulge. Using orbital frequency analysis of star particles, we show that when an SMBH is introduced early and the bar forms around it, the bar is populated by different families of regular bar-supporting orbits than when the bar forms without an SMBH. △ Less

Submitted 16 November, 2023; v1 submitted 12 June, 2023; originally announced June 2023.

Comments: 26 pages, 15 figures. Updated figures and some arguments in Sec. 5.3,5.4, Results Unchanged. Updates reflect format and style of final published version

Journal ref: ApJ 958 119 (2023)

The Early Data Release of the Dark Energy Spectroscopic Instrument

Authors: DESI Collaboration, A. G. Adame, J. Aguilar, S. Ahlen, S. Alam, G. Aldering, D. M. Alexander, R. Alfarsy, C. Allende Prieto, M. Alvarez, O. Alves, A. Anand, F. Andrade-Oliveira, E. Armengaud, J. Asorey, S. Avila, A. Aviles, S. Bailey, A. Balaguera-Antolínez, O. Ballester, C. Baltay, A. Bault, J. Bautista, J. Behera, S. F. Beltran , et al. (240 additional authors not shown)

Abstract: The Dark Energy Spectroscopic Instrument (DESI) completed its five-month Survey Validation in May 2021. Spectra of stellar and extragalactic targets from Survey Validation constitute the first major data sample from the DESI survey. This paper describes the public release of those spectra, the catalogs of derived properties, and the intermediate data products. In total, the public release includes… ▽ More The Dark Energy Spectroscopic Instrument (DESI) completed its five-month Survey Validation in May 2021. Spectra of stellar and extragalactic targets from Survey Validation constitute the first major data sample from the DESI survey. This paper describes the public release of those spectra, the catalogs of derived properties, and the intermediate data products. In total, the public release includes good-quality spectral information from 466,447 objects targeted as part of the Milky Way Survey, 428,758 as part of the Bright Galaxy Survey, 227,318 as part of the Luminous Red Galaxy sample, 437,664 as part of the Emission Line Galaxy sample, and 76,079 as part of the Quasar sample. In addition, the release includes spectral information from 137,148 objects that expand the scope beyond the primary samples as part of a series of secondary programs. Here, we describe the spectral data, data quality, data products, Large-Scale Structure science catalogs, access to the data, and references that provide relevant background to using these spectra. △ Less

Submitted 15 June, 2023; v1 submitted 9 June, 2023; originally announced June 2023.

Comments: 43 pages, 7 figures, 17 tables, submitted to AJ, DESI EDR references added

Validation of the Scientific Program for the Dark Energy Spectroscopic Instrument

Authors: DESI Collaboration, A. G. Adame, J. Aguilar, S. Ahlen, S. Alam, G. Aldering, D. M. Alexander, R. Alfarsy, C. Allende Prieto, M. Alvarez, O. Alves, A. Anand, F. Andrade-Oliveira, E. Armengaud, J. Asorey, S. Avila, A. Aviles, S. Bailey, A. Balaguera-Antolínez, O. Ballester, C. Baltay, A. Bault, J. Bautista, J. Behera, S. F. Beltran , et al. (239 additional authors not shown)

Abstract: The Dark Energy Spectroscopic Instrument (DESI) was designed to conduct a survey covering 14,000 deg$^2$ over five years to constrain the cosmic expansion history through precise measurements of Baryon Acoustic Oscillations (BAO). The scientific program for DESI was evaluated during a five month Survey Validation (SV) campaign before beginning full operations. This program produced deep spectra of… ▽ More The Dark Energy Spectroscopic Instrument (DESI) was designed to conduct a survey covering 14,000 deg$^2$ over five years to constrain the cosmic expansion history through precise measurements of Baryon Acoustic Oscillations (BAO). The scientific program for DESI was evaluated during a five month Survey Validation (SV) campaign before beginning full operations. This program produced deep spectra of tens of thousands of objects from each of the stellar (MWS), bright galaxy (BGS), luminous red galaxy (LRG), emission line galaxy (ELG), and quasar target classes. These SV spectra were used to optimize redshift distributions, characterize exposure times, determine calibration procedures, and assess observational overheads for the five-year program. In this paper, we present the final target selection algorithms, redshift distributions, and projected cosmology constraints resulting from those studies. We also present a `One-Percent survey' conducted at the conclusion of Survey Validation covering 140 deg$^2$ using the final target selection algorithms with exposures of a depth typical of the main survey. The Survey Validation indicates that DESI will be able to complete the full 14,000 deg$^2$ program with spectroscopically-confirmed targets from the MWS, BGS, LRG, ELG, and quasar programs with total sample sizes of 7.2, 13.8, 7.46, 15.7, and 2.87 million, respectively. These samples will allow exploration of the Milky Way halo, clustering on all scales, and BAO measurements with a statistical precision of 0.28% over the redshift interval $z<1.1$, 0.39% over the redshift interval $1.1<z<1.9$, and 0.46% over the redshift interval $1.9<z<3.5$. △ Less

Submitted 12 January, 2024; v1 submitted 9 June, 2023; originally announced June 2023.

Comments: 42 pages, 18 figures, accepted by AJ

Chemodynamical Properties and Ages of Metal-Poor Stars in S-PLUS

Authors: Felipe Almeida-Fernandes, Vinicius Placco, Helio Rocha-Pinto, Marcelo Borges Fernandes, Guilherme Limberg, Leandro Beraldo e Silva, João A. S. Amarante, Hélio Perottoni, Roderik Overzier, William Schoenell, Tiago Ribeiro, Antonio Kanaan, Claudia Mendes de Oliveira

Abstract: Metal-poor stars are key to our understanding of the early stages of chemical evolution in the Universe. New multi-filter surveys, such as the Southern Photometric Local Universe Survey (S-PLUS), are greatly advancing our ability to select low-metallicity stars. In this work, we analyse the chemodynamical properties and ages of 522 metal-poor candidates selected from the S-PLUS data release 3. Abo… ▽ More Metal-poor stars are key to our understanding of the early stages of chemical evolution in the Universe. New multi-filter surveys, such as the Southern Photometric Local Universe Survey (S-PLUS), are greatly advancing our ability to select low-metallicity stars. In this work, we analyse the chemodynamical properties and ages of 522 metal-poor candidates selected from the S-PLUS data release 3. About 92% of these stars were confirmed to be metal-poor ([Fe/H] $\leq -1$) based on previous medium-resolution spectroscopy. We calculated the dynamical properties of a subsample containing 241 stars, using the astrometry from Gaia Data Release 3. Stellar ages are estimated by a Bayesian isochronal method formalized in this work. We analyse the metallicity distribution of these metal-poor candidates separated into different subgroups of total velocity, dynamical properties, and ages. Our results are used to propose further restrictions to optimize the selection of metal-poor candidates in S-PLUS. The proposed astrometric selection ($\mathrm{parallax}>0.85$ mas) is the one that returns the highest fraction of extremely metal-poor stars (16.3% have [Fe/H] $\leq -3$); the combined selection provides the highest fraction of very metal-poor stars (91.0% have [Fe/H] $\leq -2$), whereas the dynamical selection (eccentricity > 0.35 and diskness < 0.75) is better for targetting metal-poor (99.5% have [Fe/H] $\leq -1$). Using only S-PLUS photometric selections, it is possible to achieve selection fractions of 15.6%, 88.5% and 98.3% for metallicities below $-$3, $-$2 and $-$1, respectively. We also show that it is possible to use S-PLUS to target metal-poor stars in halo substructures such as Gaia-Sausage/Enceladus, Sequoia, Thamnos and the Helmi stream. △ Less

Submitted 20 May, 2023; originally announced May 2023.

Comments: 18 pages, 13 figures. To be published in MNRAS main journal (accepted 15-may-2023)

High-energy synchrotron flares powered by strongly radiative relativistic magnetic reconnection: 2D and 3D PIC simulations

Authors: K. M. Schoeffler, T. Grismayer, D. Uzdensky, L. O. Silva

Abstract: The time evolution of high-energy synchrotron radiation generated in a relativistic pair plasma energized by reconnection of strong magnetic fields is investigated with two- and three-dimensional (2D and 3D) particle-in-cell (PIC) simulations. The simulations in this 2D/3D comparison study are conducted with the radiative PIC code OSIRIS, which self-consistently accounts for the synchrotron radiat… ▽ More The time evolution of high-energy synchrotron radiation generated in a relativistic pair plasma energized by reconnection of strong magnetic fields is investigated with two- and three-dimensional (2D and 3D) particle-in-cell (PIC) simulations. The simulations in this 2D/3D comparison study are conducted with the radiative PIC code OSIRIS, which self-consistently accounts for the synchrotron radiation reaction on the emitting particles, and enables us to explore the effects of synchrotron cooling. Magnetic reconnection causes compression of the plasma and magnetic field deep inside magnetic islands (plasmoids), leading to an enhancement of the flaring emission, which may help explain some astrophysical gamma-ray flare observations. Although radiative cooling weakens the emission from plasmoid cores, it facilitates additional compression there, further amplifying the magnetic field $B$ and plasma density~$n$, and thus partially mitigating this effect. Novel simulation diagnostics utilizing 2D histograms in the $n\mbox{-}B$ space are developed and used to visualize and quantify the effects of compression. The $n\mbox{-}B$ histograms are observed to be bounded by relatively sharp power-law boundaries marking clear limits on compression. Theoretical explanations for some of these compression limits are developed, rooted in radiative resistivity or 3D kinking instabilities. Systematic parameter-space studies with respect to guide magnetic field, system size, and upstream magnetization are conducted and suggest that stronger compression, brighter high-energy radiation, and perhaps significant quantum electrodynamic (QED) effects such as pair production, may occur in environments with larger reconnection-region sizes and higher magnetization, particularly when magnetic field strengths approach the critical (Schwinger) field, as found in magnetar magnetospheres. △ Less

Submitted 29 March, 2023; originally announced March 2023.

Comments: 31 pages, 23 figures

The imprint of clump formation at high redshift. II. The chemistry of the bulge

Authors: Victor P. Debattista, David J. Liddicott, Oscar A. Gonzalez, Leandro Beraldo e Silva, Joao A. S. Amarante, Ilin Lazar, Manuela Zoccali, Elena Valenti, Deanne B. Fisher, Tigran Khachaturyants, David L. Nidever, Thomas R. Quinn, Min Du, Susan Kassin

Abstract: In Paper I we showed that clumps in high-redshift galaxies, having a high star formation rate density (Σ_SFR), produce disks with two tracks in the [Fe/H]-[α/Fe] chemical space, similar to that of the Milky Way's (MW's) thin + thick disks. Here we investigate the effect of clumps on the bulge's chemistry. The chemistry of the MW's bulge is comprised of a single track with two density peaks separat… ▽ More In Paper I we showed that clumps in high-redshift galaxies, having a high star formation rate density (Σ_SFR), produce disks with two tracks in the [Fe/H]-[α/Fe] chemical space, similar to that of the Milky Way's (MW's) thin + thick disks. Here we investigate the effect of clumps on the bulge's chemistry. The chemistry of the MW's bulge is comprised of a single track with two density peaks separated by a trough. We show that the bulge chemistry of an N-body + smoothed particle hydrodynamics clumpy simulation also has a single track. Star formation within the bulge is itself in the high-Σ_SFR clumpy mode, which ensures that the bulge's chemical track follows that of the thick disk at low [Fe/H] and then extends to high [Fe/H], where it peaks. The peak at low metallicity instead is comprised of a mixture of in-situ stars and stars accreted via clumps. As a result, the trough between the peaks occurs at the end of the thick disk track. We find that the high-metallicity peak dominates near the mid-plane and declines in relative importance with height, as in the MW. The bulge is already rapidly rotating by the end of the clump epoch, with higher rotation at low [α/Fe]. Thus clumpy star formation is able to simultaneously explain the chemodynamic trends of the MW's bulge, thin + thick disks and the Splash. △ Less

Submitted 14 March, 2023; originally announced March 2023.

Comments: 25 pages, in press at ApJ

Orbital support and evolution of flat profiles of bars (shoulders)

Authors: Leandro Beraldo e Silva, Victor P. Debattista, Stuart R. Anderson, Monica Valluri, Peter Erwin, Kathryne J. Daniel, Nathan Deg

Abstract: Many barred galaxies exhibit upturns (shoulders) in their bar major-axis density profile. Simulation studies have suggested that shoulders are supported by looped $x_1$ orbits, occur in growing bars, and can appear after bar-buckling. We investigate the orbital support and evolution of shoulders via frequency analyses of orbits in simulations. We confirm that looped orbits are shoulder-supporting,… ▽ More Many barred galaxies exhibit upturns (shoulders) in their bar major-axis density profile. Simulation studies have suggested that shoulders are supported by looped $x_1$ orbits, occur in growing bars, and can appear after bar-buckling. We investigate the orbital support and evolution of shoulders via frequency analyses of orbits in simulations. We confirm that looped orbits are shoulder-supporting, and can remain so, to a lesser extent, after being vertically thickened. We show that looped orbits appear at the resonance $(Ω_\varphi - Ω_\mathrm{P})/Ω_R=1/2$ (analogous to the classical Inner Lindblad Resonance, and here called ILR) with vertical-to-radial frequency ratios $1 \lesssimΩ_z/Ω_R \lesssim 3/2$ (vertically warm orbits). Cool orbits at the ILR (those with $Ω_z/Ω_R > 3/2$) are vertically thin and have no loops, contributing negligibly to shoulders. As bars slow and thicken, either secularly or by buckling, they populate warm orbits at the ILR. Further thickening carries these orbits towards crossing the vertical ILR [vILR, $(Ω_\varphi - Ω_\mathrm{P})/Ω_z=1/2$], where they convert in-plane to vertical motion, become chaotic, kinematically hotter and less shoulder-supporting. Hence, persistent shoulders require bars to trap new stars, consistent with the need for a growing bar. Since buckling speeds up trapping on warm orbits at the ILR, it can be followed by shoulder formation, as seen in simulations. This sequence supports the recent observational finding that shoulders likely precede the emergence of BP-bulges. The python module for the frequency analysis, naif, is made available. △ Less

Submitted 14 September, 2023; v1 submitted 8 March, 2023; originally announced March 2023.

Comments: Published in ApJ; A few additional figures in response to referee's request. Conclusions unchanged

Journal ref: ApJ 955 38 (2023)

Magnetic frame-dragging correction to the electromagnetic solution of a compact neutron star

Authors: R. Torres, T. Grismayer, F. Cruz, L. O. Silva

Abstract: Neutron stars are usually modelled as spherical, rotating perfect conductors with a predominant intrinsic dipolar magnetic field anchored to their stellar crust. Due to their compactness, General Relativity corrections must be accounted for in Maxwell's equations, leading to modified interior and exterior electromagnetic solutions. We present analytical solutions for slowly-rotating magnetised neu… ▽ More Neutron stars are usually modelled as spherical, rotating perfect conductors with a predominant intrinsic dipolar magnetic field anchored to their stellar crust. Due to their compactness, General Relativity corrections must be accounted for in Maxwell's equations, leading to modified interior and exterior electromagnetic solutions. We present analytical solutions for slowly-rotating magnetised neutron stars taking into account the magnetic frame-dragging correction. For typical compactness values, i.e. $R_s \sim 0.5 [R_*]$, we show that the new terms lead to a percent order correction in the magnetic field orientation and strength compared to the case with no magnetic frame-dragging correction. Also, we obtain a self-consistent redistribution of the surface azimuthal current. We verify the validity of the derived solution through two-dimensional particle-in-cell simulations of an isolated neutron star. Defining the azimuthal electric and magnetic field amplitudes during the transient phase as observables, we prove that the magnetic frame-dragging correction reduces the transient wave amplitude, as expected from the analytical solution. We show that simulations are more accurate and stable when we include all first-order terms. The increased accuracy at lower spatiotemporal resolutions translates into a reduction in simulation runtimes. △ Less

Submitted 7 January, 2024; v1 submitted 20 February, 2023; originally announced February 2023.

Comments: 12 pages, 8 figures

Journal ref: MNRAS, Volume 524, Issue 3, September 2023, Pages 4116-4127

Climate and atmospheric models of rocky planets: habitability and observational properties

Authors: L. Silva, R. Bevilacqua, L. Biasiotti, E. Bisesi, S. L. Ivanovski, M. Maris, S. Monai, G. Murante, P. Simonetti, G. Vladilo

Abstract: The quest for atmospheric spectral signatures that may witness biological activity in exoplanets is focused on rocky planets. The best targets for future, challenging spectroscopic observations will be selected among potentially habitable planets. Surface habitability can be quantified and explored with climate and atmospheric models according to temperature-based criteria. The conceptual, modelli… ▽ More The quest for atmospheric spectral signatures that may witness biological activity in exoplanets is focused on rocky planets. The best targets for future, challenging spectroscopic observations will be selected among potentially habitable planets. Surface habitability can be quantified and explored with climate and atmospheric models according to temperature-based criteria. The conceptual, modellistic, technological and interpretative complexity of the problem requires to develop flexible climate and atmospheric models suited for a comprehensive exploration of observationally unconstrained parameters, and to simulate and interpret definitely non-terrestrial conditions. We present a summary and preliminary results on the work we are performing on multi-parametric explorations of the habitability and observational properties of rocky planets. △ Less

Submitted 22 January, 2023; originally announced January 2023.

Comments: to appear on MemSAIt, vol 94. Proceedings of the Hack100 Conference: Past, Present and Future of Astrophysical Spectroscopy, 6-10 June 2022, Trieste, Italy

Radiation reaction cooling as a source of anisotropic momentum distributions with inverted populations

Authors: Pablo J. Bilbao, Luis O. Silva

Abstract: Under the presence of strong electromagnetic fields and radiation reaction, plasmas develop anisotropic momentum distributions, characterized by a population inversion. This is a general property of collisionless plasmas when the radiation reaction force is taken into account. We study the case of a plasma in a strong magnetic field and demonstrate the development of ring momentum distributions. T… ▽ More Under the presence of strong electromagnetic fields and radiation reaction, plasmas develop anisotropic momentum distributions, characterized by a population inversion. This is a general property of collisionless plasmas when the radiation reaction force is taken into account. We study the case of a plasma in a strong magnetic field and demonstrate the development of ring momentum distributions. The timescales for ring formation are derived for this configuration. The analytical results for the ring properties and the timescales for ring formation are confirmed with particle-in-cell simulations. The resulting momentum distributions are kinetically unstable and are known to lead to coherent radiation emission in astrophysical plasmas and laboratory setups. △ Less

Submitted 29 March, 2023; v1 submitted 23 December, 2022; originally announced December 2022.

Comments: Accepted for publication in PRL. 8 pages, 4 figures

Journal ref: Phys. Rev. Lett. 130, 165101 (2023)

Conservative Evolution of Black Hole Perturbations with Time-Symmetric Numerical Methods

Authors: Michael F. O'Boyle, Charalampos Markakis, Lidia J. Gomes Da Silva, Rodrigo Panosso Macedo, Juan A. Valiente Kroon

Abstract: The scheduled launch of the LISA Mission in the next decade has called attention to the gravitational self-force problem. Despite an extensive body of theoretical work, long-time numerical computations of gravitational waves from extreme-mass-ratio-inspirals remain challenging. This work proposes a class of numerical evolution schemes suitable to this problem based on Hermite integration. Their mo… ▽ More The scheduled launch of the LISA Mission in the next decade has called attention to the gravitational self-force problem. Despite an extensive body of theoretical work, long-time numerical computations of gravitational waves from extreme-mass-ratio-inspirals remain challenging. This work proposes a class of numerical evolution schemes suitable to this problem based on Hermite integration. Their most important feature is time-reversal symmetry and unconditional stability, which enables these methods to preserve symplectic structure, energy, momentum and other Noether charges over long time periods. We apply Noether's theorem to the master fields of black hole perturbation theory on a hyperboloidal slice of Schwarzschild spacetime to show that there exist constants of evolution that numerical simulations must preserve. We demonstrate that time-symmetric integration schemes based on a 2-point Taylor expansion (such as Hermite integration) numerically conserve these quantities, unlike schemes based on a 1-point Taylor expansion (such as Runge-Kutta). This makes time-symmetric schemes ideal for long-time EMRI simulations. △ Less

Submitted 5 October, 2022; originally announced October 2022.

Comments: 43 pages, 8 figures

MSC Class: 49S05; 49K20; 65M22; 65M25; 65M70; 70S10; 83-10; 83C25

The pattern speeds of vertical breathing waves

Authors: Tigran Khachaturyants, Victor P. Debattista, Soumavo Ghosh, Leandro Beraldo e Silva, Kathryne J. Daniel

Abstract: We measure and compare the pattern speeds of vertical breathing, vertical bending, and spiral density waves in two isolated N-body+SPH simulations, using windowed Fourier transforms over 1 Gyr time intervals. We show that the pattern speeds of the breathing waves match those of the spirals but are different from those of the bending waves. We also observe matching pattern speeds between the bar an… ▽ More We measure and compare the pattern speeds of vertical breathing, vertical bending, and spiral density waves in two isolated N-body+SPH simulations, using windowed Fourier transforms over 1 Gyr time intervals. We show that the pattern speeds of the breathing waves match those of the spirals but are different from those of the bending waves. We also observe matching pattern speeds between the bar and breathing waves. Our results not only strengthen the case that, throughout the disc, breathing motions are driven by spirals but indeed that the breathing motions are part and parcel of the spirals. △ Less

Submitted 1 October, 2022; originally announced October 2022.

Comments: 5 pages, 4 figures, published in MNRAS Letters

Journal ref: Monthly Notices of the Royal Astronomical Society: Letters, Volume 517, Issue 1, November 2022, Pages L55-L59

DESI Observations of the Andromeda Galaxy: Revealing the Immigration History of our Nearest Neighbor

Authors: Arjun Dey, Joan R. Najita, S. E. Koposov, J. Josephy-Zack, Gabriel Maxemin, Eric F. Bell, C. Poppett, E. Patel, L. Beraldo e Silva, A. Raichoor, D. Schlegel, D. Lang, A. Meisner, Adam D. Myers, J. Aguilar, S. Ahlen, C. Allende Prieto, D. Brooks, A. P. Cooper, K. S. Dawson, A. de la Macorra, P. Doel, A. Font-Ribera, Juan Garcia-Bellido, S. Gontcho A Gontcho , et al. (23 additional authors not shown)

Abstract: We present DESI observations of the inner halo of M31, which reveal the kinematics of a recent merger - a galactic immigration event - in exquisite detail. Of the 11,416 sources studied in 3.75 hour of on-sky exposure time, 7,438 are M31 sources with well measured radial velocities. The observations reveal intricate coherent kinematic structure in the positions and velocities of individual stars:… ▽ More We present DESI observations of the inner halo of M31, which reveal the kinematics of a recent merger - a galactic immigration event - in exquisite detail. Of the 11,416 sources studied in 3.75 hour of on-sky exposure time, 7,438 are M31 sources with well measured radial velocities. The observations reveal intricate coherent kinematic structure in the positions and velocities of individual stars: streams, wedges, and chevrons. While hints of coherent structures have been previously detected in M31, this is the first time they have been seen with such detail and clarity in a galaxy beyond the Milky Way. We find clear kinematic evidence for shell structures in the Giant Stellar Stream, the Northeast Shelf and Western Shelf regions. The kinematics are remarkably similar to the predictions of dynamical models constructed to explain the spatial morphology of the inner halo. The results are consistent with the interpretation that much of the substructure in the inner halo of M31 is produced by a single galactic immigration event 1 - 2 Gyr ago. Significant numbers of metal-rich stars ([Fe/H]$>-0.5$) are present in all of the detected substructures, suggesting that the immigrating galaxy had an extended star formation history. We also investigate the ability of the shells and Giant Stellar Stream to constrain the gravitational potential of M31, and estimate the mass within a projected radius of 125 kpc to be ${\rm log_{10}}\, M_{\rm NFW}(<125\,{\rm kpc})/M_\odot = 11.80_{-0.10}^{+0.12}$. The results herald a new era in our ability to study stars on a galactic scale and the immigration histories of galaxies. △ Less

Submitted 20 January, 2023; v1 submitted 24 August, 2022; originally announced August 2022.

Comments: 45 pages, 22 figures, 8 tables; Astrophysical Journal in press; Data at https://zenodo.org/record/6977494

Overview of the DESI Milky Way Survey

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

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

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

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

EOS-ESTM: A flexible climate model for habitable exoplanets

Authors: L. Biasiotti, P. Simonetti, G. Vladilo, L. Silva, G. Murante, S. Ivanovski, M. Maris, S. Monai, E. Bisesi, J. von Hardenberg, A. Provenzale

Abstract: Rocky planets with temperate conditions provide the best chance for discovering habitable worlds and life outside the Solar System. In the last decades, new instrumental facilities and large observational campaigns have been driven by the quest for habitable worlds. Climate models aimed at studying the habitability of rocky planets are essential tools to pay off these technological and observation… ▽ More Rocky planets with temperate conditions provide the best chance for discovering habitable worlds and life outside the Solar System. In the last decades, new instrumental facilities and large observational campaigns have been driven by the quest for habitable worlds. Climate models aimed at studying the habitability of rocky planets are essential tools to pay off these technological and observational endeavours. In this context, we present EOS-ESTM, a fast and flexible model aimed at exploring the impact on habitability of multiple climate factors, including those unconstrained by observations. EOS-ESTM is built on ESTM, a seasonal-latitudinal energy balance model featuring an advanced treatment of the meridional and vertical transport. The novel features of EOS-ESTM include: (1) parameterizations for simulating the climate impact of oceans, land, ice, and clouds as a function of temperature and stellar zenith distance; (2) a procedure (EOS) for calculating the radiative transfer in atmospheres with terrestrial and non-terrestrial compositions illuminated by solar- and non-solar-type stars. By feeding EOS-ESTM with Earth's stellar, orbital and planetary parameters we derive a reference model that satisfies a large number of observational constraints of the Earth's climate system. Validation tests of non-terrestrial conditions yield predictions that are in line with comparable results obtained with a hierarchy of climate models. The application of EOS-ESTM to planetary atmospheres in maximum greenhouse conditions demonstrates the possibility of tracking the snowball transition at the outer edge of the HZ for a variety of planetary parameters, paving the road for multi-parametric studies of the HZ. △ Less

Submitted 10 June, 2022; originally announced June 2022.

Comments: 22 pages, 14 figures, 10 tables, accepted for publication in MNRAS

The secular growth of bars revealed by flat (peak + shoulders) density profiles

Authors: Stuart Robert Anderson, Victor P. Debattista, Peter Erwin, David J. Liddicott, Nathan Deg, Leandro Beraldo e Silva

Abstract: The major-axis density profiles of bars are known to be either exponential or 'flat'. We develop an automated non-parametric algorithm to detect flat profiles and apply it to a suite of simulations (with and without gas). We demonstrate that flat profiles are a manifestation of a bar's secular growth, producing a 'shoulder' region (an overdensity above an exponential) in its outskirts. Shoulders a… ▽ More The major-axis density profiles of bars are known to be either exponential or 'flat'. We develop an automated non-parametric algorithm to detect flat profiles and apply it to a suite of simulations (with and without gas). We demonstrate that flat profiles are a manifestation of a bar's secular growth, producing a 'shoulder' region (an overdensity above an exponential) in its outskirts. Shoulders are not present when bars form, but develop as the bar grows. If the bar does not grow, shoulders do not form. Shoulders are often accompanied by box/peanut bulges, but develop separately from them and are independent tracers of a bar's growth. They can be observed at a wide range of viewing orientations with only their slope varying significantly with inclination. We present evidence that shoulders are produced by looped x1 orbits. Since the growth rate of the bar moderately correlates with the growth rate of the shoulder strength, these orbits are probably recently trapped. Shoulders therefore are evidence of bar growth. The properties of the shoulders do not, however, establish the age of a bar, because secondary buckling or strong spirals may destroy shoulders, and also because shoulders do not form if the bar does not grow much. In particular, our results show that an exponential profile is not necessarily an indication of a young bar. △ Less

Submitted 13 May, 2022; originally announced May 2022.

Comments: 20 pages, 23 figures, 2 tables. Published in MNRAS

Journal ref: MNRAS 513, 1642-1661 (2022)

The impact of very massive stars on the chemical evolution of extremely metal-poor galaxies

Authors: S. Goswami, L. Silva, A. Bressan, V. Grisoni, G. Costa, P. Marigo, G. L. Granato, A. Lapi, M. Spera

Abstract: In recent observations of extremely metal-poor low-mass starburst galaxies, almost solar Fe/O ratios are reported, despite N/O ratios consistent with the low metallicity. We investigate if the peculiar Fe/O ratios can be a distinctive signature of an early enrichment produced by very massive objects dying as Pair-Instability Supernovae (PISN). We run chemical evolution models with yields that acco… ▽ More In recent observations of extremely metal-poor low-mass starburst galaxies, almost solar Fe/O ratios are reported, despite N/O ratios consistent with the low metallicity. We investigate if the peculiar Fe/O ratios can be a distinctive signature of an early enrichment produced by very massive objects dying as Pair-Instability Supernovae (PISN). We run chemical evolution models with yields that account for the contribution by PISN. We use both the recent non-rotating stellar yields from Goswami et al. 2021, and new yields from rotating very massive stars calculated on purpose in this work. We also search for the best initial mass function (IMF) that is able to reproduce the observations. We can reproduce the observations by adopting a bi-modal IMF and by including an initial burst of rotating very massive stars. Only with a burst of very massive stars can we reproduce the almost solar Fe/O ratios at the estimated young ages. We also confirm that rotation is absolutely needed to concomitantly reproduce the observed N/O ratios. These results stress the importance of very massive stars in galactic chemical evolution studies and strongly support a top-heavy initial mass function in the very early evolutionary stages of metal poor starburst galaxies. △ Less

Submitted 6 May, 2022; originally announced May 2022.

Comments: 12 pages, 6 figures, accepted for publication in A&A

Journal ref: A&A 663, A1 (2022)

GASTRO library I: the simulated chemodynamical properties of several GSE-like stellar halos

Authors: João A. S. Amarante, Victor P. Debattista, Leandro Beraldo e Silva, Chervin F. P. Laporte, Nathan Deg

Abstract: The Milky Way stellar halo contains relics of ancient mergers that tell the story of our Galaxy's formation. Some of them are identified due to their similarity in energy, actions and chemistry, referred to as the "chemodynamical space", and are often attributed to distinct merger events. It is also known that our Galaxy went through a significant merger event that shaped the local stellar halo du… ▽ More The Milky Way stellar halo contains relics of ancient mergers that tell the story of our Galaxy's formation. Some of them are identified due to their similarity in energy, actions and chemistry, referred to as the "chemodynamical space", and are often attributed to distinct merger events. It is also known that our Galaxy went through a significant merger event that shaped the local stellar halo during its first Gyr. Previous studies using $N$-body only and cosmological hydrodynamical simulations have shown that such single massive merger can produce several "signatures" in the chemodynamical space, which can potentially be misinterpreted as distinct merger events. Motivated by these, in this work we use a subset of the GASTRO, library which consists of several SPH+$N$-body models of single accretion event in a Milky Way-like galaxy. Here, we study models with orbital properties similar to the main merger event of our Galaxy and explore the implications to known stellar halo substructures. We find that: $i.$ supernova feedback efficiency influences the satellite's structure and orbital evolution, resulting in distinct chemodynamical features for models with the same initial conditions, $ii.$ very retrograde high energy stars are the most metal-poor of the accreted dwarf galaxy and could be misinterpreted as a distinct merger $iii.$ the most bound stars are more metal-rich in our models, the opposite of what is observed in the Milky Way, suggesting a secondary massive merger, and finally $iv.$ our models can reconcile other known substructures to an unique progenitor. △ Less

Submitted 17 September, 2022; v1 submitted 26 April, 2022; originally announced April 2022.

Comments: Published in ApJ

Journal ref: ApJ 937 12 (2022)

Model of pulsar pair cascades in non uniform electric fields: growth rate, density profile and screening time

Authors: Fábio Cruz, Thomas Grismayer, Simon Iteanu, Paolo Tortone, Luis O. Silva

Abstract: Time-dependent cascades of electron-positron pairs are thought to be the main source of plasma in pulsar magnetospheres and a primary ingredient to explain the nature of pulsar radio emission, a longstanding open problem in high-energy astrophysics. During these cascades - positive feedback loops of gamma-ray photon emission, via curvature radiation by TeV electrons and positrons, and pair product… ▽ More Time-dependent cascades of electron-positron pairs are thought to be the main source of plasma in pulsar magnetospheres and a primary ingredient to explain the nature of pulsar radio emission, a longstanding open problem in high-energy astrophysics. During these cascades - positive feedback loops of gamma-ray photon emission, via curvature radiation by TeV electrons and positrons, and pair production -, the plasma self-consistently develops inductive waves that couple to electromagnetic modes capable of escaping the pulsar dense plasma. In this work, we present an analytical description of pair cascades relevant in pulsars, including their onset, exponential growth and saturation stages. We study this problem in the case of a background linear electric field, relevant in pulsar polar caps, and using an heuristic model of the pair production process. The analytical results are confirmed with particle-in-cell simulations performed with OSIRIS including heuristic pair production. △ Less

Submitted 7 April, 2022; originally announced April 2022.

Comments: 9 pages, 7 figures, accepted for publication in PoP

Variability of the Sun's luminosity places constraints on the thermal equilibrium of the convection zone

Authors: L. E. A. Vieira, G. Kopp, T. Dudok de Wit, L. A. da Silva, F. Carlesso, A. Barbosa, A. Muralikrishna, R. Santos

Abstract: Luminosity, which is the total amount of radiant energy emitted by an object, is one of the most critical quantities in astrophysics for characterizing stars. Equally important is the temporal evolution of a star's luminosity because of its intimate connection with the stellar energy budget, large-scale convective motion, and heat storage in the stellar interior. Here, we model the solar luminosit… ▽ More Luminosity, which is the total amount of radiant energy emitted by an object, is one of the most critical quantities in astrophysics for characterizing stars. Equally important is the temporal evolution of a star's luminosity because of its intimate connection with the stellar energy budget, large-scale convective motion, and heat storage in the stellar interior. Here, we model the solar luminosity by extending a semi-empirical total solar irradiance (TSI) model that uses solar-surface magnetism to reconstruct solar irradiance over the entire 4π solid angle around the Sun. This model was constrained by comparing its output to the irradiance in the Earth's direction with the measured TSI. Comparing the solar luminosity to the TSI on timescales from days to for cycles 23 and 24, we find poor agreement on short timescales (< solar rotation). On longer timescales, however, we find good agreement between the luminosity model and the TSI, which suggests that the extrapolation of luminosities to multi-cycle timescales based on TSI reconstructions may be possible. We show that the solar luminosity is not constant but varies in phase with the solar cycle. This variation has an amplitude of 0.14% from minimum to maximum for solar cycle 23. Considering the energetics in the solar convection zone, it is therefore obvious that a steady-state input from the radiative zone at the solar minimum level would lead to a gradual reduction in the energy content in the convection zone over multi-century timescales. We show that the luminosity at the base of the convection zone should be approximately 0.032% higher than that at the solar surface during solar minimum to maintain net energy equilibrium through the solar cycle. These different energy-input scenarios place constraints on the long-term evolution of the total solar irradiance and its impact on the solar forcing of climate variability. △ Less

Submitted 6 April, 2022; originally announced April 2022.

Comments: Accepted for publication in ApJS. (26.March.2022)

Snowmass2021 Cosmic Frontier White Paper: Prospects for obtaining Dark Matter Constraints with DESI

Authors: Monica Valluri, Solene Chabanier, Vid Irsic, Eric Armengaud, Michael Walther, Connie Rockosi, Miguel A. Sanchez-Conde, Leandro Beraldo e Silva, Andrew P. Cooper, Elise Darragh-Ford, Kyle Dawson, Alis J. Deason, Simone Ferraro, Jaime E. Forero-Romero, Antonella Garzilli, Ting Li, Zarija Lukic, Christopher J. Manser, Nathalie Palanque-Delabrouille, Corentin Ravoux, Ting Tan, Wenting Wang, Risa Wechsler, Andreia Carrillo, Arjun Dey , et al. (7 additional authors not shown)

Abstract: Despite efforts over several decades, direct-detection experiments have not yet led to the discovery of the dark matter (DM) particle. This has led to increasing interest in alternatives to the Lambda CDM (LCDM) paradigm and alternative DM scenarios (including fuzzy DM, warm DM, self-interacting DM, etc.). In many of these scenarios, DM particles cannot be detected directly and constraints on thei… ▽ More Despite efforts over several decades, direct-detection experiments have not yet led to the discovery of the dark matter (DM) particle. This has led to increasing interest in alternatives to the Lambda CDM (LCDM) paradigm and alternative DM scenarios (including fuzzy DM, warm DM, self-interacting DM, etc.). In many of these scenarios, DM particles cannot be detected directly and constraints on their properties can ONLY be arrived at using astrophysical observations. The Dark Energy Spectroscopic Instrument (DESI) is currently one of the most powerful instruments for wide-field surveys. The synergy of DESI with ESA's Gaia satellite and future observing facilities will yield datasets of unprecedented size and coverage that will enable constraints on DM over a wide range of physical and mass scales and across redshifts. DESI will obtain spectra of the Lyman-alpha forest out to z~5 by detecting about 1 million QSO spectra that will put constraints on clustering of the low-density intergalactic gas and DM halos at high redshift. DESI will obtain radial velocities of 10 million stars in the Milky Way (MW) and Local Group satellites enabling us to constrain their global DM distributions, as well as the DM distribution on smaller scales. The paradigm of cosmological structure formation has been extensively tested with simulations. However, the majority of simulations to date have focused on collisionless CDM. Simulations with alternatives to CDM have recently been gaining ground but are still in their infancy. While there are numerous publicly available large-box and zoom-in simulations in the LCDM framework, there are no comparable publicly available WDM, SIDM, FDM simulations. DOE support for a public simulation suite will enable a more cohesive community effort to compare observations from DESI (and other surveys) with numerical predictions and will greatly impact DM science. △ Less

Submitted 1 July, 2022; v1 submitted 14 March, 2022; originally announced March 2022.

Comments: Contributed white paper to Snowmass 2021, CF03; minor revisions

Bending waves excited by irregular gas inflow along warps

Authors: Tigran Khachaturyants, Leandro Beraldo e Silva, Victor P. Debattista, Kathryne J. Daniel

Abstract: Gaia has revealed clear evidence of bending waves in the vertical kinematics of stars in the Solar Neighbourhood. We study bending waves in two simulations, one warped, with the warp due to misaligned gas inflow, and the other unwarped. We find slow, retrograde bending waves in both models, with the ones in the warped model having larger amplitudes. We also find fast, prograde bending waves. Progr… ▽ More Gaia has revealed clear evidence of bending waves in the vertical kinematics of stars in the Solar Neighbourhood. We study bending waves in two simulations, one warped, with the warp due to misaligned gas inflow, and the other unwarped. We find slow, retrograde bending waves in both models, with the ones in the warped model having larger amplitudes. We also find fast, prograde bending waves. Prograde bending waves in the unwarped model are very weak, in agreement with the expectation that these waves should decay on short, ~ crossing, timescales, due to strong winding. However, prograde bending waves are much stronger for the duration of the warped model, pointing to irregular gas inflow along the warp as a continuous source of excitation. We demonstrate that large amplitude bending waves that propagate through the Solar Neighbourhood give rise to a correlation between the mean vertical velocity and the angular momentum, with a slope consistent with that found by Gaia. The bending waves affect populations of all ages, but the sharpest features are found in the young populations, hinting that short wavelength waves are not supported by the older, kinematically hotter, populations. Our results demonstrate the importance of misaligned gas accretion as a recurrent source of vertical perturbations of disc galaxies, including in the Milky Way. △ Less

Submitted 7 March, 2022; originally announced March 2022.

Comments: 21 pages, 20 figures; Published on MNRAS

EOS: Atmospheric Radiative Transfer in Habitable Worlds with HELIOS

Authors: Paolo Simonetti, Giovanni Vladilo, Laura Silva, Michele Maris, Stavro L. Ivanovski, Lorenzo Biasiotti, Matej Malik, Jost von Hardenberg

Abstract: We present EOS, a procedure for determining the Outgoing Longwave Radiation (OLR) and top-of-atmosphere (TOA) albedo for a wide range of conditions expected to be present in the atmospheres of rocky planets with temperate conditions. EOS is based on HELIOS and HELIOS-K, which are novel and publicly available atmospheric radiative transfer (RT) codes optimized for fast calculations with GPU process… ▽ More We present EOS, a procedure for determining the Outgoing Longwave Radiation (OLR) and top-of-atmosphere (TOA) albedo for a wide range of conditions expected to be present in the atmospheres of rocky planets with temperate conditions. EOS is based on HELIOS and HELIOS-K, which are novel and publicly available atmospheric radiative transfer (RT) codes optimized for fast calculations with GPU processors. These codes were originally developed for the study of giant planets. In this paper we present an adaptation for applications to terrestrial-type, habitable planets, adding specific physical recipes for the gas opacity and vertical structure of the atmosphere. To test the reliability of the procedure we assessed the impact of changing line opacity profile, continuum opacity model, atmospheric lapse rate and tropopause position prescriptions on the OLR and the TOA albedo. The results obtained with EOS are in line with those of other RT codes running on traditional CPU processors, while being at least one order of magnitude faster. The adoption of OLR and TOA albedo data generated with EOS in a zonal and seasonal climate model correctly reproduce the fluxes of the present-day Earth measured by the CERES spacecraft. The results of this study disclose the possibility to incorporate fast RT calculations in climate models aimed at characterizing the atmospheres of habitable exoplanets. △ Less

Submitted 29 November, 2021; v1 submitted 22 October, 2021; originally announced October 2021.

Comments: Accepted for publication in ApJ

How stars formed in warps settle into (and contaminate) thick discs

Authors: Tigran Khachaturyants, Leandro Beraldo e Silva, Victor P. Debattista

Abstract: In recent years star formation has been discovered in the Milky Way's warp. These stars formed in the warp (warp stars) must eventually settle into the plane of the disc. We use an $N$-body$+$smooth particle hydrodynamics model of a warped galaxy to study how warp stars settle into the disc. By following warp stars in angular momentum space, we show that they first tilt to partially align with the… ▽ More In recent years star formation has been discovered in the Milky Way's warp. These stars formed in the warp (warp stars) must eventually settle into the plane of the disc. We use an $N$-body$+$smooth particle hydrodynamics model of a warped galaxy to study how warp stars settle into the disc. By following warp stars in angular momentum space, we show that they first tilt to partially align with the main disc in a time scale of $\sim1$ Gyr. Then, once differential precession halts this process, they phase mix into an axisymmetric distribution on a time scale of $\sim 6$ Gyr. The warp stars end up contaminating the geometric thick disc. Because the warp in our fiducial simulation is growing, the {\it warp stars} settle to a distribution with a negative vertical age gradient as younger stars settle further from the mid-plane. While vertically extended, warp star orbits are still nearly circular and they are therefore subject to radial migration, with a net movement inwards. As a result warp stars can be found throughout the disc. The density distribution of a given population of warp stars evolves from a torus to an increasingly centrally filled-in density distribution. Therefore we argue that, in the Milky Way, warp stars should be found in the Solar Neighbourhood. Moreover, settled warp stars may constitute part of the young flaring population seen in the Milky Way's outskirts. △ Less

Submitted 12 October, 2021; originally announced October 2021.

Comments: 21 pages, 25 figures; Published on MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 508, Issue 2, December 2021, Pages 2350-2369

Coherent emission from QED cascades in pulsar polar caps

Authors: Fábio Cruz, Thomas Grismayer, Alexander Y. Chen, Anatoly Spitkovsky, Luis O. Silva

Abstract: Pulsar magnetospheres are thought to be filled with electron-positron plasma generated in pair cascades. The driving mechanism of these cascades is the emission of gamma-ray photons and their conversion into pairs via Quantum Electrodynamics (QED) processes. In this work, we present 2D particle-in-cell simulations of pair cascades in pulsar polar caps with realistic magnetic field geometry that in… ▽ More Pulsar magnetospheres are thought to be filled with electron-positron plasma generated in pair cascades. The driving mechanism of these cascades is the emission of gamma-ray photons and their conversion into pairs via Quantum Electrodynamics (QED) processes. In this work, we present 2D particle-in-cell simulations of pair cascades in pulsar polar caps with realistic magnetic field geometry that include the relevant QED processes from first principles. Our results show that, due to variation of magnetic field curvature across the polar cap, pair production bursts self-consistently develop an inclination with respect to the local magnetic field that favors the generation of coherent electromagnetic modes with properties consistent with pulsar radio emission. We show that this emission is peaked along the magnetic axis and close to the polar cap edge and may thus offer an explanation for the core and conal components of pulsar radio emission. △ Less

Submitted 26 August, 2021; originally announced August 2021.

Comments: 7 pages, 4 figures, accepted for publication in ApJL

An ALMA view of 11 Dusty Star Forming Galaxies at the peak of Cosmic Star Formation History

Authors: L. Pantoni, M. Massardi, A. Lapi, D. Donevski, Q. D'Amato, M. Giulietti, F. Pozzi, M. Talia, C. Vignali, A. Cimatti, L. Silva, A. Bressan, T. Ronconi

Abstract: We present the ALMA view of 11 main-sequence DSFGs, (sub-)millimeter selected in the GOODS-S field, and spectroscopically confirmed to be at the peak of Cosmic SFH (z = 2-3). Our study combines the analysis of galaxy SED with ALMA continuum and CO spectral emission, by using ALMA Science Archive products at the highest spatial resolution currently available for our sample (< 1 arcsec). We include… ▽ More We present the ALMA view of 11 main-sequence DSFGs, (sub-)millimeter selected in the GOODS-S field, and spectroscopically confirmed to be at the peak of Cosmic SFH (z = 2-3). Our study combines the analysis of galaxy SED with ALMA continuum and CO spectral emission, by using ALMA Science Archive products at the highest spatial resolution currently available for our sample (< 1 arcsec). We include galaxy multi-band images and photometry (in the optical, radio and X-rays) to investigate the interlink between dusty, gaseous and stellar components and the eventual presence of AGN. We use multi-band sizes and morphologies to gain an insight on the processes that lead galaxy evolution, e.g. gas condensation, star formation, AGN feedback. The 11 DSFGs are very compact in the (sub-)millimeter (median r(ALMA) = 1.15 kpc), while the optical emission extends tolarger radii (median r(H)/r(ALMA) = 2.05). CO lines reveal the presence of a rotating disc of molecular gas, but we can not exclude either the presence of interactions and/or molecular outflows. Images at higher (spectral and spatial) resolution are needed to disentangle from the possible scenarios. Most of the galaxies are caught in the compaction phase, when gas cools and falls into galaxy centre, fuelling the dusty burst of star formation and the growing nucleus. We expect these DSFGs to be the high-zstar-forming counterparts of massive quiescent galaxies. Some features of CO emission in three galaxies are suggestive of forthcoming/ongoing AGN feedback, that is thought to trigger the morphological transition from star-forming disks to ETGs. △ Less

Submitted 12 August, 2021; originally announced August 2021.

Comments: 18 pages, 6 figures, 7 tables