Showing papers on "Infrared dark cloud published in 2019"

Widespread Molecular Outflows in the Infrared Dark Cloud G28.37+0.07: Indications of Orthogonal Outflow-filament Alignment

Abstract: We present ALMA CO(2-1) observations toward a massive infrared dark cloud G28.37+0.07. The ALMA data reveal numerous molecular (CO) outflows with a wide range of sizes throughout the cloud. Sixty-two 1.3 mm continuum cores were identified to be driving molecular outflows. We have determined the position angle in the plane-of-sky of 120 CO outflow lobes and studied their distribution. We find that the distribution of the plane-of-sky outflow position angles peaks at about 100 degrees, corresponding to a concentration of outflows with an approximately east-west direction. For most outflows, we have been able to estimate the plane-of-sky angle between the outflow axis and the filament that harbors the protostar that powers the outflow. Statistical tests strongly indicate that the distribution of outflow-filament orientations is consistent with most outflow axes being mostly orthogonal to their parent filament in three dimensions. Such alignment may result from filament fragmentation or continuous mass transportation from the filament to the embedded protostellar core. The latter is suggested by recent numerical studies with moderately strong magnetic fields.

Gravity, Magnetic Field, and Turbulence: Relative Importance and Impact on Fragmentation in the Infrared Dark Cloud G34.43+00.24

Abstract: We investigate the interplay between magnetic (B) field, gravity, and turbulence in the fragmentation process of cores within the filamentary infrared dark cloud G34.43+00.24. We observe the magnetic field morphology across G34.43, traced with thermal dust polarization at 350 μm with an angular resolution of 10'' (0.18 pc), and compare with the kinematics obtained from N2H+ across the filament. We derive local velocity gradients from N2H+, tracing motion in the plane of sky, and compare with the observed local B field orientations in the plane of sky. The B field orientations are found to be perpendicular to the long axis of the filament toward the MM1 and MM2 ridge, suggesting that the B field can guide material toward the filament. Toward MM3, the B field orientations appear more parallel to the filament and aligned with the elongated core of MM3, indicating a different role of the B field. In addition to a large-scale east–west velocity gradient, we find a close alignment between local B field orientations and local velocity gradients toward the MM1/MM2 ridge. This local correlation in alignment suggests that gas motions are influenced by the B field morphology or vice versa. Additionally, this alignment seems to become even closer with increasing integrated emission in N2H+, possibly indicating that a growing gravitational pull alignes the B field and gas motion more and more. We analyze and quantify B field, gravity, turbulence, and their relative importance toward the MM1, MM2, and MM3 regions with various techniques over two scales, a larger clump area at 2 pc scale and the smaller core area at 0.6 pc scale. While gravitational energy, B field, and turbulent pressure all grow systematically from large to small scale, the ratios among the three constituents clearly develop differently over scale. We propose that this varying relative importance between B field, gravity, and turbulence over scale drives and explains the different fragmentation types seen at subparsec scale (no fragmentation in MM1; aligned fragmentation in MM2; clustered fragmentation in MM3). We discuss uncertainties, subtleties, and the robustness of our conclusion, and we stress that a multiscale joint analysis is required to understand the dynamics in these systems.

Magnetized interstellar molecular clouds - II. The large-scale structure and dynamics of filamentary molecular clouds

Abstract: We perform ideal MHD high resolution AMR simulations with driven turbulence and self-gravity and find that long filamentary molecular clouds are formed at the converging locations of large-scale turbulence flows and the filaments are bounded by gravity. The magnetic field helps shape and reinforce the long filamentary structures. The main filamentary cloud has a length of ~4.4 pc. Instead of a monolithic cylindrical structure, the main cloud is shown to be a collection of fiber/web-like sub-structures similar to filamentary clouds such as L1495. Unless the line-of-sight is close to the mean field direction, the large-scale magnetic field and striations in the simulation are found roughly perpendicular to the long axis of the main cloud, similar to 1495. This provides strong support for a large-scale moderately strong magnetic field surrounding L1495. We find that the projection effect from observations can lead to incorrect interpretations of the true three-dimensional physical shape, size, and velocity structure of the clouds. Helical magnetic field structures found around filamentary clouds that are interpreted from Zeeman observations can be explained by a simple bending of the magnetic field that pierces through the cloud. We demonstrate that two dark clouds form a T-shape configuration which are strikingly similar to the Infrared dark cloud SDC13 leading to the interpretation that SDC13 results from a collision of two long filamentary clouds. We show that a moderately strong magnetic field (M_A ~ 1) is crucial for maintaining a long and slender filamentary cloud for a long period of time ~0.5 million years.

The Core Mass Function in the Infrared Dark Cloud G28.37+0.07

Abstract: In this paper, we analyze the 1.3 mm continuum ALMA data that cover the majority of the infrared dark cloud (IRDC) G28.37+0.07. With a spatial resolution of 05 (2500 au at 5 kpc), the continuum image reveals five groups of dense cores. Each core group has a projected physical scale of about 1 pc, with core masses spanning a dynamic range of about 100. We use the dendrogram method (astrodendro) and a newly developed graph method (astrograph) to identify individual cores. The core masses are estimated through the millimeter continuum flux, assuming constant temperature and using an NH3-based gas temperature. We construct core mass functions (CMFs) based on the two methods and fit a power-law relation dN/d log M ∝ M −α to the CMFs for M > 0.79 M ⊙. In the constant-temperature scenario, astrograph gives α = 0.80 ± 0.10, while astrodendro gives α = 0.71 ± 0.11, both significantly shallower than the Salpeter-type initial mass function with α = 1.35. In the scenario where the NH3 gas temperature is applied to cores, astrograph gives α = 1.37 ± 0.06, while astrodendro gives α = 0.87 ± 0.07. Regional CMF slope variation is seen between the core groups. We also compare CMFs in three different environments, including IRDC G28.37+0.07, IRDC clumps, and G286.21+0.17, using the identical dendrogram method. Results show that IRDCs have smaller α than the cluster-forming cloud G286.21+0.17.

On the size of the CO-depletion radius in the IRDC G351.77−0.51

Abstract: An estimate of the degree of CO-depletion ($f_D$) provides information on the physical conditions occurring in the innermost and densest regions of molecular clouds. A key parameter in these studies is the size of the depletion radius, i.e. the radius within which the C-bearing species, and in particular CO, are largely frozen onto dust grains. A strong depletion state (i.e. $f_D>10$, as assumed in our models) is highly favoured in the innermost regions of dark clouds, where the temperature is $<20$ K and the number density of molecular hydrogen exceeds a few $\times$10$^{4}$ cm$^{-3}$. In this work, we estimate the size of the depleted region by studying the Infrared Dark Cloud (IRDC) G351.77-0.51. Continuum observations performed with the $Herschel$ $Space$ $Observatory$ and the $LArge$ $APEX$ $BOlometer$ $CAmera$, together with APEX C$^{18}$O and C$^{17}$O J=2$\rightarrow$1 line observations, allowed us to recover the large-scale beam- and line-of-sight-averaged depletion map of the cloud. We built a simple model to investigate the depletion in the inner regions of the clumps in the filament and the filament itself. The model suggests that the depletion radius ranges from 0.02 to 0.15 pc, comparable with the typical filament width (i.e.$\sim$0.1 pc). At these radii, the number density of H$_2$ reaches values between 0.2 and 5.5$\times$10$^{5}$ cm$^{-3}$. These results provide information on the approximate spatial scales on which different chemical processes operate in high-mass star-forming regions and also suggest caution when using CO for kinematical studies in IRDCs.

Gravity, Magnetic Field, and Turbulence: Relative Importance and Impact on Fragmentation in the Infrared Dark Cloud G34.43+00.24

Abstract: (Abbreviated) We investigate the interplay between magnetic (B) field, gravity, and turbulence in the fragmentation process of cores within the filamentary infrared dark cloud G34.43+00.24. We observe the magnetic field (B) morphology across G34.43 and compare with the kinematics obtained from N2H+ across the filament. We derive local velocity gradients from N2H+, tracing motion in the plane of sky, and compare with the observed local B field orientations in the plane of sky. Besides a large-scale east-west velocity gradient, we find a close alignment between local B field orientations and local velocity gradients toward the MM1/MM2 ridge. This local correlation in alignment suggests that gas motions are influenced by the magnetic field morphology or vice versa. Additionally, this alignment seems to be getting even closer with increasing integrated emission in N2H+, possibly indicating that a growing gravitational pull is more and more aligning B field and gas motion. We analyze and quantify B field, gravity, turbulence, and their relative importance toward the MM1, MM2 and MM3 regions with various techniques over two scales, a larger clump area at 2 pc scale and the smaller core area at 0.6 pc scale. While gravitational energy, B field, and turbulent pressure all grow systematically from large to small scale, the ratios among the three constituents develop clearly differently over scale. We propose that this varying relative importance between B field, gravity, and turbulence over scale drives and explains the different fragmentation types seen at sub-pc scale (no fragmentation in MM1; aligned fragmentation in MM2; clustered fragmentation in MM3). We discuss uncertainties, subtleties, and the robustness of our conclusion, and we stress the need of a multi-scale joint analysis to understand the dynamics in these systems.

Interstellar Plunging Waves: ALMA resolves the physical structure of non-stationary MHD shocks.

Abstract: Magneto-hydrodynamic (MHD) shocks are violent events that inject large amounts of energy in the interstellar medium (ISM) dramatically modifying its physical properties and chemical composition. Indirect evidence for the presence of such shocks has been reported from the especial chemistry detected toward a variety of astrophysical shocked environments. However, the internal physical structure of these shocks remains unresolved since their expected spatial scales are too small to be measured with current instrumentation. Here we report the first detection of a fully spatially resolved, MHD shock toward the Infrared Dark Cloud (IRDC) G034.77-00.55. The shock, probed by Silicon Monoxide (SiO) and observed with the Atacama Large Millimetre/sub-millimetre Array (ALMA), is associated with the collision between the dense molecular gas of the cloud and a molecular gas flow pushed toward the IRDC by the nearby supernova remnant (SNR) W44. The interaction is occurring on sub-parsec spatial scales thanks to the enhanced magnetic field of the SNR, making the dissipation region of the MHD shock large enough to be resolved with ALMA. Our observations suggest that molecular flow-flow collisions can be triggered by stellar feedback, inducing shocked molecular gas densities compatible with those required for massive star formation.

Widespread Molecular Outflows in the Infrared Dark Cloud G28.37+0.07: Indications of Orthogonal Outflow-Filament Alignment

Abstract: We present ALMA CO(2-1) observations toward a massive infrared dark cloud G28.37+0.07. The ALMA data reveal numerous molecular (CO) outflows with a wide range of sizes throughout the cloud. Sixty-two 1.3 mm continuum cores were identified to be driving molecular outflows. We have determined the position angle in the plane-of-sky of 120 CO outflow lobes and studied their distribution. We find that the distribution of the plane-of-sky outflow position angles peaks at about 100 degree, corresponding to a concentration of outflows with an approximately east-west direction. For most outflows, we have been able to estimate the plane-of-sky angle between the outflow axis and the filament that harbors the protostar that powers the outflow. Statistical tests strongly indicate that the distribution of outflow-filament orientations is consistent with most outflow axes being mostly orthogonal to their parent filament in 3D. Such alignment may result from filament fragmentation or continuous mass transportation from filament to the embedded protostellar core. The latter is suggested by recent numerical studies with moderately strong magnetic fields.

Near-Infrared Imaging Polarimetry toward M17 SWex

Abstract: We conducted near-infrared (JHKs) imaging polarimetry toward the infrared dark cloud (IRDC) M17 SWex, including almost all of the IRDC filaments as well as its outskirts, with the polarimeter SIRPOL on the IRSF 1.4 m telescope. We revealed the magnetic fields of M17 SWex with our polarization-detected sources that were selected by some criteria based on their near-IR colors and the column densities toward them, which were derived from the Herschel data. The selected sources indicate not only that the ordered magnetic field is perpendicular to the cloud elongation as a whole, but also that at both ends of the elongated cloud the magnetic field appears to bent toward its central part, i.e., large-scale hourglass-shaped magnetic field perpendicular to the cloud elongation. In addition to this general trend, the elongations of the filamentary subregions within the dense parts of the cloud appear to be mostly perpendicular to their local magnetic fields, while the magnetic fields of the outskirts appear to follow the thin filaments that protrude from the dense parts. The magnetic strengths were estimated to be ~70-300 microG in the subregions, of which lengths and average number densities are ~3-9 pc and ~2-7x10^3 cm^{-3}, respectively, by the Davis-Chandrasekhar-Fermi method with the angular dispersion of our polarization data and the velocity dispersion derived from the C^{18}O (J=1-0) data obtained by the Nobeyama 45 m telescope. These field configurations and our magnetic stability analysis of the subregions imply that the magnetic field have controlled the formation/evolution of the M17 SWex cloud.

Extreme fragmentation and complex kinematics at the center of the L1287 cloud

Abstract: Aims. The filamentary ~10-pc-scale infrared dark cloud L1287 located at a parallax distance of ~929 pc is actively forming a dense cluster of low-mass young stellar objects (YSOs) at its inner ~0.1 pc region. To help understand the origin of this low-mass YSO cluster, the present work aims at resolving the gas structures and kinematics with high angular resolution.Methods. We performed ~1′′ angular resolution (~930 AU) observations at ~1.3 mm wavelengths using the Submillimeter Array (SMA), which simultaneously cover the dust continuum emission and various molecular line tracers for dense gas, warm gas, shocks, and outflows.Results. From a 1.3-mm continuum image with a resolution of ~2′′ we identified six dense cores, namely SMA1-6. Their gas masses are in the range of ~0.4– 4 M ⊙ . From a 1.3-mm continuum image with a resolution of ~1′′, we find a high fragmentation level, with 14 compact millimeter sources within 0.1 pc: SMA3 contains at least nine internal condensations; SMA5 and SMA6 are also resolved with two internal condensations. Intriguingly, one condensation in SMA3 and another in SMA5 appear associated with the known accretion outburst YSOs RNO 1C and RNO 1B. The dense gas tracer DCN (3–2) well traces the dust continuum emission and shows a clear velocity gradient along the NW-SE direction centered at SMA3. There is another velocity gradient with opposite direction around the most luminous YSO, IRAS 00338 + 6312.Conclusions. The fragmentation within 0.1 pc in L1287 is very high compared to other regions at the same spatial scales. The incoherent motions of dense gas flows are sometimes interpreted by being influenced by (proto)stellar feedback (e.g., outflows), which is not yet ruled out in this particular target source. On the other hand, the velocities (with respect to the systemic velocity) traced by DCN are small, and the directions of the velocity gradients traced by DCN are approximately perpendicular to those of the dominant CO outflow(s). Therefore, we alternatively hypothesize that the velocity gradients revealed by DCN trace the convergence from the ≳0.1 pc scales infalling motion towards the rotational motions around the more compact (~0.02 pc) sources. This global molecular gas converging flow may feed the formation of the dense low-mass YSO cluster. Finally, we also found that IRAS 00338 + 6312 is the most likely powering source of the dominant CO outflow. A compact blue-shifted outflow from RNO 1C is also identified.

The CARMA-NRO Orion Survey: Core Emergence and Kinematics in the Orion A Cloud

Abstract: We have investigated the formation and kinematics of submillimeter (submm) continuum cores in the Orion A molecular cloud. A comparison between submm continuum and near-infrared extinction shows a continuum core detection threshold of A_V ~ 5–10 mag. The threshold is similar to the star formation extinction threshold of A_V ~ 7 mag proposed by recent work, suggesting a universal star formation extinction threshold among clouds within 500 pc to the Sun. A comparison between the Orion A cloud and a massive infrared dark cloud G28.37+0.07 indicates that Orion A produces more dense gas within the extinction range 15 mag ≾ A V ≾ 60 mag. Using data from the CARMA-NRO Orion Survey, we find that dense cores in the integral-shaped filament (ISF) show subsonic core-to-envelope velocity dispersion that is significantly less than the local envelope line dispersion, similar to what has been found in nearby clouds. Dynamical analysis indicates that the cores are bound to the ISF. An oscillatory core-to-envelope motion is detected along the ISF. Its origin is to be further explored.

Interstellar Plunging Waves: ALMA Resolves the Physical Structure of Nonstationary MHD Shocks

Abstract: Magnetohydrodynamic (MHD) shocks are violent events that inject large amounts of energy in the interstellar medium dramatically modifying its physical properties and chemical composition. Indirect evidence for the presence of such shocks has been reported from the especial chemistry detected toward a variety of astrophysical shocked environments. However, the internal physical structure of these shocks remains unresolved since their expected spatial scales are too small to be measured with current instrumentation. Here we report the first detection of a fully spatially resolved, MHD shock toward the infrared dark cloud (IRDC) G034.77-00.55. The shock, probed by silicon monoxide (SiO) and observed with the Atacama Large Millimeter/submillimeter Array (ALMA), is associated with the collision between the dense molecular gas of the cloud and a molecular gas flow pushed toward the IRDC by the nearby supernova remnant (SNR) W44. The interaction is occurring on subparsec spatial scales thanks to the enhanced magnetic field of the SNR, making the dissipation region of the MHD shock large enough to be resolved with ALMA. Our observations suggest that molecular flow–flow collisions can be triggered by stellar feedback, inducing shocked molecular gas densities compatible with those required for massive star formation.

A survey of molecular cores in M17 SWex

Abstract: A survey of molecular cores covering the infrared dark cloud known as the M17 southwest extension (M17 SWex) has been carried out with the 45 m Nobeyama Radio Telescope. Based on the N2H+ (J= 1-0) data obtained, we have identified 46 individual cores whose masses are in the range 43 to 3026 Mo. We examined the relationship between the physical parameters of the cores and those of young stellar objects (YSOs) associated with the cores found in the literature. The comparison of the virial mass and the core mass indicates that most of the cores can be gravitationally stable if we assume a large external pressure. Among the 46 cores, we found four massive cores with YSOs. They have large mass of >~ 1000Mo and line width of >~ 2.5 km s^-1 which are similar to those of clumps forming high mass stars. However, previous studies have shown that there is no active massive star formation in this region. Recent measurements of near-infrared polarization infer that the magnetic field around M17 SWex is likely to be strong enough to support the cores against self-gravity. We therefore suggest that the magnetic field may prevent the cores from collapsing, causing the low-level of massive star formation in M17 SWex.

Massive Young Stellar Objects and Outflow in the Infrared Dark Cloud G79.3+0.3

Abstract: G79.3+0.3 is an infrared dark cloud in the Cygnus-X complex that is home to massive deeply embedded young stellar objects (YSOs). We have produced a submillimeter array (SMA) 1.3 mm continuum image and $^{12}$CO line maps of the eastern section of G79.3+0.3 in which we detect five separate YSOs. We have estimated physical parameters for these five YSOs and others in the vicinity of G79.3+0.3 by fitting existing photometry from Spitzer, Herschel, and ground-based telescopes to spectral energy distribution models. Through these model fits we find that the most massive YSOs seen in the SMA 1.3mm continuum emission have masses in the $5 - 6M_{sun}$ range. One of the SMA sources was observed to power a massive collimated $^{12}$CO outflow extending at least 0.94pc in both directions from the protostar, with a total mass of 0.83$M_{sun}$ and a dynamical timescale of 23kyr.

Millimeter and far-IR study of the IRCD SDC 341.232-0.268

Abstract: We analyze the molecular gas and dust associated with the infrared dark cloud SDC 341232-0268 in order to investigate the characteristics and parameters of the gas, determine the evolutionary status of four embedded EGO candidates and establish possible infall or outflow gas motions We based our study on $^{12}$CO(2-1), $^{13}$CO(2-1), and C$^{18}$O(2-1) data obtained with the APEX telescope, molecular data of high density tracers from the MALT90 survey and IR images from Spitzer, Herschel and ATLASGAL The study reveals two clumps at -44 km/s towards the IRDC, with densities of > 10$^4$ cm$^{-3}$ , typical of IRDCs, while high density tracers show H$_2$ densities > 10$^5$ FIR images reveals the presence of cold dust linked to the molecular clumps and EGOs A comparison of spectra of the optically thin and optically thick molecular lines towards the EGOs suggests the existence of infall and outflow motions

Massive Young Stellar Objects and Outflow in the Infrared-Dark Cloud G79.3+0.3

Abstract: G79.3+0.3 is an infrared dark cloud in the Cygnus-X complex that is home to massive deeply embedded young stellar objects (YSOs). We have produced a submillimeter array (SMA) 1.3 mm continuum image and $^{12}$CO line maps of the eastern section of G79.3+0.3 in which we detect five separate YSOs. We have estimated physical parameters for these five YSOs and others in the vicinity of G79.3+0.3 by fitting existing photometry from Spitzer, Herschel, and ground-based telescopes to spectral energy distribution models. Through these model fits we find that the most massive YSOs seen in the SMA 1.3mm continuum emission have masses in the $5 - 6M_{sun}$ range. One of the SMA sources was observed to power a massive collimated $^{12}$CO outflow extending at least 0.94pc in both directions from the protostar, with a total mass of 0.83$M_{sun}$ and a dynamical timescale of 23kyr.

The CARMA-NRO Orion Survey: Core Emergence and Kinematics in the Orion A Cloud.

Abstract: We have investigated the formation and kinematics of sub-mm continuum cores in the Orion A molecular cloud. A comparison between sub-mm continuum and near infrared extinction shows a continuum core detection threshold of $A_V\sim$ 5-10 mag. The threshold is similar to the star formation extinction threshold of $A_V\sim$ 7 mag proposed by recent work, suggesting a universal star formation extinction threshold among clouds within 500 pc to the Sun. A comparison between the Orion A cloud and a massive infrared dark cloud G28.37+0.07 indicates that Orion A produces more dense gas within the extinction range 15 mag $\lesssim A_V \lesssim$ 60 mag. Using data from the CARMA-NRO Orion Survey, we find that dense cores in the integral-shaped filament (ISF) show sub-sonic core-to-envelope velocity dispersion that is significantly less than the local envelope line dispersion, similar to what has been found in nearby clouds. Dynamical analysis indicates that the cores are bound to the ISF. An oscillatory core-to-envelope motion is detected along the ISF. Its origin is to be further explored.

Study of the 13 CO/C 18 O abundance ratio towards the filamentary infrared dark cloud IRDC 34.43 + 0.24

Abstract: Nowadays, there are several observational studies about the 13CO/C18O abundance ratio ( ) towards nearby molecular clouds. These works give observational support to the C18O selective photodissociation due to the interaction between the far ultraviolet (FUV) radiation and the molecular gas. It is necessary to increase the sample of molecular clouds located at different distances and affected in different ways by nearby or embedded H ii regions and OB associations to study the selective photodissociation. Using 12CO, 13CO, and C18O J = 1–0 data obtained from the FOREST unbiased Galactic plane imaging survey performed with the Nobeyama 45-m telescope, we analyse the filamentary infrared dark cloud IRDC located at the distance of about 3.9 kpc. This infrared dark cloud (IRDC) is related to several H ii regions and young stellar objects. Assuming local thermodynamic equilibrium, we obtain: N(13CO) cm–2 (average value cm–2), N(C18O) cm–2 (average value cm–2), and 3 30 (average ) across the whole IRDC. Larger values of were found towards portions of the cloud related to the H ii regions associated with the N61 and N62 bubbles and with the photodissociation regions, precisely the regions in which FUV photons are strongly interacting with the molecular gas. Our result represents an observational support to the C18O selectively photodissociation phenomenon occurring in a quite distant filamentary IRDC. Additionally, based on IR data from the Hi-GAL survey, the FUV radiation field was estimated in Habing units, and the dust temperature (T ) and H2 column density (N(H2)) distribution were studied. Using the average of N(H2), values in close agreement with the ‘canonical’ abundance ratios [H2]/[13CO] and [H2]/[C18O] were derived. However, the obtained ranges in the abundance ratios show that if an accurate analysis of the molecular gas is required, the use of the ‘canonical’ values may introduce some bias. Thus, it is important to consider how the gas is irradiated by the FUV photons across the molecular cloud. The analysis of is a good tool to perform that. Effects of beam dilution and clumpiness were studied.

NIKA2 observations around LBV stars: emission from stars and circumstellar material

Abstract: Luminous Blue Variable (LBV) stars are evolved massive objects, previous to core-collapse supernova. LBVs are characterized by photometric and spectroscopic variability, produced by strong and dense winds, mass-loss events and very intense UV radiation. LBVs strongly disturb their surroundings by heating and shocking, and produce important amounts of dust. The study of the circumstellar material is therefore crucial to understand how these massive stars evolve, and also to characterize their effects onto the interstellar medium. The versatility of NIKA2 is a key in providing simultaneous observations of both the stellar continuum and the extended, circumstellar contribution. The NIKA2 frequencies (150 and 260 GHz) are in the range where thermal dust and free-free emission compete, and hence NIKA2 has the capacity to provide key information about the spatial distribution of circumstellar ionized gas, warm dust and nearby dark clouds; non-thermal emission is also possible even at these high frequencies. We show the results of the first NIKA2 survey towards five LBVs. We detected emission from four stars, three of them immersed in tenuous circumstellar material. The spectral indices show a complex distribution and allowed us to separate and characterize different components. We also found nearby dark clouds, with spectral indices typical of thermal emission from dust. Spectral indices of the detected stars are negative and hard to be explained only by free-free processes. In one of the sources, G79.29+0.46, we also found a strong correlation of the 1mm and 2mm continuum emission with respect to nested molecular shells at 0.1 pc from the LBV. The spectral index in this region clearly separates four components: the LBV star, a bubble characterized by free-free emission, and a shell interacting with a nearby infrared dark cloud.