Showing 1–50 of 218 results for author: Gomez, F

First resolved stellar halo kinematics of a MW-mass galaxy outside the Local Group: A flat counter-rotating halo in NGC 4945

Authors: Camila Beltrand, Antonela Monachesi, Richard D'Souza, Eric F. Bell, Roelof S. de Jong, Facundo A. Gomez, Jeremy Bailin, In Sung Jang, Adam Smercina

Abstract: Stellar halos of galaxies, primarily formed through the accretion of smaller objects, are important to understand the hierarchical mass assembly of galaxies. However, the inner regions of stellar halos in disk galaxies are predicted to have an in-situ component that is expected to be prominent along the major axis. Kinematic information is crucial to disentangle the contribution of the in-situ com… ▽ More Stellar halos of galaxies, primarily formed through the accretion of smaller objects, are important to understand the hierarchical mass assembly of galaxies. However, the inner regions of stellar halos in disk galaxies are predicted to have an in-situ component that is expected to be prominent along the major axis. Kinematic information is crucial to disentangle the contribution of the in-situ component from the accreted stellar halos. The low surface brightness of stellar halos makes it inaccessible with traditional integrated light spectroscopy. In this work, using a novel technique, we study the kinematics of the stellar halo of the edge-on galaxy NGC 4945. We couple new deep Multi Unit Spectroscopic Explorer spectroscopic observations with existing Hubble Space Telescope imaging data to spectroscopically measure the line-of-sight (LOS) heliocentric velocity and velocity dispersion in two fields at a galactocentric distance of 12.2 kpc (outer disk field) and 34.6 kpc (stellar halo field) along NGC 4945 major axis, by stacking individual spectra of red giant branch and asymptotic giant branch stars. We obtain a LOS velocity and dispersion of 673+/-11 km/s and 73+/-14 km/s, respectively, for the outer disk field. This is consistent with the mean HI velocity of the disk at that distance. For the halo field we obtain a LOS velocity and dispersion of 519+/-12 km/s and 42+/-22 km/s. The halo fields' velocity measurement is within ~40 km/s from the systemic LOS velocity of NGC 4945, which is 563 km/s, suggesting that its stellar halo at 34.6 kpc along the major axis is counter-rotating and is of likely accretion origin. This provides the first ever kinematic measurement of the stellar halo of a Milky Way-mass galaxy outside the Local Group from its resolved stellar population, and establishes a powerful technique for measuring the velocity field of the stellar halos of nearby galaxies. △ Less

Submitted 25 June, 2024; originally announced June 2024.

Comments: Accepted to A&A. 23 pages, 12 figures, 2 appendices

Galaxy evolution in compact groups II. Witnessing the influence of major structures in their evolution

Authors: Gissel P. Montaguth, Antonela Monachesi, Sergio Torres-Flores, Facundo A. Gómez, Ciria Lima-Dias, Arianna Cortesi, Claudia Mendes de Oliveira, Eduardo Telles, Swayamtrupta Panda, Marco Grossi, Paulo A. A. Lopes, Ana Laura O'Mill, Jose A. Hernandez-Jimenez, D. E. Olave-Rojas, Ricardo Demarco, Antonio Kanaan, Tiago Ribeiro, William Schoenell

Abstract: Compact groups (CGs) of galaxies are extreme environments for morphological transformations and the cessation of star formation. Our objective is to understand the dynamics of CGs and how their surrounding environment impacts galaxy properties. We selected a sample of 340 CGs in the Stripe 82 region, totaling 1083 galaxies, and a control sample of 2281 field galaxies. We find that at least 27\% of… ▽ More Compact groups (CGs) of galaxies are extreme environments for morphological transformations and the cessation of star formation. Our objective is to understand the dynamics of CGs and how their surrounding environment impacts galaxy properties. We selected a sample of 340 CGs in the Stripe 82 region, totaling 1083 galaxies, and a control sample of 2281 field galaxies. We find that at least 27\% of our sample of CGs are part of major structures, i.e. non-isolated CGs. We find a bimodality in the effective radius ($R_e$)-Sérsic index ($n$) plane for all transition galaxies (those with $(u-r) > 2.3$ and $n<2.5$) in CGs. Additionally, transition galaxies in isolated CGs populate more densely the $R_e-n$ plane for $n < 1.75$. In contrast, transition galaxies in non-isolated CGs have smoothly increasing $n$ values, suggesting these galaxies have already suffered morphological transformation, and primarily contribute to the distribution of more compact galaxies in the $R_e-n$ plane for all transition galaxies in CGs. We also find significant differences in the specific star-formation rate (sSFR) distribution between the late-type galaxies (LTGs) ($(u-r)<2.3$ and $n< 2.5$) in non-isolated CGs and the same type of galaxies in the control sample, suggesting that the evolution of LTGs differs in non-isolated CGs. Early-type galaxies ($(u-r)>2.3$ and $n>2.5$) and transition galaxies in non-isolated CGs have lower sSFR values and a higher fraction of quenched galaxies, compared to those in isolated CGs. Based on our results, we propose an evolutionary scenario where the major structures in which the CGs are embedded accelerate the morphological transformations of their members. Our findings highlight the importance of considering the larger structures in which CGs may be located, when analysing the properties of their galaxy, as this can significantly affect the evolution of CGs and their galaxies. △ Less

Submitted 20 June, 2024; originally announced June 2024.

Comments: Submitted to A&A, 13 pages, 9 figures

Bar formation and evolution in the cosmological context: Inputs from the Auriga simulations

Authors: Francesca Fragkoudi, Robert Grand, Rüdiger Pakmor, Facundo Gómez, Federico Marinacci, Volker Springel

Abstract: Galactic bars drive the internal evolution of spiral galaxies, while their formation is tightly coupled to the properties of their host galaxy and dark matter halo. To explore what drives bar formation in the cosmological context and how these structures evolve throughout cosmic history, we use the Auriga suite of magneto-hydrodynamical cosmological zoom-in simulations. We find that bars are robus… ▽ More Galactic bars drive the internal evolution of spiral galaxies, while their formation is tightly coupled to the properties of their host galaxy and dark matter halo. To explore what drives bar formation in the cosmological context and how these structures evolve throughout cosmic history, we use the Auriga suite of magneto-hydrodynamical cosmological zoom-in simulations. We find that bars are robust and long-lived structures, and we recover a decreasing bar fraction with increasing redshift which plateaus around $\sim20\%$ at $z\sim3$. We find that bars which form at low and intermediate redshifts grow longer with time, while bars that form at high redshifts are born `saturated' in length, likely due to their merger-induced formation pathway. This leads to a larger bar-to-disc size ratio at high redshifts as compared to the local Universe. We subsequently examine the multi-dimensional parameter space thought to drive bar formation. We find that barred galaxies tend to have lower Toomre $Q$ values at the time of their formation, while we do not find a difference in the gas fraction of barred and unbarred populations when controlling for stellar mass. Barred galaxies tend to be more baryon-dominated at all redshifts, assembling their stellar mass earlier, while galaxies that are baryon-dominated but that do not host a bar, have a higher ex-situ bulge fraction. We explore the implications of the baryon-dominance of barred galaxies on the Tully-Fisher relation, finding an offset from the unbarred relation; confirming this in observations would serve as additional evidence for dark matter, as this behaviour is not readily explained in modified gravity scenarios. △ Less

Submitted 19 June, 2024; v1 submitted 12 June, 2024; originally announced June 2024.

Comments: 21 pages including Appendices, 17 figures, submitted to MNRAS

Cosmological gas accretion history onto the stellar discs of Milky Way-like galaxies in the Auriga simulations -- (II) The inside-out growth of discs

Authors: Federico G. Iza, Sebastián E. Nuza, Cecilia Scannapieco, Robert J. J. Grand, Facundo A. Gómez, Volker Springel, Rüdiger Pakmor, Federico Marinacci, Francesca Fragkoudi

Abstract: We investigate the growth of stellar discs in Milky Way-mass galaxies using the magnetohydrodynamical simulations of the Auriga Project in a full cosmological context. We focus on the gas accretion process along the discs, calculating the net, infall and outflow rates as a function of galactocentric distance, and investigate the relation between them and the star formation activity. The stellar di… ▽ More We investigate the growth of stellar discs in Milky Way-mass galaxies using the magnetohydrodynamical simulations of the Auriga Project in a full cosmological context. We focus on the gas accretion process along the discs, calculating the net, infall and outflow rates as a function of galactocentric distance, and investigate the relation between them and the star formation activity. The stellar distributions of around 70% of the simulated galaxies exhibit an ``inside-out'' pattern, with older (younger) stellar populations preferentially located in the inner (outer) disc regions. In all cases, we find a very tight correlation between the infall, outflow and net accretion rates, as well as between these three quantities and the star formation rate. This is because the amount of gas which is ultimately available for star formation in each radial ring depends not only on the infall rates, but also on the amount of gas leaving the disc in outflows, which directly relates to the local star formation level. Therefore, any of these rates can be used to identify galaxies with inside-out growth. For these galaxies, the correlation between the dominant times of accretion/star formation and disc radius is well fitted by a linear function. We also find that, when averaged over galaxies with formation histories similar to the Milky Way, the simulated accretion rates show a similar evolution (both temporally- and radially-integrated) to the usual accretion prescriptions used in chemical evolution models, although some major differences arise at early times and in the inner disc regions. △ Less

Submitted 10 April, 2024; originally announced April 2024.

Comments: 19 pages, 13 figures

Cosmological evolution of metallicity correlation functions from the Auriga simulations

Authors: Zefeng Li, Robert J. J. Grand, Emily Wisnioski, J. Trevor Mendel, Mark R. Krumholz, Yuan-Sen Ting, Ruediger Pakmor, Facundo A. Gómez, Federico Marinacci, Ioana Ciucă

Abstract: We study the cosmological evolution of the two-point correlation functions of galactic gas-phase metal distributions using the 28 simulated galaxies from the Auriga Project. Using mock observations of the $z = 0$ snapshots to mimic our past work, we show that the correlation functions of the simulated mock observations are well matched to the correlation functions measured from local galaxy survey… ▽ More We study the cosmological evolution of the two-point correlation functions of galactic gas-phase metal distributions using the 28 simulated galaxies from the Auriga Project. Using mock observations of the $z = 0$ snapshots to mimic our past work, we show that the correlation functions of the simulated mock observations are well matched to the correlation functions measured from local galaxy surveys. This comparison suggests that the simulations capture the processes important for determining metal correlation lengths, the key parameter in metallicity correlation functions. We investigate the evolution of metallicity correlations over cosmic time using the true simulation data, showing that individual galaxies undergo no significant systematic evolution in their metal correlation functions from $z\sim 3$ to today. In addition, the fluctuations in metal correlation length are correlated with but lag ahead fluctuations in star formation rate. This suggests that re-arrangement of metals within galaxies occurs at a higher cadence than star formation activity, and is more sensitive to the changes of environment, such as galaxy mergers, gas inflows / outflows, and fly-bys. △ Less

Submitted 13 February, 2024; originally announced February 2024.

Comments: 12 pages, 10 figures, 1 table, accepted for publication in MNRAS

A HITRAN-formatted UV line list of S$_2$ containing transitions involving $X\,^{3}Σ^{-}_{g}$, $B\,^{3}Σ^{-}_{u}$, and $B''\,^{3}Π_{u}$ electronic states

Authors: Frances M. Gomez, Robert J. Hargreaves, Iouli E. Gordon

Abstract: The sulfur dimer (S$_2$) is an important molecular constituent in cometary atmospheres and volcanic plumes on Jupiter's moon Io. It is also expected to play an important role in the photochemistry of exoplanets. The UV spectrum of S$_2$ contains transitions between vibronic levels above and below the dissociation limit, giving rise to a distinctive spectral signature. By using spectroscopic inform… ▽ More The sulfur dimer (S$_2$) is an important molecular constituent in cometary atmospheres and volcanic plumes on Jupiter's moon Io. It is also expected to play an important role in the photochemistry of exoplanets. The UV spectrum of S$_2$ contains transitions between vibronic levels above and below the dissociation limit, giving rise to a distinctive spectral signature. By using spectroscopic information from the literature, and the spectral simulation program PGOPHER, a UV line list of S$_2$ is provided. This line list includes the primary $B\,^{3}Σ^{-}_{u}-X\,^{3}Σ^{-}_{g}$ ($v'$=0-27, $v''$=0-10) electronic transition, where vibrational bands with $v'$$\geq$10 are predissociated. Intensities have been calculated from existing experimental and theoretical oscillator strengths, and semi-empirical strengths for the predissociated bands of S$_2$ have been derived from comparisons with experimental cross-sections. The S$_2$ line list also includes the $B''\,^{3}Π_{u}-X\,^{3}Σ^{-}_{g}$ ($v'$=0-19, $v''$=0-10) vibronic bands due to the strong interaction with the $B$ state. In summary, we present the new HITRAN-formatted S$_2$ line list and its validation against existing laboratory spectra. The extensive line list covers the spectral range 21700$-$41300~cm$^{-1}$ ($\sim$242$-$461~nm) and can be used for modeling both absorption and emission. △ Less

Submitted 19 January, 2024; originally announced January 2024.

Comments: Accepted to the Monthly Notices of the Royal Astronomical Society

Overview and public data release of the augmented Auriga Project: cosmological simulations of dwarf and Milky Way-mass galaxies

Authors: Robert J. J. Grand, Francesca Fragkoudi, Facundo A. Gómez, Adrian Jenkins, Federico Marinacci, Rüdiger Pakmor, Volker Springel

Abstract: We present an extended suite of the Auriga cosmological gravo-magnetohydrodynamical "zoom-in" simulations of 40 Milky Way-mass halos and 26 dwarf galaxy-mass halos run with the moving-mesh code Arepo. Auriga adopts the $Λ$ Cold Dark Matter ($Λ$CDM) cosmogony and includes a comprehensive galaxy formation physics model following the coupled cosmic evolution of dark matter, gas, stars, and supermassi… ▽ More We present an extended suite of the Auriga cosmological gravo-magnetohydrodynamical "zoom-in" simulations of 40 Milky Way-mass halos and 26 dwarf galaxy-mass halos run with the moving-mesh code Arepo. Auriga adopts the $Λ$ Cold Dark Matter ($Λ$CDM) cosmogony and includes a comprehensive galaxy formation physics model following the coupled cosmic evolution of dark matter, gas, stars, and supermassive black holes which has been shown to produce numerically well-converged galaxy properties for Milky Way-mass systems. We describe the first public data release of this augmented suite of Auriga simulations, which includes raw snapshots, group catalogues, merger trees, initial conditions, and supplementary data, as well as public analysis tools with worked examples of how to use the data. To demonstrate the value and robustness of the simulation predictions, we analyse a series of low-redshift global properties that compare well with many observed scaling relations, such as the Tully-Fisher relation, the star-forming main sequence, and HI gas fraction/disc thickness. Finally, we show that star-forming gas discs appear to build rotation and velocity dispersion rapidly for $z\gtrsim 3$ before they "settle" into ever-increasing rotation-dispersion ratios ($V/σ$). This evolution appears to be in rough agreement with some kinematic measurements from H$α$ observations, and demonstrates an application of how to utilise the released data. △ Less

Submitted 3 July, 2024; v1 submitted 16 January, 2024; originally announced January 2024.

Comments: Matches final MNRAS version

Stellar populations and origin of thick disks in AURIGA simulations

Authors: Francesca Pinna, Daniel Walo-Martín, Robert J. J. Grand, Marie Martig, Francesca Fragkoudi, Facundo A. Gómez, Federico Marinacci, Rüdiger Pakmor

Abstract: The origin of thick disks and their evolutionary connection with thin disks are still a matter of debate. We provide new insights into this topic by connecting the stellar populations of thick disks at redshift $z=0$ with their past formation and growth, in 24 Milky Way-mass galaxies from the AURIGA zoom-in cosmological simulations. We projected each galaxy edge on, and decomposed it morphological… ▽ More The origin of thick disks and their evolutionary connection with thin disks are still a matter of debate. We provide new insights into this topic by connecting the stellar populations of thick disks at redshift $z=0$ with their past formation and growth, in 24 Milky Way-mass galaxies from the AURIGA zoom-in cosmological simulations. We projected each galaxy edge on, and decomposed it morphologically into two disk components, in order to define geometrically the thin and the thick disks as usually done in observations. We produced age, metallicity and [Mg/Fe] edge-on maps. We quantified the impact of satellite mergers by mapping the distribution of ex-situ stars. Thick disks are on average $\sim 3$~Gyr older, $\sim 0.25$~dex more metal poor and $\sim 0.06$~dex more [Mg/Fe]-enhanced than thin disks. Their average ages range from $\sim 6$ to $\sim 9$~Gyr, metallicities from $\sim -0.15$ to $\sim 0.1$~dex, and [Mg/Fe] from $\sim 0.12$ to $\sim 0.16$~dex. These properties are the result of an early initial in-situ formation, followed by a later growth driven by the combination of direct accretion of stars, some in-situ star formation fueled by mergers, and dynamical heating of stars. The balance between these processes varies from galaxy to galaxy. Mergers play a key role in the mass assembly of thick disks, contributing an average accreted mass fraction of $\sim 22$\% in the analyzed thick-disk dominated regions. In two galaxies, about half of the geometric thick-disk mass was directly accreted. While primordial thick disks form at high redshift in all galaxies, young metal-rich thin disks, with much lower [Mg/Fe] abundances, start to form later but at different times (higher or lower redshift) depending on the galaxy. We conclude that thick disks result from the interplay of external processes with the internal evolution of the galaxy. △ Less

Submitted 2 April, 2024; v1 submitted 22 November, 2023; originally announced November 2023.

Comments: 57 pages, 43 figures, 3 tables; Published in Astronomy & Astrophysics

Journal ref: A&A, 683 (2024) A236

Galactoseismology in cosmological simulations: Vertical perturbations by dark matter, satellite galaxies and gas

Authors: B. García-Conde, T. Antoja, S. Roca-Fàbrega, F. Gómez, P. Ramos, N. Garavito-Camargo, MA. Gómez-Flechoso

Abstract: Only recently, complex models that include the global dynamics from dwarf satellite galaxies, dark matter halo structure, gas infalls, and stellar disk in a cosmological context became available to study the dynamics of disk galaxies such as the Milky Way (MW). We use a MW model from a high-resolution hydrodynamical cosmological simulation named GARROTXA to establish the relationship between the v… ▽ More Only recently, complex models that include the global dynamics from dwarf satellite galaxies, dark matter halo structure, gas infalls, and stellar disk in a cosmological context became available to study the dynamics of disk galaxies such as the Milky Way (MW). We use a MW model from a high-resolution hydrodynamical cosmological simulation named GARROTXA to establish the relationship between the vertical disturbances seen in its galactic disk and multiple perturbations, from the dark matter halo, satellites and gas. We calculate the bending modes in the galactic disk in the last 6 Gyr of evolution. To quantify the impact of dark matter and gas we compute the vertical acceleration exerted by these components onto the disk and compare them with the bending behavior with Fourier analysis. We find complex bending patterns at different radii and times, such as an inner retrograde mode with high frequency, as well as an outer slower retrograde mode excited at different times. The amplitudes of these bending modes are highest during the early stages of the thin disk formation and reach up to 8.5 km s-1 in the late disk evolution. We find that the infall of satellite galaxies leads to a tilt of the disk, and produces anisotropic gas accretion with subsequent star formation events, and supernovae, creating significant vertical accelerations onto the disk plane. The misalignment between the disk and the inner stellar/dark matter triaxial structure, formed during the ancient assembly of the galaxy, creates a strong vertical acceleration on the stars. We conclude that several agents trigger the bending of the stellar disk and its phase spirals in this simulation, including satellite galaxies, dark sub-halos, misaligned gaseous structures, and the inner dark matter profile, which coexist and influence each other, making it challenging to establish direct causality. △ Less

Submitted 13 November, 2023; originally announced November 2023.

An SiO Toroid and Wide-angle Outflow associated with the Massive Protostar W75N(B)-VLA2

Authors: Jose F. Gomez, Jose M. Torrelles, Josep M. Girart, Gabriele Surcis, Jeong-Sook Kim, Jorge Canto, Guillem Anglada, Salvador Curiel, Wouter H. T. Vlemmings, Carlos Carrasco-Gonzalez, Adriana R. Rodriguez-Kamenetzky, Soon-Wook Kim, Ciriaco Goddi, Huib J. van Langevelde, Alvaro Sanchez-Monge

Abstract: We have carried out ALMA observations of the massive star-forming region W75N(B), which contains the massive protostars VLA1, VLA2, and VLA3. Particularly, VLA2 is an enigmatic protostar associated with a wind-driven H$_2$O maser shell, which has evolved from an almost isotropic outflow to a collimated one in just 20 years. The shell expansion seemed to be halted by an obstacle located to the nort… ▽ More We have carried out ALMA observations of the massive star-forming region W75N(B), which contains the massive protostars VLA1, VLA2, and VLA3. Particularly, VLA2 is an enigmatic protostar associated with a wind-driven H$_2$O maser shell, which has evolved from an almost isotropic outflow to a collimated one in just 20 years. The shell expansion seemed to be halted by an obstacle located to the northeast of VLA2. Here we present our findings from observing the 1.3 mm continuum and H$_2$CO and SiO emission lines. Within a region of $\sim 30"$ ($\sim 39,000$ au) diameter, we have detected 40 compact mm-continuum sources, three of them coinciding with VLA1, VLA2, and VLA3. While the H$_2$CO emission is mainly distributed in a fragmented structure around the three massive protostars, but without any of the main H$_2$CO clumps spatially coinciding with them, the SiO is highly concentrated on VLA2, indicating the presence of very strong shocks generated near this protostar. The SiO emission is clearly resolved into an elongated structure ($\sim 0.6"\times0.3"$; $\sim 780$ au$\times$390 au) perpendicular to the major axis of the wind-driven maser shell. The structure and kinematics of the SiO emission are consistent with a toroid and a wide-angle outflow surrounding a central mass of $\sim 10$ M$_{\odot}$, thus supporting previous theoretical predictions regarding the evolution of the outflow. Additionally, we have identified the expected location and estimated the gas density of the obstacle that is hindering the expansion of the maser shell. △ Less

Submitted 13 October, 2023; v1 submitted 12 October, 2023; originally announced October 2023.

Comments: To be published in The Astrophysical Journal Letters. Sixteen pages, seven figures. Updated metadata

Galaxy evolution in modified gravity simulations: using passive galaxies to constrain gravity with upcoming surveys

Authors: Diego Pallero, Facundo A. Gómez, Nelson D. Padilla, Yara L. Jaffé, Carlton M. Baugh, Baojiu Li, César Hernández-Aguayo, Christian Arnold

Abstract: We present a quantitative analysis of the properties of galaxies and structures evolving in universes dominated by different modified gravitational models, including two variants of the f(R)-gravity (F) and two of the Dvali-Gabdadze-Poratti (N) braneworld model, which respectively feature the chameleon and Vainshtein screening mechanisms. Using the Simulation HYdrodynamics BeyONd Einstein (SHYBONE… ▽ More We present a quantitative analysis of the properties of galaxies and structures evolving in universes dominated by different modified gravitational models, including two variants of the f(R)-gravity (F) and two of the Dvali-Gabdadze-Poratti (N) braneworld model, which respectively feature the chameleon and Vainshtein screening mechanisms. Using the Simulation HYdrodynamics BeyONd Einstein (SHYBONE) cosmological hydrodynamical full-physics simulations suite, we study the departures in the properties of galaxies residing in different environments with respect to the standard model (GR). Using two different criteria to compare, we find that structures formed within modified gravity tend to show a denser gas density profile than their GR counterparts. Within the different modified gravity models, N1 and F5 gravity models show greater departures from the standard model, with gas density profiles $ρ_{\rm IGM} \geq 30\%$ denser in the outskirts for the N1 model, and in the inner parts for the F5 model. Additionally, we find that haloes evolving in MG universes show, in general, larger quenched fractions than GR, reaching up to $20\%$ larger quenching fractions in F5 regardless of the stellar mass of the galaxy. With respect to the other models, F6, N1 and N5 show slightly larger quenched fractions, but no strong differences can be found. These results directly impact the colour distribution of galaxies, making them in MG models redder and older than their GR counterparts. Like GR, once the environment starts to play a role, galaxies rapidly get quenched and the differences between models vanish. △ Less

Submitted 3 October, 2023; originally announced October 2023.

Comments: 17 pages, 12 figures, submitted to MNRAS

FLASHING: Project Overview

Authors: Hiroshi Imai, Yuhki Hamae, Kei Amada, Keisuke Nakashima, Ka-Yiu Shum, Rina Kasai, Jose F. Gomez, Lucero Uscanga, Daniel Tafoya, Gabor Orosz, Ross A. Burns

Abstract: This paper describes the overview of the FLASHING (Finest Legacy Acquisitions of SiO-/ H$_2$O-maser Ignitions by the Nobeyama Generation) project promoted using the 45 m telescope of Nobeyama Radio Observatory, which aims to intensively monitor H$_2$O (22 GHz) and SiO (43 GHz) masers associated with so-called "water fountain" sources. Here we show scientific results on the basis of the data taken… ▽ More This paper describes the overview of the FLASHING (Finest Legacy Acquisitions of SiO-/ H$_2$O-maser Ignitions by the Nobeyama Generation) project promoted using the 45 m telescope of Nobeyama Radio Observatory, which aims to intensively monitor H$_2$O (22 GHz) and SiO (43 GHz) masers associated with so-called "water fountain" sources. Here we show scientific results on the basis of the data taken in for the first five seasons of FLASHING, from 2018 December to 2023 April). We have found the evolution of the H$_2$O maser spectra, such as new spectral components breaking the record of the jet's top speed and/or systematic velocity drifts in the spectrum indicating acceleration or deceleration of the maser gas clumps. For the 43 GHz SiO maser emission, we have found its new detection in a source while its permanent disappearance in other source. Our finding may imply that the jets from these water fountains can be accelerated or decelerated, and show how cicumstellar envelopes have been destroyed. △ Less

Submitted 1 September, 2023; originally announced September 2023.

Comments: 31 pages, 22 figures, to be published in the Publications of the Astronomical Society of Japan

Clash of Titans: The Impact of Cluster Mergers in the Galaxy Cluster Red Sequence

Authors: Franklin Aldás, Alfredo Zenteno, Facundo Gómez, Daniel Hernandez-Lang, Eleazar R. Carrasco, Cristian A. Vega-Martínez, J. L. Nilo Castellón

Abstract: Merging of galaxy clusters are some of the most energetic events in the Universe, and they provide a unique environment to study galaxy evolution. We use a sample of 84 merging and relaxed SPT galaxy clusters candidates, observed with the Dark Energy Camera in the $0.11<z<0.88$ redshift range, to build colour-magnitude diagrams to characterize the impact of cluster mergers on the galaxy population… ▽ More Merging of galaxy clusters are some of the most energetic events in the Universe, and they provide a unique environment to study galaxy evolution. We use a sample of 84 merging and relaxed SPT galaxy clusters candidates, observed with the Dark Energy Camera in the $0.11<z<0.88$ redshift range, to build colour-magnitude diagrams to characterize the impact of cluster mergers on the galaxy population. We divided the sample between relaxed and disturbed, and in two redshifts bin at $z = 0.55$. When comparing the high-z to low-z clusters we find the high-z sample is richer in blue galaxies, independently of the cluster dynamical state. In the high-z bin we find that disturbed clusters exhibit a larger scatter in the Red Sequence, with wider distribution and an excess of bluer galaxies compared to relaxed clusters, while in the low-z bin we find a complete agreement between the relaxed and disturbed clusters. Our results support the scenario in which massive cluster halos at $z<0.55$ galaxies are quenched as satellites of another structure, i.e. outside the cluster, while at $z \geq 0.55$ the quenching is dominated by in-situ processes. △ Less

Submitted 21 July, 2023; originally announced July 2023.

Comments: 9 pages, 5 figures, paper accepted in MNRAS

ACM Class: J.2.3

Galaxy evolution in compact groups I: Revealing a transitional galaxy population through a multiwavelength approach

Authors: Gissel P. Montaguth, Sergio Torres-Flores, Antonela Monachesi, Facundo A. Gómez, Ciria Lima-Dias, Arianna Cortesi, Claudia Mendes de Oliveira, Eduardo Telles, Swayamtrupta Panda, Marco Grossi, Paulo A. A. Lopes, Jose A. Hernandez-Jimenez, Antonio Kanaan, Tiago Ribeiro, William Schoenell

Abstract: Compact groups of galaxies (CGs) show members with morphological disturbances, mainly products of galaxy-galaxy interactions, thus making them ideal systems to study galaxy evolution, in high-density environment. To understand how this environment affects the properties of galaxies, we select a sample of 340 CGs in the Stripe 82 region, for a total of 1083 galaxies, and a sample of 2281 field gala… ▽ More Compact groups of galaxies (CGs) show members with morphological disturbances, mainly products of galaxy-galaxy interactions, thus making them ideal systems to study galaxy evolution, in high-density environment. To understand how this environment affects the properties of galaxies, we select a sample of 340 CGs in the Stripe 82 region, for a total of 1083 galaxies, and a sample of 2281 field galaxies as a control sample. By performing a multi-wavelength morphological fitting process using S-PLUS data, we divide our sample into early-type (ETG), late-type (LTG), and transition galaxies using the r-band Sérsic index and the colour (u-r). We find a bimodal distribution in the plane of the effective radius-Sérsic index, where a secondary "peculiar" galaxy population of smaller and more compact galaxies is found in CGs, which is not observed in the control sample. This indicates that galaxies are undergoing a morphological transformation in CGs. In addition, we find significant statistical differences in the distribution of specific Star Formation Rate (sSFR) when we compare both environments for LTGs and ETGs. We also find a higher fraction of quenched galaxies and a lower median sSFR in CGs than in the control sample, suggesting the existence of environmental effects favoring the cessation of star formation, regardless of galaxy type. Our results support the notion that CGs promote morphological and physical transformations, highlighting their potential as ideal systems for galaxy pre-processing. △ Less

Submitted 21 July, 2023; originally announced July 2023.

Comments: 17 pages, 15 figures. Accepted for publication in MNRAS

A reproduction of the Milky Way's Faraday rotation measure map in galaxy simulations from global to local scales

Authors: Stefan Reissl, Ralf S. Klessen, Eric W. Pellegrini, Daniel Rahner, Rüdiger Pakmor, Robert Grand, Facundo Gomez, Federico Marinacci, Volker Springel

Abstract: Magnetic fields are of critical importance for our understanding of the origin and long-term evolution of the Milky Way. This is due to their decisive role in the dynamical evolution of the interstellar medium (ISM) and their influence on the star-formation process. Faraday rotation measures (RM) along many different sightlines across the Galaxy are a primary means to infer the magnetic field topo… ▽ More Magnetic fields are of critical importance for our understanding of the origin and long-term evolution of the Milky Way. This is due to their decisive role in the dynamical evolution of the interstellar medium (ISM) and their influence on the star-formation process. Faraday rotation measures (RM) along many different sightlines across the Galaxy are a primary means to infer the magnetic field topology and strength from observations. However, the interpretation of the data has been hampered by the failure of previous attempts to explain the observations in theoretical models and to synthesize a realistic multi-scale all-sky RM map. We here utilize a cosmological magnetohydrodynamic (MHD) simulation of the formation of the Milky Way, augment it with a novel star cluster population synthesis model for a more realistic structure of the local interstellar medium, and perform detailed polarized radiative transfer calculations on the resulting model. This yields a faithful first principles prediction of the Faraday sky as observed on Earth. The results reproduce the observations of the Galaxy not only on global scales, but also on local scales of individual star-forming clouds. They also imply that the Local Bubble containing our Sun dominates the RM signal over large regions of the sky. Modern cosmological MHD simulations of the Milky Way's formation, combined with a simple and plausible model for the fraction of free electrons in the ISM, explain the RM observations remarkably well, thus indicating the emergence of a firm theoretical understanding of the genesis of magnetic fields in our Universe across cosmic time. △ Less

Submitted 11 July, 2023; originally announced July 2023.

Comments: 14 pages, 3 figures

Structure, Kinematics, and Observability of the Large Magellanic Cloud's Dynamical Friction Wake in Cold vs. Fuzzy Dark Matter

Authors: Hayden R. Foote, Gurtina Besla, Philip Mocz, Nicolás Garavito-Camargo, Lachlan Lancaster, Martin Sparre, Emily C. Cunningham, Mark Vogelsberger, Facundo A. Gómez, Chervin F. P. Laporte

Abstract: The Large Magellanic Cloud (LMC) will induce a dynamical friction (DF) wake on infall to the Milky Way (MW). The MW's stellar halo will respond to the gravity of the LMC and the dark matter (DM) wake, forming a stellar counterpart to the DM wake. This provides a novel opportunity to constrain the properties of the DM particle. We present a suite of high-resolution, windtunnel-style simulations of… ▽ More The Large Magellanic Cloud (LMC) will induce a dynamical friction (DF) wake on infall to the Milky Way (MW). The MW's stellar halo will respond to the gravity of the LMC and the dark matter (DM) wake, forming a stellar counterpart to the DM wake. This provides a novel opportunity to constrain the properties of the DM particle. We present a suite of high-resolution, windtunnel-style simulations of the LMC's DF wake that compare the structure, kinematics, and stellar tracer response of the DM wake in cold DM (CDM), with and without self-gravity, vs. fuzzy DM (FDM) with $m_a = 10^{-23}$ eV. We conclude that the self-gravity of the DM wake cannot be ignored. Its inclusion raises the wake's density by $\sim 10\%$, and holds the wake together over larger distances ($\sim$ 50 kpc) than if self-gravity is ignored. The DM wake's mass is comparable to the LMC's infall mass, meaning the DM wake is a significant perturber to the dynamics of MW halo tracers. An FDM wake is more granular in structure and is $\sim 20\%$ dynamically colder than a CDM wake, but with comparable density. The granularity of an FDM wake increases the stars' kinematic response at the percent level compared to CDM, providing a possible avenue of distinguishing a CDM vs. FDM wake. This underscores the need for kinematic measurements of stars in the stellar halo at distances of 70-100 kpc. △ Less

Submitted 8 September, 2023; v1 submitted 30 June, 2023; originally announced July 2023.

Comments: 38 pages, 31 figures. Updated to published version

Journal ref: ApJ 954 163 (2023)

Lopsidedness as a tracer of early galactic assembly history

Authors: Arianna Dolfi, Facundo A. Gomez, Antonela Monachesi, Silvio Varela-Lavin, Patricia B. Tissera, Cristobal Sifon, Gaspar Galaz

Abstract: Large-scale asymmetries (i.e. lopsidedness) are a common feature in the stellar density distribution of nearby disk galaxies both in low- and high-density environments. In this work, we characterize the present-day lopsidedness in a sample of 1435 disk-like galaxies selected from the TNG50 simulation. We find that the percentage of lopsided galaxies (10%-30%) is in good agreement with observations… ▽ More Large-scale asymmetries (i.e. lopsidedness) are a common feature in the stellar density distribution of nearby disk galaxies both in low- and high-density environments. In this work, we characterize the present-day lopsidedness in a sample of 1435 disk-like galaxies selected from the TNG50 simulation. We find that the percentage of lopsided galaxies (10%-30%) is in good agreement with observations if we use similar radial ranges to the observations. However, the percentage (58%) significantly increases if we extend our measurement to larger radii. We find a mild or lack of correlation between lopsidedness amplitude and environment at z=0 and a strong correlation between lopsidedness and galaxy morphology regardless of the environment. Present-day galaxies with more extended disks, flatter inner galactic regions and lower central stellar mass density (i.e. late-type disk galaxies) are typically more lopsided than galaxies with smaller disks, rounder inner galactic regions and higher central stellar mass density (i.e. early-type disk galaxies). Interestingly, we find that lopsided galaxies have, on average, a very distinct star formation history within the last 10 Gyr, with respect to their symmetric counterparts. Symmetric galaxies have typically assembled at early times (~8-6 Gyr ago) with relatively short and intense bursts of central star formation, while lopsided galaxies have assembled on longer timescales and with milder initial bursts of star formation, continuing building up their mass until z=0. Overall, these results indicate that lopsidedness in present-day disk galaxies is connected to the specific evolutionary histories of the galaxies that shaped their distinct internal properties. △ Less

Submitted 7 June, 2023; originally announced June 2023.

Comments: Submitted to MNRAS, 23 pages (including the Appendix), 23 figures

Exploring the diversity and similarity of radially anisotropic Milky Way-like stellar haloes: implications for disrupted dwarf galaxy searches

Authors: Matthew D. A. Orkney, Chervin F. P. Laporte, Robert J. J. Grand, Facundo A. Gómez, Freeke van de Voort, Azadeh Fattahi, Federico Marinacci, Rüdiger Pakmor, Francesca Fragkoudi, Volker Springel

Abstract: We investigate the properties of mergers comparable to the Gaia-Sausage-Enceladus (GSE) using cosmological hydrodynamical simulations of Milky Way-like galaxies. The merger progenitors span an order of magnitude in their peak stellar mass ($3\times10^8<M_{\star}/\rm{M}_{\odot}<4\times10^9$) and include both rotation and pressure-supported galaxies ($0.10<D/T<0.77$). In a minority of cases, the GSE… ▽ More We investigate the properties of mergers comparable to the Gaia-Sausage-Enceladus (GSE) using cosmological hydrodynamical simulations of Milky Way-like galaxies. The merger progenitors span an order of magnitude in their peak stellar mass ($3\times10^8<M_{\star}/\rm{M}_{\odot}<4\times10^9$) and include both rotation and pressure-supported galaxies ($0.10<D/T<0.77$). In a minority of cases, the GSE-like debris is comprised of stars from more than one merger progenitor. However, there is a close similarity in their chemodynamical properties and the triaxial shapes of their debris, and so it is not always possible to distinguish them. The merger progenitors host a variety of luminous satellites ($0-8$ with $M_{\star}>10^6\,\rm{M}_{\odot}$), but most of these do not follow the merger to low orbital energies. Between $0-1$ of these satellites may survive to $z=0$, but with no clear signatures of their past association. We show that the fraction of stars originating from GSE-like mergers is reduced for lower metallicities (reaching a minimum around $\text{[Fe/H]} = -2$), and also within $5\,$kpc of the galactic centre. Whilst these central regions are dominated by in-situ stars, the ex-situ fraction trends towards a 100 per cent asymptote when considering the most metal-poor stars ($\text{[Fe/H]}\ll-2.5$). Considering this, its near proximity, and its small volume on the sky, the Galactic centre lends itself as a prime environment in the search for the stars from the earliest galaxies, whilst avoiding contamination from GSE stars. △ Less

Submitted 7 September, 2023; v1 submitted 3 March, 2023; originally announced March 2023.

Comments: 23 pages, 18 figures, accepted by MNRAS

Evolution of the outflow in the water fountain source IRAS 18043$-$2116

Authors: L. Uscanga, H. Imai, J. F. Gómez, D. Tafoya, G. Orosz, T. P. McCarthy, Y. Hamae, K. Amada

Abstract: We present the spectral and spatial evolution of H$_2$O masers associated with the water fountain source IRAS 18043$-$2116 found in the observations with the Nobeyama 45 m telescope and the Australia Telescope Compact Array. We have found new highest velocity components of the H$_2$O masers (at the red-shifted side $V_{\rm LSR}\simeq376~\mathrm{km~s}^{-1}$ and at the blue-shifted side… ▽ More We present the spectral and spatial evolution of H$_2$O masers associated with the water fountain source IRAS 18043$-$2116 found in the observations with the Nobeyama 45 m telescope and the Australia Telescope Compact Array. We have found new highest velocity components of the H$_2$O masers (at the red-shifted side $V_{\rm LSR}\simeq376~\mathrm{km~s}^{-1}$ and at the blue-shifted side $V_{\rm LSR}\simeq$ $-165~\mathrm{km~s}^{-1}$), and the resulting velocity spread of $\simeq 540~\mathrm{km~s}^{-1}$ breaks the speed record of fast jets/outflows in this type of sources. The locations of those components have offsets from the axis joining the two major maser clusters, indicating a large opening angle of the outflow ($\sim60^{\circ}$). The evolution of the maser cluster separation of $\sim$2.9 mas yr$^{-1}$ and the compact ($\sim0.''2$) CO emission source mapped with the Atacama Large Millimeter-submillimeter Array suggest a very short ($\sim$30 yr) timescale of the outflow. We also confirmed the increase in the flux density of the 22 GHz continuum source. The properties of the jet and the continuum sources and their possible evolution in the transition to the planetary nebula phase are further discussed. △ Less

Submitted 2 March, 2023; originally announced March 2023.

Comments: 14 pages, 7 figures. Accepted by ApJ

Monitoring of the polarized $H_{2}O$ maser emission around the massive protostars W75N(B)-VLA1 and W75N(B)-VLA2

Authors: G. Surcis, W. H. T. Vlemmings, C. Goddi, J. M. Torrelles, J. F. Gómez, A. Rodríguez-Kamenetzky, C. Carrasco-González, S. Curiel, S. -W. Kim, J. -S. Kim, H. J. van Langevelde

Abstract: Several radio sources have been detected in the HMSFR W75N(B), among them the massive YSOs VLA1 and VLA2 are of great interest. These are thought to be in different evolutionary stages. In particular, VLA1 is at the early stage of the photoionization and it is driving a thermal radio jet, while VLA2 is a thermal, collimated ionized wind surrounded by a dusty disk or envelope. In both sources 22 GH… ▽ More Several radio sources have been detected in the HMSFR W75N(B), among them the massive YSOs VLA1 and VLA2 are of great interest. These are thought to be in different evolutionary stages. In particular, VLA1 is at the early stage of the photoionization and it is driving a thermal radio jet, while VLA2 is a thermal, collimated ionized wind surrounded by a dusty disk or envelope. In both sources 22 GHz water masers have been detected in the past. Those around VLA1 show a persistent distribution along the radio jet and those around VLA2 have instead traced the evolution from a non-collimated to a collimated outflow over a period of 20 years. By monitoring the polarized emission of the water masers around both VLA1 and VLA2 over a period of 6 years, we aim to determine whether the maser distributions show any variation over time and whether the magnetic field behaves accordingly. The EVN was used in full polarization and phase-reference mode to measure the absolute positions of the masers and to determine both the orientation and the strength of the magnetic field. We observed four epochs separated by two years from 2014 to 2020. We detected polarized emission from the water masers around both the YSOs in all the epochs. We find that the masers around VLA1 are tracing a nondissociative shock originating from the expansion of the thermal radio jet, while the masers around VLA2 are tracing an asymmetric expansion of the gas that is halted in the northeast where the gas likely encounters a very dense medium. We also found that the magnetic field inferred from the water masers in each epoch can be considered as a portion of a quasi-static magnetic field estimated in that location rather than in that time. This allowed us to study locally the morphology of the magnetic field around both YSOs in a larger area by considering the vectors estimated in all the epochs as a whole. △ Less

Submitted 28 February, 2023; originally announced March 2023.

Comments: 32 pages, 15 figures, 17 tables, accepted by Astronomy & Astrophysics

Journal ref: A&A 673, A10 (2023)

The impact of the Large Magellanic Cloud on dark matter direct detection signals

Authors: Adam Smith-Orlik, Nima Ronaghi, Nassim Bozorgnia, Marius Cautun, Azadeh Fattahi, Gurtina Besla, Carlos S. Frenk, Nicolás Garavito-Camargo, Facundo A. Gómez, Robert J. J. Grand, Federico Marinacci, Annika H. G. Peter

Abstract: We study the effect of the Large Magellanic Cloud (LMC) on the dark matter (DM) distribution in the Solar neighborhood, utilizing the Auriga magneto-hydrodynamical simulations of Milky Way (MW) analogues that have an LMC-like system. We extract the local DM velocity distribution at different times during the orbit of the LMC around the MW in the simulations. As found in previous idealized simulati… ▽ More We study the effect of the Large Magellanic Cloud (LMC) on the dark matter (DM) distribution in the Solar neighborhood, utilizing the Auriga magneto-hydrodynamical simulations of Milky Way (MW) analogues that have an LMC-like system. We extract the local DM velocity distribution at different times during the orbit of the LMC around the MW in the simulations. As found in previous idealized simulations of the MW-LMC system, we find that the DM particles in the Solar neighborhood originating from the LMC analogue dominate the high speed tail of the local DM speed distribution. Furthermore, the native DM particles of the MW in the Solar region are boosted to higher speeds as a result of a response to the LMC's motion. We simulate the signals expected in near future xenon, germanium, and silicon direct detection experiments, considering DM interactions with target nuclei or electrons. We find that the presence of the LMC causes a considerable shift in the expected direct detection exclusion limits towards smaller cross sections and DM masses, with the effect being more prominent for low mass DM. Hence, our study shows, for the first time, that the LMC's influence on the local DM distribution is significant even in fully cosmological MW analogues. △ Less

Submitted 25 October, 2023; v1 submitted 8 February, 2023; originally announced February 2023.

Comments: 35 pages, 14 figures, 3 tables, JCAP accepted version

Lopsided Galaxies in a cosmological context: a new galaxy-halo connection

Authors: Silvio Varela-Lavin, Facundo A. Gómez, Patricia B. Tissera, Gurtina Besla, Nicolás Garavito-Camargo, Federico Marinacci, Chervin F. P. Laporte

Abstract: Disc galaxies commonly show asymmetric features in their morphology, such as warps and lopsidedness. These features can provide key information regarding the recent evolution of a given disc galaxy. In the nearby Universe, up to $\sim30$ percent of late-type galaxies display a global non-axisymmetric lopsided mass distribution. However, the origin of this perturbation is not well understood. In th… ▽ More Disc galaxies commonly show asymmetric features in their morphology, such as warps and lopsidedness. These features can provide key information regarding the recent evolution of a given disc galaxy. In the nearby Universe, up to $\sim30$ percent of late-type galaxies display a global non-axisymmetric lopsided mass distribution. However, the origin of this perturbation is not well understood. In this work, we study the origin of lopsided perturbations in simulated disc galaxies extracted from the TNG50 simulation of the IllustrisTNG project. We statistically explore different excitation mechanisms for this perturbation, such as direct satellite tidal interactions and distortions of the underlying dark matter distributions. We also characterize the main physical conditions that lead to lopsided perturbations. 50 percent of our sample galaxy have lopsided modes $m=1$ greater than $\sim 0.12$. We find a strong correlation between internal galaxy properties, such as central stellar surface density and disc radial extension with the strength of lopsided modes. The majority of lopsided galaxies have lower central surface densities and more extended discs than symmetric galaxies. As a result, such lopsided galaxies are less self-gravitationally cohesive, and their outer disc region is more susceptible to different types of external perturbations. However, we do not find strong evidence that tidal interactions with satellite galaxies are the main driving agent of lopsided modes. Lopsided galaxies tend to live in asymmetric dark matter halos with high spin, indicating strong galaxy-halo connections in late-type lopsided galaxies. △ Less

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

Comments: 16 pages, 16 figures, accepted for publication in MNRAS

An ever-present $Gaia$ snail shell triggered by a dark matter wake

Authors: Robert J. J. Grand, Rüdiger Pakmor, Francesca Fragkoudi, Facundo A. Gómez, Wilma Trick, Christine M. Simpson, Freeke van de Voort, Rebekka Bieri

Abstract: We utilize a novel numerical technique to model star formation in cosmological simulations of galaxy formation - called Superstars - to simulate a Milky Way-like galaxy with $\gtrsim10^8$ star particles to study the formation and evolution of out-of-equilibrium stellar disc structures in a full cosmological setting. In the plane defined by the coordinate and velocity perpendicular to the mid-plane… ▽ More We utilize a novel numerical technique to model star formation in cosmological simulations of galaxy formation - called Superstars - to simulate a Milky Way-like galaxy with $\gtrsim10^8$ star particles to study the formation and evolution of out-of-equilibrium stellar disc structures in a full cosmological setting. In the plane defined by the coordinate and velocity perpendicular to the mid-plane (vertical phase space, $\{Z,V_Z\}$), stars in Solar-like volumes at late times exhibit clear spirals qualitatively similar in shape and amplitude to the $Gaia$ ``Snail shell'' phase spiral. We show that the phase spiral forms at a look back time of $\sim 6$ Gyr during the pericentric passage of a $\sim10^{10}$ $\rm M_{\odot}$ satellite on a polar orbit. This satellite stimulates the formation of a resonant wake in the dark matter halo while losing mass at a rate of $\sim0.5$-$1$ dex per orbit loop. The peak magnitude of the wake-induced gravitational torque at the Solar radius is $\sim 8$ times that from the satellite, and triggers the formation of a disc warp that wraps up into a vertical phase spiral over time. As the wake decays, the phase spiral propagates several Gigayears to present-day and can be described as ``ever-present'' once stable disc evolution is established. These results suggest an alternative scenario to explain the $Gaia$ phase spiral which does not rely on a perturbation from bar buckling or a recent direct hit from a satellite. △ Less

Submitted 28 June, 2023; v1 submitted 15 November, 2022; originally announced November 2022.

Comments: accepted for pub in mnras. 11 + 5 pages. For a high-cadence animated version of figure 1, see https://wwwmpa.mpa-garching.mpg.de/auriga/movies/multi_halo_6_sf64.mp4

Cosmological gas accretion history onto the stellar discs of Milky Way-like galaxies in the Auriga simulations -- (I) Temporal dependency

Authors: Federico G. Iza, Cecilia Scannapieco, Sebastián E. Nuza, Robert J. J. Grand, Facundo A. Gómez, Volker Springel, Rüdiger Pakmor, Federico Marinacci

Abstract: We use the 30 simulations of the Auriga Project to estimate the temporal dependency of the inflow, outflow and net accretion rates onto the discs of Milky Way-like galaxies. The net accretion rates are found to be similar for all galaxies at early times, increasing rapidly up to $\sim 10~\mathrm{M}_\odot \, \mathrm{yr}^{-1}$. After $\sim 6~\mathrm{Gyr}$ of evolution, however, the net accretion rat… ▽ More We use the 30 simulations of the Auriga Project to estimate the temporal dependency of the inflow, outflow and net accretion rates onto the discs of Milky Way-like galaxies. The net accretion rates are found to be similar for all galaxies at early times, increasing rapidly up to $\sim 10~\mathrm{M}_\odot \, \mathrm{yr}^{-1}$. After $\sim 6~\mathrm{Gyr}$ of evolution, however, the net accretion rates are diverse: in most galaxies, these exhibit an exponential-like decay, but some systems instead present increasing or approximately constant levels up to the present time. An exponential fit to the net accretion rates averaged over the MW analogues yields typical decay time-scale of $7.2~\mathrm{Gyr}$. The analysis of the time-evolution of the inflow and outflow rates, and their relation to the star formation rate (SFR) in the discs, confirms the close connection between these quantities. First, the inflow$/$outflow ratio stays approximately constant, with typical values of $\dot{M}_\mathrm{out}/ \dot{M}_\mathrm{in} \sim 0.75$, indicating that the gas mass involved in outflows is of the order of 25% lower compared to that involved in inflows. A similar behaviour is found for the SFR$/$inflow rate ratio, with typical values between 0.1 and 0.3, and for the outflow rate$/$SFR which varies in the range $3.5$--$5.5$. Our results show that continuous inflow is key to the SFR levels in disc galaxies, and that the star formation activity and the subsequent feedback in the discs is able to produce mass-loaded galaxy winds in the disc-halo interface. △ Less

Submitted 6 October, 2022; originally announced October 2022.

Comments: 21 pages, 15 figures

Journal ref: Mon. Not. Roy. Astron. Soc., 517, 832-852, 2022

Colour and infall time distributions of satellite galaxies in simulated Milky-Way analogs

Authors: Yue Pan, Christine M. Simpson, Andrey Kravtsov, Facundo A. Gómez, Robert J. J. Grand, Federico Marinacci, Rüdiger Pakmor, Viraj Manwadkar, Clarke J. Esmerian

Abstract: We use the Auriga simulations to probe different satellite quenching mechanisms operating at different mass scales ($10^5 M_\odot \lesssim M_\star \lesssim 10^{11} M_\odot$) in Milky Way-like hosts. Our goal is to understand the origin of the satellite colour distribution and star-forming properties in both observations and simulations. We find that the satellite populations in the Auriga simulati… ▽ More We use the Auriga simulations to probe different satellite quenching mechanisms operating at different mass scales ($10^5 M_\odot \lesssim M_\star \lesssim 10^{11} M_\odot$) in Milky Way-like hosts. Our goal is to understand the origin of the satellite colour distribution and star-forming properties in both observations and simulations. We find that the satellite populations in the Auriga simulations, which was originally designed to model Milky Way-like host galaxies, resemble the populations in the Exploration of Local VolumE Satellites (ELVES) Survey and the Satellites Around Galactic Analogs (SAGA) survey in their luminosity function in the luminosity range $-12 \lesssim M_V \lesssim -15$ and resemble ELVES in their quenched fraction and colour--magnitude distribution in the luminosity range $-12 \lesssim M_g \lesssim -15$. We find that satellites transition from blue colours to red colours at the luminosity range $-15 \lesssim M_g \lesssim -12$ in both the simulations and observations and we show that this shift is driven by environmental effects in the simulations. We demonstrate also that the colour distribution in both simulations and observations can be decomposed into two statistically distinct populations based on their morphological type or star-forming status that are statistically distinct. In the simulations, these two populations also have statistically distinct infall time distributions. The comparison presented here seems to indicate that the tension between the quenched fraction in SAGA and simulations is resolved by the improved target selection of ELVES, but there are still tensions in understanding the colours of faint galaxies, of which ELVES appears to have a significant population of faint blue satellites not recovered in Auriga. △ Less

Submitted 30 January, 2023; v1 submitted 29 August, 2022; originally announced August 2022.

Comments: 15 pages, 10 figures. Published in MNRAS on December 15, 2022

Searching for nascent planetary nebulae: OHPNe candidates in the SPLASH survey

Authors: Roldán A. Cala, José F. Gómez, Luis F. Miranda, Lucero Uscanga, Shari L. Breen, Joanne R. Dawson, Itziar de Gregorio-Monsalvo, Hiroshi Imai, Hai-Hua Qiao, Olga Suárez

Abstract: The evolution of asymptotic giant branch stars from the spherical symmetry into the diverse shapes of planetary nebulae (PNe) is a topic of intensive research. Young PNe provide a unique opportunity to characterize the onset of this transitional phase. In particular, OH maser-emitting PNe (OHPNe) are considered nascent PNe. In fact, only 6 OHPNe have been confirmed to date. In order to identify an… ▽ More The evolution of asymptotic giant branch stars from the spherical symmetry into the diverse shapes of planetary nebulae (PNe) is a topic of intensive research. Young PNe provide a unique opportunity to characterize the onset of this transitional phase. In particular, OH maser-emitting PNe (OHPNe) are considered nascent PNe. In fact, only 6 OHPNe have been confirmed to date. In order to identify and characterize more OHPNe, we processed the unpublished continuum data of the interferometric follow-up of the Southern Parkes Large-Area Survey in Hydroxyl (SPLASH). We then matched the interferometric positions of OH maser and radio continuum emission, considering the latter as a possible tracer of free-free emission from photoionized gas, characteristic of PNe. We report 8 objects with a positive coincidence, 4 of which are classified as candidate OHPNe here for the first time (IRAS 16372-4808, IRAS 17494-2645, IRAS 18019-2216 and OH 341.6811+00.2634). Available evidence strongly indicates that they are evolved stars, while the comparison with confirmed OHPNe indicates that they are likely to be PNe. Their final confirmation as bona fide PNe, however, requires optical/infrared spectroscopy. The obtained spectral indices of the radio continuum emission (between $\simeq$ 0.4 - 1.3) are consistent with partially optically thick free-free emission from photoionized gas. Also, they cluster in the same region of a WISE colour-colour diagram as that of the confirmed OHPNe ($9.5 \lesssim [3.4]-[22] \lesssim 13.5$, and $4.0 \lesssim [4.6]-[12] \lesssim 7.0$), thus this diagram could help to identify more OHPNe candidates in the future. △ Less

Submitted 19 August, 2022; originally announced August 2022.

Comments: 17 pages, 7 figures, 5 tables. Accepted for publication in MNRAS

Machine Learning for Galactic Archaeology: A chemistry-based neural network method for identification of accreted disc stars

Authors: Thorold Tronrud, Patricia B. Tissera, Facundo A. Gómez, Robert J. J. Grand, Ruediger Pakmor, Federico Marinacci, Christine M. Simpson

Abstract: We develop a method ('Galactic Archaeology Neural Network', GANN) based on neural network models (NNMs) to identify accreted stars in galactic discs by only their chemical fingerprint and age, using a suite of simulated galaxies from the Auriga Project. We train the network on the target galaxy's own local environment defined by the stellar halo and the surviving satellites. We demonstrate that th… ▽ More We develop a method ('Galactic Archaeology Neural Network', GANN) based on neural network models (NNMs) to identify accreted stars in galactic discs by only their chemical fingerprint and age, using a suite of simulated galaxies from the Auriga Project. We train the network on the target galaxy's own local environment defined by the stellar halo and the surviving satellites. We demonstrate that this approach allows the detection of accreted stars that are spatially mixed into the disc. Two performance measures are defined - recovery fraction of accreted stars, and the probability that a star with a positive (accreted) classification is a true-positive result, P(TP). As the NNM output is akin to an assigned probability, we are able to determine positivity based on flexible threshold values that can be adjusted easily to refine the selection of presumed-accreted stars. We find that GANN identifies accreted disc stars within simulated galaxies, with high recovery fraction and/or high P(TP). We also find that stars in Gaia-Enceladus-Sausage (GES) mass systems are over 50% recovered by our NNMs in the majority (18/24) of cases. Additionally, nearly every individual source of accreted stars is detected at 10% or more of its peak stellar mass in the disc. We also demonstrate that a conglomerated NNM, trained on the halo and satellite stars from all of the Auriga galaxies provides the most consistent results, and could prove to be an intriguing future approach as our observational capabilities expand. △ Less

Submitted 19 July, 2022; v1 submitted 13 July, 2022; originally announced July 2022.

Comments: 19 pages, 12 figures

Is the core-cusp problem a matter of perspective: Jeans Anisotropic Modeling against numerical simulations

Authors: Wenting Wang, Ling Zhu, Zhaozhou Li, Yang Chen, Jiaxin Han, Feihong He, Xiaohu Yang, Yipeng Jing, Carlos Frenk, Jialu Nie, Hao Tian, Chao Liu, Yanan Cao, Xiaoqing Qiu, John Helly, Robert J. J. Grand, Facundo A. Gomez

Abstract: Mock member stars for 28 dwarf galaxies are constructed from the cosmological Auriga simulation, which reflect the dynamical status of realistic stellar tracers. The axis-symmetric Jeans Anisotropic Multi-Gaussian Expansion (JAM) modeling is applied to 6,000 star particles for each system, to recover the underlying matter distribution. The stellar or dark matter component individually is poorly re… ▽ More Mock member stars for 28 dwarf galaxies are constructed from the cosmological Auriga simulation, which reflect the dynamical status of realistic stellar tracers. The axis-symmetric Jeans Anisotropic Multi-Gaussian Expansion (JAM) modeling is applied to 6,000 star particles for each system, to recover the underlying matter distribution. The stellar or dark matter component individually is poorly recovered, but the total profile is constrained more reasonably. The mass within the half-mass radius of tracers is recovered the tightest, and the mass between 200 and 300 pc, $M(200-300\mathrm{pc})$, is constrained ensemble unbiasedly, with a scatter of 0.167 dex. If using 2,000 particles and only line-of-sight velocities with typical errors, the scatter in $M(200-300\mathrm{pc})$ is increased by $\sim$50%. Quiescent Sagittarius dSph-like systems and star-forming systems with strong outflows show distinct features, with $M(200-300\mathrm{pc})$ mostly under-estimated for the former, and likely over-estimated for the latter. The biases correlate with the dynamical status, which is a result of contraction motions due to tidal effects in quiescent systems or galactic winds in star-forming systems, driving them out of equilibrium. After including Gaia DR3 proper motion errors, we find proper motions can be as useful as line-of-sight velocities for nearby systems at $<\sim$60 kpc. By extrapolating the actual density profiles and the dynamical constraints down to scales below the resolution, we find the mass within 150 pc can be constrained ensemble unbiasedly, with a scatter of $\sim$0.255 dex. In the end, we show that the contraction of member stars in nearby systems is detectable based on Gaia DR3 proper motion errors. △ Less

Submitted 18 November, 2022; v1 submitted 24 June, 2022; originally announced June 2022.

Comments: accepted by ApJ

WiNDS: An H$α$ kinematics survey of nearby spiral galaxies -- Vertical perturbations in nearby disk-type galaxies

Authors: Catalina Urrejola-Mora, Facundo A. Gómez, Sergio Torres-Flores, Philippe Amram, Benoît Epinat, Antonela Monachesi, Federico Marinacci, Claudia Mendes de Oliveira

Abstract: We present the Waves in Nearby Disk galaxies Survey (WiNDS) consisting of 40 nearby low inclination disk galaxies observed through H$α$ high-resolution Fabry Perot interferometry. WiNDS consists of 12 new galaxy observations and 28 data archived observations obtained from different galaxy surveys. We derive two-dimensional line-of-sight velocity fields that are analyzed to identify the possible pr… ▽ More We present the Waves in Nearby Disk galaxies Survey (WiNDS) consisting of 40 nearby low inclination disk galaxies observed through H$α$ high-resolution Fabry Perot interferometry. WiNDS consists of 12 new galaxy observations and 28 data archived observations obtained from different galaxy surveys. We derive two-dimensional line-of-sight velocity fields that are analyzed to identify the possible presence of vertical velocity flows in the galactic disks of these low-inclination late-type galaxies using velocity residual maps, derived from the subtraction of an axisymmetric rotation model to rotational velocity map. Large and globally coherent flows in the line-of-sight velocity of nearly face-on galaxies can be associated with large vertical displacement of the disk with respect to its mid-plane. Our goal is to characterize how frequent vertical perturbations, such as those observed in the Milky Way, arise in the Local Universe. Our currently available data have allowed us to identify 20$\%$ of WiNDS galaxies with strong velocity perturbations that are consistent with vertically perturbed galactic disks. △ Less

Submitted 20 June, 2022; originally announced June 2022.

Comments: 37 pages, 30 figures. Accepted for publication in ApJ

The impact of two massive early accretion events in a Milky Way-like galaxy: repercussions for the buildup of the stellar disc and halo

Authors: Matthew D. A. Orkney, Chervin F. P. Laporte, Robert J. J. Grand, Facundo A. Gómez, Freeke van de Voort, Federico Marinacci, Ruediger Pakmor, Volker Springel

Abstract: We identify and characterise a Milky Way-like realisation from the Auriga simulations with two consecutive massive mergers $\sim2\,$Gyr apart at high redshift, comparable to the reported Kraken and Gaia-Sausage-Enceladus. The Kraken-like merger ($z=1.6$, $M_{\rm Tot} = 8\times10^{10}\,$M$_{\odot}$) is gas-rich, deposits most of its mass in the inner $10\,$kpc, and is largely isotropic. The Sausage… ▽ More We identify and characterise a Milky Way-like realisation from the Auriga simulations with two consecutive massive mergers $\sim2\,$Gyr apart at high redshift, comparable to the reported Kraken and Gaia-Sausage-Enceladus. The Kraken-like merger ($z=1.6$, $M_{\rm Tot} = 8\times10^{10}\,$M$_{\odot}$) is gas-rich, deposits most of its mass in the inner $10\,$kpc, and is largely isotropic. The Sausage-like merger ($z=1.14$, $M_{\rm Tot} = 1\times10^{11}\,$M$_{\odot}$) leaves a more extended mass distribution at higher energies, and has a radially anisotropic distribution. For the higher redshift merger, the stellar mass ratio of the satellite to host galaxy is 1:3. As a result, the chemistry of the remnant is indistinguishable from contemporaneous in-situ populations, making it challenging to identify this component through chemical abundances. This naturally explains why all abundance patterns attributed so far to Kraken are in fact fully consistent with the metal-poor in-situ so-called Aurora population and thick disc. However, our model makes a falsifiable prediction: if the Milky Way underwent a gas-rich double merger at high redshift, then this should be imprinted on its star formation history with bursts about $\sim2\,$Gyrs apart. This may offer constraining power on the highest-redshift major mergers. △ Less

Submitted 7 November, 2022; v1 submitted 18 June, 2022; originally announced June 2022.

Comments: 5 pages, 5 figures, published by MNRAS letters. Comments welcome

Journal ref: Volume 517, Issue 1, November 2022, Pages L138-L142

Possibilities and Limitations of Kinematically Identifying Stars from Accreted Ultra-Faint Dwarf Galaxies

Authors: Kaley Brauer, Hillary Diane Andales, Alexander P. Ji, Anna Frebel, Mohammad K. Mardini, Facundo A. Gomez, Brian W. O'Shea

Abstract: The Milky Way has accreted many ultra-faint dwarf galaxies (UFDs), and stars from these galaxies can be found throughout our Galaxy today. Studying these stars provides insight into galaxy formation and early chemical enrichment, but identifying them is difficult. Clustering stellar dynamics in 4D phase space ($E$, $L_z$, $J_r$, $J_z$) is one method of identifying accreted structure which is curre… ▽ More The Milky Way has accreted many ultra-faint dwarf galaxies (UFDs), and stars from these galaxies can be found throughout our Galaxy today. Studying these stars provides insight into galaxy formation and early chemical enrichment, but identifying them is difficult. Clustering stellar dynamics in 4D phase space ($E$, $L_z$, $J_r$, $J_z$) is one method of identifying accreted structure which is currently being utilized in the search for accreted UFDs. We produce 32 simulated stellar halos using particle tagging with the \textit{Caterpillar} simulation suite and thoroughly test the abilities of different clustering algorithms to recover tidally disrupted UFD remnants. We perform over 10,000 clustering runs, testing seven clustering algorithms, roughly twenty hyperparameter choices per algorithm, and six different types of data sets each with up to 32 simulated samples. Of the seven algorithms, HDBSCAN most consistently balances UFD recovery rates and cluster realness rates. We find that even in highly idealized cases, the vast majority of clusters found by clustering algorithms do not correspond to real accreted UFD remnants and we can generally only recover $6\%$ of UFDs remnants at best. These results focus exclusively on groups of stars from UFDs, which have weak dynamic signatures compared to the background of other stars. The recoverable UFD remnants are those that accreted recently, $z_{\text{accretion}}\lesssim 0.5$. Based on these results, we make recommendations to help guide the search for dynamically-linked clusters of UFD stars in observational data. We find that real clusters generally have higher median energy and $J_r$, providing a way to help identify real vs. fake clusters. We also recommend incorporating chemical tagging as a way to improve clustering results. △ Less

Submitted 1 August, 2022; v1 submitted 14 June, 2022; originally announced June 2022.

Comments: 21 pages, 14 figures, accepted to ApJ

The stellar halo in Local Group Hestia simulations III. Chemical abundance relations for accreted and in-situ stars

Authors: Sergey Khoperskov, Ivan Minchev, Noam Libeskind, Vasily Belokurov, Matthias Steinmetz, Facundo A. Gomez, Robert J. J. Grand, Yehuda Hoffman, Alexander Knebe, Jenny G. Sorce, Martin Sparre, Elmo Tempel, Mark Vogelsberger

Abstract: Since the chemical abundances of stars are the fossil records of the physical conditions in galaxies, they provide the key information for recovering the assembly history of galaxies. In this work, we explore the chemo-chrono-kinematics of accreted and in-situ stars, by analyzing six M31/MW analogues from the HESTIA suite of cosmological hydrodynamics zoom-in simulations of the Local Group. We fou… ▽ More Since the chemical abundances of stars are the fossil records of the physical conditions in galaxies, they provide the key information for recovering the assembly history of galaxies. In this work, we explore the chemo-chrono-kinematics of accreted and in-situ stars, by analyzing six M31/MW analogues from the HESTIA suite of cosmological hydrodynamics zoom-in simulations of the Local Group. We found that the merger debris are chemically distinct from the survived dwarf galaxies. The mergers debris have abundances expected for stars originating from dwarfs that had their star formation activity quenched at early times. Accreted stellar haloes, including individual debris, reveal abundance gradients in the ELz, where the most metal-rich stars have formed in the inner parts of the disrupted systems before the merger and mainly contribute to the central regions of the hosts. Therefore, we suggest that abundance measurements in the inner MW will allow constraining better the parameters of building blocks of the MW stellar halo. The MDFs of the individual debris show several peaks and the majority of debris have lower metallicity than the in-situ stars for Lz>0, while non-rotating and retrograde accreted stars are similar to the in-situ. Prograde accreted stars show a prominent knee in the [Fe/H]-[Mg/Fe] plane while the retrograde stars typically deposit to a high-[Mg/Fe] sequence. We found that the metal-poor stars ([Fe/H]<-1) of the HESTIA galaxies exhibit between zero to 80 km/s net rotation which is consistent with the Aurora population. At higher metallicities, we detect a sharp transition (spin-up) from the turbulent phase to a disk-like rotation. Mergers debris are similar in the [Fe/H]-[Mg/Fe] plane. However, combining a set of abundances allows to capture chemical patterns corresponding to different debris, which are the most prominent as a function of stellar age. △ Less

Submitted 12 September, 2023; v1 submitted 11 June, 2022; originally announced June 2022.

Comments: Paper III in series of III, 21 pages, 17 figures, A&A in press

Journal ref: A&A 677, A91 (2023)

The stellar halo in Local Group Hestia simulations II. The accreted component

Authors: Sergey Khoperskov, Ivan Minchev, Noam Libeskind, Misha Haywood, Paola Di Matteo, Vasily Belokurov, Matthias Steinmetz, Facundo A. Gomez, Robert J. J. Grand, Yehuda Hoffman, Alexander Knebe, Jenny G. Sorce, Martin Sparre, Elmo Tempel, Mark Vogelsberger

Abstract: In the Milky Way, recent progress in the exploration of its assembly history is driven by the tremendous amount of high-quality data delivered by Gaia, which has revealed a number of substructures potentially linked to several ancient accretion events. In this work, aiming to explore the phase-space structure of accreted stars, we analyze six M31/MW analogues from the HESTIA suite of cosmological… ▽ More In the Milky Way, recent progress in the exploration of its assembly history is driven by the tremendous amount of high-quality data delivered by Gaia, which has revealed a number of substructures potentially linked to several ancient accretion events. In this work, aiming to explore the phase-space structure of accreted stars, we analyze six M31/MW analogues from the HESTIA suite of cosmological hydrodynamics zoom-in simulations of the Local Group. We found that all the HESTIA galaxies experience a few dozen mergers but only 1-4 mergers have the stellar mass ratio >0.2 where, depending on the halo definition, the most massive merger contributes from 20% to 70% of the total stellar halo. Individual merger remnants show diverse density distributions at z=0, significantly overlapping with each other and with the in-situ stars in the ELz, UV and RVphi coordinates. The mergers debris often change their position in the ELz with time due to the galactic mass growth and the non-axisymmetry of the potential. In agreement with previous works, we show that even individual merger debris exhibit a number of distinct ELz features. In the UV plane, all HESTIA galaxies reveal radially hot, non-rotating or weakly counter-rotating, Gaia-Sausage-like features. We found an age gradient in Elz space for the individual debris, where the youngest stars, formed in the inner regions of accreting systems, deposit to the innermost regions of the host. The bulk of these stars is being formed during the last stages of accretion, making it possible to date the merger. In actions space (Jr, Jz, Jφ), the mergers debris do not appear as isolated substructures but are instead scattered over a large parameters area and overlapping with the in-situ stars. We also introduce a purely kinematic space (Jz/Jr-eccentricity), where different merger debris can be disentangled better from each other and from the in-situ stars. △ Less

Submitted 12 September, 2023; v1 submitted 9 June, 2022; originally announced June 2022.

Comments: Paper II in a series of III, 22 pages, 16 figures, A&A in press

Journal ref: A&A 677, A90 (2023)

The stellar halo in Local Group Hestia simulations I. The in-situ component and the effect of mergers

Authors: Sergey Khoperskov, Ivan Minchev, Noam Libeskind, Misha Haywood, Paola Di Matteo, Vasily Belokurov, Matthias Steinmetz, Facundo A. Gomez, Robert J. J. Grand, Yehuda Hoffman, Alexander Knebe, Jenny G. Sorce, Martin Sparre, Elmo Tempel, Mark Vogelsberger

Abstract: Theory suggests that mergers play an important role in shaping galactic discs and stellar haloes, which was observationally confirmed in the MW thanks to Gaia data. In this work, aiming to probe the contribution of mergers to the in situ stellar halo formation, we analyse six M31/MW analogues from the HESTIA suite of cosmological hydrodynamical zoom-in simulations of the LG. We found that all the… ▽ More Theory suggests that mergers play an important role in shaping galactic discs and stellar haloes, which was observationally confirmed in the MW thanks to Gaia data. In this work, aiming to probe the contribution of mergers to the in situ stellar halo formation, we analyse six M31/MW analogues from the HESTIA suite of cosmological hydrodynamical zoom-in simulations of the LG. We found that all the HESTIA galaxies experience between one to four mergers with stellar mass ratios between 0.2 and 1 relative to the host at the time of the merger. These significant mergers, with a single exception, happened 7-11Gyr ago. The overall impact of the most massive mergers in HESTIA is clearly seen as a sharp increase in the orbital eccentricity (and a corresponding decrease in the rotational velocity Vphi of pre-existing disc stars of the main progenitor, thus nicely reproducing the Splash-, Plume-like feature that was discovered in the MW. We do find a correlation between mergers and close pericentric passages of massive satellites and bursts of star formation in the in situ component. Massive mergers sharply increase the disc velocity dispersion of the in situ stars; however, the latest significant merger often heats up the disc up to the numbers when the contribution of the previous ones is less prominent in the age-velocity dispersion relation. In HESTIA galaxies, the in situ halo is an important component of the inner stellar halo where its fraction is about 30-40%, while in the outer parts it typically does not exceed ~5% beyond 15 kpc. The simulations suggest that this component of the stellar haloes continues to grow well after mergers conclude; however, the most significant contribution comes from stars that formed recently before the merger. The orbital analysis of the HESTIA galaxies suggests that wedges in Rmax-Zmax space are mainly populated by the stars born between significant mergers. △ Less

Submitted 12 September, 2023; v1 submitted 9 June, 2022; originally announced June 2022.

Comments: Paper I in a series of III, 17 pages, 15 figures, A&A in press

Journal ref: A&A 677, A89 (2023)

Smart: A program to automatically compute accelerations and variational equations

Authors: Daniel D. Carpintero, Nicolás P. Maffione, Facundo A. Gómez

Abstract: Modern astronomical potentials modeling galaxies or stellar systems can be rather involved, and deriving their first derivatives (accelerations) and second derivatives (variational equations) in order to compute orbits and their chaoticity may be a formidable task. We present here a fully automated routine, dubbed Smart, with which the accelerations and the variational equations of an arbitrary po… ▽ More Modern astronomical potentials modeling galaxies or stellar systems can be rather involved, and deriving their first derivatives (accelerations) and second derivatives (variational equations) in order to compute orbits and their chaoticity may be a formidable task. We present here a fully automated routine, dubbed Smart, with which the accelerations and the variational equations of an arbitrary potential that has been written in the Fortran 77 language can be computed. Almost any Fortran 77 statement is admitted in the potential, and the output are standard Fortran 77 routines ready to use. We validate our algorithm with a set of potentials including time-dependent, velocity-dependent and very complex potentials that even involve auxiliary routines. We also describe with some detail a realistic seven-component Galactic potential, MilkyWayHydra, which yields very involved derivatives, thus being a good test bed for Smart. △ Less

Submitted 31 May, 2022; originally announced May 2022.

Comments: 37 pages, 5 figures. To be published in Astronomy and Computing, in press

Preparing for low surface brightness science with the Vera C. Rubin Observatory: characterisation of tidal features from mock images

Authors: G. Martin, A. E. Bazkiaei, M. Spavone, E. Iodice, J. C. Mihos, M. Montes, J. A. Benavides, S. Brough, J. L. Carlin, C. A. Collins, P. A. Duc, F. A. Gómez, G. Galaz, H. M. Hernández-Toledo, R. A. Jackson, S. Kaviraj, J. H. Knapen, C. Martínez-Lombilla, S. McGee, D. O'Ryan, D. J. Prole, R. M. Rich, J. Román, E. A. Shah, T. K. Starkenburg , et al. (28 additional authors not shown)

Abstract: Tidal features in the outskirts of galaxies yield unique information about their past interactions and are a key prediction of the hierarchical structure formation paradigm. The Vera C. Rubin Observatory is poised to deliver deep observations for potentially of millions of objects with visible tidal features, but the inference of galaxy interaction histories from such features is not straightforwa… ▽ More Tidal features in the outskirts of galaxies yield unique information about their past interactions and are a key prediction of the hierarchical structure formation paradigm. The Vera C. Rubin Observatory is poised to deliver deep observations for potentially of millions of objects with visible tidal features, but the inference of galaxy interaction histories from such features is not straightforward. Utilising automated techniques and human visual classification in conjunction with realistic mock images produced using the NEWHORIZON cosmological simulation, we investigate the nature, frequency and visibility of tidal features and debris across a range of environments and stellar masses. In our simulated sample, around 80 per cent of the flux in the tidal features around Milky Way or greater mass galaxies is detected at the 10-year depth of the Legacy Survey of Space and Time (30-31 mag / sq. arcsec), falling to 60 per cent assuming a shallower final depth of 29.5 mag / sq. arcsec. The fraction of total flux found in tidal features increases towards higher masses, rising to 10 per cent for the most massive objects in our sample (M*~10^{11.5} Msun). When observed at sufficient depth, such objects frequently exhibit many distinct tidal features with complex shapes. The interpretation and characterisation of such features varies significantly with image depth and object orientation, introducing significant biases in their classification. Assuming the data reduction pipeline is properly optimised, we expect the Rubin Observatory to be capable of recovering much of the flux found in the outskirts of Milky Way mass galaxies, even at intermediate redshifts (z<0.2). △ Less

Submitted 7 May, 2022; v1 submitted 15 March, 2022; originally announced March 2022.

Comments: 29 pages, 25 figures; accepted for publication in MNRAS following minor corrections

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 513, Issue 1, June 2022, Pages 1459-1487,

SPLASH: The Southern Parkes Large-Area Survey in Hydroxyl -- Data Description & Release

Authors: J. R. Dawson, P. A. Jones, C. Purcell, A. J. Walsh, S. L. Breen, C. Brown, E. Carretti, M. R. Cunningham, J. M. Dickey, S. P. Ellingsen, S. J. Gibson, J. F. Gomez, J. A. Green, H. Imai, V. Krishnan, N. Lo, V. Lowe, M. Marquarding, N. M. McClure-Griffiths, .

Abstract: We present the full data release for the Southern Parkes Large-Area Survey in Hydroxyl (SPLASH), a sensitive, unbiased single-dish survey of the Southern Galactic Plane in all four ground-state transitions of the OH radical at 1612, 1665, 1667 and 1720 MHz. The survey covers the inner Galactic Plane, Central Molecular Zone and Galactic Centre over the range $|b|<$ 2$^{\circ}$, 332$^{\circ}$… ▽ More We present the full data release for the Southern Parkes Large-Area Survey in Hydroxyl (SPLASH), a sensitive, unbiased single-dish survey of the Southern Galactic Plane in all four ground-state transitions of the OH radical at 1612, 1665, 1667 and 1720 MHz. The survey covers the inner Galactic Plane, Central Molecular Zone and Galactic Centre over the range $|b|<$ 2$^{\circ}$, 332$^{\circ}$ $< l <$ 10$^{\circ}$, with a small extension between 2$^{\circ}$ $< b <$ 6$^{\circ}$, 358$^{\circ}$ $< l <$ 4$^{\circ}$. SPLASH is the most sensitive large-scale survey of OH to-date, reaching a characteristic root-mean-square sensitivity of $\sim15$ mK for an effective velocity resolution of $\sim0.9$ km/s. The spectral line datacubes are optimised for the analysis of extended, quasi-thermal OH, but also contain numerous maser sources, which have been confirmed interferometrically and published elsewhere. We also present radio continuum images at 1612, 1666 and 1720 MHz. Based on initial comparisons with $^{12}$CO(J=1-0), we find that OH rarely extends outside CO cloud boundaries in our data, but suggest that large variations in CO-to-OH brightness temperature ratios may reflect differences in the total gas column density traced by each. Column density estimation in the complex, continuum-bright Inner Galaxy is a challenge, and we demonstrate how failure to appropriately model sub-beam structure and the line-of-sight source distribution can lead to order-of-magnitude errors. Anomalous excitation of the 1612 and 1720 MHz satellite lines is ubiquitous in the inner Galaxy, but is disabled by line overlap in and around the Central Molecular Zone. △ Less

Submitted 9 March, 2022; originally announced March 2022.

Comments: Main paper: 22 pages, 15 figures. Online-only material (appended to ArXiv PDF): 11 pages, 3 figure sets. Accepted to MNRAS. Survey data is available from: https://docs.datacentral.org.au/splash/

Observational identification of a sample of likely recent Common-Envelope Events

Authors: Theo Khouri, Wouter Vlemmings, Daniel Tafoya, Andrés F. Pérez-Sánchez, Carmen Sánchez Contreras, José F. Gómez, Hiroshi Imai, Raghvendra Sahai

Abstract: One of the most poorly understood stellar evolutionary paths is that of binary systems undergoing common-envelope evolution, when the envelope of a giant star engulfs the orbit of a companion. Although this interaction leads to a great variety of astrophysical systems, direct empirical studies are difficult because few objects experiencing common-envelope evolution are known. We present ALMA obser… ▽ More One of the most poorly understood stellar evolutionary paths is that of binary systems undergoing common-envelope evolution, when the envelope of a giant star engulfs the orbit of a companion. Although this interaction leads to a great variety of astrophysical systems, direct empirical studies are difficult because few objects experiencing common-envelope evolution are known. We present ALMA observations towards sources known as water fountains that reveal they had low initial masses ($<4~{\rm M}_\odot$) and ejected a significant fraction of it over less than a few hundred years. The only mechanism able to explain such rapid mass ejection is common-envelope evolution. Our calculations show that the water-fountain sample accounts for a large fraction of the systems in our Galaxy which have just experienced the common-envelope phase. Since water-fountain sources show characteristic fast bipolar outflows, outflows and jets likely play an important role right before, during or immediately after the common-envelope phase. △ Less

Submitted 17 December, 2021; originally announced December 2021.

Comments: Published on Nature Astronomy on Dec 16, 2021: https://rdcu.be/cDlX8

Clash of Titans: a MUSE dynamical study of the extreme cluster merger SPT-CL J0307-6225

Authors: D. Hernández-Lang, A. Zenteno, A. Diaz-Ocampo, H. Cuevas, J. Clancy, H. Prado P., F. Aldás, D. Pallero, R. Monteiro-Oliveira, F. A. Gómez, A. Ramirez, J. Wynter, E. R. Carrasco, G. K. T. Hau, B. Stalder, M. McDonald, M. Bayliss, B. Floyd, G. Garmire, A. Katzenberger, K. J. Kim, M. Klein, G. Mahler, J. L. Nilo Castellon, A. Saro , et al. (1 additional authors not shown)

Abstract: We present VLT/MUSE spectroscopy, along with archival Gemini/GMOS spectroscopy, Magellan/Megacam imaging, and Chandra X-ray emission for SPT-CL J0305-6225, a z=0.58 major merging galaxy cluster with a large BCG-SZ centroid separation and a highly disturbed X-ray morphology. The galaxy density distribution shows two main overdensities with separations of 0.144 and 0.017 arcmin to their respective B… ▽ More We present VLT/MUSE spectroscopy, along with archival Gemini/GMOS spectroscopy, Magellan/Megacam imaging, and Chandra X-ray emission for SPT-CL J0305-6225, a z=0.58 major merging galaxy cluster with a large BCG-SZ centroid separation and a highly disturbed X-ray morphology. The galaxy density distribution shows two main overdensities with separations of 0.144 and 0.017 arcmin to their respective BCGs. We characterize the central regions of the two colliding structures, namely 0307-6225N and 0307-6225S, finding velocity derived masses of $M_{200,N}=$ 2.44 $\pm$ 1.41 $\times10^{14}$ M$_\odot$ and $M_{200,S}=$ 3.16 $\pm$ 1.88 $\times10^{14}$ M$_\odot$, with a line-of-sight velocity difference of $|Δv| = 342$ km s$^{-1}$. The total dynamically derived mass is consistent with the SZ derived mass of 7.63 h$_{70}^{-1}$ $\pm$ 1.36 $\times10^{14}$ M$_\odot$. We model the merger using the Monte Carlo Merger Analysis Code, estimating a merging angle of 36$^{+14}_{-12}$ degrees with respect to the plane of the sky. Comparing with simulations of a merging system with a mass ratio of 1:3, we find that the best scenario is that of an ongoing merger that began 0.96$^{+0.31}_{-0.18}$ Gyr ago. We also characterize the galaxy population using H$δ$ and [OII] $λ3727$ Å\ lines. We find that most of the emission-line galaxies belong to 0307-6225S, close to the X-ray peak position, with a third of them corresponding to red-cluster sequence galaxies, and the rest to blue galaxies with velocities consistent with recent periods of accretion. Moreover, we suggest that 0307-6225S suffered a previous merger, evidenced through the two equally bright BCGs at the center with a velocity difference of $\sim$674 km s$^{-1}$. △ Less

Submitted 18 January, 2023; v1 submitted 30 November, 2021; originally announced November 2021.

Comments: 23 pages, 8 tables, 20 figures, accepted for publication in MNRAS; Updated to match published version

Journal ref: Volume 517, 2022, pp.4355-4378

Exploring the outskirts of the EAGLE disc galaxies

Authors: Silvio Varela-Lavin, Patricia B. Tissera, Facundo A. Gómez, Lucas A. Bignone, Claudia del P. Lagos

Abstract: Observations show that the surface brightness of disc galaxies can be well-described by a single exponential (TI), up-bending (TIII) or down-bending (TII) profiles in the outskirts. Here we characterize the mass surface densities of simulated late-type galaxies from the EAGLE project according to their distribution of mono-age stellar populations, the star formation activity and angular momentum c… ▽ More Observations show that the surface brightness of disc galaxies can be well-described by a single exponential (TI), up-bending (TIII) or down-bending (TII) profiles in the outskirts. Here we characterize the mass surface densities of simulated late-type galaxies from the EAGLE project according to their distribution of mono-age stellar populations, the star formation activity and angular momentum content. We find a clear correlation between the inner scale-lengths and the stellar spin parameter, λ, for all three disc types with λ > 0.35. The outer scale-lengths of TII and TIII discs show a positive trend with λ, albeit weaker for the latter. TII discs prefer fast rotating galaxies. With regards to the stellar age distribution, negative and U-shape age profiles are the most common for all disc types. Positive age profiles are determined by a more significant contributions of young stars in the central regions, which decrease rapidly in the outer parts. TII discs prefer relative higher contributions of old stars compared to other mono-age populations across the discs whereas TIII discs become progressively more dominated by intermediate age (2-6 Gyrs) stars for increasing radius. The change in slope of the age profiles is located after the break of the mass surface density. We find evidence of larger flaring for the old stellar populations in TIII systems compared to TI and TII, which could indicate the action of other processes. Overall, the relative distributions of mono-age stellar populations and the dependence of the star formation activity on radius are found to shape the different disc types and age profiles. △ Less

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

Comments: 15 pages, 11 figures, 2 Tables. Accepted for the publication in MNRAS

High and low Sérsic index bulges in Milky Way- and M31-like galaxies: origin and connection to the bar with TNG50

Authors: Ignacio D. Gargiulo, Antonela Monachesi, Facundo A. Gómez, Dylan Nelson, Annalisa Pillepich, Rüdiger Pakmor, R. J. J. Grand, Francesca Fragkoudi, Lars Hernquist, Mark Lovell, Federico Marinacci

Abstract: We study bulge formation in MW/M31-like galaxies in a $Λ$-cold dark matter scenario, focusing on the origin of high- and low-Sersic index bulges. For this purpose we use TNG50, a simulation of the IllustrisTNG project that combines a resolution of $\sim 8 \times 10^4 M_{\odot}$ in stellar particles with a cosmological volume 52 cMpc in extent. We parametrize bulge surface brightness profiles by th… ▽ More We study bulge formation in MW/M31-like galaxies in a $Λ$-cold dark matter scenario, focusing on the origin of high- and low-Sersic index bulges. For this purpose we use TNG50, a simulation of the IllustrisTNG project that combines a resolution of $\sim 8 \times 10^4 M_{\odot}$ in stellar particles with a cosmological volume 52 cMpc in extent. We parametrize bulge surface brightness profiles by the Sérsic index and the bulge-to-total (B/T) ratio obtained from two-component photometric decompositions. In our sample of 287 MW/M31-like simulated galaxies, $17.1\%$ of photometric bulges exhibit high-Sérsic indices and $82.9\%$ show low-Sérsic indices. We study the impact that the environment, mergers and bars have in shaping the surface brightness profiles. We explore two different definitions for local environment and find no correlation between bulge properties and the environment where they reside. Simulated galaxies with higher Sérsic indices show, on average, a higher fraction of ex-situ stars in their kinematically selected bulges. For this bulge population the last significant merger (total mass ratio $m_{\rm sat}/m_{\rm host} > 0.1$) occurs, on average, at later times. However, a substantial fraction of low-Sérsic index bulges also experience a late significant merger. We find that bars play an important role in the development of the different types of photometric bulges. We show that the fraction of simulated galaxies with strong bars is smaller for the high- than for the low-Sérsic index population, reaching differences of $20\%$ at $z > 1$. Simulated galaxies with high fractions of ex-situ stars in the bulge do not develop strong bars. Conversely, simulated galaxies with long-lived strong bars have bulges with ex-situ fractions, $f_{\rm ex-situ} < 0.2$. △ Less

Submitted 26 November, 2021; originally announced November 2021.

Comments: 19 pages, 15 figures (+1 in the appendix) , submitted to MNRAS

The Physical Properties of the SVS 13 Protobinary System: Two Circumstellar Disks and a Spiraling Circumbinary Disk in the Making

Authors: Ana K. Diaz-Rodriguez, Guillem Anglada, Guillermo Blázquez-Calero, Mayra Osorio, José F. Gómez, Gary A. Fuller, Robert Estalella, José M. Torrelles, Sylvie Cabrit, Luis F. Rodríguez, Charlène Lefèvre, Enrique Macías, Carlos Carrasco-González, Luis A. Zapata, Itziar de Gregorio-Monsalvo, Paul T. P. Ho

Abstract: We present VLA and ALMA observations of the close (0.3" = 90 au separation) protobinary system SVS 13. We detect two small circumstellar disks (radii $\sim$12 and $\sim$9 au in dust, and $\sim$30 au in gas) with masses of $\sim$0.004-0.009 $M_{sun}$ for VLA 4A (the western component) and $\sim$0.009-0.030 $M_{sun}$ for VLA 4B (the eastern component). A circumbinary disk with prominent spiral arms… ▽ More We present VLA and ALMA observations of the close (0.3" = 90 au separation) protobinary system SVS 13. We detect two small circumstellar disks (radii $\sim$12 and $\sim$9 au in dust, and $\sim$30 au in gas) with masses of $\sim$0.004-0.009 $M_{sun}$ for VLA 4A (the western component) and $\sim$0.009-0.030 $M_{sun}$ for VLA 4B (the eastern component). A circumbinary disk with prominent spiral arms extending $\sim$500 au and a mass of $\sim$0.052 $M_{sun}$ appears to be in the earliest stages of formation. The dust emission is more compact and with a very high optical depth toward VLA 4B, while toward VLA 4A the dust column density is lower, allowing the detection of stronger molecular transitions. We infer rotational temperatures of $\sim$140 K, on scales of $\sim$30 au, across the whole source, and a rich chemistry. Molecular transitions typical of hot corinos are detected toward both protostars, being stronger toward VLA 4A, with several ethylene glycol transitions detected only toward this source. There are clear velocity gradients, that we interpret in terms of infall plus rotation of the circumbinary disk, and purely rotation of the circumstellar disk of VLA 4A. We measured orbital proper motions and determined a total stellar mass of 1 $M_{sun}$. From the molecular kinematics we infer the geometry and orientation of the system, and stellar masses of $\sim$0.26 $M_{sun}$ for VLA 4A and $\sim$0.60 $M_{sun}$ for VLA 4B. △ Less

Submitted 29 December, 2021; v1 submitted 23 November, 2021; originally announced November 2021.

Comments: 59 pages, 28 Figures, 14 Tables, accepted by The Astrophysical Journal. Includes changes after the ApJ proof corrections

GASKAP-HI Pilot Survey Science I: ASKAP Zoom Observations of HI Emission in the Small Magellanic Cloud

Authors: N. M. Pingel, J. Dempsey, N. M. McClure-Griffiths, J. M. Dickey, K. E. Jameson, H. Arce, G. Anglada, J. Bland-Hawthorn, S. L. Breen, F. Buckland-Willis, S. E. Clark, J. R. Dawson, H. Dénes, E. M. Di Teodoro, B. -Q. For, Tyler J. Foster, J. F. Gómez, H. Imai, G. Joncas, C. -G. Kim, M. -Y. Lee, C. Lynn, D. Leahy, Y. K. Ma, A. Marchal , et al. (31 additional authors not shown)

Abstract: We present the most sensitive and detailed view of the neutral hydrogen (HI) emission associated with the Small Magellanic Cloud (SMC), through the combination of data from the Australian Square Kilometre Array Pathfinder (ASKAP) and Parkes (Murriyang), as part of the Galactic Australian Square Kilometre Array Pathfinder (GASKAP) pilot survey. These GASKAP-HI pilot observations, for the first time… ▽ More We present the most sensitive and detailed view of the neutral hydrogen (HI) emission associated with the Small Magellanic Cloud (SMC), through the combination of data from the Australian Square Kilometre Array Pathfinder (ASKAP) and Parkes (Murriyang), as part of the Galactic Australian Square Kilometre Array Pathfinder (GASKAP) pilot survey. These GASKAP-HI pilot observations, for the first time, reveal HI in the SMC on similar physical scales as other important tracers of the interstellar medium, such as molecular gas and dust. The resultant image cube possesses an rms noise level of 1.1 K (1.6 mJy/beam) per 0.98 km s$^{-1}$ spectral channel with an angular resolution of 30$''$ ($\sim$10 pc). We discuss the calibration scheme and the custom imaging pipeline that utilizes a joint deconvolution approach, efficiently distributed across a computing cluster, to accurately recover the emission extending across the entire $\sim$25 deg$^2$ field-of-view. We provide an overview of the data products and characterize several aspects including the noise properties as a function of angular resolution and the represented spatial scales by deriving the global transfer function over the full spectral range. A preliminary spatial power spectrum analysis on individual spectral channels reveals that the power-law nature of the density distribution extends down to scales of 10 pc. We highlight the scientific potential of these data by comparing the properties of an outflowing high velocity cloud with previous ASKAP+Parkes HI test observations. △ Less

Submitted 10 December, 2021; v1 submitted 9 November, 2021; originally announced November 2021.

Comments: Accepted for publication in PASA, 34 pages, 18 figures, 5 tables

The effects of AGN feedback on the structural and dynamical properties of Milky Way-mass galaxies in cosmological simulations

Authors: Dimitrios Irodotou, Francesca Fragkoudi, Ruediger Pakmor, Robert J. J. Grand, Dimitri A. Gadotti, Tiago Costa, Volker Springel, Facundo A. Gómez, Federico Marinacci

Abstract: Feedback from active galactic nuclei (AGN) has become established as a fundamental process in the evolution of the most massive galaxies. Its impact on Milky Way (MW)-mass systems, however, remains comparatively unexplored. In this work, we use the Auriga simulations to probe the impact of AGN feedback on the dynamical and structural properties of galaxies, focussing on the bar, bulge, and disc. W… ▽ More Feedback from active galactic nuclei (AGN) has become established as a fundamental process in the evolution of the most massive galaxies. Its impact on Milky Way (MW)-mass systems, however, remains comparatively unexplored. In this work, we use the Auriga simulations to probe the impact of AGN feedback on the dynamical and structural properties of galaxies, focussing on the bar, bulge, and disc. We analyse three galaxies -- two strongly and one unbarred/weakly barred -- using three setups: (i) the fiducial Auriga model, which includes both radio and quasar mode feedback, (ii) a setup with no radio mode, and (iii) one with neither the radio nor the quasar mode. When removing the radio mode, gas in the circumgalactic medium cools more efficiently and subsequently settles in an extended disc, with little effect on the inner disc. Contrary to previous studies, we find that although the removal of the quasar mode results in more massive central components, these are in the form of compact discs, rather than spheroidal bulges. Therefore, galaxies without quasar mode feedback are more baryon-dominated and thus prone to forming stronger and shorter bars, which reveals an anti-correlation between the ejective nature of AGN feedback and bar strength. Hence, we report that the effect of AGN feedback (i.e. ejective or preventive) can significantly alter the dynamical properties of MW-like galaxies. Therefore, the observed dynamical and structural properties of MW-mass galaxies can be used as additional constraints for calibrating the efficiency of AGN feedback models. △ Less

Submitted 11 May, 2022; v1 submitted 21 October, 2021; originally announced October 2021.

Comments: 21 pages, 11 main (9 old and 2 new) + 5 Appendix figures. Accepted for publications by MNRAS

The Clustering of Orbital Poles Induced by the LMC: Hints for the Origin of Planes of Satellites

Authors: Nicolas Garavito-Camargo, Ekta Patel, Gurtina Besla, Adrian M. Price-Whelan, Facundo A. Gomez, Chervin F. P Laporte, Kathryn V. Johnston

Abstract: A significant fraction of Milky Way (MW) satellites exhibit phase-space properties consistent with a coherent orbital plane. Using tailored N--body simulations of a spherical MW halo that recently captured a massive (1.8$\times 10^{11}$M$\odot$) LMC-like satellite, we identify the physical mechanisms that may enhance the clustering of orbital poles of objects orbiting the MW. The LMC deviates the… ▽ More A significant fraction of Milky Way (MW) satellites exhibit phase-space properties consistent with a coherent orbital plane. Using tailored N--body simulations of a spherical MW halo that recently captured a massive (1.8$\times 10^{11}$M$\odot$) LMC-like satellite, we identify the physical mechanisms that may enhance the clustering of orbital poles of objects orbiting the MW. The LMC deviates the orbital poles of MW dark matter (DM) particles from the present-day random distribution. Instead, the orbital poles of particles beyond $R\approx 50$kpc cluster near the present-day orbital pole of the LMC along a sinusoidal pattern across the sky. The density of orbital poles is enhanced near the LMC by a factor $δρ_{max}$=30\%(50\%) with respect to underdense regions, and $δρ_{iso}$=15\%(30\%) relative to the isolated MW simulation (no LMC) between 50-150 kpc (150-300 kpc). The clustering appears after the LMC's pericenter ($\approx$ 50 Myr ago, 49 kpc) and lasts for at least 1 Gyr. Clustering occurs because of three effects: 1) the LMC shifts the velocity and position of the central density of the MW's halo and disk; 2) the DM dynamical friction wake and collective response induced by the LMC changes the kinematics of particles; 3) observations of particles selected within spatial planes suffer from a bias, such that measuring orbital poles in a great circle in the sky enhances the probability of their orbital poles being clustered. This scenario should be ubiquitous in hosts that recently captured a massive satellite (at least $\approx$ 1:10 mass ratio), causing the clustering of orbital poles of halo tracers. △ Less

Submitted 16 August, 2021; originally announced August 2021.

Comments: Submitted to ApJ. 23 pages and 16 figures. Comments are welcome

Satellites Around Milky Way Analogs: Tension in the Number and Fraction of Quiescent Satellites Seen in Observations Versus Simulations

Authors: Ananthan Karunakaran, Kristine Spekkens, Kyle A. Oman, Christine M. Simpson, Azadeh Fattahi, David J. Sand, Paul Bennet, Denija Crnojević, Carlos S. Frenk, Facundo A. Gómez, Robert J. J. Grand, Michael G. Jones, Federico Marinacci, Burçin Mutlu-Pakdil, Julio F. Navarro, Dennis Zaritsky

Abstract: We compare the star-forming properties of satellites around Milky Way (MW) analogs from the Stage~II release of the Satellites Around Galactic Analogs Survey (SAGA-II) to those from the APOSTLE and Auriga cosmological zoom-in simulation suites. We use archival GALEX UV imaging as a star-formation indicator for the SAGA-II sample and derive star-formation rates (SFRs) to compare with those from APO… ▽ More We compare the star-forming properties of satellites around Milky Way (MW) analogs from the Stage~II release of the Satellites Around Galactic Analogs Survey (SAGA-II) to those from the APOSTLE and Auriga cosmological zoom-in simulation suites. We use archival GALEX UV imaging as a star-formation indicator for the SAGA-II sample and derive star-formation rates (SFRs) to compare with those from APOSTLE and Auriga. We compare our detection rates from the NUV and FUV bands to the SAGA-II H$α$ detections and find that they are broadly consistent with over $85\%$ of observed satellites detected in all three tracers. We apply the same spatial selection criteria used around SAGA-II hosts to select satellites around the MW-like hosts in APOSTLE and Auriga. We find very good overall agreement in the derived SFRs for the star-forming satellites as well as the number of star-forming satellites per host in observed and simulated samples. However, the number and fraction of quenched satellites in the SAGA-II sample are significantly lower than those in APOSTLE and Auriga below a stellar mass of $M_*\sim10^{8}\,M_{\odot}$, even when the SAGA-II incompleteness and interloper corrections are included. This discrepancy is robust with respect to the resolution of the simulations and persists when alternative star-formation tracers are employed. We posit that this disagreement is not readily explained by vagaries in the observed or simulated samples considered here, suggesting a genuine discrepancy that may inform the physics of satellite populations around MW analogs. △ Less

Submitted 19 July, 2021; v1 submitted 19 May, 2021; originally announced May 2021.

Comments: Accepted for publication in ApJ Letters after minor changes to text

Linking the brightest stellar streams with the accretion history of Milky Way-like galaxies

Authors: Alex Vera-Casanova, Facundo A. Gómez, Antonela Monachesi, Ignacio Gargiulo, Diego Pallero, Robert J. J. Grand, Federico Marinacci, Rüdiger Pakmor, Christine M. Simpson, Carlos S. Frenk, Gustavo Morales

Abstract: According to the current galaxy formation paradigm, mergers and interactions play an important role in shaping present-day galaxies. The remnants of this merger activity can be used to constrain galaxy formation models. In this work we use a sample of thirty hydrodynamical simulations of Milky Way-mass halos, from the AURIGA project, to generate surface brightness maps and search for the brightest… ▽ More According to the current galaxy formation paradigm, mergers and interactions play an important role in shaping present-day galaxies. The remnants of this merger activity can be used to constrain galaxy formation models. In this work we use a sample of thirty hydrodynamical simulations of Milky Way-mass halos, from the AURIGA project, to generate surface brightness maps and search for the brightest stream in each halo as a function of varying limiting magnitude. We find that none of the models shows signatures of stellar streams at $μ_{r}^{lim} \leq 25$ mag arcsec$^{-2}$. The stream detection increases significantly between 27 and 28 mag arcsec$^{-2}$. Nevertheless, even at 30 mag arcsec$^{-2}$, 13 percent of our models show no detectable streams. We study the properties of the brightest streams progenitors (BSPs). We find that BSPs are accreted within a broad range of infall times, from 1.6 to 10 Gyr ago, with only 25 percent accreted within the last 5 Gyrs; thus most BSPs correspond to relatively early accretion events. We also find that 37 percent of the BSPs survive to the present day. The median infall times for surviving and disrupted BSPs are 5.6 and 6.7 Gyr, respectively. We find a clear relation between infall time and infall mass of the BSPs, such that more massive progenitors tend to be accreted at later times. However, we find that the BSPs are not, in most cases, the dominant contributor to the accreted stellar halo of each galaxy. △ Less

Submitted 9 June, 2022; v1 submitted 13 May, 2021; originally announced May 2021.

Comments: Accepted for publication in MNRAS, 15 pages, 11 figures

The reactivation of water maser emission in the planetary nebula IRAS 18061--2505 through a born-again episode

Authors: L. F. Miranda, O. Suárez, L. Olguín, R. Vázquez, L. Sabin, P. F. Guillén, J. F. Gómez, L. Uscanga, P. García-Lario, I. de Gregorio-Monsalvo, A. Aller, A. Manchado, P. Boumis, H. Riesgo, J. M. Matías

Abstract: Water maser emitting planetary nebulae (H$_2$O-PNe) are believed to be among the youngest PNe. We present new optical narrow- and broad-band images, intermediate- and high-resolution long-slit spectra, and archival optical images of the H$_2$O-PN IRAS 18061--2505. It appears a pinched-waist bipolar PN consisting of knotty lobes with some point-symmetric regions, a bow-shock near the tip of each lo… ▽ More Water maser emitting planetary nebulae (H$_2$O-PNe) are believed to be among the youngest PNe. We present new optical narrow- and broad-band images, intermediate- and high-resolution long-slit spectra, and archival optical images of the H$_2$O-PN IRAS 18061--2505. It appears a pinched-waist bipolar PN consisting of knotty lobes with some point-symmetric regions, a bow-shock near the tip of each lobe, and a very compact inner nebula where five components are identified in the spectra by their kinematic and emission properties. The water masers most probably reside in an oxygen-rich ring tracing the equatorial region of the bipolar lobes. These two structures probably result from common envelope evolution plus several bipolar and non-bipolar collimated outflows that have distorted the lobes. The bow-shocks could be related to a previous phase to that of common envelope. The inner nebula may be attributed to a late or very late thermal pulse that occurred before 1951.6 when it was not detectable in the POSSI-Blue image. Chemical abundances and other properties favour a 3--4 M$_{\odot}$ progenitor, although if the common envelope phase accelerated the evolution of the central star, masses <1.5 M$_{\odot}$ cannot be discarded. The age of the bipolar lobes is incompatible with the existence of water masers in IRAS 18061--2505, which may have been lately reactivated through shocks in the oxygen-rich ring, that are generated by the thermal pulse, implying that this PN is not extremely young. We discuss H$_2$O-PNe and possibly related objects in the light of our results for IRAS 18061--2505. △ Less

Submitted 11 May, 2021; originally announced May 2021.

Comments: Submitted to MNRAS (second version incorporating the referee's comments). 20 pages, 11 figures, 4 tables; 3 appendices. The quality of some figues has been reduced

Determining the full satellite population of a Milky Way-mass halo in a highly resolved cosmological hydrodynamic simulation

Authors: Robert J. J. Grand, Federico Marinacci, Rüdiger Pakmor, Christine M. Simpson, Ashley J. Kelly, Facundo A. Gómez, Adrian Jenkins, Volker Springel, Carlos S. Frenk, Simon D. M. White

Abstract: We investigate the formation of the satellite galaxy population of a Milky Way-mass halo in a very highly resolved magneto-hydrodynamic cosmological zoom-in simulation (baryonic mass resolution $m_b =$ 800 $\rm M_{\odot}$). We show that the properties of the central star-forming galaxy, such as the radial stellar surface density profile and star formation history, are: i) robust to stochastic vari… ▽ More We investigate the formation of the satellite galaxy population of a Milky Way-mass halo in a very highly resolved magneto-hydrodynamic cosmological zoom-in simulation (baryonic mass resolution $m_b =$ 800 $\rm M_{\odot}$). We show that the properties of the central star-forming galaxy, such as the radial stellar surface density profile and star formation history, are: i) robust to stochastic variations associated with the so-called ``Butterfly Effect''; and ii) well converged over 3.5 orders of magnitude in mass resolution. We find that there are approximately five times as many satellite galaxies at this high resolution compared to a standard ($m_b\sim 10^{4-5}\, \rm M_{\odot}$) resolution simulation of the same system. This is primarily because 2/3rds of the high resolution satellites do not form at standard resolution. A smaller fraction (1/6th) of the satellites present at high resolution form and disrupt at standard resolution; these objects are preferentially low-mass satellites on intermediate- to low-eccentricity orbits with impact parameters $\lesssim 30$ kpc. As a result, the radial distribution of satellites becomes substantially more centrally concentrated at higher resolution, in better agreement with recent observations of satellites around Milky Way-mass haloes. Finally, we show that our galaxy formation model successfully forms ultra-faint galaxies and reproduces the stellar velocity dispersion, half-light radii, and $V$-band luminosities of observed Milky Way and Local Group dwarf galaxies across 6 orders of magnitude in luminosity ($10^3$-$10^{9}$ $\rm L_{\odot}$). △ Less

Submitted 3 September, 2021; v1 submitted 10 May, 2021; originally announced May 2021.

Comments: accepted for publication in MNRAS. Main changes include new figure 4 and figure 11

A new analytic ram pressure profile for satellite galaxies

Authors: Cristian A. Vega-Martínez, Facundo A. Gómez, Sofía A. Cora, Tomás Hough

Abstract: We present a new analytic fitting profile to model the ram pressure exerted over satellite galaxies on different environments and epochs. The profile is built using the information of the gas particle distribution in hydrodynamical simulations of groups and clusters of galaxies to measure the ram pressure directly. We show that predictions obtained by a previously introduced $β$-profile model can… ▽ More We present a new analytic fitting profile to model the ram pressure exerted over satellite galaxies on different environments and epochs. The profile is built using the information of the gas particle distribution in hydrodynamical simulations of groups and clusters of galaxies to measure the ram pressure directly. We show that predictions obtained by a previously introduced $β$-profile model can not consistently reproduce the dependence of the ram pressure on halocentric distance and redshift for a given halo mass. It features a systematic underestimation of the predicted ram pressure at high redshifts ($z > 1.5$), which increases towards the central regions of the haloes and it is independent of halo mass, reaching differences larger than two decades for satellites at $r<0.4R_\mathrm{vir}$. This behaviour reverses as redshift decreases, featuring an increasing over-estimation with halocentric distance at $z=0$. As an alternative, we introduce a new universal analytic model for the profiles which can recover the ram pressure dependence on halo mass, halocentric distance and redshift. We analyse the impact of our new profile on galaxy properties by applying a semi-analytic model of galaxy formation and evolution on top of the simulations. We show that galaxies experiencing large amounts of cumulative ram pressure stripping typically have low stellar masses ($M_\star \leq 10^{9.5} \text{M}_\odot$). Besides, their specific star formation histories depend on the ram pressure modelling applied, particularly at high redshifts ($z > 1.5$). △ Less

Submitted 8 October, 2021; v1 submitted 29 January, 2021; originally announced February 2021.

Comments: 15 pages, 8 figures. Accepted for publication in MNRAS