Nearby Massive Stars Research Articles

Chemistry in Externally FUV-irradiated Disks in the Outskirts of the Orion Nebula Cluster

Most stars are born in stellar clusters, and their protoplanetary disks, which are the birthplaces of planets, can, therefore, be affected by the radiation of nearby massive stars. However, little is known about the chemistry of externally irradiated disks, including whether or not their properties are similar to the so-far better-studied isolated disks. Motivated by this question, we present ALMA Band 6 observations of two irradiated Class II protoplanetary disks in the outskirts of the Orion Nebula Cluster to explore the chemical composition of disks exposed to (external) far-ultraviolet (FUV) radiation fields: the 216-0939 disk and the binary system 253-1536A/B, which are exposed to radiation fields of 102–103 times the average interstellar radiation field. We detect lines from CO isotopologues, HCN, H2CO, and C2H toward both protoplanetary disks. Based on the observed disk-integrated line fluxes and flux ratios, we do not find significant differences between isolated and irradiated disks. The observed differences seem to be more closely related to the different stellar masses than to the external radiation field. This suggests that these disks are far enough away from the massive Trapezium stars, that their chemistry is no longer affected by external FUV radiation. Additional observations toward lower-mass disks and disks closer to the massive Trapezium stars are required to elucidate the level of external radiation required to make an impact on the chemistry of planet formation in different kinds of disks.

Atacama Large Aperture Submillimeter Telescope (AtLAST) science: Our Galaxy.

As we learn more about the multi-scale interstellar medium (ISM) of our Galaxy, we develop a greater understanding for the complex relationships between the large-scale diffuse gas and dust in Giant Molecular Clouds (GMCs), how it moves, how it is affected by the nearby massive stars, and which portions of those GMCs eventually collapse into star forming regions. The complex interactions of those gas, dust and stellar populations form what has come to be known as the ecology of our Galaxy. Because we are deeply embedded in the plane of our Galaxy, it takes up a significant fraction of the sky, with complex dust lanes scattered throughout the optically recognizable bands of the Milky Way. These bands become bright at (sub-)millimetre wavelengths, where we can study dust thermal emission and the chemical and kinematic signatures of the gas. To properly study such large-scale environments, requires deep, large area surveys that are not possible with current facilities. Moreover, where stars form, so too do planetary systems, growing from the dust and gas in circumstellar discs, to planets and planetesimal belts. Understanding the evolution of these belts requires deep imaging capable of studying belts around young stellar objects to Kuiper belt analogues around the nearest stars. Here we present a plan for observing the Galactic Plane and circumstellar environments to quantify the physical structure, the magnetic fields, the dynamics, chemistry, star formation, and planetary system evolution of the galaxy in which we live with AtLAST; a concept for a new, 50m single-dish sub-mm telescope with a large field of view which is the only type of facility that will allow us to observe our Galaxy deeply and widely enough to make a leap forward in our understanding of our local ecology.

HCN emission from translucent gas and UV-illuminated cloud edges revealed by wide-field IRAM 30 m maps of the Orion B GMC

Context. Massive stars form within dense clumps inside giant molecular clouds (GMCs). Finding appropriate chemical tracers of the dense gas (n(H2) > several 104 cm−3 or AV > 8 mag) and linking their line luminosity with the star formation rate is of critical importance. Aims. Our aim is to determine the origin and physical conditions of the HCN-emitting gas and study their relation to those of other molecules. Methods. In the context of the IRAM 30m ORION-B large program, we present 5 deg2 (~250 pc2) HCN, HNC, HCO+, and CO J =1–0 maps of the Orion B GMC, complemented with existing wide-field [C I] 492 GHz maps, as well as new pointed observations of rotationally excited HCN, HNC, H13CN, and HN13C lines. We compare the observed HCN line intensities with radiative transfer models including line overlap effects and electron excitation. Furthermore, we study the HCN/HNC isomeric abundance ratio with updated photochemical models. Results. We spectroscopically resolve the HCN J = 1–0 hyperfine structure (HFS) components (and partially resolved J = 2−1 and 3−2 components). We detect anomalous HFS line intensity (and line width) ratios almost everywhere in the cloud. About 70% of the total HCN J = 1−0 luminosity, L′(HCN J = 1−0) = 110 K km s−1 pc−2, arises from AV < 8 mag. The HCN/CO J = 1−0 line intensity ratio, widely used as a tracer of the dense gas fraction, shows a bimodal behavior with an inflection point at AV < 3 mag typical of translucent gas and illuminated cloud edges. We find that most of the HCN J = 1−0 emission arises from extended gas with n(H2) < 104 cm−3, and even lower density gas if the ionization fraction is χe ≥ 10−5 and electron excitation dominates. This result contrasts with the prevailing view of HCN J = 1−0 emission as a tracer of dense gas and explains the low-AV branch of the HCN/CO J = 1−0 intensity ratio distribution. Indeed, the highest HCN/CO ratios (~ 0.1) at AV < 3 mag correspond to regions of high [C I] 492 GHz/CO J = 1−0 intensity ratios (>1) characteristic of low-density photodissociation regions. The low surface brightness (≲ 1 K km s−1) and extended HCN and HCO+ J = 1−0 emission scale with IFIR – a proxy of the stellar far-ultraviolet (FUV) radiation field – in a similar way. Together with CO J = 1−0, these lines respond to increasing IFIR up to G0 ≃ 20. On the other hand, the bright HCN J = 1−0 emission (> 6 K km s−1) from dense gas in star-forming clumps weakly responds to IFIR once the FUV field becomes too intense (G0 > 1500). In contrast, HNC J = 1−0 and [C I] 492 GHz lines weakly respond to IFIR for all G0. The different power law scalings (produced by different chemistries, densities, and line excitation regimes) in a single but spatially resolved GMC resemble the variety of Kennicutt-Schmidt law indexes found in galaxy averages. Conclusions. Given the widespread and extended nature of the [C I] 492 GHz emission, as well as its spatial correlation with that of HCO+, HCN, and 13CO J = 1−0 lines (in this order), we argue that the edges of GMCs are porous to FUV radiation from nearby massive stars. Enhanced FUV radiation favors the formation and excitation of HCN on large scales, not only in dense star-forming clumps, and it leads to a relatively low value of the dense gas mass to total luminosity ratio, α (HCN) = 29 M⊙/(K km s−1pc2) in Orion B. As a corollary for extragalactic studies, we conclude that high HCN/CO J = 1−0 line intensity ratios do not always imply the presence of dense gas, which may be better traced by HNC than by HCN.

Dependence of Chemical Abundance on the Cosmic-Ray Ionization Rate in IC 348

Ions (e.g., , H2O+) have been used extensively to quantify the cosmic-ray ionization rate (CRIR) in diffuse sightlines. However, measurements of the CRIR in low-to-intermediate density gas environments are rare, especially when background stars are absent. In this work, we combine molecular line observations of CO, OH, CH, and HCO+ in the star-forming cloud IC 348, and chemical models to constrain the value of the CRIR and study the response of the distribution of chemical abundances. The cloud boundary is found to have an A V of approximately 4 mag. From the interior to the exterior of the cloud, the observed 13CO line intensities drop by an order of magnitude. The calculated average abundance of 12CO (assuming 12C/13C=65) is (1.2 ± 0.9) × 10−4, which decreases by a factor of 6 from the interior to the outside regions. The average abundance of CH (3.3 ± 0.7 × 10−8) is in good agreement with previous findings in diffuse and translucent clouds (A V < 5 mag). However, we did not find a decline in CH abundance in regions of high extinction (A V ≃ 8 mag) as previously reported in Taurus. By comparing the observed molecular abundances and chemical models, we find a decreasing trend of the CRIR as A V increases. The inferred CRIR of ζ cr = (4.7 ± 1.5) × 10−16 s−1 at low A V is consistent with measurements toward two nearby massive stars.

The growth and migration of massive planets under the influence of external photoevaporation

ABSTRACT The formation of gas giant planets must occur during the first few Myr of a star’s lifetime, when the protoplanetary disc still contains sufficient gas to be accreted on to the planetary core. The majority of protoplanetary discs are exposed to strong ultraviolet irradiation from nearby massive stars, which drives winds and depletes the mass budget for planet formation. It remains unclear to what degree external photoevaporation affects the formation of massive planets. In this work, we present a simple one dimensional model for the growth and migration of a massive planet under the influence of external FUV fields. We find that even moderate FUV fluxes $F_\mathrm{FUV}\gtrsim 100 \, G_0$ have a strong influence on planet mass and migration. By decreasing the local surface density and shutting off accretion on to the planet, external irradiation suppresses planet masses and halts migration early. The distribution of typical stellar birth environments can therefore produce an anticorrelation between semi-major axis and planet mass, which may explain the apparent decrease in planet occurrence rates at orbital periods Porb ≳ 103 d. Even moderate fluxes FFUV strongly suppress giant planet formation and inward migration for any initial semi-major axis if the stellar host mass $M_*\lesssim 0.5\, {\rm M}_\odot$, consistent with findings that massive planet occurrence is much lower around such stars. The outcomes of our prescription for external disc depletion show significant differences to the current approximation adopted in state-of-the-art population synthesis models, motivating future careful treatment of this important process.

An extended scattered light disk around AT Pyx

Aims. To understand how the multitude of planetary systems that have been discovered come to be, we need to study systems at different evolutionary stages, with different central stars but also in different environments. The most challenging environment for planet formation may be the harsh UV radiation field of nearby massive stars which quickly erodes disks by external photo-evaporation. We observed the AT Pyx system, located in the head of a cometary globule in the Gum Nebula, to search for signs of ongoing planet formation. Methods. We used the extreme adaptive optics imager VLT/SPHERE in Dual Beam Polarization Imaging Mode in H-band as well as in IRDIFS Extended mode (K12-band imaging and Y-H integral field spectroscopy) to observe AT Pyx in polarized light and total intensity. Additionally, we employed VLT/NACO to observe the system in the L-band. Results. We resolve the disk around AT Pyx for the first time in scattered light across multiple wavelengths in polarized light and total intensity. We find an extended (≥126 au) disk, with an intermediate inclination of between 35° and 42°. The disk shows a complex substructure and we identify two or possibly three spiral-like features. Depending on the precise geometry of the disk (which we cannot unambiguously infer from our data), the disk may be eccentric with an eccentricity of ~0.16 or partially self-shadowed. The spiral features and possible eccentricity are both consistent with signatures of an embedded gas giant planet with a mass of ~1 MJup. Our own observations can rule out brown dwarf companions embedded in the resolved disk, but are nevertheless not sensitive enough to confirm or rule out the presence of a gas giant. Conclusions. AT Pyx is the first disk to be spatially resolved in a cometary globule in the Gum Nebula. By comparison with disks in the Orion Nebula Cluster we note that the extension of the disk may be exceptional for this environment if the external UV radiation field is indeed comparable to other cometary globules in the region. The signposts of ongoing planet formation are intriguing and need to be followed up with either higher sensitivity or at different wavelengths.

Diagnosing Triggered Star Formation in the Galactic H ii region Sh 2-142

Abstract Stars are formed by gravitational collapse, spontaneously or, in some cases under the constructive influence of nearby massive stars, out of molecular cloud cores. Here we present an observational diagnosis of such triggered formation processes in the prominent H ii region Sh 2-142, which is associated with the young star cluster NGC 7380, and with some bright-rimmed clouds as the signpost of photoionization of molecular cloud surfaces. Using near- (2MASS) and mid-infrared (WISE) colors, we identified candidate young stars at different evolutionary stages, including embedded infrared sources having spectral energy distributions indicative of active accretion. We have also used data from our optical observations to be used in SEDs, and from Gaia EDR3 to study the kinematics of young objects. With this young stellar sample, together with the latest CO line emission data (spectral resolution ∼0.16 km s−1, sensitivity ∼0.5 K), a positional and aging sequence relative to the neighboring cloud complex, and to the bright-rimmed clouds, is inferred. The propagating stellar birth may be responsible, at least partially, for the formation of the cluster a few million years ago, and for the ongoing activity now witnessed in the cloud complex.

W51 North: A protocluster emerging out of a thermally inhibited fragmenting cloud

Context. The fragmentation process in massive star-forming regions is one of the contemporary problems in astrophysics, and several physical processes have been proposed to control the fragmentation including turbulence, magnetic field, rotation, stellar feedback, and gravity. However, the fragmentation process has been poorly studied at small spatial scales well below 1000 AU. Aims. W51 IRS2 is a well-known massive star-forming region whose fragmentation properties have not been thoroughly investigated yet. We aim to use the Atacama Large Millimeter and Submillimeter Array (ALMA) high angular resolution data to identify the fragments in W51 IRS2 and to study the fragmentation properties on a spatial scale of 0.04′′ (200 AU). Methods. We used ALMA data of W51 IRS2 from three projects, which give an angular resolution of 0.028′′ (144 AU) at millimeter wavelengths. The continuum images reveal a significant substructure in an east-west ridge, where both W51 North and W51d2 are embedded. A spectral index map has been generated from the 3 and 1.3 mm high-resolution continuum images. We identified compact fragments by using uv-range constrained 1.3 mm continuum data. A mean surface density of companions (MSDC) analysis has been performed to study the separations between fragments. Results. A total number of 33 continuum sources are identified and 29 out of them are defined as fragments in the surveyed region.The MSDC analysis reveals two breaks corresponding to spatial scales of 1845 and 7346 AU, indicative of a two-level clustering phenomenon, along with a linear regime below 1845 AU, mostly associated with W51 North, whose slope is consistent with the slope for the clustering regime of other cluster-like regions in the Galaxy. Conclusions. The typical masses and separations of the fragments as well as the relation between the density and number of fragments can be explained through a thermal Jeans process operating at high temperatures of 200–400 K, consistent with previous measurements of the temperature in the region, and produced by the nearby massive stars. Therefore, although W51 IRS2 seems to be undergoing a thermally inhibited fragmentation phase, this does not seem to prevent the formation of a protocluster associated with W51 North.

Small Protoplanetary Disks in the Orion Nebula Cluster and OMC1 with ALMA

Abstract The Orion Nebula Cluster (ONC) is the nearest dense star-forming region at ∼400 pc away, making it an ideal target to study the impact of high stellar density and proximity to massive stars (the Trapezium) on protoplanetary disk evolution. The OMC1 molecular cloud is a region of high extinction situated behind the Trapezium in which actively forming stars are shielded from the Trapezium’s strong radiation. In this work, we survey disks at high resolution with Atacama Large Millimeter/submillimeter Array at three wavelengths with resolutions of 0.″095 (3 mm; Band 3), 0.″048 (1.3 mm; Band 6), and 0.″030 (0.85 mm; Band 7) centered on radio Source I. We detect 127 sources, including 15 new sources that have not previously been detected at any wavelength. 72 sources are spatially resolved at 3 mm, with sizes from ∼8–100 au. We classify 76 infrared-detected sources as foreground ONC disks and the remainder as embedded OMC1 disks. The two samples have similar disk sizes, but the OMC1 sources have a dense and centrally concentrated spatial distribution, indicating they may constitute a spatially distinct subcluster. We find smaller disk sizes and a lack of large (&gt;75 au) disks in both our samples compared to other nearby star-forming regions, indicating that environmental disk truncation processes are significant. While photoevaporation from nearby massive Trapezium stars may account for the smaller disks in the ONC, the embedded sources in OMC1 are hidden from this radiation and thus must truncated by some other mechanism, possibly dynamical truncation or accretion-driven contraction.

A slim disc approach to external photoevaporation of discs

ABSTRACT The photoevaporation of protoplanetary discs by nearby massive stars present in their birth cluster plays a vital role in their evolution. Previous modelling assumes that the disc behaves like a classical Keplerian accretion disc out to a radius where the photoevaporative outflow is launched. There is then an abrupt change in the angular velocity profile, and the outflow is modelled by forcing the fluid parcels to conserve their specific angular momenta. Instead, we model externally photoevaporating discs using the slim disc formalism. The slim disc approach self-consistently includes the advection of radial and angular momentum as well as angular momentum redistribution by internal viscous torques. Our resulting models produce a smooth transition from a rotationally supported Keplerian disc to a photoevaporative driven outflow, where this transition typically occurs over ∼4–5 scale heights. The penetration of ultraviolet photons predominately sets the radius of the transition and the viscosity’s strength plays a minor role. By studying the entrainment of dust particles in the outflow, we find a rapid change in the dust size and surface density distribution in the transition region due to the steep gas density gradients present. This rapid change in the dust properties leaves a potentially observable signature in the continuum spectral index of the disc at mm wavelengths. Using the slim disc formalism in future evolutionary calculations will reveal how both the gas and dust evolve in their outer regions and the observable imprints of the external photoevaporation process.

Probing protoplanetary disk evolution in the Chamaeleon II region

Context. Characterizing the evolution of protoplanetary disks is necessary to improve our understanding of planet formation. Constraints on both dust and gas are needed to determine the dominant disk dissipation mechanisms. Aims. We aim to compare the disk dust masses in the Chamaeleon II (Cha II) star-forming region with other regions with ages between 1 and 10 Myr. Methods. We use ALMA band 6 observations (1.3 mm) to survey 29 protoplanetary disks in Cha II. Dust mass estimates are derived from the continuum data. Results. Out of our initial sample of 29 disks, we detect 22 sources in the continuum, 10 in 12CO, 3 in 13CO, and none in C18O (J = 2−1). Additionally, we detect two companion candidates in the continuum and 12CO emission. Most disk dust masses are lower than 10 M⊕, assuming thermal emission from optically thin dust. Including non-detections, we derive a median dust mass of 4.5 ± 1.5 M⊕ from survival analysis. We compare consistent estimations of the distributions of the disk dust mass and the disk-to-stellar mass ratios in Cha II with six other low mass and isolated star-forming regions in the age range of 1–10 Myr: Upper Sco, CrA, IC 348, Cha I, Lupus, and Taurus. When comparing the dust-to-stellar mass ratio, we find that the masses of disks in Cha II are statistically different from those in Upper Sco and Taurus, and we confirm that disks in Upper Sco, the oldest region of the sample, are statistically less massive than in all other regions. Performing a second statistical test of the dust mass distributions from similar mass bins, we find no statistical differences between these regions and Cha II. Conclusions. We interpret these trends, most simply, as a sign of decline in the disk dust masses with time or dust evolution. Different global initial conditions in star-forming regions may also play a role, but their impact on the properties of a disk population is difficult to isolate in star-forming regions lacking nearby massive stars.

Cosmic symplectite recorded irradiation by nearby massive stars in the solar system's parent molecular cloud

The Sun’s astrophysical birth environment affected the formation and composition of the Solar System. Primitive meteorites display mass-independent oxygen isotope anomalies that were likely caused by ultraviolet (UV) photochemistry of CO gas-phase molecules, either (i) in the outer solar nebula by light from the young Sun or (ii) in the parent molecular cloud by light from nearby stars. However, measurements of oxygen isotopes alone cannot unambiguously constrain the UV spectrum of the source responsible for the photochemistry. Sulfur, with four stable isotopes, can be used as a more direct probe of the astrophysical environment of mass-independent photochemistry. Here, we report the in situ isotopic analysis of paired oxygen and sulfur isotope systematics in cosmic symplectite (COS), magnetite-pentlandite intergrowths, in the primitive ungrouped carbonaceous chondrite Acfer 094. We show that COS grains contain mass-independent sulfur isotope anomalies (weighted means of Δ33S = +3.84 ± 0.72‰ and Δ36S = −6.05 ± 2.25‰, 2SE) consistent with H2S photochemistry by UV from massive O and B stars close to the Solar System’s parent molecular cloud, and inconsistent with UV from the protosun. The presence of coupled mass-independent sulfur and oxygen (Δ17O = 86 ± 6‰, 2SE) isotope anomalies in COS imply that these anomalies originated in the same astrophysical environment. We propose that this environment is the photodissociation region (PDR) of the Solar System’s parent molecular cloud, where nearby massive stars irradiated the edge of the cloud. We conclude that the Sun’s stellar neighbors, likely O and B stars in a massive-star-forming region, affected the composition of the Solar System’s primordial building blocks.

A Spectroscopic View on Cosmic PAH Emission

ConspectusPolycyclic aromatic hydrocarbon molecules (PAHs) are ubiquitously present at high abundances in the Universe. They are detected through their infrared (IR) fluorescence UV pumped by nearby massive stars. Hence, their infrared emission is used to determine the star formation rate in galaxies, one of the key indicators for understanding the evolution of galaxies. Together with fullerenes, PAHs are the largest molecules found in space. They significantly partake in a variety of physical and chemical processes in space, influencing star and planet formation as well as galaxy evolution.Since the IR features from PAHs originate from chemical bonds involving only nearest neighbor atoms, they have only a weak dependence on the size and structure of the molecule, and it is therefore not possible to identify the individual PAH molecules that make up the cosmic PAH family. This strongly hampers the interpretation of their astronomical fingerprints. Despite the lack of identification, constraints can be set on the characteristics of the cosmic PAH family thanks to a joint effort of astronomers, physicists, and chemists.This Account presents the spectroscopic properties of the cosmic PAH emission as well as the intrinsic spectroscopic properties of PAHs and astronomical modeling of the PAH evolution required for the interpretation of the cosmic PAH characteristics. We discuss the observed spectral signatures tracing PAH properties such as charge, size, and structure and highlight the related challenges. We discuss the recent success of anharmonic calculations of PAH infrared absorption and emission spectra and outline the path forward. Finally, we illustrate the importance of models on PAH processing for the interpretation of the astronomical data in terms of the charge balance and PAH destruction.Throughout this Account, we emphasize that huge progress is on the horizon on the astronomical front. Indeed, the world is eagerly awaiting the launch of the James Webb Space Telescope (JWST). With its incredible improvement in spatial resolution, combined with its complete spectral coverage of the PAH infrared emission bands at medium spectral resolution and superb sensitivity, the JWST will revolutionize PAH research. Previous observations could only present spectra averaged over regions with vastly different properties, thus greatly confusing their interpretation. The amazing spatial resolution of JWST will disentangle these different regions. This will allow us to quantify precisely how PAHs are modified by the physical conditions of their host environment and thus trace how PAHs evolve across space. However, this will only be achieved when the necessary (and still missing) fundamental properties of PAHs, outlined in this Account, are known. We strongly encourage you to join this effort.

Warm millimetre dust in protoplanetary discs near massive stars

ABSTRACT Dust plays a key role in the formation of planets and its emission also provides one of our most accessible views of protoplanetary discs. If set by radiative equilibrium with the central star, the temperature of dust in the disc plateaus at around 10–20 K in the outer regions. However, sufficiently nearby massive stars can heat the outer disc to substantially higher temperatures. In this paper, we study the radiative equilibrium temperature of discs in the presence of massive external sources and gauge the effect that it has on millimetre dust mass estimates. Since millimetre grains are not entrained in any wind, we focus on geometrically simple 2D axisymmetric disc models using radiative transfer calculations with both the host star and an external source. Recent surveys have searched for evidence of massive stars influencing disc evolution using disc properties as a function of projected separation. In assuming a disc temperature of 20 K for a disc a distance D from a strong radiation source, disc masses are overestimated by a factor that scales with D−1/2 interior to the separation that external heating becomes important. This could significantly alter dust mass estimates of discs in close proximity to θ1C in the Orion Nebular Cluster (ONC). We also make an initial assessment of the effect upon snow lines. Within a parsec of an O star like θ1C a CO snow line no longer exists, though the water snow line is virtually unaffected except for very close separations of $\le 0.01\,$pc.

On the compressive nature of turbulence driven by ionizing feedback in the pillars of the Carina Nebula

ABSTRACT The ionizing radiation of massive stars sculpts the surrounding neutral gas into pillar-like structures. Direct signatures of star formation through outflows and jets are observed in these structures, typically at their tips. Recent numerical simulations have suggested that this star formation could potentially be triggered by photoionizing radiation, driving compressive modes of turbulence in the pillars. In this study, we use recent high-resolution ALMA observations of 12CO, 13CO, and C18O, J = 2 − 1 emission to test this hypothesis for pillars in the Carina Nebula. We analyse column density and intensity-weighted velocity maps, and subtract any large-scale bulk motions in the plane of the sky to isolate the turbulent motions. We then reconstruct the dominant turbulence driving mode in the pillars, by computing the turbulence driving parameter b, characterized by the relation $\sigma _{\rho /\rho _0} = b \mathcal {M}$ between the standard deviation of the density contrast $\sigma _{\rho /\rho _0}$ (with gas density ρ and its average ρ0) and the turbulent Mach number $\mathcal {M}$. We find values of b ∼ 0.7–1.0 for most of the pillars, suggesting that predominantly compressive modes of turbulence are driven in the pillars by the ionizing radiation from nearby massive stars. We find that this range of b values can produce star formation rates in the pillars that are a factor ∼3 greater than with b ∼ 0.5, a typical average value of b for spiral-arm molecular clouds. Our results provide further evidence for the potential triggering of star formation in pillars through compressive turbulent motions.

Star cluster formation and cloud dispersal by radiative feedback: dependence on metallicity and compactness

ABSTRACT We study star cluster formation in various environments with different metallicities and column densities by performing a suite of 3D radiation hydrodynamics simulations. We find that the photoionization feedback from massive stars controls the star formation efficiency (SFE) in a star-forming cloud, and its impact sensitively depends on the gas metallicity Z and initial cloud surface density Σ. At Z = 1 Z⊙, SFE increases as a power law from 0.03 at Σ = 10 M⊙ pc−2 to 0.3 at $\Sigma = 300\,\mathrm{M}_{\odot }\, {\rm pc^{-2}}$. In low-metallicity cases $10^{-2}\!-\!10^{-1}\, \mathrm{Z}_{\odot }$, star clusters form from atomic warm gases because the molecule formation time is not short enough with respect to the cooling or dynamical time. In addition, the whole cloud is disrupted more easily by expanding H ii bubbles that have higher temperature owing to less efficient cooling. With smaller dust attenuation, the ionizing radiation feedback from nearby massive stars is stronger and terminate star formation in dense clumps. These effects result in inefficient star formation in low-metallicity environments: the SFE drops by a factor of ∼3 at Z = 10−2 Z⊙ compared to the results for Z = 1 Z⊙, regardless of Σ. Newborn star clusters are also gravitationally less bound. We further develop a new semi-analytical model that can reproduce the simulation results well, particularly the observed dependencies of the SFEs on the cloud surface densities and metallicities.

Carina High-contrast Imaging Project for massive Stars (CHIPS)

Context.Massive stars like company. However, low-mass companions have remained extremely difficult to detect at angular separations (ρ) smaller than 1″ (approx. 1000–3000 au, considering the typical distance to nearby massive stars) given the large brightness contrast between the companion and the central star. Constraints on the low-mass end of the companions mass-function for massive stars are needed, however, for helping, for example, to distinguish among the various scenarios that describe the formation of massive stars.Aims.With the aim of obtaining a statistically significant constraint on the presence of low-mass companions beyond the typical detection limit of current surveys (Δmag ≲ 5 atρ ≲ 1″), we initiated a survey of O and Wolf-Rayet stars in the Carina region using the Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) coronagraphic instrument on the Very Large Telescope (VLT). In this, the first paper of the series, we aim to introduce the survey, to present the methodology and to demonstrate the capability of SPHERE for massive stars using the multiple system QZ Car.Methods.We obtained VLT-SPHERE snapshot observations in the IRDIFS_EXT mode, which combines the IFS and IRDIS sub-systems and simultaneously provides us four-dimensional (4D) data cubes in two different fields-of-view: 1.73″ × 1.73″ for IFS (39 spectral channels across theYJHbands) and 12″ × 12″ for IRDIS (two spectral channels across theKband). Angular- and spectral-differential imaging techniques as well as PSF-fitting were applied to detect and measure the relative flux of the companions in each spectral channel. The latter were then flux-calibrated using theoretical SED models of the central object and compared to a grid of ATLAS9 atmosphere model and (pre-)main-sequence evolutionary tracks, providing a first estimate of the physical properties of the detected companions.Results.Detection limits of 9 mag atρ &gt; 200 mas for IFS, and as faint as 13 mag atρ &gt; 1.​″8 for IRDIS (corresponding to sub-solar masses for potential companions), can be reached in snapshot observations of only a few minutes integration times, allowing us to detect 19 sources around the QZ Car system. All but two are reported here for the first time. With near-IR magnitude contrasts in the range of 4 to 7.5 mag, the three brightest sources (Ab, Ad, and E) are most likely to be physically bound. They have masses in the range of 2 to 12M⊙and are potentially co-eval with QZ Car central system. The remaining sources have flux contrast of 1.5 × 105to 9.5 × 106(ΔK ≈ 11 to 13 mag). Their presence can be explained by the local source density and they are, thus, likely to be chance alignments. If they were members of the Carina nebula, they would be sub-solar-mass pre-main sequence stars.Conclusions.Based on this proof of concept, we show that the VLT/SPHERE allows us to reach the sub-solar mass regime of the companion mass function. It paves the way for this type of observation with a large sample of massive stars to provide novel constraints on the multiplicity of massive stars in a region of the parameter space that has remained inaccessible so far.

Physical conditions in the gas phases of the giant H II region LMC-N 11

Context. The ambiguous origin of the [C II] 158μm line in the interstellar medium complicates its use for diagnostics concerning the star-formation rate and physical conditions in photodissociation regions. Aims. We investigate the origin of [C II] in order to measure the total molecular gas content, the fraction of CO-dark H2 gas, and how these parameters are impacted by environmental effects such as stellar feedback. Methods. We observed the giant H II region N 11 in the Large Magellanic Cloud with SOFIA/GREAT. The [C II] line is resolved in velocity and compared to H I and CO, using a Bayesian approach to decompose the line profiles. A simple model accounting for collisions in the neutral atomic and molecular gas was used in order to derive the H2 column density traced by C+. Results. The profile of [C II] most closely resembles that of CO, but the integrated [C II] line width lies between that of CO and that of H I. Using various methods, we find that [C II] mostly originates from the neutral gas. We show that [C II] mostly traces the CO-dark H2 gas but there is evidence of a weak contribution from neutral atomic gas preferentially in the faintest components (as opposed to components with low [C II]/CO or low CO column density). Most of the molecular gas is CO-dark. The CO-dark H2 gas, whose density is typically a few 100s cm−3 and thermal pressure in the range 103.5−5 K cm−3, is not always in pressure equilibrium with the neutral atomic gas. The fraction of CO-dark H2 gas decreases with increasing CO column density, with a slope that seems to depend on the impinging radiation field from nearby massive stars. Finally we extend previous measurements of the photoelectric-effect heating efficiency, which we find is constant across regions probed with Herschel, with [C II] and [O I] being the main coolants in faint and diffuse, and bright and compact regions, respectively, and with polycyclic aromatic hydrocarbon emission tracing the CO-dark H2 gas heating where [C II] and [O I] emit. Conclusions. We present an innovative spectral decomposition method that allows statistical trends to be derived for the molecular gas content using CO, [C II], and H I profiles. Our study highlights the importance of velocity-resolved photodissociation region (PDR) diagnostics and higher spatial resolution for H I observations as future steps.

Double-peaked [O i] Profile: A Likely Signature of the Gaseous Ring around KH 15D

Abstract KH 15D is a well-known spectroscopic binary because of its unique and dramatic photometric variability. The variability is explained by a circumbinary dust ring, but the ring itself was never directly detected. We present a new interpretation of the double-peaked [O I] λ6300 profiles as originating from the hot disk surface of KH 15D. By modeling these profiles, we measure emitting radii between ∼0.5 and 5 au, basically a gaseous ring very similar in radial extent to the dust ring inferred from modeling the system’s photometric variability. We also discuss the possibility that external photoevaporation driven by ultraviolet photons from the nearby massive star HD 47887 has truncated the outer edge of the disk to the observed value.

Constraining physical conditions for the PDR of Trumpler 14 in the Carina Nebula

We investigate the physical conditions of the CO gas, based on the submillimeter imaging spectroscopy from a 2′ × 7′ (1.5 × 5 pc2) area near the young star cluster, Trumpler 14 of the Carina Nebula. The observations presented in this work are taken with the Fourier Transform Spectrometer (FTS) of the Spectral and Photometric Imaging REceiver (SPIRE) onboard the Herschel Space Observatory. The newly observed spectral lines include [CI] 370 μm [CI] 609 μm, and CO transitions from J = 4−3 to J = 13−12. Our field of view covers the edge of a cavity carved by Trumpler 14 about 1 Myr ago and marks the transition from H ii regions to photo-dissociation regions. The observed CO intensities are the most prominent at the northwest region, Car I-E. With the state-of-the-art Meudon PDR code, we successfully derive the physical conditions, which include the thermal pressure (P) and the scaling factor of radiation fields (GUV), from the observed CO spectral line energy distributions (SLEDs) in the observed region. The derived GUV values generally show excellent agreement with the UV radiation fields created by nearby OB-stars and thus confirm that the main excitation source of the observed CO emission is the UV-photons provided by the massive stars. The derived thermal pressure is in the range 0.5−3 × 108 K cm-3 with the highest values found along the ionization front in Car I-E region facing Trumpler 14, hinting that the cloud structure is similar to the recent observations of the Orion Bar. We also note a discrepancy at a local position (&lt;0.17 × 0.17 pc2) between the photo-dissociation region (PDR) modeling result and the UV radiation fields estimated from nearby massive stars, which requires further investigation on nearby objects that could contribute to local heating, including outflow. Comparing the derived thermal pressure with the radiation fields, we report the first observationally derived and spatially resolved P ~ 2 × 104 GUV relationship. As direct comparisons of the modeling results to the observed 13CO, [O I] 63 μm, and [C II] 158 μm intensities are not straightforward, we urge the reader to be cautious when constraining the physical conditions of PDRs with combinations of 12CO, 13CO, [C I], [O I] 63 μm, and [C II] 158 μm observations.