Showing 1–50 of 67 results for author: Relaño, M

CAVITY, Calar Alto Void Integral-field Treasury surveY and project extension

Authors: I. Pérez, S. Verley, L. Sánchez-Menguiano, T. Ruiz-Lara, R. García-Benito, S. Duarte Puertas, A. Jiménez, J. Domínguez-Gómez, D. Espada, R. F. Peletier, J. Román, M. I. Rodríguez, P. Sánchez Alarcón, M. Argudo-Fernández, G. Torres-Ríos, B. Bidaran, M. Alcázar-Laynez, R. van de Weygaert, S. F. Sánchez, U. Lisenfeld, A. Zurita, E. Florido, J. M. van der Hulst, G. Blázquez-Calero, P. Villalba-González , et al. (36 additional authors not shown)

Abstract: We have learnt in the last decades that the majority of galaxies belong to high density regions interconnected in a sponge-like fashion. This large-scale structure is characterised by clusters, filaments, walls, where most galaxies concentrate, but also under-dense regions, called voids. The void regions and the galaxies within represent an ideal place for the study of galaxy formation and evoluti… ▽ More We have learnt in the last decades that the majority of galaxies belong to high density regions interconnected in a sponge-like fashion. This large-scale structure is characterised by clusters, filaments, walls, where most galaxies concentrate, but also under-dense regions, called voids. The void regions and the galaxies within represent an ideal place for the study of galaxy formation and evolution as they are largely unaffected by the complex physical processes that transform galaxies in high-density environments. These void galaxies can hold the key as well to answer current challenges to the $Λ$CDM paradigm. The Calar Alto Void Integral-field Treasury surveY (CAVITY) is a Legacy project approved by the Calar Alto Observatory to obtain spatially resolved spectroscopic information of $\sim300$ void galaxies in the Local Universe (0.005 < z < 0.050) covering from -17.0 to -21.5 in $\rm r$ band absolute magnitude. It officially started in January 2021 and has been awarded 110 useful dark observing nights at the 3.5 m telescope using the PMAS spectrograph. Complementary follow-up projects including deep optical imaging, integrated, as well as resolved CO data, and integrated HI spectra, have joint the PMAS observations and naturally complete the scientific aim of characterising galaxies in cosmic voids. The extension data has been denominated CAVITY+. The data will be available to the whole community in different data releases, the first of which is planned for July 2024, and it will provide the community with PMAS data cubes for around 100 void galaxies through a user friendly, and well documented, database platform. We present here the survey, sample selection, data reduction, quality control schemes, science goals, and some examples of the scientific power of the CAVITY and CAVITY+ data. △ Less

Submitted 24 May, 2024; v1 submitted 7 May, 2024; originally announced May 2024.

Comments: 19 pages, 11 figures. Accepted for publication in A&A - Replacement after A&A minor language edition

JWST Observations of Starbursts: Cold Clouds and Plumes Launching in the M82 Outflow

Authors: Deanne B. Fisher, Alberto D. Bolatto, John Chisholm, Drummond Fielding, Rebecca C. Levy, Elizabeth Tarantino, Martha L. Boyer, Serena A. Cronin, Laura A. Lopez, J. D. Smith, Danielle A. Berg, Sebastian Lopez, Sylvain Veilleux, Paul P. van der Werf, Torsten Böker, Leindert A. Boogaard, Laura Lenkić, Simon C. O. Glover, Vicente Villanueva, Divakara Mayya, Thomas S. -Y. Lai, Daniel A. Dale, Kimberly L. Emig, Fabian Walter, Monica Relaño , et al. (6 additional authors not shown)

Abstract: In this paper we study the filamentary substructure of 3.3 $μ$m PAH emission from JWST/NIRCam observations in the base of the M82 star-burst driven wind. We identify plume-like substructure within the PAH emission with widths of $\sim$50 pc. Several of the plumes extend to the edge of the field-of-view, and thus are at least 200-300 pc in length. In this region of the outflow, the vast majority (… ▽ More In this paper we study the filamentary substructure of 3.3 $μ$m PAH emission from JWST/NIRCam observations in the base of the M82 star-burst driven wind. We identify plume-like substructure within the PAH emission with widths of $\sim$50 pc. Several of the plumes extend to the edge of the field-of-view, and thus are at least 200-300 pc in length. In this region of the outflow, the vast majority ($\sim$70\%) of PAH emission is associated with the plumes. We show that those structures contain smaller scale "clouds" with widths that are $\sim$5-15 pc, and they are morphologically similar to the results of "cloud-crushing" simulations. We estimate the cloud-crushing time-scales of $\sim$0.5-3 Myr, depending on assumptions. We show this time scale is consistent with a picture in which these observed PAH clouds survived break-out from the disk rather than being destroyed by the hot wind. The PAH emission in both the midplane and the outflow is shown to tightly correlate with that of Pa$α$ emission (from HST/NICMOS data), at the scale of both plumes and clouds, though the ratio of PAH-to-Pa$α$ increases at further distances from the midplane. Finally, we show that the outflow PAH emission is suppressed in regions of the M82 wind that are bright in X-ray emission. Overall, our results are broadly consistent with a picture in which cold gas in galactic outflows is launched via hierarchically structured plumes, and those small scale clouds are more likely to survive the wind environment when collected into the larger plume structure. △ Less

Submitted 6 May, 2024; originally announced May 2024.

Comments: Submitted to MNRAS

The ALMA-CRISTAL survey: Extended [CII] emission in an interacting galaxy system at z ~ 5.5

Authors: A. Posses, M. Aravena, J. González-López, N. M. Förster Schreiber, D. Liu, L. Lee, M. Solimano, T. Díaz-Santos, R. J. Assef, L. Barcos-Muñoz, S. Bovino, R. A. A. Bowler, G. Calistro Rivera, E. da Cunha, R. L. Davies, M. Killi, I. De Looze, A. Ferrara, D. B. Fisher, R. Herrera-Camus, R. Ikeda, T. Lambert, J. Li, D. Lutz, I. Mitsuhashi , et al. (9 additional authors not shown)

Abstract: The ALMA [CII] Resolved Ism in STar-forming gALaxies (CRISTAL) survey is a Cycle 8 ALMA Large Programme that studies the cold gas component of high-redshift galaxies. Its sub-arcsecond resolution observations are key to disentangling physical mechanisms that shape galaxies during cosmic dawn. In this paper, we explore the morphology and kinematics of the cold gas, star-forming, and stellar compone… ▽ More The ALMA [CII] Resolved Ism in STar-forming gALaxies (CRISTAL) survey is a Cycle 8 ALMA Large Programme that studies the cold gas component of high-redshift galaxies. Its sub-arcsecond resolution observations are key to disentangling physical mechanisms that shape galaxies during cosmic dawn. In this paper, we explore the morphology and kinematics of the cold gas, star-forming, and stellar components in the star-forming main-sequence galaxy CRISTAL-05/HZ3, at z = 5.54. Our analysis includes 0.3" spatial resolution (~2 kpc) ALMA observations of the [CII] line. While CRISTAL-05 was previously classified as a single source, our observations reveal that the system is a close interacting pair surrounded by an extended component of carbon-enriched gas. This is imprinted in the disturbed elongated [CII] morphology and the separation of the two components in the position-velocity diagram (~100 km/s). The central region is composed of two components, named C05-NW and C05-SE, with the former being the dominant one. A significant fraction of the [CII] arises beyond the close pair up to 10 kpc, while the regions forming new massive stars and the stellar component seem compact (r_[CII] ~ 4 r_UV), as traced by rest-frame UV and optical imaging obtained with the Hubble Space Telescope and the James Webb Space Telescope. Our kinematic model, using the DYSMALpy software, yields a minor contribution of dark matter of C05-NW within a radius of ~2x Reff. Finally, we explore the resolved [CII]/FIR ratios as a proxy for shock-heating produced by this merger. We argue that the extended [CII] emission is mainly caused by the merger, which could not be discerned with lower-resolution observations. Our work emphasizes the need for high-resolution observations to fully characterize the dynamic stages of infant galaxies and the physical mechanisms that drive the metal enrichment of the circumgalactic medium. △ Less

Submitted 5 March, 2024; originally announced March 2024.

Comments: Submitted to A&A - comments are welcome! - 19 pages, 13 figures

Observational signatures of the dust size evolution in isolated galaxy simulations

Authors: Kosei Matsumoto, Hiroyuki Hirashita, Kentaro Nagamine, Stefan van der Giessen, Leonard E. C. Romano, Monica Relaño, Ilse De Looze, Maarten Baes, Angelos Nersesian, Peter Camps, Kuan-chou Hou, Yuri Oku

Abstract: We aim to provide observational signatures of the dust size evolution in the ISM, in particular exploring indicators of polycyclic aromatic hydrocarbon (PAH) mass fraction ($q_{PAH}$) defined as the mass fraction of PAHs relative to total dust grains. Additionally, we validate our dust evolution model by comparing the observational signatures from our simulations to those from observations. We mod… ▽ More We aim to provide observational signatures of the dust size evolution in the ISM, in particular exploring indicators of polycyclic aromatic hydrocarbon (PAH) mass fraction ($q_{PAH}$) defined as the mass fraction of PAHs relative to total dust grains. Additionally, we validate our dust evolution model by comparing the observational signatures from our simulations to those from observations. We model the evolution of grain size distribution of Milky Way-like and NGC 628-like galaxies representing star-forming galaxies with a hydrodynamic simulation code, GADGET4-OSAKA, which considers dust production and interstellar processing. Furthermore, we perform post-processing dust radiative transfer with SKIRT based on the simulations to predict the observational properties. We find that the intensity ratio between 8 um and 24 um correlates with $q_{PAH}$ and can be used as an indicator of PAH mass fraction. However, this ratio is influenced by the radiation field. As another indicator, we suggest the 8 um-to-total infrared intensity ratio ($νI_ν(8 μm)/I$(TIR)), which tightly correlates with $q_{PAH}$. Furthermore, we explore the spatial evolution of $q_{PAH}$ in the simulated Milky Way-like galaxy using $νI_ν(8 μm)/I$(TIR). We find that the spatially resolved $q_{PAH}$ increases with metallicity at lower metallicity (Z<0.2 Zsun) due to the interplay between accretion and shattering while it decreases with metallicity at higher metallicity (Z>0.2 Zsun) due to coagulation. Finally, we compare the above indicators in the NGC 628-like simulation with those observed in NGC 628. Consequently, our simulation underestimates the PAH mass fraction throughout the entire galaxy. This is probably because PAH is too efficiently lost by coagulation in the interstellar medium in our model, which suggests that the inhibition of coagulation of the PAHs is key to enhancing PAH formation. △ Less

Submitted 4 February, 2024; originally announced February 2024.

Comments: 16 pages, 12 figures, submit to A&A

JWST Observations of Starbursts: Polycyclic Aromatic Hydrocarbon Emission at the Base of the M 82 Galactic Wind

Authors: Alberto D. Bolatto, Rebecca C. Levy, Elizabeth Tarantino, Martha L. Boyer, Deanne B. Fisher, Adam K. Leroy, Serena A. Cronin, Ralf S. Klessen, J. D. Smith, Dannielle A. Berg, Torsten Boeker, Leindert A. Boogaard, Eve C. Ostriker, Todd A. Thompson, Juergen Ott, Laura Lenkic, Laura A. Lopez, Daniel A. Dale, Sylvain Veilleux, Paul P. van der Werf, Simon C. O. Glover, Karin M. Sandstrom, Evan D. Skillman, John Chisholm, Vicente Villanueva , et al. (15 additional authors not shown)

Abstract: We present new observations of the central 1 kpc of the M 82 starburst obtained with the James Webb Space Telescope (JWST) near-infrared camera (NIRCam) instrument at a resolution ~0.05"-0.1" (~1-2 pc). The data comprises images in three mostly continuum filters (F140M, F250M, and F360M), and filters that contain [FeII] (F164N), H2 v=1-0 (F212N), and the 3.3 um PAH feature (F335M). We find promine… ▽ More We present new observations of the central 1 kpc of the M 82 starburst obtained with the James Webb Space Telescope (JWST) near-infrared camera (NIRCam) instrument at a resolution ~0.05"-0.1" (~1-2 pc). The data comprises images in three mostly continuum filters (F140M, F250M, and F360M), and filters that contain [FeII] (F164N), H2 v=1-0 (F212N), and the 3.3 um PAH feature (F335M). We find prominent plumes of PAH emission extending outward from the central starburst region, together with a network of complex filamentary substructure and edge-brightened bubble-like features. The structure of the PAH emission closely resembles that of the ionized gas, as revealed in Paschen alpha and free-free radio emission. We discuss the origin of the structure, and suggest the PAHs are embedded in a combination of neutral, molecular, and photoionized gas. △ Less

Submitted 21 April, 2024; v1 submitted 29 January, 2024; originally announced January 2024.

Comments: Submitted to The Astrophysical Journal

The ALMA-CRISTAL survey. Discovery of a 15 kpc-long gas plume in a $z=4.54$ Lyman-$α$ blob

Authors: M. Solimano, J. González-López, M. Aravena, R. Herrera-Camus, I. De Looze, N. M. Förster Schreiber, J. Spilker, K. Tadaki, R. J. Assef, L. Barcos-Muñoz, R. L. Davies, T. Díaz-Santos, A. Ferrara, D. B. Fisher, L. Guaita, R. Ikeda, E. J. Johnston, D. Lutz, I. Mitsuhashi, C. Moya-Sierralta, M. Relaño, T. Naab, A. C. Posses, K. Telikova, H. Übler , et al. (2 additional authors not shown)

Abstract: Massive star-forming galaxies in the high-redshift universe host large reservoirs of cold gas in their circumgalactic medium (CGM). Traditionally, these reservoirs have been linked to diffuse H I Lyman-$α$ (Ly$α)$ emission extending beyond $\approx 10$ kpc scales. In recent years, millimeter/submillimeter observations are starting to identify even colder gas in the CGM through molecular and/or ato… ▽ More Massive star-forming galaxies in the high-redshift universe host large reservoirs of cold gas in their circumgalactic medium (CGM). Traditionally, these reservoirs have been linked to diffuse H I Lyman-$α$ (Ly$α)$ emission extending beyond $\approx 10$ kpc scales. In recent years, millimeter/submillimeter observations are starting to identify even colder gas in the CGM through molecular and/or atomic tracers such as the [C II] $158\,μ$m transition. In this context, we study the well-known J1000+0234 system at $z=4.54$ that hosts a massive dusty star-forming galaxy (DSFG), a UV-bright companion, and a Ly$α$ blob. We combine new ALMA [C II] line observations taken by the CRISTAL survey with data from previous programs targeting the J1000+0234 system, and achieve a deep view into a DSFG and its rich environment at a 0.2" resolution. We identify an elongated [C II]-emitting structure with a projected size of 15 kpc stemming from the bright DSFG at the center of the field, with no clear counterpart at any other wavelength. The plume is oriented $\approx 40^{\circ}$ away from the minor axis of the DSFG, and shows significant spatial variation of its spectral parameters. In particular, the [C II] emission shifts from 180 km/s to 400 km/s between the bottom and top of the plume, relative to the DSFG's systemic velocity. At the same time, the line width starts at 400-600 km/s but narrows down to 190 km/s at top end of the plume. We discuss four possible scenarios to interpret the [C II] plume: a conical outflow, a cold accretion stream, ram pressure stripping, and gravitational interactions. While we cannot strongly rule out any of these with the available data, we disfavor the ram pressure stripping scenario due to the requirement of special hydrodynamic conditions. △ Less

Submitted 9 January, 2024; originally announced January 2024.

Comments: 17 pages (14 main text, 2 for references and 1 appendix page), 7 figures and 4 tables. Submitted to A&A

The ALMA-CRISTAL survey: Widespread dust-obscured star formation in typical star-forming galaxies at z=4-6

Authors: Ikki Mitsuhashi, Ken-ichi Tadaki, Ryota Ikeda, Rodrigo Herrera-Camus, Manuel Aravena, Ilse De Looze, Natascha M. Förster Schreiber, Jorge González-López, Justin Spilker, Roberto J. Assef, Rychard Bouwens, Loreto Barcos-Munoz, Jack Birkin, Rebecca A. A. Bowler, Gabriela Calistro Rivera, Rebecca Davies, Elisabete Da Cunha, Tanio Díaz-Santos, Andrea Ferrara, Deanne Fisher, Lilian L. Lee, Juno Li, Dieter Lutz, Monica Relaño, Thorsten Naab , et al. (7 additional authors not shown)

Abstract: We present the morphological parameters and global properties of dust-obscured star formation in typical star-forming galaxies at z=4-6. Among 26 galaxies composed of 20 galaxies observed by the Cycle-8 ALMA Large Program, CRISTAL, and six galaxies from archival data, we have individually detected rest-frame 158$μ$m dust continuum emission from 19 galaxies, nine of which are reported for the first… ▽ More We present the morphological parameters and global properties of dust-obscured star formation in typical star-forming galaxies at z=4-6. Among 26 galaxies composed of 20 galaxies observed by the Cycle-8 ALMA Large Program, CRISTAL, and six galaxies from archival data, we have individually detected rest-frame 158$μ$m dust continuum emission from 19 galaxies, nine of which are reported for the first time. The derived far-infrared luminosities are in the range $\log_{10} L_{\rm IR}\,[L_{\odot}]=$10.9-12.4, an order of magnitude lower than previously detected massive dusty star-forming galaxies (DSFGs). The average relationship between the fraction of dust-obscured star formation ($f_{\rm obs}$) and the stellar mass is consistent with previous results at z=4-6 in a mass range of $\log_{10} M_{\ast}\,[M_{\odot}]\sim$9.5-11.0 and show potential evolution from z=6-9. The individual $f_{\rm obs}$ exhibits a significant diversity, and it shows a correlation with the spatial offset between the dust and the UV continuum, suggesting the inhomogeneous dust reddening may cause the source-to-source scatter in $f_{\rm obs}$. The effective radii of the dust emission are on average $\sim$1.5 kpc and are $\sim2$ times more extended than the rest-frame UV. The infrared surface densities of these galaxies ($Σ_{\rm IR}\sim2.0\times10^{10}\,L_{\odot}\,{\rm kpc}^{-2}$) are one order of magnitude lower than those of DSFGs that host compact central starbursts. On the basis of the comparable contribution of dust-obscured and dust-unobscured star formation along with their similar spatial extent, we suggest that typical star-forming galaxies at z=4-6 form stars throughout the entirety of their disks. △ Less

Submitted 29 November, 2023; originally announced November 2023.

Metal and dust evolution in ALMA REBELS galaxies: insights for future JWST observations

Authors: Marco Palla, Ilse De Looze, Monica Relaño, Stefan van der Giessen, Pratika Dayal, Andrea Ferrara, Raffaella Schneider, Luca Graziani, Hiddo S. B. Algera, Manuel Aravena, Rebecca A. A. Bowler, Alexander P. S. Hygate, Hanae Inami, Ivana van Leeuwen, Rychard Bouwens, Jacqueline Hodge, Renske Smit, Mauro Stefanon, Paul van der Werf

Abstract: ALMA observations revealed the presence of significant amounts of dust in the first Gyr of Cosmic time. However, the metal and dust buildup picture remains very uncertain due to the lack of constraints on metallicity. JWST has started to reveal the metal content of high-redshift targets, which may lead to firmer constraints on high-redshift dusty galaxies evolution. In this work, we use detailed c… ▽ More ALMA observations revealed the presence of significant amounts of dust in the first Gyr of Cosmic time. However, the metal and dust buildup picture remains very uncertain due to the lack of constraints on metallicity. JWST has started to reveal the metal content of high-redshift targets, which may lead to firmer constraints on high-redshift dusty galaxies evolution. In this work, we use detailed chemical and dust evolution models to explore the evolution of galaxies within the ALMA REBELS survey, testing different metallicity scenarios that could be inferred from JWST observations. In the models, we track the buildup of stellar mass by using non-parametric SFHs for REBELS galaxies. Different scenarios for metal and dust evolution are simulated by allowing different prescriptions for gas flows and dust processes. The model outputs are compared with measured dust scaling relations, by employing metallicity-dependent calibrations for the gas mass based on the [CII]158micron line. Independently of the galaxies metal content, we found no need for extreme dust prescriptions to explain the dust masses revealed by ALMA. However, different levels of metal enrichment will lead to different dominant dust production mechanisms, with stardust production dominant over other ISM dust processes only in the metal-poor case. This points out how metallicity measurements from JWST will significantly improve our understanding of the dust buildup in high-redshift galaxies. We also show that models struggle to reproduce observables such as dust-to-gas and dust-to-stellar ratios simultaneously, possibly indicating an overestimation of the gas mass through current calibrations, especially at high metallicities. △ Less

Submitted 27 November, 2023; originally announced November 2023.

Comments: 16 pages + appendices, 9 Figures, 1 Table. Resubmitted to MNRAS after moderate revision

Stellar mass-metallicity relation throughout the large-scale structure of the Universe: CAVITY mother sample

Authors: Jesús Domínguez-Gómez, Isabel Pérez, Tomás Ruiz-Lara, Reynier F. Peletier, Patricia Sánchez-Blázquez, Ute Lisenfeld, Bahar Bidaran, Jesús Falcón-Barroso, Manuel Alcázar-Laynez, María Argudo-Fernández, Guillermo Blázquez-Calero, Hélène Courtois, Salvador Duarte Puertas, Daniel Espada, Estrella Florido, Rubén García-Benito, Andoni Jiménez, Kathryn Kreckel, Mónica Relaño, Laura Sánchez-Menguiano, Thijs van der Hulst, Rien van de Weygaert, Simon Verley, Almudena Zurita

Abstract: Void galaxies are essential for understanding the physical processes that drive galaxy evolution because they are less affected by external factors than galaxies in denser environments, that is, in filaments, walls, and clusters. The stellar metallicity of a galaxy traces the accumulated fossil record of the star formation through the entire life of the galaxy. A comparison of the stellar metallic… ▽ More Void galaxies are essential for understanding the physical processes that drive galaxy evolution because they are less affected by external factors than galaxies in denser environments, that is, in filaments, walls, and clusters. The stellar metallicity of a galaxy traces the accumulated fossil record of the star formation through the entire life of the galaxy. A comparison of the stellar metallicity of galaxies in various environments, including voids, filaments, walls, and clusters can provide valuable insights into how the large-scale environment affects the chemical evolution of the galaxy. We present the first comparison of the relation of the total stellar mass versus central stellar metallicity between galaxies in voids, filaments, walls, and clusters with different star formation history (SFH) types, morphologies, and colours for stellar masses between $10^{8.0}$ to $10^{11.5}$ solar masses and redshift $0.01<z<0.05$. We applied non-parametric full spectral fitting techniques (pPXF and STECKMAP) to 10807 spectra from the SDSS-DR7 (987 in voids, 6463 in filaments and walls, and 3357 in clusters) and derived their central mass-weighted average stellar metallicity ($\rm [M/H]_M$). We find that galaxies in voids have slightly lower stellar metallicities on average than galaxies in filaments and walls (by~$\sim~0.1$~dex), and they are much lower than those of galaxies in clusters (by~$\sim~0.4$~dex). These differences are more significant for low-mass ($ \sim~10^{9.25}~{\rm M_\odot}$) than for high-mass galaxies, for long-timescale SFH (extended along time) galaxies than for short-timescale SFHs (concentrated at early times) galaxies, for spiral than for elliptical galaxies, and for blue than for red galaxies. △ Less

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

Comments: Accepted in A&A, 39 pages, 18 figures, 24 tables. Language corrected version. Title change and tables correction

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

Galaxies in voids assemble their stars slowly

Authors: J. Domínguez-Gómez, I. Pérez, T. Ruiz-Lara, R. F. Peletier, P. Sánchez-Blázquez, U. Lisenfeld, J. Falcón-Barroso, M. Alcázar-Laynez, M. Argudo-Fernández, G. Blázquez-Calero, H. Courtois, S. Duarte Puertas, D. Espada, E. Florido, R. García-Benito, A. Jiménez, K. Kreckel, M. Relaño, L. Sánchez-Menguiano, T. van der Hulst, R. van de Weygaert, S. Verley, A. Zurita

Abstract: Galaxies in the Universe are distributed in a web-like structure characterised by different large-scale environments: dense clusters, elongated filaments, sheetlike walls, and under-dense regions, called voids. The low density in voids is expected to affect the properties of their galaxies. Indeed, previous studies have shown that galaxies in voids are on average bluer and less massive, and have l… ▽ More Galaxies in the Universe are distributed in a web-like structure characterised by different large-scale environments: dense clusters, elongated filaments, sheetlike walls, and under-dense regions, called voids. The low density in voids is expected to affect the properties of their galaxies. Indeed, previous studies have shown that galaxies in voids are on average bluer and less massive, and have later morphologies and higher current star formation rates than galaxies in denser large-scale environments. However, it has never been observationally proved that the star formation histories (SFHs) in void galaxies are substantially different from those in filaments, walls, and clusters. Here we show that void galaxies have had, on average, slower SFHs than galaxies in denser large-scale environments. We also find two main SFH types present in all the environments: 'short-timescale' galaxies are not affected by their large-scale environment at early times but only later in their lives; 'long-timescale' galaxies have been continuously affected by their environment and stellar mass. Both types have evolved slower in voids than in filaments, walls, and clusters. △ Less

Submitted 29 June, 2023; originally announced June 2023.

Comments: Published in Nature. 27 pages, 4 Figures, 7 Extended Data Figures, and 1 Extended Data Table

SOFIA/HAWC+ observations of the Crab Nebula: dust properties from polarised emission

Authors: Jérémy Chastenet, Ilse De Looze, Brandon S. Hensley, Bert Vandenbroucke, Mike J. Barlow, Jeonghee Rho, Aravind P. Ravi, Haley L. Gomez, Florian Kirchschlager, Juan Macías-Pérez, Mikako Matsuura, Kate Pattle, Nicolas Ponthieu, Felix D. Priestley, Monica Relaño, Alessia Ritacco, Roger Wesson

Abstract: Supernova remnants (SNRs) are well-recognised dust producers, but their net dust production rate remains elusive due to uncertainties in grain properties that propagate into observed dust mass uncertainties, and determine how efficiently these grains are processed by reverse shocks. In this paper, we present a detection of polarised dust emission in the Crab pulsar wind nebula, the second SNR with… ▽ More Supernova remnants (SNRs) are well-recognised dust producers, but their net dust production rate remains elusive due to uncertainties in grain properties that propagate into observed dust mass uncertainties, and determine how efficiently these grains are processed by reverse shocks. In this paper, we present a detection of polarised dust emission in the Crab pulsar wind nebula, the second SNR with confirmed polarised dust emission after Cassiopeia A. We constrain the bulk composition of the dust with new SOFIA/HAWC+ polarimetric data in band C 89 um and band D 154 um. After correcting for synchrotron polarisation, we report dust polarisation fractions ranging between 3.7-9.6 per cent and 2.7-7.6 per cent in three individual dusty filaments at 89 and 154 um, respectively. The detected polarised signal suggests the presence of large (> 0.05-0.1 um) grains in the Crab Nebula. With the observed polarisation, and polarised and total fluxes, we constrain the temperatures and masses of carbonaceous and silicate grains. We find that the carbon-rich grain mass fraction varies between 12 and 70 per cent, demonstrating that carbonaceous and silicate grains co-exist in this SNR. Temperatures range from 40 K to 70 K and from 30 K to 50 K for carbonaceous and silicate grains, respectively. Dust masses range from 10^{-4} Msol to 10^{-2} Msol for carbonaceous grains and to 10^{-1} Msol for silicate grains, in three individual regions. △ Less

Submitted 23 August, 2022; originally announced August 2022.

Comments: Accepted for publication in MNRAS. 16 pages, 8 figures + 2 Appendices, 2 figures

Dust grain size evolution in local galaxies: a comparison between observations and simulations

Authors: M. Relano, I. De Looze, A. Saintonge, K. -C. Hou, L. Romano, K. Nagamine, H. Hirashita, S. Aoyama, I. Lamperti, U. Lisenfeld, M. Smith, J. Chastenet, T. Xiao, Y. Gao, M. Sargent, S. A. van der Giessen

Abstract: The evolution of the dust grain size distribution has been studied in recent years with great detail in cosmological hydrodynamical simulations taking into account all the channels under which dust evolves in the interstellar medium. We present a systematic analysis of the observed spectral energy distribution of a large sample of galaxies in the local universe in order to derive not only the tota… ▽ More The evolution of the dust grain size distribution has been studied in recent years with great detail in cosmological hydrodynamical simulations taking into account all the channels under which dust evolves in the interstellar medium. We present a systematic analysis of the observed spectral energy distribution of a large sample of galaxies in the local universe in order to derive not only the total dust masses but also the relative mass fraction between small and large dust grains (DS/DL). Simulations reproduce fairly well the observations except for the high stellar mass regime where dust masses tend to be overestimated. We find that ~45% of galaxies exhibit DS/DL consistent with the expectations of simulations, while there is a sub-sample of massive galaxies presenting high DS/DL (log(DS/DL)~-0.5), and deviating from the prediction in simulations. For these galaxies, which also have high molecular gas mass fractions and metallicities, coagulation is not an important mechanism affecting the dust evolution. Including diffusion, transporting large grains from dense regions to a more diffuse medium where they can be easily shattered, would explain the observed high DS/DL values in these galaxies. With this study we reinforce the use of the small-to-large grain mass ratio to study the relative importance of the different mechanisms in the dust life cycle. Multi-phase hydrodynamical simulations with detailed feedback prescriptions and more realistic subgrid models for the dense phase could help to reproduce the evolution of the dust grain size distribution traced by observations. △ Less

Submitted 26 July, 2022; originally announced July 2022.

Comments: 32 pages, 22 figures, 4 tables. Accepted in MNRAS

JINGLE -- IV. Dust, HI gas and metal scaling laws in the local Universe

Authors: I. De Looze, I. Lamperti, A. Saintonge, M. Relano, M. W. L. Smith, C. J. R. Clark, C. D. Wilson, M. Decleir, A. P. Jones, R. C. Kennicutt, G. Accurso, E. Brinks, M. Bureau, P. Cigan, D. L. Clements, P. De Vis, L Fanciullo, Y. Gao, W. K. Gear, L. C. Ho, H. S. Hwang, M. J. Michalowski, J. C. Lee, C. Li, L. Lin , et al. (7 additional authors not shown)

Abstract: Scaling laws of dust, HI gas and metal mass with stellar mass, specific star formation rate and metallicity are crucial to our understanding of the buildup of galaxies through their enrichment with metals and dust. In this work, we analyse how the dust and metal content varies with specific gas mass ($M_{\text{HI}}$/$M_{\star}$) across a diverse sample of 423 nearby galaxies. The observed trends a… ▽ More Scaling laws of dust, HI gas and metal mass with stellar mass, specific star formation rate and metallicity are crucial to our understanding of the buildup of galaxies through their enrichment with metals and dust. In this work, we analyse how the dust and metal content varies with specific gas mass ($M_{\text{HI}}$/$M_{\star}$) across a diverse sample of 423 nearby galaxies. The observed trends are interpreted with a set of Dust and Element evolUtion modelS (DEUS) - incluidng stellar dust production, grain growth, and dust destruction - within a Bayesian framework to enable a rigorous search of the multi-dimensional parameter space. We find that these scaling laws for galaxies with $-1.0\lesssim \log M_{\text{HI}}$/$M_{\star}\lesssim0$ can be reproduced using closed-box models with high fractions (37-89$\%$) of supernova dust surviving a reverse shock, relatively low grain growth efficiencies ($ε$=30-40), and long dus lifetimes (1-2\,Gyr). The models have present-day dust masses with similar contributions from stellar sources (50-80\,$\%$) and grain growth (20-50\,$\%$). Over the entire lifetime of these galaxies, the contribution from stardust ($>$90\,$\%$) outweighs the fraction of dust grown in the interstellar medium ($<$10$\%$). Our results provide an alternative for the chemical evolution models that require extremely low supernova dust production efficiencies and short grain growth timescales to reproduce local scaling laws, and could help solving the conundrum on whether or not grains can grow efficiently in the interstellar medium. △ Less

Submitted 2 June, 2020; originally announced June 2020.

Comments: 56 pages, 30 figures, accepted for publication in MNRAS

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)

Evolution of the grain size distribution in galactic discs

Authors: M. Relano, U. Lisenfeld, K. C. Hou, I. De Looze, J. M. Vilchez, R. C. Kennicutt

Abstract: Dust is formed out of stellar material and is constantly affected by different mechanisms occurring in the ISM. Dust grains behave differently under these mechanisms depending on their sizes, and therefore the dust grain size distribution also evolves as part of the dust evolution itself. Following how the grain size distribution evolves is a difficult computing task that is just recently being ov… ▽ More Dust is formed out of stellar material and is constantly affected by different mechanisms occurring in the ISM. Dust grains behave differently under these mechanisms depending on their sizes, and therefore the dust grain size distribution also evolves as part of the dust evolution itself. Following how the grain size distribution evolves is a difficult computing task that is just recently being overtaking. Smoothed particle hydrodynamic (SPH) simulations of a single galaxy as well as cosmological simulations are producing the first predictions of the evolution of the dust grain size distribution. We compare for the first time the evolution of the dust grain size distribution predicted by the SPH simulations with the results provided by the observations. We analyse how the radial distribution of the small to large grain mass ratio (D(S)/D(L)) changes over the whole discs in three galaxies: M 101, NGC 628 and M 33. We find good agreement between the observed radial distribution of D(S)/D(L) and what is obtained from the SPH simulations of a single galaxy. The central parts of NGC 628, at high metallicity and with a high molecular gas fraction, are mainly affected not only by accretion but also by coagulation of dust grains. The centre of M 33, having lower metallicity and lower molecular gas fraction, presents an increase of D(S)/D(L), showing that shattering is very effective in creating a large fraction of small grains. Observational results provided by our galaxies confirm the general relations predicted by the cosmological simulations based on the two grain size approximation. However, we present evidence that the simulations could be overestimating the amount of large grains in high massive galaxies. △ Less

Submitted 5 February, 2020; originally announced February 2020.

Comments: 11 pages, 6 Figures, to be published in A&A

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

Modeling Dust and Starlight in Galaxies Observed by Spitzer and Herschel: The KINGFISH Sample

Authors: G. Aniano, B. T. Draine, L. K. Hunt, K. Sandstrom, D. Calzetti, R. C. Kennicutt, D. A. Dale, M. Galametz, K. D. Gordon, A. K. Leroy, J. -D. T. Smith, H. Roussel, M. Sauvage, F. Walter, L. Armus, A. D. Bolatto, M. Boquien, A. Crocker, I. De Looze, J. Donovan Meyer, G. Helou, J. Hinz, B. D. Johnson, J. Koda, A. Miller , et al. (8 additional authors not shown)

Abstract: Dust and starlight are modeled for the KINGFISH project galaxies. With data from 3.6 micron to 500 micron, models are strongly constrained. For each pixel in each galaxy we estimate (1) dust surface density; (2) q_PAH, the dust mass fraction in PAHs; (3) distribution of starlight intensities heating the dust; (4) luminosity emitted by the dust; and (5) dust luminosity from regions with high starli… ▽ More Dust and starlight are modeled for the KINGFISH project galaxies. With data from 3.6 micron to 500 micron, models are strongly constrained. For each pixel in each galaxy we estimate (1) dust surface density; (2) q_PAH, the dust mass fraction in PAHs; (3) distribution of starlight intensities heating the dust; (4) luminosity emitted by the dust; and (5) dust luminosity from regions with high starlight intensity. The models successfully reproduce both global and resolved spectral energy distributions. We provide well-resolved maps for the dust properties. As in previous studies, we find q_PAH to be an increasing function of metallicity, above a threshold Z/Z_sol approx 0.15. Dust masses are obtained by summing the dust mass over the map pixels; these "resolved" dust masses are consistent with the masses inferred from model fits to the global photometry. The global dust-to-gas ratios obtained from this study correlate with galaxy metallicities. Systems with Z/Z_sol > 0.5 have most of their refractory elements locked up in dust, whereas when Z/Z_sol < 0.3 most of these elements tend to remain in the gas phase. Within galaxies, we find that q_PAH is suppressed in regions with unusually warm dust with nu L_nu(70 um) > 0.4L_dust. With knowledge of one long-wavelength flux density ratio (e.g., f_{160}/f_{500}), the minimum starlight intensity heating the dust (U_min) can be estimated to within ~50%. For the adopted dust model, dust masses can be estimated to within ~0.07 dex accuracy using the 500 micron luminosity nu L_nu(500) alone. There are additional systematic errors arising from the choice of dust model, but these are hard to estimate. These calibrated prescriptions may be useful for studies of high-redshift galaxies. △ Less

Submitted 10 December, 2019; originally announced December 2019.

Comments: ApJ, accepted

High-Resolution Radiative Transfer Modelling of M33

Authors: Thomas G. Williams, Maarten Baes, Ilse De Looze, Monica Relaño, Matthew W. L. Smith, Sam Verstocken, Sébastien Viaene

Abstract: In this work, we characterise the contributions from both ongoing star formation and the ambient radiation field in Local Group galaxy M33, as well as estimate the scale of the local dust-energy balance (i.e. the scale at which the dust is re-emitting starlight generated in that same region) in this galaxy through high-resolution radiative transfer (RT) modelling, with defined stellar and dust geo… ▽ More In this work, we characterise the contributions from both ongoing star formation and the ambient radiation field in Local Group galaxy M33, as well as estimate the scale of the local dust-energy balance (i.e. the scale at which the dust is re-emitting starlight generated in that same region) in this galaxy through high-resolution radiative transfer (RT) modelling, with defined stellar and dust geometries. We have characterised the spectral energy distribution (SED) of M33 from UV to sub-mm wavelengths, at a spatial scale of 100 pc. We constructed input maps of the various stellar and dust geometries for use in the RT modelling. By modifying our dust mix (fewer very small carbon grains and a lower silicate-to-carbon ratio as compared to the Milky Way), we can much better fit the sub-mm dust continuum. Using this new dust composition, we find that we are able to well reproduce the observed SED of M33 using our adopted model. In terms of stellar attenuation by dust, we find a reasonably strong, broad UV bump, as well as significant systematic differences in the amount of dust attenuation when compared to standard SED modelling. We also find discrepancies in the residuals of the spiral arms versus the diffuse interstellar medium (ISM), indicating a difference in properties between these two regimes. The dust emission is dominated by heating due to the young stellar populations at all wavelengths ($\sim$80% at 10 $μ$m to $\sim$50% at 1 mm). We find that the local dust-energy balance is restored at spatial scales greater than around 1.5 kpc. △ Less

Submitted 5 June, 2019; v1 submitted 23 May, 2019; originally announced May 2019.

Comments: 19 pages, 14 Figures, 2 Tables. Accepted for publication in MNRAS. Included missing reference

Metals and dust content across the galaxies M101 and NGC 628

Authors: J. M. Vilchez, M. Relano, R. Kennicutt, I. De Looze, M. Molla, M. Galametz

Abstract: We present a spatially resolved study of the relation between dust and metallicity in the nearby spiral galaxies M101 (NGC 5457) and NGC 628 (M74). We explore the relation between the chemical abundances of their gas and stars with their dust content and their chemical evolution. The empirical spatially resolved oxygen effective yield and the gas to dust mass ratio (GDR) across both disc galaxies… ▽ More We present a spatially resolved study of the relation between dust and metallicity in the nearby spiral galaxies M101 (NGC 5457) and NGC 628 (M74). We explore the relation between the chemical abundances of their gas and stars with their dust content and their chemical evolution. The empirical spatially resolved oxygen effective yield and the gas to dust mass ratio (GDR) across both disc galaxies are derived, sampling one dex in oxygen abundance. We find that the metal budget of the NGC 628 disc and most of the M101 disc appears consistent with the predictions of the simple model of chemical evolution for an oxygen yield between half and one solar, whereas the outermost region (R<0.8R25) of M101 presents deviations suggesting the presence of gas flows. The GDR-metallicity relation shows a two slopes behaviour, with a break at 12+log(O/H)~8.4, a critical metallicity predicted by theoretical dust models when stardust production equals grain growth. A relation between GDR and the fraction of molecular to total gas, Sigma(H2)/Sigma(gas) is also found. We suggest an empirical relationship between GDR and the combination of 12+log(O/H), for metallicity, and Sigma(H2)/Sigma(gas), a proxy for the molecular clouds fraction. The GDR is closely related with metallicity at low abundance and with Sigma(H2)/Sigma(gas) for higher metallicities suggesting ISM dust growth. The ratio Sigma(dust)/Sigma(star) correlates well with 12 + log(O/H) and strongly with log(N/O) in both galaxies. For abundances below the critical one, the 'stardust' production gives us a constant value suggesting a stellar dust yield similar to the oxygen yield. △ Less

Submitted 26 November, 2018; originally announced November 2018.

Comments: 16 pages, 14 Figures, accepted in MNRAS

Comprehensive comparison of models for spectral energy distributions from 0.1 micron to 1 mm of nearby star-forming galaxies

Authors: L. K. Hunt, I. De Looze, M. Boquien, R. Nikutta, A. Rossi, S. Bianchi, D. A. Dale, G. L. Granato, R. C. Kennicutt, L. Silva, L. Ciesla, M. Relano, S. Viaene, B. Brandl, D. Calzetti, K. V. Croxall, B. T. Draine, M. Galametz, K. D. Gordon, B. A. Groves, G. Helou, R. Herrera-Camus, J. L. Hinz, J. Koda, S. Salim , et al. (4 additional authors not shown)

Abstract: We have fit the far-ultraviolet (FUV) to sub-millimeter (850 micron) spectral energy distributions (SEDs) of the 61 galaxies from the "Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel" (KINGFISH). The fitting has been performed using three models: the Code for Investigating GALaxy Evolution (CIGALE), the GRAphite-SILicate approach (GRASIL), and the Multi-wavelength Analysis of… ▽ More We have fit the far-ultraviolet (FUV) to sub-millimeter (850 micron) spectral energy distributions (SEDs) of the 61 galaxies from the "Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel" (KINGFISH). The fitting has been performed using three models: the Code for Investigating GALaxy Evolution (CIGALE), the GRAphite-SILicate approach (GRASIL), and the Multi-wavelength Analysis of Galaxy PHYSical properties (MAGPHYS). We have analyzed the results of the three codes in terms of the SED shapes, and by comparing the derived quantities with simple "recipes" for stellar mass (Mstar), star-formation rate (SFR), dust mass (Mdust), and monochromatic luminosities. Although the algorithms rely on different assumptions for star-formation history, dust attenuation and dust reprocessing, they all well approximate the observed SEDs and are in generally good agreement for the associated quantities. However, the three codes show very different behavior in the mid-infrared regime, in particular between 25 and 70 micron where there are no observational constraints for the KINGFISH sample. We find that different algorithms give discordant SFR estimates for galaxies with low specific SFR, and that the standard "recipes" for calculating FUV absorption overestimate the extinction compared to the SED-fitting results. Results also suggest that assuming a "standard" constant stellar mass-to-light ratio overestimates Mstar relative to the SED fitting, and we provide new SED-based formulations for estimating Mstar from WISE W1 (3.4 micron) luminosities and colors. From a Principal Component Analysis of Mstar, SFR, Mdust, and O/H, we reproduce previous scaling relations among Mstar, SFR, and O/H, and find that Mdust can be predicted to within roughly 0.3 dex using only Mstar and SFR. △ Less

Submitted 13 November, 2018; v1 submitted 11 September, 2018; originally announced September 2018.

Comments: 42 pages, 23 figures, 8 tables, 2 appendices: updated to v2, A&A, accepted

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

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)

The Wolf-Rayet star population in the dwarf galaxy NGC 625

Authors: Ana Monreal-Ibero, Jeremy R. Walsh, Jorge Iglesias-Paramo, Christer Sandin, M. Relano, E. Perez-Montero, J. Vilchez

Abstract: ABRIGED: Quantifying the number, type and distribution of W-R stars is a key component in the context of galaxy evolution, since they put constraints on the age of the star formation bursts. Nearby galaxies (d<5 Mpc) are particularly relevant in this context since they fill the gap between studies in the Local Group, where individual stars can be resolved, and galaxies in the Local Volume and beyo… ▽ More ABRIGED: Quantifying the number, type and distribution of W-R stars is a key component in the context of galaxy evolution, since they put constraints on the age of the star formation bursts. Nearby galaxies (d<5 Mpc) are particularly relevant in this context since they fill the gap between studies in the Local Group, where individual stars can be resolved, and galaxies in the Local Volume and beyond. We intend to characterize the W-R star population in NGC625, a low-metallicity dwarf galaxy suffering a currently declining burst of star formation. Optical IFS data have been obtained with the VIMOS-IFU covering the starburst region. We estimate the number of W-R stars using a linear combination of 3 W-R templates: 1 early-type nitrogen (WN) star, 1 late-type WN star and 1 carbon-type (WC) star (or oxygen-type (WO) star). Fits using several ensembles of templates were tested. Results were confronted with: i) high spatial resolution HST photometry; ii) numbers of W-R stars in nearby galaxies; iii) model predictions. The W-R star population is spread over the main body of the galaxy, not necessarily coincident with the overall stellar distribution. Our best estimation for the number of W-R stars yields a total of 28 W-R stars in the galaxy, out of which 17 are early- type WN, 6 are late-type WN and 5 are WC stars. The width of the stellar features nicely correlates with the dominant W-R type found in each aperture. The distribution of the different types of WR in the galaxy is roughly compatible with the way star formation has propagated in the galaxy, according to previous findings using HST images. Fits using templates at the metallicity of the LMC yield more reasonable number of W-R than those using templates at the metallicity of the SMC. Given the metallicity of NGC 625, this suggests a non-linear relation between the metallicity and the luminosity of the W-R spectral features. △ Less

Submitted 19 April, 2017; originally announced April 2017.

Comments: 13 pages, 6 figures, accepted by A&A; Labels in the axes of some figures are not displayed properly after the compilation on arXiv. Another version with the proper figures can be found here: https://www.dropbox.com/s/u2bx2li5jogmfkp/manuscript_astroph.pdf?dl=0

Journal ref: A&A 603, A130 (2017)

Dust Properties in HII Regions in M33

Authors: M. Relano, R. Kennicutt, U. Lisenfeld, S. Verley, I. Hermelo, M. Boquien, M. Albrecht, C. Kramer, J. Braine, E. Perez-Montero, I. De Looze, M. Xilouris, A. Kovacs, J. Staguhn

Abstract: The conversion of the IR emission into star formation rate can be strongly dependent on the physical properties of the dust, which are affected by the environmental conditions where the dust is embedded. We study here the dust properties of a set of HII regions in the Local Group Galaxy M33 presenting different spatial configurations between the stars, gas and dust to understand the dust evolution… ▽ More The conversion of the IR emission into star formation rate can be strongly dependent on the physical properties of the dust, which are affected by the environmental conditions where the dust is embedded. We study here the dust properties of a set of HII regions in the Local Group Galaxy M33 presenting different spatial configurations between the stars, gas and dust to understand the dust evolution under different environments. We model the SED of each region using the DustEM tool and obtain the mass relative to hydrogen for Very Small Grains (YVSG), Polycyclic Aromatic Hydrocarbons (YPAH) and Big Grains (YBG). The relative mass of the VSGs (YVSG/YTOT) is a factor of 1.7 higher for HII regions classified as filled and mixed than for regions presenting a shell structure. The enhancement of VSGs within NGC 604 and NGC 595 is correlated to expansive gas structures with velocities greater than 50 km/s. The gas-to-dust ratio derived for the HII regions in our sample exhibits two regimes related to the HI-H2 transition of the ISM. Regions corresponding to the HI diffuse regime present a gas-to-dust ratio compatible with the expected value if we assume that the gas-to-dust ratio scales linearly with metallicity, while regions corresponding to a H2 molecular phase present a flatter dust-gas surface density distribution. The fraction of VSGs can be affected by the conditions of the interstellar environment: strong shocks of 50-90 km/s existing in the interior of the most luminous HII regions can lead to fragmentation of BGs into smaller ones, while the more evolved shell and clear shell objects provide a more quiescent environment where reformation of dust BG grains might occur. The gas-to-dust variations found in this analysis might imply that grain coagulation and/or gas-phase metals incorporation to the dust mass is occurring in the interior of the HII regions in M33. △ Less

Submitted 11 June, 2016; originally announced June 2016.

Comments: Accepted for publication in A&A, 29 pages, 16 figures

Journal ref: A&A 595, A43 (2016)

Towards universal hybrid star formation rate estimators

Authors: M. Boquien, R. Kennicutt, D. Calzetti, D. Dale, M. Galametz, M. Sauvage, K. Croxall, B. Draine, A. Kirkpatrick, N. Kumari, L. Hunt, I. De Looze, E. Pellegrini, M. Relano, J. -D. Smith, F. Tabatabaei

Abstract: To compute the SFR of galaxies from the rest-frame UV it is essential to take into account the obscuration by dust. To do so, one of the most popular methods consists in combining the UV with the emission from the dust itself in the IR. Yet, different studies have derived different estimators, showing that no such hybrid estimator is truly universal. In this paper we aim at understanding and quant… ▽ More To compute the SFR of galaxies from the rest-frame UV it is essential to take into account the obscuration by dust. To do so, one of the most popular methods consists in combining the UV with the emission from the dust itself in the IR. Yet, different studies have derived different estimators, showing that no such hybrid estimator is truly universal. In this paper we aim at understanding and quantifying what physical processes drive the variations between different hybrid estimators. Doing so, we aim at deriving new universal UV+IR hybrid estimators to correct the UV for dust attenuation, taking into account the intrinsic physical properties of galaxies. We use the CIGALE code to model the spatially-resolved FUV to FIR SED of eight nearby star-forming galaxies drawn from the KINGFISH sample. This allows us to determine their local physical properties, and in particular their UV attenuation, average SFR, average specific SFR (sSFR), and their stellar mass. We then examine how hybrid estimators depend on said properties. We find that hybrid UV+IR estimators strongly depend on the stellar mass surface density (in particular at 70 and 100 micron) and on the sSFR (in particular at 24 micron and the TIR). Consequently, the IR scaling coefficients for UV obscuration can vary by almost an order of magnitude. This result contrasts with other groups who found relatively constant coefficients with small deviations. We exploit these variations to construct a new class of hybrid estimators based on observed UV to near-IR colours and near-IR luminosity densities per unit area. We find that they can reliably be extended to entire galaxies. The new estimators provide better estimates of attenuation-corrected UV emission than classical hybrid estimators. Naturally taking into account the variable impact of dust heated by old stellar populations, they constitute a step towards universal estimators. △ Less

Submitted 30 March, 2016; originally announced March 2016.

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

Millimeter and Submillimeter Excess Emission in M33 revealed by Planck and LABOCA

Authors: I. Hermelo, M. Relaño, U. Lisenfeld, S. Verley, C. Kramer, T. Ruiz-Lara, M. Boquien, E. M. Xilouris, M. Albrecht

Abstract: Previous studies have shown the existence of an excess of emission at submillimeter (submm) and millimeter (mm) wavelengths in the spectral energy distribution (SED) of many low-metallicity galaxies. The goal of the present study is to model separately the emission from the star forming (SF) component and the emission from the diffuse interstellar medium (ISM) in the nearby spiral galaxy M33. We d… ▽ More Previous studies have shown the existence of an excess of emission at submillimeter (submm) and millimeter (mm) wavelengths in the spectral energy distribution (SED) of many low-metallicity galaxies. The goal of the present study is to model separately the emission from the star forming (SF) component and the emission from the diffuse interstellar medium (ISM) in the nearby spiral galaxy M33. We decomposed the observed SED of M33 into its SF and diffuse components. Mid-infrared (MIR) and far-infrared (FIR) fluxes were extracted from Spitzer and Herschel data. At submm and mm wavelengths, we used ground-based observations from APEX to measure the emission from the SF component and data from the Planck space telescope to estimate the diffuse emission. Both components were separately fitted using radiation transfer models based on standard dust properties and a realistic geometry. Both modeled SEDs were combined to build the global SED of M33. In addition, the radiation field necessary to power the dust emission in our modeling was compared with observations from GALEX, Sloan, and Spitzer. Our modeling is able to reproduce the observations at MIR and FIR wavelengths, but we found a strong excess of emission at submm and mm wavelengths, where the model expectations severely underestimate the LABOCA and Planck fluxes. We also found that the ultraviolet (UV) radiation escaping the galaxy is 70% higher than the model predictions. We determined a gas-to-dust mass ratio Gdust~100, significantly lower than the value expected from the sub-solar metallicity of M33. We discussed different hypotheses to explain the discrepancies found in our study (i.e., excess of emission at submm and mm wavelengths, deficit of UV attenuation, and abnormally low value of Gdust), concluding that different dust properties in M33 is the most plausible explanation. △ Less

Submitted 7 March, 2016; originally announced March 2016.

Comments: accepted by A&A, 22 pages

Journal ref: A&A 590, A56 (2016)

Measuring star formation with resolved observations: the test case of M33

Authors: M. Boquien, D. Calzetti, S. Aalto, A. Boselli, J. Braine, V. Buat, F. Combes, F. Israel, C. Kramer, S. Lord, M. Relano, E. Rosolowsky, G. Stacey, F. Tabatabaei, F. van der Tak, P. van der Werf, S. Verley, M. Xilouris

Abstract: Context. Measuring star formation at a local scale is important to constrain star formation laws. Yet, it is not clear whether and how the measure of star formation is affected by the spatial scale at which a galaxy is observed. Aims. We want to understand the impact of the resolution on the determination of the spatially resolved star formation rate (SFR) and other directly associated physical pa… ▽ More Context. Measuring star formation at a local scale is important to constrain star formation laws. Yet, it is not clear whether and how the measure of star formation is affected by the spatial scale at which a galaxy is observed. Aims. We want to understand the impact of the resolution on the determination of the spatially resolved star formation rate (SFR) and other directly associated physical parameters such as the attenuation. Methods. We have carried out a multi-scale, pixel-by-pixel study of the nearby galaxy M33. Assembling FUV, Halpha, 8, 24, 70, and 100 micron maps, we have systematically compared the emission in individual bands with various SFR estimators from a resolution of 33 pc to 2084 pc. Results. We have found that there are strong, scale-dependent, discrepancies up to a factor 3 between monochromatic SFR estimators and Halpha+24 micron. The scaling factors between individual IR bands and the SFR show a strong dependence on the spatial scale and on the intensity of star formation. Finally, strong variations of the differential reddening between the nebular emission and the stellar continuum are seen, depending on the specific SFR (sSFR) and on the resolution. At the finest spatial scales, there is little differential reddening at high sSFR. The differential reddening increases with decreasing sSFR. At the coarsest spatial scales the differential reddening is compatible with the canonical value found for starburst galaxies. Conclusions. Our results confirm that monochromatic estimators of the SFR are unreliable at scales smaller than 1 kpc. Furthermore, the extension of local calibrations to high redshift galaxies presents non-trivial challenges as the properties of these systems may be poorly known. △ Less

Submitted 4 February, 2015; originally announced February 2015.

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

The giant HII region NGC 588 as a benchmark for 2D photoionisation models

Authors: Enrique Perez-Montero, Ana Monreal-Ibero, Monica Relano, Jose M. Vilchez, Carolina Kehrig, Christophe Morisset

Abstract: We use optical integral field spectroscopy and 8 and 24 micron mid-IR observations of the giant HII region NGC 588 in the disc of M33 as input and constraints for two-dimensional tailor-made photoionisation models. Two different geometrical approaches are followed for the modelling structure: i) Each spatial element of the emitting gas is studied individually using models which assume that the ion… ▽ More We use optical integral field spectroscopy and 8 and 24 micron mid-IR observations of the giant HII region NGC 588 in the disc of M33 as input and constraints for two-dimensional tailor-made photoionisation models. Two different geometrical approaches are followed for the modelling structure: i) Each spatial element of the emitting gas is studied individually using models which assume that the ionisation structure is complete in each element to look for azimuthal variations across gas and dust. ii) A single model is considered, and the two-dimensional structure of the gas and the dust are assumed to be due to the projection of an emitting sphere onto the sky. The models in both assumptions reproduce the radial profiles of Hbeta surface brightness, the observed number of ionising photons, and the strong optical emission-line relative intensities. The first approach produces a constant-density matter-bounded thin shell of variable thickness and dust-to-gas ratio, while the second gives place to a radiation-bounded thick shell sphere of decreasing particle density. However, the radial profile of the 8/24 microns IR ratio, depending on the gas and dust geometry, only fits well when the thick-shell model is used. The resulting dust-to-gas mass ratio, which was obtained empirically from the derived dust mass using data from Spitzer, also has a better fit using the thick-shell solution. In both approaches, models support the chemical homogeneity, and the ionisation-parameter radial decrease, These results must be taken with care in view of the very low extinction values that are derived from the IR, as compared to that derived from the Balmer decrement. Besides, the IR can be possibly contaminated with the emission from a cloud of diffuse gas and dust above the plane of the galaxy detected at 250 micron Herschel image. △ Less

Submitted 9 April, 2014; v1 submitted 8 April, 2014; originally announced April 2014.

Comments: 16 pages, 14 figures. Accepted for publication in Astronomy & Astrophysics. Fig. 1 included

The effects of spatial resolution on Integral Field Spectrograph surveys at different redshifts. The CALIFA perspective

Authors: D. Mast, F. F. Rosales-Ortega, S. F. Sanchez, J. M. Vílchez, J. Iglesias-Paramo, C. J. Walcher, B. Husemann, I. Marquez, R. A. Marino, R. C. Kennicutt, A. Monreal-Ibero, L. Galbany, A. de Lorenzo-Caceres, J. Mendez-Abreu, C. Kehrig, A. del Olmo, M. Relano, L. Wisotzki, E. Marmol-Queralto, S. Bekeraite, P. Papaderos, V. Wild, J. A. L. Aguerri, J. Falcon-Barroso, D. J. Bomans , et al. (5 additional authors not shown)

Abstract: Over the past decade, 3D optical spectroscopy has become the preferred tool for understanding the properties of galaxies and is now increasingly used to carry out galaxy surveys. Low redshift surveys include SAURON, DiskMass, ATLAS3D, PINGS and VENGA. At redshifts above 0.7, surveys such as MASSIV, SINS, GLACE, and IMAGES have targeted the most luminous galaxies to study mainly their kinematic pro… ▽ More Over the past decade, 3D optical spectroscopy has become the preferred tool for understanding the properties of galaxies and is now increasingly used to carry out galaxy surveys. Low redshift surveys include SAURON, DiskMass, ATLAS3D, PINGS and VENGA. At redshifts above 0.7, surveys such as MASSIV, SINS, GLACE, and IMAGES have targeted the most luminous galaxies to study mainly their kinematic properties. The on-going CALIFA survey ($z\sim0.02$) is the first of a series of upcoming Integral Field Spectroscopy (IFS) surveys with large samples representative of the entire population of galaxies. Others include SAMI and MaNGA at lower redshift and the upcoming KMOS surveys at higher redshift. Given the importance of spatial scales in IFS surveys, the study of the effects of spatial resolution on the recovered parameters becomes important. We explore the capability of the CALIFA survey and a hypothetical higher redshift survey to reproduce the properties of a sample of objects observed with better spatial resolution at lower redshift. Using a sample of PINGS galaxies, we simulate observations at different redshifts. We then study the behaviour of different parameters as the spatial resolution degrades with increasing redshift. △ Less

Submitted 15 November, 2013; originally announced November 2013.

Comments: 20 pages, 16 figures. Accepted for publication in A&A

Variation in the dust emissivity index across M33 with Herschel and Spitzer (HerM33es)

Authors: F. S. Tabatabaei, J. Braine, E. M. Xilouris, C. Kramer, M. Boquien, F. Combes, C. Henkel, M. Relano, S. Verley, P. Gratier, F. Israel, M. C. Wiedner, M. Roellig, K. F. Schuster, P. van derWerf

Abstract: We study the wavelength dependence of the dust emission as a function of position and environment across the disk of M33 at a linear resolution of 160 pc using Spitzer and Herschel photometric data. Expressing the emissivity of the dust as a power law, the power-law exponent (beta) is estimated from two independent approaches designed to properly treat the degeneracy between beta and the dust temp… ▽ More We study the wavelength dependence of the dust emission as a function of position and environment across the disk of M33 at a linear resolution of 160 pc using Spitzer and Herschel photometric data. Expressing the emissivity of the dust as a power law, the power-law exponent (beta) is estimated from two independent approaches designed to properly treat the degeneracy between beta and the dust temperature. Both beta and the dust temperature are higher in the inner disk than in the outer disk, contrary to reported beta-T anti-correlations found in other sources. In the cold + warm dust model, the warm component and the ionized gas (Halpha) have a very similar distribution across the galaxy, demonstrating that the model separates the components in an appropriate fashion. The flocculent spiral arms and the dust lanes are evident in the map of the cold component. Both cold and warm dust column densities are high in star forming regions and reach their maxima toward the giant star forming complexes NGC604 and NGC595. beta declines from close to 2 in the center to about 1.3 in the outer disk. beta is positively correlated with star formation and with molecular gas column, as traced by Halpha and CO emission. The lower dust emissivity index in the outer parts of M33 is likely related to the reduced metallicity (different grain composition) and possibly different size distribution. It is not due to the decrease in stellar radiation field or temperature in a simple way because the FIR-bright regions in the outer disk also have a low beta. Like most spirals, M33 has a (decreasing) radial gradient in star formation and molecular-to-atomic gas ratio such that the regions bright in Halpha or CO tend to trace the inner disk, making it difficult to distinguish between their effects on the dust. △ Less

Submitted 15 October, 2013; originally announced October 2013.

Comments: 18 pages, 16 figures, Accepted for Publication in Astronomy & Astrophysics Journal

Gas and dust cooling along the major axis of M33 (HerM33es): ISO/LWS CII observations

Authors: C. Kramer, J. Abreu-Vicente, S. Garcia-Burillo, M. Relano, S. Aalto, M. Boquien, J. Braine, C. Buchbender, P. Gratier, F. P. Israel, T. Nikola, M. Roellig, S. Verley, P. van der Werf, E. M. Xilouris

Abstract: We aim to better understand the heating of the gas by observing the prominent gas cooling line [CII] at 158um in the low-metallicity environment of the Local Group spiral galaxy M33 at scales of 280pc. In particular, we aim at describing the variation of the photoelectric heating efficiency with galactic environment. In this unbiased study, we used ISO/LWS [CII] observations along the major axis o… ▽ More We aim to better understand the heating of the gas by observing the prominent gas cooling line [CII] at 158um in the low-metallicity environment of the Local Group spiral galaxy M33 at scales of 280pc. In particular, we aim at describing the variation of the photoelectric heating efficiency with galactic environment. In this unbiased study, we used ISO/LWS [CII] observations along the major axis of M33, in combination with Herschel PACS and SPIRE continuum maps, IRAM 30m CO 2-1 and VLA HI data to study the variation of velocity integrated intensities. The ratio of [CII] emission over the far-infrared continuum is used as a proxy for the heating efficiency, and models of photon-dominated regions are used to study the local physical densities, FUV radiation fields, and average column densities of the molecular clouds. The heating efficiency stays constant at 0.8% in the inner 4.5kpc radius of the galaxy where it starts to increase to reach values of ~3% in the outskirts at about 6kpc radial distance. The rise of efficiency is explained in the framework of PDR models by lowered volume densities and FUV fields, for optical extinctions of only a few magnitudes at constant metallicity. In view of the significant fraction of HI emission stemming from PDRs, and for typical pressures found in the Galactic cold neutral medium (CNM) traced by HI emission, the CNM contributes ~15% to the observed [CII] emission in the inner 2kpc radius of M33. The CNM contribution remains largely undetermined in the south, while positions between 2 and 7.3kpc radial distance in the north of M33 show a contribution of ~40%+-20%. △ Less

Submitted 31 March, 2013; originally announced April 2013.

Comments: (16 pages, a version with higher resolution figures may be downloaded from http://www.iram.es/IRAMES/hermesWiki/HermesPublications)

Spectral Energy Distributions of HII regions in M33 (HerM33es)

Authors: M. Relano, S. Verley, I. Perez, C. Kramer, D. Calzetti, E. M. Xilouris, M. Boquien, J. Abreu-Vicente, F. Combes, F. Israel, F. S. Tabatabaei, J. Braine, C. Buchbender, M. Gonzalez, P. Gratier, S. Lord, B. Mookerjea, G. Quintana-Lacaci, P. van der Werf

Abstract: Within the framework of the Herschel M 33 extended survey HerM33es we study the Spectral Energy Distribution (SED) of a set of HII regions in M 33 as a function of the morphology. We present a catalogue of 119 HII regions morphologically classified: 9 filled, 47 mixed, 36 shell, and 27 clear shell HII regions. For each object we extract the photometry at twelve available wavelength bands (from FUV… ▽ More Within the framework of the Herschel M 33 extended survey HerM33es we study the Spectral Energy Distribution (SED) of a set of HII regions in M 33 as a function of the morphology. We present a catalogue of 119 HII regions morphologically classified: 9 filled, 47 mixed, 36 shell, and 27 clear shell HII regions. For each object we extract the photometry at twelve available wavelength bands (from FUV-1516A to IR-250mi) and obtain the SED. We also obtain emission line profiles across the regions to study the location of the stellar, ionised gas, and dust components. We find trends for the SEDs related to the morphology, showing that the star and gas-dust configuration affects the ratios of the emission in different bands. The mixed and filled regions show higher emission at 24mi than the shells and clear shells, which could be due to the proximity of the dust to the stellar clusters in the case of mixed and filled regions. The FIR peak for shells and clear shells seems to be located towards longer wavelengths, indicating that the dust is colder for this type of objects.The logarithmic 100/70mi ratio for filled and mixed regions remains constant over one order of magnitude in Halpha and FUV surface brightness, while the shells and clear shells exhibit a wider range of values of almost two orders of magnitude. We derive dust masses and temperatures fitting the individual SEDs with dust models proposed in the literature. The derived dust mass range is between 10^2-10^4 Msun and the cold dust temperature spans T(cold)~12-27 K. The spherical geometrical model proposed for the Halpha clear shells is confirmed by the emission profile obtained from the observations and is used to infer the electron density within the envelope: the typical electron density is 0.7+-0.3 cm^-3, while filled regions can reach values two to five times higher. △ Less

Submitted 24 January, 2013; originally announced January 2013.

Comments: 21 pages, 19 figures. Accepted to A&A

The dust SED of dwarf galaxies I. The case of NGC 4214

Authors: Israel Hermelo, Ute Lisenfeld, Monica Relaño, Richard J. Tuffs, Cristina C. Popescu, Brent Groves

Abstract: The goal of the present study is to establish the physical origin of dust heating and emission based on radiation transfer models, which self-consistently connect the emission components from diffuse dust and the dust in massive star forming regions. NGC 4214 is a nearby dwarf galaxy with a large set of ancillary data, ranging from the ultraviolet (UV) to radio, including maps from SPITZER, HERSCH… ▽ More The goal of the present study is to establish the physical origin of dust heating and emission based on radiation transfer models, which self-consistently connect the emission components from diffuse dust and the dust in massive star forming regions. NGC 4214 is a nearby dwarf galaxy with a large set of ancillary data, ranging from the ultraviolet (UV) to radio, including maps from SPITZER, HERSCHEL and detections from PLANCK. We mapped this galaxy with MAMBO at 1.2 mm at the IRAM 30 m telescope. We extract separate dust emission components for the HII regions (plus their associated PDRs on pc scales) and for the diffuse dust (on kpc scales). We analyse the full UV to FIR/submm SED of the galaxy using a radiation transfer model which self-consistently treats the dust emission from diffuse and SF complexes components, considering the illumination of diffuse dust both by the distributed stellar populations, and by escaping light from the HII regions. While maintaining consistency with the framework of this model we additionally use a model that provides a detailed description of the dust emission from the HII regions and their surrounding PDRs on pc scales. Due to the large amount of available data and previous studies for NGC 4214 very few free parameters remained in the model fitting process. We achieve a satisfactory fit for the emission from HII+PDR regions on pc scales, with the exception of the emission at 8\mi, which is underpredicted by the model. For the diffuse emission we achieve a good fit if we assume that about 30-70% of the emission escaping the HII+PDR regions is able to leave the galaxy without passing through a diffuse ISM, which is not an unlikely scenario for a dwarf galaxy which has recently undergone a nuclear starburst. We determine a dust-to-gas mass ratio of 350-390 which is close to the expected value based on the metallicity. △ Less

Submitted 26 October, 2012; originally announced October 2012.

Comments: 19 pages, 7 figures, accepted by A&A

Cool and warm dust emission from M33 (HerM33es)

Authors: E. M. Xilouris, F. S. Tabatabaei, M. Boquien, C. Kramer, C. Buchbender, F. Bertoldi, S. Anderl, J. Braine, S. Verley, M. Relano, G. Quintana-Lacaci, S. Akras, R. Beck, D. Calzetti, F. Combes, M. Gonzalez, P. Gratier, C. Henkel, F. Israel, B. Koribalski, S. Lord, B. Mookerjea, E. Rosolowsky, G. Stacey, R. P. J. Tilanus , et al. (2 additional authors not shown)

Abstract: We study the far-infrared emission from the nearby spiral galaxy M33 in order to investigate the dust physical properties such as the temperature and the luminosity density across the galaxy. Taking advantage of the unique wavelength coverage (100, 160, 250, 350 and 500 micron) of the Herschel Space Observatory and complementing our dataset with Spitzer-IRAC 5.8 and 8 micron and Spitzer-MIPS 24 an… ▽ More We study the far-infrared emission from the nearby spiral galaxy M33 in order to investigate the dust physical properties such as the temperature and the luminosity density across the galaxy. Taking advantage of the unique wavelength coverage (100, 160, 250, 350 and 500 micron) of the Herschel Space Observatory and complementing our dataset with Spitzer-IRAC 5.8 and 8 micron and Spitzer-MIPS 24 and 70 micron data, we construct temperature and luminosity density maps by fitting two modified blackbodies of a fixed emissivity index of 1.5. We find that the 'cool' dust grains are heated at temperatures between 11 and 28 K with the lowest temperatures found in the outskirts of the galaxy and the highest ones in the center and in the bright HII regions. The infrared/submillimeter total luminosity (5 - 1000 micron) is estimated to be 1.9x10^9 Lsun. 59% of the total luminosity of the galaxy is produced by the 'cool' dust grains (~15 K) while the rest 41% is produced by 'warm' dust grains (~55 K). The ratio of the cool-to-warm dust luminosity is close to unity (within the computed uncertainties), throughout the galaxy, with the luminosity of the cool dust being slightly enhanced in the center of the galaxy. Decomposing the emission of the dust into two components (one emitted by the diffuse disk of the galaxy and one emitted by the spiral arms) we find that the fraction of the emission in the disk in the mid-infrared (24 micron) is 21%, while it gradually rises up to 57% in the submillimeter (500 micron). We find that the bulk of the luminosity comes from the spiral arm network that produces 70% of the total luminosity of the galaxy with the rest coming from the diffuse dust disk. The 'cool' dust inside the disk is heated at a narrow range of temperatures between 18 and 15 K (going from the center to the outer parts of the galaxy). △ Less

Submitted 7 May, 2012; originally announced May 2012.

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

On how leakage can affect the Star Formation Rate estimation using Halpha luminosity

Authors: M. Relano, R. C. Jr. Kennicutt, J. J. Eldridge, J. C. Lee, S. Verley

Abstract: We present observational evidence that leakage of ionising photons from star-forming regions can affect the quantification of the star formation rate (SFR) in galaxies. This effect could partially explain the differences between the SFR estimates using the far ultraviolet (FUV) and the Halpha emission. We find that leakage could decrease the SFR(Ha)/SFR(FUV) ratio by up to a 25 per cent. The evide… ▽ More We present observational evidence that leakage of ionising photons from star-forming regions can affect the quantification of the star formation rate (SFR) in galaxies. This effect could partially explain the differences between the SFR estimates using the far ultraviolet (FUV) and the Halpha emission. We find that leakage could decrease the SFR(Ha)/SFR(FUV) ratio by up to a 25 per cent. The evidence is based on the observation that the SFR(Ha)/SFR(FUV) ratio is lower for objects showing a shell Halpha structure than for regions exhibiting a much more compact morphology. The study has been performed on three object samples: low luminosity dwarf galaxies from the Local Volume Legacy survey and star-forming regions in the Large Magellanic Cloud and the nearby Local Group galaxy M33. For the three samples we find differences (1.1-1.4sigma) between the SFR(Ha)/SFR(FUV) for compact and shell objects. Although leakage cannot entirely explain the observed trend of SFR(Ha)/SFR(FUV) ratios for systems with low SFR, we show the mechanism can lead to different SFR estimates when using Halpha and FUV luminosities. Therefore, further study is needed to constrain the contribution of leakage to the low SFR(Ha)/SFR(FUV) ratios observed in dwarf galaxies and its impact on the Halpha flux as a SFR indicator in such objects. △ Less

Submitted 19 April, 2012; originally announced April 2012.

Comments: 8 pages, 7 figures; accepted for publication in MNRAS

Dust and gas power-spectrum in M33 (HERM33ES)

Authors: F. Combes, M. Boquien, C. Kramer, E. M. Xilouris, F. Bertoldi, J. Braine, C. Buchbender, D. Calzetti, P. Gratier, F. Israel, B. Koribalski, S. Lord, G. Quintana-Lacaci, M. Relano, M. Roellig, G. Stacey, F. S. Tabatabaei, R. P. J. Tilanus, F. van der Tak, P. van der Werf, S. Verley

Abstract: Power spectra of de-projected images of late-type galaxies in gas and/or dust emission are very useful diagnostics of the dynamics and stability of their interstellar medium. Previous studies have shown that the power spectra can be approximated as two power-laws, a shallow one at large scales (larger than 500 pc) and a steeper one at small scales, with the break between the two corresponding to t… ▽ More Power spectra of de-projected images of late-type galaxies in gas and/or dust emission are very useful diagnostics of the dynamics and stability of their interstellar medium. Previous studies have shown that the power spectra can be approximated as two power-laws, a shallow one at large scales (larger than 500 pc) and a steeper one at small scales, with the break between the two corresponding to the line-of-sight thickness of the galaxy disk. We present a thorough analysis of the power spectra of the dust and gas emission at several wavelengths in the nearby galaxy M33. In particular, we use the recently obtained images at five wavelengths by PACS and SPIRE onboard Herschel. The large dynamical range (2-3 dex in scale) of most images allow us to determine clearly the change in slopes from -1.5 to -4, with some variations with wavelength. The break scale is increasing with wavelength, from 100 pc at 24 and 100micron to 350 pc at 500micron, suggesting that the cool dust lies in a thicker disk than the warm dust, may be due to star formation more confined to the plane. The slope at small scale tends to be steeper at longer wavelength, meaning that the warmer dust is more concentrated in clumps. Numerical simulations of an isolated late-type galaxy, rich in gas and with no bulge, like M33, are carried out, in order to better interpret these observed results. Varying the star formation and feedback parameters, it is possible to obtain a range of power-spectra, with two power-law slopes and breaks, which nicely bracket the data. The small-scale power-law is indeed reflecting the 3D behaviour of the gas layer, steepening strongly while the feedback smoothes the structures, by increasing the gas turbulence. M33 appears to correspond to a fiducial model with an SFR of $\sim$ 0.7 Mo/yr, with 10% supernovae energy coupled to the gas kinematics. △ Less

Submitted 12 January, 2012; originally announced January 2012.

Comments: 11 pages, 24 figures, accepted in Astronomy & Astrophysics

Modeling the dust Spectral Energy Distribution of NGC 4214

Authors: I. Hermelo, U. Lisenfeld, M. Relano, R. Tuffs, J. Fischera, B. Groves, C. Popescu

Abstract: We have carried out a detailed modeling of the dust Spectral Energy Distribution (SED) of the nearby, starbursting dwarf galaxy NGC 4214. A key point of our modeling is that we distinguish the emission from (i) HII regions and their associated photodissociation regions (PDRs) and (ii) diffuse dust. For both components we apply templates from the literature calculated with a realistic geometry and… ▽ More We have carried out a detailed modeling of the dust Spectral Energy Distribution (SED) of the nearby, starbursting dwarf galaxy NGC 4214. A key point of our modeling is that we distinguish the emission from (i) HII regions and their associated photodissociation regions (PDRs) and (ii) diffuse dust. For both components we apply templates from the literature calculated with a realistic geometry and including radiation transfer. The large amount of existing data from the ultraviolet (UV) to the radio allows the direct measurement of most of the input parameters of the models. We achieve a good fit for the total dust SED of NGC 4214. In the present contribution we describe the available data, the data reduction and the determination of the model parameters, whereas a description of the general outline of our work is presented in the contribution of Lisenfeld et al. in this volume. △ Less

Submitted 24 November, 2011; originally announced November 2011.

Comments: 3 pages, 1 figure. To appear in 'The Spectral Energy Distribution of Galaxies' Proceedings IAU Symposium No 284, 2011

Dust in dwarf galaxies: The case of NGC 4214

Authors: U. Lisenfeld, I. Hermelo, M. Relano, R. Tuffs, J. Fischera, B. Groves, C. Popescu

Abstract: We have carried out a detailed modelling of the dust heating and emission in the nearby, starbursting dwarf galaxy NGC 4214. Due to its proximity and the great wealth of data from the UV to the millimeter range (from GALEX, HST, {\it Spitzer}, Herschel, Planck and IRAM) it is possible to separately model the emission from HII regions and their associated photodissociation regions (PDRs) and the em… ▽ More We have carried out a detailed modelling of the dust heating and emission in the nearby, starbursting dwarf galaxy NGC 4214. Due to its proximity and the great wealth of data from the UV to the millimeter range (from GALEX, HST, {\it Spitzer}, Herschel, Planck and IRAM) it is possible to separately model the emission from HII regions and their associated photodissociation regions (PDRs) and the emission from diffuse dust. Furthermore, most model parameters can be directly determined from the data leaving very few free parameters. We can fit both the emission from HII+PDR regions and the diffuse emission in NGC 4214 with these models with "normal" dust properties and realistic parameters. △ Less

Submitted 24 November, 2011; originally announced November 2011.

Comments: 4pages, 3 figures. To appear in 'The Spectral Energy Distribution of Galaxies' Proceedings IAU Symposium No 284, 2011

Spectral Energy Distributions of a set of HII regions in M33 (HerM33es)

Authors: M. Relano, S. Verley, I. Perez, C. Kramer, E. M. Xilouris, M. Boquien, J. Braine, D. Calzetti, C. Henkel, HerM33es Team

Abstract: Within the framework of the HerM33es Key Project for Herschel and in combination with multi-wavelength data, we study the Spectral Energy Distribution (SED) of a set of HII regions in the Local Group Galaxy M33. Using the Halpha emission, we perform a classification of a selected HII region sample in terms of morphology, separating the objects in filled, mixed, shell and clear shell objects. We ob… ▽ More Within the framework of the HerM33es Key Project for Herschel and in combination with multi-wavelength data, we study the Spectral Energy Distribution (SED) of a set of HII regions in the Local Group Galaxy M33. Using the Halpha emission, we perform a classification of a selected HII region sample in terms of morphology, separating the objects in filled, mixed, shell and clear shell objects. We obtain the SED for each HII region as well as a representative SED for each class of objects. We also study the emission distribution of each band within the regions. We find different trends in the SEDs for each morphological type that are related to properties of the dust and their associated stellar cluster. The emission distribution of each band within the region is different for each morphological type of object. △ Less

Submitted 23 November, 2011; originally announced November 2011.

Comments: 3pages, 4 figures. To appear in 'The Spectral Energy Distribution of Galaxies' Proceedings IAU Symposium No 284, 2011

CALIFA, the Calar Alto Legacy Integral Field Area survey: I. Survey presentation

Authors: S. F. Sanchez, R. C. Kennicutt, A. Gil de Paz, G. van de Ven, J. M. Vílchez, L. Wisotzki, C. J. Walcher, D. Mast, J. A. L. Aguerri, S. Albiol-Perez, A. Alonso-Herrero, J. Alves, J. Bakos, T. Bartakova, J. Bland-Hawthorn, A. Boselli, D. J. Bomans, A. Castillo-Morales, C. Cortijo-Ferrero, A. de Lorenzo-Caceres, A. del Olmo, R. -J. Dettmar, A. Díaz, S. Ellis, J. Falcon-Barroso , et al. (47 additional authors not shown)

Abstract: We present here the Calar Alto Legacy Integral Field Area (CALIFA) survey, which has been designed to provide a first step in this direction.We summarize the survey goals and design, including sample selection and observational strategy.We also showcase the data taken during the first observing runs (June/July 2010) and outline the reduction pipeline, quality control schemes and general characteri… ▽ More We present here the Calar Alto Legacy Integral Field Area (CALIFA) survey, which has been designed to provide a first step in this direction.We summarize the survey goals and design, including sample selection and observational strategy.We also showcase the data taken during the first observing runs (June/July 2010) and outline the reduction pipeline, quality control schemes and general characteristics of the reduced data. This survey is obtaining spatially resolved spectroscopic information of a diameter selected sample of $\sim600$ galaxies in the Local Universe (0.005< z <0.03). CALIFA has been designed to allow the building of two-dimensional maps of the following quantities: (a) stellar populations: ages and metallicities; (b) ionized gas: distribution, excitation mechanism and chemical abundances; and (c) kinematic properties: both from stellar and ionized gas components. CALIFA uses the PPAK Integral Field Unit (IFU), with a hexagonal field-of-view of $\sim1.3\sq\arcmin'$, with a 100% covering factor by adopting a three-pointing dithering scheme. The optical wavelength range is covered from 3700 to 7000 Å, using two overlapping setups (V500 and V1200), with different resolutions: R\sim850 and R\sim1650, respectively. CALIFA is a legacy survey, intended for the community. The reduced data will be released, once the quality has been guaranteed. The analyzed data fulfill the expectations of the original observing proposal, on the basis of a set of quality checks and exploratory analysis. We conclude from this first look at the data that CALIFA will be an important resource for archaeological studies of galaxies in the Local Universe. △ Less

Submitted 4 November, 2011; v1 submitted 3 November, 2011; originally announced November 2011.

Comments: 32 pages, 29 figures, Accepted for publishing in Astronomy and Astrophysics

Detection of infalling hydrogen in transfer between the interacting galaxies NGC 5426 and NGC 5427

Authors: Joan Font, John E. Beckman, Margarita Rosado, Benoît Epinat, Kambiz Fathi, Olivier Hernandez, Claude Carignan, Leonel Gutiérrez, Monica Relaño, Javier Blasco-Herrera, Isaura Fuentes-Carrera

Abstract: Using velocity tagging we have detected hydrogen from NGC 5426 falling onto its interacting partner NGC 5427. Our observations, with the GHaFaS Fabry-Perot spectrometer, produced maps of the two galaxies in Halpha surface brightness and radial velocity. We found emission with the range of velocities associated with NGC 5426 along lines of sight apparently emanating from NGC 5427, superposed on the… ▽ More Using velocity tagging we have detected hydrogen from NGC 5426 falling onto its interacting partner NGC 5427. Our observations, with the GHaFaS Fabry-Perot spectrometer, produced maps of the two galaxies in Halpha surface brightness and radial velocity. We found emission with the range of velocities associated with NGC 5426 along lines of sight apparently emanating from NGC 5427, superposed on the velocity map of the latter. After excluding instrumental effects we assign the anomalous emission to gas pulled from NGC 5426 during its passage close to NGC 5427. Its distribution, more intense between the arms and just outside the disk of NGC 5427, and weak, or absent, in the arms, suggests that the infalling gas is behind the disk., ionized by Lyman continuum photons escaping from NGC 5427. Modeling this, we estimate the distances of these gas clouds- behind the plane: a few hundred pc to a few kpc. We also estimate the mass of the infalling (ionized plus neutral) gas, finding an infall rate of 10 solar masses per year, consistent with the high measured SFR across the disk of NGC 5427 and with the detected circumnuclear galactic wind. △ Less

Submitted 1 August, 2011; originally announced August 2011.

Comments: 10 pages, 4 figures. Accepted for publication in ApJL

Dust heating sources in galaxies: the case of M33 (HERM33ES)

Authors: M. Boquien, D. Calzetti, F. Combes, C. Henkel, F. Israel, C. Kramer, M. Relaño, S. Verley, P. van der Werf, E. M. Xilouris

Abstract: Dust emission is one of the main windows to the physics of galaxies and to star formation as the radiation from young, hot stars is absorbed by the dust and reemitted at longer wavelengths. The recently launched Herschel satellite now provides a view of dust emission in the far-infrared at an unequaled resolution and quality up to 500 μm. In the context of the Herschel HERM33ES open time key proje… ▽ More Dust emission is one of the main windows to the physics of galaxies and to star formation as the radiation from young, hot stars is absorbed by the dust and reemitted at longer wavelengths. The recently launched Herschel satellite now provides a view of dust emission in the far-infrared at an unequaled resolution and quality up to 500 μm. In the context of the Herschel HERM33ES open time key project, we are studying the moderately inclined Scd local group galaxy M33 which is located only 840 kpc away. In this article, using Spitzer and Herschel data ranging from 3.6 μm to 500 μm, along with HI, Hα maps, and GALEX ultraviolet data we have studied the emission of the dust at the high spatial resolution of 150 pc. Combining Spitzer and Herschel bands, we have provided new, inclination corrected, resolved estimators of the total infrared brightness and of the star formation rate from any combination of these bands. The study of the colors of the warm and cold dust populations shows that the temperature of the former is, at high brightness, dictated by young massive stars but, at lower brightness, heating is taken over by the evolved populations. Conversely, the temperature of the cold dust is tightly driven by the evolved stellar populations. △ Less

Submitted 12 July, 2011; originally announced July 2011.

Comments: 29 pages, 12 figures, accepted for publication in AJ

On the equivalent effective temperature of massive young star clusters: The case of NGC 595

Authors: Enrique Pérez-Montero, José M. Vílchez, Mónica Relaño, Ana Monreal-Ibero

Abstract: The softness parameter is based on the relative intensity of several optical emission lines emitted by the gas ionized by young massive star clusters and can be used to derive the equivalent effective temperature (T*) in those objects whose stellar population cannot be resolved. This method has several uncertainties due to the disagreement between different synthesis model atmospheres but it is ro… ▽ More The softness parameter is based on the relative intensity of several optical emission lines emitted by the gas ionized by young massive star clusters and can be used to derive the equivalent effective temperature (T*) in those objects whose stellar population cannot be resolved. This method has several uncertainties due to the disagreement between different synthesis model atmospheres but it is robust to study the relative variations between objects. Following the 2D photoionization models of the giant HII region NGC 595 (Pérez-Montero et al. 2011) we show that the determination of T* with the softness parameter is also robust in different regions of a same object with large variations in the geometry of the gas and in the dust-to-gas ratio. △ Less

Submitted 4 July, 2011; originally announced July 2011.

Comments: Proceeding for the Workshop "Star clusters and Associations: A RIA Workshop on Gaia", held in Granada from 23rd to 27th May, 2011. 6 pages and 4 figures

The Herschel M33 extended survey (HerM33es): PACS spectroscopy of the star forming region BCLMP 302

Authors: B. Mookerjea, C. Kramer, C. Buchbender, M. Boquien, S. Verley, M. Relano, G. Quintana-Lacaci, S. Aalto, J. Braine, D. Calzetti, F. Combes, S. Garcia-Burillo, P. Gratier, C. Henkel, F. Israel, S. Lord, T. Nikola, M. Roellig, G. Stacey, F. S. Tabatabaei, F. van der Tak, P. van der Werf

Abstract: Context: The emission line of [CII] at 158 micron is one of the strongest cooling lines of the interstellar medium (ISM) in galaxies. Aims: Disentangling the relative contributions of the different ISM phases to [CII] emission, is a major topic of the HerM33es program, a Herschel key project to study the ISM in the nearby spiral galaxy M33. Methods: Using PACS, we have mapped the emission of [CII]… ▽ More Context: The emission line of [CII] at 158 micron is one of the strongest cooling lines of the interstellar medium (ISM) in galaxies. Aims: Disentangling the relative contributions of the different ISM phases to [CII] emission, is a major topic of the HerM33es program, a Herschel key project to study the ISM in the nearby spiral galaxy M33. Methods: Using PACS, we have mapped the emission of [CII] 158 micron, [OI] 63 micron, and other FIR lines in a 2'x2' region of the northern spiral arm of M33, centered on the HII region BCLMP302. At the peak of H-alpha emission, we have observed in addition a velocity resolved [CII] spectrum using HIFI. We use scatterplots to compare these data with PACS 160 micron continuum maps, and with maps of CO and HI data, at a common resolution of 12 arcsec or 50 pc. Maps of H-alpha and 24 micron emission observed with Spitzer are used to estimate the SFR. We have created maps of the [CII] and [OI] 63 micron emission and detected [NII] 122 micron and NIII 57 micron at individual positions. Results: The [CII] line observed with HIFI is significantly broader than that of CO, and slightly blue-shifted. In addition, there is little spatial correlation between [CII] observed with PACS and CO over the mapped region. There is even less spatial correlation between [CII] and the atomic gas traced by HI. Detailed comparison of the observed intensities towards the HII region with models of photo ionization and photon dominated regions, confirms that a significant fraction, 20--30%, of the observed [CII] emission stems from the ionized gas and not from the molecular cloud. The gas heating efficiency, using the ratio between [CII] and the TIR as a proxy, varies between 0.07 and 1.5%, with the largest variations found outside the HII region. △ Less

Submitted 27 June, 2011; originally announced June 2011.

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

Star formation in M33 (HerM33es)

Authors: C. Kramer, M. Boquien, J. Braine, C. Buchbender, D. Calzetti, P. Gratier, B. Mookerjea, M. Relano, S. Verley

Abstract: Within the key project "Herschel M33 extended survey" (HerM33es), we are studying the physical and chemical processes driving star formation and galactic evolution in the nearby galaxy M33, combining the study of local conditions affecting individual star formation with properties only becoming apparent on global scales. Here, we present recent results obtained by the HerM33es team. Combining Spit… ▽ More Within the key project "Herschel M33 extended survey" (HerM33es), we are studying the physical and chemical processes driving star formation and galactic evolution in the nearby galaxy M33, combining the study of local conditions affecting individual star formation with properties only becoming apparent on global scales. Here, we present recent results obtained by the HerM33es team. Combining Spitzer and Herschel data ranging from 3.6um to 500um, along with HI, Halpha, and GALEX UV data, we have studied the dust at high spatial resolutions of 150pc, providing estimators of the total infrared (TIR) brightness and of the star formation rate. While the temperature of the warm dust at high brightness is driven by young massive stars, evolved stellar populations appear to drive the temperature of the cold dust. Plane-parallel models of photon dominated regions (PDRs) fail to reproduce fully the [CII], [OI], and CO maps obtained in a first spectroscopic study of one 2'x2' subregion of M33, located on the inner, northern spiral arm and encompassing the HII region BCLMP302. △ Less

Submitted 26 May, 2011; originally announced May 2011.

Comments: 6 pages, to appear in the proceedings of the 5th Zermatt ISM Symposium "Conditions and impact of star formation: New results with Herschel and beyond"

A 2D multiwavelength study of the ionized gas and stellar population in the Giant HII Region NGC 588

Authors: A. Monreal-Ibero, M. Relaño, C. Kehrig, E. Pérez-Montero, J. M. Vílchez, A. Kelz, M. M. Roth, O. Streicher

Abstract: (ABRIDGED) We present an analysis of NGC588 based on IFS data with PMAS, together with Spitzer images at 8 mi and 24 mi. The extinction distribution in the optical shows complex structure, with maxima correlating in position with those of the emission at 24 mi and 8 mi. The Ha luminosity absorbed by the dust within the GHIIR reproduces the structure observed in the 24 mi image, supporting the use… ▽ More (ABRIDGED) We present an analysis of NGC588 based on IFS data with PMAS, together with Spitzer images at 8 mi and 24 mi. The extinction distribution in the optical shows complex structure, with maxima correlating in position with those of the emission at 24 mi and 8 mi. The Ha luminosity absorbed by the dust within the GHIIR reproduces the structure observed in the 24 mi image, supporting the use of this band as a tracer of recent star formation. A velocity difference of ~50 km/s was measured between the areas of high and low surface brightness, which would be expected if NGC588 were an evolved GHIIR. Line ratios used in the BPT diagnostic diagrams show a larger range of variation in the low surface brightness areas. The ranges are ~0.5 to 1.2 dex for [NII]/Ha, 0.7 to 1.7 dex for [SII]/Ha, and 0.3 to 0.5 dex for [OIII]/Hb. Ratios corresponding to large ionization parameter (U) are found between the peak of the emission in Hb and the main ionizing source decreasing radially outwards within the region. Differences between the integrated and local values of the U tracers can be as high as ~0.8 dex. [OII]/Hb and [OIII]/[OII] yield similar local values for U and consistent with those expected from the integrated spectrum of an HII region ionized by a single star. The ratio [SII]/Ha departs significantly from the range predicted by this scenario, indicating the complex ionization structure in GHIIRs. There is a significant scatter in derivations of Z using strong line tracers as a function of position, caused by variations in the degree of ionization. The scatter is smaller for N2O3 which points to this tracer as a better Z tracer than N2. The comparison between integrated and local line ratio values indicates that measurements of the line ratios of GHIIR in galaxies at distances >~25 Mpc may be dominated by the ionization conditions in their low surface brightness areas. △ Less

Submitted 10 January, 2011; originally announced January 2011.

Comments: accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal, 2010 January 5; 18 pages; 15 figures

A photoionization model of the spatial distribution of the optical and mid-IR properties in NGC595

Authors: Enrique Perez-Montero, Monica Relano, Jose M. Vilchez, Ana Monreal-Ibero

Abstract: We present a set of photoionization models that reproduce simultaneously the observed optical and mid-infrared spatial distribution of the HII region NGC595 in the disk of M33 using the code CLOUDY. Both optical (PMAS-Integral Field Spectroscopy) and mid-infrared (8 mi and 24 mi bands from Spitzer) data provide enough spatial resolution to model in a novel approach the inner structure of the HII r… ▽ More We present a set of photoionization models that reproduce simultaneously the observed optical and mid-infrared spatial distribution of the HII region NGC595 in the disk of M33 using the code CLOUDY. Both optical (PMAS-Integral Field Spectroscopy) and mid-infrared (8 mi and 24 mi bands from Spitzer) data provide enough spatial resolution to model in a novel approach the inner structure of the HII region. We define a set of elliptical annular regions around the central ionizing cluster with an uniformity in their observed properties and consider each annulus as an independent thin shell structure. For the first time our models fit the relative surface brightness profiles in both the optical (Halpha, [OII], [OIII]) and the mid-infrared emissions (8 mi and 24 mi), under the assumption of a uniform metallicity (12+log(O/H) = 8.45; Esteban et al. 2009) and an age for the stellar cluster of 4.5 Myr (Malumuth et al. 1996). Our models also reproduce the observed uniformity of the R23 parameter and the increase of the [OII]/[OIII] ratio due to the decrease of the ionization parameter. The variation of the Halpha profile is explained in terms of the differences of the occupied volume (the product of filling factor and total volume of the shell) in a matter-bounded geometry, which also allows to reproduce the observed pattern of the extinction. The 8 mi/24 mi ratio is low (ranging between 0.04 and 0.4) because it is dominated by the surviving of small dust grains in the HII region, while the PAHs emit more weakly because they cannot be formed in these thin HII gas shells. The ratio is also well fitted in our models by assuming a dust-to-gas ratio in each annulus compatible with the integrated estimate for the whole HII region after the 70 mi, and 160 mi Spitzer observations. △ Less

Submitted 28 October, 2010; originally announced October 2010.

Comments: Accepted for publication in MNRAS, 9 pages, 17 figures

The red supergiants & Wolf-Rayet stars of NGC 604

Authors: John J. Eldridge, Monica Relano

Abstract: We study the post-main sequence stars in NGC 604, the most luminous HII region in M33. Previously, a number of Wolf-Rayet (WR) stars and one red supergiant (RSG) have been discovered. Based on broadband photometry of the region, we present evidence that is consistent with the presence of this RSG and with three more RSG candidates. Using SED fitting based on HST UVIJHK photometry we estimate the a… ▽ More We study the post-main sequence stars in NGC 604, the most luminous HII region in M33. Previously, a number of Wolf-Rayet (WR) stars and one red supergiant (RSG) have been discovered. Based on broadband photometry of the region, we present evidence that is consistent with the presence of this RSG and with three more RSG candidates. Using SED fitting based on HST UVIJHK photometry we estimate the ages of the WR stars and RSGs finding that the two populations are from distinct formation episodes with ages 3.2$\pm%1.0Myrs and 12.4$\pm$2.1Myrs, respectively. The RSGs have greater extinctions towards their line of sight than the WR stars consistent with the RSGs producing large amount of dust. Using the WR and RSG populations and similar SED fits to the most massive O stars we estimate that the total stellar mass is (3.8 $\pm$ 0.6) x 10^5Msun. We find a large discrepancy between the expected Hα flux from such a massive cluster and that one observed. This suggests that 49 (+16,-19) percent of the ionizing photons produced by massive stars in NGC 604 is leaking from this HII region. We also suggest that the implications of an old RSG population mean that if NGC 604 was more distant and only observed in the infrared (IR) it would be difficult to study the youngest burst of star formation due to the contamination of RSGs. △ Less

Submitted 9 September, 2010; originally announced September 2010.

Comments: Accepted for publication in MNRAS, 13 pages, 10 figures, 3 tables

Cool gas and dust in M33: Results from the Herschel M33 extended survey (HERM33ES)

Authors: J. Braine, P. Gratier, C. Kramer, E. M. Xilouris, E. Rosolowsky, C. Buchbender, M. Boquien, D. Calzetti, G. Quintana-Lacaci, F. Tabatabaei, S. Verley, F. Israel, F. van der Tak, S. Aalto, F. Combes, S. Garcia-Burillo, M. Gonzalez, C. Henkel, B. Koribalski, B. Mookerjea, M. Roellig, K. F. Schuster, M. Relano, F. Bertoldi, P. van der Werf , et al. (1 additional authors not shown)

Abstract: We present an analysis of the first space-based far-IR-submm observations of M 33, which measure the emission from the cool dust and resolve the giant molecular cloud complexes. With roughly half-solar abundances, M33 is a first step towards young low-metallicity galaxies where the submm may be able to provide an alternative to CO mapping to measure their H$_2$ content. In this Letter, we measure… ▽ More We present an analysis of the first space-based far-IR-submm observations of M 33, which measure the emission from the cool dust and resolve the giant molecular cloud complexes. With roughly half-solar abundances, M33 is a first step towards young low-metallicity galaxies where the submm may be able to provide an alternative to CO mapping to measure their H$_2$ content. In this Letter, we measure the dust emission cross-section $σ$ using SPIRE and recent CO and \HI\ observations; a variation in $σ$ is present from a near-solar neighborhood cross-section to about half-solar with the maximum being south of the nucleus. Calculating the total H column density from the measured dust temperature and cross-section, and then subtracting the \HI\ column, yields a morphology similar to that observed in CO. The H$_2$/\HI\ mass ratio decreases from about unity to well below 10% and is about 15% averaged over the optical disk. The single most important observation to reduce the potentially large systematic errors is to complete the CO mapping of M 33. △ Less

Submitted 19 May, 2010; originally announced May 2010.

Comments: 5 pages, 5 figures Accepted for publication in Astronomy and Astrophysics

PACS and SPIRE photometer maps of M33: First results of the Herschel M33 extended survey (HERM33ES)

Authors: C. Kramer, C. Buchbender, E. M. Xilouris, M. Boquien, J. Braine, D. Calzetti, S. Lord, B. Mookerjea, G. Quintana-Lacaci, M. Relano, G. Stacey, F. S. Tabatabaei, S. Verley, S. Aalto, S. Akras, M. Albrecht, S. Anderl, R. Beck, F. Bertoldi, F. Combes, M. Dumke, S. Garcia-Burillo, M. Gonzalez, P. Gratier, R. Gueusten , et al. (15 additional authors not shown)

Abstract: Within the framework of the HERM33ES key project, we are studying the star forming interstellar medium in the nearby, metal-poor spiral galaxy M33, exploiting the high resolution and sensitivity of Herschel. We use PACS and SPIRE maps at 100, 160, 250, 350, and 500 micron wavelength, to study the variation of the spectral energy distributions (SEDs) with galacto-centric distance. Detailed SED mode… ▽ More Within the framework of the HERM33ES key project, we are studying the star forming interstellar medium in the nearby, metal-poor spiral galaxy M33, exploiting the high resolution and sensitivity of Herschel. We use PACS and SPIRE maps at 100, 160, 250, 350, and 500 micron wavelength, to study the variation of the spectral energy distributions (SEDs) with galacto-centric distance. Detailed SED modeling is performed using azimuthally averaged fluxes in elliptical rings of 2 kpc width, out to 8 kpc galacto-centric distance. Simple isothermal and two-component grey body models, with fixed dust emissivity index, are fitted to the SEDs between 24 and 500 micron using also MIPS/Spitzer data, to derive first estimates of the dust physical conditions. The far-infrared and submillimeter maps reveal the branched, knotted spiral structure of M33. An underlying diffuse disk is seen in all SPIRE maps (250-500 micron). Two component fits to the SEDs agree better than isothermal models with the observed, total and radially averaged flux densities. The two component model, with beta fixed at 1.5, best fits the global and the radial SEDs. The cold dust component clearly dominates; the relative mass of the warm component is less than 0.3% for all the fits. The temperature of the warm component is not well constrained and is found to be about 60K plus/minus 10K. The temperature of the cold component drops significantly from about 24K in the inner 2 kpc radius to 13K beyond 6 kpc radial distance, for the best fitting model. The gas-to-dust ratio for beta=1.5, averaged over the galaxy, is higher than the solar value by a factor of 1.5 and is roughly in agreement with the subsolar metallicity of M33. △ Less

Submitted 14 May, 2010; originally announced May 2010.

Comments: 5 pages, 3 figures, accepted for publication in the A&A Herschel Special Issue

Properties of compact 250 μm emission and HII regions in M33 (HERM33ES)

Authors: S. Verley, M. Relaño, C. Kramer, E. M. Xilouris, M. Boquien, D. Calzetti, F. Combes, C. Buchbender, J. Braine, G. Quintana-Lacaci, F. S. Tabatabaei, S. Lord, F. Israel, G. Stacey, P. van der Werf

Abstract: Within the framework of the HERM33ES Key Project, using the high resolution and sensitivity of the Herschel photometric data, we study the compact emission in the Local Group spiral galaxy M33 to investigate the nature of the compact SPIRE emission sources. We extracted a catalogue of sources at 250um in order to investigate the nature of this compact emission. Taking advantage of the unprecedente… ▽ More Within the framework of the HERM33ES Key Project, using the high resolution and sensitivity of the Herschel photometric data, we study the compact emission in the Local Group spiral galaxy M33 to investigate the nature of the compact SPIRE emission sources. We extracted a catalogue of sources at 250um in order to investigate the nature of this compact emission. Taking advantage of the unprecedented Herschel resolution at these wavelengths, we also focus on a more precise study of some striking Halpha shells in the northern part of the galaxy. We present a catalogue of 159 compact emission sources in M33 identified by SExtractor in the 250um SPIRE band that is the one that provides the best spatial resolution. We also measured fluxes at 24um and Halpha for those 159 extracted sources. The morphological study of the shells also benefits from a multiwavelength approach including Halpha, far-UV from GALEX, and infrared from both Spitzer IRAC 8um and MIPS 24um in order to make comparisons. For the 159 compact sources selected at 250um, we find a very strong Pearson correlation coefficient with the MIPS 24um emission (r24 = 0.94) and a rather strong correlation with the Halpha emission, although with more scatter (rHa = 0.83). The morphological study of the Halpha shells shows a displacement between far-ultraviolet, Halpha, and the SPIRE bands. The cool dust emission from SPIRE clearly delineates the Halpha shell structures. The very strong link between the 250um compact emission and the 24um and Halpha emissions, by recovering the star formation rate from standard recipes for HII regions, allows us to provide star formation rate calibrations based on the 250um compact emission alone. The different locations of the Halpha and far-ultraviolet emissions with respect to the SPIRE cool dust emission leads to a dynamical age of a few Myr for the Halpha shells and the associated cool dust. △ Less

Submitted 13 May, 2010; originally announced May 2010.

Comments: 4 pages, 3 figures, Accpeted for publication in the A&A Herschel Special Issue

100 mum and 160 mum emission as resolved star-formation rate estimators in M33 (HERM33ES)

Authors: M. Boquien, D. Calzetti, C. Kramer, E. M. Xilouris, F. Bertoldi, J. Braine, C. Buchbender, F. Combes, F. Israel, B. Koribalski, S. Lord, G. Quintana-Lacaci, M. Relaño, M. Röllig, G. Stacey, F. S. Tabatabaei, R. P. J. Tilanus, F. van der Tak, P. van der Werf, S. Verley

Abstract: Over the past few years several studies have provided estimates of the SFR (star-formation rate) or the total infrared luminosity from just one infrared band. However these relations are generally derived for entire galaxies, which are known to contain a large scale diffuse emission that is not necessarily related to the latest star-formation episode. We provide new relations to estimate the SFR f… ▽ More Over the past few years several studies have provided estimates of the SFR (star-formation rate) or the total infrared luminosity from just one infrared band. However these relations are generally derived for entire galaxies, which are known to contain a large scale diffuse emission that is not necessarily related to the latest star-formation episode. We provide new relations to estimate the SFR from resolved star-forming regions at 100 mum and 160 mum. We select individual star-forming regions in the nearby (840 kpc) galaxy M33. We estimate the SFR combining the emission in Halpha and at 24 mum to calibrate the emission at 100 mum and 160 mum as SFR estimators, as mapped with PACS/Herschel. The data are obtained in the framework of the HERM33ES open time key project. There is less emission in the HII regions at 160 mum than at 100 mum. Over a dynamic range of almost 2 dex in Sigma(SFR) we find that the 100 mum emission is a nearly linear estimator of the SFR, whereas that at 160 mum is slightly superlinear. The behaviour of individual star-forming regions is surprisingly similar to that of entire galaxies. At high Sigma(SFR), star formation drives the dust temperature, whereas uncertainties and variations in radiation-transfer and dust-heated processes dominate at low Sigma(SFR). Detailed modelling of both galaxies and individual star forming regions will be needed to interpret similarities and differences between the two and assess the fraction of diffuse emission in galaxies. △ Less

Submitted 12 May, 2010; originally announced May 2010.

Comments: 5 pages, 3 figures, accepted for publication in the A&A Herschel special issue