Showing 1–50 of 100 results for author: Lada, C J

The Molecular Clouds of M31

Authors: Charles J. Lada, Jan Forbrich, Glen Petitpas, Sebastien Viaene

Abstract: Deep interferometric observations of CO and dust continuum emission are obtained with the Sub-Millimeter Array (SMA) at 230 GHz to investigate the physical nature of the giant molecular cloud (GMC) population in the Andromeda galaxy (M31). We use J = 2-1 $^{12}$CO and $^{13}$CO emission to derive the masses, sizes and velocity dispersions of 162 spatially resolved GMCs. We perform a detailed study… ▽ More Deep interferometric observations of CO and dust continuum emission are obtained with the Sub-Millimeter Array (SMA) at 230 GHz to investigate the physical nature of the giant molecular cloud (GMC) population in the Andromeda galaxy (M31). We use J = 2-1 $^{12}$CO and $^{13}$CO emission to derive the masses, sizes and velocity dispersions of 162 spatially resolved GMCs. We perform a detailed study of a subset of 117 GMCs that exhibit simple, single component line profile shapes. Examining the Larson scaling relations for these GMCs we find: 1- a highly correlated mass-size relation in both $^{12}$CO and $^{13}$CO emission; 2- a weakly correlated $^{12}$CO linewidth-size (LWS) relation along with a weaker, almost non-existent, $^{13}$CO LWS relation, suggesting a possible dependence of the LWS relation on spatial scale; and 3-that only 43\% of these GMCs are gravitationally bound. We identify two classes of GMCs based on the strength and extent of their $^{13}$CO emission. Examination of the Larson relations finds that both classes are individually characterized by strong $^{12}$CO mass-size relations and much weaker $^{12}$CO and $^{13}$CO LWS relations. The majority (73\%) of strong $^{13}$CO emitting GMCs are found to be gravitationally bound. However, only 25\% of the weak $^{13}$CO emitting GMCs are bound. The resulting breakdown in the Larson relations in the weak $^{13}$CO emitting population decouples the mass-size and LWS relations demonstrating that independent physical causes are required to understand the origin of each. Finally, in nearly every aspect, the physical properties of the M31 GMCs are found to be very similar to those of local Milky Way clouds. △ Less

Submitted 4 June, 2024; v1 submitted 28 March, 2024; originally announced March 2024.

Comments: Accepted for publication in the Astrophysical Journal. Corrected caption to Figure 13. Updated references

Journal ref: The Astrophysical Journal, 2024, Vol 966, 193

Tracing Dense Gas in Six Resolved GMCs of the Andromeda Galaxy

Authors: Jan Forbrich, Charles J. Lada, Jérôme Pety, Glen Petitpas

Abstract: We present dense-gas--tracing molecular observations of six resolved Giant Molecular Clouds (GMCs) in the Andromeda Galaxy (M31). Using the NOEMA interferometer, we observed the transitions of HCN(1-0), HCO$^+$(1-0), and HNC(1-0), as well as $^{13}$CO(1-0) and 100 GHz continuum emission. This complements our earlier work with the Submillimeter Array (SMA), including resolved dust continuum detecti… ▽ More We present dense-gas--tracing molecular observations of six resolved Giant Molecular Clouds (GMCs) in the Andromeda Galaxy (M31). Using the NOEMA interferometer, we observed the transitions of HCN(1-0), HCO$^+$(1-0), and HNC(1-0), as well as $^{13}$CO(1-0) and 100 GHz continuum emission. This complements our earlier work with the Submillimeter Array (SMA), including resolved dust continuum detections of these clouds at 230 GHz. In this work, we first compare different continuum measurements to conclude that the average free-free contamination of the observed flux is 71% at 3 mm but only 13% at 1 mm, confirming that emission at 3 mm is less reliable than that at 1 mm for calculating dust masses of star-forming clouds. While the $^{13}$CO emission is more extended than both HCN and HCO$^+$ emission, which in turn is more extended than HNC emission, we find that both HCN and HCO$^+$ are spatially coincident with, and similarly extended as, the 230 GHz dust emission. This suggests that both the 230 GHz dust continuum and most importantly the HCN emission traces the dense gas component of these GMCs. From comparison of the molecular emission with dust masses derived from the 230 GHz continuum emission, we obtain the first direct measurements of the dust-mass-to-light ratios ($α^\prime_{HCN}$ and $α^\prime_{HCO^+}$) in GMCs of an external galaxy. For HCN, the result is broadly similar to a measurement in the local Perseus cloud suggesting that these are indeed dense gas conversion factors. A larger cloud sample will be required to assess whether HCN is tracing comparable cloud-scale density regimes across the environments of M31. △ Less

Submitted 24 August, 2023; originally announced August 2023.

Comments: MNRAS, in press

A Complete HCN Survey of the Perseus Molecular Cloud

Authors: T. M. Dame, Charles J. Lada

Abstract: We present a survey of the Perseus molecular cloud in the J $=$ 1$\rightarrow$0 transition of HCN, a widely used tracer of dense molecular gas. The survey was conducted with the CfA 1.2 m telescope, which at 89 GHz has a beam width of 11' and a spectral resolution of 0.85 km s$^{-1}$. A total of 8.1 deg$^2$ was surveyed on a uniform 10' grid to a sensitivity of 14 mK per channel. The survey was co… ▽ More We present a survey of the Perseus molecular cloud in the J $=$ 1$\rightarrow$0 transition of HCN, a widely used tracer of dense molecular gas. The survey was conducted with the CfA 1.2 m telescope, which at 89 GHz has a beam width of 11' and a spectral resolution of 0.85 km s$^{-1}$. A total of 8.1 deg$^2$ was surveyed on a uniform 10' grid to a sensitivity of 14 mK per channel. The survey was compared with similar surveys of CO and dust in order to study and calibrate the HCN line as a dense gas tracer. We find the HCN emission to extend over a considerable fraction of the cloud. We show that the HCN intensity remains linear with H$_2$ column density well into the regime where the CO line saturates. We use radiative transfer modeling to show that this likely results from subthermal excitation of HCN in a cloud where the column and volume densities of H$_2$ are positively correlated. To match our HCN observations the model requires an exponential decrease in HCN abundance with increasing extinction, consistent with HCN depletion onto grains. The modeling also reveals that the mean volume density of H$_2$ in the HCN emitting regions is $\sim$ 10$^4$ cm$^{-3}$, well below the HCN critical density. For the first time, we obtain a direct measurement of the ratio of dense gas mass to HCN luminosity for an entire nearby molecular cloud: $α$(HCN) $=$ 92 M$_\odot$/(K km s$^{-1}$ pc$^2$). △ Less

Submitted 18 January, 2023; originally announced January 2023.

Comments: 14 pages, 7 figures, Accepted for the Astrophysical Journal, Jan 16, 2023

Systematic Investigation of Dust and Gaseous CO in 12 Nearby Molecular Clouds

Authors: John Arban Lewis, Charles J. Lada, Thomas Dame

Abstract: We report the first uniform and systematic study of dust and molecular gas in nearby molecular clouds. We use surveys of dust extinction and emission to determine the opacity and map the distribution of the dust within a dozen local clouds in order to derive a uniform set of basic cloud properties. We find: 1) the average dust opacity $\langleκ_{d,353}\rangle = 0.8\ {\rm cm^{2}\, g^{-1}}$ with var… ▽ More We report the first uniform and systematic study of dust and molecular gas in nearby molecular clouds. We use surveys of dust extinction and emission to determine the opacity and map the distribution of the dust within a dozen local clouds in order to derive a uniform set of basic cloud properties. We find: 1) the average dust opacity $\langleκ_{d,353}\rangle = 0.8\ {\rm cm^{2}\, g^{-1}}$ with variations of a factor of $\sim$ 2 between clouds, 2) cloud PDFs are exquisitely described by steeply falling power-laws with a narrow range of slope, and 3) a tight $M_{\rm GMC} \sim R_{\rm GMC}^2$ scaling relation for the cloud sample, indicative of a cloud population with an exactingly constant average surface density above a common fixed boundary. We compare these results to uniformly analyzed CO surveys. We measure the CO mass conversion factors and assess the efficacy of CO for tracing the physical properties of molecular clouds. We find $\langle α_{\rm CO}\rangle = 4.31 \pm 0.67$ M$_\odot$ (K km s$^{-1}$ pc$^2$)$^{-1}$ (corresponding to $X_{\rm CO}$ = 1.97 $\times$ 10$^{20}$ cm$^{-2}$(K km s$^{-1}$)$^{-1}$). We demonstrate that CO observations are a poor tracer of column density and structure on sub-cloud spatial scales. On cloud scales, CO observations can provide measurements consistent with those of the dust, provided data are analyzed in a similar, self-consistent fashion. Measurements of average GMC surface density are sensitive to choice of cloud boundary. Care must be exercised to adopt common fixed boundaries when comparing surface densities for cloud populations within and between galaxies. △ Less

Submitted 30 March, 2022; v1 submitted 14 March, 2022; originally announced March 2022.

Comments: 27 pages, 19 figures, 5 tables. Accepted ApJL. v2 Corrected bad coordinates on Appendix C maps. Brought text into alignment with final ApJL version

Simultaneous Deep Measurements of CO isotopologues and Dust Emission in Giant Molecular Clouds in the Andromeda Galaxy

Authors: Sébastien Viaene, Jan Forbrich, Charles J. Lada, Glen Petitpas, Christopher Faesi

Abstract: We present simultaneous measurements of emission from dust continuum at 230 GHz and the J=2-1 $^{12}$CO, $^{13}$CO and C$^{18}$O isotopologues at $\sim$ 15 pc resolution from individual Giant Molecular Clouds (GMCs) in the Andromeda galaxy (M31). These observations were obtained in an ongoing survey of this galaxy being conducted with the Submillimeter Array (SMA). Initial results describing the c… ▽ More We present simultaneous measurements of emission from dust continuum at 230 GHz and the J=2-1 $^{12}$CO, $^{13}$CO and C$^{18}$O isotopologues at $\sim$ 15 pc resolution from individual Giant Molecular Clouds (GMCs) in the Andromeda galaxy (M31). These observations were obtained in an ongoing survey of this galaxy being conducted with the Submillimeter Array (SMA). Initial results describing the continuum and $^{12}$CO emission were published earlier. Here we primarily analyze the observations of $^{13}$CO and C$^{18}$O emission and compare them to the measurements of dust continuum and $^{12}$CO emission. We also report additional dust continuum and CO measurements from newly added GMCs to the M31 sample. We detect spatially resolved $^{13}$CO emission with high signal-to-noise in 31 objects. We find the extent of the $^{13}$CO emission to be nearly comparable to that of $^{12}$CO, typically covering 75\% of the area of the $^{12}$CO emission. We derive $^{13}$CO and C$^{18}$O abundances of 2.9 $\times 10^{-6}$ and 4.4 $\times 10^{-7}$ relative to H$_2$, respectively, by comparison with hydrogen column densities of the same regions derived from the dust continuum observations assuming a Milky Way gas-to-dust ratio. We find the isotopic abundance ratio [$^{13}$CO]/[C$^{18}$O] = 6.7$\pm$2.9 to be consistent with the Milky Way value (8.1). Finally, we derive the mass-to-light conversion factors for all three CO species to be $α_{12} = 8.7 \pm 3.9$, $α_{13} = 48.9 \pm 20.4$ and $α_{18} = 345^{+25}_{-31}$ M$_\odot$ (K km s$^{-1}$pc$^2$)$^{-1}$ for the J=2-1 transitions of $^{12}$CO, $^{13}$CO and C$^{18}$O, respectively. △ Less

Submitted 14 February, 2021; originally announced February 2021.

Comments: 24 pages, 11 figures, accepted for publication in The Astrophysical Journal

The Mass-Size Relation and the Constancy of GMC Surface Densities in the Milky Way

Authors: C. J. Lada, T. M. Dame

Abstract: We use two existing molecular cloud catalogs derived from the same CO survey and two catalogs derived from local dust extinction surveys to investigate the nature of the GMC mass-size relation in the Galaxy. We find that the four surveys are well described by $M_{GMC} \sim R^2$ implying a constant mean surface density, $Σ_{GMC}$, for the cataloged clouds. However, the scaling coefficients and scat… ▽ More We use two existing molecular cloud catalogs derived from the same CO survey and two catalogs derived from local dust extinction surveys to investigate the nature of the GMC mass-size relation in the Galaxy. We find that the four surveys are well described by $M_{GMC} \sim R^2$ implying a constant mean surface density, $Σ_{GMC}$, for the cataloged clouds. However, the scaling coefficients and scatter differ significantly between the CO and extinction derived relations. We find that the additional scatter seen in the CO relations is due to a systematic variation in $Σ_{GMC}$ with Galactic radius that is unobservable in the local extinction data. We decompose this radial variation of $Σ_{GMC}$ into two components, a linear negative gradient with Galactic radius and a broad peak coincident with the molecular ring and superposed on the linear gradient. We show that the former may be due to a radial dependence of X$_{CO}$ on metallicity while the latter likely results from a combination of increased surface densities of individual GMCs and a systematic upward bias in the measurements of $Σ_{GMC}$ due to cloud blending in the molecular ring. We attribute the difference in scaling coefficients between the CO and extinction data to an underestimate of X$_{CO}$. We recalibrate the CO observations of nearby GMCs using extinction measurements to find that locally X$_{CO}$ $=$ 3.6$\pm$0.3 $\times$ 10$^{20}$ cm$^{-2}$ (K-km/s)$^{-1}$. We conclude that outside the molecular ring the GMC population of the Galaxy can be described to relatively good precision by a constant $Σ_{GMC}$ of 35 M$_\odot$ pc$^{-2}$. △ Less

Submitted 16 June, 2020; v1 submitted 15 June, 2020; originally announced June 2020.

Comments: 27 pages, 7 figures, accepted for publication in the ApJ, corrected typos in earlier version

First Resolved Dust Continuum Measurements of Individual Giant Molecular Clouds in the Andromeda Galaxy

Authors: Jan Forbrich, Charles J. Lada, Sébastien Viaene, Glen Petitpas

Abstract: In our local Galactic neighborhood, molecular clouds are best studied using a combination of dust measurements, to determine robust masses, sizes and internal structures of the clouds, and molecular-line observations to determine cloud kinematics and chemistry. We present here the first results of a program designed to extend such studies to nearby galaxies beyond the Magellanic Clouds. Utilizing… ▽ More In our local Galactic neighborhood, molecular clouds are best studied using a combination of dust measurements, to determine robust masses, sizes and internal structures of the clouds, and molecular-line observations to determine cloud kinematics and chemistry. We present here the first results of a program designed to extend such studies to nearby galaxies beyond the Magellanic Clouds. Utilizing the wideband upgrade of the Submillimeter Array (SMA) at 230 GHz we have obtained the first continuum detections of the thermal dust emission on sub-GMC scales ($\sim$ 15 pc) within the Andromeda galaxy (M31). These include the first resolved continuum detections of dust emission from individual GMCs beyond the Magellanic Clouds. Utilizing a powerful capability of the SMA, we simultaneously recorded CO(2-1) emission with identical $(u,\,v)$ coverage, astrometry and calibration, enabling the first measurements of the CO conversion factor, $α_{\rm\,CO(2-1)}$, toward individual GMCs across an external galaxy. Our direct measurement yields an average CO--to--dust mass conversion factor of $α^\prime_{\rm CO-dust} = 0.042\pm0.018$ $M_\odot$ (K km s$^{-1}$ pc$^2$)$^{-1}$ for the $J= 2-1$ transition. This value does not appear to vary with galactocentric radius. Assuming a constant gas-to-dust ratio of 136, the resulting $α_{\rm CO}$ $=$ 5.7 $\pm$ 2.4 $M_\odot$ (K km s$^{-1}$ pc$^2$)$^{-1}$ for the 2-1 transition is in excellent agreement with that of Milky Way GMCs, given the uncertainties. Finally, using the same analysis techniques, we compare our results with observations of the local Orion molecular clouds, placed at the distance of M31 and simulated to appear as they would if observed by the SMA. △ Less

Submitted 9 January, 2020; originally announced January 2020.

Comments: ApJ, in press

VISION - Vienna survey in Orion. III. Young stellar objects in Orion A

Authors: Josefa E. Großschedl, João Alves, Paula S. Teixeira, Hervé Bouy, Jan Forbrich, Charles J. Lada, Stefan Meingast, Álvaro Hacar, Joana Ascenso, Christine Ackerl, Birgit Hasenberger, Rainer Köhler, Karolina Kubiak, Irati Larreina, Lorenz Linhardt, Marco Lombardi, Torsten Möller

Abstract: We have extended and refined the existing young stellar object (YSO) catalogs for the Orion A molecular cloud, the closest massive star-forming region to Earth. This updated catalog is driven by the large spatial coverage (18.3 deg$^2$, $\sim$950 pc$^2$), seeing limited resolution ($\sim$0.7$"$), and sensitivity ($K_s<19$ mag) of the ESO-VISTA near-infrared survey of the Orion A cloud (VISION). Co… ▽ More We have extended and refined the existing young stellar object (YSO) catalogs for the Orion A molecular cloud, the closest massive star-forming region to Earth. This updated catalog is driven by the large spatial coverage (18.3 deg$^2$, $\sim$950 pc$^2$), seeing limited resolution ($\sim$0.7$"$), and sensitivity ($K_s<19$ mag) of the ESO-VISTA near-infrared survey of the Orion A cloud (VISION). Combined with archival mid- to far-infrared data, the VISTA data allow for a refined and more robust source selection. We estimate that among previously known protostars and pre-main-sequence stars with disks, source contamination levels (false positives) are at least $\sim$6.4% and $\sim$2.3%, respectively, mostly due to background galaxies and nebulosities. We identify 274 new YSO candidates using VISTA/Spitzer based selections within previously analyzed regions, and VISTA/WISE based selections to add sources in the surroundings, beyond previously analyzed regions. The WISE selection method recovers about 59% of the known YSOs in Orion A's low-mass star-forming part L1641, which shows what can be achieved by the all-sky WISE survey in combination with deep near-infrared data in regions without the influence of massive stars. The new catalog contains 2980 YSOs, which were classified based on the de-reddened mid-infrared spectral index into 188 protostars, 185 flat-spectrum sources, and 2607 pre-main-sequence stars with circumstellar disks. We find a statistically significant difference in the spatial distribution of the three evolutionary classes with respect to regions of high dust column-density, confirming that flat-spectrum sources are at a younger evolutionary phase compared to Class IIs, and are not a sub-sample seen at particular viewing angles. △ Less

Submitted 19 December, 2018; v1 submitted 1 October, 2018; originally announced October 2018.

Comments: 39 pages, 25 figures, Accepted for publication by A&A

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

A global correlation linking young stars, clouds, and galaxies. Towards a unified view of star formation

Authors: I. Mendigutía, C. J. Lada, R. D. Oudmaijer

Abstract: (abridged) The star formation rate (SFR) linearly correlates with the amount of dense gas mass (Mdg) involved in the formation of stars both for distant galaxies and clouds in our Galaxy. Similarly, the mass accretion rate (Macc) and the disk mass (Mdisk) of young, Class II stars are also linearly correlated. We plotted the corresponding observational data together, finding a statistically signifi… ▽ More (abridged) The star formation rate (SFR) linearly correlates with the amount of dense gas mass (Mdg) involved in the formation of stars both for distant galaxies and clouds in our Galaxy. Similarly, the mass accretion rate (Macc) and the disk mass (Mdisk) of young, Class II stars are also linearly correlated. We plotted the corresponding observational data together, finding a statistically significant correlation that spans ~ 16 orders of magnitude. This probably represents one of the widest ranges of any empirical correlation known, encompassing galaxies that are several kpc in size, pc-size star-forming clouds within our Galaxy, down to young, pre-main sequence stars with au-size protoplanetary disks. We propose a bottom-up hypothesis suggesting that a relation between Macc and the total circumstellar mass surrounding Class 0/I sources (Mcs; disk+envelope) drives the correlation in protostellar-hosting clouds and cloud-hosting galaxies. This is consistent with the fact that the SFRs derived for clouds over a timescale of 2 Myr can be roughly recovered from the sum of instantaneous Macc values of the protostars embedded within them, implying that galactic SFRs averaged over ~ 10-100 Myr should be constant over this period too. Moreover, the sum of the Mcs values directly participating in the formation of the protostellar population in a cloud likely represents a non-negligible fraction of the Mdg within the cloud. If such fraction is ~ 1-35% of the Mdg associated with star-forming clouds and galaxies, then the global correlation for all scales has a near unity slope and an intercept consistent with the (proto-)stellar accretion timescale, Mcs/Macc. Therefore, an additional critical test of our hypothesis is that the Macc-Mdisk correlation for Class II stars should also be observed between Macc and Mcs for Class 0/I sources with similar slope and intercept. △ Less

Submitted 22 August, 2018; originally announced August 2018.

Comments: Accepted in A&A. 5 figures, 4 Appendices

Journal ref: A&A 618, A119 (2018)

3D shape of Orion A from Gaia DR2

Authors: Josefa E. Grossschedl, Joao Alves, Stefan Meingast, Christine Ackerl, Joana Ascenso, Herve Bouy, Andreas Burkert, Jan Forbrich, Verena Fuernkranz, Alyssa Goodman, Alvaro Hacar, Gabor Herbst-Kiss, Charles J. Lada, Irati Larreina, Kieran Leschinski, Marco Lombardi, Andre Moitinho, Daniel Mortimer, Eleonora Zari

Abstract: We use the $\mathit{Gaia}$ DR2 distances of about 700 mid-infrared selected young stellar objects in the benchmark giant molecular cloud Orion A to infer its 3D shape and orientation. We find that Orion A is not the fairly straight filamentary cloud that we see in (2D) projection, but instead a cometary-like cloud oriented toward the Galactic plane, with two distinct components: a denser and enhan… ▽ More We use the $\mathit{Gaia}$ DR2 distances of about 700 mid-infrared selected young stellar objects in the benchmark giant molecular cloud Orion A to infer its 3D shape and orientation. We find that Orion A is not the fairly straight filamentary cloud that we see in (2D) projection, but instead a cometary-like cloud oriented toward the Galactic plane, with two distinct components: a denser and enhanced star-forming (bent) Head, and a lower density and star-formation quieter $\sim$75 pc long Tail. The true extent of Orion A is not the projected $\sim$40 pc but $\sim$90 pc, making it by far the largest molecular cloud in the local neighborhood. Its aspect ratio ($\sim$30:1) and high column-density fraction ($\sim45\%$) make it similar to large-scale Milky Way filaments ("bones"), despite its distance to the galactic mid-plane being an order of magnitude larger than typically found for these structures. △ Less

Submitted 17 August, 2018; originally announced August 2018.

Comments: 9 pages, 5 figures, 4 tables, accepted for publication in A&A

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

The HP2 Survey - IV. The Pipe nebula: Effective dust temperatures in dense cores

Authors: Birgit Hasenberger, Marco Lombardi, João Alves, Jan Forbrich, Alvaro Hacar, Charles J. Lada

Abstract: Multi-wavelength observations in the sub-mm regime provide information on the distribution of both the dust column density and the effective dust temperature in molecular clouds. In this study, we created high-resolution and high-dynamic-range maps of the Pipe nebula region and explored the value of dust-temperature measurements in particular towards the dense cores embedded in the cloud. The maps… ▽ More Multi-wavelength observations in the sub-mm regime provide information on the distribution of both the dust column density and the effective dust temperature in molecular clouds. In this study, we created high-resolution and high-dynamic-range maps of the Pipe nebula region and explored the value of dust-temperature measurements in particular towards the dense cores embedded in the cloud. The maps are based on data from the Herschel and Planck satellites, and calibrated with a near-infrared extinction map based on 2MASS observations. We have considered a sample of previously defined cores and found that the majority of core regions contain at least one local temperature minimum. Moreover, we observed an anti-correlation between column density and temperature. The slope of this anti-correlation is dependent on the region boundaries and can be used as a metric to distinguish dense from diffuse areas in the cloud if systematic effects are addressed appropriately. Employing dust-temperature data thus allows us to draw conclusions on the thermodynamically dominant processes in this sample of cores: external heating by the interstellar radiation field and shielding by the surrounding medium. In addition, we have taken a first step towards a physically motivated core definition by recognising that the column-density-temperature anti-correlation is sensitive to the core boundaries. Dust-temperature maps therefore clearly contain valuable information about the physical state of the observed medium. △ Less

Submitted 11 July, 2018; originally announced July 2018.

Comments: 14 pages, 22 figures. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 620, A24 (2018)

The ALMA view of GMCs in NGC 300: Physical Properties and Scaling Relations at 10 pc Resolution

Authors: Christopher M. Faesi, Charles J. Lada, Jan Forbrich

Abstract: We have conducted a $^{12}$CO(2-1) survey of several molecular gas complexes in the vicinity of H II regions within the spiral galaxy NGC 300 using the Atacama Large Millimeter Array. Our observations attain a resolution of 10 pc and 1 km s$^{-1}$, sufficient to fully resolve Giant Molecular Clouds (GMCs), and are the highest to date obtained beyond the Local Group. We use the CPROPS algorithm to… ▽ More We have conducted a $^{12}$CO(2-1) survey of several molecular gas complexes in the vicinity of H II regions within the spiral galaxy NGC 300 using the Atacama Large Millimeter Array. Our observations attain a resolution of 10 pc and 1 km s$^{-1}$, sufficient to fully resolve Giant Molecular Clouds (GMCs), and are the highest to date obtained beyond the Local Group. We use the CPROPS algorithm to identify and characterize 250 GMCs across the observed regions. GMCs in NGC 300 appear qualitatively and quantitatively similar to those in the Milky Way disk: they show an identical scaling relationship between size $R$ and linewidth $ΔV$ ($ΔV \propto R^{0.48\pm0.05}$), appear to be mostly in virial equilibrium, and are consistent with having a constant surface density of $60$ M$_{\odot}$ pc$^{-2}$. The GMC mass spectrum is similar to those in the inner disks of spiral galaxies (including the Milky Way). Our results suggest that global galactic properties such as total stellar mass, morphology, and average metallicity may not play a major role in setting GMC properties, at least within the disks of galaxies on the star-forming main sequence. Instead, GMC properties may be more strongly influenced by local environmental factors such as the mid-plane disk pressure. In particular, in the inner disk of NGC 300 we find this pressure to be similar to that in the local Milky Way but markedly lower than that in the disk of M51 where GMCs are characterized by systematically higher surface densities and a higher coefficient for the size-linewidth relation. △ Less

Submitted 18 January, 2018; originally announced January 2018.

Comments: 40 pages, 26 figures, accepted for publication in ApJ

Gathering dust: A galaxy-wide study of dust emission from cloud complexes in NGC 300

Authors: M. Riener, C. M. Faesi, J. Forbrich, C. J. Lada

Abstract: We used multi-band observations by the Herschel Space Observatory to study the dust emission properties of the nearby spiral galaxy NGC 300. We compiled a first catalogue of the population of giant dust clouds (GDCs) in NGC 300 and give an estimate of the total dust mass of the galaxy. We carried out source detection with the multiwavelength source extraction algorithm getsources and calculated ph… ▽ More We used multi-band observations by the Herschel Space Observatory to study the dust emission properties of the nearby spiral galaxy NGC 300. We compiled a first catalogue of the population of giant dust clouds (GDCs) in NGC 300 and give an estimate of the total dust mass of the galaxy. We carried out source detection with the multiwavelength source extraction algorithm getsources and calculated physical properties of the GDCs, including mass and temperature, from five-band Herschel PACS and SPIRE observations from 100-500 $μ$m; the final size and mass estimates are based on the observations at 250 $μ$m that have an effective spatial resolution of $\sim$170 pc. We correlated our final catalogue of GDCs to pre-existing catalogues of HII regions to infer the number of GDCs associated with high-mass star formation and determined the H$α$ emission of the GDCs. Our final catalogue of GDCs includes 146 sources, 90 of which are associated with known HII regions. We find that the dust masses of the GDCs are completely dominated by the cold dust component and range from $\sim$1.1$\cdot$10$^{3}$ - 1.4$\cdot$10$^{4}$ M$_{\odot}$. The GDCs have effective temperatures of $\sim$13-23 K and show a distinct cold dust effective temperature gradient from the centre towards the outer parts of the stellar disk. We find that the population of GDCs in our catalogue constitutes $\sim$16% of the total dust mass of NGC 300, which we estimate to be about 5.4$\cdot$10$^{6}$ M$_{\odot}$. At least about 87% of our GDCs have a high enough average dust mass surface density to provide sufficient shielding to harbour molecular clouds. We compare our results to previous pointed molecular gas observations in NGC 300 and results from other nearby galaxies and also conclude that it is very likely that most of our GDCs are associated with complexes of giant molecular clouds. △ Less

Submitted 28 November, 2017; originally announced November 2017.

Comments: 18 pages, 19 figures, appendix with 2 figures

Journal ref: A&A 612, A81 (2018)

HP2 survey: III The California Molecular Cloud--A Sleeping Giant Revisited

Authors: Charles J. Lada, John A. Lewis, Marco Lombardi, João Alves

Abstract: We present new high resolution and dynamic range dust column density and temperature maps of the California Molecular Cloud derived from a combination of Planck and Herschel dust-emission maps, and 2MASS NIR dust-extinction maps. We used these data to determine the ratio of the 2.2 micron extinction coefficient to the 850 micron opacity and found the value to be close to that found in similar stud… ▽ More We present new high resolution and dynamic range dust column density and temperature maps of the California Molecular Cloud derived from a combination of Planck and Herschel dust-emission maps, and 2MASS NIR dust-extinction maps. We used these data to determine the ratio of the 2.2 micron extinction coefficient to the 850 micron opacity and found the value to be close to that found in similar studies of the Orion B and Perseus clouds but higher than that characterizing the Orion A cloud, indicating that variations in the fundamental optical properties of dust may exist between local clouds. We show that over a wide range of extinction, the column density probability distribution function (PDF$_N$) of the cloud can be well described by a simple power law with an index that represents a steeper decline with column density than found in similar studies of the Orion and Perseus clouds. Using only the protostellar population of the cloud and our extinction maps we investigate the Schmidt relation within the cloud. We show that the protostellar surface density, $Σ_*$, is directly proportional to the ratio of the protostellar and cloud pdfs. We use the cumulative distribution of protostars to infer the functional forms for both $Σ_*$ and PDF$_*$. We find that $Σ_*$ is best described by two power-law functions with steeper indicies than found in other local GMCs. We find that the protostellar pdf is a declining function of extinction also best described by two power-laws whose behavior mirrors that of $Σ_*$. Our observations suggest that variations both in the slope of the Schmidt relation and in the sizes of the protostellar populations between GMCs are largely driven by variations in the slope of the cloud pdf. This confirms earlier studies suggesting that cloud structure plays a major role in setting the global star formation rates in GMCs. △ Less

Submitted 26 September, 2017; v1 submitted 25 August, 2017; originally announced August 2017.

Comments: Accepted for publication in Astronomy and Astrophysics. Corrected typos in source coordinates in table A.1

Journal ref: A&A 606, A100 (2017)

A new method to unveil embedded stellar clusters

Authors: Marco Lombardi, Charles J. Lada, Joao Alves

Abstract: In this paper we present a novel method to identify and characterize stellar clusters deeply embedded in a dark molecular cloud. The method is based on measuring stellar surface density in wide-field infrared images using star counting techniques. It takes advantage of the differing $H$-band luminosity functions (HLFs) of field stars and young stellar populations and is able to statistically assoc… ▽ More In this paper we present a novel method to identify and characterize stellar clusters deeply embedded in a dark molecular cloud. The method is based on measuring stellar surface density in wide-field infrared images using star counting techniques. It takes advantage of the differing $H$-band luminosity functions (HLFs) of field stars and young stellar populations and is able to statistically associate each star in an image as a member of either the background stellar population or a young stellar population projected on or near the cloud. Moreover, the technique corrects for the effects of differential extinction toward each individual star. We have tested this method against simulations as well as observations. In particular, we have applied the method to 2MASS point sources observed in the Orion A and B complexes, and the results obtained compare very well with those obtained from deep Spitzer and Chandra observations where presence of infrared excess or X-ray emission directly determines membership status for every star. Additionally, our method also identifies unobscured clusters and a low resolution version of the Orion stellar surface density map shows clearly the relatively unobscured and diffuse OB 1a and 1b sub-groups and provides useful insights on their spatial distribution. △ Less

Submitted 26 July, 2017; originally announced July 2017.

Comments: A&A, in press; 13 pages, multi-layer figures can be displayed with Adobe Acrobat Reader

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

Protostars at Low Extinction in Orion A

Authors: John Arban Lewis, Charles J Lada

Abstract: In the list of young stellar objects compiled by Megeath et al. (2012) for the Orion A molecular cloud, only 44 out of 1208 sources found projected onto low extinction (Ak<0.8 mag) gas are identified as protostars. These objects are puzzling because protostars are not typically expected to be associated with extended low extinction material. Here, we use high resolution extinction maps generated f… ▽ More In the list of young stellar objects compiled by Megeath et al. (2012) for the Orion A molecular cloud, only 44 out of 1208 sources found projected onto low extinction (Ak<0.8 mag) gas are identified as protostars. These objects are puzzling because protostars are not typically expected to be associated with extended low extinction material. Here, we use high resolution extinction maps generated from Herschel data, optical/infrared and Spitzer Space Telescope photometry and spectroscopy of the low extinction protostellar candidate sources to determine if they are likely true protostellar sources or contaminants. Out of 44 candidate objects, we determine that 10 sources are likely protostars, with the rest being more evolved young stellar objects (18), galaxies (4), false detections of nebulosity and cloud edges (9), or real sources for which more data are required to ascertain their nature (3). We find none of the confirmed protostars to be associated with recognizable dense cores and we briefly discuss possible origins for these orphaned objects. △ Less

Submitted 11 May, 2016; v1 submitted 9 May, 2016; originally announced May 2016.

Comments: 26 pages, 59 figures. Accepted in the ApJ

Herschel-Planck dust optical depth and column density maps - II. Perseus

Authors: E. Zari, M. Lombardi, J. Alves, C. J. Lada, H. Bouy

Abstract: We present optical depth and temperature maps of the Perseus molecular cloud, obtained combining dust emission data from the Herschel and Planck satellites and 2MASS/NIR dust extinction maps. The maps have a resolution of 36 arcsec in the Herschel regions, and of 5 arcmin elsewhere. The dynamic range of the optical depth map ranges from $1\times10^{-2}\, \mathrm{mag}$ up to $20 \,\mathrm{mag}$ in… ▽ More We present optical depth and temperature maps of the Perseus molecular cloud, obtained combining dust emission data from the Herschel and Planck satellites and 2MASS/NIR dust extinction maps. The maps have a resolution of 36 arcsec in the Herschel regions, and of 5 arcmin elsewhere. The dynamic range of the optical depth map ranges from $1\times10^{-2}\, \mathrm{mag}$ up to $20 \,\mathrm{mag}$ in the equivalent K band extinction. We also evaluate the ratio between the $2.2 \,\mathrm{μm}$ extinction coefficient and the $850 \,\mathrm{μm}$ opacity. The value we obtain is close to the one found in the Orion B molecular cloud. We show that the cumulative and the differential area function of the data (which is proportional to the probability distribution function of the cloud column density) follow power laws with index respectively $\simeq -2$, and $\simeq -3$. We use WISE data to improve current YSO catalogues based mostly on \emph{Spitzer} data and we build an up-to-date selection of Class~I/0 objects. Using this selection, we evaluate the local Schmidt law, $Σ_{\mathrm{YSO}} \propto Σ_{\mathrm{gas}}^β$, showing that $β=2.4 \pm 0.6$. Finally, we show that the area-extinction relation is important for determining the star formation rate in the cloud, which is in agreement with other recent works. △ Less

Submitted 26 November, 2015; originally announced November 2015.

Comments: 15 pages, 19 figures. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 587, A106 (2016)

Star Formation in the Local Milky Way

Authors: Charles J. Lada

Abstract: Studies of molecular clouds and young stars near the sun have provided invaluable insights into the process of star formation. Indeed, much of our physical understanding of this topic has been derived from such studies. Perhaps the two most fundamental problems confronting star formation research today are: 1) determining the origin of stellar mass and 2) deciphering the nature of the physical pro… ▽ More Studies of molecular clouds and young stars near the sun have provided invaluable insights into the process of star formation. Indeed, much of our physical understanding of this topic has been derived from such studies. Perhaps the two most fundamental problems confronting star formation research today are: 1) determining the origin of stellar mass and 2) deciphering the nature of the physical processes that control the star formation rate in molecular gas. As I will briefly outline here, observations and studies of local star forming regions are making particularly significant contributions toward the solution of both these important problems. △ Less

Submitted 11 August, 2015; originally announced August 2015.

Comments: To appear in "Young Stars and Planets Near the Sun", Proceedings of IAU Symposium No. 314, J.H. Kastner, B. Stelzer, S.A. Metchev, eds

Smoke in the Pipe Nebula: dust emission and grain growth in the starless core FeSt 1-457

Authors: Jan Forbrich, Charles J. Lada, Marco Lombardi, Carlos Roman-Zuñiga, João Alves

Abstract: (abridged) Methods: We derive maps of submillimeter dust optical depth and effective dust temperature from Herschel data that were calibrated against Planck. After calibration, we then fit a modified blackbody to the long-wavelength Herschel data, using the Planck-derived dust opacity spectral index beta, derived on scales of 30' (or ~1 pc). We use this model to make predictions of the submillimet… ▽ More (abridged) Methods: We derive maps of submillimeter dust optical depth and effective dust temperature from Herschel data that were calibrated against Planck. After calibration, we then fit a modified blackbody to the long-wavelength Herschel data, using the Planck-derived dust opacity spectral index beta, derived on scales of 30' (or ~1 pc). We use this model to make predictions of the submillimeter flux density at 850 micron, and we compare these in turn with APEX-Laboca observations. Results: A comparison of the submillimeter dust optical depth and near-infrared extinction data reveals evidence for an increased submillimeter dust opacity at high column densities, interpreted as an indication of grain growth in the inner parts of the core. Additionally, a comparison of the Herschel dust model and the Laboca data reveals that the frequency dependence of the submillimeter opacity, described by the spectral index beta, does not change. A single beta that is only slightly different from the Planck-derived value is sufficient to describe the data, beta=1.53+/-0.07. We apply a similar analysis to Barnard 68, a core with significantly lower column densities than FeSt 1-457, and we do not find evidence for grain growth but also a single beta. Conclusions: While we find evidence for grain growth from the dust opacity in FeSt 1-457, we find no evidence for significant variations in the dust opacity spectral index beta on scales 0.02<x<1 pc (or 36"<x<30'). The correction to the Planck-derived dust beta that we find in both cases is on the order of the measurement error, not including any systematic errors, and it would thus be reasonable to directly apply the dust beta from the Planck all-sky dust model. As a corollary, reliable effective temperature maps can be derived which would be otherwise affected by beta variations. △ Less

Submitted 22 May, 2015; originally announced May 2015.

Comments: A&A, in press

Journal ref: A&A 580, A114 (2015)

The Relationship Between the Dust and Gas-Phase CO Across the California Molecular Cloud

Authors: S. Kong, C. J. Lada, E. A. Lada, C. Román-Zúñiga, J. H. Bieging, M. Lombardi, J. Forbrich, J. F. Alves

Abstract: A deep, wide-field, near-infrared imaging survey was used to construct an extinction map of the southeastern part of the California Molecular Cloud (CMC) with $\sim$ 0.5 arc min resolution. The same region was also surveyed in the $^{12}$CO(2-1), $^{13}$CO(2-1), C$^{18}$O(2-1) emission lines at the same angular resolution. Strong spatial variations in the abundances of $^{13}$CO and C$^{18}$O were… ▽ More A deep, wide-field, near-infrared imaging survey was used to construct an extinction map of the southeastern part of the California Molecular Cloud (CMC) with $\sim$ 0.5 arc min resolution. The same region was also surveyed in the $^{12}$CO(2-1), $^{13}$CO(2-1), C$^{18}$O(2-1) emission lines at the same angular resolution. Strong spatial variations in the abundances of $^{13}$CO and C$^{18}$O were found to be correlated with variations in gas temperature, consistent with temperature dependent CO depletion/desorption on dust grains. The $^{13}$CO to C$^{18}$O abundance ratio was found to increase with decreasing extinction, suggesting selective photodissociation of C$^{18}$O by the ambient UV radiation field. The cloud averaged X-factor is found to be $<$X$_{\rm CO}$$>$ $=$ 2.53 $\times$ 10$^{20}$ ${\rm cm}^{-2}~({\rm K~km~s}^{-1})^{-1}$, somewhat higher than the Milky Way average. On sub-parsec scales we find no single empirical value of the X-factor that can characterize the molecular gas in cold (T$_{\rm k}$ $\lesssim$ 15 K) regions, with X$_{\rm CO}$ $\propto$ A$_{\rm V}$$^{0.74}$ for A$_{\rm V}$ $\gtrsim$ 3 magnitudes. However in regions containing relatively hot (T$_{\rm ex}$ $\gtrsim$ 25 K) gas we find a clear correlation between W($^{12}$CO) and A$_{\rm V}$ over a large (3 $\lesssim$ A$_{\rm V}$ $\lesssim$ 25 mag) extinction range. This suggests a constant X$_{\rm CO}$ $=$ 1.5 $\times$ 10$^{20}$ ${\rm cm}^{-2}~({\rm K~km~s}^{-1})^{-1}$ for the hot gas, a lower value than either the average for the CMC or Milky Way. We find a correlation between X$_{\rm CO}$ and T$_{\rm ex}$ with X$_{\rm CO}$ $\propto$ T$_{\rm ex}$$^{-0.7}$ suggesting that the global X-factor of a cloud may depend on the relative amounts of hot gas within it. △ Less

Submitted 10 April, 2015; v1 submitted 11 March, 2015; originally announced March 2015.

Comments: 21 pages, 17 figures, accepted by ApJ

Molecular clouds have power-law probability distribution functions

Authors: Marco Lombardi, João Alves, Charles J. Lada

Abstract: In this Letter we investigate the shape of the probability distribution of column densities (PDF) in molecular clouds. Through the use of low-noise, extinction-calibrated \textit{Herschel}/\textit{Planck} emission data for eight molecular clouds, we demonstrate that, contrary to common belief, the PDFs of molecular clouds are not described well by log-normal functions, but are instead power laws w… ▽ More In this Letter we investigate the shape of the probability distribution of column densities (PDF) in molecular clouds. Through the use of low-noise, extinction-calibrated \textit{Herschel}/\textit{Planck} emission data for eight molecular clouds, we demonstrate that, contrary to common belief, the PDFs of molecular clouds are not described well by log-normal functions, but are instead power laws with exponents close to two and with breaks between $A_K \simeq 0.1$ and $0.2\,\mathrm{mag}$, so close to the CO self-shielding limit and not far from the transition between molecular and atomic gas. Additionally, we argue that the intrinsic functional form of the PDF cannot be securely determined below $A_K \simeq 0.1\,\mathrm{mag}$, limiting our ability to investigate more complex models for the shape of the cloud PDF. △ Less

Submitted 12 February, 2015; originally announced February 2015.

Comments: Letter to the Editor, to appear in A&A

On Schmidt's Conjecture and Star Formation Scaling Laws

Authors: Charles J. Lada

Abstract: Ever since the pioneering work of Schmidt a half-century ago there has been great interest in finding an appropriate empirical relation that would directly link some property of interstellar gas with the process of star formation within it. Schmidt conjectured that this might take the form of a power-law relation between the rate of star formation (SFR) and the surface density of interstellar gas.… ▽ More Ever since the pioneering work of Schmidt a half-century ago there has been great interest in finding an appropriate empirical relation that would directly link some property of interstellar gas with the process of star formation within it. Schmidt conjectured that this might take the form of a power-law relation between the rate of star formation (SFR) and the surface density of interstellar gas. However, recent observations suggest that a linear scaling relation between the total SFR and the amount of dense gas within molecular clouds appears to be the underlying physical relation that most directly connects star formation with interstellar gas from scales of individual GMCs to those encompassing entire galaxies both near and far. Although Schmidt relations are found to exist within local GMCs, there is no Schmidt relation observed between GMCs. The implications of these results for interpreting and understanding the Kennicutt-Schmidt scaling law for galaxies are discussed. △ Less

Submitted 4 November, 2014; originally announced November 2014.

Comments: To appear in Proceedings of IAU Symposium 309: Galaxies in 3D across the Universe eds. B.L. Ziegler, F. Combes, H. Dannerbauer, M. Verdugo

Some like it cold: molecular emission and effective dust temperatures of dense cores in the Pipe Nebula

Authors: Jan Forbrich, Karin Öberg, Charles J. Lada, Marco Lombardi, Alvaro Hacar, João Alves, Jill M. Rathborne

Abstract: (abridged) [...] Methods: In a continued study of the molecular core population of the Pipe Nebula, we present a molecular-line survey of 52 cores. Previous research has shown a variety of different chemical evolutionary stages among the cores. Using the Mopra radio telescope, we observed the ground rotational transitions of HCO+, H13CO+, HCN, H13CN, HNC, and N2H+. These data are complemented with… ▽ More (abridged) [...] Methods: In a continued study of the molecular core population of the Pipe Nebula, we present a molecular-line survey of 52 cores. Previous research has shown a variety of different chemical evolutionary stages among the cores. Using the Mopra radio telescope, we observed the ground rotational transitions of HCO+, H13CO+, HCN, H13CN, HNC, and N2H+. These data are complemented with near-infrared extinction maps to constrain the column densities, effective dust temperatures derived from Herschel data, and NH3-based gas kinetic temperatures. Results: The target cores are located across the nebula, span visual extinctions between 5 and 67 mag, and effective dust temperatures (averaged along the lines of sight) between 13 and 19 K. The extinction-normalized integrated line intensities, a proxy for the abundance in constant excitation conditions of optically thin lines, vary within an order of magnitude for a given molecule. The effective dust temperatures and gas kinetic temperatures are correlated, but the effective dust temperatures are consistently higher than the gas kinetic temperatures. Combining the molecular line and temperature data, we find that N2H+ is only detected toward the coldest and densest cores while other lines show no correlation with these core properties. Conclusions: Within this large sample, N2H+ is the only species to exclusively trace the coldest and densest cores, in agreement with chemical considerations. In contrast, the common high-density tracers HCN and HNC are present in a majority of cores, demonstrating the utility of these molecules to characterize cores over a large range of extinctions. The correlation between the effective dust temperatures and the gas kinetic temperatures suggests that the former are dominated by dust that is both dense and thermodynamically coupled to the dense gas traced by NH3. [...] △ Less

Submitted 2 June, 2014; originally announced June 2014.

Comments: A&A, in press

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

Molecular Cloud-scale Star Formation in NGC 300

Authors: Christopher M. Faesi, Charles J. Lada, Jan Forbrich, Karl Menten, Herve Bouy

Abstract: We present the results of a galaxy-wide study of molecular gas and star formation in a sample of 76 HII regions in the nearby spiral galaxy NGC 300. We have measured the molecular gas at 250 pc scales using pointed CO(J=2-1) observations with the APEX telescope. We detect CO in 42 of our targets, deriving molecular gas masses ranging from our sensitivity limit of ~10^5 Msun to 7x10^5 Msun. We find… ▽ More We present the results of a galaxy-wide study of molecular gas and star formation in a sample of 76 HII regions in the nearby spiral galaxy NGC 300. We have measured the molecular gas at 250 pc scales using pointed CO(J=2-1) observations with the APEX telescope. We detect CO in 42 of our targets, deriving molecular gas masses ranging from our sensitivity limit of ~10^5 Msun to 7x10^5 Msun. We find a clear decline in the CO detection rate with galactocentric distance, which we attribute primarily to the decreasing radial metallicity gradient in NGC 300. We combine GALEX FUV, Spitzer 24 micron, and H-alpha narrowband imaging to measure the star formation activity in our sample. We have developed a new direct modeling approach for computing star formation rates that utilizes these data and population synthesis models to derive the masses and ages of the young stellar clusters associated with each of our HII region targets. We find a characteristic gas depletion time of 230 Myr at 250 pc scales in NGC 300, more similar to the results obtained for Milky Way Giant Molecular Clouds than the longer (>2 Gyr) global depletion times derived for entire galaxies and kpc-sized regions within them. This difference is partially due to the fact that our study accounts for only the gas and stars within the youngest star forming regions. We also note a large scatter in the NGC 300 SFR-molecular gas mass scaling relation that is furthermore consistent with the Milky Way cloud results. This scatter likely represents real differences in giant molecular cloud physical properties such as the dense gas fraction. △ Less

Submitted 9 May, 2014; originally announced May 2014.

Comments: 42 pages, 14 figures, 7 tables. Includes a complete image atlas of our HII region targets. ASCII versions of tables will be available electronically after paper is published. Accepted for publication in the Astrophysical Journal

Herschel-Planck dust optical-depth and column-density maps: I. Method description and results for Orion

Authors: Marco Lombardi, Herve Bouy, Joao Alves, Charles J. Lada

Abstract: We present high-resolution, high dynamic range column-density and color-temperature maps of the Orion complex using a combination of Planck dust-emission maps, Herschel dust-emission maps, and 2MASS NIR dust-extinction maps. The column-density maps combine the robustness of the 2MASS NIR extinction maps with the resolution and coverage of the Herschel and Planck dust-emission maps and constitute t… ▽ More We present high-resolution, high dynamic range column-density and color-temperature maps of the Orion complex using a combination of Planck dust-emission maps, Herschel dust-emission maps, and 2MASS NIR dust-extinction maps. The column-density maps combine the robustness of the 2MASS NIR extinction maps with the resolution and coverage of the Herschel and Planck dust-emission maps and constitute the highest dynamic range column-density maps ever constructed for the entire Orion complex, covering $0.01 \, \mathrm{mag} < A_K < 30 \,\mathrm{mag}$, or $2 \times 10^{20} \, \mathrm{cm}^{-2} < N < 5 \times 10^{23} \,\mathrm{cm}^{-2}$. We determined the ratio of the 2.2 microns extinction coefficient to the 850 microns opacity and found that the values obtained for both Orion A and B are significantly lower than the predictions of standard dust models, but agree with newer models that incorporate icy silicate-graphite conglomerates for the grain population. We show that the cloud projected pdf, over a large range of column densities, can be well fitted by a simple power law. Moreover, we considered the local Schmidt-law for star formation, and confirm earlier results, showing that the protostar surface density $Σ_*$ follows a simple law $Σ_* \propto Σ_{gas}^β$, with $β\sim 2$. △ Less

Submitted 2 April, 2014; v1 submitted 31 March, 2014; originally announced April 2014.

Comments: 19 pages, A&A in press. High resolution figures available at http://www.marcolombardi.org/research/recent-papers/orion

Fitting density models to observational data - The local Schmidt law in molecular clouds

Authors: Marco Lombardi, Charles J. Lada, João Alves

Abstract: We consider the general problem of fitting a parametric density model to discrete observations, taken to follow a non-homogeneous Poisson point process. This class of models is very common, and can be used to describe many astrophysical processes, including the distribution of protostars in molecular clouds. We give the expression for the likelihood of a given spatial density distribution of proto… ▽ More We consider the general problem of fitting a parametric density model to discrete observations, taken to follow a non-homogeneous Poisson point process. This class of models is very common, and can be used to describe many astrophysical processes, including the distribution of protostars in molecular clouds. We give the expression for the likelihood of a given spatial density distribution of protostars and apply it to infer the most probable dependence of the protostellar surface density on the gas surface density. Finally, we apply this general technique to model the distribution of protostars in the Orion molecular cloud and robustly derive the local star formation scaling (Schmidt) law for a molecular cloud. We find that in this cloud the protostellar surface density, $Σ_\mathrm{YSO}$, is directly proportional to the square gas column density, here expressed as infrared extinction in the $K$-band, $A_K$: more precisely, $Σ_\mathrm{YSO} = (1.65 \pm 0.19) A_K^{(2.03 \pm 0.15)}$ stars pc$^{-2}$. △ Less

Submitted 15 October, 2013; v1 submitted 14 October, 2013; originally announced October 2013.

Comments: 6 pages, A&A in press. References fixed

The mid-infrared extinction law in the darkest cores of the Pipe Nebula

Authors: J. Ascenso, C. J. Lada, J. Alves, C. G. Román-Zúñiga, M. Lombardi

Abstract: Context. The properties of dust grains, in particular their size distribution, are expected to differ from the interstellar medium to the high-density regions within molecular clouds. Aims. We measure the mid-infrared extinction law produced by dense material in molecular cloud cores. Since the extinction at these wavelengths is caused by dust, the extinction law in cores should depart from that f… ▽ More Context. The properties of dust grains, in particular their size distribution, are expected to differ from the interstellar medium to the high-density regions within molecular clouds. Aims. We measure the mid-infrared extinction law produced by dense material in molecular cloud cores. Since the extinction at these wavelengths is caused by dust, the extinction law in cores should depart from that found in low-density environments if the dust grains have different properties. Methods. We use the unbiased LINES method to measure the slope of the reddening vectors in color-color diagrams. We derive the mid-infrared extinction law toward the dense cores B59 and FeSt 1-457 in the Pipe Nebula over a range of visual extinction between 10 and 50 magnitudes, using a combination of Spitzer/IRAC, and ESO NTT/VLT data. Results. The mid-infrared extinction law in both cores departs significantly from a power-law between 3.6 and 8 micron, suggesting that these cores contain dust with a considerable fraction of large dust grains. We find no evidence for a dependence of the extinction law with column density up to 50 magnitudes of visual extinction in these cores, and no evidence for a variation between our result and those for other clouds at lower column densities reported elsewhere in the literature. This suggests that either large grains are present even in low column density regions, or that the existing dust models need to be revised at mid-infrared wavelengths. We find a small but significant difference in the extinction law of the two cores, that we tentatively associate with the onset of star formation in B59. △ Less

Submitted 28 November, 2012; originally announced November 2012.

Comments: 8 pages, 6 figures. Accepted to A&A

The extinction law from photometric data: linear regression methods

Authors: Joana Ascenso, Marco Lombardi, Charles J. Lada, João Alves

Abstract: Context. The properties of dust grains, in particular their size distribution, are expected to differ from the interstellar medium to the high-density regions within molecular clouds. Since the extinction at near-infrared wavelengths is caused by dust, the extinction law in cores should depart from that found in low-density environments if the dust grains have different properties. Aims. We explor… ▽ More Context. The properties of dust grains, in particular their size distribution, are expected to differ from the interstellar medium to the high-density regions within molecular clouds. Since the extinction at near-infrared wavelengths is caused by dust, the extinction law in cores should depart from that found in low-density environments if the dust grains have different properties. Aims. We explore methods to measure the near-infrared extinction law produced by dense material in molecular cloud cores from photometric data. Methods. Using controlled sets of synthetic and semi-synthetic data, we test several methods for linear regression applied to the specific problem of deriving the extinction law from photometric data. We cover the parameter space appropriate to this type of observations. Results. We find that many of the common linear-regression methods produce biased results when applied to the extinction law from photometric colors. We propose and validate a new method, LinES, as the most reliable for this effect. We explore the use of this method to detect whether or not the extinction law of a given reddened population has a break at some value of extinction. △ Less

Submitted 23 March, 2012; originally announced March 2012.

Comments: 15 pages, 18 figures, accepted to A&A, in press

Spitzer observations of NGC2264: The nature of the disk population

Authors: P. S. Teixeira, C. J. Lada, M. Marengo, E. A. Lada

Abstract: NGC2264 is a young cluster with a rich circumstellar disk population which makes it an ideal target for studying the evolution of stellar clusters. Our goal is to study its star formation history and to analyse the primordial disk evolution of its members. The study presented is based on data obtained with Spitzer IRAC and MIPS, combined with deep NIR ground-based FLAMINGOS imaging and previously… ▽ More NGC2264 is a young cluster with a rich circumstellar disk population which makes it an ideal target for studying the evolution of stellar clusters. Our goal is to study its star formation history and to analyse the primordial disk evolution of its members. The study presented is based on data obtained with Spitzer IRAC and MIPS, combined with deep NIR ground-based FLAMINGOS imaging and previously published optical data. We build NIR dust extinction maps of the molecular cloud associated with the cluster, and determine it to have a mass of 2.1x10^3Msun above an Av of 7mag. Using a differential K_s-band luminosity function of the cluster, we estimate the size of its population to be 1436$\pm$242 members. The star formation efficiency is ~25%. We identify the disk population: (i) optically thick inner disks, (ii) anaemic inner disks, and (iii) disks with inner holes, or transition disks. We analyse the spatial distribution of these sources and find that sources with thick disks segregate into sub-clusterings, whereas sources with anaemic disks do not. Furthermore, sources with anaemic disks are found to be unembedded (Av<3mag), whereas the clustered sources with thick disks are still embedded within the parental cloud. NGC2264 has undergone more than one star-forming event, where the anaemic and extincted thick disk population appear to have formed in separate episodes. We also find tentative evidence of triggered star-formation in the Fox Fur Nebula. In terms of disk evolution, our findings support the emerging disk evolution paradigm of two distinct evolutionary paths for primordial optically thick disks: a homologous one where the disk emission decreases uniformly at NIR and MIR wavelengths, and a radially differential one where the emission from the inner region of the disk decreases more rapidly than from the outer region (forming transition disks). △ Less

Submitted 16 March, 2012; originally announced March 2012.

Comments: accepted for publishing in A&A

Star Formation Rates in Molecular Clouds and the Nature of the Extragalactic Scaling Relations

Authors: Charles J. Lada, Jan Forbrich, Marco Lombardi, Joao F. Alves

Abstract: In this paper we investigate scaling relations between star formation rates and molecular gas masses for both local Galactic clouds and a sample of external galaxies. We specifically consider relations between the star formation rates and measurements of dense, as well as total, molecular gas masses. We argue that there is a fundamental empirical scaling relation that directly connects the local s… ▽ More In this paper we investigate scaling relations between star formation rates and molecular gas masses for both local Galactic clouds and a sample of external galaxies. We specifically consider relations between the star formation rates and measurements of dense, as well as total, molecular gas masses. We argue that there is a fundamental empirical scaling relation that directly connects the local star formation process with that operating globally within galaxies. Specifically, the total star formation rate in a molecular cloud or galaxy is linearly proportional to the mass of dense gas within the cloud or galaxy. This simple relation, first documented in previous studies, holds over a span of mass covering nearly nine orders of magnitude and indicates that the rate of star formation is directly controlled by the amount of dense molecular gas that can be assembled within a star formation complex. We further show that the star formation rates and total molecular masses, characterizing both local clouds and galaxies, are correlated over similarly large scales of mass and can be described by a family of linear star formation scaling laws, parameterized by $f_{DG}$, the fraction of dense gas contained within the clouds or galaxies. That is, the underlying star formation scaling law is always linear for clouds and galaxies with the same dense gas fraction. These considerations provide a single unified framework for understanding the relation between the standard (non-linear) extragalactic Schmidt-Kennicutt scaling law, that is typically derived from CO observations of the gas, and the linear star formation scaling law derived from HCN observations of the dense gas. △ Less

Submitted 19 December, 2011; originally announced December 2011.

Comments: 14 pages + 2 figures. Accepted for publication in ApJ 16 December 2011

Ices in the Quiescent IC 5146 Dense Cloud

Authors: J. E. Chiar, Y. J. Pendleton, L. J. Allamandola, A. C. A. Boogert, K. Ennico, T. P. Greene, T. R. Geballe, J. V. Keane, C. J. Lada, R. E. Mason, T. L. Roellig, S. A. Sandford, A. G. G. M. Tielens, M. W. Werner, D. C. B. Whittet, L. Decin, K. Eriksson

Abstract: This paper presents spectra in the 2 to 20 micron range of quiescent cloud material located in the IC 5146 cloud complex. The spectra were obtained with NASA's Infrared Telescope Facility (IRTF) SpeX instrument and the Spitzer Space Telescope's Infrared Spectrometer. We use these spectra to investigate dust and ice absorption features in pristine regions of the cloud that are unaltered by embedded… ▽ More This paper presents spectra in the 2 to 20 micron range of quiescent cloud material located in the IC 5146 cloud complex. The spectra were obtained with NASA's Infrared Telescope Facility (IRTF) SpeX instrument and the Spitzer Space Telescope's Infrared Spectrometer. We use these spectra to investigate dust and ice absorption features in pristine regions of the cloud that are unaltered by embedded stars. We find that the H2O-ice threshold extinction is 4.03+/-0.05 mag. Once foreground extinction is taken into account, however, the threshold drops to 3.2 mag, equivalent to that found for the Taurus dark cloud, generally assumed to be the touchstone quiescent cloud against which all other dense cloud and embedded young stellar object observations are compared. Substructure in the trough of the silicate band for two sources is attributed to CH3OH and NH3 in the ices, present at the ~2% and ~5% levels, respectively, relative to H2O-ice. The correlation of the silicate feature with the E(J-K) color excess is found to follow a much shallower slope relative to lines of sight that probe diffuse clouds, supporting the previous results by Chiar et al. (2007). △ Less

Submitted 12 February, 2011; originally announced February 2011.

Comments: 13 pages, 13 figures with multiple parts, accepted for publication in Astrophysical Journal, Feb. 2011

Deep Near-Infrared Survey of the Pipe Nebula II: Data, Methods, and Dust Extinction Maps

Authors: Carlos Gerardo Román-Zúñiga, João F. Alves, Charles J. Lada, Marco Lombardi

Abstract: We present a new set of high resolution dust extinction maps of the nearby and essentially starless Pipe Nebula molecular cloud. The maps were constructed from a concerted deep near-infrared imaging survey with the ESO-VLT, ESO-NTT, CAHA 3.5m telescopes, and 2MASS data. The new maps have a resolution three times higher than the previous extinction map of this cloud by Lombardi et al. (2006) and ar… ▽ More We present a new set of high resolution dust extinction maps of the nearby and essentially starless Pipe Nebula molecular cloud. The maps were constructed from a concerted deep near-infrared imaging survey with the ESO-VLT, ESO-NTT, CAHA 3.5m telescopes, and 2MASS data. The new maps have a resolution three times higher than the previous extinction map of this cloud by Lombardi et al. (2006) and are able to resolve structure down to 2600 AU. We detect 244 significant extinction peaks across the cloud. These peaks have masses between 0.1 and 18.4 M_sun, diameters between 1.2 and 5.7e4 AU (0.06 and 0.28 pc), and mean densities of about 1e4 cm${^-3}$, all in good agreement with previous results. From the analysis of the Mean Surface Density of Companions we find a well defined scale near 1.4e4 AU below which we detect a significant decrease in structure of the cloud. This scale is smaller than the Jeans Length calculated from the mean density of the peaks. The surface density of peaks is not uniform but instead it displays clustering. Extinction peaks in the Pipe Nebula appear to have a spatial distribution similar to the stars in Taurus, suggesting that the spatial distribution of stars evolves directly from the primordial spatial distribution of high density material. △ Less

Submitted 5 November, 2010; originally announced November 2010.

Comments: Accepted for Publication in The Astrophysical Journal, 10/10. 49 pages, 16 figures

On the Star Formation Rates in Molecular Clouds

Authors: Charles J. Lada, Marco Lombardi, João F. Alves

Abstract: In this paper we investigate the level of star formation activity within nearby molecular clouds. We employ a uniform set of infrared extinction maps to provide accurate assessments of cloud mass and structure and compare these with inventories of young stellar objects within the clouds. We present evidence indicating that both the yield and rate of star formation can vary considerably in local cl… ▽ More In this paper we investigate the level of star formation activity within nearby molecular clouds. We employ a uniform set of infrared extinction maps to provide accurate assessments of cloud mass and structure and compare these with inventories of young stellar objects within the clouds. We present evidence indicating that both the yield and rate of star formation can vary considerably in local clouds, independent of their mass and size. We find that the surface density structure of such clouds appears to be important in controlling both these factors. In particular, we find that the star formation rate (SFR) in molecular clouds is linearly proportional to the cloud mass (M_{0.8}) above an extinction threshold of A_K approximately equal to 0.8 magnitudes, corresponding to a gas surface density threshold of approximaely 116 solar masses per square pc. We argue that this surface density threshold corresponds to a gas volume density threshold which we estimate to be n(H_2) approximately equal to 10^4\cc. Specifically we find SFR (solar masses per yr) = 4.6 +/- 2.6 x 10^{-8} M_{0.8} (solar masses) for the clouds in our sample. This relation between the rate of star formation and the amount of dense gas in molecular clouds appears to be in excellent agreement with previous observations of both galactic and extragalactic star forming activity. It is likely the underlying physical relationship or empirical law that most directly connects star formation activity with interstellar gas over many spatial scales within and between individual galaxies. These results suggest that the key to obtaining a predictive understanding of the star formation rates in molecular clouds and galaxies is to understand those physical factors which give rise to the dense components of these clouds. △ Less

Submitted 15 September, 2010; originally announced September 2010.

Comments: accepted for publicaton in the Astrophysical Journal; 22 pages, 4 figures

Larson's third law and the universality of molecular cloud structure

Authors: Marco Lombardi, Joao Alves, Charles J. Lada

Abstract: Larson (1981) first noted a scaling relation between masses and sizes in molecular clouds that implies that these objects have approximately constant column densities. This original claim, based upon millimeter observations of carbon monoxide lines, has been challenged by many theorists, arguing that the apparent constant column density observed is merely the result of the limited dynamic range of… ▽ More Larson (1981) first noted a scaling relation between masses and sizes in molecular clouds that implies that these objects have approximately constant column densities. This original claim, based upon millimeter observations of carbon monoxide lines, has been challenged by many theorists, arguing that the apparent constant column density observed is merely the result of the limited dynamic range of observations, and that in reality clouds have column density variations over two orders of magnitudes. In this letter we investigate a set of nearby molecular clouds with near-infrared excess methods, which guarantee very large dynamic ranges and robust column density measurements, to test the validity of Larson's third law. We verify that different clouds have almost identical average column densities above a given extinction threshold; this holds regardless of the extinction threshold, but the actual average surface mass density is a function of the specific threshold used. We show that a second version of Larson's third law, involving the mass-radius relation for single clouds and cores, does not hold in our sample, indicating that individual clouds are not objects that can be described by constant column density. Our results instead indicate that molecular clouds are characterized by a universal structure. Finally we point out that this universal structure can be linked to the log-normal nature of cloud column density distributions. △ Less

Submitted 25 August, 2010; originally announced August 2010.

Comments: 5 pages, 4 figures, A&A in press (letter)

The Age, Stellar Content and Star Formation Timescale of the B59 Dense Core

Authors: Kevin R. Covey, Charles J. Lada, Carlos Roman-Zuniga, August A. Muench, Jan Forbrich, Joana Ascenso

Abstract: We have used moderate resolution, near-infrared spectra from the SpeX spectrograph on the NASA Infrared Telescope facility to characterize the stellar content of Barnard 59 (B59), the most active star-forming core in the Pipe Nebula. Measuring luminosity and temperature sensitive features in the spectra of 20 candidate YSOs, we identified likely background giant stars and measured each star's spec… ▽ More We have used moderate resolution, near-infrared spectra from the SpeX spectrograph on the NASA Infrared Telescope facility to characterize the stellar content of Barnard 59 (B59), the most active star-forming core in the Pipe Nebula. Measuring luminosity and temperature sensitive features in the spectra of 20 candidate YSOs, we identified likely background giant stars and measured each star's spectral type, extinction, and NIR continuum excess. We find that B59 is composed of late type (K4-M6) low-mass (0.9--0.1 M_sun) YSOs whose median stellar age is comparable to, if not slightly older than, that of YSOs within the Rho Oph, Taurus, and Chameleon star forming regions. Deriving absolute age estimates from pre-main sequence models computed by D'Antona et al., and accounting only for statistical uncertainties, we measure B59's median stellar age to be 2.6+/-0.8 Myrs. Including potential systematic effects increases the error budget for B59's median (DM98) stellar age to 2.6+4.1/-2.6 Myrs. We also find that the relative age orderings implied by pre-main sequence evolutionary tracks depend on the range of stellar masses sampled, as model isochrones possess significantly different mass dependencies. The maximum likelihood median stellar age we measure for B59, and the region's observed gas properties, suggest that the B59 dense core has been stable against global collapse for roughly 6 dynamical timescales, and is actively forming stars with a star formation efficiency per dynamical time of ~6%. This maximum likelihood value agrees well with recent star formation simulations that incorporate various forms of support against collapse, such as sub-critical magnetic fields, outflows, and radiative feedback from protostellar heating. [abridged] △ Less

Submitted 15 July, 2010; v1 submitted 13 July, 2010; originally announced July 2010.

Comments: 20 pages, 12 figures, accepted for publication in the Astrophysical Journal; updated to amend acknowledgments

Nothing to hide: An X-ray survey for young stellar objects in the Pipe Nebula

Authors: Jan Forbrich, Bettina Posselt, Kevin R. Covey, Charles J. Lada

Abstract: We have previously analyzed sensitive mid-infrared observations to establish that the Pipe Nebula has a very low star-formation efficiency. That study focused on YSOs with excess infrared emission (i.e, protostars and pre-main sequence stars with disks), however, and could have missed a population of more evolved pre-main sequence stars or Class III objects (i.e., young stars with dissipated disks… ▽ More We have previously analyzed sensitive mid-infrared observations to establish that the Pipe Nebula has a very low star-formation efficiency. That study focused on YSOs with excess infrared emission (i.e, protostars and pre-main sequence stars with disks), however, and could have missed a population of more evolved pre-main sequence stars or Class III objects (i.e., young stars with dissipated disks that no longer show excess infrared emission). Evolved pre-main sequence stars are X-ray bright, so we have used ROSAT All-Sky Survey data to search for diskless pre-main sequence stars throughout the Pipe Nebula. We have also analyzed archival XMM-Newton observations of three prominent areas within the Pipe: Barnard 59, containing a known cluster of young stellar objects; Barnard 68, a dense core that has yet to form stars; and the Pipe molecular ring, a high-extinction region in the bowl of the Pipe. We additionally characterize the X-ray properties of YSOs in Barnard 59. The ROSAT and XMM-Newton data provide no indication of a significant population of more evolved pre-main sequence stars within the Pipe, reinforcing our previous measurement of the Pipe's very low star formation efficiency. △ Less

Submitted 17 June, 2010; originally announced June 2010.

Comments: Accepted for publication in ApJ

Disentangling protostellar evolutionary stages in clustered environments using Spitzer-IRS spectra and comprehensive SED modeling

Authors: Jan Forbrich, Achim Tappe, Thomas Robitaille, August A. Muench, Paula S. Teixeira, Elizabeth A. Lada, Andrea Stolte, Charles J. Lada

Abstract: When studying the evolutionary stages of protostars that form in clusters, the role of any intracluster medium cannot be neglected. High foreground extinction can lead to situations where young stellar objects (YSOs) appear to be in earlier evolutionary stages than they actually are, particularly when using simple criteria like spectral indices. To address this issue, we have assembled detailed SE… ▽ More When studying the evolutionary stages of protostars that form in clusters, the role of any intracluster medium cannot be neglected. High foreground extinction can lead to situations where young stellar objects (YSOs) appear to be in earlier evolutionary stages than they actually are, particularly when using simple criteria like spectral indices. To address this issue, we have assembled detailed SED characterizations of a sample of 56 Spitzer-identified candidate YSOs in the clusters NGC 2264 and IC 348. For these, we use spectra obtained with the Infrared Spectrograph onboard the Spitzer Space Telescope and ancillary multi-wavelength photometry. The primary aim is twofold: 1) to discuss the role of spectral features, particularly those due to ices and silicates, in determining a YSO's evolutionary stage, and 2) to perform comprehensive modeling of spectral energy distributions (SEDs) enhanced by the IRS data. The SEDs consist of ancillary optical-to-submillimeter multi-wavelength data as well as an accurate description of the 9.7 micron silicate feature and of the mid-infrared continuum derived from line-free parts of the IRS spectra. We find that using this approach, we can distinguish genuine protostars in the cluster from T Tauri stars masquerading as protostars due to external foreground extinction. Our results underline the importance of photometric data in the far-infrared/submillimeter wavelength range, at sufficiently high angular resolution to more accurately classify cluster members. Such observations are becoming possible now with the advent of the Herschel Space Observatory. △ Less

Submitted 11 May, 2010; originally announced May 2010.

Comments: Accepted for publication in ApJ

The physics and modes of star cluster formation: observations

Authors: Charles J. Lada

Abstract: Stellar clusters are born in cold and dusty molecular clouds and the youngest clusters are embedded to various degrees in dusty dark molecular material. Such embedded clusters can be considered protocluster systems. The most deeply buried examples are so heavily obscured by dust that they are only visible at infrared wavelengths. These embedded protoclusters constitute the nearest laboratories f… ▽ More Stellar clusters are born in cold and dusty molecular clouds and the youngest clusters are embedded to various degrees in dusty dark molecular material. Such embedded clusters can be considered protocluster systems. The most deeply buried examples are so heavily obscured by dust that they are only visible at infrared wavelengths. These embedded protoclusters constitute the nearest laboratories for direct astronomical investigation of the physical processes of cluster formation and early evolution. I review the present state of empirical knowledge concerning embedded cluster systems and discuss the implications for understanding their formation and subsequent evolution to produce bound stellar clusters. △ Less

Submitted 4 November, 2009; originally announced November 2009.

Comments: 18 pages, 7 figures. To appear as invited review article in a special issue of the Phil. Trans. Royal Soc. A: Ch. 2 "Star clusters as tracers of galactic star-formation histories" (ed. R. de Grijs). Fully peer reviewed. LaTeX, requires rspublic.cls style file

A Spitzer Census of Star Formation Activity in the Pipe Nebula

Authors: Jan Forbrich, Charles J. Lada, August A. Muench, João Alves, Marco Lombardi

Abstract: The Pipe Nebula, a large nearby molecular cloud lacks obvious signposts of star formation in all but one of more than 130 dust extinction cores that have been identified within it. In order to quantitatively determine the current level of star formation activity in the Pipe Nebula, we analyzed 13 square degrees of sensitive mid-infrared maps of the entire cloud, obtained with the Multiband Imagi… ▽ More The Pipe Nebula, a large nearby molecular cloud lacks obvious signposts of star formation in all but one of more than 130 dust extinction cores that have been identified within it. In order to quantitatively determine the current level of star formation activity in the Pipe Nebula, we analyzed 13 square degrees of sensitive mid-infrared maps of the entire cloud, obtained with the Multiband Imaging Photometer for Spitzer (MIPS) at wavelengths of 24 micron and 70 micron to search for candidate Young Stellar Objects (YSOs) in the high-extinction regions. We argue that our search is complete for class I and typical class II YSOs with luminosities of Lbol~0.2 Lo and greater. We find only 18 candidate YSOs in the high-extinction regions of the entire Pipe cloud. Twelve of these sources are previously known members of a small cluster associated with Barnard 59, the largest and most massive dense core in the cloud. With only six candidate class I and class II YSOs detected towards extinction cores outside of this cluster, our findings emphatically confirm the notion of an extremely low level of star formation activity in the Pipe Nebula. The resulting star formation efficiency for the entire cloud mass is only ~0.06 %. △ Less

Submitted 27 August, 2009; originally announced August 2009.

Comments: 23 pages, 8 figures, accepted for publication in ApJ

High Resolution Near-Infrared Survey of the Pipe Nebula I: A Deep Infrared Extinction Map of Barnard 59

Authors: Carlos G. Román-Zúñiga, Charles J. Lada, João F. Alves

Abstract: We present our analysis of a fully sampled, high resolution dust extinction map of the Barnard 59 complex in the Pipe Nebula. The map was constructed with the infrared color excess technique applied to a photometric catalog that combines data from both ground and space based observations. The map resolves for the first time the high density center of the main core in the complex, that is associa… ▽ More We present our analysis of a fully sampled, high resolution dust extinction map of the Barnard 59 complex in the Pipe Nebula. The map was constructed with the infrared color excess technique applied to a photometric catalog that combines data from both ground and space based observations. The map resolves for the first time the high density center of the main core in the complex, that is associated with the formation of a small cluster of stars. We found that the central core in Barnard 59 shows an unexpected lack of significant substructure consisting of only two significant fragments. Overall, the material appears to be consistent with being a single, large core with a density profile that can be well fit by a King model. A series of NH$_3$ pointed observations towards the high column density center of the core appear to show that the core is still thermally dominated, with sub-sonic non-thermal motions. The stars in the cluster could be providing feedback to support the core against collapse, but the relatively narrow radio lines suggest that an additional source of support, for example a magnetic field, may be required to stabilize the core. Outside the central core our observations reveal the structure of peripheral cores and resolve an extended filament into a handful of significant substructures whose spacing and masses appear to be consistent with Jeans fragmentation. △ Less

Submitted 25 August, 2009; originally announced August 2009.

Comments: Accepted to ApJ, DOP October 2009. 14 pages, 11 figures, 3 tables

The California Molecular Cloud

Authors: Charles J. Lada, Marco Lombardi, Joao F. Alves

Abstract: We present an analysis of wide-field infrared extinction maps of a region in Perseus just north of the Taurus-Auriga dark cloud complex. From this analysis we have identified a massive, nearby, but previously unrecognized, giant molecular cloud (GMC). From comparison of foreground star counts with Galactic models we derive a distance of 450 +/- 23 parsecs to the cloud. At this distance the cloud… ▽ More We present an analysis of wide-field infrared extinction maps of a region in Perseus just north of the Taurus-Auriga dark cloud complex. From this analysis we have identified a massive, nearby, but previously unrecognized, giant molecular cloud (GMC). From comparison of foreground star counts with Galactic models we derive a distance of 450 +/- 23 parsecs to the cloud. At this distance the cloud extends over roughly 80 pc and has a mass of approximately 10^5 solar masses, rivaling the Orion (A) Molecular Cloud as the largest and most massive GMC in the solar neighborhood. Although surprisingly similar in mass and size to the more famous Orion Molecular Cloud (OMC) the newly recognized cloud displays significantly less star formation activity with more than an order of magnitude fewer young stellar objects than found in the OMC, suggesting that both the level of star formation and perhaps the star formation rate in this cloud are an order of magnitude or more lower than in the OMC. Analysis of extinction maps of both clouds shows that the new cloud contains only 10% the amount of high extinction (A_K > 1.0 mag) material as is found in the OMC. This, in turn, suggests that the level of star formation activity and perhaps the star formation rate in these two clouds may be directly proportional to the total amount of high extinction material and presumably high density gas within them and that there might be a density threshold for star formation on the order of a few times 10^4 hydrogen molecules per cc. △ Less

Submitted 5 August, 2009; originally announced August 2009.

Comments: Accepted for publication in the ApJ. This submission has low resolution figures. Preprint with full resolution figures can be found at http://www.cfa.harvard.edu/~clada/pubs_html/California.html

IRS Characterization of a Debris Disk around an M-type star in NGC2547

Authors: Paula S. Teixeira, Charles J. Lada, Kenneth Wood, Thomas P. Robitaille, Kevin L. Luhman

Abstract: We present 5 to 15 micron Spitzer Infrared Spectrograph (IRS) low resolution spectral data of a candidate debris disk around an M4.5 star identified as a likely member of the ~40 Myr old cluster NGC2547. The IRS spectrum shows a silicate emission feature, indicating the presence of warm, small, (sub)micron-sized dust grains in the disk. Of the fifteen previously known candidate debris disks arou… ▽ More We present 5 to 15 micron Spitzer Infrared Spectrograph (IRS) low resolution spectral data of a candidate debris disk around an M4.5 star identified as a likely member of the ~40 Myr old cluster NGC2547. The IRS spectrum shows a silicate emission feature, indicating the presence of warm, small, (sub)micron-sized dust grains in the disk. Of the fifteen previously known candidate debris disks around M-type stars, the one we discuss in this paper is the first to have an observed mid-infrared spectrum and is also the first to have measured silicate emission. We combined the IRS data with ancillary data (optical, JHKs, and Spitzer InfraRed Array Camera and 24 micron data) to build the spectral energy distribution (SED) of the source. Monte Carlo radiation transfer modeling of the SED characterized the dust disk as being very flat (h100=2AU) and extending inward within at least 0.13AU of the central star. Our analysis shows that the disk is collisionally dominated and is likely a debris disk. △ Less

Submitted 15 May, 2009; originally announced May 2009.

Comments: Accepted for publication in ApJ. The paper has 19 pages, 4 figures, and one table

Journal ref: Astrophys.J.700:454-459,2009

Dense cores in the Pipe Nebula: An improved core mass function

Authors: J. M. Rathborne, C. J. Lada, A. A. Muench, J. F. Alves, J. Kainulainen, M. Lombardi

Abstract: In this paper we derive an improved core mass function (CMF) for the Pipe Nebula from a detailed comparison between measurements of visual extinction and molecular-line emission. We have compiled a refined sample of 201 dense cores toward the Pipe Nebula using a 2-dimensional threshold identification algorithm informed by recent simulations of dense core populations. Measurements of radial veloc… ▽ More In this paper we derive an improved core mass function (CMF) for the Pipe Nebula from a detailed comparison between measurements of visual extinction and molecular-line emission. We have compiled a refined sample of 201 dense cores toward the Pipe Nebula using a 2-dimensional threshold identification algorithm informed by recent simulations of dense core populations. Measurements of radial velocities using complimentary C18O (1-0) observations enable us to cull out from this sample those 43 extinction peaks that are either not associated with dense gas or are not physically associated with the Pipe Nebula. Moreover, we use the derived C18O, central velocities to differentiate between single cores with internal structure and blends of two or more physically distinct cores, superposed along the same line-of-sight. We then are able to produce a more robust dense core sample for future follow-up studies and a more reliable CMF than was possible previously. We confirm earlier indications that the CMF for the Pipe Nebula departs from a single power-law like form with a break or knee at M ~ 2.7 +/- 1.3 Msun. Moreover, we also confirm that the CMF exhibits a similar shape to the stellar IMF, but is scaled to higher masses by a factor of ~4.5. We interpret this difference in scaling to be a measure of the star formation efficiency (22 +/- 8%). This supports earlier suggestions that the stellar IMF may originate more or less directly from the CMF. △ Less

Submitted 27 April, 2009; originally announced April 2009.

Comments: 30 pages. 5 figures. Accepted by ApJ

Journal ref: Astrophys.J.699:742-753,2009

The Last Gasp of Gas Giant Planet Formation: A Spitzer Study of the 5 Myr-old Cluster NGC 2362

Authors: Thayne Currie, Charles J. Lada, Peter Plavchan, Thomas Robitaille, Jonathan Irwin, Scott J. Kenyon

Abstract: (Abridged) We describe Spitzer IRAC and MIPS observations of the populous, 5 Myr-old open cluster NGC 2362. Early/intermediate-type confirmed/candidate cluster members either have photospheric mid-IR emission or weak, optically-thin infrared excess emission at < 24 microns consistent with debris disks. Few late-type, solar/subsolar-mass stars have primordial disks. The disk population around lat… ▽ More (Abridged) We describe Spitzer IRAC and MIPS observations of the populous, 5 Myr-old open cluster NGC 2362. Early/intermediate-type confirmed/candidate cluster members either have photospheric mid-IR emission or weak, optically-thin infrared excess emission at < 24 microns consistent with debris disks. Few late-type, solar/subsolar-mass stars have primordial disks. The disk population around late-type stars is dominated by disks with inner holes (canonical 'transition disks') and 'homologously depleted' disks. Both types of disks represent an intermediate stage between primordial disks and debris disks. Thus, we find that multiple paths for the primordial-to-debris disk transition exist. Our results undermine standard arguments in favor of a ~ 0.01 Myr year timescale for the transition based on data from Taurus-Auriga and rule out standard UV photoevaporation scenarios as the primary mechanism to explain the transition. Combining our data with other Spitzer surveys, we investigate the evolution of debris disks around high/intermediate-mass stars and investigate timescales for giant planet formation. If the gas and dust in disks evolve on similar timescales, the formation timescale for gas giant planets surrounding early-type, high/intermediate-mass stars is likely 1--5 Myr. Most solar/subsolar-mass stars detected by Spitzer have SEDs that indicate their disks may be actively leaving the primordial disk phase. Thus, gas giant planet formation may also occur by 5 Myr around solar/subsolar-mass stars as well. △ Less

Submitted 23 May, 2009; v1 submitted 15 March, 2009; originally announced March 2009.

Comments: 60 pages, 19 figures, five tables. Published in The Astrophysical Journal (2009, 698, 1)

Journal ref: Astrophys. J.698:1-27, 2009

On the fidelity of the core mass functions derived from dust column density data

Authors: J. Kainulainen, C. J. Lada, J. M. Rathborne, J. F. Alves

Abstract: Aims: We examine the recoverability and completeness limits of the dense core mass functions (CMFs) derived for a molecular cloud using extinction data and a core identification scheme based on two-dimensional thresholding. Methods: We performed simulations where a population of artificial cores was embedded into the variable background extinction field of the Pipe nebula. We extracted the cor… ▽ More Aims: We examine the recoverability and completeness limits of the dense core mass functions (CMFs) derived for a molecular cloud using extinction data and a core identification scheme based on two-dimensional thresholding. Methods: We performed simulations where a population of artificial cores was embedded into the variable background extinction field of the Pipe nebula. We extracted the cores from the simulated extinction maps, constructed the CMFs, and compared them to the input CMFs. The simulations were repeated using a variety of extraction parameters and several core populations with differing input mass functions and differing degrees of crowding. Results: The fidelity of the observed CMF depends on the parameters selected for the core extraction algorithm for our background. More importantly, it depends on how crowded the core population is. We find that the observed CMF recovers the true CMF reliably when the mean separation of cores is larger than their mean diameter (f>1). If this condition holds, the derived CMF is accurate and complete above M > 0.8-1.5 Msun, depending on the parameters used for the core extraction. In the simulations, the best fidelity was achieved with the detection threshold of 1 or 2 times the rms-noise of the extinction data, and with the contour level spacings of 3 times the rms-noise. Choosing larger threshold and wider level spacings increases the limiting mass. The simulations show that when f>1.5, the masses of individual cores are recovered with a typical uncertainty of 25-30 %. When f=1 the uncertainty is ~60 %. In very crowded cases where f<1 the core identification algorithm is unable to recover the masses of the cores adequately. For the cores of the Pipe nebula f~2.0 and therefore the use of the method in that region is justified. △ Less

Submitted 27 January, 2009; originally announced January 2009.

Comments: 9 pages, 6 figures, accepted for publication in A&A

2MASS wide field extinction maps: II. The Ophiuchus and the Lupus cloud complexe

Authors: Marco Lombardi, Charles J. Lada, Joao Alves

Abstract: We present an extinction map of a ~1,700 deg sq region that encloses the Ophiuchus, the Lupus, and the Pipe dark complexes using 42 million stars from the Two Micron All Sky Survey (2MASS) point source catalog. The use of a robust and optimal near-infrared method to map dust column density (Nicer, described in Lombardi & Alves 2001) allow us to detect extinction as low as A_K = 0.05 mag with a 2… ▽ More We present an extinction map of a ~1,700 deg sq region that encloses the Ophiuchus, the Lupus, and the Pipe dark complexes using 42 million stars from the Two Micron All Sky Survey (2MASS) point source catalog. The use of a robust and optimal near-infrared method to map dust column density (Nicer, described in Lombardi & Alves 2001) allow us to detect extinction as low as A_K = 0.05 mag with a 2-sigma significance, and still to have a resolution of 3 arcmin on our map. We also present a novel, statistically sound method to characterize the small-scale inhomogeneities in molecular clouds. Finally, we investigate the cloud structure function, and show that significant deviations from the results predicted by turbulent models are observed. △ Less

Submitted 22 September, 2008; originally announced September 2008.

Comments: 16 pages, A&A in press

The dense ring in the Coalsack: the merging of two subsonic flows

Authors: J. M. Rathborne, C. J. Lada, W. Walsh, M. Saul, H. M. Butner

Abstract: A recent high angular resolution extinction map toward the most opaque molecular globule, Globule 2, in the Coalsack Nebula revealed that it contains a strong central ring of dust column density. This ring represents a region of high density and pressure that is likely a transient and possibly turbulent structure. Dynamical models suggest that the ring has formed as a result of a sudden increase… ▽ More A recent high angular resolution extinction map toward the most opaque molecular globule, Globule 2, in the Coalsack Nebula revealed that it contains a strong central ring of dust column density. This ring represents a region of high density and pressure that is likely a transient and possibly turbulent structure. Dynamical models suggest that the ring has formed as a result of a sudden increase in external pressure which is driving a compression wave into the Globule. Here we combine the extinction measurements with a detailed study of the C18O (1-0) molecular line profiles toward Globule 2 in order to investigate the overall kinematics and, in doing so, test this dynamical model. We find that the ring corresponds to an enhancement in the C18O non-thermal velocity dispersion and non-thermal pressure. We observe a velocity gradient across the Globule that appears to trace two distinct systematic subsonic velocity flows that happen to converge within the ring. We suggest, therefore, that the ring has formed as two subsonic flows of turbulent gas merge within the Globule. The fact that the outer layers of the Globule appear stable against collapse yet there is no centrally condensed core, suggests that the Globule may be evolving from the outside in and has yet to stabilize, confirming its youth. △ Less

Submitted 11 September, 2008; originally announced September 2008.

Comments: 33 pages. 11 figures. Accepted for publication in ApJ

Journal ref: Astrophys.J.690:1659-1671,2009

New M dwarf debris disk candidates in NGC 2547

Authors: Jan Forbrich, Charles J. Lada, August A. Muench, Paula S. Teixeira

Abstract: With only six known examples, M-dwarf debris disks are rare, even though M dwarfs constitute the majority of stars in the Galaxy. After finding a new M dwarf debris disk in a shallow mid-infrared observation of NGC 2547, we present a considerably deeper Spitzer-MIPS image of the region, with a maximum exposure time of 15 minutes per pixel. Among sources selected from a previously published membe… ▽ More With only six known examples, M-dwarf debris disks are rare, even though M dwarfs constitute the majority of stars in the Galaxy. After finding a new M dwarf debris disk in a shallow mid-infrared observation of NGC 2547, we present a considerably deeper Spitzer-MIPS image of the region, with a maximum exposure time of 15 minutes per pixel. Among sources selected from a previously published membership list, we identify nine new M dwarfs with excess emission at 24 micron tracing warm material close to the snow line of these stars, at orbital radii of less than 1 AU. We argue that these are likely debris disks, suggesting that planet formation is under way in these systems. Interestingly, the estimated excess fraction of M stars appears to be higher than that of G and K stars in our sample. △ Less

Submitted 23 July, 2008; originally announced July 2008.

Comments: 16 pages, 8 figures; accepted for publication in ApJ

Discovery of superthermal hydroxyl (OH) in the HH211 outflow

Authors: A. Tappe, C. J. Lada, J. H. Black, A. A. Muench

Abstract: We present a 5-37 micron infrared spectrum obtained with the Spitzer Space Telescope toward the southeastern lobe of the young protostellar outflow HH211. The spectrum shows an extraordinary sequence of OH emission lines arising in highly excited rotational levels up to an energy E/k~28200K above the ground level. This is, to our knowledge, by far the highest rotational excitation of OH observed… ▽ More We present a 5-37 micron infrared spectrum obtained with the Spitzer Space Telescope toward the southeastern lobe of the young protostellar outflow HH211. The spectrum shows an extraordinary sequence of OH emission lines arising in highly excited rotational levels up to an energy E/k~28200K above the ground level. This is, to our knowledge, by far the highest rotational excitation of OH observed outside Earth. The spectrum also contains several pure rotational transitions of H2O (v=0), H2 (v=0) S(0) to S(7), HD (v=0) R(3) to R(6), and atomic fine-structure lines of [Fe II], [Si II], [Ne II], [S I], and [Cl I]. The origin of the highly excited OH emission is most likely the photodissociation of H2O by the UV radiation generated in the terminal outflow shock of HH211. △ Less

Submitted 12 May, 2008; originally announced May 2008.

Comments: ApJ Letters, in press

Embedded Oscillating Starless Cores

Authors: Avery E. Broderick, Ramesh Narayan, Eric Keto, Charles J. Lada

Abstract: In a previous paper we demonstrated that non-radial hydrodynamic oscillations of a thermally-supported (Bonnor-Ebert) sphere embedded in a low-density, high-temperature medium persist for many periods. The predicted column density variations and molecular spectral line profiles are similar to those observed in the Bok globule B68 suggesting that the motions in some starless cores may be oscillat… ▽ More In a previous paper we demonstrated that non-radial hydrodynamic oscillations of a thermally-supported (Bonnor-Ebert) sphere embedded in a low-density, high-temperature medium persist for many periods. The predicted column density variations and molecular spectral line profiles are similar to those observed in the Bok globule B68 suggesting that the motions in some starless cores may be oscillating perturbations on a thermally supported equilibrium structure. Such oscillations can produce molecular line maps which mimic rotation, collapse or expansion, and thus could make determining the dynamical state from such observations alone difficult. However, while B68 is embedded in a very hot, low-density medium, many starless cores are not, having interior/exterior density contrasts closer to unity. In this paper we investigate the oscillation damping rate as a function of the exterior density. For concreteness we use the same interior model employed in Broderick et al. (2007), with varying models for the exterior gas. We also develop a simple analytical formalism, based upon the linear perturbation analysis of the oscillations, which predicts the contribution to the damping rates due to the excitation of sound waves in the external medium. We find that the damping rate of oscillations on globules in dense molecular environments is always many periods, corresponding to hundreds of thousands of years, and persisting over the inferred lifetimes of the globules. △ Less

Submitted 10 April, 2008; originally announced April 2008.

Comments: 8 pages, 3 figures; accepted by ApJ