Showing 1–50 of 142 results for author: Braine, J

ALMA-IMF XIV: Free-Free Templates Derived from H$41α$ and Ionized Gas Content in Fifteen Massive Protoclusters

Authors: Roberto Galván-Madrid, Daniel J. Díaz-González, Frédérique Motte, Adam Ginsburg, Nichol Cunningham, Karl M. Menten, Mélanie Armante, Mélisse Bonfand, Jonathan Braine, Timea Csengeri, Pierre Dell'Ova, Fabien Louvet, Thomas Nony, Rudy Rivera-Soto, Patricio Sanhueza, Amelia M. Stutz, Friedrich Wyrowski, Rodrigo H. Álvarez-Gutiérrez, Tapas Baug, Sylvain Bontemps, Leonardo Bronfman, Manuel Fernández-López, Antoine Gusdorf, Atanu Koley, Hong-Li Liu , et al. (3 additional authors not shown)

Abstract: We use the H$41α$ recombination line to create templates of the millimeter free-free emission in the ALMA-IMF continuum maps, which allows to separate it from dust emission. This method complements spectral-index information and extrapolation from centimeter wavelength maps. We use the derived maps to estimate the properties of up to 34 HII regions across the ALMA-IMF protoclusters. The hydrogen i… ▽ More We use the H$41α$ recombination line to create templates of the millimeter free-free emission in the ALMA-IMF continuum maps, which allows to separate it from dust emission. This method complements spectral-index information and extrapolation from centimeter wavelength maps. We use the derived maps to estimate the properties of up to 34 HII regions across the ALMA-IMF protoclusters. The hydrogen ionizing-photon rate $Q_0$ and spectral types follow the evolutionary trend proposed by Motte et al. The youngest protoclusters lack detectable ionized gas, followed by protoclusters with increasing numbers of OB stars. The total $Q_0$ increases from $\sim 10^{45}$ s$^{-1}$ to $> 10^{49}$ s$^{-1}$. We used the adjacent He$41α$ line to measure the relative number abundances of helium, finding values consistent with the Galactic interstellar medium, although a few outliers are discussed. A search for sites of maser amplification of the H$41α$ line returned negative results. We looked for possible correlations between the electron densities ($n_e$), emission measures (EM), and $Q_0$ with HII region size $D$. The latter are the better correlated, with $Q_0 \propto D^{2.49\pm0.18}$. This favors interpretations where smaller ultracompact HII regions are not necessarily the less dynamically evolved versions of larger ones, but rather are ionized by less massive stars. Moderate correlations were found between dynamical width $ΔV_\mathrm{dyn}$ with $D$ and $Q_0$. $ΔV_\mathrm{dyn}$ increases from about one to two times the ionized-gas sound speed. Finally, an outlier HII region south of W43-MM2 is discussed. We suggest that this source could harbor an embedded stellar or disk wind. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Comments: Accepted to AAS Journals (ApJS). Data available at https://dataverse.harvard.edu/dataverse/h41a-freefree

Discovery of a dormant 33 solar-mass black hole in pre-release Gaia astrometry

Authors: Gaia Collaboration, P. Panuzzo, T. Mazeh, F. Arenou, B. Holl, E. Caffau, A. Jorissen, C. Babusiaux, P. Gavras, J. Sahlmann, U. Bastian, Ł. Wyrzykowski, L. Eyer, N. Leclerc, N. Bauchet, A. Bombrun, N. Mowlavi, G. M. Seabroke, D. Teyssier, E. Balbinot, A. Helmi, A. G. A. Brown, A. Vallenari, T. Prusti, J. H. J. de Bruijne , et al. (390 additional authors not shown)

Abstract: Gravitational waves from black-hole merging events have revealed a population of extra-galactic BHs residing in short-period binaries with masses that are higher than expected based on most stellar evolution models - and also higher than known stellar-origin black holes in our Galaxy. It has been proposed that those high-mass BHs are the remnants of massive metal-poor stars. Gaia astrometry is exp… ▽ More Gravitational waves from black-hole merging events have revealed a population of extra-galactic BHs residing in short-period binaries with masses that are higher than expected based on most stellar evolution models - and also higher than known stellar-origin black holes in our Galaxy. It has been proposed that those high-mass BHs are the remnants of massive metal-poor stars. Gaia astrometry is expected to uncover many Galactic wide-binary systems containing dormant BHs, which may not have been detected before. The study of this population will provide new information on the BH-mass distribution in binaries and shed light on their formation mechanisms and progenitors. As part of the validation efforts in preparation for the fourth Gaia data release (DR4), we analysed the preliminary astrometric binary solutions, obtained by the Gaia Non-Single Star pipeline, to verify their significance and to minimise false-detection rates in high-mass-function orbital solutions. The astrometric binary solution of one source, Gaia BH3, implies the presence of a 32.70 \pm 0.82 M\odot BH in a binary system with a period of 11.6 yr. Gaia radial velocities independently validate the astrometric orbit. Broad-band photometric and spectroscopic data show that the visible component is an old, very metal-poor giant of the Galactic halo, at a distance of 590 pc. The BH in the Gaia BH3 system is more massive than any other Galactic stellar-origin BH known thus far. The low metallicity of the star companion supports the scenario that metal-poor massive stars are progenitors of the high-mass BHs detected by gravitational-wave telescopes. The Galactic orbit of the system and its metallicity indicate that it might belong to the Sequoia halo substructure. Alternatively, and more plausibly, it could belong to the ED-2 stream, which likely originated from a globular cluster that had been disrupted by the Milky Way. △ Less

Submitted 19 April, 2024; v1 submitted 16 April, 2024; originally announced April 2024.

Comments: 23 pages, accepted fro publication in A&A Letters. New version with small fixes

Gaia Focused Product Release: Sources from Service Interface Function image analysis -- Half a million new sources in omega Centauri

Authors: Gaia Collaboration, K. Weingrill, A. Mints, J. Castañeda, Z. Kostrzewa-Rutkowska, M. Davidson, F. De Angeli, J. Hernández, F. Torra, M. Ramos-Lerate, C. Babusiaux, M. Biermann, C. Crowley, D. W. Evans, L. Lindegren, J. M. Martín-Fleitas, L. Palaversa, D. Ruz Mieres, K. Tisanić, A. G. A. Brown, A. Vallenari, T. Prusti, J. H. J. de Bruijne, F. Arenou, A. Barbier , et al. (378 additional authors not shown)

Abstract: Gaia's readout window strategy is challenged by very dense fields in the sky. Therefore, in addition to standard Gaia observations, full Sky Mapper (SM) images were recorded for nine selected regions in the sky. A new software pipeline exploits these Service Interface Function (SIF) images of crowded fields (CFs), making use of the availability of the full two-dimensional (2D) information. This ne… ▽ More Gaia's readout window strategy is challenged by very dense fields in the sky. Therefore, in addition to standard Gaia observations, full Sky Mapper (SM) images were recorded for nine selected regions in the sky. A new software pipeline exploits these Service Interface Function (SIF) images of crowded fields (CFs), making use of the availability of the full two-dimensional (2D) information. This new pipeline produced half a million additional Gaia sources in the region of the omega Centauri ($ω$ Cen) cluster, which are published with this Focused Product Release. We discuss the dedicated SIF CF data reduction pipeline, validate its data products, and introduce their Gaia archive table. Our aim is to improve the completeness of the {\it Gaia} source inventory in a very dense region in the sky, $ω$ Cen. An adapted version of {\it Gaia}'s Source Detection and Image Parameter Determination software located sources in the 2D SIF CF images. We validated the results by comparing them to the public {\it Gaia} DR3 catalogue and external Hubble Space Telescope data. With this Focused Product Release, 526\,587 new sources have been added to the {\it Gaia} catalogue in $ω$ Cen. Apart from positions and brightnesses, the additional catalogue contains parallaxes and proper motions, but no meaningful colour information. While SIF CF source parameters generally have a lower precision than nominal {\it Gaia} sources, in the cluster centre they increase the depth of the combined catalogue by three magnitudes and improve the source density by a factor of ten. This first SIF CF data publication already adds great value to the {\it Gaia} catalogue. It demonstrates what to expect for the fourth {\it Gaia} catalogue, which will contain additional sources for all nine SIF CF regions. △ Less

Submitted 8 November, 2023; v1 submitted 10 October, 2023; originally announced October 2023.

Journal ref: A&A 680, A35 (2023)

Gaia Focused Product Release: A catalogue of sources around quasars to search for strongly lensed quasars

Authors: Gaia Collaboration, A. Krone-Martins, C. Ducourant, L. Galluccio, L. Delchambre, I. Oreshina-Slezak, R. Teixeira, J. Braine, J. -F. Le Campion, F. Mignard, W. Roux, A. Blazere, L. Pegoraro, A. G. A. Brown, A. Vallenari, T. Prusti, J. H. J. de Bruijne, F. Arenou, C. Babusiaux, A. Barbier, M. Biermann, O. L. Creevey, D. W. Evans, L. Eyer, R. Guerra , et al. (376 additional authors not shown)

Abstract: Context. Strongly lensed quasars are fundamental sources for cosmology. The Gaia space mission covers the entire sky with the unprecedented resolution of $0.18$" in the optical, making it an ideal instrument to search for gravitational lenses down to the limiting magnitude of 21. Nevertheless, the previous Gaia Data Releases are known to be incomplete for small angular separations such as those ex… ▽ More Context. Strongly lensed quasars are fundamental sources for cosmology. The Gaia space mission covers the entire sky with the unprecedented resolution of $0.18$" in the optical, making it an ideal instrument to search for gravitational lenses down to the limiting magnitude of 21. Nevertheless, the previous Gaia Data Releases are known to be incomplete for small angular separations such as those expected for most lenses. Aims. We present the Data Processing and Analysis Consortium GravLens pipeline, which was built to analyse all Gaia detections around quasars and to cluster them into sources, thus producing a catalogue of secondary sources around each quasar. We analysed the resulting catalogue to produce scores that indicate source configurations that are compatible with strongly lensed quasars. Methods. GravLens uses the DBSCAN unsupervised clustering algorithm to detect sources around quasars. The resulting catalogue of multiplets is then analysed with several methods to identify potential gravitational lenses. We developed and applied an outlier scoring method, a comparison between the average BP and RP spectra of the components, and we also used an extremely randomised tree algorithm. These methods produce scores to identify the most probable configurations and to establish a list of lens candidates. Results. We analysed the environment of 3 760 032 quasars. A total of 4 760 920 sources, including the quasars, were found within 6" of the quasar positions. This list is given in the Gaia archive. In 87\% of cases, the quasar remains a single source, and in 501 385 cases neighbouring sources were detected. We propose a list of 381 lensed candidates, of which we identified 49 as the most promising. Beyond these candidates, the associate tables in this Focused Product Release allow the entire community to explore the unique Gaia data for strong lensing studies further. △ Less

Submitted 10 October, 2023; originally announced October 2023.

Comments: 35 pages, 60 figures, accepted for publication by Astronomy and Astrophysics

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

Gaia Focused Product Release: Radial velocity time series of long-period variables

Authors: Gaia Collaboration, Gaia Collaboration, M. Trabucchi, N. Mowlavi, T. Lebzelter, I. Lecoeur-Taibi, M. Audard, L. Eyer, P. García-Lario, P. Gavras, B. Holl, G. Jevardat de Fombelle, K. Nienartowicz, L. Rimoldini, P. Sartoretti, R. Blomme, Y. Frémat, O. Marchal, Y. Damerdji, A. G. A. Brown, A. Guerrier, P. Panuzzo, D. Katz, G. M. Seabroke, K. Benson , et al. (382 additional authors not shown)

Abstract: The third Gaia Data Release (DR3) provided photometric time series of more than 2 million long-period variable (LPV) candidates. Anticipating the publication of full radial-velocity (RV) in DR4, this Focused Product Release (FPR) provides RV time series for a selection of LPVs with high-quality observations. We describe the production and content of the Gaia catalog of LPV RV time series, and the… ▽ More The third Gaia Data Release (DR3) provided photometric time series of more than 2 million long-period variable (LPV) candidates. Anticipating the publication of full radial-velocity (RV) in DR4, this Focused Product Release (FPR) provides RV time series for a selection of LPVs with high-quality observations. We describe the production and content of the Gaia catalog of LPV RV time series, and the methods used to compute variability parameters published in the Gaia FPR. Starting from the DR3 LPVs catalog, we applied filters to construct a sample of sources with high-quality RV measurements. We modeled their RV and photometric time series to derive their periods and amplitudes, and further refined the sample by requiring compatibility between the RV period and at least one of the $G$, $G_{\rm BP}$, or $G_{\rm RP}$ photometric periods. The catalog includes RV time series and variability parameters for 9\,614 sources in the magnitude range $6\lesssim G/{\rm mag}\lesssim 14$, including a flagged top-quality subsample of 6\,093 stars whose RV periods are fully compatible with the values derived from the $G$, $G_{\rm BP}$, and $G_{\rm RP}$ photometric time series. The RV time series contain a mean of 24 measurements per source taken unevenly over a duration of about three years. We identify the great most sources (88%) as genuine LPVs, with about half of them showing a pulsation period and the other half displaying a long secondary period. The remaining 12% consists of candidate ellipsoidal binaries. Quality checks against RVs available in the literature show excellent agreement. We provide illustrative examples and cautionary remarks. The publication of RV time series for almost 10\,000 LPVs constitutes, by far, the largest such database available to date in the literature. The availability of simultaneous photometric measurements gives a unique added value to the Gaia catalog (abridged) △ Less

Submitted 9 October, 2023; originally announced October 2023.

Comments: 36 pages, 38 figures

Molecular and Ionized Gas in Tidal Dwarf Galaxies: The Spatially Resolved Star-Formation Relation

Authors: Navyasree Kovakkuni, Federico Lelli, Pierre-alain Duc, Médéric Boquien, Jonathan Braine, Elias Brinks, Vassilis Charmandaris, Francoise Combes, Jeremy Fensch, Ute Lisenfeld, Stacy McGaugh, J. Chris Mihos, Marcel. S. Pawlowski, Yves. Revaz, Peter. M. Weilbacher

Abstract: Tidal dwarf galaxies (TDGs) are low-mass objects that form within tidal and/or collisional debris ejected from more massive interacting galaxies. We use CO($1-0$) observations from ALMA and integral-field spectroscopy from MUSE to study molecular and ionized gas in three TDGs: two around the collisional galaxy NGC 5291 and one in the late-stage merger NGC 7252. The CO and H$α$ emission is more com… ▽ More Tidal dwarf galaxies (TDGs) are low-mass objects that form within tidal and/or collisional debris ejected from more massive interacting galaxies. We use CO($1-0$) observations from ALMA and integral-field spectroscopy from MUSE to study molecular and ionized gas in three TDGs: two around the collisional galaxy NGC 5291 and one in the late-stage merger NGC 7252. The CO and H$α$ emission is more compact than the HI emission and displaced from the HI dynamical center, so these gas phases cannot be used to study the internal dynamics of TDGs. We use CO, HI, and H$α$ data to measure the surface densities of molecular gas ($Σ_{\rm mol}$), atomic gas ($Σ_{\rm atom}$) and star-formation rate ($Σ_{\rm SFR}$), respectively. We confirm that TDGs follow the same spatially integrated $Σ_{\rm SFR}-Σ_{\rm gas}$ relation of regular galaxies, where $Σ_{\rm gas} = Σ_{\rm mol} + Σ_{\rm atom}$, even though they are HI dominated. We find a more complex behaviour in terms of the spatially resolved $Σ_{\rm SFR}-Σ_{\rm mol}$ relation on sub-kpc scales. The majority ($\sim$60$\%$) of SF regions in TDGs lie on the same $Σ_{\rm SFR}-Σ_{\rm mol}$ relation of normal spiral galaxies but show a higher dispersion around the mean. The remaining fraction of SF regions ($\sim$40$\%$) lie in the starburst region and are associated with the formation of massive super star clusters, as shown by Hubble Space Telescope images. We conclude that the local SF activity in TDGs proceeds in a hybrid fashion, with some regions comparable to normal spiral galaxies and others to extreme starbursts. △ Less

Submitted 12 September, 2023; originally announced September 2023.

Comments: 11 pages, 6 figures, 3 tables, Accepted for publication in MNRAS

Modeling the molecular gas content and CO-to-H2 conversion factors in low-metallicity star-forming dwarf galaxies

Authors: L. Ramambason, V. Lebouteiller, S. C. Madden, F. Galliano, C. T. Richardson, A. Saintonge, I. De Looze, M. Chevance, N. P. Abel, S. Hernandez, J. Braine

Abstract: Low-metallicity dwarf galaxies often show no or little CO emission, despite the intense star formation observed in local samples. Both simulations and resolved observations indicate that molecular gas in low-metallicity galaxies may reside in small dense clumps, surrounded by a substantial amount of more diffuse gas, not traced by CO. Constraining the relative importance of CO-bright versus CO-dar… ▽ More Low-metallicity dwarf galaxies often show no or little CO emission, despite the intense star formation observed in local samples. Both simulations and resolved observations indicate that molecular gas in low-metallicity galaxies may reside in small dense clumps, surrounded by a substantial amount of more diffuse gas, not traced by CO. Constraining the relative importance of CO-bright versus CO-dark H2 star-forming reservoirs is crucial to understand how star formation proceeds at low metallicity. We put to the test classically used single component radiative transfer models and compare their results to those obtained assuming an increasingly complex structure of the interstellar gas, mimicking an inhomogeneous distribution of clouds with various physical properties. We compute representative models of the interstellar medium as combinations of several gas components, each with a specific set of physical parameters. We introduce physically-motivated models assuming power-law distributions for the density, ionization parameter, and the depth of molecular clouds. We confirm the presence of a predominantly CO-dark molecular reservoir in low-metallicity galaxies. The predicted total H2 mass is best traced by [C II]158um and, to a lesser extent, by [CI] 609um, rather than by CO(1-0). We examine the CO-to-H2 conversion factor vs. metallicity relation and find that its dispersion increases significantly when different geometries of the gas are considered. We define a clumpiness parameter that anti-correlates with [CII]/CO and explains the dispersion of the CO-to-H2 conversion factor vs. metallicity relation. We find that low-metallicity galaxies with high clumpiness may have CO-to-H2 conversion factor as low as the Galactic value. We identify the clumpiness of molecular gas as a key parameter to understand variations of geometry-sensitive quantities, such as CO-to-H2 conversion factor. △ Less

Submitted 17 October, 2023; v1 submitted 26 June, 2023; originally announced June 2023.

Comments: Accepted for publication in A&A, 24 pages, 14 figures

ALMA-IMF. VII. First release of the full spectral line cubes: Core kinematics traced by DCN J=(3-2)

Authors: N. Cunningham, A. Ginsburg, R. Galván-Madrid, F. Motte, T. Csengeri, A. M. Stutz, M. Fernández-López, R. H. Álvarez-Gutiérrez, M. Armante, T. Baug, M. Bonfand, S. Bontemps, J. Braine, N. Brouillet, G. Busquet, D. J. Díaz-González, J. Di Francesco, A. Gusdorf, F. Herpin, H. Liu, A. López-Sepulcre, F. Louvet, X. Lu, L. Maud, T. Nony , et al. (8 additional authors not shown)

Abstract: ALMA-IMF is an Atacama Large Millimeter/submillimeter Array (ALMA) Large Program designed to measure the core mass function (CMF) of 15 protoclusters chosen to span their early evolutionary stages. It further aims to understand their kinematics, chemistry, and the impact of gas inflow, accretion, and dynamics on the CMF. We present here the first release of the ALMA-IMF line data cubes (DR1), prod… ▽ More ALMA-IMF is an Atacama Large Millimeter/submillimeter Array (ALMA) Large Program designed to measure the core mass function (CMF) of 15 protoclusters chosen to span their early evolutionary stages. It further aims to understand their kinematics, chemistry, and the impact of gas inflow, accretion, and dynamics on the CMF. We present here the first release of the ALMA-IMF line data cubes (DR1), produced from the combination of two ALMA 12m-array configurations. The data include 12 spectral windows, with eight at 1.3mm and four at 3mm. The broad spectral coverage of ALMA-IMF (~6.7 GHz bandwidth coverage per field) hosts a wealth of simple atomic, molecular, ionised, and complex organic molecular lines. We describe the line cube calibration done by ALMA and the subsequent calibration and imaging we performed. We discuss our choice of calibration parameters and optimisation of the cleaning parameters, and we demonstrate the utility and necessity of additional processing compared to the ALMA archive pipeline. As a demonstration of the scientific potential of these data, we present a first analysis of the DCN (3-2) line. We find that DCN traces a diversity of morphologies and complex velocity structures, which tend to be more filamentary and widespread in evolved regions and are more compact in the young and intermediate-stage protoclusters. Furthermore, we used the DCN (3-2) emission as a tracer of the gas associated with 595 continuum cores across the 15 protoclusters, providing the first estimates of the core systemic velocities and linewidths within the sample. We find that DCN (3-2) is detected towards a higher percentage of cores in evolved regions than the young and intermediate-stage protoclusters and is likely a more complete tracer of the core population in more evolved protoclusters. The full ALMA 12m-array cubes for the ALMA-IMF Large Program are provided with this DR1 release. △ Less

Submitted 26 June, 2023; originally announced June 2023.

Comments: 75 pages (21 main body; 54 appendix), 37 figures. The ALMA-IMF DR1 line release is hosted at https://dataverse.harvard.edu/dataverse/alma-imf

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

ViCTORIA project: MeerKAT HI observations of the ram pressure stripped galaxy NGC 4523

Authors: A. Boselli, P. Serra, F. de Gasperin, B. Vollmer, P. Amram, H. W. Edler, M. Fossati, G. Consolandi, P. Cote, J. C. Cuillandre, L. Ferrarese, S. Gwyn, J. Postma, M. Boquien, J. Braine, F. Combes, G. Gavazzi, G. Hensler, M. A. Miville-Deschenes, M. Murgia, J. Roediger, Y. Roehlly, R. Smith, H. X. Zhang, N. Zabel

Abstract: We present the first results of a 21 cm HI line pilot observation carried out with MeerKAT in preparation for the ViCTORIA project, an untargeted survey of the Virgo galaxy cluster. The extraordinary quality of the data in terms of sensitivity and angular resolution (rms~0.65 mJy beam^-1 at ~27"x39" and 11 km/s resolution) allowed us to detect an extended (~10 kpc projected length) low column dens… ▽ More We present the first results of a 21 cm HI line pilot observation carried out with MeerKAT in preparation for the ViCTORIA project, an untargeted survey of the Virgo galaxy cluster. The extraordinary quality of the data in terms of sensitivity and angular resolution (rms~0.65 mJy beam^-1 at ~27"x39" and 11 km/s resolution) allowed us to detect an extended (~10 kpc projected length) low column density (N(HI) < 2.5x10^20 cm^-2) HI gas tail associated with the dwarf irregular galaxy NGC4523 at the northern edge of the cluster. The morphology of the tail and of the stellar disc suggest that the galaxy is suffering a hydrodynamic interaction with the surrounding hot intracluster medium (ICM; ram pressure stripping). The orientation of the trailing tail, the gradient in the HI gas column density at the interface between the cold ISM and the hot ICM, the velocity of the galaxy with respect to that of the cluster, and its position indicate that NGC4523 is infalling for the first time into Virgo from the NNW background of the cluster. Using a grid of hydrodynamic simulations we derive the impact parameters with the surrounding ICM, and estimate that the galaxy will be at pericentre (D~500-600 kpc) in ~1 Gyr, where ram pressure stripping will be able to remove most, if not all, of its gas. The galaxy is located on the star formation main sequence when its star formation rate is derived using Halpha images obtained during the VESTIGE survey, suggesting that NGC4523 is only at the beginning of its interaction with the surrounding environment. A few HII regions are detected in the Halpha images within the HI gas tail outside the stellar disc. Their ages, derived by comparing their Halpha, FUV, NUV, and optical colours with the predictions of SED fitting models, are <30 Myr, and suggest that these HII regions have formed within the stripped gas. △ Less

Submitted 22 June, 2023; originally announced June 2023.

Comments: Accepted for publication in A&A

Journal ref: A&A 676, A92 (2023)

Asymmetries in random motions of neutral Hydrogen gas in spiral galaxies

Authors: P. Adamczyk, P. Amram, L. Chemin, B. Epinat, J. Braine, F. Combes, W. G. J. de Blok

Abstract: (Abridged). It has been recently shown that random motions of the neutral Hydrogen gas of the Triangulum galaxy (M33) exhibit a bisymmetric perturbation which is aligned with the minor axis of the galaxy, suggesting a projection effect. To investigate if perturbations in the velocity dispersion of nearby discs are comparable to those of M33, the sample is extended to 32 galaxies from The HI Nearby… ▽ More (Abridged). It has been recently shown that random motions of the neutral Hydrogen gas of the Triangulum galaxy (M33) exhibit a bisymmetric perturbation which is aligned with the minor axis of the galaxy, suggesting a projection effect. To investigate if perturbations in the velocity dispersion of nearby discs are comparable to those of M33, the sample is extended to 32 galaxies from The HI Nearby Galaxy Survey and the Westerbork HI Survey of Spiral and Irregular Galaxies. We study velocity asymmetries in the disc planes by performing Fourier transforms of high-resolution HI velocity dispersion maps corrected for beam smearing effects, and measure the amplitudes and phase angles of the Fourier harmonics. We find strong perturbations of first, second and fourth orders. The strongest asymmetry is the bisymmetry, which is predominantly associated with the presence of spiral arms. The first order asymmetry is generally oriented close to the disc major axis, and the second and fourth order asymmetries are preferentially oriented along intermediate directions between the major and minor axes of the discs. These results are evidence that strong projection effects shape the HI velocity dispersion maps. The most likely source of systematic orientations is the anisotropy of velocities, through the projection of streaming motions stronger along one of the planar directions in the discs. Moreover, systematic phase angles of asymmetries in the HI velocity dispersion could arise from tilted velocity ellipsoids. We expect a larger incidence of correlation between the radial and tangential velocities of HI gas. Our methodology is a powerful tool to constrain the dominant direction of streaming motions and thus the shape of the velocity ellipsoid of HI gas, which is de facto anisotropic at the angular scales probed by the observations. △ Less

Submitted 12 June, 2023; originally announced June 2023.

Comments: 40 pages, 33 figures. Accepted for publication in Astronomy & Astrophysics. Full resolution version available upon request

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

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

Clusters, Clouds, and Correlations: Relating Young Clusters to Giant Molecular Clouds in M33 and M31

Authors: Joshua Peltonen, Erik Rosolowsky, L. Clifton Johnson, Anil C. Seth, Julianne Dalcanton, Eric F. Bell, Jonathan Braine, Eric W. Koch, Margaret Lazzarini, Adam K. Leroy, Evan D. Skillman, Adam Smercina, Tobin Wainer, Benjamin F. Williams

Abstract: We use young clusters and giant molecular clouds (GMCs) in the galaxies M33 and M31 to constrain temporal and spatial scales in the star formation process. In M33, we compare the PHATTER catalogue of 1214 clusters with ages measured via colour-magnitude diagram (CMD) fitting to 444 GMCs identified from a new 35 pc resolution ALMA $^{12}$CO(2-1) survey. In M31, we compare the PHAT catalogue of 1249… ▽ More We use young clusters and giant molecular clouds (GMCs) in the galaxies M33 and M31 to constrain temporal and spatial scales in the star formation process. In M33, we compare the PHATTER catalogue of 1214 clusters with ages measured via colour-magnitude diagram (CMD) fitting to 444 GMCs identified from a new 35 pc resolution ALMA $^{12}$CO(2-1) survey. In M31, we compare the PHAT catalogue of 1249 clusters to 251 GMCs measured from a CARMA $^{12}$CO(1-0) survey with 20 pc resolution. Through two-point correlation analysis, we find that young clusters have a high probability of being near other young clusters, but correlation between GMCs is suppressed by the cloud identification algorithm. By comparing the positions, we find that younger clusters are closer to GMCs than older clusters. Through cross-correlation analysis of the M33 cluster data, we find that clusters are statistically associated when they are $\leq$10 Myr old. Utilizing the high precision ages of the clusters, we find that clusters older than $\approx 18$ Myr are uncorrelated with the molecular ISM. Using the spatial coincidence of the youngest clusters and GMCs in M33, we estimate that clusters spend $\approx$4-6 Myr inside their parent GMC. Through similar analysis, we find that the GMCs in M33 have a total lifetime of $\approx 11$-15 Myr. We also develop a drift model and show that the above correlations can be explained if the clusters in M33 have a 5-10 km s$^{-1}$ velocity dispersion relative to the molecular ISM. △ Less

Submitted 5 May, 2023; originally announced May 2023.

Comments: 14 pages, 13 figures, 1 tables, accepted for publication at MNRAS

A Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE).XIV. The main sequence relation in a rich environment down to M_star ~ 10^6 Mo

Authors: A. Boselli, M. Fossati, J. Roediger, M. Boquien, M. Fumagalli, M. Balogh, S. Boissier, J. Braine, L. Ciesla, P. Côté, J. C. Cuillandre, L. Ferrarese, G. Gavazzi, S. Gwyn, Junais, G. Hensler, A. Longobardi, M. Sun

Abstract: Using a compilation of Halpha fluxes for 384 star forming galaxies detected during the VESTIGE survey, we study several important scaling relations for a complete sample of galaxies in a rich environment. The extraordinary sensitivity of the data allows us to sample the whole dynamic range of the Halpha luminosity function, from massive (M*~10^11 Mo) to dwarf systems (M*~10^6 Mo). This extends pre… ▽ More Using a compilation of Halpha fluxes for 384 star forming galaxies detected during the VESTIGE survey, we study several important scaling relations for a complete sample of galaxies in a rich environment. The extraordinary sensitivity of the data allows us to sample the whole dynamic range of the Halpha luminosity function, from massive (M*~10^11 Mo) to dwarf systems (M*~10^6 Mo). This extends previous works to a dynamic range in stellar mass and star formation rate (10^-4<SFR<10 Mo yr^-1) never explored so far. The main sequence (MS) relation derived for all star forming galaxies within one virial radius of the Virgo cluster has a slope comparable to that observed in other nearby samples of isolated objects, but has a dispersion ~3 times larger. The dispersion is tightly connected to the available amount of HI gas, with gas-poor systems located far below objects of similar stellar mass but with a normal HI content. When measured on unperturbed galaxies with a normal HI gas content, the relation has a slope a=0.92, an intercept b=-1.57, and a scatter ~0.40. We compare these observational results to the prediction of models. The observed scatter in the MS relation can be reproduced only after a violent and active stripping process such as ram-pressure that removes gas from the disc and quenches star formation on short (<1 Gyr) timescales. This rules out milder processes such as starvation. This interpretation is also consistent with the position of galaxies of different star formation activity and gas content within the phase-space diagram. We also show that the star forming regions formed in the stripped material outside perturbed galaxies are located well above the MS relation drawn by unperturbed systems. These HII regions, which might be at the origin of compact sources typical in rich environments, are living a starburst phase lasting only <50 Myr, later becoming quiescent systems. △ Less

Submitted 3 November, 2022; originally announced November 2022.

Comments: Accepted for publication on A&A

Journal ref: A&A 669, A73 (2023)

Cloud-scale Radio Surveys of Star Formation and Feedback in Triangulum Galaxy M 33: VLA Observations

Authors: F. S. Tabatabaei, W. Cotton, E. Schinnerer, R. Beck, A. Brunthaler, K. M. Menten, J. Braine, E. Corbelli, C. Kramer, J. E. Beckman, J. H. Knapen, R. Paladino, E. Koch, A. Camps Farina

Abstract: Studying the interplay between massive star formation and the interstellar medium (ISM) is paramount to understand the evolution of galaxies. Radio continuum (RC) emission serves as an extinction-free tracer of both massive star formation and the energetic components of the interstellar medium. We present a multi-band radio continuum survey of the local group galaxy M 33 down to ~30 pc linear reso… ▽ More Studying the interplay between massive star formation and the interstellar medium (ISM) is paramount to understand the evolution of galaxies. Radio continuum (RC) emission serves as an extinction-free tracer of both massive star formation and the energetic components of the interstellar medium. We present a multi-band radio continuum survey of the local group galaxy M 33 down to ~30 pc linear resolution observed with the Karl G. Jansky Very Large Array (VLA). We calibrate the star-formation rate surface density and investigate the impact of diffuse emission on this calibration using a structural decomposition. Separating the thermal and nonthermal emission components, the correlation between different phases of the interstellar medium and the impact of massive star formation are also investigated. Radio sources with sizes <~ 200 pc constitute about 36% (46%) of the total RC emission at 1.5 GHz (6.3 GHz) in the inner 18' x 18' (or 4kpc x 4kpc) disk of M 33. The nonthermal spectral index becomes flatter with increasing star-formation rate surface density, indicating the escape of cosmic ray electrons {from their birth places}. The magnetic field strength also increases with star-formation rate following a bi-modal relation, indicating that the small-scale turbulent dynamo acts more efficiently at higher luminosities and star-formation rates. Although the correlations are tighter in star-forming regions, the nonthermal emission is correlated also with the more quiescent molecular gas in the ISM. An almost linear molecular star-formation law exists in M 33 when excluding diffuse structures. Massive star formation amplifies the magnetic field and increases the number of high-energy cosmic ray electrons, which can help the onset of winds and outflows. △ Less

Submitted 3 September, 2022; originally announced September 2022.

Comments: 20 pages, 12 figures, accepted for publication in MNRAS journal

Predicting HCN, HCO + , multi-transition CO, and dust emission of star-forming galaxies: Constraining the properties of resolved gas and dust disks of local spiral galaxies

Authors: T. Lizée, B. Vollmer, J. Braine, P. Gratier, F. Bigiel

Abstract: The ISM is a turbulent, multi-phase, and multi-scale medium following scaling relations. Analytical models of galactic gaseous disks need to take into account the multi-scale and multi-phase nature of the interstellar medium. They can be described as clumpy star-forming accretion disks in vertical hydrostatic equilibrium, with the mid-plane pressure balancing the gravity of the gaseous and stellar… ▽ More The ISM is a turbulent, multi-phase, and multi-scale medium following scaling relations. Analytical models of galactic gaseous disks need to take into account the multi-scale and multi-phase nature of the interstellar medium. They can be described as clumpy star-forming accretion disks in vertical hydrostatic equilibrium, with the mid-plane pressure balancing the gravity of the gaseous and stellar disk. ISM turbulence is taken into account by applying Galactic scaling relations to the cold atomic and molecular gas phases. Turbulence is maintained through energy injection by supernovae. With the determination of the gas mass fraction at a given spatial scale, the equilibrium gas temperature between turbulent heating and line cooling, the molecular abundances, and the molecular line emission can be calculated. The resulting model radial profiles of IR, H{\sc i}, CO, HCN, and HCO$^+$ emission are compared to THINGS, HERACLES, EMPIRE, SINGS, and GALEX observations of 17 local spiral galaxies. The Toomre parameter, which measures the stability against star formation (cloud collapse), exceeds unity in the inner disk of a significant number of galaxies. In two galaxies it also exceeds unity in the outer disk. Therefore, in spiral galaxies $Q_{\rm tot}=1$ is not ubiquitous. The model gas velocity dispersion is consistent with the observed H{\sc i} velocity dispersion where available. Within our model HCN and HCO$^+$ is already detectable in relatively low-density gas ($\sim 1000$~cm$^{-3}$). CO and HCN conversion factors and molecular gas depletion time were derived. Both conversion factors are consistent with values found in the literature. Whereas in the massive galaxies the viscous timescale greatly exceeds the star formation timescale, the viscous timescale is smaller than the star formation timescale within $\rm{R}~\sim~2~\rm{R}_{\rm d}$, the disk scale length, in the low-mass galaxies. △ Less

Submitted 15 May, 2022; originally announced May 2022.

Comments: 22 pages, 10 figures (without appendix) 44 pages, 27 figures

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

The CO emission in the Taffy Galaxies (UGC 12914/5) at 60pc resolution-I: The battle for star formation in the turbulent Taffy Bridge

Authors: P. N. Appleton, B. Emonts, U. Lisenfeld, E. Falgarone, P. Guillard, F. Boulanger, J. Braine, P. Ogle, C. Struck, B. Vollmer, T. Yeager

Abstract: We present ALMA observations at a spatial resolution of 0.2 arcsec (60 pc) of CO emission from the Taffy galaxies (UGC 12914/5). The observations are compared with narrow-band Pa$α$, mid-IR, radio continuum and X-ray imaging, plus optical spectroscopy. The galaxies have undergone a recent head-on collision, creating a massive gaseous bridge which is known to be highly turbulent. The bridge contain… ▽ More We present ALMA observations at a spatial resolution of 0.2 arcsec (60 pc) of CO emission from the Taffy galaxies (UGC 12914/5). The observations are compared with narrow-band Pa$α$, mid-IR, radio continuum and X-ray imaging, plus optical spectroscopy. The galaxies have undergone a recent head-on collision, creating a massive gaseous bridge which is known to be highly turbulent. The bridge contains a complex web of narrow molecular filaments and clumps. The majority of the filaments are devoid of star formation, and fall significantly below the Kennicutt-Schmidt relationship for normal galaxies, especially for the numerous regions undetected in Pa$α$ emission. Within the loosely connected filaments and clumps of gas we find regions of high velocity dispersion which appear gravitationally unbound for a wide range of likely values of $X_{\rm CO}$. Like the "Firecracker" region in the Antennae system, they would require extremely high external dynamical or thermal pressure to stop them dissipating rapidly on short crossing timescales of 2-5~Myrs. We suggest that the clouds may be transient structures within a highly turbulent multi-phase medium which is strongly suppressing star formation. Despite the overall turbulence in the system, stars seem to have formed in compact hotspots within a kpc-sized extragalactic HII region, where the molecular gas has a lower velocity dispersion than elsewhere, and shows evidence for a collision with an ionized gas cloud. Like the shocked gas in the Stephan's Quintet group, the conditions in the Taffy bridge shows how difficult it is to form stars within a turbulent, multi-phase, gas. △ Less

Submitted 31 March, 2022; originally announced March 2022.

Comments: 23 pages, 1 Appendix, 17 figures

ALMA-IMF III -- Investigating the origin of stellar masses: Top-heavy core mass function in the W43-MM2&MM3 mini-starburst

Authors: Y. Pouteau, F. Motte, T. Nony, R. Galván-Madrid, A. Men'shchikov, S. Bontemps, J. -F. Robitaille, F. Louvet, A. Ginsburg, F. Herpin, A. López-Sepulcre, P. Dell'Ova, A. Gusdorf, P. Sanhueza, A. M. Stutz, N. Brouillet, B. Thomasson, M. Armante, T. Baug, G. Busquet, T. Csengeri, N. Cunningham, M. Fernández-López, H. -L. Liu, F. Olguin , et al. (13 additional authors not shown)

Abstract: The ALMA-IMF Large Program observed the W43-MM2-MM3 ridge, whose 1.3mm and 3mm ALMA 12m array continuum images reach a 2500au spatial resolution. We used both the best-sensitivity and the line-free ALMA-IMF images, reduced the noise with the multi-resolution segmentation technique MnGSeg, and derived the most complete and most robust core catalog possible. Using two different extraction software p… ▽ More The ALMA-IMF Large Program observed the W43-MM2-MM3 ridge, whose 1.3mm and 3mm ALMA 12m array continuum images reach a 2500au spatial resolution. We used both the best-sensitivity and the line-free ALMA-IMF images, reduced the noise with the multi-resolution segmentation technique MnGSeg, and derived the most complete and most robust core catalog possible. Using two different extraction software packages, getsf and GExt2D, we identified 200 compact sources, whose 100 common sources have on average fluxes consistent to within 30%. We filtered sources with non-negligible free-free contamination and corrected fluxes from line contamination, resulting in a W43-MM2-MM3 catalog of 205 getsf cores. With a median deconvolved FWHM size of 3400au, core masses range from 0.1Msun to 70Msun and the getsf catalog is 90% complete down to 0.8Msun. The high-mass end of the core mass function (CMF) of W43-MM2-MM3 is top-heavy compared to the canonical IMF. Fitting the cumulative CMF with a single power law of the form N(>logM)\propto M^a, we measured a=-0.95\pm0.04, compared to the canonical a=-1.35 Salpeter IMF slope. The slope of the CMF is robust with respect to map processing, extraction software package, and reasonable variations in the assumptions taken to estimate core masses. We explore several assumptions on how cores transfer their mass to stars and sub-fragment to predict the IMF resulting from the W43-MM2-MM3 CMF. In stark contrast to the commonly accepted paradigm, our result argues against the universality of the CMF shape. More robust functions of the star-formation efficiency and core sub-fragmentation are required to better predict the resulting IMF, here suggested to remain top-heavy at the end of the star-formation phase. If confirmed, the IMFs emerging from starburst events could inherit their top-heavy shape from their parental CMFs, challenging the IMF universality. △ Less

Submitted 28 April, 2022; v1 submitted 7 March, 2022; originally announced March 2022.

Comments: 28 pages, 16 figures, 6 tables, A&A accepted on 02/08/22, revised following the language referee report

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

ALMA-IMF II -- investigating the origin of stellar masses: Continuum Images and Data Processing

Authors: A. Ginsburg, T. Csengeri, R. Galván-Madrid, N. Cunningham, R. H. Álvarez-Gutiérrez, T. Baug, M. Bonfand, S. Bontemps, G. Busquet, D. J. Díaz-González, M. Fernández-López, A. Guzmán, F. Herpin, H. Liu, A. López-Sepulcre, F. Louvet, L. Maud, F. Motte, F. Nakamura, T. Nony, F. A. Olguin, Y. Pouteau, P. Sanhueza, A. M. Stutz, A. P. M. Towner , et al. (27 additional authors not shown)

Abstract: We present the first data release of the ALMA-IMF Large Program, which covers the 12m-array continuum calibration and imaging. The ALMA-IMF Large Program is a survey of fifteen dense molecular cloud regions spanning a range of evolutionary stages that aims to measure the core mass function (CMF). We describe the data acquisition and calibration done by the Atacama Large Millimeter/submillimeter Ar… ▽ More We present the first data release of the ALMA-IMF Large Program, which covers the 12m-array continuum calibration and imaging. The ALMA-IMF Large Program is a survey of fifteen dense molecular cloud regions spanning a range of evolutionary stages that aims to measure the core mass function (CMF). We describe the data acquisition and calibration done by the Atacama Large Millimeter/submillimeter Array (ALMA) observatory and the subsequent calibration and imaging we performed. The image products are combinations of multiple 12m array configurations created from a selection of the observed bandwidth using multi-term, multi-frequency synthesis imaging and deconvolution. The data products are self-calibrated and exhibit substantial noise improvements over the images produced from the delivered data. We compare different choices of continuum selection, calibration parameters, and image weighting parameters, demonstrating the utility and necessity of our additional processing work. Two variants of continuum selection are used and will be distributed: the "best-sensitivity" data, which include the full bandwidth, including bright emission lines that contaminate the continuum, and "cleanest", which select portions of the spectrum that are unaffected by line emission. We present a preliminary analysis of the spectral indices of the continuum data, showing that the ALMA products are able to clearly distinguish free-free emission from dust emission, and that in some cases we are able to identify optically thick emission sources. The data products are made public with this release. △ Less

Submitted 13 May, 2023; v1 submitted 15 December, 2021; originally announced December 2021.

Comments: Data released on Zenodo at https://zenodo.org/record/5702966. Accepted for publication in A&A

Journal ref: A&A 662, A9 (2022)

ALMA-IMF I -- Investigating the origin of stellar masses: Introduction to the Large Program and first results

Authors: F. Motte, S. Bontemps, T. Csengeri, Y. Pouteau, F. Louvet, A. M. Stutz, N. Cunningham, A. López-Sepulcre, N. Brouillet, R. Galván-Madrid, A. Ginsburg, L. Maud, A. Men'shchikov, F. Nakamura, T. Nony, P. Sanhueza, R. H. Álvarez-Gutiérrez, M. Armante, T. Baug, M. Bonfand, G. Busquet, E. Chapillon, D. Díaz-González, M. Fernández-López, A. E. Guzmán , et al. (39 additional authors not shown)

Abstract: The ALMA-IMF Large Program imaged a total noncontiguous area of 53pc2, covering 15 extreme, nearby protoclusters of the Milky Way. They were selected to span relevant early protocluster evolutionary stages. Our 1.3mm and 3mm observations provide continuum images that are homogeneously sensitive to point-like cores with masses of 0.2 and 0.6Msun, respectively, with a matched spatial resolution of 2… ▽ More The ALMA-IMF Large Program imaged a total noncontiguous area of 53pc2, covering 15 extreme, nearby protoclusters of the Milky Way. They were selected to span relevant early protocluster evolutionary stages. Our 1.3mm and 3mm observations provide continuum images that are homogeneously sensitive to point-like cores with masses of 0.2 and 0.6Msun, respectively, with a matched spatial resolution of 2000au. We also detect lines that probe the protocluster structure, kinematics, chemistry, and feedback over scales from clouds to filaments to cores. We classify ALMA-IMF protoclusters as Young, Intermediate, or Evolved based on the amount of dense gas in the cloud that has potentially been impacted by HII regions. The ALMA-IMF catalog contains 700 cores that span a mass range of 0.15-250Msun at a typical size of 2100au. We show that this core sample has no significant distance bias and can be used to build core mass functions at similar physical scales. Significant gas motions, which we highlight here in the G353.41 region, are traced down to core scales and can be used to look for inflowing gas streamers and to quantify the impact of the possible associated core mass growth on the shape of the CMF with time. Our first analysis does not reveal any significant evolution of the matter concentration from clouds to cores or from the youngest to more evolved protoclusters, indicating that cloud dynamical evolution and stellar feedback have for the moment only had a slight effect on the structure of high-density gas in our sample. Furthermore, the first-look analysis of the line richness toward bright cores indicates that the survey encompasses several tens of hot cores, of which we highlight the most massive in the G351.77 cloud. Their homogeneous characterization can be used to constrain the emerging molecular complexity in protostars of high to intermediate masses. △ Less

Submitted 15 December, 2021; originally announced December 2021.

Comments: 27 pages, 9 figures, accepted for publication by A&A

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

The structure and characteristic scales of the HI gas in galactic disks

Authors: Sami Dib, Jonathan Braine, Maheswar Gopinathan, Maritza A. Lara-López, Valery V. Kravtsov, Archana Soam, Ekta Sharma, Svitlana Zhukovska, Charles Aouad, José-Antonio Belinchón, George Helou, Di Li

Abstract: The spatial distribution of the HI gas in galaxies holds important clues on the physical processes that shape the structure and dynamics of the interstellar medium (ISM). In this work, we quantify the structure of the HI gas in a sample of 33 nearby galaxies taken from the THINGS Survey using the delta-variance spectrum. The THINGS galaxies display a large diversity in their spectra, however, ther… ▽ More The spatial distribution of the HI gas in galaxies holds important clues on the physical processes that shape the structure and dynamics of the interstellar medium (ISM). In this work, we quantify the structure of the HI gas in a sample of 33 nearby galaxies taken from the THINGS Survey using the delta-variance spectrum. The THINGS galaxies display a large diversity in their spectra, however, there are a number of recurrent features. In many galaxies, we observe a bump in the spectrum on scales of a few to several hundred pc. We find the characteristic scales associated with the bump to be correlated with galactic SFR for values of the SFR > 0.5 M$_{sol}$ yr$^{-1}$ and also with the median size of the HI shells detected in those galaxies. On larger scales, we observe the existence of two self-similar regimes. The first one, on intermediate scales is shallow and the power law that describes this regime has an exponent in the range [0.1-1] with a mean value of 0.55 which is compatible with the density field being generated by supersonic turbulence in the cold phase of the HI gas. The second power law is steeper, with a range of exponents between [0.5-1.5] and a mean value of 1.5. These values are associated with subsonic turbulence which is characteristic of the warm phase of the HI gas. The spatial scale at which the transition between the two regimes occurs is found to be $\approx 0.5 R_{25}$ which is similar to the size of the molecular disk in the THINGS galaxies. Overall, our results suggest that on scales < $0.5 R_{25}$, the structure of the ISM is affected by the effects of supernova explosions. On larger scales (> 0.5 $R_{25}$), stellar feedback has no significant impact, and the structure of the ISM is determined by large scale processes that govern the dynamics of the gas in the warm neutral medium such as the flaring of the HI disk and the effects of ram pressure stripping. △ Less

Submitted 8 October, 2021; v1 submitted 28 July, 2021; originally announced July 2021.

Comments: Accepted to A&A. Matches the published version

Journal ref: A&A 655, A101 (2021)

Low star formation efficiency due to turbulent adiabatic compression in the Taffy bridge

Authors: B. Vollmer, J. Braine, B. Mazzilli-Ciraulo, B. Schneider

Abstract: The Taffy system (UGC 12914/15) consists of two massive spiral galaxies which had a head-on collision about 20 Myr ago. New sensitive, high-resolution CO(1-0) observations of the Taffy system with the IRAM PdBI are presented. About 25% of the total interferometric CO luminosity stems from the bridge region. Assuming a Galactic N(H2)/ICO conversion factor for the galactic disks and a third of this… ▽ More The Taffy system (UGC 12914/15) consists of two massive spiral galaxies which had a head-on collision about 20 Myr ago. New sensitive, high-resolution CO(1-0) observations of the Taffy system with the IRAM PdBI are presented. About 25% of the total interferometric CO luminosity stems from the bridge region. Assuming a Galactic N(H2)/ICO conversion factor for the galactic disks and a third of this value for the bridge gas, about 10% of the molecular gas mass is located in the bridge region. The giant HII region close to UGC 12915 is located at the northern edge of the high-surface brightness giant molecular cloud association (GMA), which has the highest velocity dispersion among the bridge GMAs. The bridge GMAs are clearly not virialized because of their high velocity dispersion. Three dynamical models are presented and while no single model reproduces all of the observed features, they are all present in at least one of the models. Most of the bridge gas detected in CO does not form stars. We suggest that turbulent adiabatic compression is responsible for the exceptionally high velocity dispersion of the molecular ISM and the suppression of star formation in the Taffy bridge. In this scenario the turbulent velocity dispersion of the largest eddies and turbulent substructures/clouds increase such that giant molecular clouds are no longer in global virial equilibrium. The increase of the virial parameter leads to a decrease of the star formation efficiency. Most of the low-surface density, CO-emitting gas will disperse without forming stars but some of the high-density gas will probably collapse and form dense star clusters, such as the luminous HII region close to UGC 12915. We suggest that globular clusters and super star clusters formed and still form through the gravitational collapse of gas previously compressed by turbulent adiabatic compression during galaxy interactions. △ Less

Submitted 22 January, 2021; v1 submitted 18 January, 2021; originally announced January 2021.

Comments: accepted for publication in A&A, corrected Observations section

Gas compression and stellar feedback in the tidally interacting and ram-pressure stripped Virgo spiral galaxy NGC 4654

Authors: T. Lizée, B. Vollmer, J. Braine, F. Nehlig

Abstract: NGC 4654 is a Virgo galaxy seen almost face-on, which undergoes nearly edge-on gas ram pressure stripping and a fly-by gravitational interaction with another massive galaxy, NGC 4639. NGC 4654 shows a strongly compressed gas region near the outer edge of the optical disk, with HI surface densities (HSDR), exceeding the canonical value of 10-15 Msun/pc2. New IRAM 30m HERA CO(2-1) data of NGC 4654 a… ▽ More NGC 4654 is a Virgo galaxy seen almost face-on, which undergoes nearly edge-on gas ram pressure stripping and a fly-by gravitational interaction with another massive galaxy, NGC 4639. NGC 4654 shows a strongly compressed gas region near the outer edge of the optical disk, with HI surface densities (HSDR), exceeding the canonical value of 10-15 Msun/pc2. New IRAM 30m HERA CO(2-1) data of NGC 4654 are used to study the physical conditions of the ISM. The CO-to-H$_2$ conversion factor was estimated and found to be one to two times the Galactic value with significant decrease in the ratio between the molecular fraction and the total ISM pressure in the HSDR, self-gravitating gas, a Toomre parameter below $Q=1$ and star-formation efficiency 1.5-2 times higher. Analytical models were used to reproduce radial profiles of the SFR and the atomic and molecular surface densities. A Toomre parameter of $\rm Q \sim 0.8$ combined with an increase in the velocity dispersion of 5 km/s are necessary conditions to simultaneously reproduce the gas surface densities and the SFR. A dynamical model was used to reproduce the gas distribution of NGC 4654. The comparison between the velocity dispersion given by the moment 2 map and the intrinsic 3D velocity dispersion from the model were used to discriminate between regions of broader linewidths caused by a real increase in the velocity dispersion and those caused by an unresolved velocity gradient only. We found that the 5 km/s increase in the intrinsic velocity dispersion is compatible with observations. During a period of gas compression through external interactions, the gas surface density is enhanced, leading to an increased SFR and stellar feedback. Under the influence of stellar feedback, the gas density increases only moderately. The stellar feedback acts as a regulator of star-formation, increasing the turbulent velocity within the region. △ Less

Submitted 20 November, 2020; originally announced November 2020.

Comments: 27 pages, 48 figures, 4 tables, accepted for publication in A&A

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

A Virgo Environmental Survey Tracing Ionised Gas Emission. VESTIGE VIII. Bridging the cluster-ICM-galaxy evolution at small scales

Authors: A. Longobardi, A. Boselli, M. Fossati, J. A. Villa-Vélez, S. Bianchi, V. Casasola, E. Sarpa, F. Combes, G. Hensler8, D. Burgarella, C. Schimd, A. Nanni, P. Côté, V. Buat1, P. Amram, L. Ferrarese, J. Braine, G. Trinchieri, S. Boissier, M. Boquien, P. Andreani, S. Gwyn, J. C. Cuillandre

Abstract: We measure FIR emission from tails of stripped dust following the ionised and atomic gas components in galaxies undergoing ram pressure stripping. We study the dust-to-gas relative distribution and mass ratio in the stripped interstellar medium and relate them to those of the intra-cluster medium, thus linking the cluster-ICM-galaxy evolution at small-scales. The galaxy sample consists of three Sc… ▽ More We measure FIR emission from tails of stripped dust following the ionised and atomic gas components in galaxies undergoing ram pressure stripping. We study the dust-to-gas relative distribution and mass ratio in the stripped interstellar medium and relate them to those of the intra-cluster medium, thus linking the cluster-ICM-galaxy evolution at small-scales. The galaxy sample consists of three Scd Virgo galaxies with stellar masses in the range $10^9\lesssim \mathrm{M_{*}} \lesssim 10^{10}\, \mathrm{M_{\odot}}$, and within 1 Mpc from the cluster centre, namely NGC 4330, NGC 4522, and NGC 4654. Through the analysis of VESTIGE H$α$, $Herschel$ SPIRE far-infrared, and VIVA HI data, we trace the spatial distribution of the tails and infer the dust and gas masses from the measured far-infrared 250 $μ$m and HI flux densities. Dust-to-gas mass ratios (DGRs) in the tails are analysed as a function of the galaxy mass, metallicity, and dust temperature. Along the stripped component, the dust distribution closely follows the HI and H$α$ emitting gas, all extending beyond the optical disc. In these regions, the DGRs are $2.0\pm0.6\times10^{-3}$, $0.7\pm0.1\times10^{-3}$, and $0.4\pm0.03\times10^{-3}$, for NGC 4330, NGC 4522, and NGC 4654, respectively, i.e. up to a factor of 15 less than the values measured in the main body of nearby galaxies. We also find a negative trend in the DGR as a function of the metallicity that can be explained in terms of a dust component more centrally concentrated in more metal-rich systems. Together with the finding that the stripped dust is cold, $T_{d} \lesssim 25\, K$, our results support an outside-in stripping scenario of the galaxy interstellar medium. This study shows that ram pressure stripping is a key mechanism in the building up of the Virgo intra-cluster component injecting dust grains into the ICM, thus contributing to its metal enrichment. △ Less

Submitted 5 October, 2020; originally announced October 2020.

Comments: 15 pages, 7 images, 4 tables. Accepted on A\&A

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

Anisotropy of random motions of gas in Messier 33

Authors: Laurent Chemin, Jonathan Braine, Françoise Combes, Zacharie S. Kam, Claude Carignan

Abstract: (Abridged) We study the properties of anisotropic and axisymmetric velocity ellipsoids from maps of the gas velocity dispersion in nearby galaxies. This data allow us to measure the azimuthal-to-radial axis ratio of gas velocity ellipsoids, which is a useful tool to study the structure of gaseous orbits in the disk. We also present the first estimates of perturbations in gas velocity dispersion ma… ▽ More (Abridged) We study the properties of anisotropic and axisymmetric velocity ellipsoids from maps of the gas velocity dispersion in nearby galaxies. This data allow us to measure the azimuthal-to-radial axis ratio of gas velocity ellipsoids, which is a useful tool to study the structure of gaseous orbits in the disk. We also present the first estimates of perturbations in gas velocity dispersion maps by applying an alternative model that considers isotropic and asymmetric random motions. High-quality velocity dispersion maps of the atomic medium at various angular resolutions of the nearby spiral galaxy Messier 33, are used to test the anisotropic and isotropic velocity models. The velocity dispersions of hundreds of individual molecular clouds are also analyzed. The HI velocity dispersion of M33 is systematically larger along the minor axis, and lower along the major axis. Isotropy is only possible if asymmetric motions are considered. Fourier transforms of the HI velocity dispersions reveal a bisymmetric mode which is mostly stronger than other asymmetric motions and aligned with the minor axis of the galaxy. Within the anisotropic and axisymmetric velocity model, the stronger bisymmetry is explained by a radial component that is larger than the azimuthal component of the ellipsoid of random motions, thus by gaseous orbits that are dominantly radial. The azimuthal anisotropy parameter is not strongly dependent on the choice of the vertical dispersion. The velocity anisotropy parameter of the molecular clouds is observed highly scattered. Perturbations such as HI spiral-like arms could be at the origin of the gas velocity anisotropy in M33. Further work is necessary to assess whether anisotropic velocity ellispsoids can also be invoked to explain the asymmetric gas random motions of other galaxies. △ Less

Submitted 12 June, 2020; originally announced June 2020.

Comments: 10 pages, 9 figures, published in Astronomy & Astrophysics

Radio emission during the formation of stellar clusters in M33

Authors: Edvige Corbelli, Jonathan Braine, Fatemeh S. Tabatabaei

Abstract: We investigate thermal and non-thermal radio continuum associated with the early formation and evolution of Young Stellar Clusters (YSCs) selected by their MIR emission in M33. For the first time in an external galaxy it has been possible to identify radio counterparts to more than 300 star forming regions. We proof the nature of candidate YSCs fully embedded in molecular clouds, by recovering the… ▽ More We investigate thermal and non-thermal radio continuum associated with the early formation and evolution of Young Stellar Clusters (YSCs) selected by their MIR emission in M33. For the first time in an external galaxy it has been possible to identify radio counterparts to more than 300 star forming regions. We proof the nature of candidate YSCs fully embedded in molecular clouds, by recovering their associated faint radio continuum luminosities. Using the Halpha line to identify free-free radio emission at 5 GHz in the more evolved, partially exposed YSCs, we retrieve information on the relevance of magnetic fields and cosmic rays across the M33 disk at 25 pc spatial scales. A cross-correlation of MIR and radio continuum luminosities is established from bright to very faint YSCs, with MIR-to-radio emission ratio showing a gradual decline towards the outer disk, while the magnetic field is pervasive at all radii. We establish and discuss the tight relation between radio continuum and other star formation indicators, such as Halpha. This relation holds for individual YSCs over four orders of magnitude as well as for molecular clouds hosting YSCs. On average about half of radio emission at 5 GHz in YSCs is non-thermal. For exposed but compact YSCs the non-thermal radio fraction increases with source brightness, while for large HII regions the fraction is lower and shows no clear trend. This has been found for YSCs with and without identified SNRs and underlines the possible role of massive stars in triggering particle acceleration through winds and shocks: these particles diffuse throughout the native molecular cloud prior to cloud dispersal. △ Less

Submitted 12 May, 2020; originally announced May 2020.

Comments: 14 pages, 10 figures, accepted for publication in A&A

Journal ref: A&A 639, A27 (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)

A Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE).VI. Environmental quenching on HII region scales

Authors: A. Boselli, M. Fossati, A. Longobardi, S. Boissier, M. Boquien, J. Braine, P. Cote, J. C. Cuillandre, B. Epinat, L. Ferrarese, G. Gavazzi, S. Gwyn, G. Hensler, H. Plana, Y. Roehlly, C. Schimd, M. Sun, G. Trinchieri

Abstract: The Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE) is a blind narrow-band Halpha+[NII] imaging survey of the Virgo cluster carried out with MegaCam at the Canada-French-Hawaii telescope (CFHT). We use a new set of data extracted from VESTIGE to study the impact of the hostile cluster environment on the star formation process down to the scale of HII regions (~ 50 pc). HII region… ▽ More The Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE) is a blind narrow-band Halpha+[NII] imaging survey of the Virgo cluster carried out with MegaCam at the Canada-French-Hawaii telescope (CFHT). We use a new set of data extracted from VESTIGE to study the impact of the hostile cluster environment on the star formation process down to the scale of HII regions (~ 50 pc). HII regions are identified and their parameters measured using the HIIphot code on a sample of 114 late-type galaxies spanning a wide range in morphological type (Sa-Sd, Im, BCD), stellar mass (10^6.5 <= M_star <= 10^11 Mo), and star formation activity (10^-3 <= SFR <= 10 Mo yr^-1). Owing to the exquisite average resolution of the VESTIGE data (0.65 arcsec), we detect 11302 HII regions with an Halpha luminosity L(Halpha) >= 10^37 erg s^-1. We show that the typical number of HII regions in gas-stripped objects is significantly lower than in healthy late-types of similar stellar mass. We also show that in these gas-stripped galaxies the number of HII regions significantly drops outside the effective radius, suggesting that the quenching process occurs outside-in, in agreement with other multifrequency observations. These new results consistently confirm that the main mechanism responsible for the decrease of the star formation activity observed in cluster galaxies is ram pressure, allowing us to discard other milder processes such as starvation or strangulation unable to reproduce the observed radially truncated profiles. △ Less

Submitted 17 January, 2020; originally announced January 2020.

Comments: Accepted for publication on A&A Letter

Journal ref: A&A 634, L1 (2020)

Rotation of molecular clouds in M~51

Authors: J. Braine, A. Hughes, E. Rosolowsky, P. Gratier, D. Colombo, S. Meidt, E. Schinnerer

Abstract: The grand-design spiral galaxy M~51 was observed at 40pc resolution in CO(1--0) by the PAWS project. A large number of molecular clouds were identified and we search for velocity gradients in two high signal-to-noise subsamples, containing 682 and 376 clouds. The velocity gradients are found to be systematically prograde oriented, as was previously found for the rather flocculent spiral M~33. This… ▽ More The grand-design spiral galaxy M~51 was observed at 40pc resolution in CO(1--0) by the PAWS project. A large number of molecular clouds were identified and we search for velocity gradients in two high signal-to-noise subsamples, containing 682 and 376 clouds. The velocity gradients are found to be systematically prograde oriented, as was previously found for the rather flocculent spiral M~33. This strongly supports the idea that the velocity gradients reflect cloud rotation, rather than more random dynamical forces, such as turbulence. Not only are the gradients prograde, but their $\frac{\partial v}{\partial x}$ and $\frac{\partial v}{\partial y}$ coefficients follow galactic shear in sign, although with a lower amplitude. No link is found between the orientation of the gradient and the orientation of the cloud. The values of the cloud angular momenta appear to be an extension of the values noted for galactic clouds despite the orders of magnitude difference in cloud mass. Roughly 30\% of the clouds show retrograde velocity gradients. For a strictly rising rotation curve, as in M~51, gravitational contraction would be expected to yield strictly prograde rotators within an axisymmetric potential. In M~51, the fraction of retrograde rotators is found to be higher in the spiral arms than in the disk as a whole. Along the leading edge of the spiral arms, a majority of the clouds are retrograde rotators. While this work should be continued on other nearby galaxies, the M~33 and M~51 studies have shown that clouds rotate and that they rotate mostly prograde, although the amplitudes are not such that rotational energy is a significant support mechanism against gravitation. In this work, we show that retrograde rotation is linked to the presence of a spiral gravitational potential. △ Less

Submitted 20 November, 2019; originally announced November 2019.

Comments: Accepted by Astronomy and Astrophysics

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

Physical conditions in the gas phases of the giant HII region LMC-N11: II. Origin of [CII] and fraction of CO-dark gas

Authors: V. Lebouteiller, D. Cormier, S. C. Madden, M. Galametz, S. Hony, Frédéric Galliano, M. Chevance, Min-Young Lee, J. Braine, F. L. Polles, M. Angel Requeña-Torres, R. Indebetouw, A. Hughes, N. Abel

Abstract: (abridged) The ambiguous origin of [CII] 158um in the interstellar medium complicates its use for diagnostics concerning the star-formation rate and physical conditions in photodissociation regions (PDRs). We observed the giant HII region N11 in the Large Magellanic Cloud with SOFIA/GREAT in order to investigate the origin of [CII] to obtain the total H2 gas content, the fraction of CO-dark H2 gas… ▽ More (abridged) The ambiguous origin of [CII] 158um in the interstellar medium complicates its use for diagnostics concerning the star-formation rate and physical conditions in photodissociation regions (PDRs). We observed the giant HII region N11 in the Large Magellanic Cloud with SOFIA/GREAT in order to investigate the origin of [CII] to obtain the total H2 gas content, the fraction of CO-dark H2 gas, and the influence of environmental effects such as stellar feedback. We present an innovative spectral decomposition method that allows statistical trends to be derived. The [CII] line is resolved in velocity and compared to HI and CO, using a Bayesian approach to decompose the profiles. A simple model accounting for collisions in the neutral atomic and molecular gas was used in order to derive the H2 column density traced by C+. The profile of [CII] most closely resembles that of CO, but the integrated [CII] line width lies between that of CO and that of HI. Using various methods, we find that [CII] mostly originates from the neutral gas. We show that [CII] mostly traces the CO-dark H2 gas but there is evidence of a weak contribution from neutral atomic gas preferentially in the faintest components. Most of the molecular gas is CO-dark. The fraction of CO-dark H2 gas decreases with increasing CO column density, with a slope that seems to depend on the impinging radiation field from nearby massive stars. Finally we extend previous measurements of the photoelectric-effect heating efficiency, which we find is constant across regions probed with Herschel, with [CII] and [OI] being the main coolants in faint and diffuse, and bright and compact regions, respectively, and with PAH emission tracing the CO-dark H2 gas heating where [CII] and [OI] emit. Our study highlights the importance of velocity-resolved PDR diagnostics and higher spatial resolution for HI observations. △ Less

Submitted 8 November, 2019; originally announced November 2019.

Comments: Accepted for publication in A&A

Journal ref: A&A 632, A106 (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

Warm H$_2$ as a probe of massive accretion and feedback through shocks and turbulence across cosmic time

Authors: Philip Appleton, Lee Armus, Francois Boulanger, Charles M. Bradford, Jonathan Braine, Volker Bromm, Peter Capak, Michelle Cluver, Asantha Cooray, Tanio Diaz-Santos, Eiichi Egami, Bjorn Emonts, Pierre Guillard, George Helou, Lauranne Lanz, Susanne Madden, Anne Medling, Ewan O'Sullivan, Patrick Ogle, Alexandra Pope, Guillaume Pineau des Forêts, J. Michael Shull, John-David Smith, Aditya Togi, C. Kevin Xu

Abstract: Galaxy formation depends on a complex interplay between gravitational collapse, gas accretion, merging, and feedback processes. Yet, after many decades of investigation, these concepts are poorly understood. This paper presents the argument that warm H$_2$ can be used as a tool to unlock some of these mysteries. Turbulence, shocks and outflows, driven by star formation, AGN activity or inflows, ma… ▽ More Galaxy formation depends on a complex interplay between gravitational collapse, gas accretion, merging, and feedback processes. Yet, after many decades of investigation, these concepts are poorly understood. This paper presents the argument that warm H$_2$ can be used as a tool to unlock some of these mysteries. Turbulence, shocks and outflows, driven by star formation, AGN activity or inflows, may prevent the rapid buildup of star formation in galaxies. Central to our understanding of how gas is converted into stars is the process by which gas can dissipate its mechanical energy through turbulence and shocks in order to cool. H$_2$ lines provide direct quantitative measurements of kinetic energy dissipation in molecular gas in galaxies throughout the Universe. Based on the detection of very powerful H$_2$ lines from z = 2 galaxies and proto-clusters at the detection limits of {\it Spitzer}, we are confident that future far-IR and UV H$_2$ observations will provide a wealth of new information and insight into galaxy evolution to high-z. Finally, at the very earliest epoch of star and galaxy formation, warm H$_2$ may also provide a unique glimpse of molecular gas collapse at 7 $<$ z $<$ 12 in massive dark matter (DM) halos on their way to forming the very first galaxies. Such measurements are beyond the reach of existing and planned observatories. △ Less

Submitted 15 March, 2019; originally announced March 2019.

Comments: Submitted as a science White Paper to the Astronomy and Astrophysics Astro 2020 Decadal Survey call issued by the National Academies of Sciences, Engineering and Medicine (March 11 2019)

Comparing the Properties of GMCs in M33 from Simulations and Observations

Authors: C. L. Dobbs, E. Rosolowsky, A. R. Pettitt, J. Braine, E. Corbelli, J. Sun

Abstract: We compare the properties of clouds in simulated M33 galaxies to those observed in the real M33. We apply a friends of friends algorithm and CPROPS to identify clouds, as well as a pixel by pixel analysis. We obtain very good agreement between the number of clouds, and maximum mass of clouds. Both are lower than occurs for a Milky Way-type galaxy and thus are a function of the surface density, siz… ▽ More We compare the properties of clouds in simulated M33 galaxies to those observed in the real M33. We apply a friends of friends algorithm and CPROPS to identify clouds, as well as a pixel by pixel analysis. We obtain very good agreement between the number of clouds, and maximum mass of clouds. Both are lower than occurs for a Milky Way-type galaxy and thus are a function of the surface density, size and galactic potential of M33. We reproduce the observed dependence of molecular cloud properties on radius in the simulations, and find this is due to the variation in gas surface density with radius. The cloud spectra also show good agreement between the simulations and observations, but the exact slope and shape of the spectra depends on the algorithm used to find clouds, and the range of cloud masses included when fitting the slope. Properties such as cloud angular momentum, velocity dispersions and virial relation are also in good agreement between the simulations and observations, but do not necessarily distinguish between simulations of M33 and other galaxy simulations. Our results are not strongly dependent on the level of feedback used here (10 and 20%) although they suggest that 15% feedback efficiency may be optimal. Overall our results suggest that the molecular cloud properties are primarily dependent on the gas and mass surface density, and less dependent on the localised physics such as the details of stellar feedback, or the numerical code used. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: 15 pages, 12 figures, accepted by MNRAS

Relationship between the Line Width of the Atomic and Molecular ISM in M33

Authors: Eric W. Koch, Erik W. Rosolowsky, Andreas Schruba, Adam Leroy, Amanda Kepley, Jonathan Braine, Julianne Dalcanton, Megan C. Johnson

Abstract: We investigate how the spectral properties of atomic (HI) and molecular (H$_2$) gas, traced by CO(2-1), are related in M33 on $80$ pc scales. We find the HI and CO(2-1) velocity at peak intensity to be highly correlated, consistent with previous studies. By stacking spectra aligned to the velocity of HI peak intensity, we find that the CO line width ($σ_{\rm HWHM}=4.6\pm0.9$ km s$^{-1}$;… ▽ More We investigate how the spectral properties of atomic (HI) and molecular (H$_2$) gas, traced by CO(2-1), are related in M33 on $80$ pc scales. We find the HI and CO(2-1) velocity at peak intensity to be highly correlated, consistent with previous studies. By stacking spectra aligned to the velocity of HI peak intensity, we find that the CO line width ($σ_{\rm HWHM}=4.6\pm0.9$ km s$^{-1}$; $σ_{\rm HWHM}$ is the effective Gaussian width) is consistently smaller than the HI line width ($σ_{\rm HWHM}=6.6\pm0.1$ km s$^{-1}$), with a ratio of ${\sim}0.7$, in agreement with Druard et al. (2014). The ratio of the line widths remains less than unity when the data are smoothed to a coarser spatial resolution. In other nearby galaxies, this line width ratio is close to unity which has been used as evidence for a thick, diffuse molecular disk that is distinct from the thin molecular disk dominated by molecular clouds. The smaller line width ratio found here suggests that M33 has a marginal thick molecular disk. From modelling individual lines-of-sight, we recover a strong correlation between HI and CO line widths when only the HI located closest to the CO component is considered. The median line width ratio of the line-of-sight line widths is $0.56\pm0.01$. There is substantial scatter in the HI--CO(2-1) line width relation, larger than the uncertainties, that results from regional variations on $<500$ pc scales, and there is no significant trend in the line widths, or their ratios, with galactocentric radius. These regional line width variations may be a useful probe of changes in the local cloud environment or the evolutionary state of molecular clouds. △ Less

Submitted 13 February, 2019; originally announced February 2019.

Comments: Accepted for publication in MNRAS; 20 pages, 15 figures

Rise and fall of molecular clouds across the M33 disk

Authors: Edvige Corbelli, Jonathan Braine, Carlo Giovanardi

Abstract: We carried out deep searches for CO line emission in the outer disk of M33, at R>7 kpc, and examined the dynamical conditions that can explain variations in the mass distribution of the molecular cloud throughout the disk of M33. We used the IRAM-30~m telescope to search for CO lines in the outer disk toward 12 faint mid-infrared (MIR) selected sources and in an area of the southern outer disk hos… ▽ More We carried out deep searches for CO line emission in the outer disk of M33, at R>7 kpc, and examined the dynamical conditions that can explain variations in the mass distribution of the molecular cloud throughout the disk of M33. We used the IRAM-30~m telescope to search for CO lines in the outer disk toward 12 faint mid-infrared (MIR) selected sources and in an area of the southern outer disk hosting MA1, a bright HII region. We detect narrow CO lines at the location of two MIR sources at galactocentric distances of about 8 kpc that are associated with low-mass young stellar clusters, and at four locations in the proximity of MA1. The paucity of CO lines at the location of weak MIR-selected sources probably arises because most of them are not star-forming sites in M33, but background sources. Although very uncertain, the total molecular mass of the detected clouds around MA1 is lower than expected given the stellar mass of the cluster, because dispersal of the molecular gas is taking place as the HII region expands. The mean mass of the giant molecular clouds (GMCs) in M33 decreases radially by a factor 2 from the center out to 4 kpc, then it stays constant until it drops at R>7 kpc. We suggest that GMCs become more massive toward the center because of the fast rotation of the disk, which drives mass growth by coalescence of smaller condensations as they cross the arms. The analysis of both HI and CO spectral data gives the consistent result that corotation of the two main arms in this galaxy is at a radius of 4.7+-0.3 kpc, and spiral shock waves become subsonic beyond 3.9 kpc. Perturbations are quenched beyond 6.5 kpc, where CO lines have been detected only around sporadic condensations associated with UV and MIR emission. △ Less

Submitted 17 January, 2019; originally announced January 2019.

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

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

A Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE).V. Properties of the ionised gas filament of M87

Authors: A. Boselli, M. Fossati, A. Longobardi, G. Consolandi, P. Amram, M. Sun, P. Andreani, M. Boquien, J. Braine, F. Combes, P. Cote, J. C. Cuillandre, P. A. Duc, E. Emsellem, L. Ferrarese, G. Gavazzi, S. Gwyn, G. Hensler, E. W. Peng, H. Plana, J. Roediger, R. Sanchez-Janssen, M. Sarzi, P. Serra, G. Trinchieri

Abstract: We have observed the giant elliptical galaxy M87 during the Virgo Environmental Survey Tracing Galaxy Evolution (VESTIGE), a blind narrow-band Halpha+[NII] imaging survey of the Virgo cluster carried out with MegaCam at the Canada French Hawaii Telescope (CFHT). The deep narrow-band image confirmed the presence of a filament of ionised gas extending up to ~ 3 kpc in the north-western direction and… ▽ More We have observed the giant elliptical galaxy M87 during the Virgo Environmental Survey Tracing Galaxy Evolution (VESTIGE), a blind narrow-band Halpha+[NII] imaging survey of the Virgo cluster carried out with MegaCam at the Canada French Hawaii Telescope (CFHT). The deep narrow-band image confirmed the presence of a filament of ionised gas extending up to ~ 3 kpc in the north-western direction and ~ 8 kpc to the south-east, with a couple of plumes of ionised gas, the weakest of which, at ~ 18 kpc from the nucleus, was previously unknown. The analysis of deep optical images taken from the NGVS survey confirms that this gas filament is associated with dust seen in absorption which is now detected up to ~ 2.4 kpc from the nucleus. We have also analysed the physical and kinematical properties of the ionised gas filament using deep IFU MUSE data covering the central 4.8 x 4.8 kpc^2 of the galaxy. The spectroscopic data confirms a perturbed kinematics of the ionised gas, with differences in velocity of ~ 700-800 km s^-1 on scales of < 1 kpc. The analysis of 2D diagnostic diagrams and the observed relationship between the shock-sensitive [OI]/Halpha line ratio and the velocity dispersion of the gas suggest that the gas is shock-ionised. △ Less

Submitted 4 December, 2018; v1 submitted 23 October, 2018; originally announced October 2018.

Comments: Accepted for publication to Astronomy & Astrophysics

Journal ref: A&A 623, A52 (2019)

The Probability Distribution Function of Gas Surface Density in M33

Authors: Edvige Corbelli, Bruce G. Elmegreen, Jonathan Braine, David Thilker

Abstract: The probability distribution functions (PDFs) for atomic, molecular, and total gas surface densities of M33 are determined at a resolution of about 50~pc over regions that share coherent morphological properties to unveil fingerprints of self-gravity across the star-forming disk. Most of the total gas PDFs from the central region to the edge of the star-forming disk are well-fitted by log-normal f… ▽ More The probability distribution functions (PDFs) for atomic, molecular, and total gas surface densities of M33 are determined at a resolution of about 50~pc over regions that share coherent morphological properties to unveil fingerprints of self-gravity across the star-forming disk. Most of the total gas PDFs from the central region to the edge of the star-forming disk are well-fitted by log-normal functions whose width decreases radially outwards. Because the HI velocity dispersion is approximately constant across the disk, the decrease of the PDF width is consistent with a lower Mach number for the turbulent ISM at large galactocentric radii where a higher fraction of HI is in the warm phase. The atomic gas is found mostly at face-on column densities below N$_{H}^{lim}$=2.5 10$^{21}$~cm$^{-2}$, with small radial variations of N$_{H}^{lim}$. The molecular gas PDFs do not show strong deviations from log-normal functions in the central region where molecular fractions are high. Here the high pressure and rate of star formation shapes the PDF as a log-normal function dispersing self-gravitating complexes with intense feedback at all column densities that are spatially resolved. Power law PDFs for the molecules are found near and above N$_H^{lim}$, in the well defined southern spiral arm and in a continuous dense filament extending at larger galactocentric radii; this is evident in cloud samples at different evolutionary stages along the star formation cycle. In the filament nearly half of the molecular gas departs from a log-normal PDF and power laws are also observed in pre-star forming molecular complexes. The slope of the power law is between -1 and -2. This slope, combined with maps showing where the different parts of the power law PDFs come from, suggest a power-law stratification of density within molecular cloud complexes, which is consistent with the dominance of self-gravity. △ Less

Submitted 30 June, 2018; originally announced July 2018.

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

Journal ref: A&A 617, A125 (2018)

Kinematics of the Atomic ISM in M33 on 80 pc scales

Authors: Eric W. Koch, Erik W. Rosolowsky, Felix J. Lockman, Amanda A. Kepley, Adam Leroy, Andreas Schruba, Jonathan Braine, Julianne Dalcanton, Megan C. Johnson, Snežana Stanimirović

Abstract: We present new L-band ($1\mbox{--}2$ GHz) observations of the nearby spiral galaxy M33 with 80 pc resolution obtained with the Karl G. Jansky Very Large Array. The HI observations, combined with HI measurements from the Green Bank Telescope, improve the spectral resolution and sensitivity ($2.8$ K rms noise in a $0.2$ km s$^{-1}$ channel) compared to previous observations. We find individual profi… ▽ More We present new L-band ($1\mbox{--}2$ GHz) observations of the nearby spiral galaxy M33 with 80 pc resolution obtained with the Karl G. Jansky Very Large Array. The HI observations, combined with HI measurements from the Green Bank Telescope, improve the spectral resolution and sensitivity ($2.8$ K rms noise in a $0.2$ km s$^{-1}$ channel) compared to previous observations. We find individual profiles are usually non-Gaussian, harbouring line wings, multiple components, and asymmetries. Given this spectral complexity, we quantify the motions in the atomic ISM through moment analysis of the spectra and fits to aligned, stacked profiles. The measured value of the HI line width depends strongly on the method used, with the velocity stacked profiles aligned to the peak velocity giving the minimum value of $σ= 7$ km s$^{-1}$ and all other methods giving higher values ($σ\sim10$ km s$^{-1}$). All measurements of the line width show a shallow radial trend, with $σ$ decreasing by $\sim2$ km s$^{-1}$ from $R_{\rm gal}=0$ to $R_{\rm gal}=8$ kpc. We consider a number of energy sources that might maintain the line width against turbulent dissipation, but no single source is adequate. We find excess emission relative to a Gaussian in the stacked profile line wings, ranging from 9% to 26% depending on how the spectra are aligned. By splitting the line wings into symmetric and asymmetric components, we find that the lagging rotational disk accounts for one-third of the line wing flux. We also find emission far from the rotation axis of the galaxy in multiple discrete HI clouds, including a filament with a projected length of $\sim8$ kpc. △ Less

Submitted 8 August, 2018; v1 submitted 28 June, 2018; originally announced June 2018.

Comments: MNRAS, in press. 30 pages, 24 figures

A Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE).I. Introduction to the Survey

Authors: A. Boselli, M. Fossati, L. Ferrarese, S. Boissier, G. Consolandi, A. Longobardi, P. Amram, M. Balogh, P. Barmby, M. Boquien, F. Boulanger, J. Braine, V. Buat, D. Burgarella, F. Combes, T. Contini, L. Cortese, P. Cote, S. Cote, J. C. Cuillandre, L. Drissen, B. Epinat, M. Fumagalli, S. Gallagher, G. Gavazzi , et al. (34 additional authors not shown)

Abstract: The Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE) is a blind narrow-band Halpha+[NII] imaging survey carried out with MegaCam at the Canada-France-Hawaii Telescope. The survey covers the whole Virgo cluster region from its core to one virial radius (104 deg^2). The sensitivity of the survey is of f(Halpha) ~ 4 x 10^-17 erg sec-1 cm^-2 (5 sigma detection limit) for point sources… ▽ More The Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE) is a blind narrow-band Halpha+[NII] imaging survey carried out with MegaCam at the Canada-France-Hawaii Telescope. The survey covers the whole Virgo cluster region from its core to one virial radius (104 deg^2). The sensitivity of the survey is of f(Halpha) ~ 4 x 10^-17 erg sec-1 cm^-2 (5 sigma detection limit) for point sources and Sigma (Halpha) ~ 2 x 10^-18 erg sec^-1 cm^-2 arcsec^-2 (1 sigma detection limit at 3 arcsec resolution) for extended sources, making VESTIGE the deepest and largest blind narrow-band survey of a nearby cluster. This paper presents the survey in all its technical aspects, including the survey design, the observing strategy, the achieved sensitivity in both the narrow-band Halpha+[NII] and in the broad-band r filter used for the stellar continuum subtraction, the data reduction, calibration, and products, as well as its status after the first observing semester. We briefly describe the Halpha properties of galaxies located in a 4x1 deg^2 strip in the core of the cluster north of M87, where several extended tails of ionised gas are detected. This paper also lists the main scientific motivations of VESTIGE, which include the study of the effects of the environment on galaxy evolution, the fate of the stripped gas in cluster objects, the star formation process in nearby galaxies of different type and stellar mass, the determination of the Halpha luminosity function and of the Halpha scaling relations down to ~ 10^6 Mo stellar mass objects, and the reconstruction of the dynamical structure of the Virgo cluster. This unique set of data will also be used to study the HII luminosity function in hundreds of galaxies, the diffuse Halpha+[NII] emission of the Milky Way at high Galactic latitude, and the properties of emission line galaxies at high redshift. △ Less

Submitted 8 February, 2018; originally announced February 2018.

Comments: Accepted for publication on Astronomy & Astrophysics

Journal ref: A&A 614, A56 (2018)

Spatially resolving the dust properties and submillimetre excess in M 33

Authors: M. Relaño, I. De Looze, R. C. Kennicutt, U. Lisenfeld, A. Dariush, S. Verley, J. Braine, F. Tabatabaei, C. Kramer, M. Boquien, M. Xilouris, P. Gratier

Abstract: The relative abundance of the dust grain types in the interstellar medium (ISM) is directly linked to physical quantities that trace the evolution of galaxies. We study the dust properties of the whole disc of M33 at spatial scales of ~170 pc. This analysis allows us to infer how the relative dust grain abundance changes with the conditions of the ISM, study the existence of a submillimetre excess… ▽ More The relative abundance of the dust grain types in the interstellar medium (ISM) is directly linked to physical quantities that trace the evolution of galaxies. We study the dust properties of the whole disc of M33 at spatial scales of ~170 pc. This analysis allows us to infer how the relative dust grain abundance changes with the conditions of the ISM, study the existence of a submillimetre excess and look for trends of the gas-to-dust mass ratio (GDR) with other physical properties of the galaxy. For each pixel in the disc of M33 we fit the infrared SED using a physically motivated dust model that assumes an emissivity index beta close to 2. We derive the relative amount of the different dust grains in the model, the total dust mass, and the strength of the interstellar radiation field (ISRF) heating the dust at each spatial location. The relative abundance of very small grains tends to increase, and for big grains to decrease, at high values of Halpha luminosity. This shows that the dust grains are modified inside the star-forming regions, in agreement with a theoretical framework of dust evolution under different physical conditions. The radial dependence of the GDR is consistent with the shallow metallicity gradient observed in this galaxy. The strength of the ISRF derived in our model correlates with the star formation rate in the galaxy in a pixel by pixel basis. Although this is expected it is the first time that a correlation between both quantities is reported. We produce a map of submillimetre excess in the 500 microns SPIRE band for the disc of M33. The excess can be as high as 50% and increases at large galactocentric distances. We further study the relation of the excess with other physical properties of the galaxy and find that the excess is prominent in zones of diffuse ISM outside the main star-forming regions, where the molecular gas and dust surface density are low. △ Less

Submitted 15 January, 2018; originally announced January 2018.

Comments: Accepted to Astronomy and Astrophysics, 19 pages, 18 figures plus appendix

Journal ref: A&A 613, A43 (2018)

Properties and rotation of molecular clouds in M 33

Authors: Jonathan Braine, Erik Rosolowsky, Pierre Gratier, Edvige Corbelli, Karl Schuster

Abstract: The sample of 566 molecular clouds identified in the CO(2--1) IRAM survey covering the disk of M~33 is explored in detail.The clouds were found using CPROPS and were subsequently catalogued in terms of their star-forming properties as non-star-forming (A), with embedded star formation (B), or with exposed star formation C.We find that the size-linewidth relation among the M~33 clouds is quite weak… ▽ More The sample of 566 molecular clouds identified in the CO(2--1) IRAM survey covering the disk of M~33 is explored in detail.The clouds were found using CPROPS and were subsequently catalogued in terms of their star-forming properties as non-star-forming (A), with embedded star formation (B), or with exposed star formation C.We find that the size-linewidth relation among the M~33 clouds is quite weak but, when comparing with clouds in other nearby galaxies, the linewidth scales with average metallicity.The linewidth and particularly the line brightness decrease with galactocentric distance.The large number of clouds makes it possible to calculate well-sampled cloud mass spectra and mass spectra of subsamples.As noted earlier, but considerably better defined here, the mass spectrum steepens (i.e. higher fraction of small clouds) with galactocentric distance.A new finding is that the mass spectrum of A clouds is much steeper than that of the star-forming clouds.Further dividing the sample, this difference is strong at both large and small galactocentric distances and the A vs C difference is a stronger effect than the inner/outer disk difference in mass spectra.Velocity gradients are identified in the clouds using standard techniques.The gradients are weak and are dominated by prograde rotation; the effect is stronger for the high signal-to-noise clouds.A discussion of the uncertainties is presented.The angular momenta are low but compatible with at least some simulations.The cloud and galactic gradients are similar; the cloud rotation periods are much longer than cloud lifetimes and comparable to the galactic rotation period.The rotational kinetic energy is 1-2\% of the gravitational potential energy and the cloud edge velocity is well below the escape velocity, such that cloud-scale rotation probably has little influence on the evolution of molecular clouds. △ Less

Submitted 12 January, 2018; originally announced January 2018.

Comments: Accepted for publication in Astronomy and Astrophysics, 11 pages

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

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

Galaxy evolution studies with the SPace IR telescope for Cosmology and Astrophysics (SPICA): the power of IR spectroscopy

Authors: L. Spinoglio, A. Alonso-Herrero, L. Armus, M. Baes, J. Bernard-Salas, S. Bianchi, M. Bocchio, A. Bolatto, C. M. Bradford, J. Braine, F. J. Carrera, L. Ciesla, D. L. Clements, H. Dannerbauer, Y. Doi, A. Efstathiou, E. Egami, J. A. Fernandez-Ontiveros, A. Ferrara, J. Fischer, A. Franceschini, S. Gallerani, M. Giard, E. Gonzalez-Alfonso, C. Gruppioni , et al. (43 additional authors not shown)

Abstract: IR spectroscopy in the range 12-230 micron with the SPace IR telescope for Cosmology and Astrophysics (SPICA) will reveal the physical processes that govern the formation and evolution of galaxies and black holes through cosmic time, bridging the gap between the James Webb Space Telescope (JWST) and the new generation of Extremely Large Telescopes (ELTs) at shorter wavelengths and the Atacama Larg… ▽ More IR spectroscopy in the range 12-230 micron with the SPace IR telescope for Cosmology and Astrophysics (SPICA) will reveal the physical processes that govern the formation and evolution of galaxies and black holes through cosmic time, bridging the gap between the James Webb Space Telescope (JWST) and the new generation of Extremely Large Telescopes (ELTs) at shorter wavelengths and the Atacama Large Millimeter Array (ALMA) at longer wavelengths. SPICA, with its 2.5-m telescope actively-cooled to below 8K, will obtain the first spectroscopic determination, in the mid-IR rest-frame, of both the star-formation rate and black hole accretion rate histories of galaxies, reaching lookback times of 12 Gyr, for large statistically significant samples. Densities, temperatures, radiation fields and gas-phase metallicities will be measured in dust-obscured galaxies and active galactic nuclei (AGN), sampling a large range in mass and luminosity, from faint local dwarf galaxies to luminous quasars in the distant Universe. AGN and starburst feedback and feeding mechanisms in distant galaxies will be uncovered through detailed measurements of molecular and atomic line profiles. SPICA's large-area deep spectrophotometric surveys will provide mid-IR spectra and continuum fluxes for unbiased samples of tens of thousands of galaxies, out to redshifts of z~6. Furthermore, SPICA spectroscopy will uncover the most luminous galaxies in the first few hundred million years of the Universe, through their characteristic dust and molecular hydrogen features. △ Less

Submitted 5 October, 2017; originally announced October 2017.

Comments: Paper accepted for publication on PASA on 21st September 2017, as part of the SPICA Special Issue

SPICA and the Chemical Evolution of Galaxies: The Rise of Metals and Dust

Authors: J. A. Fernández-Ontiveros, L. Armus, M. Baes, J. Bernard-Salas, A. D. Bolatto, J. Braine, L. Ciesla, I. De Looze, E. Egami, J. Fischer, M. Giard, E. González-Alfonso, G. L. Granato, C. Gruppioni, M. Imanishi, D. Ishihara, H. Kaneda, S. Madden, M. Malkan, H. Matsuhara, M. Matsuura, T. Nagao, F. Najarro, T. Nakagawa, T. Onaka , et al. (11 additional authors not shown)

Abstract: The physical processes driving the chemical evolution of galaxies in the last $\sim 11\, \rm{Gyr}$ cannot be understood without directly probing the dust-obscured phase of star-forming galaxies and active galactic nuclei. This phase, hidden to optical tracers, represents the bulk of star formation and black hole accretion activity in galaxies at $1 < z < 3$. Spectroscopic observations with a cryog… ▽ More The physical processes driving the chemical evolution of galaxies in the last $\sim 11\, \rm{Gyr}$ cannot be understood without directly probing the dust-obscured phase of star-forming galaxies and active galactic nuclei. This phase, hidden to optical tracers, represents the bulk of star formation and black hole accretion activity in galaxies at $1 < z < 3$. Spectroscopic observations with a cryogenic infrared (IR) observatory like SPICA will be sensitive enough to peer through the dust-obscured regions of galaxies and access the rest-frame mid- to far-IR range in galaxies at high-$z$. This wavelength range contains a unique suite of spectral lines and dust features that serve as proxies for the abundances of heavy elements and the dust composition, providing tracers with a feeble response to both extinction and temperature. In this work, we investigate how SPICA observations could be exploited to understand key aspects in the chemical evolution of galaxies: the assembly of nearby galaxies based on the spatial distribution of heavy element abundances, the global content of metals in galaxies reaching the knee of the luminosity function up to $z \sim 3$, and the dust composition of galaxies at high-$z$. Possible synergies with facilities available in the late 2020s are also discussed. △ Less

Submitted 5 October, 2017; originally announced October 2017.

Comments: In press. This paper belongs to the "SPICA Special Issue" on PASA

Data processing pipeline for Herschel HIFI

Authors: R. F. Shipman, S. F. Beaulieu, D. Teyssier, P. Morris, M. Rengel, C. McCoey, K. Edwards, D. Kester, A. Lorenzani, O. Coeur-Joly, M. Melchior, J. Xie, E. Sanchez, P. Zaal, I. Avruch, C. Borys, J. Braine, C. Comito1, B. Delforge, F. Herpin1, A. Hoac, W. Kwon, S. D. Lord, A. Marston, M. Mueller , et al. (4 additional authors not shown)

Abstract: {Context}. The HIFI instrument on the Herschel Space Observatory performed over 9100 astronomical observations, almost 900 of which were calibration observations in the course of the nearly four-year Herschel mission. The data from each observation had to be converted from raw telemetry into calibrated products and were included in the Herschel Science Archive. {Aims}. The HIFI pipeline was design… ▽ More {Context}. The HIFI instrument on the Herschel Space Observatory performed over 9100 astronomical observations, almost 900 of which were calibration observations in the course of the nearly four-year Herschel mission. The data from each observation had to be converted from raw telemetry into calibrated products and were included in the Herschel Science Archive. {Aims}. The HIFI pipeline was designed to provide robust conversion from raw telemetry into calibrated data throughout all phases of the HIFI missions. Pre-launch laboratory testing was supported as were routine mission operations. {Methods}. A modular software design allowed components to be easily added, removed, amended and/or extended as the understanding of the HIFI data developed during and after mission operations. {Results}. The HIFI pipeline processed data from all HIFI observing modes within the Herschel automated processing environment as well as within an interactive environment. The same software can be used by the general astronomical community to reprocess any standard HIFI observation. The pipeline also recorded the consistency of processing results and provided automated quality reports. Many pipeline modules were in use since the HIFI pre-launch instrument level testing. {Conclusions}. Processing in steps facilitated data analysis to discover and address instrument artefacts and uncertainties. The availability of the same pipeline components from pre-launch throughout the mission made for well-understood, tested, and stable processing. A smooth transition from one phase to the next significantly enhanced processing reliability and robustness. △ Less

Submitted 13 September, 2017; originally announced September 2017.

Comments: 22 pages, 12 figures

Journal ref: A&A 608, A49 (2017)

Testing the universality of the star formation efficiency in dense molecular gas

Authors: Y. Shimajiri, Ph. André, J. Braine, V. Könyves, N. Schneider, S. Bontemps, B. Ladjelate, A. Roy, Y. Gao, H. Chen

Abstract: Recent studies have suggested that star formation (SF) in dense gas may be governed by essentially the same law in Galactic clouds and external galaxies. This conclusion remains controversial, however, because different tracers have been used to probe the dense gas mass in Galactic and extragalactic studies. We conducted observations of Aquila, Oph., and Orion B in HCN(1-0), HCO$^+$(1-0), and thei… ▽ More Recent studies have suggested that star formation (SF) in dense gas may be governed by essentially the same law in Galactic clouds and external galaxies. This conclusion remains controversial, however, because different tracers have been used to probe the dense gas mass in Galactic and extragalactic studies. We conducted observations of Aquila, Oph., and Orion B in HCN(1-0), HCO$^+$(1-0), and their isotopomers to calibrate the HCN and HCO$^+$ lines used as dense gas tracers in extragalactic studies and to test the possible universality of the star formation efficiency in dense gas, SFE$_{dense}$. H$^{13}$CO$^+$ and H$^{13}$CN were observed to be good tracers of the filaments detected with Herschel. Comparing the luminosities $L_{HCN}$ and $L_{HCO^+}$ with the reference masses $M_{Herschel}^{Av>8}$, the empirical conversion factors $α_{Herschel-HCN}$ and $α_{Herschel-HCO^+}$ were found to be anti-correlated with the local FUV strength. In agreement with Pety et al. (2017), HCN and HCO$^+$ were also found to trace gas down to Av>2. As a result, published extragalactic HCN studies must be tracing all of the moderate density gas. Estimating the contribution of this moderate density gas from the typical column density PDFs in nearby clouds, we obtained the following $G_0$-dependent HCN conversion factor for external galaxies:$α_{Herschel-HCN}^{fit'}=64\times G_0^{-0.34}$. Re-estimating the dense gas masses in external galaxies with $α_{Herschel-HCN}^{fit'}$, we found that SFE$_{dense}$ is remarkably constant over 8 orders of magnitude in dense gas mass. Our results confirm that SFE$_{dense}$ of galaxies is quasi-universal on a wide range of scales from ~1-10pc to >10kpc. Based on the tight link between SF and filamentary structure found in Herschel studies of nearby clouds, we argue that SFE$_{dense}$ is primarily set by the microphysics of SF in filaments. △ Less

Submitted 29 April, 2017; originally announced May 2017.

Comments: A&A accepted for publication

Journal ref: A&A 604, A74 (2017)

Predicting HCN, HCO+, multi-transition CO, and dust emission of star-forming galaxies From local spiral and ultraluminous infrared galaxies to high-z star-forming and submillimeter galaxies

Authors: B. Vollmer, P. Gratier, J. Braine, C. Bot

Abstract: (abridged) In this work we have a closer look at the gas content or fraction and the associated star formation rate in main sequence and starburst galaxies at z=0 and z~1-2 by applying an analytical model of galactic clumpy gas disks to samples of local spiral galaxies, ULIRGs, submillimeter (smm), and high-z starforming galaxies. The model gas and dust temperatures are determined by the heating a… ▽ More (abridged) In this work we have a closer look at the gas content or fraction and the associated star formation rate in main sequence and starburst galaxies at z=0 and z~1-2 by applying an analytical model of galactic clumpy gas disks to samples of local spiral galaxies, ULIRGs, submillimeter (smm), and high-z starforming galaxies. The model gas and dust temperatures are determined by the heating and cooling equilibrium. Dense clouds are heated by turbulent mechanical and cosmic ray heating. The molecular abundances of individual gas clouds are determined by a detailed chemical network involving the cloud lifetime, density, and temperature. Molecular line emission is calculated with an escape probability formalism. The model calculates simultaneously the total gas mass, HI/H_2 mass, the gas velocity dispersion, IR luminosity, IR spectral energy distribution, CO spectral line energy distribution (SLED), HCN(1-0), and HCO+(1-0) emission of a galaxy given its size, integrated star formation rate, stellar mass radial profile, rotation curve, and Toomre Q parameter. The model reproduces the observed CO luminosities and SLEDs of all sample galaxies within the model uncertainties (~0.3 dex). Whereas the CO emission is robust against the variation of model parameters, the HCN and HCO+ emission is sensitive to the chemistry of the interstellar medium. The CO and HCN mass-to-light conversion factors including CO-dark H_2 are given and compared to the values found in the literature. Both, the HCN and HCO+ emission trace the dense molecular gas to a factor of ~2 for the local spiral galaxies, ULIRGs and smm-galaxies. About 80% of the molecular line emission of compact starburst galaxies originates in non-selfgravitating gas clouds. The integrated Kennicutt-Schmidt law has a slope of ~1 for the local spirals, ULIRGs, and smm-galaxies, whereas the slope is 1.7 for high-z starforming galaxies. △ Less

Submitted 25 January, 2017; v1 submitted 20 December, 2016; originally announced December 2016.

Comments: accepted for publication in A&A; corrected version The smm and high-z starforming galaxies were erroneously treated as if they were local galaxies. This lead to differences in the CO fluxes of about 10% and in the HCN and HCO+ fluxes of about 25-30%. Given that the model uncertainties are of the order of a factor of 2 (Sect. 4), this does not change the discussion nor the conclusions

Dense gas in low-metallicity galaxies

Authors: J. Braine, Y. Shimajiri, P. André, S. Bontemps, Yu Gao, Hao Chen, C. Kramer

Abstract: Stars form out of the densest parts of molecular clouds. Far-IR emission can be used to estimate the Star Formation Rate (SFR) and high dipole moment molecules, typically HCN, trace the dense gas. A strong correlation exists between HCN and Far-IR emission, with the ratio being nearly constant, over a large range of physical scales. A few recent observations have found HCN to be weak with respect… ▽ More Stars form out of the densest parts of molecular clouds. Far-IR emission can be used to estimate the Star Formation Rate (SFR) and high dipole moment molecules, typically HCN, trace the dense gas. A strong correlation exists between HCN and Far-IR emission, with the ratio being nearly constant, over a large range of physical scales. A few recent observations have found HCN to be weak with respect to the Far-IR and CO in subsolar metallicity (low-Z) objects. We present observations of the Local Group galaxies M33, IC10, and NGC6822 with the IRAM 30meter and NRO 45m telescopes, greatly improving the sample of low-Z galaxies observed. HCN, HCO$^+$, CS, C$_2$H, and HNC have been detected. Compared to solar metallicity galaxies, the Nitrogen-bearing species are weak (HCN, HNC) or not detected (CN, HNCO, N$_2$H$^+$) relative to Far-IR or CO emission. HCO$^+$ and C$_2$H emission is normal with respect to CO and Far-IR. While $^{13}$CO is the usual factor 10 weaker than $^{12}$CO, C$^{18}$O emission was not detected down to very low levels. Including earlier data, we find that the HCN/HCO$^+$ ratio varies with metallicity (O/H) and attribute this to the sharply decreasing Nitrogen abundance. The dense gas fraction, traced by the HCN/CO and HCO$^+$/CO ratios, follows the SFR but in the low-Z objects the HCO$^+$ is much easier to measure. Combined with larger and smaller scale measurements, the HCO$^+$ line appears to be an excellent tracer of dense gas and varies linearly with the SFR for both low and high metallicities. △ Less

Submitted 7 December, 2016; originally announced December 2016.

Comments: 9 pages, 13 figures, A&A in press

Journal ref: A&A 2016

Dense Gas in the Outer Spiral Arm of M51

Authors: Hao Chen, Jonathan Braine, Yu Gao, Jin Koda, Qiusheng Gu

Abstract: There is a linear relation between the mass of dense gas, traced by the HCN(1-0) luminosity, and the star formation rate (SFR), traced by the far-infrared luminosity. Recent observations of galactic disks have shown some systematic variations. In order to explore the SFR-dense gas link at high resolution ($\sim 4"$, $\sim 150$ pc) in the outer disk of an external galaxy, we have mapped a region ab… ▽ More There is a linear relation between the mass of dense gas, traced by the HCN(1-0) luminosity, and the star formation rate (SFR), traced by the far-infrared luminosity. Recent observations of galactic disks have shown some systematic variations. In order to explore the SFR-dense gas link at high resolution ($\sim 4"$, $\sim 150$ pc) in the outer disk of an external galaxy, we have mapped a region about 5 kpc from the center along the northern spiral arm of M51 in the HCN(1-0), HCO$^+$(1-0) and HNC(1-0) emission lines using the Northern Extended Millimeter Array (NOEMA) interferometer. The HCN and HCO$^+$ lines were detected in 6 giant molecular associations (GMAs) while HNC emission was only detected in the two brightest GMAs. One of the GMAs hosts a powerful HII region and HCN is stronger than HCO$^+$ there. Comparing with observations of GMAs in the disks of M31 and M33 at similar angular resolution ($\sim 100$ pc), we find that GMAs in the outer disk of M51 are brighter in both HCN and HCO$^+$ lines by a factor of 3 on average. However, the $I_{HCN}/I_{CO}$ and $I_{HCO^+}/I_{CO}$ ratios are similar to the ratios in nearby galactic disks and the Galactic plane. Using the Herschel 70 $μ$m data to trace the total IR luminosity at the resolution of the GMAs, we find that both the L$_{IR}$-L$_{HCN}$ and L$_{IR}$-L$_{HCO^+}$ relations in the outer disk GMAs are consistent with the proportionality between the L$_{IR}$ and the dense gas mass established globally in galaxies within the scatter. The IR/HCN and IR/HCO$^+$ ratios of the GMAs vary by a factor of 3, probably depending on whether massive stars are forming or not. △ Less

Submitted 1 December, 2016; originally announced December 2016.

Comments: 8 pages, 5 figures, accepted by ApJ