Showing 1–50 of 225 results for author: van der Tak, F

Single Aperture Large Telescope for Universe Studies (SALTUS): Science Overview

Authors: Gordon Chin, Carrie M. Anderson, Jennifer Bergner, Nicolas Biver, Gordon L. Bjoraker, Thibault Cavalie, Michael DiSanti, Jian-Rong Gao, Paul Hartogh, Leon K. Harding, Qing Hu, Daewook Kim, Craig Kulesa, Gert de Lange, David T. Leisawitz, Rebecca C. Levy, Arthur Lichtenberger, Daniel P. Marronh, Joan Najita, Trent Newswander, George H. Rieke, Dimitra Rigopoulou, Peter Roefsema, Nathan X. Roth, Kamber Schwarz , et al. (11 additional authors not shown)

Abstract: The SALTUS Probe mission will provide a powerful far-infrared (far-IR) pointed space observatory to explore our cosmic origins and the possibility of life elsewhere. The observatory employs an innovative deployable 14-m aperture, with a sunshield that will radiatively cool the off-axis primary to <45K. This cooled primary reflector works in tandem with cryogenic coherent and incoherent instruments… ▽ More The SALTUS Probe mission will provide a powerful far-infrared (far-IR) pointed space observatory to explore our cosmic origins and the possibility of life elsewhere. The observatory employs an innovative deployable 14-m aperture, with a sunshield that will radiatively cool the off-axis primary to <45K. This cooled primary reflector works in tandem with cryogenic coherent and incoherent instruments that span the 34 to 660 micron far-IR range at both high and moderate spectral resolutions. △ Less

Submitted 21 May, 2024; originally announced May 2024.

Comments: 49 pages, 10 figures, 3 tables, submitted to SPIE JATIS

Likelihood and appearance of life beyond the Earth: An astronomical perspective

Authors: Floris van der Tak

Abstract: As of 2023, over 5500 planets are known to orbit stars other than our Sun. We can measure their sizes and orbital periods, infer their masses and temperatures, and constrain their compositions. Based on these data, about 1% of extrasolar planets are potentially habitable for life as we know it, implying that of the billions of planets in our Galaxy, some may actually be inhabited, at least by micr… ▽ More As of 2023, over 5500 planets are known to orbit stars other than our Sun. We can measure their sizes and orbital periods, infer their masses and temperatures, and constrain their compositions. Based on these data, about 1% of extrasolar planets are potentially habitable for life as we know it, implying that of the billions of planets in our Galaxy, some may actually be inhabited, at least by microbes. However, recognizing signs of alien life forms is a major challenge for current technology, because of the wide range of conditions on extrasolar planets, and because of the wide range of forms that life may take. This chapter reviews observations of exoplanets and discusses astrobiological definitions of habitability and the likelihood of finding life beyond the Earth, both within and outside the Solar system. △ Less

Submitted 8 May, 2024; originally announced May 2024.

Comments: Book chapter, to appear in First Contact: Aliens and Humans in Contemporary Science Fiction, eds. M. van Dijk, F. Bosman, K. Glimmerveen (Springer). Aimed at audience without professional background in natural sciences

Correlation between accretion rate and free-free emission in protoplanetary disks -- A multi-wavelength analysis of central mm/cm emission in transition disks

Authors: Alessia A. Rota, Jurrian D. Meijerhof, Nienke van der Marel, Logan Francis, Floris S. van der Tak, Andrew D. Sellek

Abstract: The inner regions of protoplanetary disks are the locations where most of planets are thought to form and where processes that influence the global evolution of the disk, such as MHD-winds and photoevaporation, originate. Transition disks (TDs) with large inner dust cavities are the ideal targets to study the inner tens of au of disks, as the central emission can be fully disentangled from the out… ▽ More The inner regions of protoplanetary disks are the locations where most of planets are thought to form and where processes that influence the global evolution of the disk, such as MHD-winds and photoevaporation, originate. Transition disks (TDs) with large inner dust cavities are the ideal targets to study the inner tens of au of disks, as the central emission can be fully disentangled from the outer disk emission. We present a homogeneous multi-wavelength analysis of the continuum emission in a sample of 11 TDs. We investigate the nature of the central emission close to the star, distinguishing between thermal dust and free-free emission. Spatially resolved measurements of continuum emission from archival ALMA data are combined with literature cm-wave observations to study the spectral indices of the inner and outer disks separately. While the emission from the outer disks is consistent with thermal dust emission, 10/11 of the spectral indices estimated for the central emission close to the star suggest that this emission is free-free emission, likely associated with an ionized jet or a disk wind. No correlation between the free-free luminosity and the accretion luminosity or the X-ray luminosity is found, arguing against the photoevaporative wind origin. A sub-linear correlation between the ionized mass loss rate and the accretion rate onto the star is observed, suggesting an origin in an ionized jet. The relative lack of mm-dust grains in the majority of inner disks in transition disks suggests that either such dust grains have drifted quickly towards the central star, grain growth is less efficient in the inner disk, or grains grow rapidly to planetesimal sizes in the inner disk. The observed correlation between the ionized mass loss rate and the accretion rate suggests the outflow is strictly connected with the stellar accretion and that accretion in these disks is driven by a jet. △ Less

Submitted 31 January, 2024; v1 submitted 11 January, 2024; originally announced January 2024.

Comments: Accepted for publication in A&A

Tracing Evolution in Massive Protostellar Objects (TEMPO) -- I: Fragmentation and emission properties of massive star-forming clumps in a luminosity limited ALMA sample

Authors: A. Avison, G. A. Fuller, N. Asabre Frimpong, S. Etoka, M. Hoare, B. M. Jones, N. Peretto, A. Traficante, F. van der Tak, J. E. Pineda, M. Beltrán, F. Wyrowski, M. Thompson, S. Lumsden, Z. Nagy, T. Hill, S. Viti, F. Fontani, P. Schilke

Abstract: The role of massive ($\geq$ 8M$_{\odot}$) stars in defining the energy budget and chemical enrichment of the interstellar medium in their host galaxy is significant. In this first paper from the Tracing Evolution in Massive Protostellar Objects (TEMPO) project we introduce a colour-luminosity selected (L$_*$ $\sim$ 3$\times10^3$ to 1$\times10^5$ L$_{\odot}$) sample of 38 massive star forming regio… ▽ More The role of massive ($\geq$ 8M$_{\odot}$) stars in defining the energy budget and chemical enrichment of the interstellar medium in their host galaxy is significant. In this first paper from the Tracing Evolution in Massive Protostellar Objects (TEMPO) project we introduce a colour-luminosity selected (L$_*$ $\sim$ 3$\times10^3$ to 1$\times10^5$ L$_{\odot}$) sample of 38 massive star forming regions observed with ALMA at 1.3mm and explore the fragmentation, clustering and flux density properties of the sample. The TEMPO sample fields are each found to contain multiple fragments (between 2-15 per field). The flux density budget is split evenly (53%-47%) between fields where emission is dominated by a single high flux density fragment and those in which the combined flux density of fainter objects dominates. The fragmentation scales observed in most fields are not comparable with the thermal Jeans length, $λ_J$, being larger in the majority of cases, suggestive of some non-thermal mechanism. A tentative evolutionary trend is seen between luminosity of the clump and the `spectral line richness' of the TEMPO fields; with 6.7GHz maser associated fields found to be lower luminosity and more line rich. This work also describes a method of line-free continuum channel selection within ALMA data and a generalised approach used to distinguishing sources which are potentially star-forming from those which are not, utilising interferometric visibility properties. △ Less

Submitted 11 September, 2023; originally announced September 2023.

Comments: 42 pages, 17 figures, 4 Tables including 3 Appendices (with additional tables and figures). Accepted for publication in MNRAS

Dense gas and star formation in the Outer Milky Way

Authors: Jonathan Braine, Yan Sun, Yoshito Shimajiri, Floris F. S. van der Tak, Min Fang, Philippe André, Hao Chen, Yu Gao

Abstract: We present maps and spectra of the HCN(1-0) and HCO$^+$(1-0) lines in the extreme outer Galaxy, at galactocentric radii between 14 and 22 kpc, with the 13.7 meter Delingha telescope. The 9 molecular clouds were selected from a CO/$^{13}$CO survey of the outer quadrants. The goal is to better understand the structure of molecular clouds in these poorly studied subsolar metallicity regions and the r… ▽ More We present maps and spectra of the HCN(1-0) and HCO$^+$(1-0) lines in the extreme outer Galaxy, at galactocentric radii between 14 and 22 kpc, with the 13.7 meter Delingha telescope. The 9 molecular clouds were selected from a CO/$^{13}$CO survey of the outer quadrants. The goal is to better understand the structure of molecular clouds in these poorly studied subsolar metallicity regions and the relation with star formation. The lines are all narrow, less than 2km/s at half power, enabling detection of the HCN hyperfine structure in the stronger sources and allowing us to observationally test hyperfine collision rates. The hyperfine line ratios show that the HCN emission is optically thin with column densities estimated at N(HCN)~$3x10^{12}$\scm. The HCO$^+$ emission is approximately twice as strong as the HCN (taken as the sum of all components), in contrast with the inner Galaxy and nearby galaxies where they are similarly strong. For an abundance ratio $χ_{HCN}/χ_{HCO^+} = 3$, this requires a relatively low density solution for the dense gas, with n(H2) $\sim 10^3 - 10^4$\ccm. The $^{12}$CO/$^{13}$CO line ratios are similar to solar neighborhood values, roughly 7.5, despite the low $^{13}$CO abundance expected at such large radii. The HCO$^+$/CO and HCO$^+$/$^{13}$CO integrated intensity ratios are also standard at about 1/35 and 1/5 respectively. HCN is weak compared to the CO emission, with HCN/CO $\sim 1/70$ even after summing all hyperfine components. At the parsec scales observed here, the correlation between star formation, as traced by 24~$μ$m emission as is standard in extragalactic work, and dense gas via the HCN or HCO$^+$ emission, is poor, perhaps due to the lack of dynamic range. We find that the lowest dense gas fractions are in the sources at high galactic latitude (b>2, h>300pc above the plane), possibly due to lower pressure. △ Less

Submitted 8 June, 2023; originally announced June 2023.

Comments: 17 pages, 14 figures

Star & planet formation: Upcoming opportunities in the space-based infrared

Authors: Floris van der Tak

Abstract: While ALMA and JWST are revolutionizing our view of star and planet formation with their unprecedented sensitivity and resolution at submillimeter and near-IR wavelengths, many outstanding questions can only be answered with observations in the thermal (mid- and far-) infrared domain. Many of these questions require space-based observations, to achieve the necessary sensitivity and/or wavelength c… ▽ More While ALMA and JWST are revolutionizing our view of star and planet formation with their unprecedented sensitivity and resolution at submillimeter and near-IR wavelengths, many outstanding questions can only be answered with observations in the thermal (mid- and far-) infrared domain. Many of these questions require space-based observations, to achieve the necessary sensitivity and/or wavelength coverage. In particular, how do interstellar clouds develop filamentary structures and dense cores? What are the masses and luminosities of objects at the earliest stages of star formation? What are the gas masses of planet-forming disks, and how do these disks disperse during planet formation? How is refractory and volatile material distributed within the disks, and how does this evolve with time? This article reviews how upcoming and planned balloon-borne and space-based telescopes for the mid- and far-infrared will address these questions, and outlines which further missions will be needed beyond 2030, when the ELTs will be in full operation. △ Less

Submitted 14 December, 2022; v1 submitted 28 November, 2022; originally announced November 2022.

Comments: Contribution to the proceedings of the 7th Chile-Cologne-Bonn symposium: Physics and Chemistry of Star Formation; 5 pages. Version 2 includes useful comments on the original version; version 3 includes further comments

An accurate set of H$_3$O$^+ -$ H$_2$ collisional rate coefficients for non-LTE modelling of warm interstellar clouds

Authors: Sándor Demes, François Lique, Alexandre Faure, Floris F. S. van der Tak

Abstract: Hydronium (H$_3$O$^+$) was first detected in 1986 in interstellar molecular clouds. It was reported in many galactic diffuse and dense regions, as well as in extragalactic sources. H$_3$O$^+$ plays a major role both in interstellar oxygen and water chemistry. However, despite the large number of H$_3$O$^+$ observations, its collisional excitation was investigated only partially. In the present wor… ▽ More Hydronium (H$_3$O$^+$) was first detected in 1986 in interstellar molecular clouds. It was reported in many galactic diffuse and dense regions, as well as in extragalactic sources. H$_3$O$^+$ plays a major role both in interstellar oxygen and water chemistry. However, despite the large number of H$_3$O$^+$ observations, its collisional excitation was investigated only partially. In the present work we study the state-to-state rotational de-excitation of $ortho$- and $para$-H$_3$O$^+$ in collisions both with $ortho$- and $para$-H$_2$. The cross sections are calculated within the close-coupling formalism using a highly accurate potential energy surface developed for this system. The rate coefficients are computed up to $300$ K kinetic temperature. Transitions between the lowest 21 rotation-inversion states were studied for $para$-H$_3$O$^+$, and the lowest 11 states for $ortho$-H$_3$O$^+$, i.e. all levels with rotational energies below 430 K ($\sim 300$ cm$^{-1}$) are considered (up to $j\leq5$). In order to estimate the impact of the new rate coefficients on the astrophysical models for H$_3$O$^+$, radiative transfer calculations were also carried out. We have examined how the new collisional data affect the line intensities with respect to older data previously used for the interpretation of observations. By analysing all detected transitions we find that our new, accurate rate coefficients have a significant impact (typically within a factor of 2) on radiation temperatures, allowing more accurate estimation of column densities and relative abundances of hydronium, especially in warm molecular clouds, paving the path towards better interpretation of interstellar water and oxygen chemistry. △ Less

Submitted 4 November, 2022; originally announced November 2022.

Comments: 13 pages, 7 figures

MSC Class: 85-10; 85-08

Diagnosing the interstellar medium of galaxies with far-infrared emission lines II. [C II], [O I], [O III], [N II] and [N III] up to z=6

Authors: Andrés Felipe Ramos Padilla, Lingyu Wang, F. F. S. van der Tak, Scott Trager

Abstract: Gas cooling processes in the interstellar medium (ISM) are key to understanding how star-formation processes occur in galaxies. Far-infrared (FIR) fine-structure emission lines can be used as a tool to understand the gas conditions and trace the different phases of the ISM. We model the most important far-infrared (FIR) emission lines throughout cosmic time back to $z=6$ with cosmological hydrodyn… ▽ More Gas cooling processes in the interstellar medium (ISM) are key to understanding how star-formation processes occur in galaxies. Far-infrared (FIR) fine-structure emission lines can be used as a tool to understand the gas conditions and trace the different phases of the ISM. We model the most important far-infrared (FIR) emission lines throughout cosmic time back to $z=6$ with cosmological hydrodynamical simulations. We study how different physical parameters, such as the interstellar radiation field (ISRF) and metallicity, impact the ISM phases traced by FIR line luminosities and connect those with the star-formation rate (SFR). We implement a physically motivated multi-phase model of the ISM by post-processing EAGLE cosmological simulation with Cloudy look-up tables. In this model, we assume four phases of the ISM: dense molecular gas, neutral atomic gas, diffuse ionised gas (DIG) and HII regions. Our model shows good agreement with the observed luminosity-SFR relation up to $z=6$ in the FIR emission lines analysed and we also provide linear fits. Our predictions also agree with observations in terms of diagnostic diagrams involving various line ratios. We find that [C II] is the best SFR tracer of the FIR lines even though it traces multiple ISM phases, while [O III] and [N II] can be used to understand the DIG-HII balance in the ionised phase. In addition, line ratios like [C II]/[O III] and [N II/[O I] are useful to trace parameters such as ISRF, metallicity and specific star-formation rate. These results help to interpret observations of FIR line emission from the local Universe to high-$z$ galaxies. △ Less

Submitted 24 May, 2022; originally announced May 2022.

Comments: 22 pages, 18 figures, 3 tables (plus appendix). Submitted to A&A

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

Exogenous delivery of water to Mercury

Authors: Kateryna Frantseva, David Nesvorný, Michael Mueller, Floris F. S. van der Tak, Inge Loes ten Kate, Petr Pokorný

Abstract: Radar and spacecraft observations show the permanently shadowed regions around Mercury's North Pole to contain water ice and complex organic material. One possible source of this material are impacts by interplanetary dust particles (IDPs), asteroids, and comets. We have performed numerical simulations of the dynamical evolution of asteroids and comets over the few Myr and checked for their impa… ▽ More Radar and spacecraft observations show the permanently shadowed regions around Mercury's North Pole to contain water ice and complex organic material. One possible source of this material are impacts by interplanetary dust particles (IDPs), asteroids, and comets. We have performed numerical simulations of the dynamical evolution of asteroids and comets over the few Myr and checked for their impacts with Mercury. We use the N-body integrator RMVS/Swifter to propagate the Sun and the eight planets from their current positions. We add comets and asteroids to the simulations as massless test particles, based on their current orbital distributions. Asteroid impactors are assigned a probability of being water-rich (C-class) based on the measured distribution of taxonomic types. For comets, we assume a constant water fraction. For IDPs, we use a dynamical meteoroid model to compute the dust flux on Mercury. Relative to previous work on asteroid and comet impacts (Moses et al. 1999), we leverage 20 years of progress in minor body surveys. Immediate post-impact ejection of impactor material into outer space is taken into account as is the migration efficiency of water across Mercury's surface to the polar cold traps. We find that asteroids deliver $\sim 1 \times 10^{3}$ kg/yr of water to Mercury, comets deliver $\sim 1 \times 10^{3}$ kg/yr and IDPs deliver $\sim 16 \times 10^{3}$ kg/yr within a factor of several. Over a timescale of $\sim 1$ Gyr, this is enough to deliver the minimum amount of water required by the radar and MESSENGER observations. While other sources of water on Mercury are not ruled out by our analysis, we show that they are not required to explain the currently available observational lower limits. △ Less

Submitted 25 April, 2022; originally announced April 2022.

Comments: 38 pages, 10 figures, published in Icarus

A target list for searching for habitable exomoons

Authors: Vera Dobos, András Haris, Inga E. E. Kamp, Floris F. S. van der Tak

Abstract: We investigate the habitability of hypothethical moons orbiting known exoplanets. This study focuses on big, rocky exomoons that are capable of maintaining a significant atmosphere. To determine their habitability, we calculate the incident stellar radiation and the tidal heating flux arising in the moons as the two main contributors to the energy budget. We use the runaway greenhouse and the maxi… ▽ More We investigate the habitability of hypothethical moons orbiting known exoplanets. This study focuses on big, rocky exomoons that are capable of maintaining a significant atmosphere. To determine their habitability, we calculate the incident stellar radiation and the tidal heating flux arising in the moons as the two main contributors to the energy budget. We use the runaway greenhouse and the maximum greenhouse flux limits as a definition of habitability. For each exoplanet we run our calculations for plausible ranges of physical and orbital parameters for the moons and the planet using a Monte Carlo approach. We calculate the moon habitability probability for each planet which is the fraction of the investigated cases that lead to habitable conditions. Based on our results, we provide a target list for observations of known exoplanets of which the top 10 planets have more than 50\% chance for hosting habitable moons on stable orbits. Two especially promising candidates are Kepler-62~f and Kepler-16~b, both of them with known masses and radii. Our target list can help to detect the first habitable exomoon. △ Less

Submitted 25 April, 2022; originally announced April 2022.

Comments: accepted for publication in MNRAS, 10 pages

Dents in the Veil: Protostellar feedback in Orion

Authors: U. Kavak, J. Bally, J. R. Goicoechea, C. H. M. Pabst, F. F. S. van der Tak, A. G. G. M. Tielens

Abstract: Interest in stellar feedback has recently increased because new studies suggest that radiative and mechanical feedback from young massive stars regulate the physical and chemical composition of the interstellar medium (ISM) significantly. Recent SOFIA [CII] 158 micron observations of the Orion Veil revealed that the expanding bubble is powered by stellar winds and influenced by previously active m… ▽ More Interest in stellar feedback has recently increased because new studies suggest that radiative and mechanical feedback from young massive stars regulate the physical and chemical composition of the interstellar medium (ISM) significantly. Recent SOFIA [CII] 158 micron observations of the Orion Veil revealed that the expanding bubble is powered by stellar winds and influenced by previously active molecular outflows of ionizing massive stars. We aim to investigate the mechanical feedback on the whole Veil shell by searching for jets/outflows interacting with the Veil shell and determining the origin/driving mechanisms of these collisions. In the light of these findings, as well as the momenta of the dents and their dynamical timescales, we propose that the dents are created by the interaction of collimated jets/outflows from protostars with luminosities ranging from 10$^3$ to 10$^4$ $L_\odot$ indicating B-type stars in the Orion star-forming cloud with the surrounding Veil shell. However, it is challenging to pinpoint the driving stars as they may have moved from the original ejection points of the jets/outflows. We conclude that the dynamics of the expanding Veil shell is influenced not just by the O-type stars in the Trapezium cluster, but also by less massive stars, especially B-type, in the Orion Nebula. Mechanical feedback from protostars with a range of masses appears to play an important role in determining the morphology of [HII] regions and injecting turbulence into the medium. △ Less

Submitted 22 March, 2022; originally announced March 2022.

Comments: 18 pages, 18 figures, Accepted by A&A on March 22, 2022

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

Breaking Orion's Veil bubble with fossil outflows

Authors: U Kavak, J. R. Goicoechea, C. H. M. Pabst, J. Bally, F. F. S. van der Tak, A. G. G. M. Tielens

Abstract: The role of feedback in the self-regulation of star formation is a fundamental question in astrophysics. The Orion Nebula is the nearest site of ongoing and recent massive star formation. It is a unique laboratory for the study of stellar feedback. Recent SOFIA [CII] 158 $μ$m observations revealed an expanding bubble, the Veil shell, being powered by stellar winds and ionization feedback. We have… ▽ More The role of feedback in the self-regulation of star formation is a fundamental question in astrophysics. The Orion Nebula is the nearest site of ongoing and recent massive star formation. It is a unique laboratory for the study of stellar feedback. Recent SOFIA [CII] 158 $μ$m observations revealed an expanding bubble, the Veil shell, being powered by stellar winds and ionization feedback. We have identified a protrusion-like substructure in the Northwest portion of the Orion Veil Shell that may indicate additional feedback mechanisms that are highly directional. Our goal is to investigate the origin of the protrusion by quantifying its possible driving mechanisms. We use the [CII] 158 $μ$m map of the Orion Nebula obtained with the upGREAT instrument onboard SOFIA. The spectral and spatial resolution of the observations are 0.3 km/s and 16 arcsec, respectively. We consider three possible origins for this protrusion: Fossil outflow cavities created by jets/outflows during the protostellar accretion phase, pre-existing clumpiness in the OMC-1 core, and the stellar wind during the main sequence phase. Based on the energetics and the morphology, we conclude that the northwestern part of the pre-existing cloud was locally perturbed by outflows ejected from massive protostars in the Trapezium cluster. This suggests that the protrusion of the Veil is the result of mechanical rather than radiative feedback. Furthermore, we argue that the location of the protrusion is a suitable place to break the Orion Veil owing to the photo-ablation from the walls of the protrusion. We conclude that the outflows of massive protostars can influence the morphology of the future \hii\,region and even cause breakages in the ionization front. Specifically, the interaction of stellar winds of main-sequence stars with the molecular core pre-processed by the protostellar jet is important. △ Less

Submitted 9 February, 2022; originally announced February 2022.

Comments: Accepted by A&A on Feb 9, 2022 (20 pages, 19 figures)

Journal ref: A&A 660, A109 (2022)

The sharp ALMA view of infall and outflow in the massive protocluster G31.41+0.31

Authors: M. T. Beltrán, V. M. Rivilla, R. Cesaroni, D. Galli, L. Moscadelli, A. Ahmadi, H. Beuther, S. Etoka, C. Goddi, P. D. Klaassen, R. Kuiper, M. S. N. Kumar, A. Lorenzani, T. Peters, Á. Sánchez-Monge, P. Schilke, F. van der Tak, S. Vig

Abstract: Context. To better understand the formation of high-mass stars, it is fundamental to investigate how matter accretes onto young massive stars, how it is ejected, and how all this differs from the low-mass case. The massive protocluster G31.41+0.31 is the ideal target to study all these processes because observations at millimeter and centimeter wavelengths have resolved the emission of the Main co… ▽ More Context. To better understand the formation of high-mass stars, it is fundamental to investigate how matter accretes onto young massive stars, how it is ejected, and how all this differs from the low-mass case. The massive protocluster G31.41+0.31 is the ideal target to study all these processes because observations at millimeter and centimeter wavelengths have resolved the emission of the Main core into at least four massive dust continuum sources, named A, B, C, and D, within 1" or 0.018 pc, and have identified signatures of infall and several outflows associated with the core. Aims. We study the interplay between infall and outflow in G31.41+0.31 by investigating at a spatial resolution of a few 100 au their properties and their possible impact on the core. Methods. We carried out molecular line observations of typical high-density tracers, such as CH3CN or H2CO, and shock and outflow tracers, such as SiO, with ALMA at 1.4 mm that achieved an angular resolution of 0.09" (340 au). Results. The observations have revealed inverse P-Cygni profiles in CH3CN and H2CO toward the four sources embedded in the Main core, suggesting that all of them are undergoing collapse. The infall rates, estimated from the red-shifted absorption are on the order of 1E-2 Msun/yr. The individual infall rates imply that the accretion timescale of the Main core is an order of magnitude smaller than its rotation timescale. This confirms that rotating toroids such as the G31 Main core are non-equilibrium, transient collapsing structures that need to be constantly replenished with fresh material from a large-scale reservoir. For sources B, C, and D, the infall could be accelerating inside the sources, while for source A, the presence of a second emission component complicates the interpretation. The SiO observations have revealed the presence of at least six outflows in the G31.41+0.31 star-forming region, ... △ Less

Submitted 25 January, 2022; originally announced January 2022.

Comments: 17 pages, 9 figures, 2 tables. Accepted by A&A

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

Interrelations between Astrochemistry and Galactic Dynamics

Authors: E. Mendoza, N. Duronea, D. Ronsó, L. C. Corazza, F. van der Tak, S. Paron, L. -Å. Nyman

Abstract: This paper presents a review of ideas that interconnect Astrochemistry and Galactic Dynamics. Since these two areas are vast and not recent, each one has already been covered separately by several reviews. After a general historical introduction, and a needed quick review of processes like the stellar nucleosynthesis which gives the base to understand the interstellar formation of simple chemical… ▽ More This paper presents a review of ideas that interconnect Astrochemistry and Galactic Dynamics. Since these two areas are vast and not recent, each one has already been covered separately by several reviews. After a general historical introduction, and a needed quick review of processes like the stellar nucleosynthesis which gives the base to understand the interstellar formation of simple chemical compounds (H2, CO, NH3 and H2O), we focus on a number of topics which are at the crossing of the two areas, Dynamics and Astrochemistry. Astrochemistry is a flourishing field which intends to study the presence and formation of molecules as well as the influence of them into the structure, evolution and dynamics of astronomical objects. The progress in the knowledge on the existence of new complex molecules and of their process of formation originates from the observational, experimental and theoretical areas which compose the field. The interfacing areas include star formation, protoplanetary disks, the role of the spiral arms and the chemical abundance gradients in the galactic disk. It often happens that the physical conditions in some regions of the ISM are only revealed by means of molecular observations. To organise a classification of chemical evolution processes, we discuss about how astrochemistry can act in three different contexts: i. the chemistry of the early universe, including external galaxies, ii. star forming regions, and iii. AGB stars and circumstellar envelopes. We mention that our research is stimulated by plans for instruments and projects, such as the on-going LLAMA, which consists in the construction of a 12m sub-mm radio telescope in the Andes. Thus, modern and new facilities can play a key role in new discoveries not only in astrochemistry but also in radio astronomy and related areas. Furthermore, the research of the origin of life is also a stimulating perspective. △ Less

Submitted 21 May, 2021; originally announced May 2021.

Comments: 29 pages, 10 figures. Article type: Review. Manuscript accepted for publication in Frontiers in Astronomy and Space Sciences, section Fundamental Astronomy

The role of planetary interior in the long-term evolution of atmospheric CO2 on Earth-like exoplanets

Authors: M. Oosterloo, D. Höning, I. E. E. Kamp, F. F. S. van der Tak

Abstract: Context: The long-term carbonate-silicate cycle plays an important role in the evolution of Earth's climate and, therefore, may also be an important mechanism in the evolution of the climates of Earth-like exoplanets. Aims: We investigate the effects of radiogenic mantle heating, core size, and planetary mass on the evolution of the atmospheric partial $CO_2$ pressure, and the ability of a long-te… ▽ More Context: The long-term carbonate-silicate cycle plays an important role in the evolution of Earth's climate and, therefore, may also be an important mechanism in the evolution of the climates of Earth-like exoplanets. Aims: We investigate the effects of radiogenic mantle heating, core size, and planetary mass on the evolution of the atmospheric partial $CO_2$ pressure, and the ability of a long-term carbon cycle driven by plate tectonics to control the atmospheric $CO_2$ pressure. Methods: We developed a box-model which connects carbon cycling to parametrized mantle convection. The carbon cycle was coupled to the thermal evolution via the plate speed, which depends on the global Rayleigh number. Results: We find decreasing atmospheric $CO_2$ pressure with time, up to an order of magnitude over 10 Gyr. Higher abundances of radioactive isotopes result in higher $CO_2$ pressures. We find a decreasing Rayleigh number and plate speed toward planets with larger core mass fractions $f_c$, which leads to lower atmospheric $CO_2$ pressure. More massive planets may favor the development of more $CO_2$ rich atmospheres due to hotter interiors. Conclusions: The dependence of plate tectonics on mantle cooling has a significant effect on the long-term evolution of the atmospheric $CO_2$ pressure. Carbon cycling mediated by plate tectonics is efficient in regulating planetary climates for a wide range of mantle radioactive isotope abundances, planet masses and core sizes. More efficient carbon cycling on planets with a high mantle abundance of thorium or uranium highlights the importance of mapping the abundances of these elements in host stars of potentially habitable exoplanets. Inefficient carbon recycling on planets with a large core mass fraction ($f_c\gtrsim 0.8$) emphasizes the importance of precise mass-radius measurements of Earth-sized exoplanets. △ Less

Submitted 17 March, 2021; originally announced March 2021.

Comments: 20 pages, 16 figures, to appear in Astronomy & Astrophysics

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

Fragmentation in the massive G31.41+0.31 protocluster

Authors: M. T. Beltrán, V. M. Rivilla, R. Cesaroni, L. T. Maud, D. Galli, L. Moscadelli, A. Lorenzani, A. Ahmadi, H. Beuther, T. Csengeri, S. Etoka, C. Goddi, P. D. Klaassen, R. Kuiper, M. S. N. Kumar, T. Peters, Á. Sánchez-Monge, P. Schilke, F. van der Tak, S. Vig, H. Zinnecker

Abstract: Context. ALMA observations at 1.4 mm and 0.2'' (750au) angular resolution of the Main core in the high-mass star forming region G31.41+0.31 have revealed a puzzling scenario: on the one hand, the continuum emission looks very homogeneous and the core appears to undergo solid-body rotation, suggesting a monolithic core stabilized by the magnetic field; on the other hand, rotation and infall speed u… ▽ More Context. ALMA observations at 1.4 mm and 0.2'' (750au) angular resolution of the Main core in the high-mass star forming region G31.41+0.31 have revealed a puzzling scenario: on the one hand, the continuum emission looks very homogeneous and the core appears to undergo solid-body rotation, suggesting a monolithic core stabilized by the magnetic field; on the other hand, rotation and infall speed up toward the core center, where two massive embedded free-free continuum sources have been detected, pointing to an unstable core having undergone fragmentation. Aims. To establish whether the Main core is indeed monolithic or its homogeneous appearance is due to a combination of large dust opacity and low angular resolution, we carried out millimeter observations at higher angular resolution and different wavelengths. Methods. We carried out ALMA observations at 1.4 mm and 3.5 mm that achieved angular resolutions of 0.1''(375 au) and 0.075'' (280 au), respectively. VLA observations at 7 mm and 1.3 cm at even higher angular resolution, 0.05'' (190 au) and 0.07'' (260 au), respectively, were also carried out to better study the nature of the free-free continuum sources detected in the core. Results. The millimeter continuum emission of the Main core has been clearly resolved into at least four sources, A, B, C, and D, within 1'', indicating that the core is not monolithic. The deconvolved radii of the dust emission of the sources, estimated at 3.5 mm, are 400-500au, their masses range from 15 to 26 Msun, and their number densities are several 1E9 cm-3. Sources A and B, located closer to the center of the core and separated by 750 au, are clearly associated with two free-free continuum sources, likely thermal radio jets, and are the brightest in the core. The spectral energy distribution of these two sources and their masses and sizes are similar and suggest a common origin. △ Less

Submitted 8 March, 2021; originally announced March 2021.

Comments: 18 pages,9 figures, 5 tables. Accepted by A&A

Journal ref: A&A 648, A100 (2021)

A search for radio jets from massive young stellar objects. Association of radio jets with H2O and CH3OH masers

Authors: U. Kavak, A. Sanchez-Monge, A. Lopez-Sepulcre, R. Cesaroni, F. F. S. van der Tak, L. Moscadelli, M. T. Beltran, P. Schilke

Abstract: Recent theoretical and observational studies debate the similarities between the formation process of high-mass (>8 Msun) and low-mass stars. The formation of low-mass star formation is directly associated with the presence of disks and jets. According to this scenario, radio jets are expected to be common in high-mass star-forming regions. We aim to increase the number of known radio jets in high… ▽ More Recent theoretical and observational studies debate the similarities between the formation process of high-mass (>8 Msun) and low-mass stars. The formation of low-mass star formation is directly associated with the presence of disks and jets. According to this scenario, radio jets are expected to be common in high-mass star-forming regions. We aim to increase the number of known radio jets in high-mass star forming regions by searching for radio jet candidates at radio continuum wavelengths. We have used the Karl G. Jansky Very Large Array (VLA) to observe 18 high-mass star-forming regions in the C band (6 cm, ~1.0 arcsec resolution) and K band (1.3 cm, ~0.3 arcsec resolution). We have searched for radio jet candidates by studying the association of radio continuum sources with shock activity signposts. We have identified 7 as the most probable radio jets. The radio luminosity of the radio jet candidates is correlated with the bolometric luminosity and the outflow momentum rate. About 7-36% of the radio jet candidates are associated with non-thermal emission. The radio jet candidates associated with 6.7 GHz CH3OH maser emission are preferentially thermal winds and jets, while a considerable fraction of radio jet candidates associated with H2O masers show non-thermal emission, likely due to strong shocks. Our sample of 18 regions is divided in 8 less evolved, infrared-dark regions and 10 more evolved, infrared-bright regions. We have found that ~71% of the identified radio jet candidates are located in the more evolved regions. Similarly, 25% of the less evolved regions harbor one of the most probable radio jets, while up to 50% of the more evolved regions contain one of these radio jet candidates. This suggests that the detection of radio jets in high-mass star forming regions is larger in slightly more evolved regions. △ Less

Submitted 30 November, 2020; originally announced November 2020.

Comments: 40 pages, 165 figures, Accepted by A&A on November 27, 2020

Journal ref: A&A 645, A29 (2021)

Diagnosing the interstellar medium of galaxies with far-infrared emission lines I. The [C II] 158 microns line at z~0

Authors: Andrés Felipe Ramos Padilla, Lingyu Wang, Sylvia Ploeckinger, F. F. S. van der Tak, Scott Trager

Abstract: Atomic fine structure lines have been detected in the local Universe and at high redshifts over the past decades. The [C II] emission line at 158 $μ$m is an important observable as it provides constraints on the interstellar medium (ISM) cooling processes. We develop a physically motivated framework to simulate the production of far-infrared line emission from galaxies in a cosmological context. T… ▽ More Atomic fine structure lines have been detected in the local Universe and at high redshifts over the past decades. The [C II] emission line at 158 $μ$m is an important observable as it provides constraints on the interstellar medium (ISM) cooling processes. We develop a physically motivated framework to simulate the production of far-infrared line emission from galaxies in a cosmological context. This first paper sets out our methodology and describes its first application, simulating the [C II] 158 $μ$m line emission in the local Universe. We combine the output from EAGLE cosmological hydrodynamical simulations with a multi-phase model of the ISM. Gas particles are divided into three phases: dense molecular gas, neutral atomic gas and diffuse ionised gas (DIG). We estimate the [C II] line emission from the three phases using a set of Cloudy cooling tables. Our results agree with previous findings regarding the contribution of these three ISM phases to the [C II] emission. Our model shows good agreement with the observed ${\rm L_{[C II]}}$-star formation rate (SFR) relation in the local Universe within 0.4 dex scatter. The fractional contribution to the [C II] line from different ISM phases depends on the total SFR and metallicity. The neutral gas phase dominates the [C II] emission in galaxies with $\rm{SFR}\sim0.01$-$1\,\rm{M_{\odot}}\,\rm{yr^{-1}}$, but the ionised phase dominates at lower SFRs. Galaxies above solar metallicity exhibit lower ${\rm L_{[C II]}}$/SFR ratios for the neutral phase. In comparison, the ${\rm L_{[C II]}}$/SFR ratio in the DIG is stable when metallicity varies. We suggest that the reduced size of the neutral clouds, caused by increased SFRs, is the likely cause for the ${\rm L_{[C II]}}$ deficit at high infrared luminosities, although EAGLE simulations do not reach these luminosities at $z=0$. △ Less

Submitted 26 November, 2020; originally announced November 2020.

Comments: Accepted for publication in A&A; 16 pages, 11 figures, 4 tables (plus appendix)

Journal ref: A&A 645, A133 (2021)

Evidence for chromium hydride in the atmosphere of Hot Jupiter WASP-31b

Authors: Marrick Braam, Floris F. S. van der Tak, Katy L. Chubb, Michiel Min

Abstract: Context. The characterisation of exoplanet atmospheres has shown a wide diversity in compositions. Hot Jupiters have the appropriate temperatures to host metallic compounds, which should be detectable through transmission spectroscopy. Aims. We aim to detect exotic species in transmission spectra of hot Jupiters, specifically WASP-31b, by testing a variety of chemical species to explain the spectr… ▽ More Context. The characterisation of exoplanet atmospheres has shown a wide diversity in compositions. Hot Jupiters have the appropriate temperatures to host metallic compounds, which should be detectable through transmission spectroscopy. Aims. We aim to detect exotic species in transmission spectra of hot Jupiters, specifically WASP-31b, by testing a variety of chemical species to explain the spectrum. Methods. We conduct a re-analysis of publicly available transmission data of WASP-31b using the Bayesian retrieval framework TauREx II. We retrieve various combinations of the opacities of 25 species to determine the minimum set needed to fit the observed spectrum. Results. We report evidence for the spectroscopic signatures of chromium hydride (CrH), H2O and K in WASP-31b. Compared to a flat model without any signatures, a CrH-only model is preferred with a statistical significance of ~3.9 sigma. A model consisting of both CrH+H2O is found with ~2.6 sigma and ~3 sigma confidence over a CrH-only and H2O-only model respectively. Weak evidence for the addition of K is found at ~2.2 sigma over the H2O+CrH model, although the fidelity of the data point associated with this signature was questioned in earlier studies. Finally, the inclusion of collision-induced absorption and a Rayleigh scattering slope is found with ~3.5 sigma confidence over a flat model. This analysis presents the first evidence for signatures of CrH in a hot Jupiter atmosphere. At a retrieved temperature of ~1481 K, the atmosphere of WASP-31b is hot enough to host gaseous Cr-bearing species. The retrieved abundances agree well with predictions from thermal equilibrium chemistry and with the abundance found in an L-type brown dwarf atmosphere. However, additional retrievals using VLT FORS2 data lead to a non-detection of CrH. Future coverage with JWST has the potential to confirm the detection and discover other CrH features. △ Less

Submitted 22 December, 2020; v1 submitted 20 November, 2020; originally announced November 2020.

Comments: Accepted for publication in A&A (Nov 24, 2020). 12 pages, 6 figures, 9 tables

Journal ref: A&A 646, A17 (2021)

Re-analysis of the 267-GHz ALMA observations of Venus: No statistically significant detection of phosphine

Authors: I. A. G. Snellen, L. Guzman-Ramirez, M. R. Hogerheijde, A. P. S. Hygate, F. F. S. van der Tak

Abstract: Context: ALMA observations of Venus at 267 GHz have been presented in the literature that show the apparent presence of phosphine (PH3) in its atmosphere. Phosphine has currently no evident production routes on the planet's surface or in its atmosphere. Aims: The aim of this work is to assess the statistical reliability of the line detection by independent re-analysis of the ALMA data. Methods… ▽ More Context: ALMA observations of Venus at 267 GHz have been presented in the literature that show the apparent presence of phosphine (PH3) in its atmosphere. Phosphine has currently no evident production routes on the planet's surface or in its atmosphere. Aims: The aim of this work is to assess the statistical reliability of the line detection by independent re-analysis of the ALMA data. Methods: The ALMA data were reduced as in the published study, following the provided scripts. First the spectral analysis presented in the study was reproduced and assessed. Subsequently, the spectrum was statistically evaluated, including its dependence on selected ALMA baselines. Results: We find that the 12th-order polynomial fit to the spectral passband utilised in the published study leads to spurious results. Following their recipe, five other >10 sigma lines can be produced in absorption or emission within 60 km/s from the PH3 1-0 transition frequency by suppressing the surrounding noise. Our independent analysis shows a feature near the PH3 frequency at a ~2 sigma level, below the common threshold for statistical significance. Since the spectral data have a non-Gaussian distribution, we consider a feature at such level as statistically unreliable that cannot be linked to a false positive probability. Conclusions: We find that the published 267-GHz ALMA data provide no statistical evidence for phosphine in the atmosphere of Venus. △ Less

Submitted 16 November, 2020; v1 submitted 19 October, 2020; originally announced October 2020.

Comments: Accepted for publication in A&A after minor revisions

Mapping the H$_{2}$D$^{+}$ and N$_{2}$H$^{+}$ emission towards prestellar cores. Testing dynamical models of the collapse using gas tracers

Authors: E. Koumpia, L. Evans, J. Di Francesco, F. F. S. van der Tak, R. D. Oudmaijer

Abstract: The study of prestellar cores is critical as they set the initial conditions in star formation and determine the final mass of the stellar object. To date, several hypotheses are describing their gravitational collapse. We perform detailed line analysis and modelling of H$_{2}$D$^{+}$ 110 -111 and N$_{2}$H$^{+}$ 4-3 emission at 372 GHz, using 2'x2' maps (JCMT). Our goal is to test the most promine… ▽ More The study of prestellar cores is critical as they set the initial conditions in star formation and determine the final mass of the stellar object. To date, several hypotheses are describing their gravitational collapse. We perform detailed line analysis and modelling of H$_{2}$D$^{+}$ 110 -111 and N$_{2}$H$^{+}$ 4-3 emission at 372 GHz, using 2'x2' maps (JCMT). Our goal is to test the most prominent dynamical models by comparing the modelled gas kinematics and spatial distribution (H$_{2}$D$^{+}$ and N$_{2}$H$^{+}$) with observations towards four prestellar (L1544, L183, L694-2, L1517B) and one protostellar core (L1521f). We perform a detailed non-LTE radiative transfer modelling using RATRAN, where we compare the predicted spatial distribution and line profiles of H$_{2}$D$^{+}$ and N$_{2}$H$^{+}$ with observations towards all cores. To do so, we adopt the physical structure for each core predicted by three different dynamical models taken from literature: Quasi-Equilibrium Bonnor-Ebert Sphere (QE-BES), Singular Isothermal Sphere (SIS), and Larson-Penston (LP) flow. Our analysis provides an updated picture of the physical structure of prestellar cores. We find that the SIS model can be clearly excluded in explaining the gas emission towards the cores, but a larger sample is required to differentiate clearly between the LP flow, the QE-BES and the static models. All models of collapse underestimate the intensity of the gas emission by up to several factors towards the only protostellar core in our sample, indicating that different dynamics take place in different evolutionary core stages. If the LP model is confirmed towards a larger sample of prestellar cores, it would indicate that they may form by compression or accretion of gas from larger scales. If the QE-BES model is confirmed, it means that quasi hydrostatic cores can exist within turbulent ISM. △ Less

Submitted 18 September, 2020; originally announced September 2020.

Comments: 32 pages, 24 figures, 5 tables (Accepted for publication in A&A)

Journal ref: A&A 643, A61 (2020)

Enrichment of the HR 8799 planets by minor bodies and dust

Authors: K. Frantseva, M. Mueller, P. Pokorný, F. F. S. van der Tak, I. L. ten Kate

Abstract: In the Solar System, minor bodies and dust deliver various materials to planetary surfaces. Several exoplanetary systems are known to host inner and outer belts, analogues of the main asteroid belt and the Kuiper belt. We study the possibility that exominor bodies and exodust deliver volatiles and refractories to the exoplanets in the system HR8799 by performing N-body simulations. The model consi… ▽ More In the Solar System, minor bodies and dust deliver various materials to planetary surfaces. Several exoplanetary systems are known to host inner and outer belts, analogues of the main asteroid belt and the Kuiper belt. We study the possibility that exominor bodies and exodust deliver volatiles and refractories to the exoplanets in the system HR8799 by performing N-body simulations. The model consists of the host star, four giant planets (HR8799 e, d, c, and b), 650000 test particles representing the inner belt, and 1450000 test particles representing the outer belt. Moreover we modelled dust populations that originate from both belts. Within a million years, the two belts evolve towards the expected dynamical structure (also derived in other works), where mean-motion resonances with the planets carve the analogues of Kirkwood gaps. We find that, after this point, the planets suffer impacts by objects from the inner and outer belt at rates that are essentially constant with time, while dust populations do not contribute significantly to the delivery process. We convert the impact rates to volatile and refractory delivery rates using our best estimates of the total mass contained in the belts and their volatile and refractory content. Over their lifetime, the four giant planets receive between $10^{-4}$ and $10^{-3}M_\bigoplus$ of material from both belts. The total amount of delivered volatiles and refractories, ${5\times10^{-3}\textrm{M}_\bigoplus}$, is small compared to the total mass of the planets, $11\times10^{3}\textrm{M}_\bigoplus$. However, if the planets were formed to be volatile-rich, their exogenous enrichment in refractory material may well be significant and observable, for example with JWST-MIRI. If terrestrial planets exist within the snow line of the system, volatile delivery would be an important astrobiological mechanism and may be observable as atmospheric trace gases. △ Less

Submitted 27 May, 2020; originally announced May 2020.

Comments: 11 pages, 8 figures, accepted for publication in Astronomy&Astrophysics

Journal ref: A&A 638, A50 (2020)

Gas and dust cooling along the major axis of M33 (HerM33es) -- Herschel/PACS [CII] and [OI] observations

Authors: Carsten Kramer, Thomas Nikola, Sibylle Anderl, Frank Bertoldi, Mederic Boquien, Jonathan Braine, Christof Buchbender, Francoise Combes, Christian Henkel, Israel Hermelo, Frank Israel, Monica Relano, Markus Roellig, Karl Schuster, Fatemeh Tabatabaei, Floris van der Tak, Simon Verley, Paul van der Werf, Martina Wiedner, Emmanuel Xilouris

Abstract: M33 is a gas rich spiral galaxy of the Local Group. We investigate the relationship between the two major gas cooling lines and the total infrared (TIR) dust continuum. We mapped the emission of gas and dust in M33 using the far-infrared lines of [CII] and [OI](63um) and the TIR. The line maps were observed with Herschel/PACS. These maps have 50pc resolution and form a ~370pc wide stripe along its… ▽ More M33 is a gas rich spiral galaxy of the Local Group. We investigate the relationship between the two major gas cooling lines and the total infrared (TIR) dust continuum. We mapped the emission of gas and dust in M33 using the far-infrared lines of [CII] and [OI](63um) and the TIR. The line maps were observed with Herschel/PACS. These maps have 50pc resolution and form a ~370pc wide stripe along its major axis covering the sites of bright HII regions, but also more quiescent arm and inter-arm regions from the southern arm at 2kpc galacto-centric distance to the south out to 5.7kpc distance to the north. Full-galaxy maps of the continuum emission at 24um from Spitzer/MIPS, and at 70um, 100um, and 160um from PACS were combined to obtain a map of the TIR. TIR and [CII] intensities are correlated over more than two orders of magnitude. The range of TIR translates to a range of far ultraviolet (FUV) emission of G0,obs~2 to 200 in units of the average Galactic radiation field. The binned [CII]/TIR ratio drops with rising TIR, with large, but decreasing scatter. Fits of modified black bodies (MBBs) to the continuum emission were used to estimate dust mass surface densities and total gas column densities. A correction for possible foreground absorption by cold gas was applied to the [OI] data before comparing it with models of photon dominated regions (PDRs). Most of the ratios of [CII]/[OI] and ([CII]+[OI])/TIR are consistent with two model solutions. The median ratios are consistent with one solution at n~2x10^2 cm-3, G0~60, and and a second low-FUV solution at n~10^4 cm-3, G0~1.5. The bulk of the gas along the lines-of-sight is represented by a low-density, high-FUV phase with low beam filling factors ~1. A fraction of the gas may, however, be represented by the second solution. △ Less

Submitted 12 May, 2020; v1 submitted 7 May, 2020; originally announced May 2020.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 639, A61 (2020)

The Leiden Atomic and Molecular Database (LAMDA): Current Status, Recent Updates, and Future Plans

Authors: Floris van der Tak, François Lique, Alex Faure, John Black, Ewine van Dishoeck

Abstract: The Leiden Atomic and Molecular Database (LAMDA) collects spectroscopic information and collisional rate coefficients for molecules, atoms, and ions of astrophysical and astrochemical interest. We describe the developments of the database since its inception in 2005, and outline our plans for the near future. Such a database is constrained both by the nature of its uses and by the availability of… ▽ More The Leiden Atomic and Molecular Database (LAMDA) collects spectroscopic information and collisional rate coefficients for molecules, atoms, and ions of astrophysical and astrochemical interest. We describe the developments of the database since its inception in 2005, and outline our plans for the near future. Such a database is constrained both by the nature of its uses and by the availability of accurate data: we suggest ways to improve the synergies among users and suppliers of data. We summarize some recent developments in computation of collisional cross sections and rate coefficients. We consider atomic and molecular data that are needed to support astrophysics and astrochemistry with upcoming instruments that operate in the mid- and far-infrared parts of the spectrum. △ Less

Submitted 23 April, 2020; originally announced April 2020.

Comments: Accepted by Atoms for the special issue on Atomic and molecular databases; 26 pages, 6 tables

The subsurface habitability of small, icy exomoons

Authors: J. Tjoa, M. Mueller, F. F. S. van der Tak

Abstract: Assuming our Solar System as typical, exomoons may outnumber exoplanets. If their habitability fraction is similar, they would thus constitute the largest portion of habitable real estate in the Universe. Icy moons in our Solar System, such as Europa and Enceladus, have already been shown to possess liquid water, a prerequisite for life on Earth. We intend to investigate under what circumstances s… ▽ More Assuming our Solar System as typical, exomoons may outnumber exoplanets. If their habitability fraction is similar, they would thus constitute the largest portion of habitable real estate in the Universe. Icy moons in our Solar System, such as Europa and Enceladus, have already been shown to possess liquid water, a prerequisite for life on Earth. We intend to investigate under what circumstances small, icy moons may sustain subsurface oceans and thus be "subsurface habitable". We pay specific attention to tidal heating. We made use of a phenomenological approach to tidal heating. We computed the orbit averaged flux from both stellar and planetary (both thermal and reflected stellar) illumination. We then calculated subsurface temperatures depending on illumination and thermal conduction to the surface through the ice shell and an insulating layer of regolith. We adopted a conduction only model, ignoring volcanism and ice shell convection as an outlet for internal heat. In doing so, we determined at which depth, if any, ice melts and a subsurface ocean forms. We find an analytical expression between the moon's physical and orbital characteristics and the melting depth. Since this expression directly relates icy moon observables to the melting depth, it allows us to swiftly put an upper limit on the melting depth for any given moon. We reproduce the existence of Enceladus' subsurface ocean; we also find that the two largest moons of Uranus (Titania & Oberon) could well sustain them. Our model predicts that Rhea does not have liquid water. Habitable exomoon environments may be found across an exoplanetary system, largely irrespective of the distance to the host star. Small, icy subsurface habitable moons may exist anywhere beyond the snow line. This may, in future observations, expand the search area for extraterrestrial habitable environments beyond the circumstellar habitable zone. △ Less

Submitted 20 March, 2020; originally announced March 2020.

Comments: 18 pages, 12 figures. To be published in Astronomy & Astrophysics, section 10. Planets and planetary systems

Journal ref: A&A 636, A50 (2020)

Bringing high spatial resolution to the Far-infrared -- A giant leap for astrophysics

Authors: Hendrik Linz, Henrik Beuther, Maryvonne Gerin, Javier R. Goicoechea, Frank Helmich, Oliver Krause, Yao Liu, Sergio Molinari, Volker Ossenkopf-Okada, Jorge Pineda, Marc Sauvage, Eva Schinnerer, Floris van der Tak, Martina Wiedner

Abstract: The far-infrared (FIR) regime is one of the few wavelength ranges where no astronomical data with sub-arcsecond spatial resolution exist. Neither of the medium-term satellite projects like SPICA, Millimetron nor O.S.T. will resolve this malady. For many research areas, however, information at high spatial and spectral resolution in the FIR, taken from atomic fine-structure lines, from highly excit… ▽ More The far-infrared (FIR) regime is one of the few wavelength ranges where no astronomical data with sub-arcsecond spatial resolution exist. Neither of the medium-term satellite projects like SPICA, Millimetron nor O.S.T. will resolve this malady. For many research areas, however, information at high spatial and spectral resolution in the FIR, taken from atomic fine-structure lines, from highly excited carbon monoxide (CO), light hydrids, and especially from water lines would open the door for transformative science. A main theme will be to trace the role of water in proto-planetary disks, to observationally advance our understanding of the planet formation process and, intimately related to that, the pathways to habitable planets and the emergence of life. Furthermore, key observations will zoom into the physics and chemistry of the star-formation process in our own Galaxy, as well as in external galaxies. The FIR provides unique tools to investigate in particular the energetics of heating, cooling and shocks. The velocity-resolved data in these tracers will reveal the detailed dynamics engrained in these processes in a spatially resolved fashion, and will deliver the perfect synergy with ground-based molecular line data for the colder dense gas. △ Less

Submitted 16 February, 2020; originally announced February 2020.

Comments: White Paper, submitted to ESA in August 2019 in response to the "Voyage 2050" process; 24 pages, 12 figures, 1 table; three references updated from preprint to journal status, one reference added in comparison to the version submitted to ESA

Exploring the formation pathways of formamide near young O-type stars

Authors: Veronica Allen, Floris F. S. van der Tak, Ana López-Sepulcre, Álvaro Sánchez-Monge, Victor Rivilla, Ricardo Cesaroni, .

Abstract: Context. As a building block for amino acids, formamide (NH$_2$CHO) is an important molecule in astrobiology and astrochemistry, but its formation path in the interstellar medium is not understood well. Aims. We aim to find empirical evidence to support the chemical relationships of formamide to HNCO and H$_2$CO. Methods. We examine high angular resolution (~0.2") Atacama Large Millimeter/subm… ▽ More Context. As a building block for amino acids, formamide (NH$_2$CHO) is an important molecule in astrobiology and astrochemistry, but its formation path in the interstellar medium is not understood well. Aims. We aim to find empirical evidence to support the chemical relationships of formamide to HNCO and H$_2$CO. Methods. We examine high angular resolution (~0.2") Atacama Large Millimeter/submillimeter Array (ALMA) maps of six sources in three high-mass star-forming regions and compare the spatial extent, integrated emission peak position, and velocity structure of HNCO and H$_2$CO line emission with that of NH$_2$CHO by using moment maps. Through spectral modeling, we compare the abundances of these three species. Results. In these sources, the emission peak separation and velocity dispersion of formamide emission is most often similar to HNCO emission, while the velocity structure is generally just as similar to H$_2$CO and HNCO (within errors). From the spectral modeling, we see that the abundances between all three of our focus species are correlated, and the relationship between NH$_2$CHO and HNCO reproduces the previously demonstrated abundance relationship. Conclusions. In this first interferometric study, which compares two potential parent species to NH$_2$CHO, we find that all moment maps for HNCO are more similar to NH$_2$CHO than H$_2$CO in one of our six sources (G24 A1). For the other five sources, the relationship between NH$_2$CHO, HNCO, and H$_2$CO is unclear as the different moment maps for each source are not consistently more similar to one species as opposed to the other. △ Less

Submitted 27 December, 2019; originally announced December 2019.

Comments: Main text: 15 pages, 17 figures, 10 tables. Appendix: 21 pages, 28 figures

ALMA and ROSINA detections of phosphorus-bearing molecules: the interstellar thread between star-forming regions and comets

Authors: V. M. Rivilla, M. N. Drozdovskaya, K. Altwegg, P. Caselli, M. T. Beltrán, F. Fontani, F. F. S. van der Tak, R. Cesaroni, A. Vasyunin, M. Rubin, F. Lique, S. Marinakis, L. Testi, the ROSINA team

Abstract: To understand how Phosphorus-bearing molecules are formed in star-forming regions, we have analysed ALMA observations of PN and PO towards the massive star-forming region AFGL 5142, combined with a new analysis of the data of the comet 67P/Churyumov-Gerasimenko taken with the ROSINA instrument onboard Rosetta. The ALMA maps show that the emission of PN and PO arises from several spots associated w… ▽ More To understand how Phosphorus-bearing molecules are formed in star-forming regions, we have analysed ALMA observations of PN and PO towards the massive star-forming region AFGL 5142, combined with a new analysis of the data of the comet 67P/Churyumov-Gerasimenko taken with the ROSINA instrument onboard Rosetta. The ALMA maps show that the emission of PN and PO arises from several spots associated with low-velocity gas with narrow linewidths in the cavity walls of a bipolar outflow. PO is more abundant than PN in most of the spots, with the PO/PN ratio increasing as a function of the distance to the protostar. Our data favor a formation scenario in which shocks sputter phosphorus from the surface of dust grains, and gas-phase photochemistry induced by UV photons from the protostar allows efficient formation of the two species in the cavity walls. Our analysis of the ROSINA data has revealed that PO is the main carrier of P in the comet, with PO/PN>10. Since comets may have delivered a significant amount of prebiotic material to the early Earth, this finding suggests that PO could contribute significantly to the phosphorus reservoir during the dawn of our planet. There is evidence that PO was already in the cometary ices prior to the birth of the Sun, so the chemical budget of the comet might be inherited from the natal environment of the Solar System, which is thought to be a stellar cluster including also massive stars. △ Less

Submitted 26 November, 2019; originally announced November 2019.

Comments: Accepted in MNRAS (20 pages, 12 figures)

Origin of hydrogen fluoride emission in the Orion Bar: An excellent tracer for CO-dark H2 gas clouds

Authors: U. Kavak, F. F. S. van der Tak, A. G. G. M. Tielens, R. F. Shipman

Abstract: The hydrogen fluoride (HF) molecule is seen in absorption in the interstellar medium (ISM) along many lines of sight. Surprisingly, it is observed in emission toward the Orion Bar, which is an interface between the ionized region around the Orion Trapezium stars and the Orion molecular cloud. We aim to understand the origin of HF emission in the Orion Bar by comparing its spatial distribution with… ▽ More The hydrogen fluoride (HF) molecule is seen in absorption in the interstellar medium (ISM) along many lines of sight. Surprisingly, it is observed in emission toward the Orion Bar, which is an interface between the ionized region around the Orion Trapezium stars and the Orion molecular cloud. We aim to understand the origin of HF emission in the Orion Bar by comparing its spatial distribution with other tracers. We examine three mechanisms to explain the HF emission: thermal excitation, radiative dust pumping, and chemical pumping. We used a Herschel/HIFI strip map of the HF J = 1-0 line, covering 0.5' by 1.5' that is oriented perpendicular to the Orion Bar. We used the RADEX non-local thermodynamic equilibrium (non-LTE) code to construct the HF column density map. We use the Meudon PDR code to explain the morphology of HF. The bulk of the HF emission at 10 km s^-1 emerges from the CO-dark molecular gas that separates the ionization front from the molecular gas that is deeper in the Orion Bar. The excitation of HF is caused mainly by collisions with H2 at a density of 10^5 cm^-3 together with a small contribution of electrons in the interclump gas of the Orion Bar. Infrared pumping and chemical pumping are not important. We conclude that the HF J = 1-0 line traces CO-dark molecular gas. Similarly, bright photodissociation regions associated with massive star formation may be responsible for the HF emission observed toward active galactic nuclei. △ Less

Submitted 1 October, 2019; originally announced October 2019.

Comments: 12 pages, 12 figures, Accepted by A&A on September 23, 2019

Journal ref: A&A 631, A117 (2019)

Origin of the PN molecule in star-forming regions: the enlarged sample

Authors: F. Fontani, V. M. Rivilla, F. F. S. van der Tak, C. Mininni, M. T. Beltrán, P. Caselli

Abstract: Phosphorus nitride (PN) is the P-bearing species with the highest number of detections in star-forming regions. Multi-line studies of the molecule have shown that the excitation temperature of PN is usually lower than the gas kinetic temperature, suggesting that PN is likely in conditions of sub-thermal excitation. We present an analysis of PN which takes the possible sub-thermal excitation condit… ▽ More Phosphorus nitride (PN) is the P-bearing species with the highest number of detections in star-forming regions. Multi-line studies of the molecule have shown that the excitation temperature of PN is usually lower than the gas kinetic temperature, suggesting that PN is likely in conditions of sub-thermal excitation. We present an analysis of PN which takes the possible sub-thermal excitation conditions into account in a sample of 24 massive star-forming regions. We observed PN (2-1), (3-2), (4-3), and (6-5) with the IRAM-30m and APEX telescopes and detected PN lines in 15 of them. Together with 9 similar sources detected in PN in previous works, we have analysed the largest sample of star-forming regions to date, made of 33 sources with 24 detections in total (among which 13 are new detections). Hence, we have increased the number of star-forming regions detected in PN by more than a factor 2. Our analysis indicates that the PN lines are indeed sub-thermally excited, but well described by a single excitation temperature. We have compared line profiles and fractional abundances of PN and SiO, a typical shock tracer, and found that almost all objects detected in PN have high-velocity SiO wings. Moreover, the SiO and PN abundances with respect to H2 are correlated over several orders of magnitude, and uncorrelated with gas temperature. This clearly shows that the production of PN is strongly linked to the presence of shocked gas, and rules out alternative scenarios based on thermal evaporation from iced grain mantles. △ Less

Submitted 29 August, 2019; originally announced August 2019.

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

Effect of galaxy mergers on star formation rates

Authors: W. J. Pearson, L. Wang, M. Alpaslan, I. Baldry, M. Bilicki, M. J. I. Brown, M. W. Grootes, B. W. Holwerda, T. D. Kitching, S. Kruk, F. F. S. van der Tak

Abstract: Galaxy mergers and interactions are an integral part of our basic understanding of how galaxies grow and evolve over time. However, the effect that galaxy mergers have on star formation rates (SFR) is contested, with observations of galaxy mergers showing reduced, enhanced and highly enhanced star formation. We aim to determine the effect of galaxy mergers on the SFR of galaxies using statisticall… ▽ More Galaxy mergers and interactions are an integral part of our basic understanding of how galaxies grow and evolve over time. However, the effect that galaxy mergers have on star formation rates (SFR) is contested, with observations of galaxy mergers showing reduced, enhanced and highly enhanced star formation. We aim to determine the effect of galaxy mergers on the SFR of galaxies using statistically large samples of galaxies, totalling over 200\,000, over a large redshift range, 0.0 to 4.0. We train and use convolutional neural networks to create binary merger identifications (merger or non-merger) in the SDSS, KiDS and CANDELS imaging surveys. We then compare the galaxy main sequence subtracted SFR of the merging and non-merging galaxies to determine what effect, if any, a galaxy merger has on SFR. We find that the SFR of merging galaxies are not significantly different from the SFR of non-merging systems. The changes in the average SFR seen in the star forming population when a galaxy is merging are small, of the order of a factor of 1.2. However, the higher the SFR above the galaxy main sequence, the higher the fraction of galaxy mergers. Galaxy mergers have little effect on the SFR of the majority of merging galaxies compared to the non-merging galaxies. The typical change in SFR is less than 0.1~dex in either direction. Larger changes in SFR can be seen but are less common. The increase in merger fraction as the distance above the galaxy main sequence increases demonstrates that galaxy mergers can induce starbursts. △ Less

Submitted 27 August, 2019; originally announced August 2019.

Comments: 20 pages, 17 figures, 6 tables, 3 appendices, accepted for publication in A&A

Journal ref: A&A 631, A51 (2019)

Substructures in the Keplerian disc around the O-type (proto)star G17.64+0.16

Authors: L. T. Maud, R. Cesaroni, M. S. N. Kumar, V. M. Rivilla, A. Ginsburg, P. D. Klaassen, D. Harsono, A. Sanchez-Monge, A. Ahmadi, V. Allen, M. T. Beltran, H. Beuther, R. Galvan-Madrid, C. Goddi, M. G. Hoare, M. R. Hogerheijde, K. G. Johnston, R. Kuiper, L. Moscadelli, T. Peters, L. Testi, F. F. S. van der Tak, W. J. de Wit

Abstract: We present the highest angular resolution (20x15mas - 44x33au) Atacama Large Millimeter/sub-millimeter Array (ALMA) observations currently possible of the proto-O-star G17.64+0.16 in Band 6. The Cycle 5 observations with baselines out to 16km probes scales <50au and reveal the rotating disc around G17.64+0.16, a massive forming O-type star. The disc has a ring-like enhancement in the dust emission… ▽ More We present the highest angular resolution (20x15mas - 44x33au) Atacama Large Millimeter/sub-millimeter Array (ALMA) observations currently possible of the proto-O-star G17.64+0.16 in Band 6. The Cycle 5 observations with baselines out to 16km probes scales <50au and reveal the rotating disc around G17.64+0.16, a massive forming O-type star. The disc has a ring-like enhancement in the dust emission, especially visible as arc structures to the north and south. The Keplerian kinematics are most prominently seen in the vibrationally excited water line, H2O (Eu=3461.9K). The mass of the central source found by modelling the Keplerian rotation is consistent with 45+/-10Mo. The H30alpha (231.9GHz) radio-recombination line and the SiO (5-4) molecular line were detected at up to the 10 sigma$ level. The estimated disc mass is 0.6-2.6Mo under the optically thin assumption. Analysis of the Toomre Q parameter, in the optically thin regime, indicates that the disc stability is highly dependent on temperature. The disc currently appears stable for temperatures >150K, this does not preclude that the substructures formed earlier through disc fragmentation. △ Less

Submitted 15 June, 2019; originally announced June 2019.

Comments: 9 Total, inc references, appendix, long author list, 4 figures (inc appendix)

Journal ref: A&A 627, L6 (2019)

Probing the cold magnetized Universe with SPICA-POL (B-BOP)

Authors: Ph. André, A. Hughes, V. Guillet, F. Boulanger, A. Bracco, E. Ntormousi, D. Arzoumanian, A. J. Maury, J. -Ph. Bernard, S. Bontemps, I. Ristorcelli, J. M. Girart, F. Motte, K. Tassis, E. Pantin, T. Montmerle, D. Johnstone, S. Gabici, A. Efstathiou, Shantanu Basu, M. Béthermin, H. Beuther, J. Braine, J. Di Francesco, E. Falgarone , et al. (31 additional authors not shown)

Abstract: SPICA, the cryogenic infrared space telescope recently pre-selected for a `Phase A' concept study as one of the three remaining candidates for ESA's fifth medium class (M5) mission, is foreseen to include a far-infrared polarimetric imager (SPICA-POL, now called B-BOP), which would offer a unique opportunity to resolve major issues in our understanding of the nearby, cold magnetized Universe. This… ▽ More SPICA, the cryogenic infrared space telescope recently pre-selected for a `Phase A' concept study as one of the three remaining candidates for ESA's fifth medium class (M5) mission, is foreseen to include a far-infrared polarimetric imager (SPICA-POL, now called B-BOP), which would offer a unique opportunity to resolve major issues in our understanding of the nearby, cold magnetized Universe. This paper presents an overview of the main science drivers for B-BOP, including high dynamic range polarimetric imaging of the cold interstellar medium (ISM) in both our Milky Way and nearby galaxies. Thanks to a cooled telescope, B-BOP will deliver wide-field 100-350 micron images of linearly polarized dust emission in Stokes Q and U with a resolution, signal-to-noise ratio, and both intensity and spatial dynamic ranges comparable to those achieved by Herschel images of the cold ISM in total intensity (Stokes I). The B-BOP 200 micron images will also have a factor ~30 higher resolution than Planck polarization data. This will make B-BOP a unique tool for characterizing the statistical properties of the magnetized interstellar medium and probing the role of magnetic fields in the formation and evolution of the interstellar web of dusty molecular filaments giving birth to most stars in our Galaxy. B-BOP will also be a powerful instrument for studying the magnetism of nearby galaxies and testing galactic dynamo models, constraining the physics of dust grain alignment, informing the problem of the interaction of cosmic rays with molecular clouds, tracing magnetic fields in the inner layers of protoplanetary disks, and monitoring accretion bursts in embedded protostars. △ Less

Submitted 9 May, 2019; originally announced May 2019.

Comments: Accepted for publication in PASA; 37 pages, 16 figures

Journal ref: PASA, 36, e029 (2019)

Multi-line Herschel/HIFI observations of water reveal infall motions and chemical segregation around high-mass protostars

Authors: Floris van der Tak, Russ Shipman, Thierry Jacq, Fabrice Herpin, Jonathan Braine, Friedrich Wyrowski

Abstract: (Abridged) We use HIFI maps of the 987 GHz H2O 2(02)-1(11) emission to measure the sizes and shapes of 19 high-mass protostellar envelopes. To identify infall, we use HIFI spectra of the optically thin C18O 9-8 and H2O-18 1(11)-0(00) lines. The high-J C18O line traces the warm central material and redshifted H2O-18 1(11)-0(00) absorption indicates material falling onto the warm core. We probe smal… ▽ More (Abridged) We use HIFI maps of the 987 GHz H2O 2(02)-1(11) emission to measure the sizes and shapes of 19 high-mass protostellar envelopes. To identify infall, we use HIFI spectra of the optically thin C18O 9-8 and H2O-18 1(11)-0(00) lines. The high-J C18O line traces the warm central material and redshifted H2O-18 1(11)-0(00) absorption indicates material falling onto the warm core. We probe small-scale chemical differentiation by comparing H2O 752 and 987 GHz spectra with those of H2O-18. Our measured radii of the central part of the H2O 2(02)-1(11) emission are 30-40% larger than the predictions from spherical envelope models, and axis ratios are <2, which we consider good agreement. For 11 of the 19 sources, we find a significant redshift of the H2O-18 1(11)-0(00) line relative to C18O 9-8. The inferred infall velocities are 0.6-3.2 km/s, and estimated mass inflow rates range from 7e-5 to 2e-2 M0/yr, with the highest mass inflow rates occurring toward the sources with the highest masses, and possibly the youngest ages. The other sources show either expanding motions or H2O-18 lines in emission. The H2O-18 1(11)-0(00) line profiles are remarkably similar to the "differences" between the H2O 2(02)-1(11) and 2(11)-2(02) profiles, suggesting that the H2O-18 line and the H2O 2(02)-1(11) absorption originate just inside the radius where water evaporates from grains, typically 1000-5000 au from the center. In some sources, the H2O-18 line is detectable in the outflow, where no C18O emission is seen. Together, the H2O-18 absorption and C18O emission profiles show that the water abundance around high-mass protostars has at least three levels: low in the cool outer envelope, high within the 100 K radius, and very high in the outflowing gas. Thus, despite the small regions, the combination of lines presented here reveals systematic inflows and chemical information about the outflows. △ Less

Submitted 27 March, 2019; originally announced March 2019.

Comments: Accepted for publication in Astronomy & Astrophysics; 10 pages body + 10 pages appendix

Probing ISM Structure in Trumpler 14 & Carina I Using The Stratospheric Terahertz Observatory 2

Authors: Young Min Seo, Paul F. Goldsmith, Chris Walker, David J. Hollenbach, Mark G. Wolfire, Craig Kulesa, Volker Tolls, Pietro N. Bernasconi, Umit Kavak, Floris F. S. van der Tak, Russ Shipman, Jian Rong Gao, Alexander Tielens, Michael G. Burton, Harold Yorke, Erick Young, William L. Peters, Abram Young, Christopher Groppi, Kristina Davis, Jorge L. Pineda, William D. Langer, Jonathan H. Kawamura, Antony Stark, Gary Melnick , et al. (4 additional authors not shown)

Abstract: We present observations of the Trumpler 14/Carina I region carried out using the Stratospheric Terahertz Observatory 2 (STO2). The Trumpler 14/Carina I region is in the west part of the Carina Nebula Complex, which is one of the most extreme star-forming regions in the Milky Way. We observed Trumpler 14/Carina I in the 158 $μ$m transition of [C\,{\sc ii}] with a spatial resolution of 48$''$ and a… ▽ More We present observations of the Trumpler 14/Carina I region carried out using the Stratospheric Terahertz Observatory 2 (STO2). The Trumpler 14/Carina I region is in the west part of the Carina Nebula Complex, which is one of the most extreme star-forming regions in the Milky Way. We observed Trumpler 14/Carina I in the 158 $μ$m transition of [C\,{\sc ii}] with a spatial resolution of 48$''$ and a velocity resolution of 0.17 km s$^{-1}$. The observations cover a 0.25$^\circ$ by 0.28$^\circ$ area with central position {\it l} = 297.34$^\circ$, {\it b} = -0.60$^\circ$. The kinematics show that bright [C\,{\sc ii}] structures are spatially and spectrally correlated with the surfaces of CO clouds, tracing the photodissociation region and ionization front of each molecular cloud. Along 7 lines of sight that traverse Tr 14 into the dark ridge to the southwest, we find that the [C\,{\sc ii}] luminosity from the HII region is 3.7 times that from the PDR. In same los we find in the PDRs an average ratio of 1:4.1:5.6 for the mass in atomic gas:dark-CO gas: molecular gas traced by CO. Comparing multiple gas tracers including HI 21cm, [C\,{\sc ii}], CO, and radio recombination lines, we find that the HII regions of the Carina Nebula Complex are well-described as HII regions with one-side freely expanding towards us, consistent with the champagne model of ionized gas evolution. The dispersal of the GMC in this region is dominated by EUV photoevaporation; the dispersal timescale is 20-30 Myr. △ Less

Submitted 22 May, 2019; v1 submitted 22 March, 2019; originally announced March 2019.

Comments: ApJ accepted

Identifying Galaxy Mergers in Observations and Simulations with Deep Learning

Authors: W. J. Pearson, L. Wang, J. W. Trayford, C. E. Petrillo, F. F. S. van der Tak

Abstract: Mergers are an important aspect of galaxy formation and evolution. We aim to test whether deep learning techniques can be used to reproduce visual classification of observations, physical classification of simulations and highlight any differences between these two classifications. With one of the main difficulties of merger studies being the lack of a truth sample, we can use our method to test b… ▽ More Mergers are an important aspect of galaxy formation and evolution. We aim to test whether deep learning techniques can be used to reproduce visual classification of observations, physical classification of simulations and highlight any differences between these two classifications. With one of the main difficulties of merger studies being the lack of a truth sample, we can use our method to test biases in visually identified merger catalogues. A convolutional neural network architecture was developed and trained in two ways: one with observations from SDSS and one with simulated galaxies from EAGLE, processed to mimic the SDSS observations. The SDSS images were also classified by the simulation trained network and the EAGLE images classified by the observation trained network. The observationally trained network achieves an accuracy of 91.5% while the simulation trained network achieves 65.2% on the visually classified SDSS and physically classified EAGLE images respectively. Classifying the SDSS images with the simulation trained network was less successful, only achieving an accuracy of 64.6%, while classifying the EAGLE images with the observation network was very poor, achieving an accuracy of only 53.0% with preferential assignment to the non-merger classification. This suggests that most of the simulated mergers do not have conspicuous merger features and visually identified merger catalogues from observations are incomplete and biased towards certain merger types. The networks trained and tested with the same data perform the best, with observations performing better than simulations, a result of the observational sample being biased towards conspicuous mergers. Classifying SDSS observations with the simulation trained network has proven to work, providing tantalizing prospects for using simulation trained networks for galaxy identification in large surveys. △ Less

Submitted 10 May, 2019; v1 submitted 27 February, 2019; originally announced February 2019.

Comments: Submitted to A&A, revised after first referee report. 20 pages, 22 figures, 14 tables, 1 appendix

Journal ref: A&A 626, A49 (2019)

Deep Learning for Galaxy Mergers in the Galaxy Main Sequence

Authors: William J. Pearson, Lingyu Wang, James Trayford, Carlo E. Petrillo, Floris F. S. van der Tak

Abstract: Starburst galaxies are often found to be the result of galaxy mergers. As a result, galaxy mergers are often believed to lie above the galaxy main sequence: the tight correlation between stellar mass and star formation rate. Here, we aim to test this claim. Deep learning techniques are applied to images from the Sloan Digital Sky Survey to provide visual-like classifications for over 340 000 objec… ▽ More Starburst galaxies are often found to be the result of galaxy mergers. As a result, galaxy mergers are often believed to lie above the galaxy main sequence: the tight correlation between stellar mass and star formation rate. Here, we aim to test this claim. Deep learning techniques are applied to images from the Sloan Digital Sky Survey to provide visual-like classifications for over 340 000 objects between redshifts of 0.005 and 0.1. The aim of this classification is to split the galaxy population into merger and non-merger systems and we are currently achieving an accuracy of 91.5%. Stellar masses and star formation rates are also estimated using panchromatic data for the entire galaxy population. With these preliminary data, the mergers are placed onto the full galaxy main sequence, where we find that merging systems lie across the entire star formation rate - stellar mass plane. △ Less

Submitted 22 January, 2019; originally announced January 2019.

Comments: 4 pages, 1 figure. For Proceedings IAU Symposium No. 341

Journal ref: Proc. IAU 15 (2019) 104-108

Chasing discs around O-type (proto)stars - ALMA evidence for an SiO disc and disc wind from G17.64+0.16

Authors: L. T. Maud, R. Cesaroni, M. S. N. Kumar, F. F. S. van der Tak, V. Allen, M. G. Hoare, P. D. Klaassen, D. Harsono, M. R. Hogerheijde, Á. Sánchez-Monge, P. Schilke, A. Ahmadi, M. T. Beltrán, H. Beuther, T. Csengeri, S. Etoka, G. Fuller, R. Galván-Madrid, C. Goddi, Th. Henning, K. G. Johnston, R. Kuiper, S. Lumsden, L. Moscadelli, J. C. Mottram , et al. (6 additional authors not shown)

Abstract: We present high angular resolution 0.2 arcsec continuum and molecular emission line Atacama Large Millimeter/sub-millimeter Array (ALMA) observations of G17.64+0.16 in Band 6 (220GHz) taken as part of a campaign in search of circumstellar discs around (proto)-O-stars. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peak… ▽ More We present high angular resolution 0.2 arcsec continuum and molecular emission line Atacama Large Millimeter/sub-millimeter Array (ALMA) observations of G17.64+0.16 in Band 6 (220GHz) taken as part of a campaign in search of circumstellar discs around (proto)-O-stars. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. We detect SiO (5-4) emission that is marginally resolved and elongated in a direction perpendicular to the large-scale outflow seen in the 13CO (2-1) line using the main ALMA array in conjunction with the Atacama Compact Array (ACA). Morphologically, the SiO appears to represent a disc-like structure. Using parametric models we show that the position-velocity profile of the SiO is consistent with the Keplerian rotation of a disc around an object between 10-30Mo in mass, only if there is also radial expansion from a separate structure. The radial motion component can be interpreted as a disc wind from the disc surface. Models with a central stellar object mass between 20 and 30Mo are the most consistent with the stellar luminosity (100000 Lo) and indicative of an O-type star. The H30a millimetre recombination line (231.9GHz) is also detected, but spatially unresolved, and is indicative of a very compact, hot, ionised region co-spatial with the dust continuum core. Accounting for all observables, we suggest that G17.64 is consistent with a O-type young stellar object in the final stages of protostellar assembly, driving a wind, but that has not yet developed into a compact HII region. The existance and detection of the disc in G17.64 is likely related to its isolated and possibly more evolved nature, traits which may underpin discs in similar sources. △ Less

Submitted 9 October, 2018; originally announced October 2018.

Comments: 23 pages, 13 fig (inc 3 Appendix figures)

Journal ref: A&A 620, A31 (2018)

Probing the High-Redshift Universe with SPICA: Toward the Epoch of Reionization and Beyond

Authors: E. Egami, S. Gallerani, R. Schneider, A. Pallottini, L. Vallini, E. Sobacchi, A. Ferrara, S. Bianchi, M. Bocchio, S. Marassi, L. Armus, L. Spinoglio, A. W. Blain, M. Bradford, D. L. Clements, H. Dannerbauer, J. A. Fernández-Ontiveros, E. González-Alfonso, M. J. Griffin, C. Gruppioni, H. Kaneda, K. Kohno, S. C. Madden, H. Matsuhara, P. Najarro , et al. (14 additional authors not shown)

Abstract: With the recent discovery of a dozen dusty star-forming galaxies and around 30 quasars at z>5 that are hyper-luminous in the infrared ($μ$$L_{\rm IR}>10^{13}$ L$_{\odot}$, where $μ$ is a lensing magnification factor), the possibility has opened up for SPICA, the proposed ESA M5 mid-/far-infrared mission, to extend its spectroscopic studies toward the epoch of reionization and beyond. In this paper… ▽ More With the recent discovery of a dozen dusty star-forming galaxies and around 30 quasars at z>5 that are hyper-luminous in the infrared ($μ$$L_{\rm IR}>10^{13}$ L$_{\odot}$, where $μ$ is a lensing magnification factor), the possibility has opened up for SPICA, the proposed ESA M5 mid-/far-infrared mission, to extend its spectroscopic studies toward the epoch of reionization and beyond. In this paper, we examine the feasibility and scientific potential of such observations with SPICA's far-infrared spectrometer SAFARI, which will probe a spectral range (35-230 $μ$m) that will be unexplored by ALMA and JWST. Our simulations show that SAFARI is capable of delivering good-quality spectra for hyper-luminous infrared galaxies (HyLIRGs) at z=5-10, allowing us to sample spectral features in the rest-frame mid-infrared and to investigate a host of key scientific issues, such as the relative importance of star formation versus AGN, the hardness of the radiation field, the level of chemical enrichment, and the properties of the molecular gas. From a broader perspective, SAFARI offers the potential to open up a new frontier in the study of the early Universe, providing access to uniquely powerful spectral features for probing first-generation objects, such as the key cooling lines of low-metallicity or metal-free forming galaxies (fine-structure and H2 lines) and emission features of solid compounds freshly synthesized by Population III supernovae. Ultimately, SAFARI's ability to explore the high-redshift Universe will be determined by the availability of sufficiently bright targets (whether intrinsically luminous or gravitationally lensed). With its launch expected around 2030, SPICA is ideally positioned to take full advantage of upcoming wide-field surveys such as LSST, SKA, Euclid, and WFIRST, which are likely to provide extraordinary targets for SAFARI. △ Less

Submitted 25 September, 2018; v1 submitted 20 September, 2018; originally announced September 2018.

Comments: 24 pages, 8 figures, accepted for publication in PASA (SPICA special issue)

Highly Ordered and Pinched Magnetic Fields in the Class 0 Proto-Binary System L1448 IRS 2

Authors: Woojin Kwon, Ian Stephens, John Tobin, Leslie Looney, Zhi-Yun Li, Floris van der Tak, Richard Crutcher

Abstract: We have carried out polarimetric observations with the Atacama Large Millimeter/submillimeter Array (ALMA) toward the Class 0 protostellar system L1448 IRS 2, which is a proto-binary embedded within a flattened, rotating structure, and for which a hint of a central disk has been suggested, but whose magnetic fields are aligned with the bipolar outflow on the cloud core scale. Our high sensitivity… ▽ More We have carried out polarimetric observations with the Atacama Large Millimeter/submillimeter Array (ALMA) toward the Class 0 protostellar system L1448 IRS 2, which is a proto-binary embedded within a flattened, rotating structure, and for which a hint of a central disk has been suggested, but whose magnetic fields are aligned with the bipolar outflow on the cloud core scale. Our high sensitivity and high resolution ($\sim 100$ au) observations show a clear hourglass magnetic field morphology centered on the protostellar system, but the central pattern is consistent with a toroidal field indicative of a circumstellar disk, although other interpretations are also possible, including field lines dragged by an equatorial accretion flow into a configuration parallel to the midplane. If a relatively large disk does exist, it would suggest that the magnetic braking catastrophe is averted in this system, not through a large misalignment between the magnetic and rotation axes, but rather through some other mechanisms, such as non-ideal magneto-hydrodynamic effects and/or turbulence. We have also found a relationship of decreasing polarization fractions with intensities and the various slopes of this relationship can be understood as multiple polarization mechanisms and/or depolarization from a changing field morphology. In addition, we found a prominent clumpy depolarization strip crossing the center perpendicular to the bipolar outflow. Moreover, a rough estimate of the magnetic field strength indicates that the field is strong enough to hinder formation of a rotationally supported disk, which is inconsistent with the feature of a central toroidal field. △ Less

Submitted 28 May, 2019; v1 submitted 18 May, 2018; originally announced May 2018.

Comments: 14 pages, accepted for publication in ApJ

Main sequence of star forming galaxies beyond the Herschel confusion limit

Authors: W. J. Pearson, L. Wang, P. D. Hurley, K. Małek, V. Buat, D. Burgarella, D. Farrah, S. J. Oliver, D. J. B. Smith, F. F. S. van der Tak

Abstract: Deep far-infrared (FIR) cosmological surveys are known to be affected by source confusion, causing issues when examining the main sequence (MS) of star forming galaxies. This has typically been partially tackled by the use of stacking. However, stacking only provides the average properties of the objects in the stack. This work aims to trace the MS over $0.2\leq z<6.0$ using the latest de-blended… ▽ More Deep far-infrared (FIR) cosmological surveys are known to be affected by source confusion, causing issues when examining the main sequence (MS) of star forming galaxies. This has typically been partially tackled by the use of stacking. However, stacking only provides the average properties of the objects in the stack. This work aims to trace the MS over $0.2\leq z<6.0$ using the latest de-blended Herschel photometry, which reaches $\approx10$ times deeper than the 5$σ$ confusion limit in SPIRE. This provides more reliable star formation rates (SFRs), especially for the fainter galaxies, and hence a more reliable MS. We built a pipeline that uses the spectral energy distribution (SED) modelling and fitting tool CIGALE to generate flux density priors in the Herschel SPIRE bands. These priors were then fed into the de-blending tool XID+ to extract flux densities from the SPIRE maps. Multi-wavelength data were combined with the extracted SPIRE flux densities to constrain SEDs and provide stellar mass (M$_{\star}$) and SFRs. These M$_{\star}$ and SFRs were then used to populate the SFR-M$_{\star}$ plane over $0.2\leq z<6.0$. No significant evidence of a high-mass turn-over was found; the best fit is thus a simple two-parameter power law of the form log(SFR)$=α[$log(M$_{\star})-10.5]+β$. The normalisation of the power law increases with redshift, rapidly at $z\lesssim1.8$, from $0.58\pm0.09$ at $z\approx0.37$ to $1.31\pm0.08$ at $z\approx1.8$. The slope is also found to increase with redshift, perhaps with an excess around $1.8\leq z<2.9$. The increasing slope indicates that galaxies become more self-similar as redshift increases, implying that the specific SFR of high-mass galaxies increases over $z=0.2$ to $z=6.0$, becoming closer to that of low-mass galaxies. The excess in the slope at $1.8\leq z<2.9$, if present, coincides with the peak of the cosmic star formation history. △ Less

Submitted 14 May, 2018; v1 submitted 10 April, 2018; originally announced April 2018.

Comments: 15 pages, 14 figures, 3 tables, 3 appendices, accepted for publication in A&A

Journal ref: A&A 615, A146 (2018)

Complex cyanides as chemical clocks in hot cores

Authors: V. Allen, F. F. S. van der Tak, C. Walsh

Abstract: In the high-mass star-forming region G35.20-0.74N, small scale (about 800 AU) chemical segregation has been observed in which complex organic molecules containing the CN group are located in a small location. We aim to determine the physical origin of the large abundance difference (about 4 orders of magnitude) in complex cyanides within G35.20-0.74 B, and we explore variations in age, gas and dus… ▽ More In the high-mass star-forming region G35.20-0.74N, small scale (about 800 AU) chemical segregation has been observed in which complex organic molecules containing the CN group are located in a small location. We aim to determine the physical origin of the large abundance difference (about 4 orders of magnitude) in complex cyanides within G35.20-0.74 B, and we explore variations in age, gas and dust temperature, and gas density. We performed gas-grain astrochemical modeling experiments with exponentially increasing (coupled) gas and dust temperature rising from 10 to 500 K at constant H$_2$ densities of 10$^7$, 10$^8$, and 10$^9$ cm$^{-3}$. We tested the effect of varying the initial ice composition, cosmic-ray ionization rate, warm-up time (over 50, 200, and 1000 kyr), and initial (10, 15, and 25 K) and final temperatures (300 and 500 K). Varying the initial ice compositions within the observed and expected ranges does not noticeably affect the modeled abundances indicating that the chemical make-up of hot cores is determined in the warm-up stage. Complex cyanides vinyl and ethyl cyanide (CH$_2$CHCN and C$_2$H$_5$CN, respectively) cannot be produced in abundances (versus H$_2$) greater than 5x10$^{-10}$ for CH$_2$CHCN and 2x10$^{-10}$ for C$_2$H$_5$CN with a fast warm-up time (52 kyr), while the lower limit for the observed abundance of C$_2$H$_5$CN toward source B3 is 3.4x10$^{-10}$. Complex cyanide abundances are reduced at higher initial temperatures and increased at higher cosmic-ray ionization rates. Reproducing the observed abundances toward G35.20-0.74 Core B3 requires a fast warm-up at a high cosmic-ray ionization rate (1x10$^{-16}$ s$^{-1}$) at a high gas density (>10$^9$ cm$^{-3}$). G35.20-0.74 source B3 only needs to be about 2000 years older than B1/B2 for the observed chemical difference to be present. (This abstract has been shortened) △ Less

Submitted 9 April, 2018; v1 submitted 22 March, 2018; originally announced April 2018.

Comments: 14 pages (plus 11 pages appendices), 12 figures main text

SPICA - a large cryogenic infrared space telescope Unveiling the obscured Universe

Authors: P. R. Roelfsema, H. Shibai, L. Armus, D. Arrazola, M. Audard, M. D. Audley, C. M. Bradford, I. Charles, P. Dieleman, Y. Doi, L. Duband, M. Eggens, J. Evers, I. Funaki, J. R. Gao, M. Giard, A. di Giorgio L. M. González Fernández, M. Griffin, F. P. Helmich, R. Hijmering, R. Huisman, D. Ishihara, N. Isobe, B. Jackson, H. Jacobs , et al. (44 additional authors not shown)

Abstract: Measurements in the infrared wavelength domain allow us to assess directly the physical state and energy balance of cool matter in space, thus enabling the detailed study of the various processes that govern the formation and early evolution of stars and planetary systems in galaxies over cosmic time. Previous infrared missions, from IRAS to Herschel, have revealed a great deal about the obscured… ▽ More Measurements in the infrared wavelength domain allow us to assess directly the physical state and energy balance of cool matter in space, thus enabling the detailed study of the various processes that govern the formation and early evolution of stars and planetary systems in galaxies over cosmic time. Previous infrared missions, from IRAS to Herschel, have revealed a great deal about the obscured Universe, but sensitivity has been limited because up to now it has not been possible to fly a telescope that is both large and cold. SPICA is a mission concept aimed at taking the next step in mid- and far-infrared observational capability by combining a large and cold telescope with instruments employing state-of-the-art ultra-sensitive detectors. The mission concept foresees a 2.5-meter diameter telescope cooled to below 8 K. With cooling provided by mechanical coolers instead of depending on a limited cryogen supply, the mission lifetime can extend significantly beyond the required three years. SPICA offers instrumentation with spectral resolving powers ranging from R ~50 through 11000 in the 17-230 $μ$m domain as well as R~28.000 spectroscopy between 12 and 18 $μ$m. Additionally SPICA will provide efficient 30-37 $μ$m broad band mapping, and polarimetric imaging in the 100-350 $μ$m range. SPICA will provide unprecedented spectroscopic sensitivity of ~5 x $10^{-20}$ W/m$^2$ (5$σ$/1hr) - at least two orders of magnitude improvement over what has been attained to date. With this exceptional leap in performance, new domains in infrared astronomy will become accessible, allowing us, for example, to unravel definitively galaxy evolution and metal production over cosmic time, to study dust formation and evolution from very early epochs onwards, and to trace the formation history of planetary systems. △ Less

Submitted 28 March, 2018; originally announced March 2018.

Comments: 34 pages, 22 figures, paper accepted for publication in PASA on 2nd February 2018, part of the PASA SPICA Special Issue

Accelerating infall and rotational spin-up in the hot molecular core G31.41+0.31

Authors: M. T. Beltrán, R. Cesaroni, V. M. Rivilla, Á. Sánchez-Monge, L. Moscadelli, A. Ahmadi, V. Allen, H. Beuther, S. Etoka, D. Galli, R. Galván-Madrid, C. Goddi, K. G. Johnston, A. Kölligan, R. Kuiper, M. S. N. Kumar, L. T. Maud, J. C. Mottram, T. Peters, P. Schilke, L. Testi, F. van der Tak, C. M. Walmsley

Abstract: As part of our effort to search for circumstellar disks around high-mass stellar objects, we observed the well-known core G31.41+0.31 with ALMA at 1.4 mm with an angular resolution of~0.22" (~1700 au). The dust continuum emission has been resolved into two cores namely Main and NE. The Main core, which has the stronger emission and is the more chemically rich, has a diameter of ~5300 au, and is as… ▽ More As part of our effort to search for circumstellar disks around high-mass stellar objects, we observed the well-known core G31.41+0.31 with ALMA at 1.4 mm with an angular resolution of~0.22" (~1700 au). The dust continuum emission has been resolved into two cores namely Main and NE. The Main core, which has the stronger emission and is the more chemically rich, has a diameter of ~5300 au, and is associated with two free-free continuum sources. The Main core looks featureless and homogeneous in dust continuum emission and does not present any hint of fragmentation. Each transition of CH3CN and CH3OCHO, both ground and vibrationally excited, as well as those of CH3CN isotopologues, shows a clear velocity gradient along the NE-SW direction, with velocity linearly increasing with distance from the center, consistent with solid-body rotation. However, when comparing the velocity field of transitions with different upper level energies, the rotation velocity increases with increasing energy of the transition, which suggests that the rotation speeds up towards the center. Spectral lines towards the dust continuum peak show an inverse P-Cygni profile that supports the existence of infall in the core. The infall velocity increases with the energy of the transition suggesting that the infall is accelerating towards the center of the core, consistent with gravitational collapse. Despite the monolithic appearance of the Main core, the presence of red-shifted absorption, the existence of two embedded free-free sources at the center, and the rotational spin-up are consistent with an unstable core undergoing fragmentation with infall and differential rotation due to conservation of angular momentum. Therefore, the most likely explanation for the monolithic morphology is that the large opacity of the dust emission prevents the detection of any inhomogeneity in the core. △ Less

Submitted 21 June, 2018; v1 submitted 14 March, 2018; originally announced March 2018.

Comments: 28 pages, 16 figures, and 7 tables. Accepted for publication in A&A

Delivery of organics to Mars through asteroid and comet impacts

Authors: Kateryna Frantseva, Michael Mueller, Inge Loes ten Kate, Floris F. S. van der Tak, Sarah Greenstreet

Abstract: Given rapid photodissociation and photodegradation, the recently discovered organics in the Martian subsurface and atmosphere were probably delivered in geologically recent times. Possible parent bodies are C-type asteroids, comets, and interplanetary dust particles (IDPs). The dust infall rate was estimated, using different methods, to be between $0.71$ and $2.96 \times 10^6$ kg/yr (Nesvorny et a… ▽ More Given rapid photodissociation and photodegradation, the recently discovered organics in the Martian subsurface and atmosphere were probably delivered in geologically recent times. Possible parent bodies are C-type asteroids, comets, and interplanetary dust particles (IDPs). The dust infall rate was estimated, using different methods, to be between $0.71$ and $2.96 \times 10^6$ kg/yr (Nesvorny et al., 2011, Borin et al., 2017, Crismani et al., 2017); assuming a carbon content of 10% (Flynn, 1996), this implies an IDP carbon flux of $0.07 - 0.3 \times 10^6$ kg/yr. We calculate for the first time the carbon flux from impacts of asteroids and comets. To this end, we perform dynamical simulations of impact rates on Mars. We use the N-body integrator RMVS/Swifter to propagate the Sun and the eight planets from their current positions. We separately add comets and asteroids to the simulations as massless test particles, based on their current orbital elements, yielding Mars impact rates of $4.34\times10^{-3}$ comets/Myr and 3.3 asteroids/Myr. We estimate the global carbon flux on Mars from cometary impacts to be $\sim 0.013 \times 10^{6}$~kg/yr within an order of magnitude, while asteroids deliver $\sim 0.05 \times 10^6$~kg/yr. These values correspond to $\sim 4-19 \%$ and $\sim 17-71 \%$, respectively, of the IDP-borne carbon flux estimated by Nesvorny et al. 2011, Borin et al. 2017 and Crismani et al. 2017. Unlike the spatially homogeneous IDP infall, impact ejecta are distributed locally, concentrated around the impact site. We find organics from asteroids and comets to dominate over IDP-borne organics at distances up to 150~km from the crater center. Our results may be important for the interpretation of in situ detections of organics on Mars. △ Less

Submitted 8 March, 2018; originally announced March 2018.

Comments: 30 pages, 5 figures, accepted for publication in Icarus

The evolution of young HII regions

Authors: P. D. Klaassen, K. G. Johnston, J. S. Urquhart, J. C. Mottram, T. Peters, R. Kuiper, H. Beuther, F. F. S. van der Tak, C. Goddi

Abstract: High-mass stars form in much richer environments than those associated with isolated low-mass stars, and once they reach a certain mass, produce ionised (HII) regions. The formation of these pockets of ionised gas are unique to the formation of high-mass stars (M $>8$ M$_\odot$), and present an excellent opportunity to study the final stages of accretion, which could include accretion through the… ▽ More High-mass stars form in much richer environments than those associated with isolated low-mass stars, and once they reach a certain mass, produce ionised (HII) regions. The formation of these pockets of ionised gas are unique to the formation of high-mass stars (M $>8$ M$_\odot$), and present an excellent opportunity to study the final stages of accretion, which could include accretion through the HII region itself. This study of the dynamics of the gas on both sides of these ionisation boundaries in very young HII regions aims to quantify the relationship between the HII regions and their immediate environments.We present high-resolution ($\sim$ 0.5$"$) ALMA observations of nine HII regions selected from the Red MSX Source (RMS) survey with compact radio emission and bolometric luminosities greater than 10$^4$ L$_\odot$. We focus on the initial presentation of the data, including initial results from the radio recombination line H29$α$, some complementary molecules, and the 256 GHz continuum emission. Of the six (out of nine) regions with H29$α$ detections, two appear to have cometary morphologies with velocity gradients across them, and two appear more spherical with velocity gradients suggestive of infalling ionised gas. The remaining two were either observed at low resolution or had signals that were too weak to draw robust conclusions. We also present a description of the interactions between the ionised and molecular gas (as traced by CS (J=5-4)), often (but not always) finding theHII region had cleared its immediate vicinity of molecules. Of our sample of nine, the observations of the two clusters expected to have the youngest HII regions (from previous radio observations) are suggestive of having infalling motions in the H29$α$ emission, which could be indicative of late stage accretion onto the stars despite the presence of an HII region. △ Less

Submitted 13 December, 2017; originally announced December 2017.

Comments: Accepted for publication in A&A. 25 pages, 15 figures, 9 tables, with Table A.2 made available electronically via the CDS

Journal ref: A&A 611, A99 (2018)

Probing the baryon cycle of galaxies with SPICA mid- and far-infrared observations

Authors: F. F. S. van der Tak, S. C. Madden, P. Roelfsema, L. Armus, M. Baes, J. Bernard-Salas, A. Bolatto, S. Bontemps, C. Bot, C. M. Bradford, J. Braine, L. Ciesla, D. Clements, D. Cormier, J. A. Fernández-Ontiveros, F. Galliano, M. Giard, H. Gomez, E. González-Alfonso, F. Herpin, D. Johnstone, A. Jones, H. Kaneda, F. Kemper, V. Lebouteiller , et al. (7 additional authors not shown)

Abstract: The SPICA mid and far-infrared telescope will address fundamental issues in our understanding of star formation and ISM physics in galaxies. A particular hallmark of SPICA is the outstanding sensitivity enabled by the cold telescope, optimized detectors, and wide instantaneous bandwidth throughout the mid- and far-infrared. The spectroscopic, imaging and polarimetric observations that SPICA will b… ▽ More The SPICA mid and far-infrared telescope will address fundamental issues in our understanding of star formation and ISM physics in galaxies. A particular hallmark of SPICA is the outstanding sensitivity enabled by the cold telescope, optimized detectors, and wide instantaneous bandwidth throughout the mid- and far-infrared. The spectroscopic, imaging and polarimetric observations that SPICA will be able to collect will help in clarifying the complex physical mechanisms which underlie the baryon cycle of galaxies. In particular: (i) The access to a large suite of atomic and ionic fine-structure lines for large samples of galaxies will shed light on the origin of the observed spread in star formation rates within and between galaxies. (ii) Observations of HD rotational lines (out to $\sim$10 Mpc) and fine structure lines such as [CII] 158 $μ$m (out to $\sim$100 Mpc) will clarify the main reservoirs of interstellar matter in galaxies, including phases where CO does not emit. (iii) Far-infrared spectroscopy of dust and ice features will address uncertainties in the mass and composition of dust in galaxies, and the contributions of supernovae to the interstellar dust budget will be quantified by photometry and monitoring of supernova remnants in nearby galaxies. (iv) Observations of far-infrared cooling lines such as [OI] 63 $μ$m from star-forming molecular clouds in our Galaxy will evaluate the importance of shocks to dissipate turbulent energy. The paper concludes with requirements for the telescope and instruments, and recommendations for the observing strategy. △ Less

Submitted 30 November, 2017; originally announced November 2017.

Comments: Accepted for publication in PASA; 19 pages, 8 figures, 1 table

Unbiased large spectroscopic surveys of galaxies selected by SPICA using dust bands

Authors: H. Kaneda, D. Ishihara, S. Oyabu, M. Yamagishi, T. Wada, L. Armus, M. Baes, V. Charmandaris, B. Czerny, A. Efstathiou, J. A. Fern'andez-Ontiveros, A. Ferrara, E. Gonz'alez-Alfonso, M. Griffin, C. Gruppioni, E. Hatziminaoglou, M. Imanishi, K. Kohno, J. Kwon, T. Nakagawa, T. Onaka, F. Pozzi, D. Scott, J. -D. T. Smith, L. Spinoglio , et al. (5 additional authors not shown)

Abstract: The mid-infrared (IR) range contains many spectral features associated with large molecules and dust grains such as polycyclic aromatic hydrocarbons (PAHs) and silicates. These are usually very strong compared to fine-structure gas lines, and thus valuable in studying the spectral properties of faint distant galaxies. In this paper, we evaluate the capability of low-resolution mid-IR spectroscopic… ▽ More The mid-infrared (IR) range contains many spectral features associated with large molecules and dust grains such as polycyclic aromatic hydrocarbons (PAHs) and silicates. These are usually very strong compared to fine-structure gas lines, and thus valuable in studying the spectral properties of faint distant galaxies. In this paper, we evaluate the capability of low-resolution mid-IR spectroscopic surveys of galaxies that could be performed by SPICA. The surveys are designed to address the question how star formation and black hole accretion activities evolved over cosmic time through spectral diagnostics of the physical conditions of the interstellar/circumnuclear media in galaxies. On the basis of results obtained with Herschel far-IR photometric surveys of distant galaxies and Spitzer and AKARI near- to mid-IR spectroscopic observations of nearby galaxies, we estimate the numbers of the galaxies at redshift z > 0.5, which are expected to be detected in the PAH features or dust continuum by a wide (10 deg^2) or deep (1 deg^2) blind survey, both for a given observation time of 600 hours. As by-products of the wide blind survey, we also expect to detect debris disks, through the mid-IR excess above the photospheric emission of nearby main-sequence stars, and we estimate their number. We demonstrate that the SPICA mid-IR surveys will efficiently provide us with unprecedentedly large spectral samples, which can be studied further in the far-IR with SPICA. △ Less

Submitted 19 October, 2017; originally announced October 2017.

Comments: Paper accepted for publication on PASA on 18th October 2017, as part of the SPICA Special Issue

Feedback and feeding in the context of galaxy evolution with SPICA: direct characterization of molecular outflows and inflows

Authors: E. González-Alfonso, L. Armus, F. J. Carrera, V. Charmandaris, A. Efstathiou, E. Egami, J. A. Fernández-Ontiveros, J. Fischer, G. L. Granato, C. Gruppioni, E. Hatziminaoglou, M. Imanishi, N. Isobe, H. Kaneda, D. Koziel-Wierzbowska, M. A. Malkan, J. Martin-Pintado, S. Mateos, H. Matsuhara, G. Miniutti, T. Nakagawa, F. Pozzi, F. Rico-Villas, G. Rodighiero, P. Roelfsema , et al. (6 additional authors not shown)

Abstract: A far-infrared observatory such as the {\it SPace Infrared telescope for Cosmology and Astrophysics} ({\it SPICA}), with its unprecedented spectroscopic sensitivity, would unveil the role of feedback in galaxy evolution during the last $\sim10$ Gyr of the Universe ($z=1.5-2$), through the use of far- and mid-infrared molecular and ionic fine structure lines that trace outflowing and infalling gas.… ▽ More A far-infrared observatory such as the {\it SPace Infrared telescope for Cosmology and Astrophysics} ({\it SPICA}), with its unprecedented spectroscopic sensitivity, would unveil the role of feedback in galaxy evolution during the last $\sim10$ Gyr of the Universe ($z=1.5-2$), through the use of far- and mid-infrared molecular and ionic fine structure lines that trace outflowing and infalling gas. Outflowing gas is identified in the far-infrared through P-Cygni line shapes and absorption blueshifted wings in molecular lines with high dipolar moments, and through emission line wings of fine-structure lines of ionized gas. We quantify the detectability of galaxy-scale massive molecular and ionized outflows as a function of redshift in AGN-dominated, starburst-dominated, and main-sequence galaxies, explore the detectability of metal-rich inflows in the local Universe, and describe the most significant synergies with other current and future observatories that will measure feedback in galaxies via complementary tracers at other wavelengths. △ Less

Submitted 6 October, 2017; originally announced October 2017.

Comments: This paper belongs to the SPICA Special Issue on PASA. Accepted for publication in PASA

Tracing the evolution of dust obscured star-formation and accretion back to the reionisation epoch with SPICA

Authors: C. Gruppioni, L. Ciesla, E. Hatziminaoglou, F. Pozzi, G. Rodighiero, P. Santini, L. Armus, M. Baes, J. Braine, V. Charmandaris, D. L. Clements, N. Christopher, H. Dannerbauer, A. Efstathiou, E. Egami, J. A. Fernandez-Ontiveros, F. Fontanot, A. Franceschini, E. Gonzalez-Alfonso, M. Griffin, H. Kaneda, L. Marchetti, P. Monaco, T. Nakagawa, T. Onaka , et al. (13 additional authors not shown)

Abstract: Our current knowledge of star formation and accretion luminosity at high-redshift (z>3-4), as well as the possible connections between them, relies mostly on observations in the rest-frame ultraviolet (UV), which are strongly affected by dust obscuration. Due to the lack of sensitivity of past and current infrared (IR) instrumentation, so far it has not been possible to get a glimpse into the earl… ▽ More Our current knowledge of star formation and accretion luminosity at high-redshift (z>3-4), as well as the possible connections between them, relies mostly on observations in the rest-frame ultraviolet (UV), which are strongly affected by dust obscuration. Due to the lack of sensitivity of past and current infrared (IR) instrumentation, so far it has not been possible to get a glimpse into the early phases of the dust-obscured Universe. Among the next generation of IR observatories, SPICA, observing in the 12-350 micron range, will be the only facility that can enable us to make the required leap forward in understanding the obscured star-formation rate and black-hole accretion rate densities (SFRD and BHARD, respectively) with respect to what Spitzer and Herschel achieved in the mid- and far-IR at z<3. In particular, SPICA will have the unique ability to trace the evolution of the obscured SFRD and BHARD over cosmic time, from the peak of their activity back to the reionisation epoch (i.e., 3<z<6-7), where its predecessors had severe limitations. Here we discuss the potential of both deep and shallow photometric surveys performed with the SPICA mid-IR instrument (SMI), enabled by the very low level of impact of dust obscuration in a band centred at 34 micron. These unique unbiased photometric surveys that SPICA will perform will be followed up by observations both with the SPICA spectrometers and with other facilities at shorter and longer wavelengths, with the aim to fully characterise the evolution of AGNs and star-forming galaxies after re-ionisation. △ Less

Submitted 6 October, 2017; originally announced October 2017.

Comments: This paper belongs to the SPICA Special Issue on PASA