Showing 1–38 of 38 results for author: Ntormousi, E

Gas-phase metallicity gradients in galaxies at $z \sim 6-8$

Authors: G. Venturi, S. Carniani, E. Parlanti, M. Kohandel, M. Curti, A. Pallottini, L. Vallini, S. Arribas, A. J. Bunker, A. J. Cameron, M. Castellano, A. Ferrara, A. Fontana, S. Gallerani, V. Gelli, R. Maiolino, E. Ntormousi, C. Pacifici, L. Pentericci, S. Salvadori, E. Vanzella

Abstract: The study of gas-phase metallicity and its spatial distribution at high redshift is crucial to understand the processes that shaped the growth and evolution of galaxies in the early Universe. Here we study the spatially resolved metallicity in three systems at $z\sim6-8$, namely A2744-YD4, BDF-3299, and COSMOS24108, with JWST NIRSpec IFU low-resolution ($R$ $\sim$ 100) spectroscopic observations.… ▽ More The study of gas-phase metallicity and its spatial distribution at high redshift is crucial to understand the processes that shaped the growth and evolution of galaxies in the early Universe. Here we study the spatially resolved metallicity in three systems at $z\sim6-8$, namely A2744-YD4, BDF-3299, and COSMOS24108, with JWST NIRSpec IFU low-resolution ($R$ $\sim$ 100) spectroscopic observations. These are among the highest-$z$ sources in which metallicity gradients have been probed so far. Each of these systems hosts several spatial components in the process of merging within a few kiloparsecs, identified from the rest-frame UV and optical stellar continuum and ionised gas emission line maps. The sources have heterogeneous properties, with stellar masses log($M_*/M_\odot) \sim 7.6-9.3$, star formation rates (SFRs) $\sim1-15$ $M_\odot$ yr$^{-1}$, and gas-phase metallicities 12+log(O/H) $\sim 7.7-8.3$, which exhibit a large scatter within each system. Their properties are generally consistent with those of the highest-redshift samples to date ($z\sim3-10$), though the sources in A2744-YD4 and COSMOS24108 are at the high end of the mass-metallicity relation (MZR) defined by the $z\sim3-10$ sources. Moreover, the targets in this work follow the predicted slope of the MZR at $z\sim 6-8$ from most cosmological simulations. The gas-phase metallicity gradients are consistent with being flat in the main sources of each system. Flat metallicity gradients are thought to arise from gas mixing processes on galaxy scales, such as mergers or galactic outflows and SN winds driven by intense stellar feedback, which wash out any gradient formed in the galaxy. The existence of flat gradients at $z\sim6-8$ sets also important constraints on cosmological simulations, whose predictions on the cosmic evolution of metallicity gradients often differ significantly, especially at high redshift. △ Less

Submitted 6 March, 2024; originally announced March 2024.

Comments: 18 pages (14 main body; 4 appendix); 15 figures (7 main body; 8 appendix); submitted to A&A

Time evolution of the galactic $B- ρ$ relation: the impact of the magnetic field morphology

Authors: A. Konstantinou, E. Ntormousi, K. Tassis, A. Pallottini

Abstract: One of the most frequently used indicators to characterize the magnetic field's influence on star formation is the relation between magnetic field strength and gas density ($B-ρ$ relation), usually expressed as $B \propto ρ^κ$. The value of $κ$ is an indication of the dynamical importance of the magnetic field during gas compression. Investigating the global magnetic field's impact on this relatio… ▽ More One of the most frequently used indicators to characterize the magnetic field's influence on star formation is the relation between magnetic field strength and gas density ($B-ρ$ relation), usually expressed as $B \propto ρ^κ$. The value of $κ$ is an indication of the dynamical importance of the magnetic field during gas compression. Investigating the global magnetic field's impact on this relation and its evolution, we conduct MHD simulations of Milky-Way-like galaxies including gravity, star formation, and supernova feedback along with non-equilibrium chemistry up to $H_2$ formation fueling star formation. Two initial magnetic field morphologies are studied: one completely ordered (toroidal) and the other completely random. In these models, we study the dynamical importance of the magnetic field through the plasma $β$ and the $B-ρ$ relation. For both magnetic morphologies, low-density regions are thermally supported, while high-density regions are magnetically dominated. Equipartition is reached earlier and at lower densities in the toroidal model. However, the $B-ρ$ relation is not unique even within the same galaxy, as it consistently includes two different branches for a given density, with $κ$ ranging from about 0.2 to 0.8. The mean value of $κ$ for each model also displays significant variations over time, which supersede the differences between the two models. While our findings suggest that the magnetic field morphology does influence the galactic $B-ρ$ relation, its impact is transient, since time-averaged differences between the models fall within the large temporal scatter. The context and time-dependent nature of the $B-ρ$ relation underscore the need for comprehensive research and observations to understand the intricate role of magnetic fields in star formation processes across diverse galactic environments. △ Less

Submitted 19 February, 2024; v1 submitted 15 February, 2024; originally announced February 2024.

Comments: The paper has been accepted for publication at A&A

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

Extragalactic Magnetism with SOFIA (SALSA Legacy Program). VII. A Tomographic View of Far-infrared and Radio Polarimetric Observations through MHD Simulations of Galaxies

Authors: Sergio Martin-Alvarez, Enrique Lopez-Rodriguez, Tara Dacunha, Susan E. Clark, Alejandro S. Borlaff, Rainer Beck, Francisco Rodríguez Montero, S. Lyla Jung, Julien Devriendt, Adrianne Slyz, Julia Roman-Duval, Evangelia Ntormousi, Mehrnoosh Tahani, Kandaswamy Subramanian, Daniel A. Dale, Pamela M. Marcum, Konstantinos Tassis, Ignacio del Moral-Castro, Le Ngoc Tram, Matt J. Jarvis

Abstract: The structure of magnetic fields in galaxies remains poorly constrained, despite the importance of magnetism in the evolution of galaxies. Radio synchrotron and far-infrared (FIR) polarization and polarimetric observations are the best methods to measure galactic scale properties of magnetic fields in galaxies beyond the Milky Way. We use synthetic polarimetric observations of a simulated galaxy t… ▽ More The structure of magnetic fields in galaxies remains poorly constrained, despite the importance of magnetism in the evolution of galaxies. Radio synchrotron and far-infrared (FIR) polarization and polarimetric observations are the best methods to measure galactic scale properties of magnetic fields in galaxies beyond the Milky Way. We use synthetic polarimetric observations of a simulated galaxy to identify and quantify the regions, scales, and interstellar medium (ISM) phases probed at FIR and radio wavelengths. Our studied suite of magnetohydrodynamical cosmological zoom-in simulations features high-resolutions (10 pc full-cell size) and multiple magnetization models. Our synthetic observations have a striking resemblance to those of observed galaxies. We find that the total and polarized radio emission extends to approximately double the altitude above the galactic disk (half-intensity disk thickness of $h_\text{I radio} \sim h_\text{PI radio} = 0.23 \pm 0.03$ kpc) relative to the total FIR and polarized emission that are concentrated in the disk midplane ($h_\text{I FIR} \sim h_\text{PI FIR} = 0.11 \pm 0.01$ kpc). Radio emission traces magnetic fields at scales of $\gtrsim 300$ pc, whereas FIR emission probes magnetic fields at the smallest scales of our simulations. These scales are comparable to our spatial resolution and well below the spatial resolution ($<300$ pc) of existing FIR polarimetric measurements. Finally, we confirm that synchrotron emission traces a combination of the warm neutral and cold neutral gas phases, whereas FIR emission follows the densest gas in the cold neutral phase in the simulation. These results are independent of the ISM magnetic field strength. The complementarity we measure between radio and FIR wavelengths motivates future multiwavelength polarimetric observations to advance our knowledge of extragalactic magnetism. △ Less

Submitted 24 March, 2024; v1 submitted 10 November, 2023; originally announced November 2023.

Comments: Accepted for publication in the ApJ. Replaced to match final version. 35 pages, 17 figures, 4 tables

Bright-Moon Sky as a Wide-Field Linear Polarimetric Flat Source for Calibration

Authors: S. Maharana, S. Kiehlmann, D. Blinov, V. Pelgrims, V. Pavlidou, K. Tassis, J. A. Kypriotakis, A. N. Ramaprakash, R. M. Anche, A. Basyrov, K. Deka, H. K. Eriksen, T. Ghosh, E. Gjerløw, N. Mandarakas, E. Ntormousi, G. V. Panopoulou, A. Papadaki, T. Pearson, S. B. Potter, A. C. S. Readhead, R. Skalidis, I. K. Wehus

Abstract: Next-generation wide-field optical polarimeters like the Wide-Area Linear Optical Polarimeters (WALOPs) have a field of view (FoV) of tens of arcminutes. For efficient and accurate calibration of these instruments, wide-field polarimetric flat sources will be essential. Currently, no established wide-field polarimetric standard or flat sources exist. This paper tests the feasibility of using the p… ▽ More Next-generation wide-field optical polarimeters like the Wide-Area Linear Optical Polarimeters (WALOPs) have a field of view (FoV) of tens of arcminutes. For efficient and accurate calibration of these instruments, wide-field polarimetric flat sources will be essential. Currently, no established wide-field polarimetric standard or flat sources exist. This paper tests the feasibility of using the polarized sky patches of the size of around ten-by-ten arcminutes, at a distance of up to 20 degrees from the Moon, on bright-Moon nights as a wide-field linear polarimetric flat source. We observed 19 patches of the sky adjacent to the bright-Moon with the RoboPol instrument in the SDSS-r broadband filter. These were observed on five nights within two days of the full-Moon across two RoboPol observing seasons. We find that for 18 of the 19 patches, the uniformity in the measured normalized Stokes parameters $q$ and $u$ is within 0.2 %, with 12 patches exhibiting uniformity within 0.07 % or better for both $q$ and $u$ simultaneously, making them reliable and stable wide-field linear polarization flats. We demonstrate that the sky on bright-Moon nights is an excellent wide-field linear polarization flat source. Various combinations of the normalized Stokes parameters $q$ and $u$ can be obtained by choosing suitable locations of the sky patch with respect to the Moon △ Less

Submitted 7 May, 2023; originally announced May 2023.

Comments: 8 pages including appendix, 6 figures and 3 tables. Submitted to Astronomy and Astrophysics for review. Comments are welcome

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

Extragalactic magnetism with SOFIA (SALSA Legacy Program) -- V: First results on the magnetic field orientation of galaxies

Authors: Alejandro S. Borlaff, Enrique Lopez-Rodriguez, Rainer Beck, Susan E. Clark, Evangelia Ntormousi, Konstantinos Tassis, Sergio Martin-Alvarez, Mehrnoosh Tahani, Daniel A. Dale, Ignacio del Moral Castro, Julia Roman-Duval, Pamela M. Marcum, John E. Beckman, Kandaswamy Subramanian, Sarah Eftekharzadeh, Leslie Proudfit

Abstract: We present the analysis of the magnetic field ($B$-field) structure of galaxies measured with far-infrared (FIR) and radio (3 and 6 cm) polarimetric observations. We use the first data release of the Survey on extragALactic magnetiSm with SOFIA (SALSA) of 14 nearby ($<20$ Mpc) galaxies with resolved (5 arcsec-18 arcsec; $90$ pc--$1$ kpc) imaging polarimetric observations using HAWC+/SOFIA from… ▽ More We present the analysis of the magnetic field ($B$-field) structure of galaxies measured with far-infrared (FIR) and radio (3 and 6 cm) polarimetric observations. We use the first data release of the Survey on extragALactic magnetiSm with SOFIA (SALSA) of 14 nearby ($<20$ Mpc) galaxies with resolved (5 arcsec-18 arcsec; $90$ pc--$1$ kpc) imaging polarimetric observations using HAWC+/SOFIA from $53$ to $214$ \um. We compute the magnetic pitch angle ($Ψ_{B}$) profiles as a function of the galactrocentric radius. We introduce a new magnetic alignment parameter ($ζ$) to estimate the disordered-to-ordered $B$-field ratio in spiral $B$-fields. We find FIR and radio wavelengths to not generally trace the same $B$-field morphology in galaxies. The $Ψ_{B}$ profiles tend to be more ordered with galactocentric radius in radio ($ζ_{\rm{6cm}} = 0.93\pm0.03$) than in FIR ($ζ_{\rm{154μm}} = 0.84\pm0.14$). For spiral galaxies, FIR $B$-fields are $2-75$\% more turbulent than the radio $B$-fields. For starburst galaxies, we find that FIR polarization is a better tracer of the $B$-fields along the galactic outflows than radio polarization. Our results suggest that the $B$-fields associated with dense, dusty, turbulent star-forming regions, those traced at FIR, are less ordered than warmer, less-dense regions, those traced at radio, of the interstellar medium. The FIR $B$-fields seem to be more sensitive to the activity of the star-forming regions and the morphology of the molecular clouds within a vertical height of few hundred pc in the disk of spiral galaxies than the radio $B$-fields. △ Less

Submitted 13 June, 2023; v1 submitted 23 March, 2023; originally announced March 2023.

Comments: 28 pages, 14 figures. Accepted for publication in The Astrophysical Journal

Witnessing the fragmentation of a filament into prestellar cores in Orion B/NGC 2024

Authors: Y. Shimajiri, Ph. André, N. Peretto, D. Arzoumanian, E. Ntormousi, V. Könyves

Abstract: Recent Herschel observations of nearby clouds have shown that filamentary structures are ubiquitous and that most prestellar cores form in filaments. Probing the density ($n$) and velocity ($V$) structure of filaments is crucial for the understanding of the star formation process. To characterize both the $n$ and the $V$ field of a fragmenting filament, we mapped NGC2024. 13CO, C18O, and H13CO+ tr… ▽ More Recent Herschel observations of nearby clouds have shown that filamentary structures are ubiquitous and that most prestellar cores form in filaments. Probing the density ($n$) and velocity ($V$) structure of filaments is crucial for the understanding of the star formation process. To characterize both the $n$ and the $V$ field of a fragmenting filament, we mapped NGC2024. 13CO, C18O, and H13CO+ trace the filament seen in the $N_{H_2}$ data. The radial profile from the $N_{H_2}$ data shows $D_{HP}$~0.081 pc, which is similar to the Herschel findings. The $D_{HP}$ from 13CO and C18O are broader, while the $D_{HP}$ from H13CO+ is narrower, than $D_{HP}$ from Herschel. These results suggest that 13CO and C18O trace only the outer part of the filament and H13CO+ only the inner part. The H13CO+ $V_{centroid}$ map reveals $V$ gradients along both filament axis, as well as $V$ oscillations with a period $λ$~0.2 pc along the major axis. Comparison between the $V$ and the $n$ distribution shows a tentative $λ$/4 shift in H13CO+ or C18O. This $λ$/4 shift is not simultaneously observed for all cores in any single tracer but is tentatively seen in either H13CO+ or C18O. We produced a toy model taking into account a transverse $V$ gradient, a longitudinal $V$ gradient, and a longitudinal oscillation mode caused by fragmentation. Examination of synthetic data shows that the oscillation component produces an oscillation pattern in the velocity structure function (VSF) of the model. The H13CO+ VSF shows an oscillation pattern, suggesting that our observations are partly tracing core-forming motions and fragmentation. We also found that the mean $M_{core}$ corresponds to the effective $M_{BE}$ in the filament. This is consistent with a scenario in which higher-mass cores form in higher line-mass filaments. △ Less

Submitted 12 February, 2023; v1 submitted 25 January, 2023; originally announced January 2023.

Comments: accepted in A&A

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

A closer look at supernovae as seeds for galactic magnetization

Authors: Evangelia Ntormousi, Fabio Del Sordo, Matteo Cantiello, Andrea Ferrara

Abstract: Explaining the currently observed magnetic fields in galaxies requires relatively strong seeding in the early Universe. One theory proposes that magnetic fields of the order of $μ$G were expelled by supernova (SN) explosions after primordial, nG or weaker fields were amplified in stellar interiors. In this work, we calculate the maximum magnetic energy that can be injected in the interstellar medi… ▽ More Explaining the currently observed magnetic fields in galaxies requires relatively strong seeding in the early Universe. One theory proposes that magnetic fields of the order of $μ$G were expelled by supernova (SN) explosions after primordial, nG or weaker fields were amplified in stellar interiors. In this work, we calculate the maximum magnetic energy that can be injected in the interstellar medium by a stellar cluster of mass $M_{cl}$ based on what is currently known about stellar magnetism. We consider early-type stars and adopt either a Salpeter or a top-heavy IMF. For their magnetic fields, we adopt either a Gaussian or a bimodal distribution. The Gaussian model assumes that all massive stars are magnetized with $10^3 < \langle B_* \rangle < 10^4$ G, while the bimodal, consistent with observations of Milky Way stars, assumes only 5-10 per cent of OB stars have $10^3 < \langle B_* \rangle < 10^4$ G, while the rest have $10 < \langle B_* \rangle < 10^2$ G. We find that the maximum magnetic energy that can be injected by a stellar population is between $10^{-10}-10^{-7}$ times the total SN energy. The highest end of these estimates is about five orders of magnitude lower than what is usually employed in cosmological simulations, where about $10^{-2}$ of the SN energy is injected as magnetic. Pure advection of the stellar magnetic field by SN explosions is a good candidate for seeding a dynamo, but not enough to magnetize galaxies. Assuming SNe as main mechanism for galactic magnetization, the magnetic field cannot exceed an intensity of $10^{-7}$ G in the best-case scenario for a population of $10^{5}$ solar masses in a superbubble of 300 pc radius, while more typical values are between $10^{-10}-10^{-9}$~G. Therefore, other scenarios for galactic magnetization at high redshift need to be explored. △ Less

Submitted 22 November, 2022; originally announced November 2022.

Comments: Accepted for publication in A&A letters

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

Galactic Dynamos

Authors: Axel Brandenburg, Evangelia Ntormousi

Abstract: Spiral galaxies, including the Milky Way, have large-scale magnetic fields with significant energy densities. The dominant theory attributes these magnetic fields to a large-scale dynamo. We review the current status of dynamo theory and discuss various numerical simulations designed to explain either particular aspects of the problem or to reproduce galactic magnetic fields globally. Our main con… ▽ More Spiral galaxies, including the Milky Way, have large-scale magnetic fields with significant energy densities. The dominant theory attributes these magnetic fields to a large-scale dynamo. We review the current status of dynamo theory and discuss various numerical simulations designed to explain either particular aspects of the problem or to reproduce galactic magnetic fields globally. Our main conclusions can be summarized as follows. Idealized direct numerical simulations produce mean magnetic fields, whose saturation energy density tends to decline with increasing magnetic Reynolds number. This is still an unsolved problem. Large-scale galactic magnetic fields of microgauss strengths can probably only be explained if helical magnetic fields of small or moderate length scales can rapidly be ejected or destroyed. Small-scale dynamos are important throughout a galaxy's life, and probably provide strong seed fields at early stages. The circumgalactic medium (CGM) may play an important role in driving dynamo action at small and large length scales. These interactions between the galactic disk and the CGM may provide important insights into our understanding of galactic dynamos. We expect future research in galactic dynamos to focus on the cosmological history of galaxies and the interaction with the CGM as means of replacing the idealized boundary conditions used in earlier work. △ Less

Submitted 30 August, 2023; v1 submitted 7 November, 2022; originally announced November 2022.

Comments: 47 pages, 15 figures, 4 tables

Report number: NORDITA-2022-077

Journal ref: Annu. Rev. Astron. Astrophys. 61, 561-606 (2023)

Extragalactic magnetism with SOFIA (SALSA Legacy Program). VI. The magnetic fields in the multi-phase interstellar medium of the Antennae galaxies

Authors: Enrique Lopez-Rodriguez, Alejandro S. Borlaff, Rainer Beck, William T. Reach, Sui Ann Mao, Evangelia Ntormousi, Konstantinos Tassis, Sergio Martin-Alvarez, Susan E. Clark, Daniel A. Dale, Ignacio del Moral-Castro

Abstract: Mergers are thought to be a fundamental channel for galaxy growth, perturbing the gas dynamics and the magnetic fields (B-fields) in the interstellar medium (ISM). However, the mechanisms that amplify and dissipate B-fields during a merger remain unclear. We characterize the morphology of the ordered B-fields in the multi-phase ISM of the closest merger of two spiral galaxies, the Antennae galaxie… ▽ More Mergers are thought to be a fundamental channel for galaxy growth, perturbing the gas dynamics and the magnetic fields (B-fields) in the interstellar medium (ISM). However, the mechanisms that amplify and dissipate B-fields during a merger remain unclear. We characterize the morphology of the ordered B-fields in the multi-phase ISM of the closest merger of two spiral galaxies, the Antennae galaxies. We compare the inferred B-fields using $154~μ$m thermal dust and $11$ cm radio synchrotron emission polarimetric observations. We find that the $154~μ$m B-fields are more ordered across the Antennae galaxies than the $11$ cm B-fields. The turbulent-to-ordered $154~μ$m B-field increases at the galaxy cores and star-forming regions. The relic spiral arm has an ordered spiral $154~μ$m B-field, while the $11$ cm B-field is radial. The $154~μ$m B-field may be dominated by turbulent dynamos with high $^{12}$CO(1-0) velocity dispersion driven by star-forming regions, while the $11$ cm B-field is cospatial with high HI velocity dispersion driven by galaxy interaction. This result shows the dissociation between the warm gas mainly disturbed by the merger, and the dense gas still following the dynamics of the relic spiral arm. We find a $\sim8.9$ kpc scale ordered B-field connecting the two galaxies. The base of the tidal tail is cospatial with the HI and $^{12}$CO(1-0) emission and has compressed and/or sheared $154~μ$m and $11$ cm B-fields driven by the merger. We suggest that amplify B-fields, with respect to the rest of the system and other spiral galaxies, may be supporting the gas flow between both galaxies and the tidal tail. △ Less

Submitted 9 December, 2022; v1 submitted 31 October, 2022; originally announced November 2022.

Comments: 11 pages, 5 figures, Accepted for publication in ApJ Letters

Starlight-polarization-based tomography of the magnetized interstellar medium: PASIPHAE's line-of-sight inversion method

Authors: V. Pelgrims, G. V. Panopoulou, K. Tassis, V. Pavlidou, A. Basyrov, D. Blinov, E. Gjerløw, S. Kiehlmann, N. Mandarakas, A. Papadaki, R. Skalidis, A. Tsouros, R. M. Anche, H. K. Eriksen, T. Ghosh, J. A. Kypriotakis, S. Maharana, E. Ntormousi, T. J. Pearson, S. B. Potter, A. N. Ramaprakash, A. C. S. Readhead, I. K. Wehus

Abstract: We present the first Bayesian method for tomographic decomposition of the plane-of-sky orientation of the magnetic field with the use of stellar polarimetry and distance. This standalone tomographic inversion method presents an important step forward in reconstructing the magnetized interstellar medium (ISM) in 3D within dusty regions. We develop a model in which the polarization signal from the m… ▽ More We present the first Bayesian method for tomographic decomposition of the plane-of-sky orientation of the magnetic field with the use of stellar polarimetry and distance. This standalone tomographic inversion method presents an important step forward in reconstructing the magnetized interstellar medium (ISM) in 3D within dusty regions. We develop a model in which the polarization signal from the magnetized and dusty ISM is described by thin layers at various distances. Our modeling makes it possible to infer the mean polarization (amplitude and orientation) induced by individual dusty clouds and to account for the turbulence-induced scatter in a generic way. We present a likelihood function that explicitly accounts for uncertainties in polarization and parallax. We develop a framework for reconstructing the magnetized ISM through the maximization of the log-likelihood using a nested sampling method. We test our Bayesian inversion method on mock data taking into account realistic uncertainties from Gaia and as expected for the optical polarization survey PASIPHAE according to the currently planned observing strategy. We demonstrate that our method is effective at recovering the cloud properties as soon as the polarization induced by a cloud to its background stars is higher than $\sim 0.1\%$ for the adopted survey exposure time and level of systematic uncertainty. Our method makes it possible to recover not only the mean polarization properties but also to characterize the intrinsic scatter, thus creating new ways to characterize ISM turbulence and the magnetic field strength. Finally, we apply our method to an existing data set of starlight polarization with known line-of-sight decomposition, demonstrating agreement with previous results and an improved quantification of uncertainties in cloud properties. △ Less

Submitted 28 February, 2023; v1 submitted 3 August, 2022; originally announced August 2022.

Comments: Closely matches the published version: A&A 670, A164, 2023 The accompanying code (BISP-1) is available at https://github.com/vpelgrims/Bisp_1

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

Extragalactic magnetism with SOFIA (SALSA Legacy Program) -- IV: Program overview and first results on the polarization fraction

Authors: Enrique Lopez-Rodriguez, Sui Ann Mao, Rainer Beck, Alejandro S. Borlaff, Evangelia Ntormousi, Konstantinos Tassis, Daniel A. Dale, Julia Roman-Duval, Kandaswamy Subramanian, Sergio Martin-Alvarez, Pamela M. Marcum, Susan E. Clark, William T. Reach, Doyal A. Harper, Ellen G. Zweibel

Abstract: We present the first data release of the Survey on extragALactic magnetiSm with SOFIA (SALSA Legacy Program) with a set of 14 nearby ($<20$ Mpc) galaxies with resolved imaging polarimetric observations using HAWC+ from $53$ to $214$ $μ$m at a resolution of $5-18$" ($90$ pc $-$ $1$ kpc). We introduce the definitions and background on extragalactic magnetism, and present the scientific motivation an… ▽ More We present the first data release of the Survey on extragALactic magnetiSm with SOFIA (SALSA Legacy Program) with a set of 14 nearby ($<20$ Mpc) galaxies with resolved imaging polarimetric observations using HAWC+ from $53$ to $214$ $μ$m at a resolution of $5-18$" ($90$ pc $-$ $1$ kpc). We introduce the definitions and background on extragalactic magnetism, and present the scientific motivation and sample selection of the program. Here, we focus on the general trends in the emissive polarization fraction. Far-infrared polarimetric observations trace the thermal polarized emission of magnetically aligned dust grains across the galaxy disks with polarization fractions of $P=0-15$% in the cold, $T_{\rm d} = [19,48]$ K, and dense, $\log_{10}(N_{\rm HI+H_{2}}) = [19.96,22.91]$, interstellar medium. The spiral galaxies show a median $\langle P_{154μm} \rangle = 3.3\pm0.9 $% across the disks. We report the first polarized spectrum of starburst galaxies showing a minimum within $89-154$ $μ$m. The falling $53-154$ $μ$m polarized spectrum may be due to a decrease in the dust grain alignment efficiency produced by variations in dust temperatures along the line-of-sight in the galactic outflow. We find that the starburst galaxies and the star-forming regions within normal galaxies have the lowest polarization fractions. We find that 50% (7 out of 14) of the galaxies require a broken power-law in the $P-N_{HI+H_{2}}$ and $P-T_{d}$ relations with three different trends. Group 1 has a relative increase of anisotropic random B-fields produced by compression or shear of B-fields in the galactic outflows, starburst rings, and inner-bar of galaxies; and Groups 2 and 3 have a relative increase of isotropic random B-fields driven by star-forming regions in the spiral arms, and/or an increase of dust grain alignment efficiency caused by shock-driven regions or evolutionary stages of a galaxy. △ Less

Submitted 7 July, 2022; v1 submitted 2 May, 2022; originally announced May 2022.

Comments: 43 pages, 28 figures. Accepted for publication to ApJ

Extragalactic magnetism with SOFIA (SALSA Legacy Program) -- III: First data release and on-the-fly polarization mapping characterization

Authors: Enrique Lopez-Rodriguez, Melanie Clarke, Sachin Shenoy, William Vacca, Simon Coude, Ryan Arneson, Peter Ashton, Sarah Eftekharzadeh, Rainer Beck, John E. Beckman, Alejandro S. Borlaff, Susan E. Clark, Daniel A. Dale, Sergio Martin-Alvarez, Evangelia Ntormousi, William T. Reach, Julia Roman-Duval, Konstantinos Tassis, Doyal A. Harper, Pamela M. Marcum

Abstract: We describe the data processing of the Survey on extragALactic magnetiSm with SOFIA (SALSA Legacy Program). This first data release presents 33% (51.34h out of 155.7h, including overheads) of the total awarded time taken from January 2020 to December 2021. Our observations were performed using the newly implemented on-the-fly mapping (OTFMAP) technique in the polarimetric mode. We present the pipe… ▽ More We describe the data processing of the Survey on extragALactic magnetiSm with SOFIA (SALSA Legacy Program). This first data release presents 33% (51.34h out of 155.7h, including overheads) of the total awarded time taken from January 2020 to December 2021. Our observations were performed using the newly implemented on-the-fly mapping (OTFMAP) technique in the polarimetric mode. We present the pipeline steps to obtain homogeneously reduced high-level data products of polarimetric maps of galaxies for use in scientific analysis. Our approach has a general design and can be applied to sources smaller than the field-of-view of the HAWC+ array in any given band. We estimate that the OTFMAP polarimetric mode offers a reduction of observing overheads by a factor 2.34, and an improvement in sensitivity by a factor 1.80 when compared to previously obtained polarimetric observations using the chopping and nodding mode. The OTFMAP is a significant optimization of the polarimetric mode of HAWC+ as it ultimately reduces the cost of operations of SOFIA/HAWC+ by increasing the science collected per hour of observation up to an overall factor of 2.49. The OTFMAP polarimetric mode is the standard observing strategy of SALSA. The results and quantitative analysis of this first data release are presented in Papers IV and V of the series. △ Less

Submitted 28 April, 2022; originally announced April 2022.

Comments: 26 pages, 16 figures, 5 tables. Submitted to ApJ (Comments welcome)

First look at the multiphase interstellar medium with synthetic observations of low-frequency Faraday tomography

Authors: A. Bracco, E. Ntormousi, V. Jelić, M. Padovani, B. Šiljeg, A. Erceg, L. Turić, L. Ceraj, I. Šnidarić

Abstract: Faraday tomography of radio polarimetric data below 200 MHz from LOFAR are providing us with a new perspective on the diffuse and magnetized interstellar medium (ISM). Of particular interest is the discovery of Faraday-rotated synchrotron polarization associated with neutral gas, as traced by atomic hydrogen (HI) and dust. Here we present the first in-depth numerical study of these LOFAR results.… ▽ More Faraday tomography of radio polarimetric data below 200 MHz from LOFAR are providing us with a new perspective on the diffuse and magnetized interstellar medium (ISM). Of particular interest is the discovery of Faraday-rotated synchrotron polarization associated with neutral gas, as traced by atomic hydrogen (HI) and dust. Here we present the first in-depth numerical study of these LOFAR results. We produce and analyze comprehensive synthetic observations of low-frequency synchrotron polarization from MHD simulations of colliding super shells in the multiphase ISM. We define five distinct gas phases over more than four orders of magnitude in gas temperature and density, ranging from hot, and warm, fully ionized gas to cold neutral medium. We focus on the contribution of each gas phase to synthetic observations of both rotation measure and synchrotron polarized intensity below 200 MHz. We also investigate the link between the latter and synthetic observations of optically thin HI gas. We find that, not only the fully ionized gas but also the warm partially ionized and neutral phases strongly contribute to the total rotation measure and polarized intensity. However, the contribution of each phase to the observables strongly depends on the choice of integration axis and the orientation of the mean magnetic field with respect to the shell-collision axis. Strong correlation between HI synthetic data and synchrotron polarized intensity, reminiscent of LOFAR results, is obtained with lines of sight perpendicular to the mean magnetic field direction. Our study suggests that multiphase modelling of MHD processes is needed in order to interpret observations of the radio sky at low frequency. This work is a first step toward understanding the complexity of low-frequency synchrotron emission that will be soon revolutionized by large-scale surveys with LOFAR and the SKA. △ Less

Submitted 15 April, 2022; v1 submitted 6 April, 2022; originally announced April 2022.

Comments: 17 pages, 20 figures, accepted by A&A, in press

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

Magnetic fields in star formation: from clouds to cores

Authors: Kate Pattle, Laura Fissel, Mehrnoosh Tahani, Tie Liu, Evangelia Ntormousi

Abstract: In this chapter we review recent advances in understanding the roles that magnetic fields play throughout the star formation process, gained through observations and simulations of molecular clouds, the dense, star-forming phase of the magnetised, turbulent interstellar medium (ISM). Recent results broadly support a picture in which the magnetic fields of molecular clouds transition from being gra… ▽ More In this chapter we review recent advances in understanding the roles that magnetic fields play throughout the star formation process, gained through observations and simulations of molecular clouds, the dense, star-forming phase of the magnetised, turbulent interstellar medium (ISM). Recent results broadly support a picture in which the magnetic fields of molecular clouds transition from being gravitationally sub-critical and near equipartition with turbulence in low-density cloud envelopes, to being energetically sub-dominant in dense, gravitationally unstable star-forming cores. Magnetic fields appear to play an important role in the formation of cloud substructure by setting preferred directions for large-scale gas flows in molecular clouds, and can direct the accretion of material onto star-forming filaments and hubs. Low-mass star formation may proceed in environments close to magnetic criticality; high-mass star formation remains less well-understood, but may proceed in more supercritical environments. The interaction between magnetic fields and (proto)stellar feedback may be particularly important in setting star formation efficiency. We also review a range of widely-used techniques for quantifying the dynamic importance of magnetic fields, concluding that better-calibrated diagnostics are required in order to use the spectacular range of forthcoming observations and simulations to quantify our emerging understanding of how magnetic fields influence the outcome of the star formation process. △ Less

Submitted 21 March, 2022; originally announced March 2022.

Comments: 34 pages, 10 figures, to appear in Protostars and Planets VII conference proceedings. Editors: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura. All figures reproduced with permission of the authors. Comments welcome before 30/06/2022

The effect of cosmic variance on the characteristics of dust polarization power spectra

Authors: V. Pelgrims, E. Ntormousi, K. Tassis

Abstract: In the context of cosmic microwave background polarization studies and the characterization of the Galactic foregrounds, the power spectrum analysis of the thermal dust polarization sky has led to intriguing evidence of an E/B asymmetry and a positive TE correlation. In this work, we produce synthesized dust polarization maps from a set of global magneto-hydrodynamic (MHD) simulations of Milky-Way… ▽ More In the context of cosmic microwave background polarization studies and the characterization of the Galactic foregrounds, the power spectrum analysis of the thermal dust polarization sky has led to intriguing evidence of an E/B asymmetry and a positive TE correlation. In this work, we produce synthesized dust polarization maps from a set of global magneto-hydrodynamic (MHD) simulations of Milky-Way-sized galaxies, and analyze their power spectra at intermediate angular scales (angular multipoles $\ell \in \left[60 ,\, 140\right]$). We study the role of the initial configuration of the large-scale magnetic field, its strength, and the feedback on the power spectrum characteristics. Using full-galaxy MHD simulations, we were able to estimate the variance induced by the peculiar location of the observer in the galaxy. We find that the polarization power spectra sensitively depend on the observer's location, impeding a distinction between different simulation setups. There is a clear statistical difference between the power spectra measured from within the spiral arms and those measured from the inter-arm regions. Also, power spectra from within supernova-driven bubbles share common characteristics, regardless of the underlying model. However, no correlation was found between the properties of the polarization power spectra and the local (with respect to the observer) mean values of physical quantities such as the density and the strength of the magnetic field. Finally, we find indications that the global strength of the magnetic field may play a role in shaping the power spectrum characteristics; as the global magnetic field strength increases, the E/B asymmetry and the TE correlation increase, whereas the viewpoint-induced variance decreases. However, we find no direct correlation with the strength of the local magnetic field that permeates the mapped region of the interstellar medium. △ Less

Submitted 10 January, 2022; originally announced January 2022.

Comments: 20 pages including 3 appendices. Accepted version at A&A

Journal ref: A&A 658, A134 (2022)

Dynamo effect in unstirred self-gravitating turbulence

Authors: Axel Brandenburg, Evangelia Ntormousi

Abstract: In many astrophysical environments, self-gravity can generate kinetic energy, which, in principle, is available for driving dynamo action. Using direct numerical simulations, we show that in unstirred self-gravitating subsonic turbulence with helicity and a magnetic Prandtl number of unity, there is a critical magnetic Reynolds number of about 25 above which the work done against the Lorentz force… ▽ More In many astrophysical environments, self-gravity can generate kinetic energy, which, in principle, is available for driving dynamo action. Using direct numerical simulations, we show that in unstirred self-gravitating subsonic turbulence with helicity and a magnetic Prandtl number of unity, there is a critical magnetic Reynolds number of about 25 above which the work done against the Lorentz force exceeds the Ohmic dissipation. The collapse itself drives predominantly irrotational motions that cannot be responsible for dynamo action. We find that, with a weak magnetic field, one-third of the work done by the gravitational force goes into compressional heating and the remaining two-thirds go first into kinetic energy of the turbulence before a fraction of it is converted further into magnetic and finally thermal energies. Close to the collapse, however, these fractions change toward 1/4 and 3/4 for compressional heating and kinetic energy, respectively. When the magnetic field is strong, the compressional heating fraction is unchanged. Out of the remaining kinetic energy, one quarter goes directly into magnetic energy via work against the Lorentz force. The fraction of vortical motions diminishes in favor of compressive motions that are almost exclusively driven by the Jeans instability. For an initially uniform magnetic field, field amplification at scales larger than those of the initial turbulence are driven by tangling. △ Less

Submitted 14 March, 2022; v1 submitted 7 December, 2021; originally announced December 2021.

Comments: 17 pages, 19 figures, 1 table, resubmitted to MNRAS

Report number: NORDITA 2021-107

The width of Herschel filaments varies with distance

Authors: G. V. Panopoulou, S. E. Clark, A. Hacar, F. Heitsch, J. Kainulainen, E. Ntormousi, D. Seifried, R. J. Smith

Abstract: Context: Filamentary structures in nearby molecular clouds have been found to exhibit a characteristic width of 0.1 pc, as observed in dust emission. Understanding the origin of this universal width has become a topic of central importance in the study of molecular cloud structure and the early stages of star formation. Aims: We investigate how the recovered widths of filaments depend on the dista… ▽ More Context: Filamentary structures in nearby molecular clouds have been found to exhibit a characteristic width of 0.1 pc, as observed in dust emission. Understanding the origin of this universal width has become a topic of central importance in the study of molecular cloud structure and the early stages of star formation. Aims: We investigate how the recovered widths of filaments depend on the distance from the observer, by using previously published results from the Herschel Gould Belt Survey. Methods: We obtained updated estimates on the distances to nearby molecular clouds observed with Herschel by using recent results based on 3D dust extinction mapping and Gaia. We examined the widths of filaments from individual clouds separately, as opposed to treating them as a single population. We used these per-cloud filament widths to search for signs of variation amongst the clouds of the previously published study. Results: We find a significant dependence of the mean per-cloud filament width with distance. The distribution of mean filament widths for nearby clouds is incompatible with that of farther away clouds. The mean per-cloud widths scale with distance approximately as 4--5 times the beam size. We examine the effects of resolution by performing a convergence study of a filament profile in the Herschel image of the Taurus Molecular Cloud. We find that resolution can severely affect the shapes of radial profiles over the observed range of distances. Conclusions: We conclude that the data are inconsistent with 0.1 pc being the universal characteristic width of filaments. △ Less

Submitted 8 August, 2022; v1 submitted 15 November, 2021; originally announced November 2021.

Comments: An erratum to this manuscript has been published here: https://www.aanda.org/articles/aa/full_html/2022/07/aa42281e-21/aa42281e-21.html

Journal ref: A&A 657, L13 (2022)

Extragalactic magnetism with SOFIA (Legacy Program) -- II: A magnetically-driven flow in the starburst ring of NGC 1097

Authors: Enrique Lopez-Rodriguez, Rainer Beck, Susan E. Clark, Annie Hughes, Alejandro S. Borlaff, Evangelia Ntormousi, Lucas Grosset, Konstantinos Tassis, John E. Beckman, Kandaswamy Subramanian, Daniel Dale, Tanio Diaz-Santos, Legacy Team

Abstract: Galactic bars are frequent in disk galaxies and they may support the transfer of matter towards the central engine of active nuclei. The barred galaxy NGC 1097 has magnetic forces controlling the gas flow at several kpc scales, which suggest that magnetic fields (B-fields) are dynamically important along the bar and nuclear ring. However, the effect of the B-field on the gas flows in the central k… ▽ More Galactic bars are frequent in disk galaxies and they may support the transfer of matter towards the central engine of active nuclei. The barred galaxy NGC 1097 has magnetic forces controlling the gas flow at several kpc scales, which suggest that magnetic fields (B-fields) are dynamically important along the bar and nuclear ring. However, the effect of the B-field on the gas flows in the central kpc scale has not been characterized. Using thermal polarized emission at $89$ $μ$m with HAWC+/SOFIA, here, we measure that the polarized flux is spatially located at the contact regions of the outer-bar with the starburst ring. The linear polarization decomposition analysis shows that the $89$ $μ$m and radio ($3.5$ and $6.2$ cm) polarization traces two different modes, $m$, of the B-field: a constant B-field orientation and dominated by $m=0$ at $89$ $μ$m, and a spiral B-field dominated by $m=2$ at radio. We show that the B-field at 89 $μ$m is concentrated in the warmest region of a shock driven by the galactic-bar dynamics in the contact regions between the outer-bar with the starburst ring. Radio polarization traces a superposition of the spiral B-field outside and within the starburst ring. According to Faraday rotation measures between $3.5$ and $6.2$ cm, the radial component of the B-field along the contact regions points toward the galaxy's center on both sides. We conclude that gas streams outside and within the starburst ring follow the B-field, which feeds the black hole with matter from the host galaxy. △ Less

Submitted 7 October, 2021; v1 submitted 19 July, 2021; originally announced July 2021.

Comments: 17 pages, 10 figures, Accepted for publication to ApJ

SMILE: Search for MIlli-LEnses

Authors: C. Casadio, D. Blinov, A. C. S. Readhead, I. W. A. Browne, P. N. Wilkinson, T. Hovatta, N. Mandarakas, V. Pavlidou, K. Tassis, H. K. Vedantham, J. A. Zensus, V. Diamantopoulos, K. E. Dolapsaki, K. Gkimisi, G. Kalaitzidakis, M. Mastorakis, K. Nikolaou, E. Ntormousi, V. Pelgrims, K. Psarras

Abstract: Dark Matter (DM) halos with masses below $\sim10^{8}$ $M_{\odot}$, which would help to discriminate between DM models, may be detected through their gravitational effect on distant sources. The same applies to primordial black holes, considered as an alternative scenario to DM particle models. However, there is still no evidence for the existence of such objects. With the aim of finding compact ob… ▽ More Dark Matter (DM) halos with masses below $\sim10^{8}$ $M_{\odot}$, which would help to discriminate between DM models, may be detected through their gravitational effect on distant sources. The same applies to primordial black holes, considered as an alternative scenario to DM particle models. However, there is still no evidence for the existence of such objects. With the aim of finding compact objects in the mass range $\sim$ 10$^{6}$ -- 10$^{9}$$M_{\odot}$, we search for strong gravitational lenses on milli (mas)-arcseconds scales (< 150 mas). For our search, we used the Astrogeo VLBI FITS image database -- the largest publicly available database, containing multi-frequency VLBI data of 13828 individual sources. We used the citizen science approach to visually inspect all sources in all available frequencies in search for images with multiple compact components on mas-scales. At the final stage, sources were excluded based on the surface brightness preservation criterion. We obtained a sample of 40 sources that passed all steps and therefore are judged to be milli-arcsecond lens candidates. These sources are currently followed-up with on-going European VLBI Network (EVN) observations at 5 and 22 GHz. Based on spectral index measurements, we suggest that two of our candidates have a higher probability to be associated with gravitational lenses. △ Less

Submitted 14 July, 2021; originally announced July 2021.

Comments: 6 pages, 1 figure; accepted for publication

Strong dependence of the physical properties of cores on spatial resolution in observations and simulations

Authors: Fabien Louvet, Patrick Hennebelle, Alexander Men'shchikov, Pierre Didelon, Evangelia Ntormousi, Pierre Didelon, Frédérique Motte

Abstract: During the last decade in star formation research, many studies have targeted low- and high-mass star formation regions located at different distances, with different telescopes having specific angular resolution capabilities. We present a systematic investigation of the angular resolution effects, with special attention being paid to the derived masses of sources as well as the shape of the resul… ▽ More During the last decade in star formation research, many studies have targeted low- and high-mass star formation regions located at different distances, with different telescopes having specific angular resolution capabilities. We present a systematic investigation of the angular resolution effects, with special attention being paid to the derived masses of sources as well as the shape of the resulting source mass functions (SMFs). We tested the impact of angular resolution, from 0.6 down to 0.02 pc, in two star-forming regions observed with Herschel (NGC6334 and Aquila), and three (magneto)-hydrodynamical simulations. We detected and measured sources at each resolution using getsf and we analysed the derived masses and sizes of the sources. We find that the number of sources does not converge from 0.6 to 0.05 pc. It increases by about two when the angular resolution increases with a similar factor. Below 0.05 pc, the number of source still increases by about 1.3 when the angular resolution increases by two, suggesting that we are close to, but not yet at, convergence. We find that the measured sizes and masses of sources linearly depend on the angular resolution with no sign of convergence to a resolution-independent value. The corresponding SMF peak also shifts with angular resolution, while the slope of the high-mass tail of the SMFs remains almost invariant. If prestellar cores, physically distinct from their background, exist in cluster-forming molecular clouds, we conclude that their mass must be lower than reported so far in the literature. We discuss various implications for the studies of star formation: the problem of determining the mass reservoirs involved in the star-formation process; the inapplicability of the Gaussian beam deconvolution to infer source sizes; and the impossibility to determine the efficiency of the mass conversion from the cores to the stars. △ Less

Submitted 12 July, 2021; originally announced July 2021.

Comments: Accepted by A&A

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

Extragalactic Magnetism with SOFIA (Legacy Program) -- I: The magnetic field in the multi-phase interstellar medium of M51

Authors: Alejandro S. Borlaff, Enrique Lopez-Rodriguez, Rainer Beck, Rodion Stepanov, Eva Ntormousi, Annie Hughes, Konstantinos Tassis, Pamela M. Marcum, Lucas Grosset, John E. Beckman, Leslie Proudfit, Susan E. Clark, Tanio Díaz-Santos, Sui Ann Mao, William T. Reach, Julia Roman-Duval, Kandaswamy Subramanian, Le Ngoc Tram, Ellen G. Zweibel, SOFIA Legacy Team

Abstract: The recent availability of high-resolution far-infrared (FIR) polarization observations of galaxies using HAWC+/SOFIA has facilitated studies of extragalactic magnetic fields in the cold and dense molecular disks.We investigate if any significant structural differences are detectable in the kpc-scale magnetic field of the grand design face-on spiral galaxy M51 when traced within the diffuse (radio… ▽ More The recent availability of high-resolution far-infrared (FIR) polarization observations of galaxies using HAWC+/SOFIA has facilitated studies of extragalactic magnetic fields in the cold and dense molecular disks.We investigate if any significant structural differences are detectable in the kpc-scale magnetic field of the grand design face-on spiral galaxy M51 when traced within the diffuse (radio) and the dense and cold (FIR) interstellar medium (ISM). Our analysis reveals a complex scenario where radio and FIR polarization observations do not necessarily trace the same magnetic field structure. We find that the magnetic field in the arms is wrapped tighter at 154um than at 3 and 6 cm; statistically significant lower values for the magnetic pitch angle are measured at FIR in the outskirts (R > 7 kpc) of the galaxy. This difference is not detected in the interarm region. We find strong correlations of the polarization fraction and total intensity at FIR and radio with the gas column density and 12CO(1-0) velocity dispersion. We conclude that the arms show a relative increase of small-scale turbulent B-fields at regions with increasing column density and dispersion velocities of the molecular gas. No correlations are found with HI neutral gas. The star formation rate shows a clear correlation with the radio polarized intensity, which is not found in FIR, pointing to a small-scale dynamo-driven B-field amplification scenario. This work shows that multi-wavelength polarization observations are key to disentangling the interlocked relation between star formation, magnetic fields, and gas kinematics in the multi-phase ISM. △ Less

Submitted 8 September, 2021; v1 submitted 19 May, 2021; originally announced May 2021.

Comments: 35 pages, 25 figures, 5 tables. Accepted for publication in ApJ

A dynamo amplifies the magnetic field of a Milky-Way-like galaxy

Authors: Evangelia Ntormousi, Konstantinos Tassis, Fabio Del Sordo, Francesca Fragkoudi, Rüdiger Pakmor

Abstract: The magnetic fields of spiral galaxies are so strong that they cannot be primordial. Their typical values are over one billion times higher than any value predicted for the early Universe. Explaining this immense growth and incorporating it in galaxy evolution theories is one of the long-standing challenges in astrophysics. So far, the most successful theory for the sustained growth of the galacti… ▽ More The magnetic fields of spiral galaxies are so strong that they cannot be primordial. Their typical values are over one billion times higher than any value predicted for the early Universe. Explaining this immense growth and incorporating it in galaxy evolution theories is one of the long-standing challenges in astrophysics. So far, the most successful theory for the sustained growth of the galactic magnetic field is the alpha-omega dynamo. This theory predicts a characteristic dipolar or quadrupolar morphology for the galactic magnetic field, which has been observed in external galaxies. However, so far, there has been no direct demonstration of a mean-field dynamo operating in direct, multi-physics simulations of spiral galaxies. We do so in this work. We employ numerical models of isolated, star-forming spiral galaxies that include a magnetized gaseous disk, a dark matter halo, stars, and stellar feedback. Naturally, the resulting magnetic field has a complex morphology that includes a strong random component. Using a smoothing of the magnetic field on small scales, we are able to separate the mean from the turbulent component and analyze them individually. We find that a mean-field dynamo naturally occurs as a result of the dynamical evolution of the galaxy and amplifies the magnetic field by an order of magnitude over half a Gyr. Despite the highly dynamical nature of these models, the morphology of the mean component of the field is identical to analytical predictions. This result underlines the importance of the mean-field dynamo in galactic evolution. Moreover, by demonstrating the natural growth of the magnetic field in a complex galactic environment, it brings us a step closer to understanding the cosmic origin of magnetic fields. △ Less

Submitted 22 June, 2020; originally announced June 2020.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 641, A165 (2020)

Turnaround radius of galaxy clusters in N-body simulations

Authors: Giorgos Korkidis, Vasiliki Pavlidou, Konstantinos Tassis, Evangelia Ntormousi, Theodore N. Tomaras, Konstantinos Kovlakas

Abstract: We use N-body simulations to examine whether a characteristic turnaround radius, as predicted from the spherical collapse model in a $\rm {ΛCDM}$ Universe, can be meaningfully identified for galaxy clusters, in the presence of full three-dimensional effects. We use The Dark Sky Simulations and Illustris-TNG dark-matter--only cosmological runs to calculate radial velocity profiles around collapsed… ▽ More We use N-body simulations to examine whether a characteristic turnaround radius, as predicted from the spherical collapse model in a $\rm {ΛCDM}$ Universe, can be meaningfully identified for galaxy clusters, in the presence of full three-dimensional effects. We use The Dark Sky Simulations and Illustris-TNG dark-matter--only cosmological runs to calculate radial velocity profiles around collapsed structures, extending out to many times the virial radius $R_{200}$. There, the turnaround radius can be unambiguously identified as the largest non-expanding scale around a center of gravity. We find that: (a) Indeed, a single turnaround scale can meaningfully describe strongly non-spherical structures. (b) For halos of masses $M_{200}>10^{13}M_\odot$, the turnaround radius $R_{ta}$ scales with the enclosed mass $M_{ta}$ as $M_{ta}^{1/3}$, as predicted by the spherical collapse model. (c) The deviation of $R_{ta}$ in simulated halos from the spherical collapse model prediction is insensitive to halo asphericity. Rather, it is sensitive to the tidal forces due to massive neighbors when such are present. (d) Halos exhibit a characteristic average density within the turnaround scale. This characteristic density is dependent on cosmology and redshift. For the present cosmic epoch and for concordance cosmological parameters ($Ω_m \sim 0.7$; $Ω_Λ\sim 0.3$) turnaround structures exhibit an average matter density contrast with the background Universe of $δ\sim 11$. Thus $R_{ta}$ is equivalent to $R_{11}$ -- in a way analogous to defining the "virial" radius as $R_{200}$ -- with the advantage that $R_{11}$ is shown in this work to correspond to a kinematically relevant scale in N-body simulations. △ Less

Submitted 17 December, 2019; originally announced December 2019.

Comments: Submitted to A&A, nine pages, 10 figures

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

Probing fragmentation and velocity sub-structure in the massive NGC 6334 filament with ALMA

Authors: Y. Shimajiri, Ph. Andre, E. Ntormousi, A. Men'shchikov, D. Arzoumanian, P. Palmeirim

Abstract: Herschel surveys of Galactic clouds support a paradigm for low-mass star formation in which dense filaments play a crucial role. The detailed fragmentation properties of star-forming filaments remain poorly understood, however, and the validity of the filament paradigm in the high-mass regime is still unclear. To investigate the density/velocity structure of the filament in the high-mass star-form… ▽ More Herschel surveys of Galactic clouds support a paradigm for low-mass star formation in which dense filaments play a crucial role. The detailed fragmentation properties of star-forming filaments remain poorly understood, however, and the validity of the filament paradigm in the high-mass regime is still unclear. To investigate the density/velocity structure of the filament in the high-mass star-forming region NGC6334, we conducted ALMA observations in the 3mm continuum and the N2H+(1-0) line at ~3arcsec resolution. The filament was detected in both tracers. We identified 26 cores at 3mm and 5 velocity-coherent fiber-like features in N2H+ within the filament. The typical length of, and velocity difference between, the fiber-like features of the NGC6334 filament are reminiscent of the properties for the fibers of the low-mass star-forming filament B211/B213. Only 2 or 3 of the 5 velocity-coherent features are well aligned with the filament and may represent genuine, fiber sub-structures. The core mass distribution has a peak at ~10Msun. They can be divided into 7 groups of cores, closely associated with ArTeMiS clumps. The projected separation between cores and the projected spacing between clumps are roughly consistent with the effective Jeans length in the filament and a physical scale of about 4 times the filament width, respectively, suggesting a bimodal filament fragmentation process. Despite being one order of magnitude denser and more massive than the B211/B213 filament, the NGC6334 filament has a similar density/velocity structure. The difference is that the cores in NGC6334 appear to be an order of magnitude denser and more massive than the cores in Taurus. This suggests that dense filaments may evolve and fragment in a similar manner in low- and high-mass star-forming regions, and that the filament paradigm may hold in the intermediate-mass (if not high-mass) star formation regime. △ Less

Submitted 10 November, 2019; v1 submitted 7 October, 2019; originally announced October 2019.

Comments: 21 pages, 18 figures, accepted for publication in A&A

Journal ref: A&A 632, A83 (2019)

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

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

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

Submitted 9 May, 2019; originally announced May 2019.

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

Journal ref: PASA, 36, e029 (2019)

The core and stellar mass functions in massive collapsing filaments

Authors: Evangelia Ntormousi, Patrick Hennebelle

Abstract: The connection between the pre-stellar core mass function (CMF) and the stellar initial mass function (IMF) lies at the heart of all star formation theories. In this paper, we study the earliest phases of star formation with a series of high-resolution numerical simulations that include the formation of sinks. In particular, we focus on the transition from cores to sinks within a massive molecular… ▽ More The connection between the pre-stellar core mass function (CMF) and the stellar initial mass function (IMF) lies at the heart of all star formation theories. In this paper, we study the earliest phases of star formation with a series of high-resolution numerical simulations that include the formation of sinks. In particular, we focus on the transition from cores to sinks within a massive molecular filament. We compare the CMF and IMF between magnetized and unmagnetized simulations, and between different resolutions. We find that selecting cores based on their kinematic virial parameter excludes collapsing objects because they host large velocity dispersions. Selecting only the thermally unstable magnetized cores, we observe that their mass-to-flux ratio spans almost two orders of magnitude for a given mass. We also see that, when magnetic fields are included, the CMF peaks at higher core mass values with respect to pure hydrodynamical simulations. Nonetheless, all models produce sink mass functions with a high-mass slope consistent with Salpeter. Finally, we examine the effects of resolution and find that, in isothermal simulations, even models with very high dynamical range fail to converge in the mass function. Our main conclusion is that, although the resulting CMFs and IMFs have similar slopes in all simulations, the cores have slightly different sizes and kinematical properties when a magnetic field is included. However, a core selection based on the mass-to-flux ratio alone is not enough to alter the shape of the CMF, if we do not take thermal stability into account. Finally, we conclude that extreme care should be given to resolution issues when studying sink formation with an isothermal equation of state. △ Less

Submitted 15 February, 2019; originally announced February 2019.

Comments: Accepted for publication by A&A

Journal ref: A&A 625, A82 (2019)

Probing accretion of ambient cloud material into the Taurus B211/B213 filament

Authors: Y. Shimajiri, Ph. Andre, P. Palmeirim, D. Arzoumanian, A. Bracco, V. Konyves, E. Ntormousi, B. Ladjelate

Abstract: Herschel observations have emphasized the role of molecular filaments in star formation. However, the origin and evolution of these filaments are not yet well understood, partly because of the lack of kinematic information. To examine whether the B211/B213 filament is accreting background gas due to its gravitational potential, we produced a toy accretion model and compared its predictions to the… ▽ More Herschel observations have emphasized the role of molecular filaments in star formation. However, the origin and evolution of these filaments are not yet well understood, partly because of the lack of kinematic information. To examine whether the B211/B213 filament is accreting background gas due to its gravitational potential, we produced a toy accretion model and compared its predictions to the 12CO(1--0) and 13CO(1--0) velocity patterns. We also examined the spatial distributions of Halpha, 857 GHz continuum, and HI emission to search for evidence of large-scale external effects. We estimated the depth of the cloud around the B211/B213 filament to be 0.3--0.7 pc under the assumption that the density of the gas is the same as the 13CO critical density. Compared to a linear extent of >10 pc in the plane of the sky, this suggests that the 3D morphology of the cloud is sheet-like. 12CO and 13CO PV diagrams perpendicular to the filament axis show that the emission from the gas surrounding B211/B213 is redshifted to the northeast of the filament and blueshifted to the southwest, respectively, and that the velocities of both components approach the filament velocity as the line of sight approaches the filament crest. The PV diagrams predicted by our accretion model are in good agreement with the observed 12CO and 13CO PV diagrams, supporting the scenario of mass accretion into the filament proposed by Palmeirim et al. Moreover, inspection of the distribution of the Halpha and 857 GHz emission in the Taurus-California-Perseus region suggests that the B211/B213 filament may have formed as a result of an expanding supershell generated by the Per OB2 association. Based on these results, we propose a scenario in which the B211/B213 filament was initially formed by large-scale compression of HI gas and then is now growing in mass due to the gravitational accretion of ambient cloud molecular gas. △ Less

Submitted 15 November, 2018; originally announced November 2018.

Comments: 16 pages, 17 figures, accepted in A&A

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

Magnetic fields in massive spirals: The role of feedback and initial conditions

Authors: Evangelia Ntormousi

Abstract: Magnetic fields play a very important role in the evolution of galaxies through their direct impact on star formation and stellar feedback-induced turbulence. However, their co-evolution with these processes has still not been thoroughly investigated, and the possible effect of the initial conditions is largely unknown. This letter presents the first results from a series of high-resolution numeri… ▽ More Magnetic fields play a very important role in the evolution of galaxies through their direct impact on star formation and stellar feedback-induced turbulence. However, their co-evolution with these processes has still not been thoroughly investigated, and the possible effect of the initial conditions is largely unknown. This letter presents the first results from a series of high-resolution numerical models, aimed at deciphering the effect of the initial conditions and of stellar feedback on the evolution of the galactic magnetic field in isolated, Milky-Way-like galaxies. The models start with an ordered, either poloidal or toroidal, magnetic field of varying strength, and are evolved with and without supernova feedback. They include a dark matter halo, a stellar and a gaseous disk, as well as the appropriate cooling and heating processes for the interstellar medium. Independently of the initial conditions, the galaxies develop a turbulent velocity field and a random magnetic field component in under 15 Myrs. Supernova feedback is extremely efficient in building a random magnetic field component up to large galactic heights. However, a random magnetic field emerges even in runs without feedback, which points to an inherent instability of the ordered component. Supernova feedback greatly affects the velocity field of the galaxy up to large galactic heights, and helps restructure the magnetic field up to 10 kpc above the disk, independently of the initial magnetic field morphology. On the other hand, the initial morphology of the magnetic field can accelerate the development of a random component at large heights. These effects have important implications for the study of the magnetic field evolution in galaxy simulations. △ Less

Submitted 19 October, 2018; originally announced October 2018.

Comments: A&A Letters, accepted

Journal ref: A&A 619, L5 (2018)

PASIPHAE: A high-Galactic-latitude, high-accuracy optopolarimetric survey

Authors: Konstantinos Tassis, Anamparambu N. Ramaprakash, Anthony C. S. Readhead, Stephen B. Potter, Ingunn K. Wehus, Georgia V. Panopoulou, Dmitry Blinov, Hans Kristian Eriksen, Brandon Hensley, Ata Karakci, John A. Kypriotakis, Siddharth Maharana, Evangelia Ntormousi, Vasiliki Pavlidou, Timothy J. Pearson, Raphael Skalidis

Abstract: PASIPHAE (the Polar-Areas Stellar Imaging in Polarization High-Accuracy Experiment) is an optopolarimetric survey aiming to measure the linear polarization from millions of stars, and use these to create a three-dimensional tomographic map of the magnetic field threading dust clouds within the Milky Way. This map will provide invaluable information for future CMB B-mode experiments searching for i… ▽ More PASIPHAE (the Polar-Areas Stellar Imaging in Polarization High-Accuracy Experiment) is an optopolarimetric survey aiming to measure the linear polarization from millions of stars, and use these to create a three-dimensional tomographic map of the magnetic field threading dust clouds within the Milky Way. This map will provide invaluable information for future CMB B-mode experiments searching for inflationary gravitational waves, providing unique information regarding line-of-sight integration effects. Optical polarization observations of a large number of stars at known distances, tracing the same dust that emits polarized microwaves, can map the magnetic field between them. The Gaia mission is measuring distances to a billion stars, providing an opportunity to produce a tomographic map of Galactic magnetic field directions, using optical polarization of starlight. Such a map will not only boost CMB polarization foreground removal, but it will also have a profound impact in a wide range of astrophysical research, including interstellar medium physics, high-energy astrophysics, and evolution of the Galaxy. Taking advantage of the novel technology implemented in our high-accuracy Wide-Area Linear Optical Polarimeters (WALOPs) currently under construction at IUCAA, India, we will engage in a large-scale optopolarimetric program that can meet this challenge: a survey of both northern and southern Galactic polar regions targeted by CMB experiments, covering over 10,000 square degrees, which will measure linear optical polarization of over 360 stars per square degree (over 3.5 million stars, a 1000-fold increase over the state of the art). The survey will be conducted concurrently from the South African Astronomical Observatory in Sutherland, South Africa in the southern hemisphere, and the Skinakas Observatory in Crete, Greece, in the north. △ Less

Submitted 12 October, 2018; originally announced October 2018.

Comments: 16 pages, 6 figures

The role of magnetic fields in the structure and interaction of supershells

Authors: Evangelia Ntormousi, Joanne R. Dawson, Patrick Hennebelle, Katharina Fierlinger

Abstract: Large-scale shocks formed by clustered feedback of young OB stars are considered an important source of mechanical energy for the ISM and a trigger of molecular cloud formation. Their interaction sites are locations where kinetic energy and magnetic fields are redistributed between ISM phases. In this work we study the effect of the magnetic field on the expansion and fragmentation of supershells… ▽ More Large-scale shocks formed by clustered feedback of young OB stars are considered an important source of mechanical energy for the ISM and a trigger of molecular cloud formation. Their interaction sites are locations where kinetic energy and magnetic fields are redistributed between ISM phases. In this work we study the effect of the magnetic field on the expansion and fragmentation of supershells and look for the signatures of supershell collisions on dense structures and on the kinetic and magnetic energy distribution of the ISM. We performed a series of high-resolution, three-dimensional simulations of colliding supershells. The shocks are created by time-dependent feedback and evolve in a diffuse turbulent environment that is either unmagnetized or has different initial magnetic field configurations. In the hydrodynamical situation, the expansion law of the superbubbles is consistent with the radius-time relation that is theoretically predicted for wind-blown bubbles. The supershells fragment over their entire surface into small dense clumps that carry more than half of the total kinetic energy in the volume. However, this is not the case when a magnetic field is introduced, either in the direction of the collision or perpendicular to the collision. In magnetized situations, the shell surfaces are more stable to dynamical instabilities. When the magnetic field opposes the collision, the expansion law of the supershells also becomes significantly flatter than in the hydrodynamical case. Although a two-phase medium arises in all cases, in the MHD simulations the cold phase is limited to lower densities. △ Less

Submitted 13 January, 2017; originally announced January 2017.

Comments: to appear in A&A Section 6: Interstellar and circumstellar matter

Journal ref: A&A 599, A94 (2017)

The effect of ambipolar diffusion on low-density molecular ISM filaments

Authors: Evangelia Ntormousi, Patrick Hennebelle, Philippe André, Jacques Masson

Abstract: The filamentary structure of the molecular interstellar medium and the potential link of this morphology to star formation have been brought into focus recently by high resolution observational surveys. An especially puzzling matter is that local interstellar filaments appear to have the same thickness, independent of their column density. This requires a theoretical understanding of their formati… ▽ More The filamentary structure of the molecular interstellar medium and the potential link of this morphology to star formation have been brought into focus recently by high resolution observational surveys. An especially puzzling matter is that local interstellar filaments appear to have the same thickness, independent of their column density. This requires a theoretical understanding of their formation process and the physics that governs their evolution. In this work we explore a scenario in which filaments are dissipative structures of the large-scale interstellar turbulence cascade and ion-neutral friction (also called ambipolar diffusion) is affecting their sizes by preventing small-scale compressions. We employ high-resolution, 3D MHD simulations, performed with the grid code RAMSES, to investigate non-ideal MHD turbulence as a filament formation mechanism. We focus the analysis on the mass and thickness distributions of the resulting filamentary structures. Simulations of both driven and decaying MHD turbulence show that the morphologies of the density and the magnetic field are different when ambipolar diffusion is included in the models. In particular, the densest structures are broader and more massive as an effect of ion-neutral friction and the power spectra of both the velocity and the density steepen at a smaller wavenumber. The comparison between ideal and non-ideal MHD simulations shows that ambipolar diffusion causes a shift of the filament thickness distribution towards higher values. However, none of the distributions exhibit the pronounced peak found in the observed local filaments. Limitations in dynamical range and the absence of self-gravity in these numerical experiments do not allow us to conclude at this time whether this is due to the different filament selection or due to the physics inherent of the filament formation. △ Less

Submitted 15 February, 2016; originally announced February 2016.

Comments: A&A accepted

Stellar feedback efficiencies: supernovae versus stellar winds

Authors: Katharina M. Fierlinger, Andreas Burkert, Evangelia Ntormousi, Peter Fierlinger, Marc Schartmann, Alessandro Ballone, Martin G. H. Krause, Roland Diehl

Abstract: Stellar winds and supernova (SN) explosions of massive stars ("stellar feedback") create bubbles in the interstellar medium (ISM) and insert newly produced heavy elements and kinetic energy into their surroundings, possibly driving turbulence. Most of this energy is thermalized and immediately removed from the ISM by radiative cooling. The rest is available for driving ISM dynamics. In this work w… ▽ More Stellar winds and supernova (SN) explosions of massive stars ("stellar feedback") create bubbles in the interstellar medium (ISM) and insert newly produced heavy elements and kinetic energy into their surroundings, possibly driving turbulence. Most of this energy is thermalized and immediately removed from the ISM by radiative cooling. The rest is available for driving ISM dynamics. In this work we estimate the amount of feedback energy retained as kinetic energy when the bubble walls have decelerated to the sound speed of the ambient medium. We show that the feedback of the most massive star outweighs the feedback from less massive stars. For a giant molecular cloud (GMC) mass of 1e5 solar masses (as e.g. found in the Orion GMCs) and a star formation efficiency of 8% the initial mass function predicts a most massive star of approximately 60 solar masses. For this stellar evolution model we test the dependence of the retained kinetic energy of the cold GMC gas on the inclusion of stellar winds. In our model winds insert 2.34 times the energy of a SN and create stellar wind bubbles serving as pressure reservoirs. We find that during the pressure driven phases of the bubble evolution radiative losses peak near the contact discontinuity (CD), and thus, the retained energy depends critically on the scales of the mixing processes across the CD. Taking into account the winds of massive stars increases the amount of kinetic energy deposited in the cold ISM from 0.1% to a few percent of the feedback energy. △ Less

Submitted 18 November, 2015; v1 submitted 16 November, 2015; originally announced November 2015.

Comments: The manuscript MN-15-2374-MJ.R2 has been accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal on 2015 November 16

Shell instability of a collapsing dense core

Authors: Evangelia Ntormousi, Patrick Hennebelle

Abstract: Understanding the formation of binary and multiple stellar systems largely comes down to studying the circumstances for the fragmentation of a condensing core during the first stages of the collapse. However, the probability of fragmentation and the number of fragments seem to be determined to a large degree by the initial conditions. In this work we study the fate of the linear perturbations of a… ▽ More Understanding the formation of binary and multiple stellar systems largely comes down to studying the circumstances for the fragmentation of a condensing core during the first stages of the collapse. However, the probability of fragmentation and the number of fragments seem to be determined to a large degree by the initial conditions. In this work we study the fate of the linear perturbations of a homogeneous gas sphere both analytically and numerically. In particular, we investigate the stability of the well-known homologous solution that describes the collapse of a uniform spherical cloud. The difficulty of the mathematical singularity in the perturbation equations is surpassed here by explicitly introducing a weak shock next to the sonic point. In parallel, we perform adaptive mesh refinement (AMR) numerical simulations of the linear stages of the collapse and compared the growth rates obtained by each method. With this combination of analytical and numerical tools, we explore the behavior of both spherically symmetric and non-axisymmetric perturbations. The numerical experiments provide the linear growth rates as a function of the core's initial virial parameter and as a function of the azimuthal wave number of the perturbation. The overlapping regime of the numerical experiments and the analytical predictions is the situation of a cold and large cloud, and in this regime the analytically calculated growth rates agree very well with the ones obtained from the simulations. The use of a weak shock as part of the perturbation allows us to find a physically acceptable solution to the equations for a continuous range of growth rates. The numerical simulations agree very well with the analytical prediction for the most unstable cores, while they impose a limit of a virial parameter of 0.1 for core fragmentation in the absence of rotation. △ Less

Submitted 12 November, 2014; originally announced November 2014.

Comments: Accepted by A&A

Journal ref: A&A 574, A130 (2015)

A Young GMC Formed at the Interface of Two Colliding Supershells: Observations Meet Simulations

Authors: J. R. Dawson, E. Ntormousi, Y. Fukui, T. Hayakawa, K. Fierlinger

Abstract: Dense, star-forming gas is believed to form at the stagnation points of large-scale ISM flows, but observational examples of this process in action are rare. We here present a giant molecular cloud (GMC) sandwiched between two colliding Milky Way supershells, which we argue shows strong evidence of having formed from material accumulated at the collision zone. Combining 12CO, 13CO and C18O(J=1-0)… ▽ More Dense, star-forming gas is believed to form at the stagnation points of large-scale ISM flows, but observational examples of this process in action are rare. We here present a giant molecular cloud (GMC) sandwiched between two colliding Milky Way supershells, which we argue shows strong evidence of having formed from material accumulated at the collision zone. Combining 12CO, 13CO and C18O(J=1-0) data with new high-resolution, 3D hydrodynamical simulations of colliding supershells, we discuss the origin and nature of the GMC (G288.5+1.5), favoring a scenario in which the cloud was partially seeded by pre-existing denser material, but assembled into its current form by the action of the shells. This assembly includes the production of some new molecular gas. The GMC is well interpreted as non-self-gravitating, despite its high mass (MH2 ~ 1.7 x 10^5 Msol), and is likely pressure confined by the colliding flows, implying that self-gravity was not a necessary ingredient for its formation. Much of the molecular gas is relatively diffuse, and the cloud as a whole shows little evidence of star formation activity, supporting a scenario in which it is young and recently formed. Drip-like formations along its lower edge may be explained by fluid dynamical instabilities in the cooled gas. △ Less

Submitted 11 November, 2014; originally announced November 2014.

Comments: 13 pages, 9 figures, accepted for publication in ApJ

Age, size, and position of H ii regions in the Galaxy. Expansion of ionized gas in turbulent molecular clouds

Authors: P. Tremblin, L. D. Anderson, P. Didelon, A. C. Raga, V. Minier, E. Ntormousi, A. Pettitt, C. Pinto, M. Samal, N. Schneider, A. Zavagno

Abstract: This work aims at improving the current understanding of the interaction between H ii regions and turbulent molecular clouds. We propose a new method to determine the age of a large sample of OB associations by investigating the development of their associated H ii regions in the surrounding turbulent medium. Using analytical solutions, one-dimensional (1D), and three-dimensional (3D) simulations,… ▽ More This work aims at improving the current understanding of the interaction between H ii regions and turbulent molecular clouds. We propose a new method to determine the age of a large sample of OB associations by investigating the development of their associated H ii regions in the surrounding turbulent medium. Using analytical solutions, one-dimensional (1D), and three-dimensional (3D) simulations, we constrained the expansion of the ionized bubble depending on the turbulent level of the parent molecular cloud. A grid of 1D simulations was then computed in order to build isochrone curves for H ii regions in a pressure-size diagram. This grid of models allowed to date large sample of OB associations and was used on the H ii Region Discovery Survey (HRDS). Analytical solutions and numerical simulations showed that the expansion of H ii regions is slowed down by the turbulence up to the point where the pressure of the ionized gas is in a quasi-equilibrium with the turbulent ram pressure. Based on this result, we built a grid of 1D models of the expansion of H ii regions in a profile based on Larson laws. The 3D turbulence is taken into account by an effective 1D temperature profile. The ages estimated by the isochrones of this grid agree well with literature values of well-known regions such as Rosette, RCW 36, RCW 79, and M16. We thus propose that this method can be used to give ages of young OB associations through the Galaxy such as the HRDS survey and also in nearby extra-galactic sources. △ Less

Submitted 6 June, 2014; originally announced June 2014.

Comments: Accepted in A&A

Journal ref: A&A 568, A4 (2014)

Delayed metal recycling in galaxies: the inefficiency of cold gas enrichment in colliding supershell simulations

Authors: Evangelia Ntormousi, Andreas Burkert

Abstract: The fate of metals ejected by young OB associations into the Interstellar Medium (ISM) is investigated numerically. In particular, we study the enrichment of the cold gas phase, which is the material that forms molecular clouds. Following previous work, the expansion and collision of two supershells in a diffuse ISM is simulated, in this case also introducing an advected quantity which represents… ▽ More The fate of metals ejected by young OB associations into the Interstellar Medium (ISM) is investigated numerically. In particular, we study the enrichment of the cold gas phase, which is the material that forms molecular clouds. Following previous work, the expansion and collision of two supershells in a diffuse ISM is simulated, in this case also introducing an advected quantity which represents the metals expelled by the young stars. We adopt the simplest possible approach, not differentiating between metals coming from stellar winds and those coming from supernovae. Even though the hot, diffuse phase of the ISM receives a significant amount of metals from the stars, the cold phase is efficiently shielded, with very little metal enrichment. Significant enrichment of the cold ISM will therefore be delayed by at least the cooling time of this hot phase. No variations in cloud metallicity with distance from the OB association or with direction are found, which means that the shell collision does little to enhance the metallicity of the cold clumps. We conclude that the stellar generation that forms out of molecular structures, triggered by shell collisions cannot be significantly enriched. △ Less

Submitted 8 November, 2011; originally announced November 2011.

Comments: 11 pages, 16 figures, submitted to the Astrophysical Journal

Formation of Cold Filamentary Structure from Wind Blown Superbubbles

Authors: Evangelia Ntormousi, Andreas Burkert, Katharina Fierlinger, Fabian Heitsch

Abstract: The expansion and collision of two wind-blown superbubbles is investigated numerically. Our models go beyond previous simulations of molecular cloud formation from converging gas flows by exploring this process with realistic flow parameters, sizes and timescales. The superbubbles are blown by time-dependent winds and supernova explosions, calculated from population synthesis models. They expand i… ▽ More The expansion and collision of two wind-blown superbubbles is investigated numerically. Our models go beyond previous simulations of molecular cloud formation from converging gas flows by exploring this process with realistic flow parameters, sizes and timescales. The superbubbles are blown by time-dependent winds and supernova explosions, calculated from population synthesis models. They expand into a uniform or turbulent diffuse medium. We find that dense, cold gas clumps and filaments form naturally in the compressed collision zone of the two superbubbles. Their shapes resemble the elongated, irregular structure of observed cold, molecular gas filaments and clumps. At the end of the simulations, between 65 and 80 percent of the total gas mass in our simulation box is contained in these structures. The clumps are found in a variety of physical states, ranging from pressure equilibrium with the surrounding medium to highly under-pressured clumps with large irregular internal motions and structures which are rotationally supported. △ Less

Submitted 26 November, 2010; originally announced November 2010.

Comments: Submitted to ApJ

Hot gas haloes around disc galaxies: O VII column densities from galaxy formation simulations

Authors: Evangelia Ntormousi, Jesper Sommer-Larsen

Abstract: Numerical models of disc galaxy formation predict the existence of extended, hot ( T~10^6 K) gas haloes around present day spirals. The X-ray luminosity of these haloes is predicted to increase strongly with galaxy mass. However, searches for their X-ray emission have not been successful so far. We calculate the all sky O VII column density distributions for the haloes of three Milky Way like disc… ▽ More Numerical models of disc galaxy formation predict the existence of extended, hot ( T~10^6 K) gas haloes around present day spirals. The X-ray luminosity of these haloes is predicted to increase strongly with galaxy mass. However, searches for their X-ray emission have not been successful so far. We calculate the all sky O VII column density distributions for the haloes of three Milky Way like disc galaxies, resulting from cosmological high-resolution, N-body/gasdynamical simulations. We perform calculations both including the disc gas and without it, so the disc contribution to the column density is quantified. It is found that the column densities estimated for Milky Way-like galaxies are just below the observational upper limit, making a test of the hot halo paradigm likely within observational reach. △ Less

Submitted 23 July, 2010; originally announced July 2010.

Comments: accepted for publication in MNRAS