Search

VISIONS is an ESO public survey of five nearby (d < 500 pc) star-forming molecular cloud complexes that are canonically associated with the constellations of Chamaeleon, Corona Australis, Lupus, Ophiuchus, and Orion. The survey was carried out with VISTA, using VIRCAM, and collected data in the near-infrared passbands J, H, and Ks. With a total on-sky exposure time of 49.4 h VISIONS covers an area of 650 deg$^2$, and it was designed to build an infrared legacy archive similar to that of 2MASS. Taking place between April 2017 and March 2022, the observations yielded approximately 1.15 million images, which comprise 19 TB of raw data. The observations are grouped into three different subsurveys: The wide subsurvey comprises shallow, large-scale observations and has visited the star-forming complexes six times over the course of its execution. The deep subsurvey of dedicated high-sensitivity observations has collected data on the areas with the largest amounts of dust extinction. The control subsurvey includes observations of areas of low-to-negligible dust extinction. Using this strategy, the VISIONS survey offers multi-epoch position measurements, is able to access deeply embedded objects, and provides a baseline for statistical comparisons and sample completeness. In particular, VISIONS is designed to measure the proper motions of point sources with a precision of 1 mas/yr or better, when complemented with data from VHS. Hence, VISIONS can provide proper motions for sources inaccessible to Gaia. VISIONS will enable addressing a range of topics, including the 3D distribution and motion of embedded stars and the nearby interstellar medium, the identification and characterization of young stellar objects, the formation and evolution of embedded stellar clusters and their initial mass function, as well as the characteristics of interstellar dust and the reddening law., Accepted for publication in Astronomy & Astrophysics on 19 January 2023

We present millimeter, optical, and soft X-ray observations of a stellar flare with an energy squarely in the regime of typical X1 solar flares. The flare was observed from Proxima Cen on 2019 May 6 as part of a larger multi-wavelength flare monitoring campaign and was captured by Chandra, LCOGT, du Pont, and ALMA. Millimeter emission appears to be a common occurrence in small stellar flares that had gone undetected until recently, making it difficult to interpret these events within the current multi-wavelength picture of the flaring process. The May 6 event is the smallest stellar millimeter flare detected to date. We compare the relationship between the soft X-ray and millimeter emission to that observed in solar flares. The X-ray and optical flare energies of 10$^{30.3\pm0.2}$ and 10$^{28.9\pm0.1}$ erg, respectively, the coronal temperature of T=11.0$\pm$2.1 MK, and the emission measure of 9.5$\pm$2.2 X 10$^{49}$ cm$^{-3}$ are consistent with M-X class solar flares. We find the soft X-ray and millimeter emission during quiescence are consistent with the Gudel-Benz Relation, but not during the flare. The millimeter luminosity is >100X higher than that of an equivalent X1 solar flare and lasts only seconds instead of minutes as seen for solar flares., 12 pages, 5 figures. Accepted to ApJ

We present the four-year survey results of monthly submillimeter monitoring of eight nearby ($< 500 $pc) star-forming regions by the JCMT Transient Survey. We apply the Lomb-Scargle Periodogram technique to search for and characterize variability on 295 submillimeter peaks brighter than 0.14 Jy beam$^{-1}$, including 22 disk sources (Class II), 83 protostars (Class 0/I), and 190 starless sources. We uncover 18 secular variables, all of them protostars. No single-epoch burst or drop events and no inherently stochastic sources are observed. We classify the secular variables by their timescales into three groups: Periodic, Curved, and Linear. For the Curved and Periodic cases, the detectable fractional amplitude, with respect to mean peak brightness, is $\sim4$ % for sources brighter than $\sim$ 0.5 Jy beam$^{-1}$. Limiting our sample to only these bright sources, the observed variable fraction is 37 % (16 out of 43). Considering source evolution, we find a similar fraction of bright variables for both Class 0 and Class I. Using an empirically motivated conversion from submillimeter variability to variation in mass accretion rate, six sources (7 % of our full sample) are predicted to have years-long accretion events during which the excess mass accreted reaches more than 40 % above the total quiescently accreted mass: two previously known eruptive Class I sources, V1647 Ori and EC 53 (V371 Ser), and four Class 0 sources, HOPS 356, HOPS 373, HOPS 383, and West 40. Considering the full protostellar ensemble, the importance of episodic accretion on few years timescale is negligible, only a few percent of the assembled mass. However, given that this accretion is dominated by events of order the observing time-window, it remains uncertain as to whether the importance of episodic events will continue to rise with decades-long monitoring., Accepted for publication in the Astrophysical Journal

We present a newly enlarged census of the compact radio population towards the Orion Nebula Cluster (ONC) using high-sensitivity continuum maps (3-10 $��$Jy bm$^{-1}$) from a total of $\sim30$ h centimeter-wavelength observations over an area of $\sim$20$'\times20'$ obtained in the C-band (4$-$8 GHz) with the Karl G. Jansky Very Large Array (VLA) in its high-resolution A-configuration. We thus complement our previous deep survey of the innermost areas of the ONC, now covering the field of view of the Chandra Orion Ultra-deep Project (COUP). Our catalog contains 521 compact radio sources of which 198 are new detections. Overall, we find that 17% of the (mostly stellar) COUP sources have radio counterparts, while 53% of the radio sources have COUP counterparts. Most notably, the radio detection fraction of X-ray sources is higher in the inner cluster and almost constant for $r>3'$ (0.36 pc) from $��^1$ Ori C suggesting a correlation between the radio emission mechanism of these sources and their distance from the most massive stars at the center of the cluster, for example due to increased photoionisation of circumstellar disks. The combination with our previous observations four years prior lead to the discovery of fast proper motions of up to $\sim$373 km s$^{-1}$ from faint radio sources associated with ejecta of the OMC1 explosion. Finally, we search for strong radio variability. We found changes in flux density by a factor of $\lesssim$5 within our observations and a few sources with changes by a factor $>$10 on long timescales of a few years., Accepted for publication in MNRAS, 17 pages, 12 figures

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 $\alpha_{12} = 8.7 \pm 3.9$, $\alpha_{13} = 48.9 \pm 20.4$ and $\alpha_{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., Comment: 24 pages, 11 figures, accepted for publication in The Astrophysical Journal

Spectroscopic observations of white dwarfs reveal that many of them are polluted by exoplanetary material, whose bulk composition can be uniquely probed this way. We present a spectroscopic and photometric analysis of the DA white dwarf WDJ181417.84$-$735459.83, an object originally identified to have a strong infrared excess in the 2MASS and WISE catalogues that we confirmed to be intrinsic to the white dwarf, and likely corresponding to the emission of a dusty disc around the star. The finding of Ca, Fe and Mg absorption lines in two X-SHOOTER spectra of the white dwarf, taken 8 years apart, is further evidence of accretion from a dusty disc. We do not report variability in the absorption lines between these two spectra. Fitting a blackbody model to the infrared excess gives a temperature of 910$\pm50$ K. We have estimated a total accretion flux from the spectroscopic metal lines of $|\dot{\rm M}| = 1.784 \times 10^{9}\, $g s$^{-1}$., Accepted for publication in MNRAS. 11 pages, 10 figures, 6 tables

Young and magnetically active low-mass stars often exhibit non-thermal coronal radio emission due to the gyration of electrons in their magnetized chromospheres. This emission is easily detectable at centimeter wavelengths with the current sensitivity of large radio interferometers like the VLA. With the aim of identifying nearby stars adequate for future accurate radio astrometric monitoring using Very Long Baseline Interferometry (VLBI), we have used the VLA in its B\, configuration to search for radio emission at \mbox{$\nu \simeq 6$\,GHz} \mbox{($\lambda \simeq 5$\,cm)} toward a sample of 170 nearby ($5000$\,K. The binarity fraction among the radio-bright stars is at least twice as high as among the radio-quiet stars. The radio-bright nearby young stars identified here provide an interesting sample for future astrometric studies using VLBI arrays aimed at searching for hitherto unknown tight binary components or even exoplanets., Comment: 24 pages, 15 figures plus electronic figures with 31 figures, 3 tables, accepted by ApJ

From Very Long Baseline Array (VLBA) observations we previously identified a population of 123 young stellar systems with nonthermal radio emission toward the core of the Orion Nebula Cluster (ONC). We find optical sources in the Gaia DR2 catalog for 34 of them within 0.2" of the radio positions. Most of the radio sources are likely companions of Gaia detections. However, there are 11 stars whose VLBA position differs from Gaia by $, Comment: Accepted for publication in ApJ. 16 pages, 4 figures, 7 tables

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, $\alpha_{\rm\,CO(2-1)}$, toward individual GMCs across an external galaxy. Our direct measurement yields an average CO--to--dust mass conversion factor of $\alpha^\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 $\alpha_{\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., Comment: ApJ, in press

We present first results of a four-epoch VLBA survey for nonthermal emission toward all 556 compact radio sources previously identified in a deep VLA survey of the Orion Nebula Cluster (ONC). We identify VLBA counterparts toward an unprecedented 123 sources. Of these, 41 do not have X-ray counterparts, of which 34 also do not display near-infrared counterparts. Since these cannot be explained by extragalactic background sources, this suggests a component of the ONC population of young stellar objects that are too deeply embedded for even X-rays to be detectable. We find pervasive variability and detect even most of the highest-S/N sources in only one out of four epochs. Neither a negative spectral index nor extreme variability in the VLA data is a good predictor of a VLBA detection., Comment: ApJ, in press

A search of the first Data Release of the VISTA Variables in the Via Lactea (VVV) Survey discovered the exceptionally red transient VVV-WIT-01 ($H-K_s=5.2$). It peaked before March 2010, then faded by $\sim$9.5 mag over the following two years. The 1.6--22 $\mu$m spectral energy distribution in March 2010 was well fit by a highly obscured black body with $T \sim 1000$ K and $A_{K_s} \sim 6.6$ mag. The source is projected against the Infrared Dark Cloud (IRDC) SDC G331.062$-$0.294. The chance projection probability is small for any single event ($p \approx 0.01$ to 0.02) which suggests a physical association, e.g. a collision between low mass protostars. However, black body emission at $T \sim 1000$ K is common in classical novae (especially CO novae) at the infrared peak in the light curve, due to condensation of dust $\sim$30--60 days after the explosion. Radio follow up with the Australia Telescope Compact Array (ATCA) detected a fading continuum source with properties consistent with a classical nova but probably inconsistent with colliding protostars. Considering all VVV transients that could have been projected against a catalogued IRDC raises the probability of a chance association to $p=0.13$ to 0.24. After weighing several options, it appears likely that VVV-WIT-01 was a classical nova event located behind an IRDC., Comment: 13 pages, 7 figures. Accepted by MNRAS

The connection between low- and high-mass filaments is a matter of strong debate. In order to bridge these two filamentary regimes, we have investigated the internal structure of the Integral Filament in Orion using ALMA observations of the N2H+ (1-0) emission line in Band 3 in combination with previous single-dish data. Our ALMA mosaics, the largest of its kind carried out so-far in local clouds, reveal the presence of multiple sonic-like fibers inside this massive filament. In combination with the identification of fibers in regions such as Taurus, Musca, and Perseus, the first unambiguous detection of fibers in Orion highlights the importance of these gas substructures as the fundamental building blocks of both low- and high-mass filaments.

We have compiled the most complete compact and ultracompact HII region catalogue to date via multi-wavelength inspection of survey data. We utilise data from the recently available SASSy 850$��$m survey to identify massive star forming clumps in the outer Galaxy ($R_{\rm{GC}}>8.5$ kpc) and cross-match with infrared and radio data of known UC HII regions from the RMS database. For the inner Galaxy sample ($R_{\rm{GC}}, updated to correct typos. 15 pages

Observations of young stellar objects (YSOs) in centimeter bands can probe the continuum emission from growing dust grains, ionized winds, and magnetospheric activity, which are intimately connected to the evolution of protoplanetary disks and the formation of planets. We have carried out sensitive continuum observations toward the Ophiuchus A star-forming region using the Karl G. Jansky Very Large Array (VLA) at 10 GHz over a field-of-view of 6$'$ with a spatial resolution of $\theta_{maj}$ $\times$ $\theta_{min}$ $\sim$ 0.4$''$ $\times$ 0.2$''$. We achieved a 5 $\mu$Jy beam$^{-1}$ root-mean-square noise level at the center of our mosaic field of view. Among the eighteen sources we detected, sixteen are YSOs (three Class 0, five Class I, six Class II, and two Class III) and two are extragalactic candidates. We find that thermal dust emission generally contributes less that 30% of the emission at 10 GHz. The radio emission is dominated by other types of emission such as gyro-synchrotron radiation from active magnetospheres, free-free emission from thermal jets, free-free emission from the outflowing photo-evaporated disk material, and/or synchrotron emission from accelerated cosmic-rays in jet or protostellar surface shocks. These different types of emission could not be clearly disentangled. Our non-detections towards Class II/III disks suggest that extreme UV-driven photoevaporation is insufficient to explain the disk dispersal, assuming that the contribution of UV photoevaporating stellar winds to radio flux does not evolve with time. The sensitivity of our data cannot exclude photoevaporation due to X-ray photons as an efficient mechanism for disk dispersal. Deeper surveys with the Square Kilometre Array will be able to provide strong constraints on disk photoevaporation., Comment: Accepted in A&A

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, 39 pages, 25 figures, Accepted for publication by A&A

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., 9 pages, 5 figures, 4 tables, accepted for publication in A&A

Studies both of entire galaxies and of local Galactic star formation indicate a dependency of a molecular cloud's star formation rate (SFR) on its dense gas mass. In external galaxies, such measurements are derived from HCN(1-0) observations, usually encompassing many Giant Molecular Clouds (GMCs) at once. The Andromeda galaxy (M31) is a unique laboratory to study the relation of the SFR and HCN emission down to GMC scales at solar-like metallicities. In this work, we correlate our composite SFR determinations with archival HCN, HCO$^+$, and CO observations, resulting in a sample of nine reasonably representative GMCs. We find that, at the scale of individual clouds, it is important to take into account both obscured and unobscured star formation to determine the SFR. When correlated against the dense-gas mass from HCN, we find that the SFR is low, in spite of these refinements. We nevertheless retrieve an SFR - dense-gas mass correlation, confirming that these SFR tracers are still meaningful on GMC scales. The correlation improves markedly when we consider the HCN/CO ratio instead of HCN by itself. This nominally indicates a dependency of the SFR on the dense-gas fraction, in contradiction to local studies. However, we hypothesize that this partly reflects the limited dynamic range in dense-gas mass, and partly that the ratio of single-pointing HCN and CO measurements may be less prone to systematics like sidelobes. In this case, the HCN/CO ratio would importantly be a better empirical measure of the dense-gas content itself., 9 pages, 6 figures, accepted for publication in MNRAS

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 $\Delta V$ ($\Delta 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., Comment: 40 pages, 26 figures, accepted for publication in ApJ

In previous contributions, we have presented an analytical model describing the evolution of molecular clouds (MCs) undergoing hierarchical gravitational contraction. The cloud's evolution is characterized by an initial increase in its mass, density, and star formation rate (SFR) and efficiency (SFE) as it contracts, followed by a decrease of these quantities as newly formed massive stars begin to disrupt the cloud. The main parameter of the model is the maximum mass reached by the cloud during its evolution. Thus, specifying the instantaneous mass and some other variable completely determines the cloud's evolutionary stage. We apply the model to interpret the observed scatter in SFEs of the cloud sample compiled by Lada et al.\ as an evolutionary effect so that, although clouds such as California and Orion A have similar masses, they are in very different evolutionary stages, causing their very different observed SFRs and SFEs. The model predicts that the California cloud will eventually reach a significantly larger total mass than the Orion A cloud. Next, we apply the model to derive estimated ages of the clouds since the time when approximately 25\% of their mass had become molecular. We find ages from $\sim 1.5$ to 27 Myr, with the most inactive clouds being the youngest. Further predictions of the model are that clouds with very low SFEs should have massive atomic envelopes constituting the majority of their gravitational mass, and that low-mass clouds ($M \sim 10^3$-$10^4 \, M_\odot$) end their lives with a mini-burst of star formation, reaching SFRs $\sim 300$-$500\, M_\odot$ Myr$^{-1}$. By this time, they have contracted to become compact ($\sim 1$ pc) massive star-forming clumps, in general embedded within larger GMCs., Comment: 11 pages, 3 figures, and 2 tables. Published in MNRAS. Journal-ref: MNRAS, 2018, Volume 479, p.3254-3263. ADS link: https://ui.adsabs.harvard.edu/abs/2018MNRAS.479.3254V/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 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., 14 pages, 22 figures. Accepted for publication in Astronomy & Astrophysics

Aims: We aim to use the high spatial resolution of the Atacama Large Millimeter/submillimeter Array (ALMA) to map the flow pattern of molecular gas near DG Tau and its disk, a young stellar object driving a jet and a molecular outflow. Methods: We use observations from ALMA in the J = 2 - 1 transition of 12CO, 13CO, and C18O to study the Keplerian disk of DG Tau and outflows that may be related to the disk and the jet. Results: We find a new wind component flowing radially at a steep angle (~25 deg from the vertical) above the disk with a velocity of ~ 3.1 km/s. It continues the trend of decreasing velocity for increasing distance from the jet axis ("onion-like velocity structure"). Conclusions: The new component is located close to the protostellar disk surface and may be related to photoevaporative winds., Comment: Astronomy and Astrophysics Letters, in press, 8 pages

Most protostars have luminosities that are fainter than expected from steady accretion over the protostellar lifetime. The solution to this problem may lie in episodic mass accretion -- prolonged periods of very low accretion punctuated by short bursts of rapid accretion. However, the timescale and amplitude for variability at the protostellar phase is almost entirely unconstrained. In "A JCMT/SCUBA-2 Transient Survey of Protostars in Nearby Star Forming Regions", we are monitoring monthly with SCUBA-2 the sub-mm emission in eight fields within nearby (, ApJ, accepted

Young stellar objects are known to exhibit strong radio variability on timescales of weeks to months, and a few reports have documented extreme radio flares, with at least an order of magnitude change in flux density on timescales of hours to days. However, there have been few constraints on the occurrence rate of such radio flares or on the correlation with pre-main sequence X-ray flares, although such correlations are known for the Sun and nearby active stars. Here we report simultaneous deep VLA radio and Chandra X-ray observations of the Orion Nebula Cluster, targeting hundreds of sources to look for the occurrence rate of extreme radio variability and potential correlation with the most extreme X-ray variability. We identify 13 radio sources with extreme radio variability, with some showing an order of magnitude change in flux density in less than 30 minutes. All of these sources show X-ray emission and variability, but only on timescales, 18 pages, 13 sources, 7 figures. ApJ, in press

Using multi-epoch VLA observations, covering a time baseline of 29.1 years, we have measured the proper motions of 88 young stars with compact radio emission in the core of the Orion Nebula Cluster (ONC) and the neighboring BN/KL region. Our work increases the number of young stars with measured proper motion at radio frequencies by a factor of 2.5 and enables us to perform a better statistical analysis of the kinematics of the region than was previously possible. Most stars (79 out of 88) have proper motions consistent with a Gaussian distribution centered on $\overline{\mu_{\alpha}\cos{\delta}}=1.07\pm0.09\quad{\rm mas\,yr^{-1}}$, and $\overline{\mu_{\delta}}=-0.84\pm0.16\quad{\rm mas\,yr^{-1}}$, with velocity dispersions of $\sigma_{\alpha}=1.08\pm0.07\quad{\rm mas\,\,yr^{-1}},$ $\sigma_{\delta}=1.27\pm0.15\quad{\rm mas\,\,yr^{-1}}$. We looked for organized movements of these stars but found no clear indication of radial expansion/contraction or rotation. The remaining nine stars in our sample show peculiar proper motions that differ from the mean proper motions of the ONC by more than 3-$\sigma$. One of these stars, V 1326 Ori, could have been expelled from the Orion Trapezium 7,000 years ago. Two could be related to the multi-stellar disintegration in the BN/KL region, in addition to the previously known sources BN, I and n. The others either have high uncertainties (so their anomalous proper motions are not firmly established) or could be foreground objects., Comment: 15 pages, 5 figures, 2 tables, 1 appendix, accepted in The Astrophysical Journal

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 $\mu$m; the final size and mass estimates are based on the observations at 250 $\mu$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$\alpha$ 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., Comment: 18 pages, 19 figures, appendix with 2 figures

Young stellar objects (YSOs) may undergo periods of active accretion (outbursts), during which the protostellar accretion rate is temporarily enhanced by a few orders of magnitude. Whether or not these accretion outburst YSOs possess similar dust/gas reservoirs to each other, and whether or not their dust/gas reservoirs are similar as quiescent YSOs, are issues not yet clarified. The aim of this work is to characterize the millimeter thermal dust emission properties of a statistically significant sample of long and short duration accretion outburst YSOs (i.e., FUors and EXors) and the spectroscopically identified candidates of accretion outbursting YSOs (i.e., FUor-like objects). We have carried out extensive Submillimeter Array (SMA) observations mostly at $\sim$225 GHz (1.33 mm) and $\sim$272 GHz (1.10 mm), from 2008 to 2017. We covered accretion outburst YSOs located at $$3-$\sigma$ significance. Detected sources except for the two cases of V883 Ori and NGC 2071 MM3 were observed with $\sim$1$"$ angular resolution. Overall our observed targets show a systematically higher millimeter luminosity distribution than those of the $M_{*}>$0.3 $M_{\odot}$ Class II YSOs in the nearby ($\lesssim$400 pc) low-mass star-forming molecular clouds (e.g., Taurus, Lupus, Upp Scorpio, and Chameleon I). In addition, at 1 mm our observed confirmed binaries or triple-system sources are systematically fainter than the rest of the sources even though their 1 mm fluxes are broadly distributed. We may have detected $\sim$30-60\% millimeter flux variability from V2494 Cyg and V2495 Cyg, from the observations separated by $\sim$1 year., Comment: 23 pages, 17 figures, accepted to A&A

We present the Coordinated Synoptic Investigation of NGC 2264, a continuous 30-day multi-wavelength photometric monitoring campaign on more than 1000 young cluster members using 16 telescopes. The unprecedented combination of multi-wavelength, high-precision, high-cadence, and long-duration data opens a new window into the time domain behavior of young stellar objects. Here we provide an overview of the observations, focusing on results from Spitzer and CoRoT. The highlight of this work is detailed analysis of 162 classical T Tauri stars for which we can probe optical and mid-infrared flux variations to 1% amplitudes and sub-hour timescales. We present a morphological variability census and then use metrics of periodicity, stochasticity, and symmetry to statistically separate the light curves into seven distinct classes, which we suggest represent different physical processes and geometric effects. We provide distributions of the characteristic timescales and amplitudes, and assess the fractional representation within each class. The largest category (>20%) are optical "dippers" having discrete fading events lasting ~1-5 days. The degree of correlation between the optical and infrared light curves is positive but weak; notably, the independently assigned optical and infrared morphology classes tend to be different for the same object. Assessment of flux variation behavior with respect to (circum)stellar properties reveals correlations of variability parameters with H$\alpha$ emission and with effective temperature. Overall, our results point to multiple origins of young star variability, including circumstellar obscuration events, hot spots on the star and/or disk, accretion bursts, and rapid structural changes in the inner disk., Comment: Published in AJ. 59 pages; 4 tables; 49 figures, most of which are highly degraded to fit size limits. Author name typo corrected. For a better resolution version, please visit http://web.ipac.caltech.edu/staff/amc/codyetal2014.pdf

We present multi-epoch astrometric radio observations with the Very Long Baseline Array (VLBA) of the young ultracool-dwarf binary LSPM J1314+1320AB . The radio emission comes from the secondary star. Combining the VLBA data with Keck near-infrared adaptive-optics observations of both components, a full astrometric fit of parallax ($\pi_{\rm abs}=57.975\pm0.045$ mas, corresponding to a distance of $d=17.249\pm0.013$ pc), proper motion ($\mu_{\rm \alpha cos \delta}=-247.99\pm0.10$ mas yr$^{-1}$, $\mu_{\delta}=-183.58\pm0.22$ mas yr$^{-1}$), and orbital motion is obtained. Despite the fact that the two components have nearly identical masses to within $\pm2$%, the secondary's radio emission exceeds that of the primary by a factor of $\gtrsim$30, suggesting a difference in stellar rotation history, which could result in different magnetic field configurations. Alternatively, the emission could be anisotropic and beamed toward us for the secondary but not for the primary. Using only reflex motion, we exclude planets of mass 0.7 to 10 $M_{\rm jup}$ with orbital periods of 600 to 10 days, respectively. Additionally, we use the full orbital solution of the binary to derive an upper limit for the semi-major axis of 0.23 AU for stable planetary orbits within this system. These limits cover a parameter space that is inaccessible with, and complementary to, near-infrared radial velocity surveys of ultracool dwarfs. Our absolute astrometry will constitute an important test for the astrometric calibration of Gaia., Comment: ApJ, in press

We present the first high angular resolution study of giant molecular clouds (GMCs) in the nearby spiral galaxy NGC 300, based on observations from the Submillimeter Array (SMA). We target eleven 500 pc-sized regions of active star formation within the galaxy in the CO(J=2-1) line at 40 pc spatial and 1 km/s spectral resolution and identify 45 individual GMCs. We characterize the physical properties of these GMCs, and find that they are similar to GMCs in the disks of the Milky Way and other nearby spiral galaxies. For example, the GMC mass spectrum in our sample has a slope of 1.80+/-0.07. Twelve clouds are spatially resolved by our observations, of which ten have virial mass estimates that agree to within a factor of two with mass estimates derived directly from CO integrated intensity, suggesting that the majority of these GMCs are bound. The resolved clouds show consistency with Larson's fundamental relations between size, linewidth, and mass observed in the Milky Way. We find that the linewidth scales with the size as DeltaV ~ R^(0.52+/-0.20), and the median surface density in the subsample is 54 Msun/pc^(-2). We detect 13CO in four GMCs and find a mean 12CO/13CO flux ratio of 6.2. Our interferometric observations recover between 30% and 100% of the integrated intensity from the APEX single dish CO observations of Faesi et al. 2014, suggesting the presence of low-mass GMCs and/or diffuse gas below our sensitivity limit. The fraction of APEX emission recovered increases with the SMA total intensity as well as with the star formation rate., 16 pages, 10 figures, accepted to ApJ

We present a deep centimeter-wavelength catalog of the Orion Nebula Cluster (ONC), based on a 30h single-pointing observation with the Karl G. Jansky Very Large Array in its high-resolution A-configuration using two 1 GHz bands centered at 4.7 GHz and 7.3 GHz. A total of 556 compact sources were detected in a map with a nominal rms noise of 3 muJy/bm, limited by complex source structure and the primary beam response. Compared to previous catalogs, our detections increase the sample of known compact radio sources in the ONC by more than a factor of seven. The new data show complex emission on a wide range of spatial scales. Following a preliminary correction for the wideband primary-beam response, we determine radio spectral indices for 170 sources whose index uncertainties are less than +/-0.5. We compare the radio to the X-ray and near-infrared point-source populations, noting similarities and differences., Comment: ApJ, in press

Orion A hosts the nearest massive star factory, thus offering a unique opportunity to resolve the processes connected with the formation of both low- and high-mass stars. Here we present the most detailed and sensitive near-infrared (NIR) observations of the entire molecular cloud to date. With the unique combination of high image quality, survey coverage, and sensitivity, our NIR survey of Orion A aims at establishing a solid empirical foundation for further studies of this important cloud. In this first paper we present the observations, data reduction, and source catalog generation. To demonstrate the data quality, we present a first application of our catalog to estimate the number of stars currently forming inside Orion A and to verify the existence of a more evolved young foreground population. We used the European Southern Observatory's (ESO) Visible and Infrared Survey Telescope for Astronomy (VISTA) to survey the entire Orion A molecular cloud in the NIR $J, H$, and $K_S$ bands, covering a total of $\sim$18.3 deg$^2$. We implemented all data reduction recipes independently of the ESO pipeline. Estimates of the young populations toward Orion A are derived via the $K_S$-band luminosity function. Our catalog (799995 sources) increases the source counts compared to the Two Micron All Sky Survey by about an order of magnitude. The 90% completeness limits are 20.4, 19.9, and 19.0 mag in $J, H$, and $K_S$, respectively. The reduced images have 20% better resolution on average compared to pipeline products. We find between 2300 and 3000 embedded objects in Orion A and confirm that there is an extended foreground population above the Galactic field, in agreement with previous work. The Orion A VISTA catalog represents the most detailed NIR view of the nearest massive star-forming region and provides a fundamental basis for future studies of star formation processes toward Orion., Comment: Accepted for publication in A&A, high-quality version available at http://homepage.univie.ac.at/stefan.meingast/vision.pdf

As part of the Young Stellar Object VARiability (YSOVAR) program, we monitored NGC 1333 for ~35 days at 3.6 and 4.5 um using the Spitzer Space Telescope. We report here on the mid-infrared variability of the point sources in the ~10x~20arcmin area centered on 03:29:06, +31:19:30 (J2000). Out of 701 light curves in either channel, we find 78 variables over the YSOVAR campaign. About half of the members are variable. The variable fraction for the most embedded SEDs (Class I, flat) is higher than that for less embedded SEDs (Class II), which is in turn higher than the star-like SEDs (Class III). A few objects have amplitudes (10-90th percentile brightness) in [3.6] or [4.5]>0.2 mag; a more typical amplitude is 0.1-0.15 mag. The largest color change is >0.2 mag. There are 24 periodic objects, with 40% of them being flat SED class. This may mean that the periodic signal is primarily from the disk, not the photosphere, in those cases. We find 9 variables likely to be 'dippers', where texture in the disk occults the central star, and 11 likely to be 'bursters', where accretion instabilities create brightness bursts. There are 39 objects that have significant trends in [3.6]-[4.5] color over the campaign, about evenly divided between redder-when-fainter (consistent with extinction variations) and bluer-when-fainter. About a third of the 17 Class 0 and/or jet-driving sources from the literature are variable over the YSOVAR campaign, and a larger fraction (~half) are variable between the YSOVAR campaign and the cryogenic-era Spitzer observations (6-7 years), perhaps because it takes time for the envelope to respond to changes in the central source. The NGC 1333 brown dwarfs do not stand out from the stellar light curves in any way except there is a much larger fraction of periodic objects (~60% of variable brown dwarfs are periodic, compared to ~30% of the variables overall)., 84 pages, 43 figures. Accepted by AJ

We present a time-variability study of young stellar objects in the cluster IRAS 20050+2720, performed at 3.6 and 4.5 micron with the Spitzer Space Telescope; this study is part of the Young Stellar Object VARiability project (YSOVAR). We have collected light curves for 181 cluster members over 40 days. We find a high variability fraction among embedded cluster members of ca. 70%, whereas young stars without a detectable disk display variability less often (in ca. 50% of the cases) and with lower amplitudes. We detect periodic variability for 33 sources with periods primarily in the range of 2-6 days. Practically all embedded periodic sources display additional variability on top of their periodicity. Furthermore, we analyze the slopes of the tracks that our sources span in the color-magnitude diagram (CMD). We find that sources with long variability time scales tend to display CMD slopes that are at least partially influenced by accretion processes, while sources with short variability time scales tend to display extinction-dominated slopes. We find a tentative trend of X-ray detected cluster members to vary on longer time scales than the X-ray undetected members., Accepted for publication in The Astronomical Journal. Data from Table 2 is given as a machine-readable file in the source of this arXiv submission

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 $$ $=$ 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., Comment: 21 pages, 17 figures, accepted by ApJ

We used the Very Large Array to carry out a multi-epoch radio continuum monitoring of the Orion Nebula Cluster and Orion Molecular Cloud. Our observations reveal the presence of 19 sources. With the exception of the sources BN and C the sources show variability between the different epochs. We have found tentative evidence of variability in the massive object related with source I. Our observations also confirm radio flux density variations of a factor >2 on timescales of hours to days in 5 sources. One of these flaring sources, OHC-E, has been detected for the first time. We conclude that the radio emission arises from: i) highly-variable non-thermal gyrosynchrotron emission produced by electrons accelerated in the magnetospheres of pre-main sequence stars; ii) thermal emission from ionized gas and/or heated dust around massive objects and proplyds. Combining our sample with other radio monitoring and a X-ray catalog, we have studied the properties of 51 radio/X-ray stars. We have found severals hints of a direct relation between the X-ray activity and the mechanisms responsible for (at least some fraction of) the radio emission. We have estimated a radio flaring rate of 0.14 flares day-1 in the densest stellar cluster of the region, suggesting that radio flares from young stars are more common events than previously thought., 20 pages, 18 figures. Accepted in The Astrophysical Journal

This paper continues our study of the foreground population to the Orion molecular clouds. The goal is to characterize the foreground population north of NGC 1981 and to investigate the star formation history in the large Orion star-forming region. We focus on a region covering about 25 square degrees, centered on the $\epsilon$ Orionis supergiant (HD 37128, B0\,Ia) and covering the Orion Belt asterism. We used a combination of optical (SDSS) and near-infrared (2MASS) data, informed by X-ray (\textit{XMM-Newton}) and mid-infrared (WISE) data, to construct a suite of color-color and color-magnitude diagrams for all available sources. We then applied a new statistical multiband technique to isolate a previously unknown stellar population in this region. We identify a rich and well-defined stellar population in the surveyed region that has about 2\,000 objects that are mostly M stars. We infer the age for this new population to be at least 5\, Myr and likely $\sim10$\,Myr and estimate a total of about 2\,500 members, assuming a normal IMF. This new population, which we call the Orion Belt population, is essentially extinction-free, disk-free, and its spatial distribution is roughly centered near $\epsilon$ Ori, although substructure is clearly present. The Orion Belt population is likely the low-mass counterpart to the Ori OB Ib subgroup. Although our results do not rule out Blaauw's sequential star formation scenario for Orion, we argue that the recently proposed blue streams scenario provides a better framework on which one can explain the Orion star formation region as a whole. We speculate that the Orion Belt population could represent the evolved counterpart of an Orion nebula-like cluster., Comment: 13 pages, 10 figures. Accepted for publication to A&A

The YSOVAR (Young Stellar Object VARiability) Spitzer Space Telescope observing program obtained the first extensive mid-infrared (3.6 & 4.5 um) time-series photometry of the Orion Nebula Cluster plus smaller footprints in eleven other star-forming cores (AFGL490, NGC1333, MonR2, GGD 12-15, NGC2264, L1688, Serpens Main, Serpens South, IRAS 20050+2720, IC1396A, and Ceph C). There are ~29,000 unique objects with light curves in either or both IRAC channels in the YSOVAR data set. We present the data collection and reduction for the Spitzer and ancillary data, and define the "standard sample" on which we calculate statistics, consisting of fast cadence data, with epochs about twice per day for ~40d. We also define a "standard sample of members", consisting of all the IR-selected members and X-ray selected members. We characterize the standard sample in terms of other properties, such as spectral energy distribution shape. We use three mechanisms to identify variables in the fast cadence data--the Stetson index, a chi^2 fit to a flat light curve, and significant periodicity. We also identified variables on the longest timescales possible of ~6 years, by comparing measurements taken early in the Spitzer mission with the mean from our YSOVAR campaign. The fraction of members in each cluster that are variable on these longest timescales is a function of the ratio of Class I/total members in each cluster, such that clusters with a higher fraction of Class I objects also have a higher fraction of long-term variables. For objects with a YSOVAR-determined period and a [3.6]-[8] color, we find that a star with a longer period is more likely than those with shorter periods to have an IR excess. We do not find any evidence for variability that causes [3.6]-[4.5] excesses to appear or vanish within our data; out of members and field objects combined, at most 0.02% may have transient IR excesses., Accepted to AJ; 38 figures, 93 pages

We present novel tests of pre$-$main-sequence models based on individual dynamical masses for the M7 binary LSPM J1314+1320AB. Joint analysis of our Keck adaptive optics astrometric monitoring along with Very Long Baseline Array radio data from a companion paper yield component masses of $0.0885\pm0.0006$ $M_{\odot}$ and $0.0875\pm0.0010$ $M_{\odot}$ and a parallactic distance of $17.249\pm0.013$ pc. We also derive component luminosities that are consistent with the system being coeval at an age of $80.8\pm2.5$ Myr, according to BHAC15 evolutionary models. The presence of lithium is consistent with model predictions, marking the first time the theoretical lithium depletion boundary has been tested with ultracool dwarfs of known mass. However, we find that the average evolutionary model-derived effective temperature ($2950\pm5$ K) is 180 K hotter than we derive from a spectral type$-$$T_{\rm eff}$ relation based on BT-Settl models ($2770\pm100$ K). We suggest that the dominant source of this discrepancy is model radii being too small by $\approx$13%. In a test that mimics the typical application of evolutionary models by observers, we derive masses on the H-R diagram using the luminosity and BT-Settl temperature. The estimated masses are $46^{+16}_{-19}$% (2.0$\sigma$) lower than we measure dynamically and would imply that this is a system of $\approx$50 $M_{\rm Jup}$ brown dwarfs, highlighting the large systematic errors possible when inferring masses from the H-R diagram. This is first time masses have been measured for ultracool ($\geq$M6) dwarfs displaying spectral signatures of low gravity. Based on features in the infrared, LSPM J1314+1320AB appears higher gravity than typical Pleiades and AB Dor members, opposite the expectation given its younger age. The components of LSPM J1314+1320AB are now the nearest, lowest mass pre$-$main-sequence stars with direct mass measurements., Comment: ApJ, in press

We compare APOGEE radial velocities (RVs) of young stars in the Orion A cloud with CO line gas emission and find a correlation between the two at large-scales, in agreement with previous studies. However, at smaller scales we find evidence for the presence of substructure in the stellar velocity field. Using a Friends-of-Friends approach we identify 37 stellar groups with almost identical RVs. These groups are not randomly distributed but form elongated chains or strings of stars with five or more members with low velocity dispersion, across lengths of 1-1.5~pc. The similarity between the kinematic properties of the APOGEE strings and the internal velocity field of the chains of dense cores and fibers recently identified in the dense ISM is striking and suggests that for most of the Orion A cloud, young stars keep memory of the parental gas substructure where they originated., 9 pages, 7 figures, Accepted for publication in A&A

LRLL 54361 is an infrared source located in the star forming region IC 348 SW. Remarkably, its infrared luminosity increases by a factor of 10 during roughly one week every 25.34 days. To understand the origin of these remarkable periodic variations, we obtained sensitive 3.3 cm JVLA radio continuum observations of LRLL 54361 and its surroundings in six different epochs: three of them during the IR-on state and three during the IR-off state. The radio source associated with LRLL 54361 remained steady and did not show a correlation with the IR variations. We suggest that the IR is tracing the results of fast (with a timescale of days) pulsed accretion from an unseen binary companion, while the radio traces an ionized outflow with an extent of $\sim$100 AU that smooths out the variability over a period of order a year. The average flux density measured in these 2014 observations, 27$\pm$5 $\mu$Jy, is about a factor of two less than that measured about 1.5 years before, $53\pm$11 $\mu$Jy, suggesting that variability in the radio is present, but over larger timescales than in the IR. We discuss other sources in the field, in particular two infrared/X-ray stars that show rapidly varying gyrosynchrotron emission., Comment: 12 pages, 3 figures, to appear in The Astrophysical Journal

(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, Comment: A&A, in press