Showing 1–50 of 107 results for author: Bruzual, G

Chemical abundances and ionizing mechanisms in the star-forming double-ring of AM 0644-741 using MUSE data

Authors: V. M. A. Gómez-González, Y. D. Mayya, J. Zaragoza-Cardiel, G. Bruzual, S. Charlot, G. Ramos-Larios, L. M. Oskinova, A. A. C. Sander, S. Reyero Serantes

Abstract: We present the analysis of archival Very Large Telescope (VLT) Multi-Unit Spectroscopic Explorer (MUSE) observations of 179 HII regions in the star-forming double-ring collisional galaxy AM 0644-741 at 98.6 Mpc. We determined ionic abundances of He, N, O and Fe using the direct method for the brightest H II region (ID 39); we report $\log\rm{(\frac{N}{O})}=-1.3\pm0.2$ and… ▽ More We present the analysis of archival Very Large Telescope (VLT) Multi-Unit Spectroscopic Explorer (MUSE) observations of 179 HII regions in the star-forming double-ring collisional galaxy AM 0644-741 at 98.6 Mpc. We determined ionic abundances of He, N, O and Fe using the direct method for the brightest H II region (ID 39); we report $\log\rm{(\frac{N}{O})}=-1.3\pm0.2$ and $12+\log\rm{(\frac{O}{H})}=8.9\pm0.2$. We also find the so-called `blue-bump', broad He II $\lambda4686$, in the spectrum of this knot of massive star-formation; its luminosity being consistent with the presence of $\sim430$ Wolf-Rayet (WR) stars of the Nitrogen late-type. We determined the O abundances for 137 HII regions using the strong-line method; we report a median value of $12+\log\rm{(\frac{O}{H})}=8.5\pm0.8$. The location of three objects, including the WR complex, coincide with that of an Ultra Luminous X-ray source. Nebular He II is not detected in any H II region. We investigate the physical mechanisms responsible for the observed spectral lines using appropriate diagnostic diagrams and ionization models. We find that the H II regions are being photoionized by star clusters with ages $\sim2.5-20$ Myr and ionization potential $-3.5<$$\log\langle U\rangle$$<-3.0$. In these diagrams, a binary population is needed to reproduce the observables considered in this work. △ Less

Submitted 20 February, 2024; originally announced February 2024.

Comments: 20 pages. Accepted in MNRAS

Quantifying the escape of Ly$α$ at $z\approx 5-6$: a census of Ly$α$ escape fraction with H$α$ emitting galaxies spectroscopically confirmed by JWST and VLT/MUSE

Authors: Xiaojing Lin, Zheng Cai, Yunjing Wu, Zihao Li, Fengwu Sun, Xiaohui Fan, Zuyi Chen, Mingyu Li, Fuyan Bian, Yuanhang Ning, Linhua Jiang, Gustavo Bruzual, Stephane Charlot, Jacopo Chevallard

Abstract: JWST provides an unprecedented opportunity for unbiased surveys of H$α$-emitting galaxies at $z>4$ with the NIRCam wide-field slitless spectroscopy (WFSS). In this work, we present a census of Ly$α$ escape fraction ($f_{esc, Lyα}$) of 165 star-forming galaxies at $z=4.9-6.3$ using their H$α$ emission directly measured from FRESCO NIRCam/WFSS data. We search for Ly$α$ emission of each H$α$-emitting… ▽ More JWST provides an unprecedented opportunity for unbiased surveys of H$α$-emitting galaxies at $z>4$ with the NIRCam wide-field slitless spectroscopy (WFSS). In this work, we present a census of Ly$α$ escape fraction ($f_{esc, Lyα}$) of 165 star-forming galaxies at $z=4.9-6.3$ using their H$α$ emission directly measured from FRESCO NIRCam/WFSS data. We search for Ly$α$ emission of each H$α$-emitting galaxy in VLT/MUSE data. The overall $f_{esc, Lyα}$ measured by stacking is $f_{esc, Lyα}$ is $0.090\pm0.006$. We find that $f_{esc, Lyα}$ displays a strong dependence on the observed UV slope ($β_{\rm obs}$) and E(B-V), such that the bluest galaxies ($β_{\rm obs}\sim-2.5$) have the largest escape fractions ($f_{\rm esc, Lyα}\approx0.6$), indicative of the crucial role of dust and gas in modulating the escape of Ly$α$ photons. $f_{esc, Lyα}$ is less well related to other parameters, including the UV luminosity and stellar mass, and the variation in $f_{esc, Lyα}$ with them can be explained by their underlying coupling with E(B-V) or $β_{\rm obs}$. Our results suggest a tentative decline in $f_{esc, Lyα}$ at $z\gtrsim 5$, implying increasing intergalactic medium attenuation towards higher redshift. Furthermore, the dependence of $f_{esc, Lyα}$ on $β_{\rm obs}$ is proportional to that of the ionizing photon escape fraction ($f_{\rm esc, LyC}$), indicating the escape of Ly$α$ and ionizing photon may be regulated by similar physical processes. With $f_{esc, Lyα}$ as a proxy to $f_{\rm esc, LyC}$, we infer that UV-faint ($M_{\rm UV}>-16$) galaxies contribute $>70\%$ of the total ionizing emissivity at $z=5-6$. If these relations hold during the epoch of reionization, UV-faint galaxies can contribute the majority of UV photon budget to reionize the Universe. △ Less

Submitted 18 April, 2024; v1 submitted 17 January, 2024; originally announced January 2024.

Comments: 17 pages, 12 figures; Accepted by the ApJS

Ages and metallicities of globular clusters in M81 using GTC/OSIRIS spectra

Authors: Luis Lomelí-Núñez, Y. D. Mayya, L. H. Rodríguez-Merino, P. A. Ovando, Jairo A. Alzate, D. Rosa-González, B. Cuevas-Otahola, Gustavo Bruzual, Arianna Cortesi, V. M. A Gómez-González, Carlos G. Escudero

Abstract: We here present the results of an analysis of the optical spectroscopy of 42 globular cluster (GC) candidates in the nearby spiral galaxy M81 (3.61~Mpc). The spectra were obtained using the long-slit and MOS modes of the OSIRIS instrument at the 10.4~m Gran Telescopio Canarias (GTC) at a spectral resolution of $\sim$1000. We used the classical H$β$ vs [MgFe]$'$ index diagram to separate genuine ol… ▽ More We here present the results of an analysis of the optical spectroscopy of 42 globular cluster (GC) candidates in the nearby spiral galaxy M81 (3.61~Mpc). The spectra were obtained using the long-slit and MOS modes of the OSIRIS instrument at the 10.4~m Gran Telescopio Canarias (GTC) at a spectral resolution of $\sim$1000. We used the classical H$β$ vs [MgFe]$'$ index diagram to separate genuine old GCs from clusters younger than 3 Gyr. Of the 30 spectra with continuum signal-to-noise ratio $>10$, we confirm 17 objects to be classical GCs (age $>10$~Gyr, $-1.4<$[Fe/H]$<-$0.4), with the remaining 13 being intermediate-age clusters (1-7.5~Gyr). We combined age and metallicity data of other nearby spiral galaxies ($\lesssim18$~Mpc) obtained using similar methodology like the one we have used here to understand the origin of GCs in spiral galaxies in the cosmological context. We find that the metal-poor ([Fe/H]<$-$1) GCs continued to form up to 6~Gyr after the first GCs were formed, with all younger systems (age $<8$~Gyr) being metal-rich. △ Less

Submitted 4 January, 2024; originally announced January 2024.

Comments: 16 pages, 13 figures, 8 tables; Accepted for publication in MNRAS

Nebular emission from young stellar populations including binary stars

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

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

Submitted 13 December, 2023; originally announced December 2023.

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

Variation of the stellar initial mass function in semi-analytical models III: testing the cosmic ray regulated integrated galaxy-wide initial mass function

Authors: Fabio Fontanot, Francesco La Barbera, Gabriella De Lucia, Rachele Cecchi, Lizhi Xie, Michaela Hirschmann, Gustavo Bruzual, Stéphane Charlot, Alexandre Vazdekis

Abstract: In our previous work, we derive the CR-IGIMF: a new scenario for a variable stellar initial mass function (IMF), which combines numerical results on the role played by cosmic rays in setting the thermal state of star forming gas, with the analytical approach of the integrated galaxy-wide IMF. In this work, we study the implications of this scenario for the properties of local Early-Type galaxies (… ▽ More In our previous work, we derive the CR-IGIMF: a new scenario for a variable stellar initial mass function (IMF), which combines numerical results on the role played by cosmic rays in setting the thermal state of star forming gas, with the analytical approach of the integrated galaxy-wide IMF. In this work, we study the implications of this scenario for the properties of local Early-Type galaxies (ETG), as inferred from dynamical, photometric and spectroscopic studies. We implement a library of CR-IGIMF shapes in the framework of the Galaxy Evolution and Assembly (GAEA) model. Our realization includes a derivation of synthetic spectral energy distribution for each model galaxy, allowing a direct derivation of the mass fraction in the mean IMF of low-mass stars (i.e. the dwarf-to-giant ratio - $\rm f_{dg}$) and a comparison with IMF sensitive spectral features. The predictions of the GAEA model implementing the CR-IGIMF confirm our previous findings: it correctly reproduces both the observed excess of z$\sim$0 dynamical mass (mass-to-light ratios) with respect to spectroscopic (photometric) estimates assuming a universal, MW-like, IMF, and the observed increase of [$α$/Fe] ratios with stellar mass in spheroidal galaxies. Moreover, this realization reproduces the increasing trends of $\rm f_{dg}$, and IMF-sensitive line-strengths with velocity dispersion, although the predicted relations are significantly shallower than the observed ones. Our results show that the CR-IGIMF is a promising scenario that reproduces at the same time dynamical, photometric and spectroscopic indications of a varying IMF in local ETGs. The shallow relations found for spectral indices suggest that either a stronger variability as a function of galaxy properties or additional dependences (e.g. as a function of star forming gas metallicity) might be required to match the strength of the observed trends. △ Less

Submitted 18 April, 2024; v1 submitted 21 November, 2023; originally announced November 2023.

Comments: 13 pages, 8 figures, 2 tables, A&A accepted

Journal ref: A&A 686, A302 (2024)

HST FUV Spectroscopy of Super Star Cluster A in the Green Pea Analog Mrk 71: Revealing the Presence of Very Massive Stars

Authors: Linda J. Smith, M. S. Oey, Svea Hernandez, Jenna Ryon, Claus Leitherer, Stephane Charlot, Gustavo Bruzual, Daniela Calzetti, You-Hua Chu, Matthew J. Hayes, Bethan L. James, Anne. E. Jaskot, Goeran Oestlin

Abstract: Mrk 71 is a low metallicity (Z = 0.16 Z_sun) starburst region in the local dwarf galaxy NGC 2366, hosting two super star clusters (SSCs A and B), and is recognized as a Green Pea (GP) analog with SSC A responsible for the GP properties. We present STIS and FOS far-ultraviolet (FUV) spectra of the embedded SSC Mrk 71-A obtained with the Hubble Space Telescope (HST). The STIS FUV spectrum shows the… ▽ More Mrk 71 is a low metallicity (Z = 0.16 Z_sun) starburst region in the local dwarf galaxy NGC 2366, hosting two super star clusters (SSCs A and B), and is recognized as a Green Pea (GP) analog with SSC A responsible for the GP properties. We present STIS and FOS far-ultraviolet (FUV) spectra of the embedded SSC Mrk 71-A obtained with the Hubble Space Telescope (HST). The STIS FUV spectrum shows the characteristic features of very massive stars (VMS, masses > 100 M_sun) and we derive an age of 1+/-1 Myr by comparison with the Charlot & Bruzual suite of spectral population synthesis models with upper mass limits of 300 and 600 M_sun. We compare the STIS spectrum with all known SSC spectra exhibiting VMS signatures: NGC 5253-5, R136a, NGC 3125-A1 and the z = 2.37 Sunburst cluster. We find that the cluster mass-loss rates and wind velocities, as characterized by the C IV P Cygni profiles and the He II emission line strengths, are very similar over Z = 0.16 to 0.4 Z_sun. This agrees with predictions that the optically thick winds of VMS will be enhanced near the Eddington limit and show little metallicity dependence. We find very strong damped Lyman-alpha absorption with log N(H I) = 22.2 cm-2 associated with Mrk 71-A. We discuss the natal environment of this young SSC in terms of radiatively-driven winds, catastrophic cooling and recent models where the cluster is surrounded by highly pressurized clouds with large neutral columns. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: 15 pages accepted for publication in the Astrophysical Journal

A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE): JWST Discovers an Overdensity around a Metal Absorption-selected Galaxy at $z\sim5.5$

Authors: Yunjing Wu, Feige Wang, Zheng Cai, Xiaohui Fan, Kristian Finlator, Jinyi Yang, Joseph F. Hennawi, Fengwu Sun, Jaclyn B. Champagne, Xiaojing Lin, Zihao Li, Zuyi Chen, Eduardo Bañados, George D. Becker, Sarah E. I. Bosman, Gstavo Bruzual, Stephane Charlot, Hsiao-Wen Chen, Jacopo Chevallard, Anna-Christina Eilers, Emanuele Paolo Farina, Xiangyu Jin, Hyunsung D. Jun, Koki Kakiichi, Mingyu Li , et al. (5 additional authors not shown)

Abstract: The launch of ${\it JWST}$ opens a new window for studying the connection between metal-line absorbers and galaxies at the end of the Epoch of Reionization (EoR). Previous studies have detected absorber-galaxy pairs in limited quantities through ground-based observations. To enhance our understanding of the relationship between absorbers and their host galaxies at $z>5$, we utilized the NIRCam Wid… ▽ More The launch of ${\it JWST}$ opens a new window for studying the connection between metal-line absorbers and galaxies at the end of the Epoch of Reionization (EoR). Previous studies have detected absorber-galaxy pairs in limited quantities through ground-based observations. To enhance our understanding of the relationship between absorbers and their host galaxies at $z>5$, we utilized the NIRCam Wide Field Slitless Spectroscopy (WFSS) to search for absorber-associated galaxies by detecting their rest-frame optical emission lines (e.g., [OIII] + H$β$). We report the discovery of a MgII-associated galaxy at $z=5.428$ using data from the ${\it JWST}$ ASPIRE program. The MgII absorber is detected on the spectrum of quasar J0305--3150 with a rest-frame equivalent width of 0.74$\mathring{A}$. The associated galaxy has an [OIII] luminosity of $10^{42.5}\ {\rm erg\ s^{-1}}$ with an impact parameter of 24.9 proper kiloparsecs (pkpc). The joint ${\it HST}$-${\it JWST}$ spectral energy distribution (SED) implies a stellar mass and star-formation rate of ${\rm M_* \approx 10^{8.8}}$ ${\rm M_{\odot}}$, ${\rm SFR}\approx 10\ {\rm M_{\odot}\ yr^{-1}}$. Its [OIII] equivalent width and stellar mass are typical of [OIII] emitters at this redshift. Furthermore, connecting the outflow starting time to the SED-derived stellar age, the outflow velocity of this galaxy is $\sim300\ {\rm km\ s^{-1}}$, consistent with theoretical expectations. We identified six additional [OIII] emitters with impact parameters of up to $\sim300$ pkpc at similar redshifts ($|dv|<1000\ {\rm km\ s^{-1}}$). The observed number is consistent with that in cosmological simulations. This pilot study suggests that systematically investigating the absorber-galaxy connection within the ASPIRE program will provide insights into the metal-enrichment history in the early universe. △ Less

Submitted 8 November, 2023; v1 submitted 28 September, 2023; originally announced September 2023.

Comments: Accepted for publication in ApJL. Main text 8 pages, 4 figures. For more information of the JWST ASPIRE program please check https://aspire-quasar.github.io/index.html

The cluster initial mass function of the M82 disk Super Star Clusters

Authors: Bolivia Cuevas-Otahola, Divakara Mayya, Jesus Arriaga-Hernandez, Ivanio Puerari, Gustavo Bruzual

Abstract: The presence of a population of a large number ($\sim$400) of almost coeval (100--300 Myr) super star clusters (SSCs) in the disk of M82 offers an opportunity to construct the Cluster Initial Mass Function (CIMF) from the observed present-day Cluster Mass Function (CMF). We carry out the dynamical and photometric evolution of the CMF assuming the clusters move in circular orbits under the gravitat… ▽ More The presence of a population of a large number ($\sim$400) of almost coeval (100--300 Myr) super star clusters (SSCs) in the disk of M82 offers an opportunity to construct the Cluster Initial Mass Function (CIMF) from the observed present-day Cluster Mass Function (CMF). We carry out the dynamical and photometric evolution of the CMF assuming the clusters move in circular orbits under the gravitational potential of the host galaxy using the semi-analytical simulation code EMACSS. We explore power-law and log-normal functions for the CIMFs, and populate the clusters in the disk assuming uniform, power-law, and exponential radial distribution functions. We find that the observed CMF is best produced by a CIMF that is power-law in form with an index of 1.8, for a power-law radial distribution function. More importantly, we establish that the observed turn-over in the present-day CMF is the result of observational incompleteness rather than due to dynamically induced effects, or an intrinsically log-normal CIMF, as was proposed for the fossil starburst region B of this galaxy. Our simulations naturally reproduce the mass-radius relation observed for a sub-sample of M82 SSCs. △ Less

Submitted 1 September, 2023; originally announced September 2023.

Comments: 17 pages, 11 figures, accepted to be published on MNRAS

MSC Class: 85 ACM Class: J.2

On the origin of optical and IR emission lines in star forming galaxies

Authors: Mariela Martínez-Paredes, Gustavo Bruzual, Christophe Morisset, Minsun Kim, Marcio Meléndez, Luc Binette

Abstract: Combining the {\sc Cloudy} photoionization code with updated stellar population synthesis results, we simultaneously model the MIR $\neiii/\neii$ vs. $\oiv/\neiii$, the MIR-FIR $\neiii/\neii$ vs. $\oiv/\oiii$ and the classical BPT diagnostic diagrams. We focus on the properties of optically classified \hii\,galaxies that lie in the normal star forming zone in the MIR diagnostic diagram. We find th… ▽ More Combining the {\sc Cloudy} photoionization code with updated stellar population synthesis results, we simultaneously model the MIR $\neiii/\neii$ vs. $\oiv/\neiii$, the MIR-FIR $\neiii/\neii$ vs. $\oiv/\oiii$ and the classical BPT diagnostic diagrams. We focus on the properties of optically classified \hii\,galaxies that lie in the normal star forming zone in the MIR diagnostic diagram. We find that a small fraction of our models lie in this zone, but most of them correspond to the lowest explored metallicity, \zstar\,=\,0.0002, at age $\sim1$ Gyr. This value of \zstar\,is, by far, lower than the values derived for these galaxies from optical emission lines, suggesting that the far-UV emission produced by post-AGB stars (a.k.a. HOLMES, hot low-mass evolved stars) is NOT the source of ionization. Instead, shock models can easily reproduce this part of the MIR diagram. We suggest that it is likely that some of these galaxies have been misclassified and that in them, shocks, produced by a weak AGN-outflow, could be an important source of ionizaton. Using a subset of our models, we derive a new demarcation line for the maximal contribution of retired galaxies in the BPT diagram. This demarcation line allows for a larger contamination from the neighbouring AGN-dominated region. Considering the importance of disentangling the different ionising mechanisms in weak or deeply obscured systems, new observational efforts to classify galaxies both in the optical and IR are required to better constrain this kind of models and understand their evolutionary paths. △ Less

Submitted 10 August, 2023; originally announced August 2023.

Comments: Accepted for publication in MNRAS

Extreme broad He\2 emission at high and low redshifts: the dominant role of VMS in NGC 3125-A1 and CDFS131717

Authors: Aida Wofford, Andrés Sixtos, Stephane Charlot, Gustavo Bruzual, Fergus Cullen, Thomas M. Stanton, Svea Hernández, Linda J. Smith, Matthew Hayes

Abstract: Super star cluster (SSC) A1 (3.1E5 Msun) in NGC 3125 has one of the strongest (EW = 4.6 +/- 0.5 Ang) broad (FWHM = 1131 +\- 40 km/s) He II 1640 emission lines in the nearby Universe and constitutes an important template for interpreting observations of extreme He II emitters out to redshifts of z = 2-3. We use Cosmic Origins Spectrograph (COS) observations of A1 to show that there is no significan… ▽ More Super star cluster (SSC) A1 (3.1E5 Msun) in NGC 3125 has one of the strongest (EW = 4.6 +/- 0.5 Ang) broad (FWHM = 1131 +\- 40 km/s) He II 1640 emission lines in the nearby Universe and constitutes an important template for interpreting observations of extreme He II emitters out to redshifts of z = 2-3. We use Cosmic Origins Spectrograph (COS) observations of A1 to show that there is no significant contamination of the He II line with nebular emission and that the line is redshifted by 121 +/-17 km/s relative to ISM lines. We compare the COS G130M + G160M observations of A1 to recent binary BPASS and single-star Charlot & Bruzual (C&B) simple stellar population (SSP) models with Very Massive Stars (VMS) of up to 300 Msun. We suggest why BPASS models fail to reproduce A1's He II emission. On the other hand, a C&B model with Z = 0.008, age = 2.2 Myr, and VMS approaching the Eddington limit provides an excellent fit to the He II emission and fits reasonably well C III 1175, N V 1238,1241, and C IV 1548, 1551. We present O V 1371 line-profile predictions showing that this line constitutes an important tracer of youth and VMS in galaxies. Finally, we discuss the presence of VMS in CDFS131717, a highly star-forming low-metallicity galaxy located at z = 3.071, which has a tentative detection of O V absorption and strong broad He II emission. These features are rare and hint to the presence of short-lived VMS in the galaxy. Our results show the effect of the latest developments of stellar wind theory and the importance of accounting for VMS in models. △ Less

Submitted 23 May, 2023; originally announced May 2023.

Comments: 17 pages, 17 figures, accepted in MNRAS

Constraints on star formation in Orion A from Gaia

Authors: Jairo A. Alzate, Gustavo Bruzual, Marina Kounkel, Gladis Magris, Lee Hartmann, Nuria Calvet, Lyra Cao

Abstract: We develop statistical methods within a Bayesian framework to infer the star formation history from photometric surveys of pre-main sequence populations. Our procedures include correcting for biases due to extinction in magnitude-limited surveys, and using distributions from subsets of stars with individual extinction measurements. We also make modest corrections for unresolved binaries. We apply… ▽ More We develop statistical methods within a Bayesian framework to infer the star formation history from photometric surveys of pre-main sequence populations. Our procedures include correcting for biases due to extinction in magnitude-limited surveys, and using distributions from subsets of stars with individual extinction measurements. We also make modest corrections for unresolved binaries. We apply our methods to samples of populations with Gaia photometry in the Orion A molecular cloud. Using two well-established sets of evolutionary tracks, we find that, although our sample is incomplete at youngest ages due to extinction, star formation has proceeded in Orion A at a relatively constant rate between ages of about 0.3 and 5 Myr, in contrast to other studies suggesting multiple epochs of star formation. Similar results are obtained for a set of tracks that attempt to take the effects of strong magnetic fields into account. We also find no evidence for a well-constrained "birthline" that would result from low-mass stars appearing first along the deuterium-burning main sequence, especially using the magnetic evolutionary tracks. While our methods have been developed to deal with Gaia data, they may be useful for analyzing other photometric surveys of star-forming regions. △ Less

Submitted 19 May, 2023; originally announced May 2023.

Report number: 2305.1182

Simulated observations of star formation regions: infrared evolution of globally collapsing clouds

Authors: Jesús M. Jáquez-Domínguez, Roberto Galván-Madrid, Jacopo Fritz, Manuel Zamora-Avilés, Peter Camps, Gustavo Bruzual, Maarten Baes, Yuxin Lin, Enrique Vázquez-Semadeni

Abstract: The direct comparison between hydrodynamical simulations and observations is needed to improve the physics included in the former and test biases in the latter. Post-processing radiative transfer and synthetic observations are now the standard way to do this. We report on the first application of the \texttt{SKIRT} radiative transfer code to simulations of a star-forming cloud. The synthetic obser… ▽ More The direct comparison between hydrodynamical simulations and observations is needed to improve the physics included in the former and test biases in the latter. Post-processing radiative transfer and synthetic observations are now the standard way to do this. We report on the first application of the \texttt{SKIRT} radiative transfer code to simulations of a star-forming cloud. The synthetic observations are then analyzed following traditional observational workflows. We find that in the early stages of the simulation, stellar radiation is inefficient in heating dust to the temperatures observed in Galactic clouds, thus the addition of an interstellar radiation field is necessary. The spectral energy distribution of the cloud settles rather quickly after $\sim3$ Myr of evolution from the onset of star formation, but its morphology continues to evolve for $\sim8$ Myr due to the expansion of \textsc{Hii} regions and the respective creation of cavities, filaments, and ridges. Modeling synthetic \textit{Herschel} fluxes with 1- or 2-component modified black bodies underestimates total dust masses by a factor of $\sim2$. Spatially-resolved fitting recovers up to about $70\%$ of the intrinsic value. This ``missing mass'' is located in a very cold dust component with temperatures below $10$ K, which does not contribute appreciably to the far-infrared flux. This effect could bias real observations if such dust exists in large amounts. Finally, we tested observational calibrations of the SFR based on infrared fluxes and concluded that they are in agreement when compared to the intrinsic SFR of the simulation averaged over $\sim100$ Myr. △ Less

Submitted 10 April, 2023; originally announced April 2023.

The relation between morphology, star formation history, and environment in local universe galaxies

Authors: D. Pérez-Millán, Jacopo Fritz, Rosa A. González-Lópezlira, Alessia Moretti, Bernardo Cervantes Sodi, Benedetta Vulcani, Marco Gullieuszik, Gustavo Bruzual, Stéphane Charlot, Daniela Bettoni

Abstract: The observed properties of galaxies are strongly dependent on both their total stellar mass and their morphology. Furthermore, the environment is known to play a strong role in shaping them. The galaxy population in the local universe that is located in virialized clusters is found to be red, poorly star-forming, and mostly composed of early morphological types. Towards a holistic understanding of… ▽ More The observed properties of galaxies are strongly dependent on both their total stellar mass and their morphology. Furthermore, the environment is known to play a strong role in shaping them. The galaxy population in the local universe that is located in virialized clusters is found to be red, poorly star-forming, and mostly composed of early morphological types. Towards a holistic understanding of the mechanisms that drive galaxy evolution, we exploit the spectrophotometric data from the WINGS and OmegaWINGS local galaxy cluster surveys, and study the role of both the local and the large-scale environments. We attempt to disentangle their effects from the intrinsic characteristics of the galaxies, in shaping the star formation activity at fixed morphological type and stellar mass. Using a sample of field galaxies from the same surveys for comparison, we analyse the effects of the environment, embodied by the local density, clustercentric distance, and close neighbours, respectively, on the star formation histories of cluster galaxies. We find that local effects have a more relevant impact on galaxy stellar properties than the large-scale environment, and that morphology needs to be taken into account to pinpoint the mechanisms that are driving the influence of clusters in galaxy evolution. △ Less

Submitted 13 February, 2023; originally announced February 2023.

Comments: Accepted by MNRAS. 24 pages, 20 figures, including 3 appendices

Journal ref: MNRAS 2023

Detection of He++ ion in the star-forming ring of the Cartwheel using MUSE data and ionizing mechanisms

Authors: Y. D. Mayya, A. Plat, V. M. A. Gómez-González, J. Zaragoza-Cardiel, S. Charlot, G. Bruzual

Abstract: We here report the detection of the nebular HeII4686 line in 32 HII regions in the metal-poor collisional ring galaxy Cartwheel using the Multi-Unit Spectroscopic Explorer (MUSE) dataset. The measured I(HeII4686)/I(Hbeta) ratio varies from 0.004 to 0.07, with a mean value of 0.010+/-0.003. Ten of these 32 HII regions are coincident with the location of an Ultra Luminous X-ray (ULX) source. We used… ▽ More We here report the detection of the nebular HeII4686 line in 32 HII regions in the metal-poor collisional ring galaxy Cartwheel using the Multi-Unit Spectroscopic Explorer (MUSE) dataset. The measured I(HeII4686)/I(Hbeta) ratio varies from 0.004 to 0.07, with a mean value of 0.010+/-0.003. Ten of these 32 HII regions are coincident with the location of an Ultra Luminous X-ray (ULX) source. We used the flux ratios of important diagnostic lines and results of photoionization by Simple Stellar Populations (SSPs) to investigate the likely physical mechanisms responsible for the ionization of He+. We find that the majority of the regions (27) are consistent with photoionization by star clusters in their Wolf-Rayet (WR) phase with initial ionization parameter -3.5<logU<-2.0. Blue Bump (BB), the characteristic feature of the WR stars, however, is not detected in any of the spectra. We demonstrate that this non-detection is due to the relatively low equivalent width (EW) of the BB in metal-poor SSPs, in spite of containing sufficient number of WR stars to reproduce the observed I(HeII4686)/I(Hbeta) ratio of 1.5% at the Cartwheel metallicity of Z=0.004. The HII regions in the WR phase that are coincident with a ULX source do not show line ratios characteristic of ionization by X-ray sources. However, the ULX sources may have a role to play in the ionization of He+ in two (#99, 144) of the five regions that are not in the WR phase. Ionization by radiative shocks along with the presence of channels for the selective leakage of ionizing photons are the likely scenarios in #17 and #148, the two regions with the highest observed I(HeII4686)/I(Hbeta) ratio. △ Less

Submitted 8 January, 2023; originally announced January 2023.

Comments: 24 pages. Accepted in MNRAS

The Mass-Metallicity Relation of Dwarf Galaxies at Cosmic Noon from JWST Observations

Authors: Mingyu Li, Zheng Cai, Fuyan Bian, Xiaojing Lin, Zihao Li, Yunjing Wu, Fengwu Sun, Shiwu Zhang, Emmet Golden-Marx, Zechang Sun, Siwei Zou, Xiaohui Fan, Eiichi Egami, Stephane Charlot, Gustavo Bruzual, Jacopo Chevallard

Abstract: We present a study of the mass-metallicity relation (MZR) of 51 dwarf galaxies ($M_\star\approx 10^{6.5} - 10^{9.5}~M_\odot$) at $z = 2-3$ from the Abell 2744 and SMACS J0723-3732 galaxy cluster fields. These dwarf galaxies are identified and confirmed by deep JWST/NIRISS imaging and slitless grism spectroscopic observations. By taking advantage of the superior performance of JWST and the gravitat… ▽ More We present a study of the mass-metallicity relation (MZR) of 51 dwarf galaxies ($M_\star\approx 10^{6.5} - 10^{9.5}~M_\odot$) at $z = 2-3$ from the Abell 2744 and SMACS J0723-3732 galaxy cluster fields. These dwarf galaxies are identified and confirmed by deep JWST/NIRISS imaging and slitless grism spectroscopic observations. By taking advantage of the superior performance of JWST and the gravitational lensing effect, we extend the previous MZR relation at $z=2-3$ to a much lower mass regime down by $\approx$ 2.5 orders of magnitude as compared with previous studies. We find that the MZR has a shallower slope at the low-mass end ($M_\star<10^{9}~M_\odot$), with a slope turnover point of $\approx$ $10^9~M_\odot$. This implies that the dominating feedback processes in dwarf galaxies may be different from that in massive galaxies. From $z=3$ to $z=2$, the metallicity of the dwarf galaxies is enhanced by $\approx0.09$ dex for a given stellar mass, consistent with the mild evolution found in galaxies with higher mass. Furthermore, we confirm the existence of a fundamental metallicity relation (FMR) between the gas-phase metallicity, stellar mass, and star formation rate in dwarf galaxies at $z=2-3$. Our derived FMR, which has no significant redshift evolution, can be used as a benchmark to understand the origin of the anti-correlation between the SFR and metallicity of dwarf galaxies in the high-$z$ Universe. △ Less

Submitted 22 September, 2023; v1 submitted 2 November, 2022; originally announced November 2022.

Comments: 19 pages, 4 figures, 3 tables, Accepted for publication in the Astrophysical Journal Letters

Journal ref: ApJL 955 L18 (2023)

The relation between globular cluster systems and supermassive black holes in spiral galaxies III. The link to the $M_\bullet-M_\ast$ correlation

Authors: Rosa A. González-Lópezlira, Luis Lomelí-Núñez, Yasna Ordenes-Briceño, Laurent Loinard, Stephen Gwyn, Karla Alamo-Martínez, Gustavo Bruzual, Ariane Lançon, Thomas H. Puzia

Abstract: We continue to explore the relationship between globular cluster total number, $N_{\rm GC}$, and central black hole mass, $M_\bullet$, in spiral galaxies. We present here results for the Sab galaxies NGC 3368, NGC 4736 (M 94) and NGC 4826 (M 64), and the Sm galaxy NGC 4395. The globular cluster (GC) candidate selection is based on the ($u^*$ - $i^\prime$) versus ($i^\prime$ - $K_s$) color-color di… ▽ More We continue to explore the relationship between globular cluster total number, $N_{\rm GC}$, and central black hole mass, $M_\bullet$, in spiral galaxies. We present here results for the Sab galaxies NGC 3368, NGC 4736 (M 94) and NGC 4826 (M 64), and the Sm galaxy NGC 4395. The globular cluster (GC) candidate selection is based on the ($u^*$ - $i^\prime$) versus ($i^\prime$ - $K_s$) color-color diagram, and $i^\prime$-band shape parameters. We determine the $M_\bullet$ versus $N_{\rm GC}$ correlation for these spirals, plus NGC 4258, NGC 253, M 104, M 81, M 31, and the Milky Way. We also redetermine the correlation for the elliptical sample in Harris, Poole, & Harris (2014), with updated galaxy types from Sahu et al. 2019b. Additionally, we derive total stellar galaxy mass, $M_\ast$, from its two-slope correlation with $N_{\rm GC}$ (Hudson, Harris, & Harris 2014), and fit $M_\bullet$ versus $M_\ast$ for both spirals and ellipticals. We obtain log $M_\bullet \propto$ (1.01 $\pm$ 0.13) log $N_{\rm GC}$ for ellipticals, and log $M_\bullet \propto$ (1.64 $\pm$ 0.24) log $N_{\rm GC}$ for late type galaxies (LTG). The linear $M_\bullet$ versus $N_{\rm GC}$ correlation in ellipticals could be due to statistical convergence through mergers, but not the much steeper correlation for LTG. However, in the $M_\bullet$ versus total stellar mass ($M_\ast$) parameter space, with $M_\ast$ derived from its correlation with $N_{\rm GC}$, $M_\bullet \propto$ (1.48 $\pm$ 0.18) log $M_\ast$ for ellipticals, and $M_\bullet \propto$ (1.21 $\pm$ 0.16) log $M_\ast$ for LTG. The observed agreement between ellipticals and LTG in this parameter space may imply that black holes and galaxies co-evolve through "calm" accretion, AGN feedback, and other secular processes. △ Less

Submitted 10 October, 2022; originally announced October 2022.

Comments: 38 pages, 18 figures, 19 tables; accepted for publication in The Astrophysical Journal. arXiv admin note: text overlap with arXiv:1612.05655

Metal-Enriched Neutral Gas Reservoir around a Strongly-lensed, Low-mass Galaxy at $z=4$ Identified by JWST/NIRISS and VLT/MUSE

Authors: Xiaojing Lin, Zheng Cai, Siwei Zou, Zihao Li, Zuyi Chen, Fuyan Bian, Fengwu Sun, Yiping Shu, Yunjing Wu, Mingyu Li, Jianan Li, Xiaohui Fan, J. Xavier Prochaska, Daniel Schaerer, Stephane Charlot, Daniel Espada, Miroslava Dessauges-Zavadsky, Eiichi Egami, Daniel Stark, Kirsten K. Knudsen, Gustavo Bruzual, Jacopo Chevallard

Abstract: Direct observations of low-mass, low-metallicity galaxies at $z\gtrsim4$ provide an indispensable opportunity for detailed inspection of the ionization radiation, gas flow, and metal enrichment in sources similar to those that reionized the Universe. Combining the James Webb Space Telescope (JWST), VLT/MUSE, and ALMA, we present detailed observations of a strongly lensed, low-mass (… ▽ More Direct observations of low-mass, low-metallicity galaxies at $z\gtrsim4$ provide an indispensable opportunity for detailed inspection of the ionization radiation, gas flow, and metal enrichment in sources similar to those that reionized the Universe. Combining the James Webb Space Telescope (JWST), VLT/MUSE, and ALMA, we present detailed observations of a strongly lensed, low-mass ($\approx 10^{7.6}$ ${\rm M}_\odot$) galaxy at $z=3.98$ (also see Vanzella et al. 2022). We identify strong narrow nebular emission, including CIV $λ\lambda1548,1550$, HeII $\lambda1640$, OIII] $λ\lambda1661,1666$, [NeIII] $\lambda3868$, [OII] $\lambda3727$, and Balmer series of Hydrogen from this galaxy, indicating a metal-poor HII region ($\lesssim 0.12\ {\rm Z}_\odot$) powered by massive stars. Further, we detect a metal-enriched damped Ly$α$ system (DLA) associated with the galaxy with the HI column density of $N_{\rm{HI}}\approx 10^{21.8}$ cm$^{-2}$. The metallicity of the associated DLA may reach the super solar metallicity (${\gtrsim Z}_\odot$). Moreover, thanks to JWST and gravitational lensing, we present the resolved UV slope ($β$) map at the spatial resolution of $\approx 100$ pc at $z=4$, with steep UV slopes reaching $β\approx -2.5$ around three star-forming clumps. Combining with low-redshift analogs, our observations suggest that low-mass, low-metallicity galaxies, which dominate reionization, could be surrounded by a high covering fraction of the metal-enriched, neutral-gaseous clouds. This implies that the metal enrichment of low-mass galaxies is highly efficient, and further support that in low-mass galaxies, only a small fraction of ionizing radiation can escape through the interstellar or circumgalactic channels with low column-density neutral gas. △ Less

Submitted 1 December, 2022; v1 submitted 7 September, 2022; originally announced September 2022.

Comments: 4 figures, 1 table; Accepted for Publication in the ApJL

SDSS-IV MaNGA: pyPipe3D analysis release for 10,000 galaxies

Authors: S. F. Sánchez, J. K. Barrera-Ballesteros, E. Lacerda, A. Mejıa-Narvaez, A. Camps-Fariña, G. Bruzual, C. Espinosa-Ponce, A. Rodrıguez-Puebla, A. R. Calette, H. Ibarra-Medel, V. Avila-Reese, H. Hernandez-Toledo, M. A. Bershady, M. Cano-Diaz, A. M. Munguia-Cordova

Abstract: We present here the analysis performed using the pyPipe3D pipeline for the final MaNGA dataset included in the SDSS seventeenth data-release. This dataset comprises more than 10,000 individual datacubes, being the integral field spectroscopy galaxy survey with the largest number of galaxies. pyPipe3D processes the IFS datacubes to extract spatially-resolved spectroscopic properties of both the ste… ▽ More We present here the analysis performed using the pyPipe3D pipeline for the final MaNGA dataset included in the SDSS seventeenth data-release. This dataset comprises more than 10,000 individual datacubes, being the integral field spectroscopy galaxy survey with the largest number of galaxies. pyPipe3D processes the IFS datacubes to extract spatially-resolved spectroscopic properties of both the stellar population and the ionized-gas emission lines. A brief summary of the properties of the sample and the characteristics of the analyzed data are included. The article provides details on (i) the performed analysis, (ii) a description of the pipeline, (iii) the adopted stellar population library, (iv) the morphological and photometric analysis, (v) the adopted datamodel for the derived spatially resolved properties and (vi) the individual integrated and characteristic galaxy properties included in a final catalog. Comparisons with results from a previous version of the pipeline for earlier data releases and from other tools using this dataset are included. A practical example on how to use of the full dataset, and final catalog illustrates how to handle the delivered product. Our full analysis can be accessed and downloaded from the webpage http://ifs.astroscu.unam.mx/MaNGA/Pipe3D_v3_1_1/. △ Less

Submitted 14 June, 2022; originally announced June 2022.

Comments: 66 pages, 29 figures, 18 Tables, submitted to ApJS

SDSS-IV MaNGA: The radial distribution of physical properties within galaxies in the nearby universe

Authors: J. K. Barrera-Ballesteros, S. F. Sánchez, C. Espinosa-Ponce, C. López-Cobá, L. Carigi, A. Z. Lugo-Aranda, E. Lacerda, G. Bruzual, H. Hernandez-Toledo, N. Boardman, N. Drory, Richard R. Lane, J. R. Brownstein

Abstract: Using the largest sample of galaxies observed with an optical integral field unit (IFU, the SDSS-IV MaNGA survey, $\sim$10000 targets), we derive the radial distribution of the physical properties obtained from the stellar continuum and the ionized-gas emission lines. Given the large sample, we are able to explore the impact of the total stellar mass and morphology by averaging those radial distri… ▽ More Using the largest sample of galaxies observed with an optical integral field unit (IFU, the SDSS-IV MaNGA survey, $\sim$10000 targets), we derive the radial distribution of the physical properties obtained from the stellar continuum and the ionized-gas emission lines. Given the large sample, we are able to explore the impact of the total stellar mass and morphology by averaging those radial distributions for different bins of both global properties. We use a piece-wise analysis to characterize the slopes of the gradients from those properties at different galactocentric distances. In general we find that most of the properties -- derived from both the stellar continuum and the ionized gas emission lines -- exhibit a negative gradient with a secondary impact by global properties such as the total stellar mass or morphology. Our results confirm the intimate interplay between the properties of the stellar component and those of the ionized gas at local (kpc) scales in order to set the observed gradients. Furthermore, the resemblance of the gradients for similar global properties (in particular for the stellar parameters) indicates statistical similar histories of star formation and chemical enrichment with an initial radial gas distribution following the potential of the galaxy. △ Less

Submitted 14 June, 2022; originally announced June 2022.

Comments: 58 pages. 59 Figures. Submitted to the RMxAA

Post-starburst galaxies in the centers of intermediate redshift clusters

Authors: Ariel Werle, Bianca Poggianti, Alessia Moretti, Callum Bellhouse, Benedetta Vulcani, Marco Gullieuszik, Mario Radovich, Jacopo Fritz, Alessandro Ignesti, Johan Richard, Geneviève Soucail, Gustavo Bruzual, Stephane Charlot, Matilde Mingozzi, Cecilia Bacchini, Neven Tomicic, Rory Smith, Andrea Kulier, Giorgia Peluso, Andrea Franchetto

Abstract: We present results from MUSE spatially-resolved spectroscopy of 21 post-starburst galaxies in the centers of 8 clusters from $z\sim0.3$ to $z\sim0.4$. We measure spatially resolved star-formation histories (SFHs), the time since quenching ($t_Q$) and the fraction of stellar mass assembled in the past 1.5 Gyr ($μ_{1.5}$). The SFHs display a clear enhancement of star-formation prior to quenching for… ▽ More We present results from MUSE spatially-resolved spectroscopy of 21 post-starburst galaxies in the centers of 8 clusters from $z\sim0.3$ to $z\sim0.4$. We measure spatially resolved star-formation histories (SFHs), the time since quenching ($t_Q$) and the fraction of stellar mass assembled in the past 1.5 Gyr ($μ_{1.5}$). The SFHs display a clear enhancement of star-formation prior to quenching for 16 out of 21 objects, with at least 10% (and up to $>50$%) of the stellar mass being assembled in the past 1.5 Gyr and $t_Q$ ranging from less than 100 Myrs to $\sim800$ Myrs. By mapping $t_Q$ and $μ_{1.5}$, we analyze the quenching patterns of the galaxies. Most galaxies in our sample have quenched their star-formation from the outside-in or show a side-to-side/irregular pattern, both consistent with quenching by ram-pressure stripping. Only three objects show an inside-out quenching pattern, all of which are at the high-mass end of our sample. At least two of them currently host an active galactic nucleus. In two post-starbursts, we identify tails of ionized gas indicating that these objects had their gas stripped by ram pressure very recently. Post-starburst features are also found in the stripped regions of galaxies undergoing ram-pressure stripping in the same clusters, confirming the link between these classes of objects. Our results point to ram-pressure stripping as the main driver of fast quenching in these environments, with active galactic nuclei playing a role at high stellar masses. △ Less

Submitted 16 March, 2022; originally announced March 2022.

Comments: Accepted for publication in ApJ

Chemical evolution history of MaNGA galaxies

Authors: Artemi Camps-Fariña, Sebastián F. Sánchez, Alfredo Mejía-Narváez, Eduardo A. D. Lacerda, Leticia Carigi, Gustavo Bruzual, Paola Alvarez-Hurtado, Niv Drory, Richard R. Lane, Nicholas Fraser Boardman, Guillermo A. Blanc

Abstract: We show the results of a study using the spectral synthesis technique study for the full MaNGA sample showing their Chemical Enrichment History (ChEH) as well as the evolution of the stellar mass-metallicity relation (MZR) over cosmic time. We find that the more massive galaxies became enriched first and the lower mass galaxies did so later, producing a change in the MZR which becomes shallower in… ▽ More We show the results of a study using the spectral synthesis technique study for the full MaNGA sample showing their Chemical Enrichment History (ChEH) as well as the evolution of the stellar mass-metallicity relation (MZR) over cosmic time. We find that the more massive galaxies became enriched first and the lower mass galaxies did so later, producing a change in the MZR which becomes shallower in time. Separating the sample into morphology and star-forming status bins some particularly interesting results appear: The mass dependency of the MZR becomes less relevant for later morphological types, to the extent that it inverts for Sd/Irr galaxies, suggesting that morphology is at least as important a factor as mass in chemical evolution. The MZR for the full sample shows a flattening at the high-mass end and another at the low-mass range, but the former only appears for retired galaxies while the latter only appears for star-forming galaxies. We also find that the average metallicity gradient is currently negative for all mass bins but for low mass galaxies it was inverted at some point in the past, before which all galaxies had a positive gradient. We also compare how diverse the ChEHs are in the different bins considered as well as what primarily drives the diversity: How much galaxies become enriched or how quickly they do so. △ Less

Submitted 2 March, 2022; originally announced March 2022.

Comments: Currently in review at ApJ

Direct Constraints on the Extremely Metal-Poor Massive Stars Underlying Nebular C IV Emission from Ultra-Deep HST/COS Ultraviolet Spectroscopy

Authors: Peter Senchyna, Daniel P. Stark, Stephane Charlot, Adele Plat, Jacopo Chevallard, Zuyi Chen, Tucker Jones, Ryan L. Sanders, Gwen C. Rudie, Thomas J. Cooper, Gustavo Bruzual

Abstract: Metal-poor nearby galaxies hosting massive stars have a fundamental role to play in our understanding of both high-redshift galaxies and low metallicity stellar populations. But while much attention has been focused on their bright nebular gas emission, the massive stars that power it remain challenging to constrain. Here we present exceptionally deep Hubble Space Telescope ultraviolet spectra tar… ▽ More Metal-poor nearby galaxies hosting massive stars have a fundamental role to play in our understanding of both high-redshift galaxies and low metallicity stellar populations. But while much attention has been focused on their bright nebular gas emission, the massive stars that power it remain challenging to constrain. Here we present exceptionally deep Hubble Space Telescope ultraviolet spectra targeting six galaxies that power strong nebular C IV emission approaching that encountered at $z>6$. We find that the strength and spectral profile of the nebular C IV in these new spectra follow a sequence evocative of resonant scattering models, indicating that the hot circumgalactic medium likely plays a key role in regulating C IV escape locally. We constrain the metallicity of the massive stars in each galaxy by fitting the forest of photospheric absorption lines, reporting measurements driven by iron that lie uniformly below 10% solar. Comparison with the gas-phase oxygen abundances reveals evidence for enhancement in O/Fe above solar across the sample, robust to assumptions about the absolute gas-phase metallicity scale. This supports the idea that these local systems are more chemically-similar to their primordial high-redshift counterparts than to the bulk of nearby galaxies. Finally, we find significant tension between the strong stellar wind profiles observed and our population synthesis models constrained by the photospheric forest in our highest-quality spectra. This reinforces the need for caution in interpreting wind lines in isolation at high-redshift, but also suggests a unique path towards validating fundamental massive star physics at extremely low metallicity with integrated ultraviolet spectra. △ Less

Submitted 22 November, 2021; originally announced November 2021.

Comments: 34 pages, 14 figures, submitted to AAS Journals

Synthetic photometry of OB star clusters with stochastically sampled IMFs: analysis of models and HST observations

Authors: Rogelio Orozco-Duarte, Aida Wofford, Alba Vidal-García, Gustavo Bruzual, Stephane Charlot, Mark R. Krumholz, Stephen Hannon, Janice Lee, Timothy Wofford, Michele Fumagalli, Daniel Dale, Matteo Messa, Eva K. Grebel, Linda Smith, Kathryn Grasha, David Cook

Abstract: We present a pilot library of synthetic NUV, U, B, V, and I photometry of star clusters with stochastically sampled IMFs and ionized gas for initial masses, $M_i=10^3$, $10^4$, and $10^5$ M$_{\odot}$; $t=1$, 3, 4, and 8 Myr; $Z=0.014$ and $Z=0.002$; and log(U$_{\rm S}$) =-2 and -3. We compare the library with predictions from deterministic models and observations of isolated low-mass ($<10^4$ M… ▽ More We present a pilot library of synthetic NUV, U, B, V, and I photometry of star clusters with stochastically sampled IMFs and ionized gas for initial masses, $M_i=10^3$, $10^4$, and $10^5$ M$_{\odot}$; $t=1$, 3, 4, and 8 Myr; $Z=0.014$ and $Z=0.002$; and log(U$_{\rm S}$) =-2 and -3. We compare the library with predictions from deterministic models and observations of isolated low-mass ($<10^4$ M$_{\odot}$) star clusters with co-spatial compact H\2 regions. The clusters are located in NGC 7793, one of the nearest galaxies observed as part of the \hst LEGUS and H$α$-LEGUS surveys. 1) For model magnitudes that only account for the stars: a) the residual |deterministic mag - median stochastic mag| can be $\ge0.5$ mag, even for $M_i=10^5$ M$_{\odot}$; and b) the largest spread in stochastic magnitudes occurs when Wolf-Rayet stars are present. 2) For $M_i=10^5$ M$_{\odot}$: a) the median stochastic mag with gas can be $>$1.0 mag more luminous than the median stochastic magnitude without gas; and b) nebular emission lines can contribute with $>50\%$ and $>30\%$ to the total emission in the V and I bands, respectively. 3) Age-dating OB-star clusters via deterministic tracks in the U-B vs. V-I plane is highly uncertain at $Z=0.014$ for $M_i\sim10^3$ M$_{\odot}$ and $Z=0.002$ for $M_i\sim10^3-10^5$ M$_{\odot}$. 4) For low-mass clusters, the V-band extinction derived with stochastic models significantly depends on the value of log(U$_{\rm S}$). 5) The youngest clusters tend to have higher extinction. 6) The majority of clusters have multi-peaked age PDFs. 7) Finally, we discuss the importance of characterising the true variance in the number of stars per mass bin in nature. △ Less

Submitted 11 October, 2021; originally announced October 2021.

J-PLUS: Detecting and studying extragalactic globular clusters -- the case of NGC 1023

Authors: Danielle de Brito Silva, Paula Coelho, Arianna Cortesi, Gustavo Bruzual, Gladis Magris C., Ana L. Chies-Santos, Jose A. Hernandez-Jimenez, Alessandro Ederoclite, Izaskun San Roman, Jesús Varela, Duncan A. Forbes, Yolanda Jiménez-Teja, Javier Cenarro, David Cristóbal-Hornillos, Carlos Hernández-Monteagudo, Carlos López-Sanjuan, Antonio Marín-Franch, Mariano Moles, Héctor Vázquez Ramió, Renato Dupke, Laerte Sodré Jr. 2, Raul E. Angulo

Abstract: Extragalactic globular clusters (GCs) are key objects for studying the history of galaxies. The arrival of wide-field surveys such as the Javalambre Photometric Local Universe Survey (J-PLUS) offers new possibilities for the study of these systems. We perform the first study of GCs in J-PLUS to recover information about the history of NGC 1023 taking advantage of wide-field images and 12 filters.… ▽ More Extragalactic globular clusters (GCs) are key objects for studying the history of galaxies. The arrival of wide-field surveys such as the Javalambre Photometric Local Universe Survey (J-PLUS) offers new possibilities for the study of these systems. We perform the first study of GCs in J-PLUS to recover information about the history of NGC 1023 taking advantage of wide-field images and 12 filters. We develop the semiautomatic pipeline GCFinder that detects GC candidates in J-PLUS images and can also be adapted to similar surveys. We study the stellar population properties of a sub-sample of GC candidates using spectral energy distribution (SED) fitting. We find 523 GC candidates in NGC 1023, of which $\sim$300 are new. We identify subpopulations of GC candidates, where age and metallicity distributions have multiple peaks. By comparing our results with simulations, we report a possible broad age-metallicity relation, evidence that NGC 1023 experienced accretion events in the past. The dominating age peak is at $10^{10}$ yr. We report a correlation between masses and ages that suggests that massive GC candidates are more likely to survive the turbulent history of the host galaxy. Modeling the light of NGC 1023, we find two spiral-like arms and detect a displacement of the galaxy's photometric center with respect to the outer isophotes and center of GC distribution ($\sim$700 pc and $\sim$1600 pc, respectively), which could be the result of ongoing interaction between NGC 1023 and NGC 1023A. By studying the GC system of NGC 1023 with J-PLUS we showcase the power of multi-band surveys for this kind of study and find evidence of a complex accretion history of the host galaxy. △ Less

Submitted 18 October, 2023; v1 submitted 8 October, 2021; originally announced October 2021.

Comments: Published in A&A, 2022

Journal ref: A&A 664, A129 (2022)

CoSHA: Code for Stellar properties Heuristic Assignment -- for the MaStar stellar library

Authors: Alfredo Mejía-Narváez, Gustavo Bruzual, Sebastian F. Sánchez, Leticia Carigi, Jorge Barrera-Ballesteros, Mabel Valerdi, Renbin Yan, Niv Drory

Abstract: We introduce \cosha{}: a Code for Stellar properties Heuristic Assignment. In order to estimate the stellar properties, \cosha{} implements a Gradient Tree Boosting algorithm to label each star across the parameter space ($T_\mathrm{eff}$, $\log{g}$, $[\mathrm{Fe}/\mathrm{H}]$, and $[α/\mathrm{Fe}]$). We use \cosha{} to estimate these stellar atmospheric parameters of $22\,$k unique stars in the M… ▽ More We introduce \cosha{}: a Code for Stellar properties Heuristic Assignment. In order to estimate the stellar properties, \cosha{} implements a Gradient Tree Boosting algorithm to label each star across the parameter space ($T_\mathrm{eff}$, $\log{g}$, $[\mathrm{Fe}/\mathrm{H}]$, and $[α/\mathrm{Fe}]$). We use \cosha{} to estimate these stellar atmospheric parameters of $22\,$k unique stars in the MaNGA Stellar Library (MaStar). To quantify the reliability of our approach, we run both internal tests using the Göttingen Stellar Library (GSL, a theoretical library) and the first data release of MaStar, and external tests by comparing the resulting distributions in the parameter space with the APOGEE estimates of the same properties. In summary, our parameter estimates span in the ranges: $T_\mathrm{eff}=[2900,12000]\,$K, $\log{g}=[-0.5,5.6]$, $[\mathrm{Fe}/\mathrm{H}]=[-3.74,0.81]$, $[α/\mathrm{Fe}]=[-0.22,1.17]$. {We report internal (external) uncertainties of the properties of $σ_{T_\mathrm{eff}}\sim43\,(240)\,$K, $σ_{\log{g}}\sim0.2\,(0.4)$, $σ_{[\mathrm{Fe}/\mathrm{H}]}\sim0.16\,(0.24)$, $σ_{[α/\mathrm{Fe}]}\sim0.09\,(0.08)$.} These uncertainties are comparable to those of other methods with similar objectives. Despite the fact that \cosha{} is not aware of the spatial distribution of these physical properties in the Milky Way, we are able to recover the main trends known in the literature. The catalog of physical properties for MaStar can be accessed in \url{http://ifs.astroscu.unam.mx/MaStar}. △ Less

Submitted 27 July, 2022; v1 submitted 3 August, 2021; originally announced August 2021.

Comments: 25 pages, 12 figures, published in ApJS

Constraints on the dust extinction law of the Galaxy with Swift/UVOT, Gaia and 2MASS

Authors: I. Ferreras, M. Tress, G. Bruzual, S. Charlot, M. Page, V. Yershov, P. Kuin, D. Kawata, M. Cropper

Abstract: We explore variations of the dust extinction law of the Milky Way by selecting stars from the Swift/UVOT Serendipitous Source Catalogue, cross-matched with Gaia DR2 and 2MASS to produce a sample of 10,452 stars out to ~4kpc with photometry covering a wide spectral window. The near ultraviolet passbands optimally encompass the 2175A bump, so that we can simultaneously fit the net extinction, quoted… ▽ More We explore variations of the dust extinction law of the Milky Way by selecting stars from the Swift/UVOT Serendipitous Source Catalogue, cross-matched with Gaia DR2 and 2MASS to produce a sample of 10,452 stars out to ~4kpc with photometry covering a wide spectral window. The near ultraviolet passbands optimally encompass the 2175A bump, so that we can simultaneously fit the net extinction, quoted in the V band (A$_V$), the steepness of the wavelength dependence ($δ$) and the bump strength (E$_b$). The methodology compares the observed magnitudes with theoretical stellar atmospheres from the models of Coelho. Significant correlations are found between these parameters, related to variations in dust composition, that are complementary to similar scaling relations found in the more complex dust attenuation law of galaxies - that also depend on the distribution of dust among the stellar populations within the galaxy. We recover the strong anticorrelation between A$_V$ and Galactic latitude, as well as a weaker bump strength at higher extinction. $δ$ is also found to correlate with latitude, with steeper laws towards the Galactic plane. Our results suggest that variations in the attenuation law of galaxies cannot be fully explained by dust geometry. △ Less

Submitted 5 May, 2021; originally announced May 2021.

Comments: 12 pages, 11 figures, 2 tables. Accepted for publication in MNRAS

Star cluster survival in dark matter halos: An old cluster in Eridanus II?

Authors: Jairo A. Alzate, Verónica Lora, Gustavo Bruzual, Luis Lomelí-Núñez, Bernardo Cervantes Sodi

Abstract: The star formation history and the internal dynamics of Milky Way satellite galaxies are often complicated. In the last years, a substantial fraction of the known faint dwarf satellites have been studied. Some of them show embedded stellar substructures, such as star clusters and even globular star clusters. In this work we study Eridanus II, a dwarf spheroidal satellite which hosts a star cluster… ▽ More The star formation history and the internal dynamics of Milky Way satellite galaxies are often complicated. In the last years, a substantial fraction of the known faint dwarf satellites have been studied. Some of them show embedded stellar substructures, such as star clusters and even globular star clusters. In this work we study Eridanus II, a dwarf spheroidal satellite which hosts a star cluster, using published and archival data from the Hubble Space Telescope Advanced Camera for Surveys. We employ a Bayesian hierarchical method to infer the star formation history of Eridanus II. We find that the bulk of the stars in Eridanus II are very old ($13.5_{-1}^{+0.5}$ Gyr) and quite metal poor ($Z$\,=\,$0.00001$). We do not find any evidence of the presence of an intermediate age or young population in Eri II. We cannot date the embedded star cluster as a separate entity, but we find it likely that the cluster has a similar age and metallicity as the bulk of the stars in Eri II. The existence of an old star cluster in a dark matter dominated old metal poor dwarf galaxy is of major importance to cast light on the dark matter distribution within dwarf galaxies. The existence of intermediate age stars is required by the recent detection of carbon stars in Eri II. Since no recent star formation is detected, {\it blue-straggler} fusions of lower mass stars are the most likely origin of the carbon star progenitors. △ Less

Submitted 4 May, 2021; originally announced May 2021.

Comments: 13 pages, 13 figures, 5 tables. To appear in MNRAS: Accepted 2021 May 3. Received 2021 April 30; in original form 2020 July 28

Modelling H$_2$ and its effects on star formation using a joint implementation of GADGET-3 and KROME

Authors: Emanuel Sillero, Patricia B. Tissera, Diego G. Lambas, Stefano Bovino, Dominik R. Schleicher, Tommaso Grassi, Gustavo Bruzual, Stéphane Charlot

Abstract: We present P-GADGET3-K, an updated version of GADGET3, that incorporates the chemistry package KROME. P-GADGET3-K follows the hydrodynamical and chemical evolution of cosmic structures, incorporating the chemistry and cooling of H$_2$ and metal cooling in non-equilibrium. We performed different runs of the same ICs to assess the impact of various physical parameters and prescriptions, namely gas m… ▽ More We present P-GADGET3-K, an updated version of GADGET3, that incorporates the chemistry package KROME. P-GADGET3-K follows the hydrodynamical and chemical evolution of cosmic structures, incorporating the chemistry and cooling of H$_2$ and metal cooling in non-equilibrium. We performed different runs of the same ICs to assess the impact of various physical parameters and prescriptions, namely gas metallicity, molecular hydrogen formation on dust, star formation recipes including or not H$_2$ dependence, and the effects of numerical resolution. We find that the characteristics of the simulated systems, both globally and at kpc-scales, are in good agreement with several observable properties of molecular gas in star-forming galaxies. The surface density profiles of SFR and H$_2$ are found to vary with the clumping factor and resolution. In agreement with previous results, the chemical enrichment of the gas component is found to be a key ingredient to model the formation and distribution of H$_2$ as a function of gas density and temperature. A SF algorithm that takes into account the H$_2$ fraction together with a treatment for the local stellar radiation field improves the agreement with observed H$_2$ abundances over a wide range of gas densities and with the molecular Kennicutt-Schmidt law, implying a more realistic modelling of the star formation process. △ Less

Submitted 9 April, 2021; originally announced April 2021.

Comments: 21 pages, 16 figures. Accepted for publication in MNRAS

The miniJPAS survey: Identification and characterization of galaxy populations with the J-PAS photometric system

Authors: R. M. González Delgado, L. A. Díaz-García, A. de Amorim, G. Bruzual, R. Cid Fernandes, E. Pérez, S. Bonoli, A. J. Cenarro, P. R. T. Coelho, A. Cortesi, R. García-Benito, R. López Fernández, G. Martínez-Solaeche, J. E. Rodríguez-Martín, G. Magris, A. Mejía-Narvaez, D. Brito-Silva, L. R. Abramo, J. M. Diego, R. A. Dupke, A. Hernán-Caballero, C. Hernández-Monteagudo, C. López-Sanjuan, A. Marín-Franch, V. Marra , et al. (17 additional authors not shown)

Abstract: J-PAS will soon start imaging 8000 deg2 of the northern sky with its unique set of 56 filters (R $\sim$ 60). Before, we observed 1 deg2 on the AEGIS field with an interim camera with all the J-PAS filters. With this data (miniJPAS), we aim at proving the scientific potential of J-PAS to identify and characterize the galaxy populations with the goal of performing galaxy evolution studies across cos… ▽ More J-PAS will soon start imaging 8000 deg2 of the northern sky with its unique set of 56 filters (R $\sim$ 60). Before, we observed 1 deg2 on the AEGIS field with an interim camera with all the J-PAS filters. With this data (miniJPAS), we aim at proving the scientific potential of J-PAS to identify and characterize the galaxy populations with the goal of performing galaxy evolution studies across cosmic time. Several SED-fitting codes are used to constrain the stellar population properties of a complete flux-limited sample (rSDSS <= 22.5 AB) of miniJPAS galaxies that extends up to z = 1. We find consistent results on the galaxy properties derived from the different codes, independently of the galaxy spectral-type or redshift. For galaxies with SNR>=10, we estimate that the J-PAS photometric system allows to derive stellar population properties with a precision that is equivalent to that obtained with spectroscopic surveys of similar SNR. By using the dust-corrected (u-r) colour-mass diagram, a powerful proxy to characterize galaxy populations, we find that the fraction of red and blue galaxies evolves with cosmic time, with red galaxies being $\sim$ 38% and $\sim$ 18% of the whole population at z = 0.1 and z = 0.5, respectively. At all redshifts, the more massive galaxies belong to the red sequence and these galaxies are typically older and more metal rich than their counterparts in the blue cloud. Our results confirm that with J-PAS data we will be able to analyze large samples of galaxies up to z $\sim$ 1, with galaxy stellar masses above of log(M$_*$/M$_{\odot}$) $\sim$ 8.9, 9.5, and 9.9 at z = 0.3, 0.5, and 0.7, respectively. The SFH of a complete sub-sample of galaxies selected at z $\sim$ 0.1 with log(M$_*$/M$_{\odot}$) > 8.3 constrain the cosmic evolution of the star formation rate density up to z $\sim$ 3 in good agreement with results from cosmological surveys. △ Less

Submitted 5 March, 2021; v1 submitted 25 February, 2021; originally announced February 2021.

Comments: 27 pages, 23 figures, 3 tables, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 649, A79 (2021)

On the Thermally Pulsing Asymptotic Giant Branch Contribution to the Light of Nearby Disk Galaxies

Authors: Eric E. Martínez-García, Gustavo Bruzual, Rosa A. González-Lópezlira, Lino H. Rodríguez-Merino

Abstract: The study of the luminosity contribution from thermally pulsing asymptotic giant branch (TP-AGB) stars to the stellar populations of galaxies is crucial to determine their physical parameters (e.g., stellar mass and age). We use a sample of 84 nearby disk galaxies to explore diverse stellar population synthesis models with different luminosity contributions from TP-AGB stars. We fit the models to… ▽ More The study of the luminosity contribution from thermally pulsing asymptotic giant branch (TP-AGB) stars to the stellar populations of galaxies is crucial to determine their physical parameters (e.g., stellar mass and age). We use a sample of 84 nearby disk galaxies to explore diverse stellar population synthesis models with different luminosity contributions from TP-AGB stars. We fit the models to optical and near-infrared (NIR) photometry, on a pixel-bypixel basis. The statistics of the fits show a preference for a low-luminosity contribution (i.e., high mass-to-light ratio in the NIR) from TP-AGB stars. Nevertheless, for 30 percent to 40 percent of the pixels in our sample a high-luminosity contribution (hence low mass-to-light ratio in the NIR) from TP-AGB stars is favored. According to our findings, the mean TP-AGB star luminosity contribution in nearby disk galaxies may vary with Hubble type. This may be a consequence of the variation of the TP-AGB mass-loss rate with metallicity, if metal-poor stars begin losing mass earlier than metal-rich stars, because of a pre-dust wind that precedes the dust-driven wind. △ Less

Submitted 17 February, 2021; v1 submitted 27 November, 2020; originally announced November 2020.

Comments: Published in 2021, ApJ, 908, 110

Star formation history of the solar neighbourhood as told by Gaia

Authors: Jairo A. Alzate, Gustavo Bruzual, Daniel J. Díaz-González

Abstract: The Gaia DR2 catalog is the best source of stellar astrometric and photometric data available today. The history of the Milky Way galaxy is written in stone in this data set. Parallaxes and photometry tell us where the stars are today, when were they formed, and with what chemical content, i.e. their star formation history (SFH). We develop a Bayesian hierarchical model suited to reconstruct the S… ▽ More The Gaia DR2 catalog is the best source of stellar astrometric and photometric data available today. The history of the Milky Way galaxy is written in stone in this data set. Parallaxes and photometry tell us where the stars are today, when were they formed, and with what chemical content, i.e. their star formation history (SFH). We develop a Bayesian hierarchical model suited to reconstruct the SFH of a resolved stellar population. We study the stars brighter than $G\,=\,15$ within 100 pc of the Sun in Gaia DR2 and derive a SFH of the solar neighbourhood in agreement with previous determinations and improving upon them because we detect chemical enrichment. Our results show a maximum of star formation activity about 10 Gyr ago, producing large numbers of stars with slightly below solar metallicity (Z=0.014), followed by a decrease in star formation up to a minimum level occurring around 8 Gyr ago. After a quiet period, star formation rises to a maximum at about 5 Gyr ago, forming stars of solar metallicity (Z=0.017). Finally, star formation has been decreasing until the present, forming stars of Z=0.03 at a residual level. We test the effects introduced in the inferred SFH by ignoring the presence of unresolved binary stars in the sample, reducing the apparent limiting magnitude, and modifying the stellar initial mass function. △ Less

Submitted 11 November, 2020; originally announced November 2020.

Comments: 33 pages, 28 figures

A quantitative demonstration that stellar feedback locally regulates galaxy growth

Authors: Javier Zaragoza-Cardiel, Jacopo Fritz, Itziar Aretxaga, Yalia D. Mayya, Daniel Rosa-González, John E. Beckman, Gustavo Bruzual, Stephane Charlot

Abstract: We have applied stellar population synthesis to 500 pc sized regions in a sample of 102 galaxy discs observed with the MUSE spectrograph. We derived the star formation history and analyse specifically the "recent" ($20\rm{Myr}$) and "past" ($570\rm{Myr}$) age bins. Using a star formation self-regulator model we can derive local mass-loading factors, $η$ for specific regions, and find that this fac… ▽ More We have applied stellar population synthesis to 500 pc sized regions in a sample of 102 galaxy discs observed with the MUSE spectrograph. We derived the star formation history and analyse specifically the "recent" ($20\rm{Myr}$) and "past" ($570\rm{Myr}$) age bins. Using a star formation self-regulator model we can derive local mass-loading factors, $η$ for specific regions, and find that this factor depends on the local stellar mass surface density, $Σ_*$, in agreement with the predictions form hydrodynamical simulations including supernova feedback. We integrate the local $η$-$Σ_*$ relation using the stellar mass surface density profiles from the Spitzer Survey of Stellar Structure in Galaxies (S4G) to derive global mass-loading factors, $η_{\rm{G}}$, as a function of stellar mass, $M_*$. The $η_{\rm{G}}$-$M_*$ relation found is in very good agreement with hydrodynamical cosmological zoom-in galaxy simulations. The method developed here offers a powerful way of testing different implementations of stellar feedback, to check on how realistic are their predictions. △ Less

Submitted 18 September, 2020; originally announced September 2020.

Comments: Accepted for publication in MNRAS

Ultraviolet spectra of extreme nearby star-forming regions: evidence for an overabundance of very massive stars

Authors: Peter Senchyna, Daniel P. Stark, Stéphane Charlot, Jacopo Chevallard, Gustavo Bruzual, Alba Vidal-García

Abstract: As deep spectroscopic campaigns extend to higher redshifts and lower stellar masses, the interpretation of galaxy spectra depends increasingly upon models for very young stellar populations. Here we present new HST/COS ultraviolet spectroscopy of seven nearby ($<120$ Mpc) star-forming regions hosting very young stellar populations ($\sim$ 4-20 Myr) with optical Wolf-Rayet stellar wind signatures,… ▽ More As deep spectroscopic campaigns extend to higher redshifts and lower stellar masses, the interpretation of galaxy spectra depends increasingly upon models for very young stellar populations. Here we present new HST/COS ultraviolet spectroscopy of seven nearby ($<120$ Mpc) star-forming regions hosting very young stellar populations ($\sim$ 4-20 Myr) with optical Wolf-Rayet stellar wind signatures, ideal laboratories in which to test these stellar models. We detect nebular C III] in all seven, but at equivalent widths uniformly $< 10$ Å. This suggests that even for very young stellar populations, the highest equivalent width C III] emission at $\geq 15$ Å is reserved for inefficiently-cooled gas at metallicities at or below that of the SMC. The spectra also reveal strong C IV P-Cygni profiles and broad He II emission formed in the winds of massive stars, including some of the most prominent He II stellar wind lines ever detected in integrated spectra. We find that the latest stellar population synthesis prescriptions with improved treatment of massive stars nearly reproduce the entire range of stellar He II wind strengths observed here. However, we find that these models cannot simultaneously match the strongest wind features alongside the optical nebular line constraints. This discrepancy can be naturally explained by an overabundance of very massive stars produced by a high incidence of binary mass transfer and mergers occurring on short $\lesssim 10$ Myr timescales, suggesting these processes may be crucial for understanding the highest-sSFR galaxies in the early Universe. Reproducing both the stellar and nebular light of young systems such as these will be a crucial benchmark for the next generation of stellar population synthesis models. △ Less

Submitted 22 August, 2020; originally announced August 2020.

Comments: 25 pages, 9 figures; submitted to MNRAS

MEGARA-IFU detection of extended HeII4686 nebular emission in the central region of NGC1569 and its ionization budget

Authors: Y. D. Mayya, E. Carrasco, V. M. A. Gomez-Gonzalez, J. Zaragoza-Cardiel, A. Gil de Paz, P. A. Ovando, M. Sanchez-Cruces, L. Lomeli-Nunez, L. Rodriguez-Merino, D. Rosa-Gonzalez, S. Silich, G. Tenorio-Tagle, G. Bruzual, S. Charlot, R. Terlevich, E. Terlevich, O. Vega, J. Gallego, J. Iglesias-Paramo, A. Castillo-Morales, M. L. Garcia-Vargas, P. Gomez-Alvarez, S. Pascual, A. Perez-Calpena

Abstract: We here report the detection of extended HeII4686 nebular emission in the central region of NGC1569 using the integral field spectrograph MEGARA at the 10.4-m Gran Telescopio Canarias. The observations cover a Field of View (FoV) of 12.5 arcsec x 11.3 arcsec at seeing-limited spatial resolution of ~15 pc and at a spectral resolution of R=6000 in the wavelength range 4330--5200 Angstrom. The emissi… ▽ More We here report the detection of extended HeII4686 nebular emission in the central region of NGC1569 using the integral field spectrograph MEGARA at the 10.4-m Gran Telescopio Canarias. The observations cover a Field of View (FoV) of 12.5 arcsec x 11.3 arcsec at seeing-limited spatial resolution of ~15 pc and at a spectral resolution of R=6000 in the wavelength range 4330--5200 Angstrom. The emission extends over a semi-circular arc of ~40 pc width and ~150 pc diameter around the super star cluster A (SSC-A). The Av derived using Balmer decrement varies from the Galactic value of 1.6 mag to a maximum of ~4.5 mag, with a mean value of 2.65+/-0.60 mag. We infer 124+/-11 Wolf-Rayet (WR) stars in SSC-A using the HeII4686 broad feature and Av=2.3 mag. The He+ ionizing photon rate from these WR stars is sufficient to explain the luminosity of the HeII4686 nebula. The observationally-determined total He+ and H0 ionizing photon rates, their ratio, and the observed number of WR stars in SSC-A are all consistent with the predictions of simple stellar population models at an age of 4.0+/-0.5 Myr, and mass of (5.5+/-0.5)x10^5 Msun. Our observations reinforce the absence of WR stars in SSC-B, the second most massive cluster in the FoV. None of the other locations in our FoV where HeII4686 emission has been reported from narrow-band imaging observations contain WR stars. △ Less

Submitted 1 August, 2020; originally announced August 2020.

Comments: 20 pages. Accepted in MNRAS

The miniJPAS survey: a preview of the Universe in 56 colours

Authors: S. Bonoli, A. Marín-Franch, J. Varela, H. Vázquez Ramió, L. R. Abramo, A. J. Cenarro, R. A. Dupke, J. M. Vílchez, D. Cristóbal-Hornillos, R. M. González Delgado, C. Hernández-Monteagudo, C. López-Sanjuan, D. J. Muniesa, T. Civera, A. Ederoclite, A. Hernán-Caballero, V. Marra, P. O. Baqui, A. Cortesi, E. S. Cypriano, S. Daflon, A. L. de Amorim, L. A. Díaz-García, J. M. Diego, G. Martínez-Solaeche , et al. (144 additional authors not shown)

Abstract: The Javalambre-Physics of the Accelerating Universe Astrophysical Survey (J-PAS) will soon start to scan thousands of square degrees of the northern extragalactic sky with a unique set of $56$ optical filters from a dedicated $2.55$m telescope, JST, at the Javalambre Astrophysical Observatory. Before the arrival of the final instrument (a 1.2 Gpixels, 4.2deg$^2$ field-of-view camera), the JST was… ▽ More The Javalambre-Physics of the Accelerating Universe Astrophysical Survey (J-PAS) will soon start to scan thousands of square degrees of the northern extragalactic sky with a unique set of $56$ optical filters from a dedicated $2.55$m telescope, JST, at the Javalambre Astrophysical Observatory. Before the arrival of the final instrument (a 1.2 Gpixels, 4.2deg$^2$ field-of-view camera), the JST was equipped with an interim camera (JPAS-Pathfinder), composed of one CCD with a 0.3deg$^2$ field-of-view and resolution of 0.23 arcsec pixel$^{-1}$. To demonstrate the scientific potential of J-PAS, with the JPAS-Pathfinder camera we carried out a survey on the AEGIS field (along the Extended Groth Strip), dubbed miniJPAS. We observed a total of $\sim 1$ deg$^2$, with the $56$ J-PAS filters, which include $54$ narrow band (NB, $\rm{FWHM} \sim 145$Angstrom) and two broader filters extending to the UV and the near-infrared, complemented by the $u,g,r,i$ SDSS broad band (BB) filters. In this paper we present the miniJPAS data set, the details of the catalogues and data access, and illustrate the scientific potential of our multi-band data. The data surpass the target depths originally planned for J-PAS, reaching $\rm{mag}_{\rm {AB}}$ between $\sim 22$ and $23.5$ for the NB filters and up to $24$ for the BB filters ($5σ$ in a $3$~arcsec aperture). The miniJPAS primary catalogue contains more than $64,000$ sources extracted in the $r$ detection band with forced photometry in all other bands. We estimate the catalogue to be complete up to $r=23.6$ for point-like sources and up to $r=22.7$ for extended sources. Photometric redshifts reach subpercent precision for all sources up to $r=22.5$, and a precision of $\sim 0.3$% for about half of the sample. (Abridged) △ Less

Submitted 9 July, 2020; v1 submitted 3 July, 2020; originally announced July 2020.

Comments: The miniJPAS data and associated value added catalogues are publicly accessible via this url: http://archive.cefca.es/catalogues/minijpas-pdr201912

Journal ref: A&A 653, A31 (2021)

Clues on the history of early-type galaxies from SDSS spectra and GALEX photometry

Authors: Ariel Werle, Roberto Cid Fernandes, Natalia Vale Asari, Paula R. T. Coelho, Gustavo Bruzual, Stephane Charlot, Reinaldo R. de Carvalho, Fábio R. Herpich, Clauda Mendes de Oliveira, Laerte Sodré Jr., Daniel Ruschel Dutra, André de Amorim, Vitor M. Sampaio

Abstract: Stellar population studies of early-type galaxies (ETGs) based on their optical stellar continuum suggest that these are quiescent systems. However, emission lines and ultraviolet photometry reveal a diverse population. We use a new version of the STARLIGHT spectral synthesis code and state-of-the-art stellar population models to simultaneously fit SDSS spectra and GALEX photometry for a sample of… ▽ More Stellar population studies of early-type galaxies (ETGs) based on their optical stellar continuum suggest that these are quiescent systems. However, emission lines and ultraviolet photometry reveal a diverse population. We use a new version of the STARLIGHT spectral synthesis code and state-of-the-art stellar population models to simultaneously fit SDSS spectra and GALEX photometry for a sample of 3453 galaxies at $z < 0.1$ with $NUV-r > 5$ that are classified as elliptical by Galaxy Zoo. We reproduce $FUV$ magnitudes for 80 per cent of UV upturn galaxies selected using criteria from the literature, suggesting that additional stellar population ingredients such as binaries and extreme horizontal branch stars may have a limited contribution to the UV upturn. The addition of ultraviolet data leads to a broadening of the distributions of mean stellar ages, metallicities and attenuation. Stellar populations younger than $1\,$Gyr are required to reproduce the ultraviolet emission in 17 per cent of our sample. These systems represent 43 per cent of the sample at $5<NUV-r<5.5$ and span the same stellar mass range as other ETGs in our sample. ETGs with young stellar components have larger $Hα$ equivalent widths ($W_{Hα}$) and larger dust attenuation. Emission line ratios and $W_{Hα}$ indicate that the ionising source in these systems is a mixture of young and old stellar populations. Their young stellar populations are metal-poor, especially for high-mass galaxies, indicating recent star formation associated with rejuvenation events triggered by external processes, such as minor mergers. △ Less

Submitted 30 July, 2020; v1 submitted 2 July, 2020; originally announced July 2020.

Comments: Accepted for publication in MNRAS

Studying the ISM at ~10 pc scale in NGC 7793 with MUSE -- I. Data description and properties of the ionised gas

Authors: Lorenza Della Bruna, Angela Adamo, Arjan Bik, Michele Fumagalli, Rene Walterbos, Göran Östlin, Gustavo Bruzual, Daniela Calzetti, Stephane Charlot, Kathryn Grasha, Linda J. Smith, David Thilker, Aida Wofford

Abstract: Using MUSE AO data, we probe the ISM in the local spiral galaxy NGC 7793 at a spatial resolution of $\sim$ 10 pc. We identify HII regions and compile a catalogue of supernova remnants (SNRs), planetary nebulae (PNe) and Wolf Rayet stars (WR). We compute electron densities and temperatures from the [SII]6716/6731 and [SIII]6312/9069 line ratios. We study the properties of the ionised gas through BP… ▽ More Using MUSE AO data, we probe the ISM in the local spiral galaxy NGC 7793 at a spatial resolution of $\sim$ 10 pc. We identify HII regions and compile a catalogue of supernova remnants (SNRs), planetary nebulae (PNe) and Wolf Rayet stars (WR). We compute electron densities and temperatures from the [SII]6716/6731 and [SIII]6312/9069 line ratios. We study the properties of the ionised gas through BPT diagrams combined with gas velocity dispersion. We spectroscopically confirm 2 previously detected WR stars and a SNR and report the discovery of 7 WR, 1 SNR, and 2 PNe. The diffuse ionized gas (DIG) fraction is between $\sim$ 27 and 42% depending on the method used to define the HII region boundaries. In agreement with previous studies, we find that the DIG exhibits enhanced [SII]/H$α$ and [NII]/H$α$ ratios and a median temperature that is $\sim$ 3000 K higher than in HII regions. We also observe an apparent inverse correlation between temperature and H$α$ surface brightness. Overall, the observed [SII]6716/6731 ratio is consistent within 1$σ$ with $n_e$ < 30 cm$^{-3}$, with an almost identical distribution for the DIG and HII regions. The velocity dispersion of the ionised gas indicates that the DIG has a higher degree of turbulence than the HII regions. Comparison with photoionisation and shock models reveals that the diffuse component can only partially be explained via shocks and that it is most likely consistent with photons leaking from density bounded HII regions or with radiation from evolved field stars. [abridged] △ Less

Submitted 20 February, 2020; originally announced February 2020.

Comments: Accepted for publication in A&A. 24 pages, 19 figures. Please see journal article for high resolution figures

Journal ref: A&A 652, C6 (2021)

High Resolution Spectral Line Indices Useful for the Analysis of Stellar Populations

Authors: Lino H. Rodriguez-Merino, Y. D. Mayya, Paula R. T. Coelho, Gustavo Bruzual, Stephane Charlot, Esperanza Carrasco, Armando Gil de Paz

Abstract: The well-known age-metallicity-attenuation degeneracy does not permit unique and good estimates of basic parameters of stars and stellar populations. The effects of dust can be avoided using spectral line indices, but current methods have not been able to break the age-metallicity degeneracy. Here we show that using at least two new spectral line indices defined and measured on high-resolution (R=… ▽ More The well-known age-metallicity-attenuation degeneracy does not permit unique and good estimates of basic parameters of stars and stellar populations. The effects of dust can be avoided using spectral line indices, but current methods have not been able to break the age-metallicity degeneracy. Here we show that using at least two new spectral line indices defined and measured on high-resolution (R= 6000) spectra of a signal-to-noise ratio (S/N) > 10 one gets unambiguous estimates of the age and metallicity of intermediate to old stellar populations. Spectroscopic data retrieved with new astronomical facilities, e.g., X-shooter, MEGARA, and MOSAIC, can be employed to infer the physical parameters of the emitting source by means of spectral line index and index--index diagram analysis. △ Less

Submitted 30 January, 2020; originally announced January 2020.

Comments: To be published in The Astrophysical Journal Letters

Spatially-resolved analysis of neutral winds, stars and ionized gas kinematics with MEGARA/GTC: new insights on the nearby galaxy UGC 10205

Authors: C. Catalán-Torrecilla, Á. Castillo-Morales, A. Gil de Paz, J. Gallego, E. Carrasco, J. Iglesias-Páramo, R. Cedazo, M. Chamorro-Cazorla, S. Pascual, M. L. García-Vargas, N. Cardiel, P. Gómez-Alvarez, A. Pérez-Calpena, I. Martínez-Delgado, B. T. Dullo, P. Coelho, G. Bruzual, S. Charlot

Abstract: We present a comprehensive analysis of the multi-phase structure of the interstellar medium (ISM) and the stellar kinematics in the edge-on nearby galaxy UGC 10205 using integral field spectroscopy (IFS) data taken with MEGARA at the GTC. We explore both the neutral and the ionized gas phases using the interstellar Na ${\small I}$ D doublet absorption (LR$-$V set-up, R $\sim$ 6000) and the H$α$ em… ▽ More We present a comprehensive analysis of the multi-phase structure of the interstellar medium (ISM) and the stellar kinematics in the edge-on nearby galaxy UGC 10205 using integral field spectroscopy (IFS) data taken with MEGARA at the GTC. We explore both the neutral and the ionized gas phases using the interstellar Na ${\small I}$ D doublet absorption (LR$-$V set-up, R $\sim$ 6000) and the H$α$ emission line (HR$-$R set-up, R $\sim$ 18000), respectively. The high-resolution data show the complexity of the H$α$ emission line profile revealing the detection of up to three kinematically distinct gaseous components. Despite of this fact, a thin disk model is able to reproduce the bulk of the ionized gas motions in the central regions of UGC 10205. The use of asymmetric drift corrections is needed to reconciliate the ionized and the stellar velocity rotation curves. We also report the detection of outflowing neutral gas material blueshifted by $\sim$ 87 km s$^{-1}$. The main physical properties that describe the observed outflow are a total mass M$_{out}$ $=$ (4.55 $\pm$ 0.06) $\times$ 10$^{7}$ M$_{\odot}$ and a cold gas mass outflow rate $\dot{M}$$_{out}$ $=$ 0.78 $\pm$ 0.03 M$_{\odot}$ yr$^{-1}$. This work points out the necessity of exploiting high-resolution IFS data to understand the multi-phase components of the ISM and the multiple kinematical components in the central regions of nearby galaxies. △ Less

Submitted 18 December, 2019; originally announced December 2019.

Comments: 16 pages, 10 figures

To use or not to use synthetic stellar spectra in population synthesis models?

Authors: Paula R. T. Coelho, Gustavo Bruzual, Stephane Charlot

Abstract: Stellar population synthesis (SPS) models are invaluable to study star clusters and galaxies. They provide means to extract stellar masses, stellar ages, star formation histories, chemical enrichment and dust content of galaxies from their integrated spectral energy distributions, colours or spectra. As most models, they contain uncertainties which can hamper our ability to model and interpret obs… ▽ More Stellar population synthesis (SPS) models are invaluable to study star clusters and galaxies. They provide means to extract stellar masses, stellar ages, star formation histories, chemical enrichment and dust content of galaxies from their integrated spectral energy distributions, colours or spectra. As most models, they contain uncertainties which can hamper our ability to model and interpret observed spectra. This work aims at studying a specific source of model uncertainty: the choice of an empirical vs. a synthetic stellar spectral library. Empirical libraries suffer from limited coverage of parameter space, while synthetic libraries suffer from modelling inaccuracies. Given our current inability to have both ideal stellar-parameter coverage with ideal stellar spectra, what should one favour: better coverage of the parameters (synthetic library) or better spectra on a star-by-star basis (empirical library)? To study this question, we build a synthetic stellar library mimicking the coverage of an empirical library, and SPS models with different choices of stellar library tailored to these investigations. Through the comparison of model predictions and the spectral fitting of a sample of nearby galaxies, we learned that: predicted colours are more affected by the coverage effect than the choice of a synthetic vs. empirical library; the effects on predicted spectral indices are multiple and defy simple conclusions; derived galaxy ages are virtually unaffected by the choice of the library, but are underestimated when SPS models with limited parameter coverage are used; metallicities are robust against limited HRD coverage, but are underestimated when using synthetic libraries. △ Less

Submitted 25 October, 2019; originally announced October 2019.

Comments: Accepted for publication in MNRAS

Constraints on the production and escape of ionizing radiation from the emission-line spectra of metal-poor star-forming galaxies

Authors: Adèle Plat, Stephane Charlot, Gustavo Bruzual, Anna Feltre, Alba Vidal-García, Christophe Morisset, Jacopo Chevallard, Helge Todt

Abstract: We explore the production and escape of ionizing photons in young galaxies by investigating the ultraviolet and optical emission-line properties of models of ionization-bounded and density-bounded HII regions, active-galactic-nucleus (AGN) narrow-line regions and radiative shocks computed all using the same physically-consistent description of element abundances and depletion on to dust grains dow… ▽ More We explore the production and escape of ionizing photons in young galaxies by investigating the ultraviolet and optical emission-line properties of models of ionization-bounded and density-bounded HII regions, active-galactic-nucleus (AGN) narrow-line regions and radiative shocks computed all using the same physically-consistent description of element abundances and depletion on to dust grains down to very low metallicities. We compare these models with a reference sample of metal-poor star-forming galaxies and Lyman-continuum (LyC) leakers at various redshifts, which allows the simultaneous exploration of more spectral diagnostics than typically available at once for individual subsamples. We confirm that current single- and binary-star population synthesis models do not produce hard-enough radiation to account for the high-ionization emission of the most metal-poor galaxies. Introducing either an AGN or radiative-shock component brings models into agreement with observations. A published model including X-ray binaries is an attractive alternative to reproduce the observed rise in HeII4686/Hbeta ratio with decreasing oxygen abundance in metal-poor star-forming galaxies, but not the high observed HeII4686/Hbeta ratios of galaxies with large EW(Hbeta). A source of harder ionizing radiation appears to be required in these extreme objects, such as an AGN or radiative-shock component, perhaps linked to an initial-mass-function bias toward massive stars at low metallicity. This would also account for the surprisingly high [OI]/[OIII] ratios of confirmed LyC leakers relative to ionization-bounded models. We find no simple by-eye diagnostic of the nature of ionizing sources and the escape of LyC photon, which require proper simultaneous fits of several lines to be discriminated against. △ Less

Submitted 16 September, 2019; originally announced September 2019.

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

Chronos: A NIR spectroscopic survey to target the most important phases of galaxy evolution across cosmic time

Authors: I. Ferreras, M. Cropper, R. Sharples, J. Bland-Hawthorn, G. Bruzual, S. Charlot, C. J. Conselice, S. Driver, J. Dunlop, A. M. Hopkins, S. Kaviraj, T. Kitching, F. La Barbera, O. Lahav, A. Pasquali, S. Serjeant, J. Silk, R. Windhorst

Abstract: (Abridged summary) Responding to ESA's Voyage 2050 call to define the long-term plan for the future space missions that will address the astrophysics science questions during the 2035-2050 cycle, we propose a dedicated, ultra-deep spectroscopic survey in the near infrared (NIR), that will target a mass-limited sample of galaxies during two of the most fundamental epochs of cosmic evolution: the fo… ▽ More (Abridged summary) Responding to ESA's Voyage 2050 call to define the long-term plan for the future space missions that will address the astrophysics science questions during the 2035-2050 cycle, we propose a dedicated, ultra-deep spectroscopic survey in the near infrared (NIR), that will target a mass-limited sample of galaxies during two of the most fundamental epochs of cosmic evolution: the formation of the first galaxies (at z>6; cosmic dawn), and at the peak of galaxy formation activity (between redshift z=1 and 3; cosmic noon). By way of NIR observations, it is possible to study the Lyman-alpha region in the former, and the optical rest-frame in the latter, allowing us to extract fundamental observables such as gas and stellar kinematics, chemical abundances, and ages, providing a unique legacy database covering these two crucial stages of cosmic evolution. A dedicated, space-based facility will overcome the challenges faced by ground-based telescopes, no matter how large the aperture, or the reduced field of view and low multiplex factor of the best space-based instrument in the near future, namely NIRSpec at the JWST. Our project (codename Chronos) aims to produce about 1 million high quality spectra, with a high S/N in the continuum, where information about the underlying stellar populations is encoded. △ Less

Submitted 22 August, 2019; originally announced August 2019.

Comments: 24 pages, 14 figures. Submitted to ESA's Voyage 2050 call to shape the future space science programme. arXiv admin note: text overlap with arXiv:1306.6333

Detection of the self-regulation of star formation in galaxy discs

Authors: Javier Zaragoza-Cardiel, Jacopo Fritz, Itziar Aretxaga, Divakara Mayya, Daniel Rosa-González, John E. Beckman, Gustavo Bruzual, Stephane Charlot, Luis Lomelí-Núñez

Abstract: Stellar feedback has a notable influence on the formation and evolution of galaxies. However, direct observational evidence is scarce. We have performed stellar population analysis using MUSE optical spectra of the spiral galaxy NGC 628 and find that current maximum star formation in spatially resolved regions is regulated according to the level of star formation in the recent past. We propose a m… ▽ More Stellar feedback has a notable influence on the formation and evolution of galaxies. However, direct observational evidence is scarce. We have performed stellar population analysis using MUSE optical spectra of the spiral galaxy NGC 628 and find that current maximum star formation in spatially resolved regions is regulated according to the level of star formation in the recent past. We propose a model based on the self-regulator or "bathtub" models, but for spatially resolved regions of the galaxy. We name it the "resolved self-regulator model" and show that the predictions of this model are in agreement with the presented observations. We observe star formation self-regulation and estimate the mass-loading factor, $η=2.5 \pm 0.5$, consistent with values predicted by galaxy formation models. The method described here will help provide better constraints on those models. △ Less

Submitted 4 June, 2019; originally announced June 2019.

Comments: Accepted for publication in MNRAS Letters. 6 pages, 3 figures

Early Low-Mass Galaxies and Star-Cluster Candidates at z~6-9 Identified by the Gravitational Lensing Technique and Deep Optical/Near-Infrared Imaging

Authors: Shotaro Kikuchihara, Masami Ouchi, Yoshiaki Ono, Ken Mawatari, Jacopo Chevallard, Yuichi Harikane, Takashi Kojima, Masamune Oguri, Gustavo Bruzual, Stephane Charlot

Abstract: We present very faint dropout galaxies at z~6-9 with a stellar mass M* down to M*~10^6Mo that are found in deep optical/near-infrared (NIR) images of the full data sets of the Hubble Frontier Fields (HFF) program in conjunction with deep ground-based and Spitzer images and gravitational lensing magnification effects. We investigate stellar populations of the HFF dropout galaxies with the optical/N… ▽ More We present very faint dropout galaxies at z~6-9 with a stellar mass M* down to M*~10^6Mo that are found in deep optical/near-infrared (NIR) images of the full data sets of the Hubble Frontier Fields (HFF) program in conjunction with deep ground-based and Spitzer images and gravitational lensing magnification effects. We investigate stellar populations of the HFF dropout galaxies with the optical/NIR photometry and BEAGLE models made of self-consistent stellar population synthesis and photoionization models, carefully including strong nebular emission impacting on the photometry. We identify 357 galaxies with M*~10^6-10^9Mo, and find that a stellar mass to UV luminosity L_UV ratio M*/L_UV is nearly constant at M*~10^6-10^9Mo. Our best-estimate M*/L_UV function is comparable to a model of star-formation duration time of 100 Myr, but 2-7 times higher than the one of 10 Myr assumed in a previous study (at the 5sigma level) that would probably underestimate M* of faint galaxies. We derive the galaxy stellar mass functions (GSMFs) at z~6-9 that agree with those obtained by previous studies with no M*/L_UV assumptions at M*>~10^8Mo, and that extends to M*~10^6Mo. Estimating the stellar mass densities rho* with the GSMFs, we find that rho* smoothly increases from log(rho*/[Mo Mpc^(-3)])=5.91 +0.75/-0.65 at z~9 to 6.21 +0.39/-0.37 at z~6-7, which is consistent with the one estimated from star-formation rate density measurements. In conjunction with the estimates of the galaxy effective radii R_e on the source plane, we have pinpointed two objects with low stellar masses (M*<=10^7Mo) and very compact morphologies (R_e<=40 physical pc) that are comparable with those of globular clusters (GCs) in the Milky Way today. These objects are candidates of star clusters that should be a part or a dominant component of high-redshift low-mass galaxy, some of which may be related to GCs today. △ Less

Submitted 16 May, 2019; originally announced May 2019.

Parameter Estimation for Scarce Stellar Populations

Authors: V. H. Ramírez-Siordia, G. Bruzual, B. Cervantes Sodi, T. Bitsakis

Abstract: We present a Bayesian method to determine simultaneously the age, metallicity, distance modulus, and interstellar reddening by dust of any resolved stellar population, by comparing the observed and synthetic color magnitude diagrams on a star by star basis, with no need to bin the data into a carefully selected magnitude grid. We test the method with mock stellar populations, and show that it work… ▽ More We present a Bayesian method to determine simultaneously the age, metallicity, distance modulus, and interstellar reddening by dust of any resolved stellar population, by comparing the observed and synthetic color magnitude diagrams on a star by star basis, with no need to bin the data into a carefully selected magnitude grid. We test the method with mock stellar populations, and show that it works correctly even for scarce stellar populations with only one or two hundred stars above the main sequence turn off. If the population is the result of two star formation bursts, we can infer the contribution of each event to the total stellar population. The code works automatically and has already been used to study massive amounts of Magellanic clouds photometric data. In this paper we analyze in detail three Large Magellanic Cloud star clusters and 6 Ultra Faint Dwarf Galaxies. For these galaxies we recover physical parameters in agreement with those quoted in the literature, age $\sim13.7$ Gyr and a very low metallicity $log\,Z\sim-4$. Searching for multiple populations in these galaxies, we find, at a very low significance level, signs of a double stellar population for Ursa Major I: a dominant old population and a younger one which contributes $\sim25$% of the stars, in agreement with independent results from other authors. △ Less

Submitted 18 April, 2019; originally announced April 2019.

Comments: 14 pages, 11 figures. Accepted for publication in MNRAS

Spectroscopy of NGC 4258 globular cluster candidates: membership confirmation and kinematics

Authors: Rosa A. Gonzaléz-Lópezlira, Y. D. Mayya, Laurent Loinard, Karla Álamo-Martínez, George Heald, Iskren Y. Georgiev, Yasna Órdenes-Briceño, Ariane Lancon, Maritza A. Lara-López, Luis Lomelí-Núñez, Gustavo Bruzual, Thomas H. Puzia

Abstract: We present multi-object spectroscopic observations of 23 globular cluster candidates (GCCs) in the prototypical megamaser galaxy NGC 4258, carried out with the OSIRIS instrument at the 10.4 m Gran Telescopio Canarias. The candidates have been selected based on the ($u^*\ -\ i^\prime$) versus ($i^\prime\ -\ K_s$) diagram, in the first application of the uiks-method to a spiral galaxy. In the spectr… ▽ More We present multi-object spectroscopic observations of 23 globular cluster candidates (GCCs) in the prototypical megamaser galaxy NGC 4258, carried out with the OSIRIS instrument at the 10.4 m Gran Telescopio Canarias. The candidates have been selected based on the ($u^*\ -\ i^\prime$) versus ($i^\prime\ -\ K_s$) diagram, in the first application of the uiks-method to a spiral galaxy. In the spectroscopy presented here, 70% of the candidates are confirmed as globular clusters. Our results validate the efficiency of the uiks-method in the sparser GC systems of spirals, and given the downward correction to the total number of GCs, the agreement of the galaxy with the correlations between black hole mass, and total number and mass of GCs is actually improved. We find that the GCs, mostly metal-poor, co-rotate with the HI disk, even at large galactocentric distances. △ Less

Submitted 9 April, 2019; v1 submitted 4 April, 2019; originally announced April 2019.

Comments: 25 pages, 15 figures; accepted for publication in The Astrophysical Journal

From kpcs to the central parsec of NGC 1097: feeding star formation and a black hole at the same time

Authors: M. Almudena Prieto, Juan A. Fernandez-Ontiveros, Gustavo Bruzual, Andreas Burkert, Marc Schartmann, Stephan Charlot

Abstract: A panchromatic view of the star forming ring and feeding process in the central kpc of the galaxy NGC 1097 is presented. The assembled IR to UV images at $\sim$10 pc resolution allow us to characterise the population of circa 250 clusters in the ring and disentangle the network of filaments of dust and gas that enshroud and feed them. The ring is a place of intermittent star bursts over the last 1… ▽ More A panchromatic view of the star forming ring and feeding process in the central kpc of the galaxy NGC 1097 is presented. The assembled IR to UV images at $\sim$10 pc resolution allow us to characterise the population of circa 250 clusters in the ring and disentangle the network of filaments of dust and gas that enshroud and feed them. The ring is a place of intermittent star bursts over the last 100 Myr. Four major episodes covering a proto-cluster phase of eleven mid-IR sources at the molecular clouds core, and two (three) previous bursts with a time separation of 20 - 30 Myr are identified. The extinction map of the inner few kpc resolves NGC1097's two major dust lanes in bundles of narrow, $<$25 pc width, filaments running along the galaxy's bar. As they approach the ring, some circularise along it, others curve to the centre to produce a nuclear spiral. We believe these are kpc-scale dust-gas streamers feeding the ring and the black-hole. The total mass in clusters formed in the ring in the last 100 Myr is $< 10^7\, \rm{M_\odot}$, i.e. $< 1\% $ of the $10^{9} M_\odot$ of molecular gas in the ring; yet, at its current star formation rate, $\sim1.8\, \rm{M_\odot \, yr^{-1}}$, an order of magnitude more in stellar mass should have been produced over that period. This means that the availability of gas in the ring is not the sole star formation driver, perhaps the rate at which dense gas accumulates in the ring is key. △ Less

Submitted 25 February, 2019; originally announced February 2019.

Comments: 16 pages, 9 color figures. To appear in MNRAS

Simultaneous analysis of SDSS spectra and GALEX photometry with STARLIGHT: Method and early results

Authors: Ariel Werle, Roberto Cid Fernandes, Natalia Vale Asari, Gustavo Bruzual, Stéphane Charlot, Rosa Gonzalez Delgado, Fabio Herpich

Abstract: We combine data from the Sloan Digital Sky Survey and the Galaxy Evolution Explorer to simultaneously analyze optical spectra and ultraviolet photometry of 231643 galaxies with the STARLIGHT spectral synthesis code using state-of-the-art stellar population models. We present a new method to estimate GALEX photometry in the SDSS spectroscopic aperture, which proves quite reliable if applied to larg… ▽ More We combine data from the Sloan Digital Sky Survey and the Galaxy Evolution Explorer to simultaneously analyze optical spectra and ultraviolet photometry of 231643 galaxies with the STARLIGHT spectral synthesis code using state-of-the-art stellar population models. We present a new method to estimate GALEX photometry in the SDSS spectroscopic aperture, which proves quite reliable if applied to large samples. In agreement with previous experiments with CALIFA, we find that adding UV constraints leads to a moderate increase on the fraction of $\sim 10^7$ - $10^8$ yr populations and a concomitant decrease of younger and older components, yielding slightly older luminosity weighted mean stellar ages. These changes are most relevant in the low-mass end of the blue cloud. An increase in dust attenuation is observed for galaxies dominated by young stars. We investigate the contribution of different stellar populations to the fraction of light in GALEX and SDSS bands across the UV-optical color-magnitude diagram. As an example application, we use this $λ$ dependence to highlight differences between retired galaxies with and without emission lines. In agreement with an independent study by Herpich et al., we find that the former show an excess of intermediate age populations when compared to the later. Finally, we test the suitability of two different prescription for dust, finding that our dataset is best fitted using the attenuation law of starburst galaxies. However, results for the Milky Way extinction curve improve with decreasing $τ_V$, especially for edge-on galaxies. △ Less

Submitted 14 December, 2018; v1 submitted 27 November, 2018; originally announced November 2018.

Comments: Accepted for publication in MNRAS

Modelling dust rings in early-type galaxies through a sequence of radiative transfer simulations and 2D image fitting

Authors: P. Bonfini, O. Gonzalez-Martin, J. Fritz, T. Bitsakis, G. Bruzual, B. Cervantes Sodi

Abstract: A large fraction of early-type galaxies (ETGs) host prominent dust features, and central dust rings are arguably the most interesting among them. We present here `Lord Of The Rings' (LOTR), a new methodology which allows to integrate the extinction by dust rings in a 2D fitting modelling of the surface brightness distribution. Our pipeline acts in two steps, first using the surface fitting softwar… ▽ More A large fraction of early-type galaxies (ETGs) host prominent dust features, and central dust rings are arguably the most interesting among them. We present here `Lord Of The Rings' (LOTR), a new methodology which allows to integrate the extinction by dust rings in a 2D fitting modelling of the surface brightness distribution. Our pipeline acts in two steps, first using the surface fitting software GALFIT to determine the unabsorbed stellar emission, and then adopting the radiative transfer code SKIRT to apply dust extinction. We apply our technique to NGC 4552 and NGC 4494, two nearby ETGs. We show that the extinction by a dust ring can mimic, in a surface brightness profile, a central point source (e.g. an unresolved nuclear stellar cluster or an active galactic nucleus; AGN) superimposed to a `core' (i.e. a central flattening of the stellar light commonly observed in massive ETGs). We discuss how properly accounting for dust features is of paramount importance to derive correct fluxes especially for low luminosity AGNs (LLAGNs). We suggest that the geometries of dust features are strictly connected with how relaxed is the gravitational potential, i.e. with the evolutionary stage of the host galaxy. Additionally, we find hints that the dust mass contained in the ring relates to the AGN activity. △ Less

Submitted 26 April, 2018; originally announced April 2018.

Comments: Accepted in MNRAS

J-PLUS: The Javalambre Photometric Local Universe Survey

Authors: A. J. Cenarro, M. Moles, D. Cristóbal-Hornillos, A. Marín-Franch, A. Ederoclite, J. Varela, C. López-Sanjuan, C. Hernández-Monteagudo, R. E. Angulo, H. Vázquez Ramió, K. Viironen, S. Bonoli, A. A. Orsi, G. Hurier, I. San Roman, N. Greisel, G. Vilella-Rojo, L. A. Díaz-García, R. Logroño-García, S. Gurung-López, D. Spinoso, D. Izquierdo-Villalba, J. A. L. Aguerri, C. Allende Prieto, C. Bonatto , et al. (97 additional authors not shown)

Abstract: J-PLUS is an ongoing 12-band photometric optical survey, observing thousands of square degrees of the Northern hemisphere from the dedicated JAST/T80 telescope at the Observatorio Astrofísico de Javalambre. T80Cam is a 2 sq.deg field-of-view camera mounted on this 83cm-diameter telescope, and is equipped with a unique system of filters spanning the entire optical range. This filter system is a com… ▽ More J-PLUS is an ongoing 12-band photometric optical survey, observing thousands of square degrees of the Northern hemisphere from the dedicated JAST/T80 telescope at the Observatorio Astrofísico de Javalambre. T80Cam is a 2 sq.deg field-of-view camera mounted on this 83cm-diameter telescope, and is equipped with a unique system of filters spanning the entire optical range. This filter system is a combination of broad, medium and narrow-band filters, optimally designed to extract the rest-frame spectral features (the 3700-4000Å Balmer break region, H$δ$, Ca H+K, the G-band, the Mgb and Ca triplets) that are key to both characterize stellar types and to deliver a low-resolution photo-spectrum for each pixel of the sky observed. With a typical depth of AB $\sim 21.25$ mag per band, this filter set thus allows for an indiscriminate and accurate characterization of the stellar population in our Galaxy, it provides an unprecedented 2D photo-spectral information for all resolved galaxies in the local universe, as well as accurate photo-z estimates ($Δ\,z\sim 0.01-0.03$) for moderately bright (up to $r\sim 20$ mag) extragalactic sources. While some narrow band filters are designed for the study of particular emission features ([OII]/$λ$3727, H$α$/$λ$6563) up to $z < 0.015$, they also provide well-defined windows for the analysis of other emission lines at higher redshifts. As a result, J-PLUS has the potential to contribute to a wide range of fields in Astrophysics, both in the nearby universe (Milky Way, 2D IFU-like studies, stellar populations of nearby and moderate redshift galaxies, clusters of galaxies) and at high redshifts (ELGs at $z\approx 0.77, 2.2$ and $4.4$, QSOs, etc). With this paper, we release $\sim 36$ sq.deg of J-PLUS data, containing about $1.5\times 10^5$ stars and $10^5$ galaxies at $r<21$ mag. △ Less

Submitted 8 April, 2018; originally announced April 2018.

Comments: Submitted to A&A

Journal ref: A&A 622, A176 (2019)