Showing 1–50 of 79 results for author: Fernández-Ontiveros, J

J-PLUS: The fraction of calcium white dwarfs along the cooling sequence

Authors: C. López-Sanjuan, P. -E. Tremblay, M. W. O'Brien, D. Spinoso, A. Ederoclite, H. Vázquez Ramió, A. J. Cenarro, A. Marín-Franch, T. Civera, J. M. Carrasco, B. T. Gänsicke, N. P. Gentile Fusillo, A. Hernán-Caballero, M. A. Hollands, A. del Pino, H. Domínguez Sánchez, J. A. Fernández-Ontiveros, F. M. Jiménez-Esteban, A. Rebassa-Mansergas, L. Schmidtobreick, R. E. Angulo, D. Cristòbal-Hornillos, R. A. Dupke, C. Hernández-Monteagudo, M. Moles , et al. (2 additional authors not shown)

Abstract: We used the Javalambre Photometric Local Universe Survey (J-PLUS) DR2 photometry in twelve optical bands over 2176 deg2 to estimate the fraction of white dwarfs with presence of CaII H+K absorption along the cooling sequence. We compared the J-PLUS photometry against metal-free theoretical models to estimate the equivalent width in the J0395 passband of 10 nm centered at 395 nm (EW_J0395), a proxy… ▽ More We used the Javalambre Photometric Local Universe Survey (J-PLUS) DR2 photometry in twelve optical bands over 2176 deg2 to estimate the fraction of white dwarfs with presence of CaII H+K absorption along the cooling sequence. We compared the J-PLUS photometry against metal-free theoretical models to estimate the equivalent width in the J0395 passband of 10 nm centered at 395 nm (EW_J0395), a proxy to detect calcium absorption. A total of 4399 white dwarfs within 30000 > Teff > 5500 K and mass M > 0.45 Msun were analyzed. Their EW_J0395 distribution was modeled using two populations, corresponding to polluted and non-polluted systems, to estimate the fraction of calcium white dwarfs (f_Ca) as a function of Teff. The probability for each individual white dwarf of presenting calcium absorption, pca, was also computed. The comparison with both the measured Ca/He abundance and the metal pollution from spectroscopy shows that EW_J0395 correlates with the presence of calcium. The fraction of calcium white dwarfs increases from f_Ca = 0 at Teff = 13500 K to f_Ca = 0.15 at Teff = 5500 K. We compare our results with the fractions derived from the 40 pc spectroscopic sample and from SDSS spectra. The trend found in J-PLUS observations is also present in the 40 pc sample, however SDSS shows a deficit of metal-polluted objects at Teff < 12000 K. Finally, we found 39 white dwarfs with pca > 0.99. Twenty of them have spectra presented in previous studies, whereas we observed six additional targets. These 26 objects were all confirmed as metal-polluted systems. The J-PLUS optical data provide a robust statistical measurement for the presence of CaII H+K absorption in white dwarfs. We find a 15 +- 3 % increase in the fraction of calcium white dwarfs from Teff = 13500 K to 5500 K, which reflects their selection function in the optical from the total population of metal-polluted systems. △ Less

Submitted 23 June, 2024; originally announced June 2024.

Comments: Submitted to Astronomy and Astrophysics. Comments are welcome

A case study of gas impacted by black-hole jets with the JWST: outflows, bow shocks, and high excitation of the gas in the galaxy IC5063

Authors: K. M. Dasyra, G. F. Paraschos, F. Combes, P. Patapis, G. Helou, M. Papachristou, J. A. Fernandez-Ontiveros, T. G. Bisbas, L. Spinoglio, L. Armus, M. Malkan

Abstract: We present James Webb Space Telescope MIRI data of the inner ~3x2 kpc^2 of the galaxy IC5063, in which the jets of a supermassive black hole interact with the gaseous disk they are crossing. Jet-driven outflows were known to be initiated along or near the jet path, and the stability conditions of clouds were known to vary because of these outflows. The MIRI data, of unprecedented resolution and se… ▽ More We present James Webb Space Telescope MIRI data of the inner ~3x2 kpc^2 of the galaxy IC5063, in which the jets of a supermassive black hole interact with the gaseous disk they are crossing. Jet-driven outflows were known to be initiated along or near the jet path, and the stability conditions of clouds were known to vary because of these outflows. The MIRI data, of unprecedented resolution and sensitivity in the infrared, now reveal that there are more than ten discrete regions with outflows, nearly doubling the number of such known regions. Outflows exist near the radio lobes, at the nucleus, in a biconical structure perpendicular to the jet, and in a bubble moving against the disk. In some of them, velocities above escape velocity are observed. Stratification is also observed, with higher ionization or excitation gas attaining higher velocities. More outflows and bow shocks, found further away from the nucleus than the radio lobes, in regions without significant radio emission, reveal the existence of past or weak radio jets that interacted with the interstellar medium. The coincidence of the bow shocks with the optical extended emission line region (EELR) suggests that the jets also contributed to the gas ionization. Maps of the H2 gas excitation temperature, T_ex, indicate that the molecular gas is most excited in regions with radio emission. There, T_ex is more than 100K higher than in the EELR interior. We argue that a combination of jet-related shocks and cosmic rays is likely responsible for this excess molecular gas excitation. △ Less

Submitted 5 June, 2024; originally announced June 2024.

Comments: Submitted to ApJ

J-PLUS: Bayesian object classification with a strum of BANNJOS

Authors: A. del Pino, C. López-Sanjuan, A. Hernán-Caballero, H. Domínguez-Sánchez, R. von Marttens, J. A. Fernández-Ontiveros, P. R. T. Coelho, A. Lumbreras-Calle, J. Vega-Ferrero, F. Jimenez-Esteban, P. Cruz, V. Marra, M. Quartin, C. A. Galarza, R. E. Angulo, A. J. Cenarro, D. Cristóbal-Hornillos, R. A. Dupke, A. Ederoclite, C. Hernández-Monteagudo, A. Marín-Franch, M. Moles, L. Sodré Jr., J. Varela, H. Vázquez Ramió

Abstract: With its 12 optical filters, the Javalambre-Photometric Local Universe Survey (J-PLUS) provides an unprecedented multicolor view of the local Universe. The third data release (DR3) covers 3,192 deg$^2$ and contains 47.4 million objects. However, the classification algorithms currently implemented in its pipeline are deterministic and based solely on the sources morphology. Our goal is classify the… ▽ More With its 12 optical filters, the Javalambre-Photometric Local Universe Survey (J-PLUS) provides an unprecedented multicolor view of the local Universe. The third data release (DR3) covers 3,192 deg$^2$ and contains 47.4 million objects. However, the classification algorithms currently implemented in its pipeline are deterministic and based solely on the sources morphology. Our goal is classify the sources identified in the J-PLUS DR3 images into stars, quasi-stellar objects (QSOs), and galaxies. For this task, we present BANNJOS, a machine learning pipeline that uses Bayesian neural networks to provide the probability distribution function (PDF) of the classification. BANNJOS is trained on photometric, astrometric, and morphological data from J-PLUS DR3, Gaia DR3, and CatWISE2020, using over 1.2 million objects with spectroscopic classification from SDSS DR18, LAMOST DR9, DESI EDR, and Gaia DR3. Results are validated using $1.4 10^5$ objects and cross-checked against theoretical model predictions. BANNJOS outperforms all previous classifiers in terms of accuracy, precision, and completeness across the entire magnitude range. It delivers over 95% accuracy for objects brighter than $r = 21.5$ mag, and ~90% accuracy for those up to $r = 22$ mag, where J-PLUS completeness is < 25%. BANNJOS is also the first object classifier to provide the full probability distribution function (PDF) of the classification, enabling precise object selection for high purity or completeness, and for identifying objects with complex features, like active galactic nuclei with resolved host galaxies. BANNJOS has effectively classified J-PLUS sources into around 20 million galaxies, 1 million QSOs, and 26 million stars, with full PDFs for each, which allow for later refinement of the sample. The upcoming J-PAS survey, with its 56 color bands, will further enhance BANNJOS's ability to detail each source's nature. △ Less

Submitted 25 April, 2024; originally announced April 2024.

Comments: Submitted to Astronomy and Astrophysics. 29 pages, 23 figures, 4 tables

Chemical abundances and deviations from the solar S/O ratio in the gas-phase ISM of galaxies based on infrared emission lines

Authors: Borja Pérez-Díaz, Enrique Pérez-Montero, Juan A. Fernández-Ontiveros, José M. Vílchez, Antonio Hernán-Caballero, Ricardo Amorín

Abstract: The infrared (IR) range is extremely useful in the context of chemical abundance studies of the gas-phase interstellar medium (ISM) due to the large variety of ionic species traced in this regime, the negligible effects from dust attenuation or temperature stratification, and the amount of data that has been and will be released in the coming years. Taking advantage of available IR emission lines,… ▽ More The infrared (IR) range is extremely useful in the context of chemical abundance studies of the gas-phase interstellar medium (ISM) due to the large variety of ionic species traced in this regime, the negligible effects from dust attenuation or temperature stratification, and the amount of data that has been and will be released in the coming years. Taking advantage of available IR emission lines, we analysed the chemical content of the gas-phase ISM in a sample of 131 Star-Forming Galaxies (SFGs) and 73 Active Galactic Nuclei (AGNs). Particularly, we derived the chemical content via their total oxygen abundance in combination with nitrogen and sulfur abundances, and with the ionisation parameter. We used a new version of the code HII-CHI-Mistry-IR v3.1 which allows us to estimate log(N/O), 12+log(O/H), log(U), and, for the first time, 12+log(S/H) from IR emission lines, which can be applied to both SFGs and AGNs. We tested that the estimations from this new version, that only considers sulfur lines for the derivation of sulfur abundances, are compatible with previous studies. While most of the SFGs and AGNs show solar log(N/O) abundances, we found a large spread in the log(S/O) relative abundances. Specifically, we found extremely low log(S/O) values (1/10th solar) in some SFGs and AGNs with solar-like oxygen abundances. This result warns against the use of optical and IR sulfur emission lines to estimate oxygen abundances when no prior estimation of log(S/O) is provided. △ Less

Submitted 5 March, 2024; originally announced March 2024.

Comments: Accepted for publication in A&A. 16 pages, 11 figures, 4 electronic tables

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

Multiwavelength exploration of Extreme Emission Line Galaxies detected in miniJPAS survey

Authors: Iris Breda, Stergios Amarantidis, José M. Vilchez, Enrique Pérez-Montero, Carolina Kehrig, Jorge Iglesias-Páramo, Antonio Arroyo-Polonio, Juan A. Fernández-Ontiveros, Rosa M. González Delgado, Luis A. Díaz-García, Raul Abramo, 5 Jailson Alcaniz, Narciso Benítez, Silvia Bonoli, Javier A. Cenarro, David Cristóbal-Hornillos, Renato Dupke, Alessandro Ederoclite, Antonio Hernán-Caballero, Carlos López-Sanjuan, Antonio Marín-Franch, Claudia Mendes de Oliveira, Mariano Moles, Laerte Sodré, Keith Taylor , et al. (2 additional authors not shown)

Abstract: Extreme Emission Line Galaxies (EELGs) stand as remarkable objects due to their extremely metal poor environment and intense star formation. Considered as local analogues of high-redshift galaxies in the peak of their star-forming activity, they offer insights into conditions prevalent during the early Universe. Assessment of their stellar and gas properties is, therefore, of critical importance,… ▽ More Extreme Emission Line Galaxies (EELGs) stand as remarkable objects due to their extremely metal poor environment and intense star formation. Considered as local analogues of high-redshift galaxies in the peak of their star-forming activity, they offer insights into conditions prevalent during the early Universe. Assessment of their stellar and gas properties is, therefore, of critical importance, which requires the assembly of a considerable sample, comprehending a broad redshift range. The Javalambre-Physics of the Accelerating Universe Astrophysical Survey (JPAS) plays a significant role in assembling such a sample, encompassing approximately 8000 deg2 and employing 54 narrow-band optical filters. The present work describes the development and subsequent application of the tools that will be employed in the forthcoming JPAS spectrophotometric data, allowing for the massive and automated characterization of EELGs that are expected to be identified. This fully automated pipeline (requiring only the object coordinates from users) constructs Spectral Energy Distributions (SEDs) by retrieving virtually all the available multi-wavelength photometric data archives, employs SED fitting tools and identifies optical emission lines. It was applied to the sample of extreme line emitters identified in the miniJPAS Survey, and its derived physical properties such as stellar mass and age, coupled with fundamental relations, mirror results obtained through spectral modeling of SDSS spectra. Thorough testing using galaxies with documented photometric measurements across different wavelengths confirmed the pipeline's accuracy, demonstrating its capability for automated analysis of sources with varying characteristics, spanning brightness, morphology, and redshifts. The modular nature of this pipeline facilitates any addition from the user. △ Less

Submitted 24 January, 2024; originally announced January 2024.

Comments: 15 pages, 9 Figures, accepted for publication in MNRAS

Exploring the hardness of the ionising radiation with the infrared softness diagram. I. Similar effective temperature scales for starbursts and (ultra)luminous infrared galaxies

Authors: Enrique Pérez-Montero, Juan A. Fernández-Ontiveros, Borja Pérez-Díaz, José M. Vílchez, Nimisha Kumari, Ricardo Amorín

Abstract: {We explored the {softness parameter} in the infrared, whose main purpose is the characterisation of the hardness of the incident ionising radiation in emission-line nebulae. This parameter is obtained from the combination of mid-infrared wavelength range transitions corresponding to consecutive ionisation stages in star-forming regions. We compiled observational data from a sample of star-forming… ▽ More {We explored the {softness parameter} in the infrared, whose main purpose is the characterisation of the hardness of the incident ionising radiation in emission-line nebulae. This parameter is obtained from the combination of mid-infrared wavelength range transitions corresponding to consecutive ionisation stages in star-forming regions. We compiled observational data from a sample of star-forming galaxies (SFGs), including luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs), to study the softness parameter and its equivalent expression in two dimensions, the softness diagram. We compared them with predictions from photoionisation models to determine the shape of the ionising continuum energy distribution in each case. We also used the measured emission-line ratios as input for HCmistry-Teff-IR, a code that performs a Bayesian-like comparison with photoionisation model predictions in order to quantify the equivalent effective temperature (T*) and the ionisation parameter. We found similar average values within the errors of the softness parameter in (U)LIRGs (-0.57) in the rest of the SFGs (-0.51), which could be interpreted as indicative of a similar incident radiation field. This result is confirmed from the analysis using HCm-Teff-IR, which simultaneously points to a slightly lower, although similar within the errors, T* scale for (U)LIRGs, even when a higher dust-to-gas mass ratio is considered in the models for these objects. These derived T* values are compatible with the ionisation from massive stars, without any need of harder ionising sources, both for (U)LIRGs and the rest of the SFGs. However, the derived T* in (U)LIRGs do not show any correlation with metallicity. This could be interpreted as a sign that their similar average T* values are due to the attenuation of the energetic incident flux from massive stars by the heated dust mixed with the gas. △ Less

Submitted 18 January, 2024; originally announced January 2024.

Comments: Accepted for publication by Astronomy & Astrophysics. 9 pages, 6 figures

The spectral energy distributions and the bolometric luminosities of local AGN: study of the complete 12 micron AGN sample

Authors: Luigi Spinoglio, Juan Antonio Fernandez-Ontiveros, Matthew A. Malkan

Abstract: We measure the bolometric luminosity of a complete and unbiased 12 micron-selected sample of active galactic nuclei (AGN) in the local Universe. For each galaxy we used a 10-band radio-to-X-ray Spectral Energy Distribution (SED) to isolate the genuine AGN continuum in each band, including sub-arcsecond measurements where available, and correcting those contaminated by the host galaxy. We derive th… ▽ More We measure the bolometric luminosity of a complete and unbiased 12 micron-selected sample of active galactic nuclei (AGN) in the local Universe. For each galaxy we used a 10-band radio-to-X-ray Spectral Energy Distribution (SED) to isolate the genuine AGN continuum in each band, including sub-arcsecond measurements where available, and correcting those contaminated by the host galaxy. We derive the median SED of Seyfert type 1 AGN, Seyferts with hidden broad-lines (HBL), Seyferts of type 2, and LINER nuclei in our sample. The median Seyfert 1 SED shows the characteristic blue bump feature in the UV, but nevertheless the largest contribution to the bolometric luminosity comes from the IR and X-ray continua. The median SEDs of both HBL and type 2 AGN are affected by starlight contamination in the optical/UV. The median SED of HBL AGN is consistent with that of Seyfert 1's, when an extinction of Av = 1.2 mag is applied. The comprehensive SEDs allowed us to measure accurate bolometric luminosities and derive robust bolometric corrections for the different tracers. The 12 micron and the K-band nuclear luminosities have good linear correlations with the bolometric luminosity, similar to those in the X-rays. We derive bolometric corrections for either continuum bands (K-band, 12 micron, 2-10 keV and 14-195 keV) and narrow emission lines (mid-IR high ionization lines of [OIV] and [NeV] and optical [OIII]5007A) as well as for combinations of IR continuum and line emission. A combination of continuum plus line emission accurately predicts the bolometric luminosity up to quasar luminosities. △ Less

Submitted 21 December, 2023; originally announced December 2023.

Comments: Accepted for publication in the Astrophysical Journal on December 20, 2023

Fast infrared winds during the radio-loud and X-ray obscured stages of the black hole transient GRS 1915+105

Authors: J. Sánchez-Sierras, T. Muñoz-Darias, S. E. Motta, R. P. Fender, A. Bahramian, C. Martínez-Sebastián, J. A. Fernández-Ontiveros, J. Casares, M. Armas Padilla, D. A. Green, D. Mata Sánchez, J. Strader, M. A. P. Torres

Abstract: The black hole transient GRS 1915+105 entered a new phase of activity in 2018, generally characterised by low X-ray and radio fluxes. This phase has been only interrupted by episodes of strong and variable radio emission, during which high levels of X-ray absorption local to the source were measured. We present 18 epochs of near-infrared spectroscopy (2018-2023) obtained with GTC/EMIR and VLT/X-sh… ▽ More The black hole transient GRS 1915+105 entered a new phase of activity in 2018, generally characterised by low X-ray and radio fluxes. This phase has been only interrupted by episodes of strong and variable radio emission, during which high levels of X-ray absorption local to the source were measured. We present 18 epochs of near-infrared spectroscopy (2018-2023) obtained with GTC/EMIR and VLT/X-shooter, spanning both radio-loud and radio-quiet periods. We demonstrate that radio-loud phases are characterised by strong P-Cygni line profiles, indicative of accretion disc winds with velocities of up to $\mathrm{\sim 3000~km~s^{-1}}$. This velocity is consistent with those measured in other black hole transients. It is also comparable to the velocity of the X-ray winds detected during the peak outburst phases in GRS 1915+105, reinforcing the idea that massive, multi-phase outflows are characteristic features of the largest and most powerful black hole accretion discs. Conversely, the evolution of the Br$γ$ line profile during the radio-quiet phases follows the expected trend for accretion disc lines in a system that is gradually decreasing its intrinsic luminosity, exhibiting weaker intensities and more pronounced double-peaks. △ Less

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

Comments: Accepted for publication in A&A Letters

The miniJPAS survey: Maximising the photo-z accuracy from multi-survey datasets with probability conflation

Authors: A. Hernán-Caballero, M. Akhlaghi, C. López-Sanjuan, H. Vázquez-Ramió, J. Laur, J. Varela, T. Civera, D. Muniesa, A. Finoguenov, J. A. Fernández-Ontiveros, H. Domínguez-Sánchez, J. Chaves-Montero, A. Fernández-Soto, A. Lumbreras-Calle, L. A. Díaz-García, A. del Pino, R. M. González Delgado, C. Hernández-Monteagudo, P. Coelho, Y. Jiménez-Teja, P. A. A. Lopes, V. Marra, E. Tempel, J. M. Vílchez, R. Abramo , et al. (13 additional authors not shown)

Abstract: We present a new method for obtaining photometric redshifts (photo-z) for sources observed by multiple photometric surveys using a combination (conflation) of the redshift probability distributions (PDZs) obtained independently from each survey. The conflation of the PDZs has several advantages over the usual method of modelling all the photometry together, including modularity, speed, and accurac… ▽ More We present a new method for obtaining photometric redshifts (photo-z) for sources observed by multiple photometric surveys using a combination (conflation) of the redshift probability distributions (PDZs) obtained independently from each survey. The conflation of the PDZs has several advantages over the usual method of modelling all the photometry together, including modularity, speed, and accuracy of the results. Using a sample of galaxies with narrow-band photometry in 56 bands from J-PAS and deeper grizy photometry from the Hyper-SuprimeCam Subaru Strategic program (HSC-SSP), we show that PDZ conflation significantly improves photo-z accuracy compared to fitting all the photometry or using a weighted average of point estimates. The improvement over J-PAS alone is particularly strong for i>22 sources, which have low signal-to-noise ratio in the J-PAS bands. For the entire i<22.5 sample, we obtain a 64% (45%) increase in the number of sources with redshift errors |Dz|<0.003, a factor 3.3 (1.9) decrease in the normalised median absolute deviation of the errors (sigma_NMAD), and a factor 3.2 (1.3) decrease in the outlier rate compared to J-PAS (HSC-SSP) alone. The photo-z accuracy gains from combining the PDZs of J-PAS with a deeper broadband survey such as HSC-SSP are equivalent to increasing the depth of J-PAS observations by ~1.2--1.5 magnitudes. These results demonstrate the potential of PDZ conflation and highlight the importance of including the full PDZs in photo-z catalogues. △ Less

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

Comments: 15 pages, 16 figures, 1 table. Accepted for publication in A&A

Journal ref: A&A 684, A61 (2024)

A plausible link between dynamically unsettled molecular gas and the radio jet in NGC 6328

Authors: Michalis Papachristou, Kalliopi Maria Dasyra, Juan Antonio Fernández-Ontiveros, Anelise Audibert, Ilaria Ruffa, Francoise Combes, Markos Polkas, Athanasia Gkogkou

Abstract: We report the detection of outflowing molecular gas at the center of the nearby radio galaxy NGC6328 (z=0.014), which has a gigahertz-peaked spectrum radio core and a compact (2 pc) young double radio lobe tracing jet. Utilizing Atacama Large Millimeter/submillimeter Array (ALMA) CO(2-1) and CO(3-2) observations, as well as a novel code developed to fit the 3D gas distribution and kinematics, to s… ▽ More We report the detection of outflowing molecular gas at the center of the nearby radio galaxy NGC6328 (z=0.014), which has a gigahertz-peaked spectrum radio core and a compact (2 pc) young double radio lobe tracing jet. Utilizing Atacama Large Millimeter/submillimeter Array (ALMA) CO(2-1) and CO(3-2) observations, as well as a novel code developed to fit the 3D gas distribution and kinematics, to study the molecular gas kinematics, we find that the bulk of the gas is situated within a highly warped disk structure, most likely the result of a past merger event. Our analysis further uncovers, within the inner regions of the gas distribution (R<300 pc) and at a position angle aligning with that of the radio jet (150 degrees), the existence of two anti-diametric molecular gas structures kinematically detached from the main disk. These structures most likely trace a jet-induced cold gas outflow with a total lower limit mass of $2\times 10^6\,\mathrm{M_\odot}$ mass, corresponding to an outflow rate of $2\,\mathrm{M_\odot\,yr^{-1}}$ and a kinetic power of $2.7\times 10^{40}\,\mathrm{erg\,s^{-1}}$. The energy required to maintain such a molecular outflow is aligned with the mechanical power of the jet. △ Less

Submitted 3 October, 2023; originally announced October 2023.

Comments: 19 pages, 13 figures, Accepted for publication in A&

Journal ref: A&A 679, A115 (2023)

The miniJPAS survey: clusters and galaxy groups detection with AMICO

Authors: M. Maturi, A. Finoguenov, P. A. A. Lopes, R. M. González Delgado, R. A. Dupke, E. S. Cypriano, E. R. Carrasco, J. M. Diego, M. Penna-Lima, J. M. Vílchez, L. Moscardini, V. Marra, S. Bonoli, J. E. Rodríguez-Martín, A. Zitrin, I. Márquez, A. Hernán-Caballero, Y. Jiménez-Teja, R. Abramo, J. Alcaniz, N. Benitez, S. Carneiro, J. Cenarro, D. Cristóbal-Hornillos, A. Ederoclite , et al. (9 additional authors not shown)

Abstract: Samples of galaxy clusters allow us to better understand the physics at play in galaxy formation and to constrain cosmological models once their mass, position (for clustering studies) and redshift are known. In this context, large optical data sets play a crucial role. We investigate the capabilities of the Javalambre-Physics of the Accelerating Universe Astrophysical Survey (J-PAS) in detecting… ▽ More Samples of galaxy clusters allow us to better understand the physics at play in galaxy formation and to constrain cosmological models once their mass, position (for clustering studies) and redshift are known. In this context, large optical data sets play a crucial role. We investigate the capabilities of the Javalambre-Physics of the Accelerating Universe Astrophysical Survey (J-PAS) in detecting and characterizing galaxy groups and clusters. We analyze the data of the miniJPAS survey, obtained with the JPAS-Pathfinder camera and covering $1$ deg$^2$ centered on the AEGIS field to the same depths and with the same 54 narrow band plus 2 broader band near-UV and near-IR filters anticipated for the full J-PAS survey. We use the Adaptive Matched Identifier of Clustered Objects (AMICO) to detect and characterize groups and clusters of galaxies down to $S/N=2.5$ in the redshift range $0.05<z<0.8$. We detect 80, 30 and 11 systems with signal-to-noise ratio larger than 2.5, 3.0 and 3.5, respectively, down to $\sim 10^{13}\,M_{\odot}/h$. We derive mass-proxy scaling relations based on Chandra and XMM-Newton X-ray data for the signal amplitude returned by AMICO, the intrinsic richness and a new proxy that incorporates the galaxies' stellar masses. The latter proxy is made possible thanks to the J-PAS filters and shows a smaller scatter with respect to the richness. We fully characterize the sample and use AMICO to derive a probabilistic membership association of galaxies to the detected groups that we test against spectroscopy. We further show how the narrow band filters of J-PAS provide a gain of up to 100% in signal-to-noise ratio in detection and an uncertainty on the redshift of clusters of only $σ_z=0.0037(1+z)$ placing J-PAS in between broadband photometric and spectroscopic surveys. The performances of AMICO and J-PAS with respect to mass sensitivity, mass-proxies quality △ Less

Submitted 12 July, 2023; originally announced July 2023.

Comments: 15 pages, 12 figures, 3 tables, submitted to A&A

Journal ref: A&A 678, A145 (2023)

ReveaLLAGN 0: First Look at JWST MIRI data of Sombrero and NGC 1052

Authors: K. Goold, A. Seth, M. Molina, D. Ohlson, J. C. Runnoe, T. Boeker, T. A. Davis, A. Dumont, M. Eracleous, J. A. Fernández-Ontiveros, E. Gallo, A. D. Goulding, J. E. Greene, L. C. Ho, S. B. Markoff, N. Neumayer, R. Plotkin, A. Prieto, S. Satyapal, G. Van De Ven, J. L. Walsh, F. Yuan, A. Feldmeier-Krause, K. Gültekin, S. Hoenig , et al. (6 additional authors not shown)

Abstract: We present the first results from the Revealing Low-Luminosity Active Galactic Nuclei (ReveaLLAGN) survey, a JWST survey of seven nearby LLAGN. We focus on two observations with the Mid-Infrared Instrument's (MIRI) Medium Resolution Spectrograph (MRS) of the nuclei of NGC 1052 and Sombrero (NGC 4594 / M104). We also compare these data to public JWST data of a higher-luminosity AGN, NGC 7319 and NG… ▽ More We present the first results from the Revealing Low-Luminosity Active Galactic Nuclei (ReveaLLAGN) survey, a JWST survey of seven nearby LLAGN. We focus on two observations with the Mid-Infrared Instrument's (MIRI) Medium Resolution Spectrograph (MRS) of the nuclei of NGC 1052 and Sombrero (NGC 4594 / M104). We also compare these data to public JWST data of a higher-luminosity AGN, NGC 7319 and NGC 7469. JWST clearly separates the AGN spectrum from the galaxy light even in Sombrero, the faintest target in our survey; the AGN components have very red spectra. We find that the emission-line widths in both NGC 1052 and Sombrero increase with increasing ionization potential, with FWHM > 1000 km/s for lines with ionization potential > 50 eV. These lines are also significantly blue-shifted in both LLAGN. The high ionization potential lines in NGC 7319 show neither broad widths or significant blue shifts. Many of the lower ionization potential emission lines in Sombrero show significant blue wings extending > 1000 km/s. These features and the emission-line maps in both galaxies are consistent with outflows along the jet direction. Sombrero has the lowest luminosity high-ionization potential lines ([Ne V] and [O IV]) ever measured in the mid-IR, but the relative strengths of these lines are consistent with higher luminosity AGN. On the other hand, the [Ne V] emission is much weaker relative to the [Ne III] and [Ne II] lines of higher-luminosity AGN. These initial results show the great promise that JWST holds for identifying and studying the physical nature of LLAGN. △ Less

Submitted 1 March, 2024; v1 submitted 3 July, 2023; originally announced July 2023.

Comments: Accepted to ApJ Feb 28, 2024

Merger-driven infall of metal-poor gas in luminous infrared galaxies: a deep dive beneath the mass-metallicity relation

Authors: Borja Pérez-Díaz, Enrique Pérez-Montero, Juan A. Fernández-Ontiveros, José M. Vílchez, Ricardo Amorín

Abstract: The build up of heavy elements and the stellar mass assembly are fundamental processes in the formation and evolution of galaxies. Although they have been extensively studied through observations and simulations, the key elements that govern these processes, such as gas accretion and outflows, are not fully understood. This is especially true for luminous and massive galaxies, which usually suffer… ▽ More The build up of heavy elements and the stellar mass assembly are fundamental processes in the formation and evolution of galaxies. Although they have been extensively studied through observations and simulations, the key elements that govern these processes, such as gas accretion and outflows, are not fully understood. This is especially true for luminous and massive galaxies, which usually suffer strong feedback in the form of massive outflows, and large-scale gas accretion triggered by galaxy interactions. For a sample of 77 luminous infrared (IR) galaxies, we derive chemical abundances using new diagnostics based on nebular IR lines, which peer through the dusty medium of these objects and allow us to include the obscured metals in our abundance determinations. In contrast to optical-based studies, our analysis reveals that most luminous IR galaxies remain close to the mass-metallicity relation. Nevertheless, four galaxies with extreme star-formation rates ($> 60$M$_{\odot }$yr$^{-1}$) in their late merger stages show heavily depressed metallicities of 12+log(O/H) $\sim 7.7$--$8.1$ along with solar-like N/O ratios, indicative of gas mixing processes affecting their chemical composition. This evidence suggests the action of a massive infall of metal-poor gas in a short phase during the late merger stages, eventually followed by a rapid enrichment. These results challenge the classical gas equilibrium scenario usually applied to main-sequence galaxies, suggesting that the chemical enrichment and stellar-mass growth in luminous IR galaxies are regulated by different processes. △ Less

Submitted 26 June, 2023; originally announced June 2023.

Comments: Submitted to Nature Astronomy. 19 pages, 4 figures

The miniJPAS survey quasar selection III: Classification with artificial neural networks and hybridisation

Authors: G. Martínez-Solaeche, Carolina Queiroz, R. M. González Delgado, Natália V. N. Rodrigues, R. García-Benito, Ignasi Pérez-Ràfols, L. Raul Abramo, Luis Díaz-García, Matthew M. Pieri, Jonás Chaves-Montero, A. Hernán-Caballero, J. E. Rodríguez-Martín, Silvia Bonoli, Sean S. Morrison, Isabel Márquez, J. M. Vílchez, C. López-Sanjuan, A. J. Cenarro, R. A. Dupke, A. Martín-Franch, J. Varel, H. Vázquez Ramió, D. Cristóbal-Hornillos, M. Moles, J. Alcaniz , et al. (6 additional authors not shown)

Abstract: This paper is part of large effort within the J-PAS collaboration that aims to classify point-like sources in miniJPAS, which were observed in 60 optical bands over $\sim$ 1 deg$^2$ in the AEGIS field. We developed two algorithms based on artificial neural networks (ANN) to classify objects into four categories: stars, galaxies, quasars at low redshift ($z < 2.1)$, and quasars at high redshift (… ▽ More This paper is part of large effort within the J-PAS collaboration that aims to classify point-like sources in miniJPAS, which were observed in 60 optical bands over $\sim$ 1 deg$^2$ in the AEGIS field. We developed two algorithms based on artificial neural networks (ANN) to classify objects into four categories: stars, galaxies, quasars at low redshift ($z < 2.1)$, and quasars at high redshift ($z \geq 2.1$). As inputs, we used miniJPAS fluxes for one of the classifiers (ANN$_1$) and colours for the other (ANN$_2$). The ANNs were trained and tested using mock data in the first place. We studied the effect of augmenting the training set by creating hybrid objects, which combines fluxes from stars, galaxies, and quasars. Nevertheless, the augmentation processing did not improve the score of the ANN. We also evaluated the performance of the classifiers in a small subset of the SDSS DR12Q superset observed by miniJPAS. In the mock test set, the f1-score for quasars at high redshift with the ANN$_1$ (ANN$_2$) are $0.99$ ($0.99$), $0.93$ ($0.92$), and $0.63$ ($0.57$) for $17 < r \leq 20$, $20 < r \leq 22.5$, and $22.5 < r \leq 23.6$, respectively, where $r$ is the J-PAS rSDSS band. In the case of low-redshift quasars, galaxies, and stars, we reached $0.97$ ($0.97$), $0.82$ ($0.79$), and $0.61$ ($0.58$); $0.94$ ($0.94$), $0.90$ ($0.89$), and $0.81$ ($0.80$); and $1.0$ ($1.0$), $0.96$ ($0.94$), and $0.70$ ($0.52$) in the same r bins. In the SDSS DR12Q superset miniJPAS sample, the weighted f1-score reaches 0.87 (0.88) for objects that are mostly within $20 < r \leq 22.5$. Finally, we estimate the number of point-like sources that are quasars, galaxies, and stars in miniJPAS. △ Less

Submitted 22 March, 2023; originally announced March 2023.

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

The miniJPAS survey: AGN & host galaxy co-evolution of X-ray selected sources

Authors: I. E. López, M. Brusa, S. Bonoli, F. Shankar, N. Acharya, B. Laloux, K. Dolag, A. Georgakakis, A. Lapi, C. Ramos Almeida, M. Salvato, J. Chaves-Montero, P. Coelho, L. A. Díaz-García, J. A. Fernández-Ontiveros, A. Hernán-Caballero, R. M. González Delgado, I. Marquez, M. Pović, R. Soria, C. Queiroz, P. T. Rahna, R. Abramo, J. Alcaniz, N. Benitez , et al. (13 additional authors not shown)

Abstract: Studies indicate strong evidence of a scaling relation in the local Universe between the supermassive black hole mass ($M_\rm{BH}$) and the stellar mass of their host galaxies ($M_\star$). They even show similar histories across cosmic times of their differential terms: star formation rate (SFR) and black hole accretion rate (BHAR). However, a clear picture of this coevolution is far from being un… ▽ More Studies indicate strong evidence of a scaling relation in the local Universe between the supermassive black hole mass ($M_\rm{BH}$) and the stellar mass of their host galaxies ($M_\star$). They even show similar histories across cosmic times of their differential terms: star formation rate (SFR) and black hole accretion rate (BHAR). However, a clear picture of this coevolution is far from being understood. We select an X-ray sample of active galactic nuclei (AGN) up to $z=2.5$ in the miniJPAS footprint. Their X-ray to infrared spectral energy distributions (SEDs) have been modeled with CIGALE, constraining the emission to 68 bands. For a final sample of 308 galaxies, we derive their physical properties (e.g., $M_\star$, $\rm{SFR}$, $\rm{SFH}$, and $L_\rm{AGN}$). We also fit their optical spectra for a subsample of 113 sources to estimate the $M_\rm{BH}$. We calculate the BHAR depending on two radiative efficiency regimes. We find that the Eddington ratios ($λ$) and its popular proxy ($L_\rm{X}$/$M_\star$) have 0.6 dex of difference, and a KS-test indicates that they come from different distributions. Our sources exhibit a considerable scatter on the $M_\rm{BH}$-$M_\star$ relation, which can explain the difference between $λ$ and its proxy. We also model three evolution scenarios to recover the integral properties at $z=0$. Using the SFR and BHAR, we show a notable diminution in the scattering between $M_\rm{BH}$-$M_\star$. For the last scenario, we consider the SFH and a simple energy budget for the AGN accretion, obtaining a relation similar to the local Universe. Our study covers $\sim 1$ deg$^2$ in the sky and is sensitive to biases in luminosity. Nevertheless, we show that, for bright sources, the link between SFR and BHAR, and their decoupling based on an energy limit is the key that leads to the local $M_\rm{BH}$-$M_\star$ scaling relation. △ Less

Submitted 2 February, 2023; originally announced February 2023.

Comments: Accepted for publication in A&A

Journal ref: A&A 672, A137 (2023)

J-PLUS: Towards an homogeneous photometric calibration using Gaia BP/RP low-resolution spectra

Authors: C. López-Sanjuan, H. Vázquez Ramió, K. Xiao, H. Yuan, J. M. Carrasco, J. Varela, D. Cristóbal-Hornillos, P. -E. Tremblay, A. Ederoclite, A. Marín-Franch, A. J. Cenarro, P. R. T. Coelho, S. Daflon, A. del Pino, H. Domínguez Sánchez, J. A. Fernández-Ontiveros, A. Hernán-Caballero, F. M. Jiménez-Esteban, J. Alcaniz, R. E. Angulo, R. A. Dupke, C. Hernández-Monteagudo, M. Moles, L. Sodré Jr

Abstract: We present the photometric calibration of the twelve optical passbands for the Javalambre Photometric Local Universe Survey (J-PLUS) third data release (DR3), comprising 1642 pointings of two square degrees each. We selected nearly 1.5 million main sequence stars with a signal-to-noise ratio larger than ten in the twelve J-PLUS passbands and available low-resolution (R = 20-80) spectrum from the b… ▽ More We present the photometric calibration of the twelve optical passbands for the Javalambre Photometric Local Universe Survey (J-PLUS) third data release (DR3), comprising 1642 pointings of two square degrees each. We selected nearly 1.5 million main sequence stars with a signal-to-noise ratio larger than ten in the twelve J-PLUS passbands and available low-resolution (R = 20-80) spectrum from the blue and red photometers (BP/RP) in Gaia DR3. We compared the synthetic photometry from BP/RP spectra with the J-PLUS instrumental magnitudes, after correcting for the magnitude and color terms between both systems, to obtain an homogeneous photometric solution for J-PLUS. To circumvent the current limitations in the absolute calibration of the BP/RP spectra, the absolute color scale was derived using the locus of 109 white dwarfs closer than 100 pc with a negligible interstellar extinction. Finally, the absolute flux scale was anchored to the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS) photometry in the r band. The precision of the J-PLUS photometric calibration, estimated from duplicated objects observed in adjacent pointings and by comparison with the spectro-photometric standard star GD 153, is ~12 mmag in u, J0378, and J0395; and ~7 mmag in J0410, J0430, g, J0515, r, J0660, i, J0861, and z. The estimated accuracy in the calibration along the surveyed area is better than 1% for all the passbands. The Gaia BP/RP spectra provide a high-quality, homogeneous photometric reference in the optical range across the full-sky, in spite of their current limitations as an absolute reference. The calibration method for J-PLUS DR3 reaches an absolute precision and accuracy of 1% in the twelve optical filters within an area of 3284 square degrees. △ Less

Submitted 29 January, 2023; originally announced January 2023.

Comments: Submitted to Astronomy and Astrophysics. 17 pages, 14 figures, 2 appendix. Comments are welcome. arXiv admin note: text overlap with arXiv:2101.12407

J-NEP: 60-band photometry and photometric redshifts for the James Webb Space Telescope North Ecliptic Pole Time-Domain Field

Authors: A. Hernán-Caballero, C. N. A. Willmer, J. Varela, C. López-Sanjuan, A. Marín-Franch, H. Vázquez Ramió, T. Civera, A. Ederoclite, D. Muniesa, J. Cenarro, S. Bonoli, R. Dupke, J. Lim, J. Chaves-Montero, J. Laur, C. Hernández-Monteagudo, J. A. Fernández-Ontiveros, A. Fernández-Soto, L. A. Díaz-García, R. M. González Delgado, C. Queiroz, J. M. Vílchez, R. Abramo, J. Alcaniz, N. Benítez , et al. (6 additional authors not shown)

Abstract: The J-PAS survey will observe ~1/3 of the northern sky with a set of 56 narrow-band filters using the dedicated 2.55 m JST telescope at the Javalambre Astrophysical Observatory. Prior to the installation of the main camera, in order to demonstrate the scientific potential of J-PAS, two small surveys were performed with the single-CCD Pathfinder camera: miniJPAS (~1 deg2 along the Extended Groth St… ▽ More The J-PAS survey will observe ~1/3 of the northern sky with a set of 56 narrow-band filters using the dedicated 2.55 m JST telescope at the Javalambre Astrophysical Observatory. Prior to the installation of the main camera, in order to demonstrate the scientific potential of J-PAS, two small surveys were performed with the single-CCD Pathfinder camera: miniJPAS (~1 deg2 along the Extended Groth Strip), and J-NEP (~0.3 deg2 around the JWST North Ecliptic Pole Time Domain Field), including all 56 J-PAS filters as well as u, g, r, and i. J-NEP is ~0.5-1.0 magnitudes deeper than miniJPAS, providing photometry for 24,618 r-band detected sources and photometric redshifts (photo-z) for the 6,662 sources with r<23. In this paper we describe the photometry and photo-z of J-NEP and demonstrate a new method for the removal of systematic offsets in the photometry based on the median colours of galaxies, dubbed "galaxy locus recalibration". This method does not require spectroscopic observations except in a few reference pointings and, unlike previous methods, is applicable to the whole J-PAS survey. We use a spectroscopic sample of 787 galaxies to test the photo-z performance for J-NEP and in comparison to miniJPAS. We find that the deeper J-NEP observations result in a factor ~1.5-2 decrease in sigma_NMAD (a robust estimate of the standard deviation of the photo-z error) and the outlier rate relative to miniJPAS for r>21.5 sources, but no improvement in brighter ones. We find the same relation between sigma_NMAD and odds in J-NEP and miniJPAS, suggesting sigma_NMAD can be predicted for any set of J-PAS sources from their odds distribution alone, with no need for additional spectroscopy to calibrate the relation. We explore the causes for photo-z outliers and find that colour-space degeneracy at low S/N, photometry artifacts, source blending, and exotic spectra are the most important factors. △ Less

Submitted 23 January, 2023; originally announced January 2023.

Comments: 16 pages, 25 figures, accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 671, A71 (2023)

J-PLUS DR3: Galaxy-Star-Quasar classification

Authors: R. von Marttens, V. Marra, M. Quartin, L. Casarini, P. O. Baqui, A. Alvarez-Candal, F. J. Galindo-Guil, J. A. Fernández-Ontiveros, Andrés del Pino, L. A. Díaz-García, C. López-Sanjuan, J. Alcaniz, R. Angulo, A. J. Cenarro, D. Cristóbal-Hornillos, R. Dupke, A. Ederoclite, C. Hernández-Monteagudo, A. Marín-Franch, M. Moles, L. Sodré, J. Varela, H. Vázquez Ramió

Abstract: The Javalambre Photometric Local Universe Survey (J-PLUS) is a 12-band photometric survey using the 83-cm JAST telescope. Data Release 3 includes 47.4 million sources. J-PLUS DR3 only provides star-galaxy classification so that quasars are not identified from the other sources. Given the size of the dataset, machine learning methods could provide a valid alternative classification and a solution t… ▽ More The Javalambre Photometric Local Universe Survey (J-PLUS) is a 12-band photometric survey using the 83-cm JAST telescope. Data Release 3 includes 47.4 million sources. J-PLUS DR3 only provides star-galaxy classification so that quasars are not identified from the other sources. Given the size of the dataset, machine learning methods could provide a valid alternative classification and a solution to the classification of quasars. Our objective is to classify J-PLUS DR3 sources into galaxies, stars and quasars, outperforming the available classifiers in each class. We use an automated machine learning tool called TPOT to find an optimized pipeline to perform the classification. The supervised machine learning algorithms are trained on the crossmatch with SDSS DR18, LAMOST DR8 and Gaia. We checked that the training set of about 660 thousand galaxies, 1.2 million stars and 270 thousand quasars is both representative and contain a minimal presence of contaminants (less than 1%). We considered 37 features: the twelve photometric bands with respective errors, six colors, four morphological parameters, galactic extinction with its error and the PSF relative to the corresponding pointing. With TPOT genetic algorithm, we found that XGBoost provides the best performance: the AUC for galaxies, stars and quasars is above 0.99 and the average precision is above 0.99 for galaxies and stars and 0.96 for quasars. XGBoost outperforms the classifiers already provided in J-PLUS DR3 and also classifies quasars. △ Less

Submitted 2 November, 2023; v1 submitted 12 December, 2022; originally announced December 2022.

Comments: 19 pages, 22 figures

Journal ref: MNRAS 527 (2), 3347-3365 (2024)

Costs and benefits of automation for astronomical facilities

Authors: A. Yanes-Díaz, S. Rueda-Teruel, R. Bello, D. Lozano-Pérez, M. Royo-Navarro, T. Civera, M. Domínguez-Martínez, N. Martínez-Olivar, S. Chueca, C. Iñiguez, A. Marin-Franch, F. Rueda-Teruel, G. López-Alegre, S. Bielsa, J. Muñoz-Maudos, H. Rueda-Asensio, A. Muñoz-Teruel, D. Garcés-Cubel, I. Soriano-Laguía, M. Almarcegui-Gracia, A. J. Cenarro, M. Moles, D. Cristobal-Hornillos, J. Varela, A. Ederoclite , et al. (35 additional authors not shown)

Abstract: The Observatorio Astrofísico de Javalambre (OAJ†1) in Spain is a young astronomical facility, conceived and developed from the beginning as a fully automated observatory with the main goal of optimizing the processes in the scientific and general operation of the Observatory. The OAJ has been particularly conceived for carrying out large sky surveys with two unprecedented telescopes of unusually l… ▽ More The Observatorio Astrofísico de Javalambre (OAJ†1) in Spain is a young astronomical facility, conceived and developed from the beginning as a fully automated observatory with the main goal of optimizing the processes in the scientific and general operation of the Observatory. The OAJ has been particularly conceived for carrying out large sky surveys with two unprecedented telescopes of unusually large fields of view (FoV): the JST/T250, a 2.55m telescope of 3deg field of view, and the JAST/T80, an 83cm telescope of 2deg field of view. The most immediate objective of the two telescopes for the next years is carrying out two unique photometric surveys of several thousands square degrees, J-PAS†2 and J-PLUS†3, each of them with a wide range of scientific applications, like e.g. large structure cosmology and Dark Energy, galaxy evolution, supernovae, Milky Way structure, exoplanets, among many others. To do that, JST and JAST are equipped with panoramic cameras under development within the J-PAS collaboration, JPCam and T80Cam respectively, which make use of large format (~ 10k x 10k) CCDs covering the entire focal plane. This paper describes in detail, from operations point of view, a comparison between the detailed cost of the global automation of the Observatory and the standard automation cost for astronomical facilities, in reference to the total investment and highlighting all benefits obtained from this approach and difficulties encountered. The paper also describes the engineering development of the overall facilities and infrastructures for the fully automated observatory and a global overview of current status, pinpointing lessons learned in order to boost observatory operations performance, achieving scientific targets, maintaining quality requirements, but also minimizing operation cost and human resources. △ Less

Submitted 24 November, 2022; v1 submitted 23 November, 2022; originally announced November 2022.

Comments: Global Observatory Control System GOCS

Journal ref: Proceedings Volume 12186, Observatory Operations: Strategies, Processes, and Systems IX; 121860Y (2022)

Compact jets dominate the continuum emission in low-luminosity active galactic nuclei

Authors: J. A. Fernández-Ontiveros, X. López-López, A. Prieto

Abstract: The disappearance of the accretion disc in low-luminosity active galactic nuclei (LLAGN) leaves behind a faint optical nuclear continuum whose nature has been largely debated, mainly due to serious observational limitations in the IR to UV range. We combine multi-wavelength sub-arcsecond resolution observations -- able to isolate the genuine nuclear continuum -- with nebular lines in the mid-IR, t… ▽ More The disappearance of the accretion disc in low-luminosity active galactic nuclei (LLAGN) leaves behind a faint optical nuclear continuum whose nature has been largely debated, mainly due to serious observational limitations in the IR to UV range. We combine multi-wavelength sub-arcsecond resolution observations -- able to isolate the genuine nuclear continuum -- with nebular lines in the mid-IR, to indirectly probe the shape of the extreme UV continuum. We found that 8 of the nearest prototype LLAGN are compatible with pure compact jet emission (self-absorbed synchrotron plus the associated self-Compton component) over more than ten orders of magnitude in frequency. When compared with typical radio galaxies, the LLAGN continua show two peculiarities: $i)$ a very steep spectral slope in the IR-to-optical/UV range ($-3.7 < α_0 < -1.3$; $F_ν\propto ν^{α_0}$); and $ii)$ a very high turnover frequency ($0.2-30\, \rm{THz}$; $1.3\,\rm{mm}-10\,\rm{μm}$). These attributes can be explained if the synchrotron continuum is mainly dominated by thermalised particles at the jet base or corona with considerably high temperatures, whereas only a small fraction of the energy ($\sim 20\%$) would be distributed along the high-energy power-law tail of accelerated particles. On the other hand, the nebular gas excitation in LLAGN is in agreement with photo-ionisation from inverse Compton radiation ($α_{\rm x} \sim -0.7$), which would dominate the nuclear continuum shortwards of $\sim 3000$ Å. Our results suggest that the LLAGN continuum can be dominated at all wavelengths by undeveloped jets, powered by a thermalised particle distribution, similar to the behaviour observed in compact jets of quiescent black hole X-ray binaries. This has important implications in the context of galaxy evolution, since LLAGN may represent a major but underestimated source of kinetic feedback in galaxies. △ Less

Submitted 17 November, 2022; originally announced November 2022.

Comments: Accepted for publication in A&A. Abstract abridged to meet arXiv requirements

Journal ref: A&A 670, A22 (2023)

The high-ionization IR fine-structure lines as bolometric indicators of the AGN power: study of the complete 12 micron AGN sample

Authors: Luigi Spinoglio, Juan Antonio Fernandez-Ontiveros, Matthew A. Malkan

Abstract: The high-ionization mid-IR lines, excited in the Narrow Line Regions (NLR) of Active Galactic Nuclei (AGN), barely affected by stellar excitation and dust extinction, trace the AGN bolometric power. We used the complete 12 micron sample of Seyfert galaxies, for which 100/116 objects have reliable 2-10keV observations. The [NeV] and [OIV] mid-IR lines linearly correlate with several AGN bolometric… ▽ More The high-ionization mid-IR lines, excited in the Narrow Line Regions (NLR) of Active Galactic Nuclei (AGN), barely affected by stellar excitation and dust extinction, trace the AGN bolometric power. We used the complete 12 micron sample of Seyfert galaxies, for which 100/116 objects have reliable 2-10keV observations. The [NeV] and [OIV] mid-IR lines linearly correlate with several AGN bolometric indicators (intrinsic 2-10keV and observed 14-195keV X-ray emission, compact nuclear 12 micron emission, [OIII] 5007A line emission), both in terms of flux and luminosity. No evidence of systematic differences in these correlations is found among the Seyfert populations, including type 1 and type 2, and Compton thick and thin AGN. Nevertheless, we find that a sequence of high-to-low Eddington ratio together with strong-to-weak line excitation (traced by the [OIV]/[Ne II] line ratio) encompasses from type 1 through type 2 AGN to Low Ionization Nuclear Emission Line Region (LINER) galaxies, showing intrinsic differences in these three AGN populations. A positive correlation between the Black Hole Accretion Rate (BHAR) and the Star Formation Rate (SFR) is found, but no correlation between the specific SFR (sSFR) and the ratio BHAR/M_(BH), simply reflecting the fact that the more massive is a galaxy, the more it is forming stars and feeding its central black hole. The JWST telescope, just beginning operations, will allow large samples of AGN to be observed in these lines in the nearby Universe (z<0.9). △ Less

Submitted 5 October, 2022; originally announced October 2022.

Comments: To be published in The Astrophysical Journal

TOPz: Photometric redshifts for J-PAS

Authors: J. Laur, E. Tempel, A. Tamm, R. Kipper, L. J. Liivamägi, A. Hernán-Caballero, M. M. Muru, J. Chaves-Montero, L. A. Díaz-García, S. Turner, T. Tuvikene, C. Queiroz, C. R. Bom, J. A. Fernández-Ontiveros, R. M. González Delgado, T. Civera, R. Abramo, J. Alcaniz, N. Benitez, S. Bonoli, S. Carneiro, J. Cenarro, D. Cristóbal-Hornillos, R. Dupke, A. Ederoclite , et al. (8 additional authors not shown)

Abstract: The importance of photometric galaxy redshift estimation is rapidly increasing with the development of specialised powerful observational facilities. We develop a new photometric redshift estimation workflow TOPz to provide reliable and efficient redshift estimations for the upcoming large-scale survey J-PAS which will observe 8500 deg2 of the northern sky through 54 narrow-band filters. TOPz reli… ▽ More The importance of photometric galaxy redshift estimation is rapidly increasing with the development of specialised powerful observational facilities. We develop a new photometric redshift estimation workflow TOPz to provide reliable and efficient redshift estimations for the upcoming large-scale survey J-PAS which will observe 8500 deg2 of the northern sky through 54 narrow-band filters. TOPz relies on template-based photo-z estimation with some added J-PAS specific features and possibilities. We present TOPz performance on data from the miniJPAS survey, a precursor to the J-PAS survey with an identical filter system. First, we generated spectral templates based on the miniJPAS sources using the synthetic galaxy spectrum generation software CIGALE. Then we applied corrections to the input photometry by minimising systematic offsets from the template flux in each filter. To assess the accuracy of the redshift estimation, we used spectroscopic redshifts from the DEEP2, DEEP3, and SDSS surveys, available for 1989 miniJPAS galaxies with r < 22 magAB. We also tested how the choice and number of input templates, photo-z priors, and photometric corrections affect the TOPz redshift accuracy. The general performance of the combination of miniJPAS data and the TOPz workflow fulfills the expectations for J-PAS redshift accuracy. Similarly to previous estimates, we find that 38.6% of galaxies with r < 22 mag reach the J-PAS redshift accuracy goal of dz/(1 + z) < 0.003. Limiting the number of spectra in the template set improves the redshift accuracy up to 5%, especially for fainter, noise-dominated sources. Further improvements will be possible once the actual J-PAS data become available. △ Less

Submitted 2 September, 2022; originally announced September 2022.

Comments: 20 pages, 24 figures

Journal ref: A&A 668, A8 (2022)

Measuring chemical abundances in AGN from infrared nebular lines: HII-CHI-Mistry-IR for AGN

Authors: Borja Pérez-Díaz, Enrique Pérez-Montero, Juan A. Fernández-Ontiveros, José M. Vílchez

Abstract: Future and on-going infrared and radio observatories such as JWST, METIS or ALMA will increase the amount of rest-frame IR spectroscopic data for galaxies by several orders of magnitude. While studies of the chemical composition of the ISM based on optical observations have been widely spread over decades for SFG and, more recently, for AGN, similar studies need to be performed using IR data. This… ▽ More Future and on-going infrared and radio observatories such as JWST, METIS or ALMA will increase the amount of rest-frame IR spectroscopic data for galaxies by several orders of magnitude. While studies of the chemical composition of the ISM based on optical observations have been widely spread over decades for SFG and, more recently, for AGN, similar studies need to be performed using IR data. This regime can be especially useful in the case of AGN given that it is less affected by temperature and dust extinction, traces higher ionic species and can also provide robust estimations of the chemical abundance ratio N/O. We present a new tool based on a bayesian-like methodology to estimate chemical abundances from IR emission lines in AGN. We use a sample of 58 AGN with IR spectroscopic data retrieved from the literature to probe the validity of our method. The estimations of the chemical abundances based on IR lines in our sample are later compared with the corresponding abundances derived from the optical emission lines in the same objects. HII-CHI-Mistry-IR takes advantage of photoionization models, characterized by the chemical abundance ratios O/H and N/O and the ionization parameter $U$, to compare their predicted emission-line fluxes with a set of observed values. Instead of matching single emission lines, the code uses some specific emission-line ratios sensitive to the above free parameters. We report mainly solar and also subsolar abundances for O/H in the nuclear region for our sample of AGN, whereas N/O clusters around solar values. We find a discrepancy between the chemical abundances derived from IR and optical emission lines, being the latter higher than the former. This discrepancy, also reported by previous studies of the composition of the ISM in AGN from IR observations, is independent from the gas density or the incident radiation field to the gas. △ Less

Submitted 18 July, 2022; originally announced July 2022.

Comments: 20 pages, 15 figures, 8 tables. Accepted for publication in A&A

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

The miniJPAS survey: The role of group environment in quenching the star formation

Authors: R. M. González Delgado, J. E. Rodríguez-Martín, L. A. Díaz-García, A. de Amorim, R. García-Benito, G. Martínez-Solaeche, P. A. A. Lopes, M. Maturi, E. Pérez, R. Cid Fernandes, A. Cortesi, A. Finoguenov, E. R. Carrasco, A. Hernán-Caballero, L. R. Abramo, J. Alcaniz, N. Benítez, S. Bonoli, A. J. Cenarro, D. Cristóbal-Hornillos, J. M. Diego, R. A. Dupke, A. Ederoclite, J. A. Fernández-Ontiveros, C. López-Sanjuan , et al. (10 additional authors not shown)

Abstract: The miniJPAS survey has observed $\sim 1$ deg$^2$ on the AEGIS field with 60 bands (spectral resolution of $R \sim 60$) in order to demonstrate the capabilities of the Javalambre-Physics of the Accelerating Universe Astrophysical Survey (J-PAS) that will map $\sim 8000$ deg$^2$ of the northern sky in the next years. This paper shows the power of J-PAS to detect low mass groups and characterise the… ▽ More The miniJPAS survey has observed $\sim 1$ deg$^2$ on the AEGIS field with 60 bands (spectral resolution of $R \sim 60$) in order to demonstrate the capabilities of the Javalambre-Physics of the Accelerating Universe Astrophysical Survey (J-PAS) that will map $\sim 8000$ deg$^2$ of the northern sky in the next years. This paper shows the power of J-PAS to detect low mass groups and characterise their galaxy populations up to $z \sim 1$. We use the spectral energy distribution fitting code BaySeAGal to derive the stellar population properties of the galaxy members in 80 groups at $z \leq 0.8$ previously detected by the AMICO code, as well as for a galaxy field sample retrieved from the whole miniJPAS sample. We identify blue, red, quiescent, and transition galaxy populations through their rest-frame (extinction corrected) colour, stellar mass ($M_\star$) and specific star formation rate. We measure their abundance as a function of $M_\star$ and environment. We find: (i) The fraction of red and quiescent galaxies in groups increases with $M_\star$ and it is always higher in groups than in the field. (ii) The quenched fraction excess (QFE) in groups strongly increases with $M_\star$, (from a few percent to higher than 60% in the mass range $10^{10} - 3 \times 10 ^{11}$ $M_\odot$. (iii) The abundance excess of transition galaxies in groups shows a modest dependence with $M_\star$ (iv) The fading time scale is very short ($<1.5$ Gyr), indicating that the star formation declines very rapidly in groups. (v) The evolution of the galaxy quenching rate in groups shows a modest but significant evolution since $z\sim0.8$, compatible with an evolution with constant $QFE=0.4$, previously measured for satellites in the nearby Universe, and consistent with a scenario where the low-mass star-forming galaxies in clusters at $z= 1-1.4$ are environmentally quenched. △ Less

Submitted 12 July, 2022; originally announced July 2022.

Comments: 21 pages, 20 figures, accepted for publication in Astronomy & Astrophysics

Linking pressure gradients with the stability of molecular clouds in galactic outflows

Authors: K. M. Dasyra, G. F. Paraschos, T. Bisbas, F. Combes, J. A. Fernandez-Ontiveros

Abstract: The jets launched by actively accreting black holes are capable of launching several of the massive (million or billion solar mass) molecular outflows observed in galaxies. These outflows could suppress or enhance star formation in galaxies. To investigate the stability of clouds capable to form stars in outflows, we modeled CO and HCO+ ALMA data of the galaxy IC5063, in which black-hole jets impa… ▽ More The jets launched by actively accreting black holes are capable of launching several of the massive (million or billion solar mass) molecular outflows observed in galaxies. These outflows could suppress or enhance star formation in galaxies. To investigate the stability of clouds capable to form stars in outflows, we modeled CO and HCO+ ALMA data of the galaxy IC5063, in which black-hole jets impact molecular clouds. Using a radiative transfer code that self-consistently performs astrochemical and thermal balance calculations based on the available gas heating sources, we found that mechanical heating and cosmic ray (CR) heating are fully capable of individually reproducing the data. In our best-fit model, CRs provide about 1/3rd of the dense gas heating at the radio lobes, emphasizing the role of this often neglected mechanism in heating the gas and potentially generating outflows. The gas temperature and density indicate that the jet passage leads to an increase of about 1 order of magnitude in the internal pressure Pi of molecular clouds (with Pi/k from 8*10^5 up to 7*10^6 K cm^-3), irrespective of the excitation mechanism. From the fluxes of [S II] and [N II] lines in VLT MUSE data, the external pressure Pe of molecular clouds increases in several regions enough to exceed Pi. This result leads us to conclude that we are observing the expansion of an ionized overpressurized cocoon that compresses molecular clouds and that could lead to their collapse. Some jet-impacted clouds, nonetheless, near pathways that the jet cleared have increased Pi and decreased Pe. They are likely to undergo evaporation of their outer layers. Part of the evaporated layers could mass load the outflow thanks to ram pressure from co-spatial ionized gas flows. The observed pressure changes thus suggest that both star formation enhancement and suppression could simultaneously occur. △ Less

Submitted 11 May, 2022; originally announced May 2022.

Comments: Submitted to Nature Astronomy

J-PLUS: Support Vector Regression to Measure Stellar Parameters

Authors: Cunshi Wang, Yu Bai, Haibo Yuan, Jifeng Liu, J. A. Fernández-Ontiveros, Paula R. T. Coelho, F. Jiménez-Esteban, Carlos Andrés Galarza, R. E. Angulo, A. J. Cenarro, D. Cristóbal-Hornillos, R. A. Dupke, A. Ederoclite, C. Hernández-Monteagudo, C. López-Sanjuan, A. Marín-Franch, M. Moles, L. Sodré Jr., H. Vázquez Ramió, J. Varela

Abstract: Context. Stellar parameters are among the most important characteristics in studies of stars, which are based on atmosphere models in traditional methods. However, time cost and brightness limits restrain the efficiency of spectral observations. The J-PLUS is an observational campaign that aims to obtain photometry in 12 bands. Owing to its characteristics, J-PLUS data have become a valuable resou… ▽ More Context. Stellar parameters are among the most important characteristics in studies of stars, which are based on atmosphere models in traditional methods. However, time cost and brightness limits restrain the efficiency of spectral observations. The J-PLUS is an observational campaign that aims to obtain photometry in 12 bands. Owing to its characteristics, J-PLUS data have become a valuable resource for studies of stars. Machine learning provides powerful tools to efficiently analyse large data sets, such as the one from J-PLUS, and enable us to expand the research domain to stellar parameters. Aims. The main goal of this study is to construct a SVR algorithm to estimate stellar parameters of the stars in the first data release of the J-PLUS observational campaign. Methods. The training data for the parameters regressions is featured with 12-waveband photometry from J-PLUS, and is cross-identified with spectrum-based catalogs. These catalogs are from the LAMOST, the APOGEE, and the SEGUE. We then label them with the stellar effective temperature, the surface gravity and the metallicity. Ten percent of the sample is held out to apply a blind test. We develop a new method, a multi-model approach in order to fully take into account the uncertainties of both the magnitudes and stellar parameters. The method utilizes more than two hundred models to apply the uncertainty analysis. Results. We present a catalog of 2,493,424 stars with the Root Mean Square Error of 160K in the effective temperature regression, 0.35 in the surface gravity regression and 0.25 in the metallicity regression. We also discuss the advantages of this multi-model approach and compare it to other machine-learning methods. △ Less

Submitted 15 August, 2022; v1 submitted 5 May, 2022; originally announced May 2022.

Comments: 18 pages, 10 figures, accepted by AA

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

The miniJPAS survey: Identification and characterization of the emission line galaxies down to $z < 0.35$ in the AEGIS field

Authors: G. Martínez-Solaeche, R. M. González Delgado, R. García-Benito, L. A. Díaz-García, J. E. Rodríguez-Martín, E. Pérez, A. de Amorim, S. Duarte Puertas, Laerte Sodré Jr., David Sobral, Jonás Chaves-Montero, J. M. Vílchez, A. Hernán-Caballero, C. López-Sanjuan, A. Cortesi, S. Bonoli, A. J. Cenarro, R. A. Dupke, A. Marín-Franch, J. Varela, H. Vázquez Ramió, L. R. Abramo, D. Cristóbal-Hornillos, M. Moles, J. Alcaniz , et al. (6 additional authors not shown)

Abstract: The Javalambre-Physics of the Accelerating Universe Astrophysical Survey (J-PAS) is expected to map thousands of square degrees of the northern sky with 56 narrowband filters in the upcoming years. This will make J-PAS a very competitive and unbiased emission line survey compared to spectroscopic or narrowband surveys with fewer filters. The miniJPAS survey covered 1 deg$^2$, and it used the same… ▽ More The Javalambre-Physics of the Accelerating Universe Astrophysical Survey (J-PAS) is expected to map thousands of square degrees of the northern sky with 56 narrowband filters in the upcoming years. This will make J-PAS a very competitive and unbiased emission line survey compared to spectroscopic or narrowband surveys with fewer filters. The miniJPAS survey covered 1 deg$^2$, and it used the same photometric system as J-PAS, but the observations were carried out with the pathfinder J-PAS camera. In this work, we identify and characterize the sample of emission line galaxies (ELGs) from miniJPAS with a redshift lower than $0.35$. Using a method based on artificial neural networks, we detect the ELG population and measure the equivalent width and flux of the $Hα$, $Hβ$, [OIII], and [NII] emission lines. We explore the ionization mechanism using the diagrams [OIII]/H$β$ versus [NII]/H$α$ (BPT) and EW(H$α$) versus [NII]/H$α$ (WHAN). We identify 1787 ELGs ($83$%) from the parent sample (2154 galaxies) in the AEGIS field. For the galaxies with reliable EW values that can be placed in the WHAN diagram (2000 galaxies in total), we obtained that $72.8 \pm 0.4$%, $17.7 \pm 0.4$% , and $9.4 \pm 0.2$% are star-forming (SF), active galactic nucleus (Seyfert), and quiescent galaxies, respectively. Based on the flux of $Hα$ we find that the star formation main sequence is described as $\log$ SFR $[M_\mathrm{\odot} \mathrm{yr}^{-1}] = 0.90^{+ 0.02}_{-0.02} \log M_{\star} [M_\mathrm{\odot}] -8.85^{+ 0.19}_{-0.20}$ and has an intrinsic scatter of $0.20^{+ 0.01}_{-0.01}$. The cosmic evolution of the SFR density ($ρ_{\text{SFR}}$) is derived at three redshift bins: $0 < z \leq 0.15$, $0.15 < z \leq 0.25$, and $0.25 < z \leq 0.35$, which agrees with previous results that were based on measurements of the $Hα$ emission line. △ Less

Submitted 4 April, 2022; originally announced April 2022.

Comments: 22 pages, 19 figures

Journal ref: A&A 661, A99 (2022)

CO in the ALMA Radio-Source Catalogue (ARC): the molecular gas content of radio galaxies as a function of redshift

Authors: A. Audibert, K. M. Dasyra, M. Papachristou, J. A. Fernández-Ontiveros, I. Ruffa, L. Bisigello, F. Combes, P. Salomé, C. Gruppioni

Abstract: To evaluate the role of radio activity in galaxy evolution, we designed a large archival CO survey of radio galaxies (RGs) to determine their molecular gas masses at different epochs. We used a sample of 120 RGs representative of the NVSS 1.4 GHz survey, when flux limited at 0.4 Jy. Of those, 66 galaxies belonged to the ALMA Radio-source Catalogue (ARC) of calibrators and had spectral window tunin… ▽ More To evaluate the role of radio activity in galaxy evolution, we designed a large archival CO survey of radio galaxies (RGs) to determine their molecular gas masses at different epochs. We used a sample of 120 RGs representative of the NVSS 1.4 GHz survey, when flux limited at 0.4 Jy. Of those, 66 galaxies belonged to the ALMA Radio-source Catalogue (ARC) of calibrators and had spectral window tunings around CO (1-0), (2-1), (3-2), or (4-3). We reduced their ALMA data, determined their H2 mass contents, and combined the results with similar results for the remaining 54 galaxies from the literature. We found that, while at all epochs the majority of RGs have undetectable reservoirs, there is a rapid increase in the H2 mass content of the CO-detected RGs with z. At 1<z<2.5, one-fourth of the RGs have at least as much molecular gas as simulations would indicate for a typical halo mass of that epoch. These galaxies plausibly have ``normal'' or even starburst hosts. Taking into account the completeness correction of the sample, we created the corresponding H$_2$ mass functions at 0.005<z<0.3 and 1<z<2.5. The local mass function reveals that the number density of low-z RGs with detectable molecular gas reservoirs is only a little lower (a factor of ~4) than that of pure (or little star-forming) type 1 and 2 AGN in simulations. At 1<z<2.5, there is a significant decrease in the number density of high-z RGs due to the rarity of bright radio galaxies. An estimate for the missing faint RGs would, nonetheless, bring populations close again. Finally, we find that the volume density of molecular gas locked up in the brightest 1/5000-1/7000 RGs is similar in the examined $z$ bins. This result likely indicates that the inflow rate on one hand and the star-formation depletion rate plus the jet-driven expulsion rate on the other hand counteract each other in the most luminous RGs of each epoch. △ Less

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

Comments: 21 pages, 10 figures

Journal ref: A&A 668, A67 (2022)

The miniJPAS survey: White dwarf science with 56 optical filters

Authors: C. López-Sanjuan, P. -E. Tremblay, A. Ederoclite, H. Vázquez Ramió, A. J. Cenarro, A. Marín-Franch, J. Varela, S. Akras, M. A. Guerrero, F. M. Jiménez-Esteban, R. Lopes de Oliveira, A. L. Chies-Santos, J. A. Fernández-Ontiveros, R. Abramo, J. Alcaniz, N. Benítez, S. Bonoli, S. Carneiro, D. Cristóbal-Hornillos, R. A. Dupke, C. Mendes de Oliveira, M. Moles, L. Sodré Jr., K. Taylor

Abstract: We analyze the white dwarf population in miniJPAS, the first square degree observed with 56 medium-band, 145 A in width optical filters by the Javalambre Physics of the accelerating Universe Astrophysical Survey (J-PAS), to provide a data-based forecast for the white dwarf science with low-resolution (R ~ 50) photo-spectra. We define the sample of the bluest point-like sources in miniJPAS with r <… ▽ More We analyze the white dwarf population in miniJPAS, the first square degree observed with 56 medium-band, 145 A in width optical filters by the Javalambre Physics of the accelerating Universe Astrophysical Survey (J-PAS), to provide a data-based forecast for the white dwarf science with low-resolution (R ~ 50) photo-spectra. We define the sample of the bluest point-like sources in miniJPAS with r < 21.5 mag, point-like probability larger than 0.5, (u-r) < 0.80 mag, and (g-i) < 0.25 mag. This sample comprises 33 sources with spectroscopic information, 11 white dwarfs and 22 QSOs. We estimate the effective temperature (Teff), the surface gravity, and the composition of the white dwarf population by a Bayesian fitting to the observed photo-spectra. The miniJPAS data permit the classification of the observed white dwarfs into H-dominated and He-dominated with 99% confidence, and the detection of calcium absorption and polluting metals down to r ~ 21.5 mag at least for sources with 7000 < Teff < 22000 K, the temperature range covered by the white dwarfs in miniJPAS. The effective temperature is estimated with a 2% uncertainty, close to the 1% from spectroscopy. A precise estimation of the surface gravity depends on the available parallax information. In addition, the white dwarf population at Teff > 7000 K can be segregated from the bluest extragalactic QSOs, providing a clean sample based on optical photometry alone. The J-PAS low-resolution photo-spectra provide precise and accurate effective temperatures and atmospheric compositions for white dwarfs, complementing the data from Gaia. J-PAS will also detect and characterize new white dwarfs beyond the Gaia magnitude limit, providing faint candidates for spectroscopic follow up. △ Less

Submitted 20 March, 2022; originally announced March 2022.

Comments: Submitted to Astronomy and Astrophyscis. 12 pages, 9 figures, 4 tables. Comments are welcome

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

Galaxy evolution through infrared and submillimeter spectroscopy: Measuring star formation and black hole accretion with JWST and ALMA

Authors: Sabrina Mordini, Luigi Spinoglio, Juan Antonio Fernández-Ontiveros

Abstract: Rest-frame mid- to far-infrared spectroscopy is a powerful tool to study how galaxies formed and evolved, because a major part of their evolution occurs in heavily dust enshrouded environments, especially at the so-called Cosmic Noon. Using the calibrations of IR lines we predict the expected fluxes of lines and features, with the aim to measure the star formation rate and the Black Hole Accretion… ▽ More Rest-frame mid- to far-infrared spectroscopy is a powerful tool to study how galaxies formed and evolved, because a major part of their evolution occurs in heavily dust enshrouded environments, especially at the so-called Cosmic Noon. Using the calibrations of IR lines we predict the expected fluxes of lines and features, with the aim to measure the star formation rate and the Black Hole Accretion rate in intermediate to high redshift galaxies. The launch of the James Webb Space Telescope will allow us a deep investigation of both the SF and the BHA obscured processes as a function of cosmic time. We assess the spectral lines and features that can be detected by JWST-MIRI in galaxies and Active Galactic Nuclei up to redshift z= 3. We confirm the fine-structure lines of [MgIV]4.49um and [ArVI]4.53um as good BHA rate tracers for the 1<z<3 range, and we propose the [NeVI]7.65um line as the best tracer for redshifts of z<1.5. We suggest the use of the [ArII]6.98um and [ArIII]8.99um lines to measure the SF rate, for z<3 and z<2. At higher redshifts, the PAH features at 6.2um and 7.7um can be observed at z<3 and z<2.7 respectively. Rest-frame far-IR spectroscopy is currently being collected in high redshift galaxies (z>3) with the Atacama Large Millimeter Array. We confirm that the [CII]158um line is a good tracer of the SF rate and can in most cases (0.9<z<2 and 3<z<9) be observed, and we propose the use of the combination of [OIII]88um and [OI]145um lines as an alternative SF rate tracer, that can be detected above z>3. We conclude, however, that the current and foreseen facilities will not be able to cover properly the peak of the obscured SF and BHA activities at the Cosmic Noon of galaxy evolution and a new IR space telescope, actively cooled to obtain very good sensitivities, covering the full IR spectral range from about 10um to 300um, will be needed. △ Less

Submitted 8 March, 2022; v1 submitted 22 December, 2021; originally announced December 2021.

Comments: Accepted by PASA

J-PLUS: Stellar Parameters, C, N, Mg, Ca and [α/Fe] Abundances for Two Million Stars from DR1

Authors: Lin Yang, Haibo Yuan, Maosheng Xiang, Fuqing Duan, Yang Huang, Jifeng Liu, Timothy C. Beers, Carlos Andrés Galarza, Simone Daflon, J. A. Fernández-Ontiveros, Javier Cenarro, David Cristóbal-Hornillos, Carlos Hernández-Monteagudo, Carlos López-Sanjuan, Antonio Marín-Franch, Mariano Moles, Jesús Varela, Héctor Vázquez Ramió, Jailson Alcaniz, Renato Dupke, Alessandro Ederoclite, Laerte Sodré Jr., Raul E. Angulo

Abstract: Context. The Javalambre Photometric Local Universe Survey (J-PLUS) has obtained precise photometry in twelve specially designed filters for large numbers of Galactic stars. Deriving their precise stellar atmospheric parameters and individual elemental abundances is crucial for studies of Galactic structure, and the assembly history and chemical evolution of our Galaxy. Aims. Our goal is to estimat… ▽ More Context. The Javalambre Photometric Local Universe Survey (J-PLUS) has obtained precise photometry in twelve specially designed filters for large numbers of Galactic stars. Deriving their precise stellar atmospheric parameters and individual elemental abundances is crucial for studies of Galactic structure, and the assembly history and chemical evolution of our Galaxy. Aims. Our goal is to estimate not only stellar parameters (effective temperature, Teff, surface gravity, log g, and metallicity, [Fe/H]), but also [α/Fe] and four elemental abundances ([C/Fe], [N/Fe], [Mg/Fe], and [Ca/Fe]) using data from J-PLUS DR1. Methods. By combining recalibrated photometric data from J-PLUS DR1, Gaia DR2, and spectroscopic labels from LAMOST, we design and train a set of cost-sensitive neural networks, the CSNet, to learn the non-linear mapping from stellar colors to their labels. Results. We have achieved precisions of δTeff {\sim}55K, δlogg{\sim}0.15dex, and δ[Fe/H]{\sim}0.07dex, respectively, over a wide range of temperature, surface gravity, and metallicity. The uncertainties of the abundance estimates for [α/Fe] and the four individual elements are in the range 0.04-0.08 dex. We compare our parameter and abundance estimates with those from other spectroscopic catalogs such as APOGEE and GALAH, and find an overall good agreement. Conclusions. Our results demonstrate the potential of well-designed, high-quality photometric data for determinations of stellar parameters as well as individual elemental abundances. Applying the method to J-PLUS DR1, we have obtained the aforementioned parameters for about two million stars, providing an outstanding data set for chemo-dynamic analyses of the Milky Way. The catalog of the estimated parameters is publicly accessible. △ Less

Submitted 14 December, 2021; originally announced December 2021.

Comments: 30 pages, 26 figures, 3 tables, accepted by A&A

Journal ref: A&A 659, A181 (2022)

J-PLUS: Uncovering a large population of extreme [OIII] emitters in the local Universe

Authors: A. Lumbreras-Calle, C. López-Sanjuan, D. Sobral, J. A. Fernández-Ontiveros, J. M. Vílchez, A. Hernán-Caballero, M. Akhlaghi, L. A. Díaz-García, J. Alcaniz, R. E. Angulo, A. J. Cenarro, D. Cristóbal-Hornillos, R. A. Dupke, A. Ederoclite, C. Hernández-Monteagudo, A. Marín-Franch, M. Moles, L. Sodré Jr., H. Vázquez Ramió, J. Varela

Abstract: Over the past decades, several studies have discovered a population of galaxies undergoing very strong star formation events, called extreme emission line galaxies (EELGs). In this work, we exploit the capabilities of the Javalambre Photometric Local Universe Survey (J-PLUS), a wide field multifilter survey, with 2000 square degrees observed. We use it to identify EELGs at low redshift by their [O… ▽ More Over the past decades, several studies have discovered a population of galaxies undergoing very strong star formation events, called extreme emission line galaxies (EELGs). In this work, we exploit the capabilities of the Javalambre Photometric Local Universe Survey (J-PLUS), a wide field multifilter survey, with 2000 square degrees observed. We use it to identify EELGs at low redshift by their [OIII]5007 emission line. We intend to provide with a more complete, deep, and less biased sample of local EELGs. We select objects with an excess of flux in the J-PLUS mediumband $J0515$ filter, which covers the [OIII] line at z$<$0.06. We remove contaminants (stars and higher redshift systems) using J-PLUS and WISE infrared data, with SDSS spectra as a benchmark. We perform spectral energy distribution fitting to estimate the properties of the galaxies: line fluxes, equivalent widths (EWs), masses, etc. We identify 466 EELGs at ${\rm z} < 0.06$ with [OIII] EW over 300 \textÅ and $r$-band mag. below 20, of which 411 were previously unknown. Most show compact morphologies, low stellar masses ($\log (M_{\star}/M_{\odot}) \sim {8.13}^{+0.61}_{-0.58}$), low dust extinction ($E(B-V)\sim{0.1}^{+0.2}_{-0.1}$), and very young bursts of star formation (${3.0}^{+2.7}_{-2.0}$ Myr). Our method is up to $\sim$ 20 times more efficient detecting EELGs per Mpc$^3$ than broadband surveys, and as complete as magnitude-limited spectroscopic surveys (and reaching fainter objects). The sample is not directly biased against strong H$α$ emitters, in contrast with broadband surveys. We demonstrate the capability of J-PLUS to identify, following a clear selection process, a large sample of previously unknown EELGs showing unique properties. A fraction of them are likely similar to the first galaxies in the Universe, but at a much lower redshift, which makes them ideal targets for follow-up studies. △ Less

Submitted 13 December, 2021; originally announced December 2021.

Comments: 22 pages, 22 figures. Submitted to Astronomy & Astrophysics

Journal ref: A&A 668, A60 (2022)

Simulating infrared spectro-photometric surveys with a SPRITZ

Authors: L. Bisigello, C. Gruppioni, F. Calura, A. Feltre, F. Pozzi, C. Vignali, L. Barchiesi, G. Rodighiero, M. Negrello, F. J. Carrera, K. M. Dasyra, J. A. Fernández-Ontiveros, M. Giard, E. Hatziminaoglou, H. Kaneda, E. Lusso, M. Pereira-Santaella, P. G. Pérez-González, C. Ricci, D. Schaerer, L. Spinoglio, L. Wang

Abstract: Mid- and far-infrared (IR) photometric and spectroscopic observations are fundamental to a full understanding of the dust-obscured Universe and the evolution of both star formation and black hole accretion in galaxies. In this work, using the specifications of the SPace Infrared telescope for Cosmology and Astrophysics (SPICA) as a baseline, we investigate the capability to study the dust-obscured… ▽ More Mid- and far-infrared (IR) photometric and spectroscopic observations are fundamental to a full understanding of the dust-obscured Universe and the evolution of both star formation and black hole accretion in galaxies. In this work, using the specifications of the SPace Infrared telescope for Cosmology and Astrophysics (SPICA) as a baseline, we investigate the capability to study the dust-obscured Universe of mid- and far-IR photometry at 34 and 70 $μ$m and low-resolution spectroscopy at 17-36 $μ$m using the state-of-the-art Spectro-Photometric Realisations of Infrared-selected Targets at all-z (SPRITZ) simulation. This investigation is also compared to the expected performance of the Origins Space Telescope and the Galaxy Evolution Probe. The photometric view of the Universe of a SPICA-like mission could cover not only bright objects (e.g. L$_{IR}$>10$^{12}$L$_\odot$) up to z=10, but also normal galaxies(L$_{IR}$<10$^{11}$L$_\odot$) up to z$\sim$4. At the same time, the spectroscopic observations of such mission could also allow us to estimate the redshifts and study the physical properties for thousands of star-forming galaxies and active galactic nuclei by observing the polycyclic aromatic hydrocarbons and a large set of IR nebular emission lines. In this way, a cold, 2.5-m size space telescope with spectro-photometric capability analogous to SPICA, could provide us with a complete three-dimensional (i.e. images and integrated spectra) view of the dust-obscured Universe and the physics governing galaxy evolution up to z$\sim$4. △ Less

Submitted 22 November, 2021; originally announced November 2021.

Comments: 22 pages, 17 figures, 4 tables, accepted for publication in PASA

SOFIA observations of far-IR fine-structure lines in galaxies to measure metallicity

Authors: Luigi Spinoglio, Juan Antonio Fernandez-Ontiveros, Matthew A. Malkan, Suyash Kumar, Miguel Pereira-Santaella, Borja Perez-Diaz, Enrique Perez-Montero, Alfred Krabbe, William Vacca, Sebastian Colditz, Christian Fischer

Abstract: We present new and archival SOFIA FIFI-LS far-IR spectroscopic observations of 25 local galaxies of either the [OIII]52um and/or the [NIII]57um lines. Including other 31 galaxies from Herschel-PACS, we discuss a local sample of 47 galaxies, including HII region, luminous IR, low-metallicity dwarf and Seyfert galaxies. Analyzing the mid- to far-IR fine-structure lines of this sample, we assess the… ▽ More We present new and archival SOFIA FIFI-LS far-IR spectroscopic observations of 25 local galaxies of either the [OIII]52um and/or the [NIII]57um lines. Including other 31 galaxies from Herschel-PACS, we discuss a local sample of 47 galaxies, including HII region, luminous IR, low-metallicity dwarf and Seyfert galaxies. Analyzing the mid- to far-IR fine-structure lines of this sample, we assess the metallicity and compare with the optical spectroscopy estimates. Using the IR, we find a similar O/H--N/O relation to that known in the optical. As opposite, we find systematically lower N/O IR abundances when compared to the optical determinations, especially at high values of N/O (log(N/O) > -0.8). We explore various hypotheses to account for this difference: (i) difference in ionization structure traced by optical (O+, N+ regions) versus IR lines (O++, N++ regions); (ii) contamination of diffuse ionized gas affecting the optical lines used to compute the N/O abundance; (iii) dust obscuration affecting the optical-based determinations. However, we have not found any correlation of the Delta(N/O)= (N/O)_OPT-(N/O)_IR with either ionization, or electron density, or optical extinction. We speculatively suggest that accretion of metal-poor gas from the circumgalactic medium could provide an explanation for this difference, because the rapid decrease of total abundances during infall is followed by a N/O ratio decrease due to primary production of young - possibly embedded - massive stars, are preferentially traced by the IR diagnostics, while optical diagnostics would better trace the secondary production, when both N/O and O/H abundance ratios will increase. △ Less

Submitted 8 November, 2021; originally announced November 2021.

Comments: Accepted in the AstrophysicalJournal on the 7th November 2021

CO kinematics unveil outflows plausibly driven by a young jet in the Gigahertz Peaked Radio Core of NGC6328

Authors: Michalis Papachristou, Kalliopi Maria Dasyra, Juan Antonio Fernández-Ontiveros, Anelise Audibert, Ilaria Ruffa, Francoise Combes

Abstract: We report the detection of outflowing molecular gas in the center of the nearby (z=0.014) massive radio galaxy NGC6328. The radio core of the galaxy, PKS B1718-649, is identified as a Gigahertz Peaked Spectrum source with a compact (2 pc) double radio lobe morphology. We used ALMA CO(2-1) and CO(3-2) observations at 100 pc resolution to study the gas kinematics up to ~5 kpc from the galaxy center.… ▽ More We report the detection of outflowing molecular gas in the center of the nearby (z=0.014) massive radio galaxy NGC6328. The radio core of the galaxy, PKS B1718-649, is identified as a Gigahertz Peaked Spectrum source with a compact (2 pc) double radio lobe morphology. We used ALMA CO(2-1) and CO(3-2) observations at 100 pc resolution to study the gas kinematics up to ~5 kpc from the galaxy center. While the bulk of the molecular gas is settled in a highly warped disk, in the inner 300 pc of the disk and along with the orientation of the radio jet, we identified high-excitation and high-velocity gas that cannot be attributed to any regular kinematic component based on our detailed 3D modeling of the ALMA data. The high-velocity dispersion in the gas also suggests that it is not part of an inflowing, shredding structure. These results suggest the presence of a molecular outflow of 3 to 8 solar masses per year. The outflow possibly originated from the interaction of the jet with the dense interstellar medium, even though the radio emission is detected closer to the center than the outflow. In this sense, this source resembles NGC1377, 4C31.04 and ESO420-G13, in which the outflows are linked to faint or past jet activity. △ Less

Submitted 31 October, 2021; originally announced November 2021.

Comments: 4 pages, 4 figures, To appear in Astronomische Nachrichten, proceedings paper from "6th Workshop on Compact Steep Spectrum and GHz-Peaked Spectrum Radio Sources"

The role of SPICA-like missions and the Origins Space Telescope in the quest for heavily obscured AGN and synergies with Athena

Authors: L. Barchiesi, F. Pozzi, C. Vignali, F. J. Carrera, F. Vito, F. Calura, L. Bisigello, G. Lanzuisi, C. Gruppioni, E. Lusso, I. Delvecchio, M. Negrello, A. Cooray, A. Feltre, J. A. Fernández-Ontiveros, S. Gallerani, H. Kaneda, S. Oyabu, M. Pereira-Santaella, E. Piconcelli, C. Ricci, G. Rodighiero, L. Spinoglio, F. Tombesi

Abstract: In the BH-galaxy co-evolution framework, most of the star-formation (SF) and the black hole (BH) accretion is expected to take place in highly obscured conditions. Thus, obscured AGN are difficult to identify in optical or X-ray bands, but shine bright in the IR. Moreover, X-ray background (XRB) synthesis models predict that a large fraction of the yet-unresolved XRB is due to the most obscured (C… ▽ More In the BH-galaxy co-evolution framework, most of the star-formation (SF) and the black hole (BH) accretion is expected to take place in highly obscured conditions. Thus, obscured AGN are difficult to identify in optical or X-ray bands, but shine bright in the IR. Moreover, X-ray background (XRB) synthesis models predict that a large fraction of the yet-unresolved XRB is due to the most obscured (Compton thick, CT) of these AGN. In this work, we investigate the synergies between putative IR missions (using SPICA, proposed for ESA/M5 but withdrawn in October 2020, and Origins Space Telescope, OST, as `templates') and the X-ray mission Athena, which should fly in early 2030s, in detecting and characterizing AGN, with a particular focus on the most obscured ones. Using an XRB synthesis model, we estimated the number of AGN and the number of those which will be detected in the X-rays. For each AGN we associated an optical-to-FIR SED from observed AGN with both X-ray data and SED decomposition, and used these SEDs to check if the AGN will be detected by SPICA-like or OST at IR wavelengths. We expect that, with the deepest Athena and SPICA-like (or OST) surveys, we will be able to detect in the IR more than $90\,\%$ of all the AGN (down to L$_{2-10\text{keV}} \sim 10^{42}\,$erg/s and up to $z \sim 10$) predicted by XRB synthesis modeling, and we will detect at least half of them in the X-rays. Athena will be extremely powerful in detecting and discerning moderate- and high-luminosity AGN. We find that the most obscured and elusive CT-AGN will be exquisitely sampled by SPICA-like mission or OST and that Athena will allow a fine characterization of the most-luminous ones. This will provide a significant step forward in the process of placing stronger constraints on the yet-unresolved XRB and investigating the BH accretion rate evolution up to very high redshift ($z \ge 4$). △ Less

Submitted 15 June, 2021; originally announced June 2021.

Comments: Accepted for publication in PASA

Dust in the central parsecs of unobscured AGN: more challenges to the torus

Authors: Almudena Prieto, Jakub Nadolny, Juan A. Fernández-Ontiveros, Mar Mezcua

Abstract: A parsec-scale dusty torus is thought to be the cause of Active Galactic Nuclei (AGN) dichotomy in the 1/2 types, narrow/broad emission lines. In a previous work, on the basis of parsec-scale resolution infrared / optical dust maps it was found that dust filaments, few parsecs wide, several hundred parsecs long, were ubiquitous features crossing the centre of type 2 AGN, their optical thickness be… ▽ More A parsec-scale dusty torus is thought to be the cause of Active Galactic Nuclei (AGN) dichotomy in the 1/2 types, narrow/broad emission lines. In a previous work, on the basis of parsec-scale resolution infrared / optical dust maps it was found that dust filaments, few parsecs wide, several hundred parsecs long, were ubiquitous features crossing the centre of type 2 AGN, their optical thickness being sufficient to fully obscure the optical nucleus. This work presents the complementary view for type 1 and intermediate-type AGN. The same type of narrow, collimated, dust filaments are equally found at the centre of these AGN. The difference now resides in their location with respect to the nucleus, next to it but not crossing it, as it is the case in type 2, and their reduced optical thickness towards the centre, $A_V \lesssim 1.5 \rm{mag}$, insufficient to obscure at UV nucleus wavelengths. It is concluded that large scale, hundred pc to kpc long, dust filaments and lanes, reminiscent of those seen in the Milky Way, are a common ingredient to the central parsec of galaxies. Their optical thickness changes along their structure, in type 2 reaching optical depths high enough to obscure the nucleus in full. Their location with respect to the nucleus and increasing gradient in optical depth towards the centre could naturally lead to the canonical type 1/2 AGN classification, making these filaments to play the role of the torus. Dust filaments and lanes show equivalent morphologies in molecular gas. Available gas kinematic indicates mass inflows at rates $ < ~ 1 M\odot~ yr^{-1}$. △ Less

Submitted 16 June, 2021; v1 submitted 14 June, 2021; originally announced June 2021.

Comments: 20 pages, 14 figures, accepted in MNRAS (07-Jun-2021)

Molecular gas kinematics in the nuclear region of nearby Seyfert galaxies with ALMA

Authors: A. Bewketu Belete, P. Andreani, J. A. Fernández-Ontiveros, E. Hatziminaoglou, F. Combes, M. Sirressi, R. Slater, C. Ricci, K. Dasyra, C. Cicone, S. Aalto, L. Spinoglio, M. Imanishi, J. R. De Medeiros

Abstract: We present the analysis of the molecular gas in the nuclear regions of NGC 4968, NGC 4845, and MCG-06-30-15, with the help of ALMA observations of the CO(2-1) emission line. The aim is to determine the kinematics of the gas in the central (~ 1 kpc) region. We use the 3D-Based Analysis of Rotating Object via Line Observations ($^{3D}$BAROLO) and DiskFit softwares. Circular motions dominate the kine… ▽ More We present the analysis of the molecular gas in the nuclear regions of NGC 4968, NGC 4845, and MCG-06-30-15, with the help of ALMA observations of the CO(2-1) emission line. The aim is to determine the kinematics of the gas in the central (~ 1 kpc) region. We use the 3D-Based Analysis of Rotating Object via Line Observations ($^{3D}$BAROLO) and DiskFit softwares. Circular motions dominate the kinematics of the gas in the central discs, mainly in NGC 4845 and MCG-06-30-15, however there is a clear evidence of non-circular motions in the central ($\sim$ 1 kpc) region of NGC 4845 and NGC 4968. The strongest non-circular motion is detected in the inner disc of NGC 4968 with velocity $\sim 115\, \rm{km\,s^{-1}}$. The bisymmetric model is found to give the best-fit for NGC 4968 and NGC 4845. If the dynamics of NGC 4968 is modeled as a corotation pattern just outside of the bar, the bar pattern speed turns out to be at $Ω_b$ = $52\, \rm{km\,s^{-1}\,kpc^{-1}}$ the corotation is set at 3.5 kpc and the inner Lindblad resonance (ILR) ring at R = 300pc corresponding to the CO emission ring. The 1.2 mm ALMA continuum is peaked and compact in NGC 4968 and MCG-06-30-15, but their CO(2-1) has an extended distribution. Allowing the CO-to-H$_{2}$ conversion factor $α_{CO}$ between 0.8 and 3.2, typical of nearby galaxies of the same type, the molecular mass M(H$_{2}$) is estimated to be $\sim 3-12\times 10^{7} ~{\rm M_\odot}$ (NGC 4968), $\sim 9-36\times 10^{7}~ {\rm M_\odot}$ (NGC 4845), and $\sim 1-4\times 10^{7}~ {\rm M_\odot}$ (MCG-06-30-15). We conclude that the observed non-circular motions in the disc of NGC 4968 and likely that seen in NGC 4845 is due to the presence of the bar in the nuclear region. At the current spectral and spatial resolution and sensitivity we cannot claim any strong evidence in these sources of the long sought feedback/feeding effect due to the AGN presence. △ Less

Submitted 14 May, 2021; originally announced May 2021.

Comments: 14 pages, 11 figures, 3 tables, accepted for publication in the Astronomy & Astrophysics (A&A) journal

Journal ref: A&A 654, A24 (2021)

Calibration of mid- to far-infrared spectral lines in galaxies

Authors: Sabrina Mordini, Luigi Spinoglio, Juan Antonio Fernández-Ontiveros

Abstract: Mid- to far-infrared (IR) lines are suited to study dust obscured regions in galaxies, because IR spectroscopy allows us to explore the most hidden regions where heavily obscured star formation as well as accretion onto supermassive black-holes occur. This is mostly important at redshifts of 1<z<3, when most of the baryonic mass in galaxies has been assembled. We provide reliable calibrations of t… ▽ More Mid- to far-infrared (IR) lines are suited to study dust obscured regions in galaxies, because IR spectroscopy allows us to explore the most hidden regions where heavily obscured star formation as well as accretion onto supermassive black-holes occur. This is mostly important at redshifts of 1<z<3, when most of the baryonic mass in galaxies has been assembled. We provide reliable calibrations of the mid- to far-IR ionic fine structure lines, the brightest H2 pure rotational lines and the Polycyclic Aromatic Hydrocarbons (PAHs) features, that will be used to analyse current and future observations in the mm/submm range from the ground, as well as mid-IR spectroscopy from the upcoming James Webb Space Telescope. We use three samples of galaxies observed in the local Universe: star forming galaxies, AGN and low-metallicity dwarf galaxies. For each population we derive different calibrations of the observed line luminosities versus the total IR luminosities. We derive spectroscopic measurements of SFR and BHAR using mid- and far-IR fine structure lines, H2 pure rotational lines and PAH features. We derive robust star-formation tracers based on the [CII]158 $μ$m line; the sum of the [OI]63$μ$m and [OIII]88$μ$m lines; a combination of the neon and sulfur mid-IR lines; the bright PAH features at 6.2 and 11.3 $μ$m, and the H2 rotational lines at 9.7, 12.3 and 17 $μ$m. We propose the [CII]158$μ$m line, the combination of two neon lines and, for solar-like metallicity galaxies that may harbor an AGN, the PAH11.3$μ$m feature as the best SFR tracers. A reliable measure of the BHAR can be obtained using the [OIV]25.9 $μ$m and the [NeV]14.3 and 24.3 $μ$m lines. For the most commonly observed fine-structure lines in the far-IR we compare our calibration with the existing ALMA observations of high redshift galaxies finding overall a good agreement with local results. △ Less

Submitted 13 May, 2021; v1 submitted 10 May, 2021; originally announced May 2021.

Comments: 30 pages, 35 figures, 3 tables, Accepted for publication by A&A

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

X-ray binary accretion states in Active Galactic Nuclei? Sensing the accretion disc of supermassive black holes with mid-infrared nebular lines

Authors: Juan A. Fernández-Ontiveros, Teo Muñoz-Darias

Abstract: Accretion states, which are universally observed in stellar-mass black holes in X-ray binaries, are also anticipated in active galactic nuclei (AGN). This is the case at low luminosities, when the jet-corona coupling dominates the energy output in both populations. Previous attempts to extend this framework to a wider AGN population have been extremely challenging due to heavy hydrogen absorption… ▽ More Accretion states, which are universally observed in stellar-mass black holes in X-ray binaries, are also anticipated in active galactic nuclei (AGN). This is the case at low luminosities, when the jet-corona coupling dominates the energy output in both populations. Previous attempts to extend this framework to a wider AGN population have been extremely challenging due to heavy hydrogen absorption of the accretion disc continuum and starlight contamination from the host galaxies. We present the luminosity-excitation diagram (LED), based on the [OIV]$_{25.9μm}$ and [NeII]$_{12.8μm}$ mid-infrared nebular line fluxes. This tool enables to probe the accretion disc contribution to the ionising continuum. When applied to a sample of 167 nearby AGN, the LED recovers the characteristic q-shaped morphology outlined by individual X-ray binaries during a typical accretion episode, allowing us to tentatively identify the main accretion states. The soft state would include broad-line Seyferts and about half of the Seyfert 2 population, showing highly excited gas and radio-quiet cores consistent with disc-dominated nuclei. The hard state mostly includes low-luminosity AGN ($\leq 10^{-3}\, \rm{L_{Edd}}$) characterised by low-excitation radio-loud nuclei and a negligible disc contribution. The remaining half of Seyfert 2 nuclei and the bright LINERs show low excitation at high accretion luminosities, and could be identified with the bright-hard and intermediate states. Their hosts show ongoing star formation in the central kiloparsecs. We discuss the above scenario, its potential links with the galaxy evolution picture and the possible presence of accretion state transitions in AGN, as suggested by the growing population of changing-look quasars. △ Less

Submitted 19 April, 2021; originally announced April 2021.

Comments: Accepted for publication in MNRAS

Measuring Chemical Abundances with Infrared Nebular Lines: HII-Chi-mistry-IR

Authors: J. A. Fernández-Ontiveros, E. Pérez-Montero, J. M. Vílchez, R. Amorín, L. Spinoglio

Abstract: We provide a new method to derive heavy element abundances based on the unique suite of nebular lines in the mid- to far-infrared (IR) range. Using grids of photo-ionisation models that cover a wide range in O/H and N/O abundances, and ionisation parameter, our code HII-Chi-mistry-IR (HCm-IR) provides model-based abundances based on extinction free and temperature insensitive tracers, two signific… ▽ More We provide a new method to derive heavy element abundances based on the unique suite of nebular lines in the mid- to far-infrared (IR) range. Using grids of photo-ionisation models that cover a wide range in O/H and N/O abundances, and ionisation parameter, our code HII-Chi-mistry-IR (HCm-IR) provides model-based abundances based on extinction free and temperature insensitive tracers, two significant advantages over optical diagnostics. The code is probed using a sample of 56 galaxies observed with $Spitzer$ and $Herschel$ covering a wide range in metallicity, $7.2 \lesssim 12+\log(O/H) \lesssim 8.9$. The IR model-based metallicities obtained are robust within a scatter of 0.03 dex when the hydrogen recombination lines, which are typically faint transitions in the IR range, are not available. When compared to the optical abundances obtained with the direct method, model-based methods, and strong-line calibrations, HCm-IR estimates show a typical dispersion of ~0.2 dex, in line with previous studies comparing IR and optical abundances, a do not introduce a noticeable systematic above $12+\log(O/H) \gtrsim 7.6$. This accuracy can be achieved using the lines [SIV]$_{10.5 μm}$, [SIII]$_{18.7,33.5 μm}$, [NeIII]$_{15.6 μm}$ and [NeII]$_{12.8 μm}$. Additionally, HCm-IR provides an independent N/O measurement when the [OIII]$_{52,88 μm}$ and [NIII]$_{57 μm}$ transitions are measured, and therefore the derived abundances in this case do not rely on particular assumptions in the N/O ratio. Large uncertainties (~0.4 dex) may affect the abundance determinations of galaxies at sub- or over-solar metallicities when a solar-like N/O ratio is adopted. Finally, the code has been applied to 8 galaxies located at $1.8 < z < 7.5$ with ground-based detections of far-IR lines redshifted in the submm range, revealing solar-like N/O and O/H abundances in agreement with recent studies. △ Less

Submitted 16 March, 2021; originally announced March 2021.

Comments: Accepted for publication in A&A. The abstract has been abridged to fit Arxiv's character limit. A script to derive chemical abundances with HII-Chi-mistry-IR is available at https://www.iaa.csic.es/~epm/HII-CHI-mistry.html

Mid-IR cosmological spectrophotometric surveys from space: Measuring AGN and star formation at the Cosmic Noon with a SPICA-like mission

Authors: Luigi Spinoglio, Sabrina Mordini, Juan Antonio Fernandez-Ontiveros, Almudena Alonso-Herrero, Lee Armus, Laura Bisigello, Francesco Calura, Francisco J. Carrera, Asantha Cooray, Helmut Dannerbauer, Roberto Decarli, Eiichi Egami, David Elbaz, Alberto Franceschini, Eduardo Gonzalez Alfonso, Luca Graziani, Carlotta Gruppioni, Evanthia Hatziminaoglou, Hidehiro Kaneda, Kotaro Kohno, Alvaro Labiano, Georgios Magdis, Matthew A. Malkan, Hideo Matsuhara, Tohru Nagao , et al. (9 additional authors not shown)

Abstract: We use the SPace Infrared telescope for Cosmology and Astrophysics (SPICA) project as a template to demonstrate how deep spectrophotometric surveys covering large cosmological volumes over extended fields (1-15 square degrees) with a mid-IR imaging spectrometer (17-36 micron) in conjunction with deep 70 micron photometry with a far-IR camera, at wavelengths which are not affected by dust extinctio… ▽ More We use the SPace Infrared telescope for Cosmology and Astrophysics (SPICA) project as a template to demonstrate how deep spectrophotometric surveys covering large cosmological volumes over extended fields (1-15 square degrees) with a mid-IR imaging spectrometer (17-36 micron) in conjunction with deep 70 micron photometry with a far-IR camera, at wavelengths which are not affected by dust extinction can answer the most crucial questions in current galaxy evolution studies. A SPICA-like mission will be able for the first time to provide an unobscured three dimensional (3-D, i.e. x, y and redshift z) view of galaxy evolution back to an age of the Universe of less than ~2 Gyrs, in the mid-IR rest-frame. This survey strategy will produce a full census of the Star formation Rate (SFR) in the Universe, using Polycyclic Aromatic Hydrocarbons (PAH) bands and fine-structure ionic lines, reaching the characteristic knee of the galaxy luminosity function, where the bulk of the population is distributed, at any redshift up to z ~3.5. Deep follow-up pointed spectroscopic observations with grating spectrometers { onboard the satellite}, across the full IR spectral range (17-210 micron), would simultaneously measure Black Hole Accretion Rate (BHAR), from high-ionization fine-structure lines, and SFR, from PAH and low- to mid-ionization lines in thousands of galaxies from solar to low metallicities, down to the knee of their luminosity functions. The analysis of the resulting atlas of IR spectra will reveal the physical processes at play in evolving galaxies across cosmic time, especially its heavily dust-embedded phase during the activity peak at the cosmic noon (z ~1-3), through IR emission lines and features that are insensitive to the dust obscuration. △ Less

Submitted 5 March, 2021; originally announced March 2021.

Comments: Paper accepted for publication on PASA on 5th March 2021, as part of the SPICA Special Issue

Optical spectroscopy of 4U 1812-12: an ultra-compact X-ray binary seen through an H II region

Authors: M. Armas Padilla, T. Muñoz-Darias, F. Jiménez-Ibarra, J. A. Fernández-Ontiveros, J. Casares, M. A. P. Torres, J. García-Rojas, V. A. Cúneo, N. Degenaar

Abstract: The persistent, low-luminosity neutron star X-ray binary 4U 1812-12 is a potential member of the scarce family of ultra-compact systems. We performed deep photometric and spectroscopic optical observations with the 10.4 m Gran Telescopio Canarias in order to investigate the chemical composition of the accreted plasma, which is a proxy for the donor star class. We detect a faint optical counterpart… ▽ More The persistent, low-luminosity neutron star X-ray binary 4U 1812-12 is a potential member of the scarce family of ultra-compact systems. We performed deep photometric and spectroscopic optical observations with the 10.4 m Gran Telescopio Canarias in order to investigate the chemical composition of the accreted plasma, which is a proxy for the donor star class. We detect a faint optical counterpart (g~25, r~23) that is located in the background of the outskirts of the Sharpless 54 H II region, whose characteristic nebular lines superimpose on the X-ray binary spectrum. Once this is corrected for, the actual source spectrum lacks hydrogen spectral features. In particular, the Halpha emission line is not detected, with an upper limit (3 sigma) on the equivalent width of <1.3 A. Helium (He I) lines are neither observed, albeit our constraints are not restrictive enough to properly test the presence of this element. We also provide stringent upper limits on the presence of emission lines from other elements, such as C and O, which are typically found in ultra-compact systems with C-O white dwarfs donors. The absence of hydrogen features, the persistent nature of the source at low luminosity, as well as the low optical to X-ray flux ratio confirm 4U 1812-12 as a compelling ultra-compact X-ray binary candidate, for which we tentatively propose a He-rich donor based on the optical spectrum and the detection of short thermonuclear X-ray bursts. In this framework, we discuss the possible orbital period of the system according to disc instability and evolutionary models. △ Less

Submitted 15 September, 2020; originally announced September 2020.

Comments: Accepted for publication in A&A

Journal ref: A&A 644, A63 (2020)

$Herschel$/PACS OH Spectroscopy of Seyfert, LINER, and Starburst Galaxies

Authors: Jordan N. Runco, Matthew A. Malkan, Juan Antonio Fernández-Ontiveros, Luigi Spinoglio, Miguel Pereira-Santaella

Abstract: We investigated the 65$μ$m, 71$μ$m, 79$μ$m, 84$μ$m, 119$μ$m, and 163$μ$m OH doublets of 178 local (0 < $z$ < 0.35) galaxies. They were observed using the $Herschel$/PACS spectrometer, including Seyfert galaxies, LINERs, and star-forming galaxies. We observe these doublets exclusively in absorption (OH71), primarily in absorption (OH65, OH84), mostly in emission (OH79), only in emission (OH163) and… ▽ More We investigated the 65$μ$m, 71$μ$m, 79$μ$m, 84$μ$m, 119$μ$m, and 163$μ$m OH doublets of 178 local (0 < $z$ < 0.35) galaxies. They were observed using the $Herschel$/PACS spectrometer, including Seyfert galaxies, LINERs, and star-forming galaxies. We observe these doublets exclusively in absorption (OH71), primarily in absorption (OH65, OH84), mostly in emission (OH79), only in emission (OH163) and an approximately even mix of the both (OH119). In 19 galaxies we find P-Cygni or reverse P-Cygni line profiles in the OH doublets. We use several galaxy observables to probe spectral classification, brightness of a central AGN/starburst component, and radiation field strength. We find that OH79, OH119, and OH163 are more likely to display strong emission for bright, unobscured AGN. For less luminous, obscured AGN and non-active galaxies, we find populations of strong absorption (OH119), weaker emission (OH163), and a mix of weak emission and weak absorption (OH79). For OH65, OH71 and OH84, we do not find significant correlations with the observables listed above. We do find relationships between OH79 and OH119 with both the 9.7$μ$m silicate feature and Balmer decrement dust extinction tracers in which more dust leads to weaker emission / stronger absorption. The origin of emission for the observed OH doublets, whether from collisional excitation, or from radiative pumping by infrared photons, is discussed. △ Less

Submitted 15 September, 2020; originally announced September 2020.

Comments: 29 pages, 10 figures, 5 machine readable tables and 6 supplemental figures available upon request

Discovery of optical outflows and inflows in the black hole candidate GRS 1716-249

Authors: V. A. Cúneo, T. Muñoz-Darias, J. Sánchez-Sierras, F. Jiménez-Ibarra, M. Armas Padilla, D. A. H. Buckley, J. Casares, P. Charles, J. M. Corral-Santana, R. Fender, J. A. Fernández-Ontiveros, D. Mata Sánchez, G. Panizo-Espinar, G. Ponti, M. A. P. Torres

Abstract: We present optical spectroscopy obtained with the GTC, VLT and SALT telescopes during the decline of the 2016-2017 outburst of the black hole candidate GRS 1716-249 (Nova Oph 1993). Our 18-epoch data set spans 6 months and reveals that the observational properties of the main emission lines are very variable, even on time scales of a few hours. Several epochs are characterised by P-Cyg (as well as… ▽ More We present optical spectroscopy obtained with the GTC, VLT and SALT telescopes during the decline of the 2016-2017 outburst of the black hole candidate GRS 1716-249 (Nova Oph 1993). Our 18-epoch data set spans 6 months and reveals that the observational properties of the main emission lines are very variable, even on time scales of a few hours. Several epochs are characterised by P-Cyg (as well as flat-top and asymmetric) profiles in the H$α$, H$β$ and He II ($λ$4686) emission lines, implying the presence of an accretion disc wind, which is likely hot and dense. The wind's terminal velocity ($\sim$2000 km s$^{-1}$) is similar to that observed in other black hole X-ray transients. These lines also show transient and sharp red-shifted absorptions, taking the form of inverted P-Cyg profiles. We argue that these profiles can be explained by the presence of infalling material at $\sim$1300 km s$^{-1}$. We propose a failed wind scenario to explain this inflow and discuss other alternatives, such as obscuration produced by an accretion-related structure (e.g. the gas stream) in a high inclination system. △ Less

Submitted 28 September, 2020; v1 submitted 27 July, 2020; originally announced July 2020.

Comments: Published in MNRAS main journal

Journal ref: Volume 498, Issue 1, pp.25-32, 2020

Low optical polarisation at the core of the optically-thin jet of M87

Authors: A. Y. Fresco, J. A. Fernández-Ontiveros, M. A. Prieto, J. A. Acosta-Pulido, A. Merloni

Abstract: We study the optical linear and circular polarisation in the optically-thin regime of the core and jet of M87. Observations were acquired two days before the Event Horizon Telescope (EHT) campaign in early April 2017. A high degree ($\sim 20$ per cent) of linear polarisation (P$_{\rm lin}$) is detected in the bright jet knots resolved at $\sim 10\, \rm{arcsec}$ to $23\, \rm{arcsec}$ ($0.8$-… ▽ More We study the optical linear and circular polarisation in the optically-thin regime of the core and jet of M87. Observations were acquired two days before the Event Horizon Telescope (EHT) campaign in early April 2017. A high degree ($\sim 20$ per cent) of linear polarisation (P$_{\rm lin}$) is detected in the bright jet knots resolved at $\sim 10\, \rm{arcsec}$ to $23\, \rm{arcsec}$ ($0.8$-$1.8\, \rm{kpc}$) from the centre, whereas the nucleus and inner jet show P$_{\rm lin} \lesssim 5$ per cent. The position angle of the linear polarisation shifts by $\sim 90$ degrees from each knot to the adjacent ones, with the core angle perpendicular to the first knot. The nucleus was in a low level of activity (P$_{\rm lin} \sim 2$-$3$ per cent), and no emission was detected from HST-1. No circular polarisation was detected either in the nucleus or the jet above a $3σ$ level of P$_{\rm circ} \leq 1.5$ per cent, discarding the conversion of P$_{\rm lin}$ into P$_{\rm circ}$. A disordered magnetic field configuration or a mix of unresolved knots polarised along axes with different orientations could explain the low P$_{\rm lin}$. The latter implies a smaller size of the core knots, in line with current interferometric observations. Polarimetry with EHT can probe this scenario in the future. A steep increase of both P$_{\rm lin}$ and P$_{\rm circ}$ with increasing frequency is expected for the optically-thin domain, above the turnover point. This work describes the methodology to recover the four Stokes parameters using a $λ/4$ wave-plate polarimeter. △ Less

Submitted 16 June, 2020; v1 submitted 29 May, 2020; originally announced May 2020.

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

The physics of Galaxy Evolution with SPICA observations

Authors: Luigi Spinoglio, Juan A. Fernandez-Ontiveros, Sabrina Mordini

Abstract: The evolution of galaxies at Cosmic Noon (redshift 1<z<3) passed through a dust-obscured phase, during which most stars formed and black holes in galactic nuclei started to shine, which cannot be seen in the optical and UV, but it needs rest frame mid-to-far IR spectroscopy to be unveiled. At these frequencies, dust extinction is minimal and a variety of atomic and molecular transitions, tracing m… ▽ More The evolution of galaxies at Cosmic Noon (redshift 1<z<3) passed through a dust-obscured phase, during which most stars formed and black holes in galactic nuclei started to shine, which cannot be seen in the optical and UV, but it needs rest frame mid-to-far IR spectroscopy to be unveiled. At these frequencies, dust extinction is minimal and a variety of atomic and molecular transitions, tracing most astrophysical domains, occur. The future IR space telescope mission, SPICA, currently under evaluation for the 5th Medium Size ESA Cosmic Vision Mission, fully redesigned with its 2.5 m mirror cooled down to T < 8K will perform such observations. SPICA will provide for the first time a 3-dimensional spectroscopic view of the hidden side of star formation and black hole accretion in all environments, from voids to cluster cores over 90% of cosmic time. Here we outline what SPICA will do in galaxy evolution studies. △ Less

Submitted 30 March, 2020; originally announced March 2020.

Comments: Contributed talk at the IAU Symp 359 Galaxy Evolution and Feedback Across Different Environments, 2020 March 2-6, Bento Concalves, Brazil. arXiv admin note: text overlap with arXiv:1911.12181

Journal ref: Proc. IAU 15 (2020) 72-77

Unveiling the physical processes that regulate Galaxy Evolution with SPICA observations

Authors: Luigi Spinoglio, Juan A. Fernández-Ontiveros, Sabrina Mordini

Abstract: To study the dust obscured phase of the galaxy evolution during the peak of the Star Formation Rate (SFR) and the Black Hole Accretion Rate (BHAR) density functions ($z = 1 - 4$), rest frame mid-to-far infrared (IR) spectroscopy is needed. At these frequencies, dust extinction is at its minimum and a variety of atomic and molecular transitions, tracing most astrophysical domains, occur. The future… ▽ More To study the dust obscured phase of the galaxy evolution during the peak of the Star Formation Rate (SFR) and the Black Hole Accretion Rate (BHAR) density functions ($z = 1 - 4$), rest frame mid-to-far infrared (IR) spectroscopy is needed. At these frequencies, dust extinction is at its minimum and a variety of atomic and molecular transitions, tracing most astrophysical domains, occur. The future IR space telescope mission, SPICA, fully redesigned with its $2.5\, \rm{m}$ mirror cooled down to $T < 8\, \rm{K}$, will be able to perform such observations. With SPICA, we will: 1) obtain a direct spectroscopic measurement of the SFR and of the BHAR histories, 2) measure the evolution of metals and dust to establish the matter cycle in galaxies, 3) uncover the feedback and feeding mechanisms in large samples of distant galaxies, either AGN- or starburst-dominated, reaching lookback times of nearly 12 Gyr. SPICA large-area deep surveys will provide low-resolution, mid-IR spectra and continuum fluxes for unbiased samples of tens of thousands of galaxies, and even the potential to uncover the youngest, most luminous galaxies in the first few hundred million years. In this talk a brief review of the scientific preparatory work that has been done in extragalactic astronomy by the SPICA Collaboration will be given. △ Less

Submitted 27 November, 2019; originally announced November 2019.

Comments: Contributed talk at the IAU Symp 356 "Nuclear Activity in Galaxies Across Cosmic Time", ed. M. Povic et al

AGN types and unification model

Authors: Luigi Spinoglio, Juan Antonio Fernández-Ontiveros

Abstract: The motivation of the "unified model" is to explain the main properties of the large zoo of active galactic nuclei with a single physical object. The discovery of broad permitted lines in the polarized spectrum of type 2 Seyfert galaxies in the mid 80's led to the idea of an obscuring torus, whose orientation with respect to our line of sight was the reason of the different optical spectra. Howeve… ▽ More The motivation of the "unified model" is to explain the main properties of the large zoo of active galactic nuclei with a single physical object. The discovery of broad permitted lines in the polarized spectrum of type 2 Seyfert galaxies in the mid 80's led to the idea of an obscuring torus, whose orientation with respect to our line of sight was the reason of the different optical spectra. However, after many years of observations with different techniques, including IR and mm interferometry, the resulting properties of the observed dust structures differ from the torus model that would be needed to explain the type 1 vs type 2 dichotomy. Moreover, in the last years, multi-frequency monitoring of active galactic nuclei has shown an increasing number of transitions from one type to the other one, which cannot be explained in terms of the simple orientation of the dusty structure surrounding the active galactic nucleus (AGN). The interrelations between the AGN and the host galaxy, as also shown in the Magorrian relation, suggest that the evolution of the host galaxy may also have an important role in the observed manifestation of the nuclei. As an example, the observed delay between the maximum star formation activity and the onset of the AGN activity, and the higher occurrence of type 2 nuclei in star forming galaxies, have suggested the possible evolutionary path from, e.g., HII $\rightarrow$ AGN2 $\rightarrow$ AGN1. In the next years the models of unification need to also consider this observational framework and not only simple orientation effects. △ Less

Submitted 27 November, 2019; originally announced November 2019.

Comments: Invited talk at IAU Symp 356 "Nuclear Activity in Galaxies Across Cosmic Time", ed. M. Povic et al

Journal ref: Proc. IAU 15 (2019) 29-43

A CO molecular gas wind 340 pc away from the Seyfert 2 nucleus in ESO 420-G13 probes an elusive radio jet

Authors: J. A. Fernández-Ontiveros, K. M. Dasyra, E. Hatziminaoglou, M. A. Malkan, M. Pereira-Santaella, M. Papachristou, L. Spinoglio, F. Combes, S. Aalto, N. Nagar, M. Imanishi, P. Andreani, C. Ricci, R. Slater

Abstract: A prominent jet-driven outflow of CO(2-1) molecular gas is found along the kinematic minor axis of the Seyfert 2 galaxy ESO 420-G13, at a distance of $340-600\, \rm{pc}$ from the nucleus. The wind morphology resembles a characteristic funnel shape, formed by a highly collimated filamentary emission at the base, likely tracing the jet propagation through a tenuous medium, until a bifurcation point… ▽ More A prominent jet-driven outflow of CO(2-1) molecular gas is found along the kinematic minor axis of the Seyfert 2 galaxy ESO 420-G13, at a distance of $340-600\, \rm{pc}$ from the nucleus. The wind morphology resembles a characteristic funnel shape, formed by a highly collimated filamentary emission at the base, likely tracing the jet propagation through a tenuous medium, until a bifurcation point at $440\, \rm{pc}$ where the jet hits a dense molecular core and shatters, dispersing the molecular gas into several clumps and filaments within the expansion cone. We also trace the jet in ionised gas within the inner $\lesssim 340\, \rm{pc}$ using the [NeII]$_{\rm 12.8 μm}$ line emission, where the molecular gas follows a circular rotation pattern. The wind outflow carries a mass of $\sim 8 \times 10^6\, \rm{M_\odot}$ at an average wind projected speed of $\sim 160\, \rm{km\,s^{-1}}$, which implies a mass outflow rate of $\sim 14\, \rm{M_\odot\,yr^{-1}}$. Based on the structure of the outflow and the budget of energy and momentum, we discard radiation pressure from the active nucleus, star formation, and supernovae as possible launching mechanisms. ESO 420-G13 is the second case after NGC 1377 where the presence of a previously unknown jet is revealed due to its interaction with the interstellar medium, suggesting that unknown jets in feeble radio nuclei might be more common than expected. Two possible jet-cloud configurations are discussed to explain the presence of an outflow at such distance from the AGN. The outflowing gas will likely not escape, thus a delay in the star formation rather than quenching is expected from this interaction, while the feedback effect would be confined within the central few hundred parsecs of the galaxy. △ Less

Submitted 31 October, 2019; originally announced November 2019.

Comments: 12 pages, 8 figures. Accepted for publication in A&A

Journal ref: A&A 633, A127 (2020)