Source author record

Jeong-Eun Lee

Jeong-Eun Lee appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

79works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

79 published item(s)

preprint2025arXiv

Early Planet Formation in Embedded Disks (eDisk). XX. Constraining the Chemical Tracers of Young Protostellar Sources

Recent studies indicate that the formation of planets in protoplanetary disks begins early in the embedded Class 0/I phases of protostellar evolution. The physical and chemical makeup of the embedded phase can provide valuable insights into the process of star and planet formation. This study aims to provide a thorough overview of the various morphologies for molecular emissions observed on disk scales toward nearby embedded sources. We present high angular resolution (0.1", 15 au) molecular line emissions for $^{12}$CO, $^{13}$CO, C$^{18}$O, SO, SiO, DCN, CH$_3$OH, H$_2$CO, and c-C$_3$H$_2$ towards 19 nearby protostellar sources in the context of the Atacama Large Millimeter/submillimeter Array (ALMA) Large Program "Early Planet Formation in Embedded Disks (eDisk)". Emissions in $^{12}$CO are seen towards all sources and primarily trace outflowing materials. A few sources also show high-velocity jets in SiO emission and high-velocity channel maps of $^{12}$CO. The $^{13}$CO and C$^{18}$O emissions are well-known tracers of high-density regions and trace the inner envelope and disk regions with clear signs of rotation seen at continuum scales. The large-scale emissions of $^{13}$CO also delineate the outflow cavity walls where the outflowing and infalling materials interact with each other, and exposure to UV radiation leads to the formation of hydrocarbons such as c-C$_3$H$_2$. Both DCN and CH$_3$OH, when detected, show compact emissions from the inner envelope and disk regions that peak at the position of the protostar. The CH$_3$OH emissions are contained within the region of DCN emissions, which suggests that CH$_3$OH traces the hot core regions. Likewise, a few sources also display emissions in CH$_3$OH towards the outflow. Both SO and H$_2$CO show complex morphology among the sources, suggesting that they are formed through multiple processes in protostellar systems.

preprint2024arXiv

The JCMT Transient Survey: Six-Year Summary of 450/850\,$μ$m Protostellar Variability and Calibration Pipeline Version 2.0

The JCMT Transient Survey has been monitoring eight Gould Belt low-mass star-forming regions since December 2015 and six somewhat more distant intermediate-mass star-forming regions since February 2020 with SCUBA-2 on the JCMT at \ShortS and \LongS and with an approximately monthly cadence. We introduce our Pipeline v2 relative calibration procedures for image alignment and flux calibration across epochs, improving on our previous Pipeline v1 by decreasing measurement uncertainties and providing additional robustness. These new techniques work at both \LongS and \ShortNS, where v1 only allowed investigation of the \LongS data. Pipeline v2 achieves better than $0.5^{\prime\prime}$ relative image alignment, less than a tenth of the submillimeter beam widths. The v2 relative flux calibration is found to be 1\% at \LongS and $<5$\% at \ShortNS. The improvement in the calibration is demonstrated by comparing the two pipelines over the first four years of the survey and recovering additional robust variables with v2. Using the full six years of the Gould Belt survey the number of robust variables increases by 50\,\%, and at \ShortS we identify four robust variables, all of which are also robust at \LongNS. The multi-wavelength light curves for these sources are investigated and found to be consistent with the variability being due to dust heating within the envelope in response to accretion luminosity changes from the central source.

preprint2022arXiv

A Gaia view on the star formation in the Monoceros OB1 and R1 associations

Stellar kinematics provides the key to understanding star formation process. In this respect, we present a kinematic study of the Monoceros OB1 (Mon OB1) and R1 (Mon R1) associations using the recent Gaia data and radial velocities of stars derived from high-resolution spectroscopy and the literature. A total of 728 members are selected using the criteria based on the intrinsic properties of young stars, parallaxes, and proper motions. The spatial distribution and kinematic properties of members show that these associations have distinct substructures. In Mon OB1, we find one northern group and two southern groups. Mon R1 is composed of three small stellar groups that are spatially and kinematically distinct. Some stars are found in a halo around these two associations. We detect patterns of expansion for most stellar groups in the associations. In addition, two stellar groups in Mon OB1 show the signature of rotation, which provides an important constraint on cluster formation. The star formation history of Mon OB1 is slightly revised. Star formation first occurred in the southern region and subsequently in the northern region. Recent star-forming events ignited deeper into the southern region, while some stars are escaping from Mon OB1, forming a halo. Mon R1 might have formed at the same epoch as the formation of the northern group in Mon OB1. Given that star formation is taking place on different scales along a large arc-like structure, Mon OB1 and Mon R1 may be the results of hierarchical star formation.

preprint2022arXiv

A VLA View of the Flared, Asymmetric Disk Around the Class 0 Protostar L1527 IRS

We present high resolution Karl G. Jansky Very Large Array (VLA) observations of the protostar L1527 IRS at 7 mm, 1.3 cm, and 2 cm wavelengths. We detect the edge-on dust disk at all three wavelengths and find that it is asymmetric, with the southern side of the disk brighter than the northern side. We confirm this asymmetry through analytic modeling and also find that the disk is flared at 7 mm. We test the data against models including gap features in the intensity profile, and though we cannot rule such models out, they do not provide a statistically significant improvement in the quality of fit to the data. From these fits, we can however place constraints on allowed properties of any gaps that could be present in the true, underlying intensity profile. The physical nature of the asymmetry is difficult to associate with physical features due to the edge-on nature of the disk, but could be related to spiral arms or asymmetries seen in other imaging of more face-on disks.

preprint2022arXiv

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): A Hot Corino Survey toward Protostellar Cores in the Orion Cloud

The presence of complex organic molecules (COMs) in the interstellar medium (ISM) is of great interest since it may link to the origin and prevalence of life in the universe. Aiming to investigate the occurrence of COMs and their possible origins, we conducted a chemical census toward a sample of protostellar cores as part of the ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP) project. We report the detection of 11 hot corino sources, which exhibit compact emissions from warm and abundant COMs, among 56 Class 0/I protostellar cores. All the hot corino sources discovered are likely Class 0 and their sizes of the warm region ($>$ 100 K) are comparable to 100 au. The luminosity of the hot corino sources exhibits positive correlations with the total number of methanol and the extent of its emissions. Such correlations are consistent with the thermal desorption picture for the presence of hot corino and suggest that the lower luminosity (Class 0) sources likely have a smaller region with COMs emissions. With the same sample selection method and detection criteria being applied, the detection rates of the warm methanol in the Orion cloud (15/37) and the Perseus cloud (28/50) are statistically similar when the cloud distances and the limited sample size are considered. Observing the same set of COM transitions will bring a more informative comparison between the cloud properties.

preprint2022arXiv

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): Deriving Inclination Angle and Velocity of the Protostellar Jets from their SiO Knots

We have selected six sources (G209.55-19.68S2, G205.46-14.56S1$_{-}$A, G203.21-11.20W2, G191.90-11.21S, G205.46-14.56S3, and G206.93-16.61W2) from the Atacama Large Millimeter/submillimeter Array Survey of Orion Planck Galactic Cold Clumps (ALMASOP), in which these sources have been mapped in the CO (J=2-1), SiO (J=5-4), and C$^{18}$O (J=2-1) lines. These sources have high-velocity SiO jets surrounded by low-velocity CO outflows. The SiO jets consist of a chain of knots. These knots have been thought to be produced by semi-periodical variations in jet velocity. Therefore, we adopt a shock-forming model, which uses such variations to estimate the inclination angle and velocity of the jets. We also derive the inclination angle of the CO outflows using the wide-angle wind-driven shell model, and find it to be broadly consistent with that of the associated SiO jets. In addition, we apply this shock-forming model to another three protostellar sources with SiO jets in the literature -- HH 211, HH 212, and L1448C(N) -- and find that their inclination angle and jet velocity are consistent with those previously estimated from proper motion and radial velocity studies.

preprint2022arXiv

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): Evidence for a Molecular Jet Launched at an Unprecedented Early Phase of Protostellar evolution

Protostellar outflows and jets play a vital role in star formation as they carry away excess angular momentum from the inner disk surface, allowing the material to be transferred toward the central protostar. Theoretically, low velocity and poorly collimated outflows appear from the beginning of the collapse, at the first hydrostatic core (FHSC) stage. With growing protostellar core mass, high-density jets are launched which entrain an outflow from the infalling envelope. Until now, molecular jets have been observed at high velocity ($\gtrsim$ 100 km/s) in early Class\,0 protostars. We, for the first time, detect a dense molecular jet in SiO emission with small-velocity ($\sim$ 4.2 km\,s$^{-1}$, deprojected $\sim$ 24 km\,s$^{-1}$) from source G208.89-20.04Walma (hereafter, G208Walma) using ALMA Band\,6 observations. This object has some characteristics of FHSCs, such as a small outflow/jet velocity, extended 1.3\,mm continuum emission, and N$_2$D$^+$ line emission. Additional characteristics, however, are typical of early protostars: collimated outflow and SiO jet. The full extent of the outflow corresponds to a dynamical time scale of $\sim$ 930$^{+200}_{-100}$ years. The spectral energy distribution also suggests a very young source having an upper limit of T$_{bol}$ $\sim$ 31 K and L$_{bol}$ $\sim$ 0.8 L$_\sun$. We conclude that G208Walma is likely in the transition phase from FHSC to protostar, and the molecular jet has been launched within a few hundred years of initial collapse. Therefore, G208Walma may be the earliest object discovered in the protostellar phase with a molecular jet.

preprint2022arXiv

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- VIII. A search for hot cores by using C$_2$H$_5$CN, CH$_3$OCHO and CH$_3$OH lines

Hot cores characterized by rich lines of complex organic molecules are considered as ideal sites for investigating the physical and chemical environments of massive star formation. We present a search for hot cores by using typical nitrogen- and oxygen-bearing complex organic molecules (C$_2$H$_5$CN, CH$_3$OCHO and CH$_3$OH), based on ALMA Three-millimeter Observations of Massive Star-forming regions (ATOMS). The angular resolutions and line sensitivities of the ALMA observations are better than 2 arcsec and 10 mJy/beam, respectively. A total of 60 hot cores are identified with 45 being newly detected, in which the complex organic molecules have high gas temperatures ($>$ 100 K) and small source sizes ($<$ 0.1 pc). So far this is the largest sample of hot cores observed with similar angular resolution and spectral coverage. The observations have also shown nitrogen and oxygen differentiation in both line emission and gas distribution in 29 hot cores. Column densities of CH$_3$OH and CH$_3$OCHO increase as rotation temperatures rise. The column density of CH$_3$OCHO correlates tightly with that of CH$_3$OH. The pathways for production of different species are discussed. Based on the spatial position difference between hot cores and UC~H{\sc ii} regions, we conclude that 24 hot cores are externally heated while the other hot cores are internally heated. The observations presented here will potentially help establish a hot core template for studying massive star formation and astrochemistry.

preprint2022arXiv

B-fields in Star-Forming Region Observations (BISTRO): Magnetic Fields in the Filamentary Structures of Serpens Main

We present 850 $μ$m polarimetric observations toward the Serpens Main molecular cloud obtained using the POL-2 polarimeter on the James Clerk Maxwell Telescope (JCMT) as part of the B-fields In STar-forming Region Observations (BISTRO) survey. These observations probe the magnetic field morphology of the Serpens Main molecular cloud on about 6000 au scales, which consists of cores and six filaments with different physical properties such as density and star formation activity. Using the histogram of relative orientation (HRO) technique, we find that magnetic fields are parallel to filaments in less dense filamentary structures where $N_{H_2} < 0.93\times 10^{22}$ cm$^{-2}$ (magnetic fields perpendicular to density gradients), while being perpendicular to filaments (magnetic fields parallel to density gradients) in dense filamentary structures with star formation activity. Moreover, applying the HRO technique to denser core regions, we find that magnetic field orientations change to become perpendicular to density gradients again at $N_{H_2} \approx 4.6 \times 10^{22}$ cm$^{-2}$. This can be interpreted as a signature of core formation. At $N_{H_2} \approx 16 \times 10^{22}$ cm$^{-2}$ magnetic fields change back to being parallel to density gradients once again, which can be understood to be due to magnetic fields being dragged in by infalling material. In addition, we estimate the magnetic field strengths of the filaments ($B_{POS} = 60-300~μ$G)) using the Davis-Chandrasekhar-Fermi method and discuss whether the filaments are gravitationally unstable based on magnetic field and turbulence energy densities.

preprint2022arXiv

Dissecting the different components of the modest accretion bursts of the very young protostar HOPS 373

Observed changes in protostellar brightness can be complicated to interpret. In our JCMT~Transient monitoring survey, we discovered that a young binary protostar, HOPS 373, is undergoing a modest $30\%$ brightness increase at 850 $μ$m, caused by a factor of 1.8$-$3.3 enhancement in the accretion rate. The initial burst occurred over a few months, with a sharp rise and then shallower decay. A second rise occurred soon after the decay, and the source is still bright one year later. The mid-IR emission, the small-scale CO outflow mapped with ALMA, and the location of variable maser emission indicate that the variability is associated with the SW component. The near-infrared and NEOWISE $W1$ and $W2$ emission is located along the blueshifted CO outflow, spatially offset by $\sim3$ to $4^{\prime\prime}$ from the SW component. The $K$-band emission imaged by UKIRT shows a compact H$_2$ emission source at the edge of the outflow, with a tail tracing the outflow back to the source. The $W1$ emission, likely dominated by scattered light, brightens by 0.7 mag, consistent with expectations based on the sub-mm lightcurve. The signal of continuum variability in $K$-band and $W2$ is masked by stable H$_2$ emission, as seen in our Gemini/GNIRS spectrum, and perhaps by CO emission. These differences in emission sources complicate infrared searches for variability of the youngest protostars.

preprint2022arXiv

Ice features of low-luminosity protostars in near-infrared spectra of AKARI/IRC

We present near-infrared spectra of three low-luminosity protostars and one background star in the Perseus molecular cloud, acquired using the Infrared Camera (IRC) onboard the \textit{AKARI} space telescope. For the comparison with different star-forming environments, we also present spectra of the massive protostar AFGL 7009S, where the protostellar envelope is heated significantly, and the low-mass protostar RNO 91, which is suspected to be undergoing an episodic burst. We detected ice absorption features of \ch{H2O}, \ch{CO2}, and \ch{CO} at all spectra around the wavelengths of 3.05, 4.27, and 4.67 $μ$m, respectively. At least two low-luminosity protostars, we also detected the \ch{XCN} ice feature at 4.62 $μ$m. The presence of the crystalline \ch{H2O} ice and \ch{XCN} ice components indicates that the low-luminosity protostars experienced a hot phase via accretion bursts during the past mass accretion process. We compared the ice abundances of low-luminosity protostars with those of the embedded low-mass protostars and the dense molecular clouds and cores, suggesting that their ice abundances reflect the strength of prior bursts and the timescale after the last burst.

preprint2022arXiv

Nobeyama Survey of Inward Motions toward Cores in Orion Identified by SCUBA-2

In this study, 36 cores (30 starless and 6 protostellar) identified in Orion were surveyed to search for inward motions. We used the Nobeyama 45 m radio telescope, and mapped the cores in the $J = 1\rightarrow0$ transitions of HCO$^+$, H$^{13}$CO$^+$, N$_2$H$^+$, HNC, and HN$^{13}$C. The asymmetry parameter $δV$, which was the ratio of the difference between the HCO$^+$ and H$^{13}$CO$^+$ peak velocities to the H$^{13}$CO$^+$ line width, was biased toward negative values, suggesting that inward motions were more dominant than outward motions. Three starless cores (10% of all starless cores surveyed) were identified as cores with blue-skewed line profiles (asymmetric profiles with more intense blue-shifted emission), and another two starless cores (7%) were identified as candidate blue-skewed line profiles. The peak velocity difference between HCO$^+$ and H$^{13}$CO$^+$ of them was up to 0.9 km s$^{-1}$, suggesting that some inward motions exceeded the speed of sound for the quiescent gas ($\sim10-17$ K). The mean of $δV$ of the five aforementioned starless cores was derived to be $-$0.5$\pm$0.3. One core, G211.16$-$19.33North3, observed using the ALMA ACA in DCO$^+$ $J = 3\rightarrow2$ exhibited blue-skewed features. Velocity offset in the blue-skewed line profile with a dip in the DCO$^+$ $J = 3\rightarrow2$ line was larger ($\sim 0.5$ km s$^{-1}$) than that in HCO$^+$ $J = 1\rightarrow0$ ($\sim 0.2$ km s$^{-1}$), which may represent gravitational acceleration of inward motions. It seems that this core is at the last stage in the starless phase, judging from the chemical evolution factor version 2.0 (CEF2.0).

preprint2021arXiv

Atomic Shocks in the Outflow of L1551 IRS 5 Identified with SOFIA-upGREAT Observations of [OI]

We present velocity resolved SOFIA/upGREAT observations of [OI] and [CII] lines toward a Class I protostar, L1551 IRS 5, and its outflows. The SOFIA observations detect [OI] emission toward only the protostar and [CII] emission toward the protostar and the red-shifted outflow. The [OI] emission has a width of $\sim$100 km s$^{-1}$ only in the blue-shifted velocity, suggesting an origin in shocked gas. The [CII] lines are narrow, consistent with an origin in a photodissociation region. Differential dust extinction from the envelope due to the inclination of the outflows is the most likely cause of the missing red-shifted [OI] emission. Fitting the [OI] line profile with two Gaussian components, we find one component at the source velocity with a width of $\sim$20 km s$^{-1}$ and another extremely broad component at -30 km s$^{-1}$ with a width of 87.5 km s$^{-1}$, the latter of which has not been seen in L1551 IRS 5. The kinematics of these two components resemble cavity shocks in molecular outflows and spot shocks in jets. Radiative transfer calculations of the [OI], high-J CO, and H$_2$O lines in the cavity shocks indicate that [OI] dominates the oxygen budget, making up more than 70% of the total gaseous oxygen abundance and suggesting [O]/[H] of $\sim$1.5$\times$10$^{-4}$. Attributing the extremely broad [OI] component to atomic winds, we estimate the intrinsic mass loss rate of (1.3$\pm$0.8)$\times$ 10$^{-6}$ M$_{\odot}$ yr$^{-1}$. The intrinsic mass loss rates derived from low-J CO, [OI], and HI are similar, supporting the model of momentum-conserving outflows, where the atomic wind carries most momentum and drives the molecular outflows.

preprint2021arXiv

Detection of a dense SiO jet in the evolved protostellar phase

Jets and outflows trace the accretion history of protostars. High-velocity molecular jets have been observed from several protostars in the early Class\,0 phase of star formation, detected with the high-density tracer SiO. Until now, no clear jet has been detected with SiO emission from isolated evolved Class\,I protostellar systems. We report a prominent dense SiO jet from a Class\,I source G205S3 (HOPS\,315: T$_{bol}$ $\sim$ 180 K, spectral index $\sim$ 0.417), with a moderately high mass-loss rate ($\sim$ 0.59 $\times$ 10$^{-6}$ M$_\odot$ yr$^{-1}$) estimated from CO emission. Together, these features suggest that G205S3 is still in a high accretion phase, similar to that expected of Class\,0 objects. We compare G205S3 to a representative Class\,0 system G206W2 (HOPS\,399) and literature Class\,0/I sources to explore the possible explanations behind the SiO emission seen at the later phase. We estimate a high inclination angle ($\sim$ 40$^\circ$) for G205S3 from CO emission, which may expose the infrared emission from the central core and mislead the spectral classification. However, the compact 1.3\,mm continuum, C$^{18}$O emission, location in the bolometric luminosity to sub-millimeter fluxes diagram, outflow force ($\sim$ 3.26 $\times$ 10$^{-5}$ M$_\odot$km s$^{-1}$/yr) are also analogous to that of Class\,I systems. We thus consider G205S3 to be at the very early phase of Class\,I, and in the late phase of ``high-accretion". The episodic ejection could be due to the presence of an unknown binary, a planetary companion, or dense clumps, where the required mass for such high accretion could be supplied by a massive circumbinary disk.

preprint2021arXiv

Molecular Cloud Cores with High Deuterium Fractions: Nobeyama Mapping Survey

We present the results of on-the-fly mapping observations of 44 fields containing 107 SCUBA-2 cores in the emission lines of molecules, N$_2$H$^+$, HC$_3$N, and CCS at 82$-$94 GHz using the Nobeyama 45-m telescope. This study aimed at investigating the physical properties of cores that show high deuterium fractions and might be close to the onset of star formation. We found that the distributions of the N$_2$H$^+$ and HC$_3$N line emissions are approximately similar to that of 850-$μ$m dust continuum emission, whereas the CCS line emission is often undetected or is distributed in a clumpy structure surrounding the peak position of the 850-$μ$m dust continuum emission. Occasionally (12%), we observe the CCS emission which is an early-type gas tracer toward the young stellar object, probably due to local high excitation. Evolution toward star formation does not immediately affect nonthermal velocity dispersion.

preprint2021arXiv

Quantifying Variability of YSOs in the Mid-IR Over Six Years with NEOWISE

Variability in young stellar objects (YSOs) can be caused by various time-dependent phenomena associated with star formation, including accretion rates, geometric changes in the circumstellar disks, stochastic hydromagnetic interactions between stellar surfaces and inner disk edges, reconnections within the stellar magnetosphere, and hot/cold spots on stellar surfaces. We uncover and characterize $\sim$1700 variables from a sample of $\sim$5400 YSOs in nearby low-mass star-forming regions using mid-IR light curves obtained from the 6.5-years NEOWISE All Sky Survey. The mid-IR variability traces a wide range of dynamical, physical, and geometrical phenomenon. We classify six types of YSO mid-IR variability based on their light curves: secular variability ($Linear, Curved, Periodic$) and stochastic variability ($Burst, Drop, Irregular$). YSOs in earlier evolutionary stages have higher fractions of variables and higher amplitudes for the variability, with the recurrence timescale of FUor-type outbursts (defined here as $Δ$W1 or $Δ$W2 $>1$ mag followed by inspection of candidates) of $\sim$1000 years in the early embedded protostellar phase. Known eruptive young stars and subluminous objects show fractions of variables similar to the fraction ($\sim55\%$) found in typical protostars, suggesting that these two distinct types are not distinct in variability over the 6.5-year timescale. Along with brightness variability, we also find a diverse range of secular color variations, which can be attributed to a competitive interplay between the variable accretion luminosity of the central source and the variable extinction by material associated with the accretion process.

preprint2021arXiv

V899 Mon: a peculiar eruptive young star close to the end of its outburst

V899 Mon is an eruptive young star showing characteristics of both FUors and EXors. It reached a peak brightness in 2010, then briefly faded in 2011, followed by a second outburst. We conducted multi-filter optical photometric monitoring, as well as optical and near-infrared spectroscopic observations of V899 Mon. The light curves and color-magnitude diagrams show that V899 Mon has been gradually fading after its second outburst peak in 2018, but smaller accretion bursts are still happening. Our spectroscopic observations taken with Gemini/IGRINS and VLT/MUSE show a number of emission lines, unlike during the outbursting stage. We used the emission line fluxes to estimate the accretion rate and found that it has significantly decreased compared to the outbursting stage. The mass loss rate is also weakening. Our 2D spectro-astrometric analysis of emission lines recovered jet and disk emission of V899 Mon. We found the emission from permitted metallic lines and the CO bandheads can be modeled well with a disk in Keplerian rotation, which also gives a tight constraint for the dynamical stellar mass of 2 ${M_{\odot}}$. After a discussion of the physical changes that led to the changes in the observed properties of V899 Mon, we suggest this object is finishing its second outburst.

preprint2020arXiv

ALMA ACA and Nobeyama observations of two Orion cores in deuterated molecular lines

We mapped two molecular cloud cores in the Orion A cloud with the ALMA ACA 7-m Array and with the Nobeyama 45-m radio telescope. These cores have bright N$_2$D$^+$ emission in single-pointing observations with the Nobeyama 45-m radio telescope, have relatively high deuterium fraction, and are thought to be close to the onset of star formation. One is a star-forming core, and the other is starless. These cores are located along filaments observed in N$_2$H$^+$, and show narrow linewidths of 0.41 km s$^{-1}$ and 0.45 km s$^{-1}$ in N$_2$D$^+$, respectively, with the Nobeyama 45-m telescope. Both cores were detected with the ALMA ACA 7m Array in the continuum and molecular lines at Band 6. The starless core G211 shows clumpy structure with several sub-cores, which in turn show chemical differences. Also, the sub-cores in G211 have internal motions that are almost purely thermal. The starless sub-core G211D, in particular, shows a hint of the inverse P Cygni profile, suggesting infall motion. The star-forming core G210 shows an interesting spatial feature of two N$_2$D$^+$ peaks of similar intensity and radial velocity located symmetrically with respect to the single dust continuum peak. One interpretation is that the two N$_2$D$^+$ peaks represent an edge-on pseudo-disk. The CO outflow lobes, however, are not directed perpendicular to the line connecting both N$_2$D$^+$ peaks.

preprint2020arXiv

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): I. Detection of New Hot Corinos with ACA

We report the detection of four new hot corino sources, G211.47-19.27S, G208.68-19.20N1, G210.49-19.79W and G192.12-11.10 from a survey study of Planck Galactic Cold Clumps in the Orion Molecular Cloud Complex with the Atacama Compact Array (ACA). Three sources had been identified as low mass Class 0 protostars in the Herschel Orion Protostar Survey (HOPS). One source in the lambda Orionis region is firstly reported as a protostellar core. We have observed abundant complex organic molecules (COMs), primarily methanol but also other oxygen-bearing COMs (in G211.47-19.27S and G208.68-19.20N1) and the molecule of prebiotic interest NH2CHO (in G211.47-19.27S), signifying the presence of hot corinos. While our spatial resolution is not sufficient for resolving most of the molecular emission structure, the large linewidth and high rotational temperature of COMs suggest that they likely reside in the hotter and innermost region immediately surrounding the protostar. In G211.47-19.27S, the D/H ratio of methanol ([CH2DOH]/[CH3OH]) and the 12C/13C ratio of methanol ([CH3OH]/[13CH3OH]) are comparable to those of other hot corinos. Hydrocarbons and long carbon-chain molecules such as c-C3H2 and HCCCN are also detected in the four sources, likely tracing the outer and cooler molecular envelopes.

preprint2020arXiv

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- I. Survey description and a first look at G9.62+0.19

The "ATOMS," standing for {\it ALMA Three-millimeter Observations of Massive Star-forming regions}, survey has observed 146 active star forming regions with ALMA Band 3, aiming to systematically investigate the spatial distribution of various dense gas tracers in a large sample of Galactic massive clumps, to study the roles of stellar feedback in star formation, and to characterize filamentary structures inside massive clumps. In this work, the observations, data analysis, and example science of the "ATOMS" survey are presented, using a case study for the G9.62+0.19 complex. Toward this source, some transitions, commonly assumed to trace dense gas, including CS $J = 2-1$, HCO$^+$ $J = 1-0$ and HCN $J = 1-0$, are found to show extended gas emission in low density regions within the clump; less than 25\% of their emission is from dense cores. SO, CH$_3$OH, H$^{13}$CN and HC$_3$N show similar morphologies in their spatial distributions and reveal well the dense cores. Widespread narrow SiO emission is present (over $\sim$1 pc), which may be caused by slow shocks from large--scale colliding flows or H{\sc ii} regions. Stellar feedback from an expanding H{\sc ii} region has greatly reshaped the natal clump, significantly changed the spatial distribution of gas, and may also account for the sequential high-mass star formation in the G9.62+0.19 complex. The ATOMS survey data can be jointly analyzed with other survey data, e.g., "MALT90", "Orion B", "EMPIRE", "ALMA\_IMF", and "ALMAGAL", to deepen our understandings of "dense gas" star formation scaling relations and massive proto-cluster formation.

preprint2020arXiv

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- II. Compact objects in ACA observations and star formation scaling relations

We report studies of the relationships between the total bolometric luminosity ($L_{\rm bol}$ or $L_{\rm TIR}$) and the molecular line luminosities of $J=1-0$ transitions of H$^{13}$CN, H$^{13}$CO$^+$, HCN, and HCO$^+$ with data obtained from ACA observations in the "ATOMS" survey of 146 active Galactic star forming regions. The correlations between $L_{\rm bol}$ and molecular line luminosities $L'_{\rm mol}$ of the four transitions all appear to be approximately linear. Line emission of isotopologues shows as large scatters in $L_{\rm bol}$-$L'_{\rm mol}$ relations as their main line emission. The log($L_{\rm bol}$/$L'_{\rm mol}$) for different molecular line tracers have similar distributions. The $L_{\rm bol}$-to-$L'_{\rm mol}$ ratios do not change with galactocentric distances ($R_{\rm GC}$) and clump masses ($M_{\rm clump}$). The molecular line luminosity ratios (HCN-to-HCO$^+$, H$^{13}$CN-to-H$^{13}$CO$^+$, HCN-to-H$^{13}$CN and HCO$^+$-to-H$^{13}$CO$^+$) all appear constant against $L_{\rm bol}$, dust temperature ($T_{\rm d}$), $M_{\rm clump}$ and $R_{\rm GC}$. Our studies suggest that both the main lines and isotopologue lines are good tracers of the total masses of dense gas in Galactic molecular clumps. The large optical depths of main lines do not affect the interpretation of the slopes in star formation relations. We find that the mean star formation efficiency (SFE) of massive Galactic clumps in the "ATOMS" survey is reasonably consistent with other measures of the SFE for dense gas, even those using very different tracers or examining very different spatial scales.

preprint2020arXiv

High-resolution spectroscopic monitoring observations of FU Orionis-type object, V960 Mon

We present the results of high-resolution (R $\ge$ 30,000) optical and near-infrared spectroscopic monitoring observations of a FU Orionis-type object, V960 Mon, which underwent an outburst in 2014 November. We have monitored this object with the Bohyunsan Optical Echelle Spectrograph (BOES) and the Immersion GRating INfrared Spectrograph (IGRINS) since 2014 December. Various features produced by a wind, disk, and outflow/jet were detected. The wind features varied over time and continually weakened after the outburst. We detected double-peaked line profiles in the optical and near-infrared, and the line widths tend to decrease with increasing wavelength, indicative of Keplerian disk rotation. The disk features in the optical and near-infrared spectra fit well with G-type and K-type stellar spectra convolved with a kernel to account for the maximum projected disk rotation velocity of about 40.3$\pm$3.8 km s$^{-1}$ and 36.3$\pm$3.9 km s$^{-1}$, respectively. We also report the detection of [S II] and H$_{2}$ emission lines, which are jet/outflow tracers and rarely found in FUors.

preprint2020arXiv

Modeling Nitrogen Fractionation in the Protoplanetary Disk around TW Hya: Model Constraints on Grain Population and Carbon-to-Oxygen Elemental Abundance Ratio

Observations conducted using the Atacama Large Millimeter/submillimeter Array on the protoplanetary disk around TW Hya show the nitrogen fractionation of HCN molecules in HC$^{14}$N/HC$^{15}$N $\sim$120 at a radius of $\sim$20 AU. In this study, we investigated the physical and chemical conditions that control this nitrogen fractionation process. To this end, a new disk model was developed, in which the isotope-selective photodissociation of N$_2$ and isotope-exchange chemical reactions have been incorporated. Our model can successfully reproduce the observed HCN column density when the elemental abundances of the gas-phase carbon and oxygen are depleted by two orders of magnitude relative to those in the interstellar medium and carbon is more abundant than oxygen ([C/O]$_{\rm elem}>$ 1). The isotope-selective photodissociation of N$_2$ is the dominant nitrogen fractionation process in our models. The observed HC$^{14}$N/HC$^{15}$N ratio, which increases outwards, can also be reproduced by the model by assuming that the small dust grains in the atmosphere of the outer disk are depleted more than those in the inner disk. This is consistent with grain evolution models, according to which small dust grains are continuously replenished in the inner disk due to fragmentation of the large dust grains that radially drift from the outer disk.

preprint2020arXiv

Molecular Cloud Cores with High Deuterium Fraction: Nobeyama Single-Pointing Survey

We present the results of a single-pointing survey of 207 dense cores embedded in Planck Galactic Cold Clumps distributed in five different environments ($λ$ Orionis, Orion A, B, Galactic plane, and high latitudes) to identify dense cores on the verge of star formation for the study of the initial conditions of star formation. We observed these cores in eight molecular lines at 76-94 GHz using the Nobeyama 45-m telescope. We find that early-type molecules (e.g., CCS) have low detection rates and that late-type molecules (e.g., N$_2$H$^+$, c-C$_3$H$_2$) and deuterated molecules (e.g., N$_2$D$^+$, DNC) have high detection rates, suggesting that most of the cores are chemically evolved. The deuterium fraction (D/H) is found to decrease with increasing distance, indicating that it suffers from differential beam dilution between the D/H pair of lines for distant cores ($>$1 kpc). For $λ$ Orionis, Orion A, and B located at similar distances, D/H is not significantly different, suggesting that there is no systematic difference in the observed chemical properties among these three regions. We identify at least eight high D/H cores in the Orion region and two at high latitudes, which are most likely to be close to the onset of star formation. There is no clear evidence of the evolutionary change in turbulence during the starless phase, suggesting that the dissipation of turbulence is not a major mechanism for the beginning of star formation as judged from observations with a beam size of 0.04 pc.

preprint2020arXiv

Radiative Transfer modeling of EC 53: An Episodically Accreting Class I Young Stellar Object

In the episodic accretion scenario, a large fraction of the protostellar mass accretes during repeated and large bursts of accretion. Since outbursts on protostars are typically identified at specific wavelengths, interpreting these outbursts requires converting this change in flux to a change in total luminosity. The Class I young stellar object EC 53 in the Serpens Main cloud has undergone repeated increases in brightness at 850 $μ$m that are likely caused by bursts of accretion. In this study, we perform two- and three-dimensional continuum radiative transfer modeling to quantify the internal luminosity rise in EC 53 that corresponds to the factor of $\sim$1.5 enhancement in flux at 850 $μ$m. We model the spectral energy distribution and radial intensity profile in both the quiescent and outburst phases. The internal luminosity in the outburst phase is $\sim 3.3$ times brighter than the luminosity in the quiescent phase. The radial intensity profile analysis demonstrates that the detected sub-mm flux variation of EC 53 comes from the heated envelope by the accretion burst. We also find that the role of external heating of the EC 53 envelope by the interstellar radiation field is insignificant.

preprint2020arXiv

The origin of a distributed stellar population in the star-forming region W4

Stellar kinematics provides the key to understanding the formation process and dynamical evolution of stellar systems. Here, we present a kinematic study of the massive star-forming region W4 in the Cassiopeia OB6 association using the Gaia Data Release 2 and high-resolution optical spectra. This star-forming region is composed of a core cluster (IC 1805) and a stellar population distributed over 20 pc, which is a typical structural feature found in many OB associations. According to a classical model, this structural feature can be understood in the context of the dynamical evolution of a star cluster. The core-extended structure exhibits internally different kinematic properties. Stars in the core have an almost isotropic motion, and they appear to reach virial equilibrium given their velocity dispersion (0.9 +/- 0.3 km/s) comparable to that in a virial state (~0.8 km/s). On the other hand, the distributed population shows a clear pattern of radial expansion. From the N-body simulation for the dynamical evolution of a model cluster in subvirial state, we reproduce the observed structure and kinematics of stars. This model cluster experiences collapse for the first 2 Myr. Some members begin to radially escape from the cluster after the initial collapse, eventually forming a distributed population. The internal structure and kinematics of the model cluster appear similar to those of W4. Our results support the idea that the stellar population distributed over 20 pc in W4 originate from the dynamical evolution of IC 1805.

preprint2020arXiv

The Relationship between Mid-Infrared and Sub-Millimetre Variability of Deeply Embedded Protostars

We study the relationship between the mid-infrared and sub-mm variability of deeply embedded protostars using the multi-epoch data from the Wide Infrared Survey Explorer ($WISE$/NEOWISE) and the ongoing James Clerk Maxwell Telescope (JCMT) transient survey. Our search for signs of stochastic (random) and/or secular (roughly monotonic in time) variability in a sample of 59 young stellar objects (YSOs) revealed that 35 are variable in at least one of the two surveys. This variability is dominated by secular changes. Of those objects with secular variability, 14 objects ($22\%$ of the sample) show correlated secular variability over mid-IR and sub-mm wavelengths. Variable accretion is the likely mechanism responsible for this type of variability. Fluxes of YSOs that vary in both wavelengths follow a relation of $\log_{10} F_{4.6}(t)=η\log_{10} F_{850}(t)$ between the mid-IR and sub-mm, with $η=5.53\pm0.29$. This relationship arises from the fact that sub-mm fluxes respond to the dust temperature in the larger envelope whereas the mid-IR emissivity is more directly proportional to the accretion luminosity. The exact scaling relation, however, depends on the structure of the envelope, the importance of viscous heating in the disc, and dust opacity laws.

preprint2020arXiv

Three dimensional projection effects on chemistry in a Planck galactic cold clump

Offsets of molecular line emission peaks from continuum peaks are very common but frequently difficult to explain with a single spherical cloud chemical model. We propose that the spatial projection effects of an irregular three dimensional (3D) cloud structure can be a solution. This work shows that the idea can be successfully applied to the Planck cold clump G224.4-0.6 by approximating it with four individual spherically symmetric cloud cores whose chemical patterns overlap with each other to produce observable line maps. With the empirical physical structures inferred from the observation data of this clump and a gas-grain chemical model, the four cores can satisfactorily reproduce its 850 $μ$m continuum map and the diverse peak offsets of CCS, HC$_3$N and N$_2$H$^+$ simultaneously at chemical ages of about $8\times 10^5\sim 3\times 10^6$ yrs. The 3D projection effects on chemistry has the potential to explain such asymmetrical distributions of chemicals in many other molecular clouds.

preprint2019arXiv

The circumstellar environment around the embedded protostar EC 53

EC53 is an embedded protostar with quasi-periodic emission in the near-IR and sub-mm. We use ALMA high-resolution observations of continuum and molecular line emission to describe the circumstellar environment of EC 53. The continuum image reveals a disk with a flux that suggests a mass of 0.075 Msun, much less than the estimated mass in the envelope, and an in-band spectral index that indicates grain growth to centimeter sizes. Molecular lines trace the outflow cavity walls, infalling and rotating envelope, and/or the Keplerian disk. The rotation profile of the C17O 3-2 line emission cannot isolate the Keplerian motion clearly although the lower limit of the protostellar mass can be calculated as 0.3 +- 0.1 Msun if the Keplerian motion is adopted. The weak CH3OH emission, which is anti-correlated with the HCO+ 4-3 line emission, indicates that the water snow line is more extended than what expected from the current luminosity, attesting to bygone outburst events. The extended snow line may persist for longer at the disk surface because the lower density increases the freeze-out timescale of methanol and water.

preprint2016arXiv

Astrochemical Properties of Planck Cold Clumps

We observed thirteen Planck cold clumps with the James Clerk Maxwell Telescope/SCUBA-2 and with the Nobeyama 45 m radio telescope. The N$_2$H$^+$ distribution obtained with the Nobeyama telescope is quite similar to SCUBA-2 dust distribution. The 82 GHz HC$_3$N, 82 GHz CCS, and 94 GHz CCS emission are often distributed differently with respect to the N$_2$H$^+$ emission. The CCS emission, which is known to be abundant in starless molecular cloud cores, is often very clumpy in the observed targets. We made deep single-pointing observations in DNC, HN$^{13}$C, N$_2$D$^+$, cyclic-C$_3$H$_2$ toward nine clumps. The detection rate of N$_2$D$^+$ is 50\%. Furthermore, we observed the NH$_3$ emission toward 15 Planck cold clumps to estimate the kinetic temperature, and confirmed that most of targets are cold ($\lesssim$ 20 K). In two of the starless clumps observe, the CCS emission is distributed as it surrounds the N$_2$H$^+$ core (chemically evolved gas), which resembles the case of L1544, a prestellar core showing collapse. In addition, we detected both DNC and N$_2$D$^+$. These two clumps are most likely on the verge of star formation. We introduce the Chemical Evolution Factor (CEF) for starless cores to describe the chemical evolutionary stage, and analyze the observed Planck cold clumps.

preprint2016arXiv

Discovery of an extremely wide-angle bipolar outflow in AFGL 5142

Most bipolar outflows are associated with individual young stellar objects and have small opening angles. Here we report the discovery of an extremely wide-angle ($\sim$180$\arcdeg$) bipolar outflow ("EWBO") in a cluster forming region AFGL 5142 from low-velocity emission of the HCN (3-2) and HCO$^{+}$ (3-2) lines. This bipolar outflow is along a north-west to south-east direction with a line-of-sight flow velocity of about 3 km~s$^{-1}$ and is spatially connected to the high-velocity jet-like outflows. It seems to be a collection of low-velocity material entrained by the high-velocity outflows due to momentum feedback. The total ejected mass and mass loss rate due to both high velocity jet-like outflows and the "EWBO" are $\sim$24.5 M$_{\sun}$ and $\sim1.7\times10^{-3}$ M$_{\sun}$~yr$^{-1}$, respectively. Global collapse of the clump is revealed by the "blue profile" in the HCO$^{+}$ (1-0) line. A hierarchical network of filaments was identified in NH$_{3}$ (1,1) emission. Clear velocity gradients of the order of 10 km~s$^{-1}$~pc$^{-1}$ are found along filaments, indicating gas inflow along the filaments. The sum of the accretion rate along filaments and mass infall rate along the line of sight is $\sim$3.1$\times10^{-3}$ M$_{\sun}$~yr$^{-1}$, which exceeds the total mass loss rate, indicating that the central cluster is probably still gaining mass. The central cluster is highly fragmented and 22 condensations are identified in 1.1 mm continuum emission. The fragmentation process seems to be determined by thermal pressure and turbulence. The magnetic field may not play an important role in fragmentation.

preprint2016arXiv

IGRINS spectroscopy of Class I sources: IRAS 03445+3242 and IRAS 04239+2436

We have detected molecular and atomic line emission from the hot and warm disks of two Class I sources, IRAS 03445+3242 and IRAS 04239+2436 using the high resolution Immersion GRating INfrared Spectrograph (IGRINS). CO overtone band transitions and near-IR lines of Na I and Ca I, all in emission, trace the hot inner disk while CO rovibrational absorption spectra of the first overtone transition trace the warm gas within the inner few AU of the disk. The emission-line profiles for both sources show evidence for Keplerian disks. A thin Keplerian disk with power-law temperature and column density profiles with a projected rotational velocity of $\sim$60--75 km s$^{-1}$ and a gas temperature of $\sim$3500 K at the innermost annulus can reproduce the CO overtone band emission. Na I and Ca I emission lines also arise from this disk, but they show complicated line features possibly affected by photospheric absorption lines. Multi-epoch observations show asymmetric variations of the line profiles on one-year (CO overtone bandhead and atomic lines for IRAS 03445+3242) or on one-day (atomic lines for IRAS 04239+2436) time scales, implying non-axisymmetric features in disks. The narrow CO rovibrational absorption spectra ($v$=0$\rightarrow$2) indicate that both warm ($>$ 150 K) and cold ($\sim$20--30 K) CO gas are present along the line of sight to the inner disk. This study demonstrates the power of IGRINS as a tool for studies of the sub-AU scale hot and AU-scale warm protoplanetary disks with its simultaneous coverage of the full H and K bands with high spectral resolution ($R$= 45,000) allowing many aspects of the sources to be investigated at once.

preprint2016arXiv

The CDF Archive: Herschel PACS and SPIRE Spectroscopic Data Pipeline and Products for Protostars and Young Stellar Objects

We present the COPS-DIGIT-FOOSH (CDF) Herschel spectroscopy data product archive, and related ancillary data products, along with data fidelity assessments, and a user-created archive in collaboration with the Herschel-PACS and SPIRE ICC groups. Our products include datacubes, contour maps, automated line fitting results, and best 1-D spectra products for all protostellar and disk sources observed with PACS in RangeScan mode for two observing programs: the DIGIT Open Time Key Program (KPOT_nevans_1 and SDP_nevans_1; PI: N. Evans), and the FOOSH Open Time Program (OT1_jgreen02_2; PI: J. Green). In addition, we provide our best SPIRE-FTS spectroscopic products for the COPS Open Time Program (OT2_jgreen02_6; PI: J. Green) and FOOSH sources. We include details of data processing, descriptions of output products, and tests of their reliability for user applications. We identify the parts of the dataset to be used with caution. The resulting absolute flux calibration has improved in almost all cases. Compared to previous reductions, the resulting rotational temperatures and numbers of CO molecules have changed substantially in some sources. On average, however, the rotational temperatures have not changed substantially (< 2%), but the number of warm (Trot ~ 300 K) CO molecules has increased by about 18%.

preprint2016arXiv

The inflow signature toward different evolutionary phases of massive star formation

We analyzed both HCN J=1-0 and HNC J=1-0 line profiles to study the inflow motions in different evolutionary stages of massive star formation: 54 infrared dark clouds (IRDCs), 69 high-mass protostellar object (HMPOs), and 54 ultra-compact HII regions (UCHIIs). The inflow asymmetry in HCN spectra seems to be prevalent throughout all the three evolutionary phases, with IRDCs showing the largest excess in blue profile. In the case of HNC spectra, the prevalence of blue sources does not appear, excepting for IRDCs. We suggest that this line is not appropriate to trace inflow motion in evolved stages of massive star formation because the abundance of HNC decreases at high temperatures. This result spotlights the importance of considering chemistry in the dynamics study of massive star-forming regions. The fact that the IRDCs show the highest blue excess in both transitions indicates that the most active inflow occurs in the early phase of star formation, i.e., the IRDC phase rather than in the later phases. However, mass is still inflowing onto some UCHIIs. We also found that the absorption dips of the HNC spectra in 6 out of 7 blue sources are red-shifted relative to their systemic velocities. These red-shifted absorption dips may indicate global collapse candidates, although mapping observations with better resolution are needed to examine this feature in more detail.

preprint2015arXiv

"Dust, Ice, and Gas In Time" (DIGIT) Herschel Observations of GSS30-IRS1 in Ophiuchus

As a part of the "Dust, Ice, and Gas In Time" (DIGIT) key program on Herschel, we observed GSS30-IRS1, a Class I protostar located in Ophiuchus (d = 120 pc), with Herschel/Photodetector Array Camera and Spectrometer (PACS). More than 70 lines were detected within a wavelength range from 50 micron to 200 micron, including CO, H2O, OH, and two atomic [O I] lines at 63 and 145 micron. The [C II] line, known as a tracer of externally heated gas by the interstellar radiation field, is also detected at 158 micron. All lines, except [O I] and [C II], are detected only at the central spaxel of 9.4" X 9.4". The [O I] emissions are extended along a NE-SW orientation, and the [C II] line is detected over all spaxels, indicative of external PDR. The total [C II] intensity around GSS30 reveals that the far-ultraviolet radiation field is in the range of 3 to 20 G0, where G0 is in units of the Habing Field, 1.6 X 10^{-3} erg cm^{-2} s^{-1}. This enhanced external radiation field heats the envelope of GSS30-IRS1, causing the continuum emission to be extended, unlike the molecular emission. The best-fit continuum model of GSS30-IRS1 with the physical structure including flared disk, envelope, and outflow shows that the internal luminosity is 10 Lsun, and the region is externally heated by a radiation field enhanced by a factor of 130 compared to the standard local interstellar radiation field.

preprint2015arXiv

An Ammonia Spectral Map of the L1495-B218 Filaments in the Taurus Molecular Cloud : I. Physical Properties of Filaments and Dense cores

We present deep NH$_3$ observations of the L1495-B218 filaments in the Taurus molecular cloud covering over a 3 degree angular range using the K-band focal plane array on the 100m Green Bank Telescope. The L1495-B218 filaments form an interconnected, nearby, large complex extending over 8 pc. We observed NH$_3$ (1,1) and (2,2) with a spectral resolution of 0.038 km/s and a spatial resolution of 31$"$. Most of the ammonia peaks coincide with intensity peaks in dust continuum maps at 350 $μ$m and 500 $μ$m. We deduced physical properties by fitting a model to the observed spectra. We find gas kinetic temperatures of 8 $-$ 15 K, velocity dispersions of 0.05 $-$ 0.25 km/s, and NH$_3$ column densities of 5$\times$10$^{12}$ $-$ 1$\times$10$^{14}$ cm$^{-2}$. The CSAR algorithm, which is a hybrid of seeded-watershed and binary dendrogram algorithms, identifies a total of 55 NH$_3$ structures including 39 leaves and 16 branches. The masses of the NH$_3$ sources range from 0.05 M$_\odot$ to 9.5 M$_\odot$. The masses of NH$_3$ leaves are mostly smaller than their corresponding virial mass estimated from their internal and gravitational energies, which suggests these leaves are gravitationally unbound structures. 9 out of 39 NH$_3$ leaves are gravitationally bound and 7 out of 9 gravitationally bound NH$_3$ leaves are associated with star formation. We also found that 12 out of 30 gravitationally unbound leaves are pressure-confined. Our data suggest that a dense core may form as a pressure-confined structure, evolve to a gravitationally bound core, and undergo collapse to form a protostar.

preprint2015arXiv

Color Variability of HBC 722 in the Post-Outburst Phases

We carried out photometric observations for HBC 722 in SDSS r, i and z bands from 2011 April to 2013 May with a Camera for Quasars in Early uNiverse attached to the 2.1m Otto Struve telescope at McDonald Observatory. The post-outburst phenomena were classified into five phases according to not only brightness but also color variations, which might be caused by physical changes in the emitting regions of optical and near-infrared bands. A series of spectral energy distribution (SED) is presented to support color variations and track the time evolution of SED in optical/near-infrared bands after the outburst. Given two years of data, possible periodicities of r, i and z bands were checked. We found out three families of signals around ~6, ~10 and ~1 days in three bands, which is broadly consistent with Green et al. (2013). We also examined short term variability (intra-day and day scales) to search for evidences of flickering by using the micro-variability method. We found clear signs of day scale variability and weak indications of intra-day scale fluctuations, which implies that the flickering event occurs in HBC 722 after outburst.

preprint2015arXiv

Detection of Infall in the Protostar B335 with ALMA

Observations of the isolated globule B335 with ALMA have yielded absorption features against the continuum that are redshifted from the systemic velocity in both HCN and HCO$^+$ lines. These features provide unambiguous evidence for infall toward a central luminosity source. Previously developed models of inside-out collapse can match the observed line profiles of HCN and HCO$^+$ averaged over the central 50 AU. At the new distance of 100 pc, the inferred infall radius is 0.012 pc, the mass infall rate is $3 \times 10^{-6}$ solar masses per year, the age is 50,000 years, and the accumulated mass in the central zone is 0.15 solar masses, most of which must be in the star or in parts of a disk that are opaque at 0.8 mm. The continuum detection indicates an optically thin mass (gas and dust) of only $7.5\times 10^{-4}$ solar masses in the central region, consistent with only a very small disk mass.

preprint2015arXiv

High Resolution Optical and NIR Spectra of HBC 722

We present the results of high resolution (R$\ge$30,000) optical and near-IR spectroscopic monitoring observations of HBC 722, a recent FU Orionis object that underwent an accretion burst in 2010. We observed HBC 722 in optical/near-IR with the BOES, HET-HRS, and IGRINS spectrographs, at various points in the outburst. We found atomic lines with strongly blueshifted absorption features or P Cygni profiles, both evidence of a wind driven by the accretion. Some lines show a broad double-peaked absorption feature, evidence of disk rotation. However, the wind-driven and disk-driven spectroscopic features are anti-correlated in time; the disk features became strong as the wind features disappeared. This anti-correlation might indicate that the rebuilding of the inner disk was interrupted by the wind pressure during the first two years. The Half-Width at Half-Depth (HWHD) of the double-peaked profiles decreases with wavelength, indicative of the Keplerian rotation; the optical spectra with the disk feature are fitted by a G5 template stellar spectrum convolved with a rotation velocity of 70 km s$^{-1}$ while the near-IR disk features are fitted by a K5 template stellar spectrum convolved with a rotation velocity of 50 km s$^{-1}$. Therefore, the optical and near-IR spectra seem to trace the disk at 39 and 76 $\textit{R}_{\odot}$, respectively. We fit a power-law temperature distribution in the disk, finding an index of 0.8, comparable to optically thick accretion disk models.

preprint2015arXiv

Infrared and Radio observations of a small group of protostellar objects in the molecular core, L1251-C

We present a multi-wavelength observational study of a low-mass star-forming region, L1251-C, with observational results at wavelengths from the near-infrared to the millimeter. Spitzer Space Telescope observations confirmed that IRAS 22343+7501 is a small group of protostellar objects. The extended emission to east-west direction with its intensity peak at the center of L1251A has been detected at 350 and 850 micron with the CSO and JCMT telescopes, tracing dense envelope materials around L1251A. The single-dish data from the KVN and TRAO telescopes show inconsistencies between the intensity peaks of several molecular line emission and that of the continuum emission, suggesting complex distributions of molecular abundances around L1251A. The SMA interferometer data, however, show intensity peaks of CO 2-1 and 13CO 2-1 located at the position of IRS 1, which is both the brightest source in IRAC image and the weakest source in the 1.3 mm dust continuum map. IRS 1 is the strongest candidate for the driving source of the newly detected compact CO 2-1 outflow. Over the whole region (14' by 14') of L125l-C, 3 Class I and 16 Class II sources have been detected, including three YSOs in L1251A. A comparison with the average projected distance among 19 YSOs in L1251-C and that among 3 YSOs in L1251A suggests L1251-C is an example of low-mass cluster formation, where protostellar objects are forming in a small group.

preprint2015arXiv

Planck cold clumps in the $λ$ Orionis complex: I. Discovery of an extremely young Class 0 protostellar object and a proto-brown dwarf candidate in a bright rimmed clump PGCC G192.32-11.88

We are performing a series of observations with ground-based telescopes toward Planck Galactic cold clumps (PGCCs) in the $λ$ Orionis complex in order to systematically investigate the effects of stellar feedback. In the particular case of PGCC G192.32-11.88, we discovered an extremely young Class 0 protostellar object (G192N) and a proto-brown dwarf candidate (G192S). G192N and G192S are located in a gravitationally bound bright-rimmed clump. The velocity and temperature gradients seen in line emission of CO isotopologues indicate that PGCC G192.32-11.88 is externally heated and compressed. G192N probably has the lowest bolometric luminosity ($\sim0.8$ L$_{\sun}$) and accretion rate (6.3$\times10^{-7}$ M$_{\sun}$~yr$^{-1}$) when compared with other young Class 0 sources (e.g. PACS Bright Red sources (PBRs)) in the Orion complex. It has slightly larger internal luminosity ($0.21\pm0.01$ L$_{\sun}$) and outflow velocity ($\sim$14 km~s$^{-1}$) than the predictions of first hydrostatic cores (FHSCs). G192N might be among the youngest Class 0 sources, which are slightly more evolved than a FHSC. Considering its low internal luminosity ($0.08\pm0.01$ L$_{\odot}$) and accretion rate (2.8$\times10^{-8}$ M$_{\sun}$~yr$^{-1}$), G192S is an ideal proto-brown dwarf candidate. The star formation efficiency ($\sim$0.3\%-0.4\%) and core formation efficiency ($\sim$1\%) in PGCC G192.32-11.88 are significantly smaller than in other giant molecular clouds or filaments, indicating that the star formation therein is greatly suppressed due to stellar feedback.

preprint2015arXiv

Radio Observations of the Star Formation Activities in the NGC 2024 FIR 4 Region

Star formation activities in the NGC 2024 FIR 4 region were studied by imaging centimeter continuum sources and water maser sources using several archival data sets from the Very Large Array. The continuum source VLA 9 is elongated in the northwest-southeast direction, consistent with the FIR 4 bipolar outflow axis, and has a flat spectrum in the 6.2-3.6 cm interval. The three water maser spots associated with FIR 4 are also distributed along the outflow axis. One of the spots is located close to VLA 9, and another one is close to an X-ray source. Examinations of the positions of compact objects in this region suggest that the FIR 4 cloud core contains a single low-mass protostar. VLA 9 is the best indicator of the protostellar position. VLA 9 may be a radio thermal jet driven by this protostar, and it is unlikely that FIR 4 contains a high-mass young stellar object (YSO). A methanol 6.7 GHz maser source is located close to VLA 9, at a distance of about 100 AU. The FIR 4 protostar must be responsible for the methanol maser action, which suggests that methanol class II masers are not necessarily excited by high-mass YSOs. Also discussed are properties of other centimeter continuum sources in the field of view and the water masers associated with FIR 6n. Some of the continuum sources are radio thermal jets, and some are magnetically active young stars.

preprint2015arXiv

Stochastic Model of the Spin Distribution of Dark Matter Halos

We employ a stochastic approach to probing the origin of the log-normal distributions of halo spin in N-body simulations. After analyzing spin evolution in halo merging trees, it was found that a spin change can be characterized by a stochastic random walk of angular momentum. Also, spin distributions generated by random walks are fairly consistent with those directly obtained from N-body simulations. We derived a stochastic differential equation from a widely used spin definition and measured the probability distributions of the derived angular momentum change from a massive set of halo merging trees. The roles of major merging and accretion are also statistically analyzed in evolving spin distributions. Several factors (local environment, halo mass, merging mass ratio, and redshift) are found to influence the angular momentum change. The spin distributions generated in the mean-field or void regions tend to shift slightly to a higher spin value compared with simulated spin distributions, which seems to be caused by the correlated random walks. We verified the assumption of randomness in the angular momentum change observed in the N-body simulation and detected several degrees of correlation between walks, which may provide a clue for the discrepancies between the simulated and generated spin distributions in the voids. However, the generated spin distributions in the group and cluster regions successfully match the simulated spin distribution. We also demonstrated that the log-normality of the spin distribution is a natural consequence of the stochastic differential equation of the halo spin, which is well described by the Geometric Brownian Motion model.

preprint2015arXiv

The HCN/HNC abundance ratio toward different evolutionary phases of massive star formation

Using the HCN and HNC J=1--0 line observations, the abundance ratio of HCN/HNC has been estimated for different evolutionary stages of massive star formation: Infrared dark clouds (IRDCs), High-mass protostellar object (HMPOs), and Ultra-compact HII regions (UCHIIs). IRDCs were divided into `quiescent IRDC cores' and `active IRDC cores', depending on star formation activity. The HCN/HNC ratio is known to be higher at active and high temperature regions related to ongoing star formation, compared to cold and quiescent regions. Our observations toward 8 quiescent IRDC cores, 16 active IRDC cores, 23 HMPOs, and 31 UCHIIs show consistent results; the ratio is 0.97~($\pm~0.10$), 2.65~($\pm~0.88$), 4.17~($\pm~1.03$) and 8.96~($\pm~3.32$) in these respective evolutionary stages, increasing from quiescent IRDC cores to UCHIIs. The change of the HCN/HNC abundance ratio, therefore, seems directly associated with the evolutionary stages of star formation, which have different temperatures. One suggested explanation for this trend is the conversion of HNC to HCN, which occurs effectively at higher temperatures. To test the explanation, we performed a simple chemical model calculation. In order to fit the observed results, the energy barrier of the conversion must be much lower than the value provided by theoretical calculations.

preprint2015arXiv

The warm CO gas along the UV-heated outflow cavity walls: a possible interpretation for the Herschel/PACS CO spectra of embedded YSOs

A fraction of the mid-$J$ ($J$= 14--13 to $J$= 24--23) CO emission detected by the \textit{Herschel}/PACS observations of embedded young stellar objects (YSOs) has been attributed to the UV-heated outflow cavity walls. We have applied our newly developed self-consistent models of Photon-Dominated Region (PDR) and non-local thermal equilibrium line Radiative transfer In general Grid (RIG) to the \textit{Herschel} FIR observations of 27 low mass YSOs and one intermediate mass YSO, NGC7129-FIRS2. When the contribution of the hot component (traced by transitions of $J> 24$) is removed, the rotational temperature of the warm component is nearly constant with $\sim250$ K. This can be reproduced by the outflow cavity wall ($n \geq 10^6\, \mathrm{cm}^{-3}$, $\log G_{0}/n \geq-4.5$, $\mathrm{log} G_0\ge 3$, $T_{\rm gas} \ge 300 $K, and X(CO)$ \ge 10^{-5}$) heated by a UV radiation field with a black body temperature of 15,000 K or 10,000 K. However, a shock model combined with an internal PDR will be required to determine the quantitative contribution of a PDR relative to a shock to the mid-$J$ CO emission.

preprint2014arXiv

A PDR model for the FIR mid-$J$ CO ladder with universal rotational temperature in star forming region

A photon dominated region (PDR) is one of the leading candidate mechanisms for the origin of the warm CO gas with near universal $\sim$300~K rotational temperature inferred from the CO emission detected towards embedded protostars by Herschel/PACS. We have developed a PDR model in general coordinates, where we can use the most adequate coordinate system for an embedded protostar having outflow cavity walls, to solve chemistry and gas energetics self-consistently for given UV radiation fields with different spectral shapes. Simple 1D tests and applications show that FIR mid-$J$ ($14\leq J\leq 24$) CO lines are emitted from near the surface of a dense region exposed to high UV fluxes. We apply our model to HH46 and find the UV-heated outflow cavity wall can reproduce the mid-$J$, CO transitions observed by Herschel/PACS. A model with UV radiation corresponding to a blackbody of 10,000~K results in the rotational temperature lower than 300~K, while models with the Draine interstellar radiation field and the 15,000 K blackbody radiation field predict the rotational temperature similar to the observed one.

preprint2014arXiv

Follow-up observations toward Planck cold clumps with ground-based radio telescopes

The physical and chemical properties of prestellar cores, especially massive ones, are still far from being well understood due to the lack of a large sample. The low dust temperature ($<$14 K) of Planck cold clumps makes them promising candidates for prestellar objects or for sources at the very initial stages of protostellar collapse. We have been conducting a series of observations toward Planck cold clumps (PCCs) with ground-based radio telescopes. In general, when compared with other star forming samples (e.g. infrared dark clouds), PCCs are more quiescent, suggesting that most of them may be in the earliest phase of star formation. However, some PCCs are associated with protostars and molecular outflows, indicating that not all PCCs are in a prestellar phase. We have identified hundreds of starless dense clumps from the mapping survey with the Purple Mountain Observatory (PMO) 13.7-m telescope. Follow-up observations suggest that these dense clumps are ideal targets to search for prestellar objects.

preprint2014arXiv

Herschel Key Program, "Dust, Ice, and Gas In Time" (DIGIT): the origin of molecular and atomic emission in low-mass protostars in Taurus

Six low-mass embedded sources (L1489, L1551-IRS5, TMR1, TMC1-A, L1527, and TMC1) in Taurus have been observed with Herschel-PACS to cover the full spectrum from 50 to 210 $μ$m as part of the Herschel key program, "Dust, Ice, and Gas In Time (DIGIT)". The relatively low intensity of the interstellar radiation field surrounding Taurus minimizes contamination of the [C II] emission associated with the sources by diffuse emission from the cloud surface, allowing study of the [C II] emission from the source. In several sources, the [C II] emission is distributed along the outflow, as is the [O I] emission. The atomic line luminosities correlate well with each other, as do the molecular lines, but the atomic and molecular lines correlate poorly. The relative contribution of CO to the total gas cooling is constant at $\sim$30 %, while the cooling fraction by H$_2$O varies from source to source, suggesting different shock properties resulting in different photodissociation levels of H$_2$O. The gas with a power-law temperature distribution with a moderately high density can reproduce the observed CO fluxes, indicative of CO close to LTE. However, H$_2$O is mostly subthermally excited. L1551-IRS5 is the most luminous source ($L_{\rm bol}=24.5$ L$_\odot$) and the [O I] 63.1 $μ$m line accounts for more than 70 % of its FIR line luminosity, suggesting complete photodissociation of H$_2$O by a J-shock. In L1551-IRS5, the central velocity shifts of the [O I] line, which exceed the wavelength calibration uncertainty ($\sim$70 km s$^{-1}$) of PACS, are consistent with the known red- and blue-shifted outflow direction.

preprint2014arXiv

Molecular Outflows Driven by Low-Mass Protostars. I. Correcting for Underestimates When Measuring Outflow Masses and Dynamical Properties

We present a survey of 28 molecular outflows driven by low-mass protostars, all of which are sufficiently isolated spatially and/or kinematically to fully separate into individual outflows. Using a combination of new and archival data from several single-dish telescopes, 17 outflows are mapped in CO (2-1) and 17 are mapped in CO (3-2), with 6 mapped in both transitions. For each outflow, we calculate and tabulate the mass, momentum, kinetic energy, mechanical luminosity, and force assuming optically thin emission in LTE at an excitation temperature of 50 K. We show that all of the calculated properties are underestimated when calculated under these assumptions. Taken together, the effects of opacity, outflow emission at low velocities confused with ambient cloud emission, and emission below the sensitivities of the observations increase outflow masses and dynamical properties by an order of magnitude, on average, and factors of 50-90 in the most extreme cases. Different (and non-uniform) excitation temperatures, inclination effects, and dissociation of molecular gas will all work to further increase outflow properties. Molecular outflows are thus almost certainly more massive and energetic than commonly reported. Additionally, outflow properties are lower, on average, by almost an order of magnitude when calculated from the CO (3-2) maps compared to the CO (2-1) maps, even after accounting for different opacities, map sensitivities, and possible excitation temperature variations. It has recently been argued in the literature that the CO (3-2) line is subthermally excited in outflows, and our results support this finding.

preprint2014arXiv

Radio Variability Survey of Very Low Luminosity Protostars

Ten very low luminosity objects were observed multiple times in the 8.5 GHz continuum in search of protostellar magnetic activities. A radio outburst of IRAM 04191+1522 IRS was detected, and the variability timescale was about 20 days or shorter. The results of this survey and archival observations suggest that IRAM 04191+1522 IRS is in active states about half the time. Archival data show that L1014 IRS and L1148 IRS were detectable previously and suggest that at least 20%-30% of very low luminosity protostars are radio variables. Considering the variability timescale and flux level of IRAM 04191+1522 IRS and the previous detection of the circular polarization of L1014 IRS, the radio outbursts of these protostars are probably caused by magnetic flares. However, IRAM 04191+1522 IRS is too young and small to develop an internal convective dynamo. If the detected radio emission is indeed coming from magnetic flares, the discovery implies that the flares may be caused by the fossil magnetic fields of interstellar origin.

preprint2014arXiv

Systematic Effects on the Genus Topology of Large Scale Structure of the Universe

Large-scale structure of the universe is a useful cosmological probe of the primordial non-Gaussianity and the expansion history of the universe because its topology does not change with time in the linear regime in the standard paradigm of structure formation. However, when the topology of iso-density contour surfaces is measured from an observational data, many systematic effects are introduced due to the finite size of pixels used to define the density field, non-linear gravitational evolution, redshift-space distortion, shot noise (discrete sampling), and bias in the distribution of the density field tracers. We study the various systematic effects on the genus curve to a great accuracy by using the Horizon Run 2 simulation of a ΛCDM cosmology. We numerically measure the genus curve from the gravitationally evolved matter and dark matter halo density fields. It is found that all the non-Gaussian deviations due to the systematic effects can be modeled by using a few low-order Hermite polynomials from H0 to H4. We compare our results with the analytic theories whenever possible, and find many new terms in the Hermite series that are making significant contributions to the non-Gaussian deviations. In particular, it is found that the amplitude drop of the genus curve due to the non-linear gravitational evolution can be accurately modeled by two terms H0 and H2 with coefficients both proportional to σ_0^2, the mean-square density fluctuation.

preprint2014arXiv

The D/H Ratio of Water Ice at Low Temperatures

We present the modeling results of deuterium fractionation of water ice, H2, and the primary deuterium isotopologues of H3+ adopting physical conditions associated with the star and planet formation process. We calculated the deuterium chemistry for a range of gas temperatures (T_gas ~ 10 - 30 K), molecular hydrogen density (n(H2)~ 10^4 - 10^7), and ortho/para ratio (opr) of H2 based on state-to-state reaction rates and explore the resulting fractionation including the formation of a water ice mantle coating grain surfaces. We find that the deuterium fractionation exhibits the expected temperature dependence of large enrichments at low gas temperature. More significantly the inclusion of water ice formation leads to large D/H ratios in water ice (>= 10^-2 at 10 K) but also alters the overall deuterium chemistry. For T < 20 K the implantation of deuterium into ices lowers the overall abundance of HD which reduces the efficiency of deuterium fractionation at high density. In agreement with an earlier study, under these conditions HD may not be the primary deuterium reservoir in the cold dense interstellar medium and H3+ will be the main charge carrier in the dense centers of pre-stellar cores and the protoplanetary disk midplane.

preprint2014arXiv

Thermal Starless Ammonia Core Surrounded by CCS in the Orion A Cloud

We imaged two starless molecular cloud cores, TUKH083 and TUKH122, in the Orion A giant molecular cloud in the CCS and ammonia (NH$_3$) emission with the Very Large Array. TUKH122 contains one NH$_3$ core "TUKH122-n," which is elongated and has a smooth oval boundary. Where observed, the CCS emission surrounds the NH$_3$ core. This configuration resembles that of the N$_2$H$^+$ and CCS distribution in the Taurus starless core L1544, a well-studied example of a dense prestellar core exhibiting infall motions. The linewidth of TUKH122-n is narrow (0.20 km s$^{-1}$) in the NH$_3$ emission line and therefore dominated by thermal motions. The smooth oval shape of the core boundary and narrow linewidth in NH$_3$ seem to imply that TUKH122-n is dynamically relaxed and quiescent. TUKH122-n is similar to L1544 in the kinetic temperature (10 K), linear size (0.03 pc), and virial mass ($\sim$ 2 $M_{\odot}$). Our results strongly suggest that TUKH122-n is on the verge of star formation. TUKH122-n is embedded in the 0.2 pc massive (virial mass $\sim$ 30 $M_{\odot}$) turbulent parent core, while the L1544 NH$_3$ core is embedded in the 0.2 pc less-massive (virial mass $\sim$ 10 $M_{\odot}$) thermal parent core. TUKH083 shows complicated distribution in NH$_3$, but was not detected in CCS. The CCS emission toward TUKH083 appears to be extended, and is resolved out in our interferometric observations.

preprint2013arXiv

An Analysis of the Environments of FU Orionis Objects with Herschel

We present Herschel-HIFI, SPIRE, and PACS 50-670 μm imaging and spectroscopy of six FU Orionis-type objects and candidates (FU Orionis, V1735 Cyg, V1515 Cyg, V1057 Cyg, V1331 Cyg, and HBC 722), ranging in outburst date from 1936-2010, from the "FOOSH" (FU Orionis Objects Surveyed with Herschel) program, as well as ancillary results from Spitzer-IRS and the Caltech Submillimeter Observatory. In their system properties (Lbol, Tbol, line emission), we find that FUors are in a variety of evolutionary states. Additionally, some FUors have features of both Class I and II sources: warm continuum consistent with Class II sources, but rotational line emission typical of Class I, far higher than Class II sources of similar mass/luminosity. Combining several classification techniques, we find an evolutionary sequence consistent with previous mid-IR indicators. We detect [O I] in every source at luminosities consistent with Class 0/I protostars, much greater than in Class II disks. We detect transitions of 13CO (J_up of 5 to 8) around two sources (V1735 Cyg and HBC 722) but attribute them to nearby protostars. Of the remaining sources, three (FU Ori, V1515 Cyg, and V1331 Cyg) exhibit only low-lying CO, but one (V1057 Cyg) shows CO up to J = 23 - 22 and evidence for H2O and OH emission, at strengths typical of protostars rather than T Tauri stars. Rotational temperatures for "cool" CO components range from 20-81 K, for ~ 10^50 total CO molecules. We detect [C I] and [N II] primarily as diffuse emission.

preprint2013arXiv

Chemical Variation in Molecular Cloud Cores in the Orion A Cloud. II

We have mapped six molecular cloud cores in the Orion A giant molecular cloud (GMC), whose kinetic temperatures range from 10 to 30 K, in CCS and N2H+ with Nobeyama 45 m radio telescope to study their chemical characteristics. We identified 31 intensity peaks in the CCS and N2H+ emission in these molecular cloud cores. It is found for cores with temperatures lower than ~ 25 K that the column density ratio of N(N2H+)/N(CCS) is low toward starless core regions while it is high toward star-forming core regions, in case that we detected both of the CCS and N2H+ emission. This is very similar to the tendency found in dark clouds (kinetic temperature ~ 10 K). The criterion found in the Orion A GMC is N(N2H+)/N(CCS) ~ 2-3. In some cases, the CCS emission is detected toward protostars as well as the N2H+ emission. Secondary late-stage CCS peak in the chemical evolution caused by CO depletion may be a possible explanation for this. We found that the chemical variation of CCS and N2H+ can also be used as a tracer of evolution in warm (10-25 K) GMC cores. On the other hand, some protostars do not accompany N2H+ intensity peaks but are associated with dust continuum emitting regions, suggesting that the N2H+ abundance might be decreased due to CO evaporation in warmer star-forming sites.

preprint2013arXiv

Embedded Protostars in the Dust, Ice, and Gas In Time (DIGIT) Key Program: Continuum SEDs, and an Inventory of Characteristic Far-Infrared Lines from PACS Spectroscopy

We present 50-210 um spectral scans of 30 Class 0/I protostellar sources, obtained with Herschel-PACS, and 0.5-1000 um SEDs, as part of the Dust, Ice, and Gas in Time (DIGIT) Key Program. Some sources exhibit up to 75 H2O lines ranging in excitation energy from 100-2000 K, 12 transitions of OH, and CO rotational lines ranging from J=14-13 up to J=40-39. [O I] is detected in all but one source in the entire sample; among the sources with detectable [O I] are two Very Low Luminosity Objects (VeLLOs). The mean 63/145 um [O I] flux ratio is 17.2 +/- 9.2. The [O I] 63 um line correlates with Lbol, but not with the time-averaged outflow rate derived from low-J CO maps. [C II] emission is in general not local to the source. The sample Lbol increased by 1.25 (1.06) and Tbol decreased to 0.96 (0.96) of mean (median) values with the inclusion of the Herschel data. Most CO rotational diagrams are characterized by two optically thin components (<N> = (0.70 +/- 1.12) x 10^49 total particles). N_CO correlates strongly with Lbol, but neither Trot nor N_CO(warm)/N_CO(hot) correlates with Lbol, suggesting that the total excited gas is related to the current source luminosity, but that the excitation is primarily determined by the physics of the interaction (e.g., UV- heating/shocks). Rotational temperatures for H2O (<Trot> = 194 +/- 85 K) and OH (<Trot> = 183 +/- 117 K) are generally lower than for CO, and much of the scatter in the observations about the best fit is attributed to differences in excitation conditions and optical depths amongst the detected lines.

preprint2013arXiv

L1448-MM observations by the Herschel Key program, "Dust, Ice, and Gas In Time" (DIGIT)

We present Herschel/PACS observations of L1448-MM, a Class 0 protostar with a prominent outflow. Numerous emission lines are detected at 55 < lambda < 210 micrometer including CO, OH, H2O, and [O I]. We investigate the spatial distribution of each transition to find that lines from low energy levels tend to distribute along the outflow direction while lines from high energy levels peak at the central spatial pixel. Spatial maps reveal that OH emission lines are formed in a relatively small area, while [O I] emission is extended. According to the rotational diagram analysis, the CO emission can be fitted by two (warm and hot) temperature components. For H2O, the ortho-to-para ratio is close to 3. The non-LTE LVG calculations suggest that CO and H2O lines could instead be formed in a high kinetic temperature (T > 1000 K) environment, indicative of a shock origin. For OH, IR-pumping processes play an important role in the level population. The molecular emission in L1448-MM is better explained with a C-shock model, but the atomic emission of PACS [O I] and Spitzer/IRS [Si II] emission is not consistent with C-shocks, suggesting multiple shocks in this region. Water is the major line coolant of L1448-MM in the PACS wavelength range, and the best-fit LVG models predict that H2O and CO emit (50-80)% of their line luminosity in the PACS wavelength range.

preprint2013arXiv

Resolved Depletion Zones and Spatial Differentiation of N2H+ and N2D+

We present a study on the spatial distribution of N2D+ and N2H+ in thirteen protostellar systems. Eight of thirteen objects observed with the IRAM 30m telescope show relative offsets between the peak N2D+ (J=2-1) and N2H+ (J=1-0) emission. We highlight the case of L1157 using interferometric observations from the Submillimeter Array and Plateau de Bure Interferometer of the N2D+ (J=3-2) and N2H+ (J=1-0) transitions respectively. Depletion of N2D+ in L1157 is clearly observed inside a radius of ~2000 AU (7") and the N2H+ emission is resolved into two peaks at radii of ~1000 AU (3.5"), inside the depletion region of N2D+. Chemical models predict a depletion zone in N2H+ and N2D+ due to destruction of H2D+ at T ~ 20 K and the evaporation of CO off dust grains at the same temperature. However, the abundance offsets of 1000 AU between the two species are not reproduced by chemical models, including a model that follows the infall of the protostellar envelope. The average abundance ratios of N2D+ to N2H+ have been shown to decrease as protostars evolve by Emprechtinger et al., but this is the first time depletion zones of N2D+ have been spatially resolved. We suggest that the difference in depletion zone radii for N2H+ and N2D+ is caused by either the CO evaporation temperature being above 20 K or an H2 ortho-to-para ratio gradient in the inner envelope.

preprint2013arXiv

Topology of Luminous Red Galaxies from the Sloan Digital Sky Survey

We present measurements of the genus topology of luminous red galaxies (LRGs) from the Sloan Digital Sky Survey (SDSS) Data Release 7 catalog, with unprecedented statistical significance. To estimate the uncertainties in the measured genus, we construct 81 mock SDSS LRG surveys along the past light cone from the Horizon Run 3, one of the largest N-body simulations to date that evolved 7210^3 particles in a 10815 Mpc/h size box. After carefully modeling and removing all known systematic effects due to finite pixel size, survey boundary, radial and angular selection functions, shot noise and galaxy biasing, we find the observed genus amplitude to reach 272 at 22 Mpc/h smoothing scale with an uncertainty of 4.2%; the estimated error fully incorporates cosmic variance. This is the most accurate constraint of the genus amplitude to date, which significantly improves on our previous results. In particular, the shape of the genus curve agrees very well with the mean topology of the SDSS LRG mock surveys in the LCDM universe. However, comparison with simulations also shows small deviations of the observed genus curve from the theoretical expectation for Gaussian initial conditions. While these discrepancies are mainly driven by known systematic effects such as those of shot noise and redshift-space distortions, they do contain important cosmological information on the physical effects connected with galaxy formation, gravitational evolution and primordial non-Gaussianity. We address here the key role played by systematics on the genus curve, and show how to accurately correct for their effects to recover the topology of the underlying matter. In a forthcoming paper, we provide an interpretation of those deviations in the context of the local model of non-Gaussianity.

preprint2013arXiv

Water and Methanol Maser Survey of Protostars in the Orion Molecular Cloud Complex

The results of a maser survey toward ninety-nine protostars in the Orion molecular cloud complex are presented. The target sources are low-mass protostars identified from infrared observations. Single-dish observations were carried out in the water maser line at 22 GHz and the methanol class I maser lines at 44, 95, and 133 GHz. Most of the detected sources were mapped to determine the source positions. Five water maser sources were detected, and they are excited by HH 1-2 VLA 3, HH 1-2 VLA 1, L1641N MM1/3, NGC 2071 IRS 1/3, and an object in the OMC 3 region. The water masers showed significant variability in intensity and velocity with time scales of a month or shorter. Four methanol emission sources were detected, and those in the OMC 2 FIR 3/4 and L1641N MM1/3 regions are probably masers. The methanol emission from the other two sources in the NGC 2071 IRS 1-3 and V380 Ori NE regions are probably thermal. For the water masers, the number of detections per protostar in the survey region is about 2%, which suggests that the water masers of low-mass protostars are rarely detectable. The methanol class I maser of low-mass protostars is an even rarer phenomenon, with a detection rate much smaller than 1%.

preprint2013arXiv

Wilson-Bappu Effect: Extended to Surface Gravity

Wilson and Bappu found a tight correlation between the stellar absolute visual magnitude (Mv) and the width of the Ca II K emission line for late-type stars in 1957. Here, we revisit the Wilson-Bappu relationship (hereafter, WBR) to claim that WBR can be an excellent indicator of stellar surface gravity of late-type stars as well as a distance indicator. We have measured the width (W) of the Ca II K emission line in high resolution spectra of 125 late-type stars, which were obtained with Bohyunsan Optical Echelle Spectrograph (BOES) and adopted from the UVES archive. Based on our measurement of the emission line width (W), we have obtained a WBR of Mv = 33.76 - 18.08 logW. In order to extend the WBR to be a surface gravity indicator, the stellar atmospheric parameters such as effective temperature (Teff), surface gravity (logg), metallicity ([Fe/H]), and micro-turbulence ($ξ_{tur}$) have been derived from the self-consistent detailed analysis using the Kurucz stellar atmospheric model and the abundance analysis code, MOOG. Using these stellar parameters and logW, we found that logg = - 5.85 logW + 9.97 logTeff - 23.48 for late-type stars.

preprint2012arXiv

CO J=2-1 Observations toward the Supernova Remnant G54.1+0.3

We present 12CO J = 2-1 line observations of G54.1+0.3, a composite supernova remnant with a mid-infrared (MIR) loop surrounding the central pulsar wind nebula (PWN). We mapped an area of 12' x 9' around the PWN and its associated MIR loop. We confirm two velocity components that had been proposed to be possibly interacting with the PWN/MIR-loop; the +53 km/s cloud that appears in contact with the eastern boundary of the PWN and the +23 km/s cloud that has CO emission coincident with the MIR loop. We have not found a direct evidence for the interaction in either of these clouds. Instead, we detected an 5'-long arc-like cloud at +15-+23 km/s with a systematic velocity gradient of ~3 km/s/arcmin and broad-line emitting CO gas having widths (FWHM) of <7 km/s in the western interior of the supernova remnant. We discuss their association with the supernova remnant.

preprint2012arXiv

CO2 Ice toward Low-luminosity, Embedded Protostars: Evidence for Episodic Mass Accretion via Chemical History

We present Spitzer IRS spectroscopy of CO2 ice bending mode spectra at 15.2 micrometer toward 19 young stellar objects with luminosity lower than 1 Lsun (3 with luminosity lower than 0.1 Lsun). Ice on dust grain surfaces can encode the history of heating because pure CO2 ice forms only at elevated temperature, T > 20 K, and thus around protostars of higher luminosity. Current internal luminosities of YSOs with L < 1 Lsun do not provide the conditions needed to produce pure CO2 ice at radii where typical envelopes begin. The presence of detectable amounts of pure CO2 ice would signify a higher past luminosity. Many of the spectra require a contribution from a pure, crystalline CO2 component, traced by the presence of a characteristic band splitting in the 15.2 micrometer bending mode. About half of the sources (9 out of 19) in the low luminosity sample have evidence for pure CO2 ice, and six of these have significant double-peaked features, which are very strong evidence of pure CO2 ice. The presence of the pure CO2 ice component indicates that the dust temperature, and hence luminosity of the central star/accretion disk system, must have been higher in the past. An episodic accretion scenario, in which mixed CO-CO2 ice is converted to pure CO2 ice during each high luminosity phase, explains the presence of pure CO2 ice, the total amount of CO2 ice, and the observed residual C18O gas.

preprint2012arXiv

Radio Imaging of the NGC 2024 FIR 5/6 Region: a Hypercompact H II Region Candidate in Orion

The NGC 2024 FIR 5/6 region was observed in the 6.9 mm continuum with an angular resolution of about 1.5 arcsec. The 6.9 mm continuum map shows four compact sources, FIR 5w, 5e, 6c, and 6n, as well as an extended structure of the ionization front associated with the optical nebulosity. FIR 6c has a source size of about 0.4 arcsec or 150 AU. The spectral energy distribution (SED) of FIR 6c is peculiar: rising steeply around 6.9 mm and flat around 1 mm. The possibility of a hypercompact H II region is explored. If the millimeter flux of FIR 6c comes from hot ionized gas heated by a single object at the center, the central object may be a B1 star of about 5800 solar luminosities and about 13 solar masses. The 6.9 mm continuum of FIR 6n may be a mixture of free-free emission and dust continuum emission. Archival data show that both FIR 6n and 6c exhibit water maser activity, suggesting the existence of shocked gas around them. The 6.9 mm continuum emission from FIR 5w has a size of about 1.8 arcsec or 760 AU. The SEDs suggest that the 6.9 mm emission of FIR 5w and 5e comes from dust, and the masses of the dense molecular gas are about 0.6 and 0.5 solar masses, respectively.

preprint2012arXiv

Water and Methanol Maser Activities in the NGC 2024 FIR 6 Region

The NGC 2024 FIR 6 region was observed in the water maser line at 22 GHz and the methanol class I maser lines at 44, 95, and 133 GHz. The water maser spectra displayed several velocity components and month-scale time variabilities. Most of the velocity components may be associated with FIR 6n, while one component was associated with FIR 4. A typical life time of the water-maser velocity-components is about 8 months. The components showed velocity fluctuations with a typical drift rate of about 0.01 km/s/day. The methanol class I masers were detected toward FIR 6. The methanol emission is confined within a narrow range around the systemic velocity of the FIR 6 cloud core. The methanol masers suggest the existence of shocks driven by either the expanding H II region of FIR 6c or the outflow of FIR 6n.

preprint2011arXiv

Disentangling the Environment of the FU Orionis Candidate HBC 722 with Herschel

We analyze the submillimeter emission surrounding the new FU Orionis-type object, HBC 722. We present the first epoch of observations of the active environs of HBC 722, with imaging and spectroscopy from PACS, SPIRE, and HIFI aboard the Herschel Space Observatory, as well as CO J= 2-1 and 350 um imaging (SHARC-II) with the Caltech Submillimeter Observatory. The primary source of submillimeter continuum emission in the region -- 2MASS 20581767+4353310 -- is located 16$\arcsec$ south-southeast of the optical flaring source while the optical and near-IR emission is dominated by HBC 722. A bipolar outflow extends over HBC 722; the most likely driver is the submillimeter source. We detect warm (100 K) and hot (246 K) CO emission in the surrounding region, evidence of outflow-driven heating in the vicinity. The region around HBC 722 itself shows little evidence of heating driven by the new outbursting source itself.

preprint2011arXiv

Medium Resolution Spectral Library of Late-Type Stellar Templates in Near-Infrared Band

We present medium resolution (R = 5000-6000) spectra in the near-infrared band, 1.4-1.8 microns, for template stars in G, K, and M types observed by the echelle spectrometer, IRCS, at the SUBARU 8.2 m telescope. The identification of lines is based on the spectra of Arcturus (K2 III) in the literature. We measured the equivalent of widths and compared our results to those of Meyer et al. (1998). We conclude that our spectral resolution (R = 6000) data can investigate more accurately the properties of lines in stellar spectra. The library of the template stellar spectra in ASCII format are available for download on the World Wide Web.

preprint2011arXiv

Radio Imaging of the NGC 1333 IRAS 4A Region: Envelope, Disks, and Outflows of a Protostellar Binary System

The NGC 1333 IRAS 4A protobinary was observed in the 1.3 cm and 6.9 mm continuum and the ammonia and SiO lines, with an angular resolution of about 0.4 arcseconds. The continuum maps show the circumstellar structures of the two protostars, A1 and A2. The A1 system is brighter and more massive than the A2 system. The ratio of mass, including dense gas and protostar, is about 6. The properties of the circumstellar disks and outflows suggest that A1 may be younger than A2. The deflected part of the northeastern jet of A2 is bright in the SiO line, and the distance between the brightest peak and deflection point suggests that the enhancement of SiO takes about 100 yr after the collision with a dense core. The ammonia maps show a small structure that seems to be a part of the obstructing core. The outflow properties were studied by comparing interferometric maps of SiO, ammonia, formaldehyde, and HCN lines. Their overall structures agree well, suggesting that these species are excited by the same mechanism. However, the intensity distributions show that SiO is chemically unique. SiO may be directly linked to the primary jet while the other species may be tracing the entrained ambient gas.

preprint2011arXiv

Radio Imaging of the NGC 1333 IRAS 4B Region

The NGC 1333 IRAS 4B region was observed in the 6.9 mm and 1.3 cm continuum with an angular resolution of about 0.4 arcseconds. IRAS 4BI was detected in both bands, and BII was detected in the 6.9 mm continuum only. The 1.3 cm source of BI seems to be a disk-like flattened structure with a size of about 50 AU. IRAS 4BI does not show any sign of multiplicity. Examinations of archival infrared images show that the dominating emission feature in this region is a bright peak in the southern outflow driven by BI, corresponding to the molecular hydrogen emission source HL 9a. Both BI and BII are undetectable in the mid-IR bands. The upper limit on the far-IR flux of IRAS 4BII suggests that it may be a very low luminosity young stellar object.

preprint2011arXiv

Submillimeter Observations of Dense Clumps in the Infrared Dark Cloud G049.40-00.01

We obtained 350 and 850 micron continuum maps of the infrared dark cloud G049.40-00.01. Twenty-one dense clumps were identified within G049.40-00.01 based on the 350 micron continuum map with an angular resolution of about 9.6". We present submillimeter continuum maps and report physical properties of the clumps. The masses of clumps range from 50 to 600 M_sun. About 70% of the clumps are associated with bright 24 micron emission sources, and they may contain protostars. The most massive two clumps show extended, enhanced 4.5 micron emission indicating vigorous star-forming activity. The clump size-mass distribution suggests that many of them are forming high mass stars. G049.40-00.01 contains numerous objects in various evolutionary stages of star formation, from pre-protostellar clumps to HII regions.

preprint2011arXiv

The Spitzer c2d Survey of Nearby Dense Cores XI : Infrared and Submillimeter Observations of CB130

We present new observations of the CB130 region, composed of three separate cores. Using the \textit{Spitzer Space Telescope} we detected a Class 0 and a Class II object in one of these, CB130-1. The observed photometric data from \textit{Spitzer} and ground-based telescopes are used to establish the physical parameters of the Class 0 object. SED fitting with a radiative transfer model shows that the luminosity of the Class 0 object is 0.14 $-$ 0.16 L$_{\odot}$, which is a low luminosity for a protostellar object. In order to constrain the chemical characteristics of the core having the low luminosity object, we compare our molecular line observations to models of lines including abundance variations. We tested both ad hoc step function abundance models and a series of self-consistent chemical evolution models. In the chemical evolution models, we consider a continuous accretion model and an episodic accretion model to explore how variable luminosity affects the chemistry. The step function abundance models can match observed lines reasonably well. The best fitting chemical evolution model requires episodic accretion and the formation of CO$_2$ ice from CO ice during the low luminosity periods. This process removes C from the gas phase, providing a much improved fit to the observed gas-phase molecular lines and the CO$_2$ ice absorption feature. Based on the chemical model result, the low luminosity of CB130-1 is explained better as a quiescent stage between episodic accretion bursts rather than being at the first hydrostatic core stage.

preprint2010arXiv

A Candidate Detection of the First Hydrostatic Core

The first hydrostatic core (FHSC) represents a very early phase in the low-mass star formation process, after collapse of the parent core has begun but before a true protostar has formed. This large (few AU), cool (100 K), pressure supported core of molecular hydrogen is expected from theory, but has yet to be observationally verified. Here we present observations of an excellent candidate for the FHSC phase: Per-Bolo 58, a dense core in Perseus that was previously believed to be starless. The 70 micron flux of 65 mJy, from new deep Spitzer MIPS observations, is consistent with that expected for the FHSC. A low signal-to-noise detection at 24 micron leaves open the possibility that Per-Bolo 58 could be a very low luminosity protostar, however. We utilize radiative transfer models to determine the best-fitting FHSC and protostar models to the spectral energy distribution and 2.9 mm visibilities of Per-Bolo 58. The source is consistent with a FHSC with some source of lower opacity through the envelope allowing 24 micron emission to escape; a small outflow cavity or a cavity in the envelope are both possible. While we are unable to rule out the presence of a protostar, if present it would be one of the lowest luminosity protostellar objects yet observed, with an internal luminosity of approximately 0.01 Lsun.

preprint2010arXiv

Effect of FIR Fluxes on Constraining Properties of YSOs

Young Stellar Objects (YSOs) in the early evolutionary stages are very embedded, and thus they emit most of their energy at long wavelengths such as far-infrared (FIR) and submillimeter (Submm). Therefore, the FIR observational data are very important to classify the accurate evolutionary stages of these embedded YSOs, and to better constrain their physical parameters in the dust continuum modeling. We selected 28 YSOs, which were detected in the AKARI Far-Infrared Surveyor (FIS), from the Spitzer c2d legacy YSO catalogs to test the effect of FIR fluxes on the classification of their evolutionary stages and on the constraining of envelope properties, internal luminosity, and UV strength of the Interstellar Radiation Field (ISRF). According to our test, one can mis-classify the evolutionary stages of YSOs, especially the very embedded ones if the FIR fluxes are not included. In addition, the total amount of heating of YSOs can be underestimated without the FIR observational data.

preprint2010arXiv

The Solar Nebula on Fire: A Solution to the Carbon Deficit in the Inner Solar System

Despite a surface dominated by carbon-based life, the bulk composition of the Earth is dramatically carbon poor when compared to the material available at formation. Bulk carbon deficiency extends into the asteroid belt representing a fossil record of the conditions under which planets are born. The initial steps of planet formation involve the growth of primitive sub-micron silicate and carbon grains in the Solar Nebula. We present a solution wherein primordial carbon grains are preferentially destroyed by oxygen atoms ignited by heating due to stellar accretion at radii < 5 AU. This solution can account for the bulk carbon deficiency in the Earth and meteorites, the compositional gradient within the asteroid belt, and for growing evidence for similar carbon deficiency in rocks surrounding other stars.

preprint2009arXiv

The Spitzer c2d Survey of Nearby Dense Cores VII: Chemistry and Dynamics in L43

We present results from the Spitzer Space Telescope and molecular line observations of 9 species toward the dark cloud L43. The Spitzer images and molecular line maps suggest it has a starless core and a Class I protostar evolving in the same environment. CO depletion is seen in both sources, and DCO+ lines are stronger toward the starless core. With a goal of testing the chemical characteristics from pre- to protostellar stages, we adopt an evolutionary chemical model to calculate the molecular abundances and compare with our observations. Among the different model parameters we tested, the best-fit model suggests a longer total timescale at the pre-protostellar stage, but with faster evolution at the later steps with higher densities.

preprint2009arXiv

The Spitzer c2d Survey of Nearby Dense Cores: Jet and Molecular Outflow Associated with a YSO in core A of L1251

A long infrared jet has been discovered by the Spitzer c2d legacy program in core A of L1251. It is associated with a very embedded Class 0 object with an accretion luminosity of about 0.9 Lsun derived by radiative transfer model fitting to the observed SED. Comparing the observed IRAC colors along the infrared jet with those calculated from a model of an admixture of gas with a power-law temperature distribution indicates that the jet is possibly created by a paraboloidal bow shock propagating into the ambient medium of n(H_2)=10^5 cm^{-3}. In addition, the variation of the power-law index along the jet suggests that the portion of hot gas decreases with distance from the jet engine. The molecular outflow in this region has been mapped for the first time using CO data. From the calculated outflow momentum flux, a very strong lower limit to the average accretion luminosity is 3.6 sin i/cos^3 i Lsun, indicative of a decrease in the accretion rate with time.

preprint2008arXiv

Centimeter Polarimetry of the R Coronae Australis Region

Circularly polarized 3.5 cm continuum emission was detected toward three radio sources in the R CrA region using the Very Large Array. The Class I protostar IRS 5b persistently showed polarized radio emission with a constant helicity over 8 yr, which suggests that its magnetosphere has a stable configuration. There is a good correlation between the Stokes I and Stokes V fluxes, and the fractional polarization is about 0.17. During active phases the fractional polarization is a weakly decreasing function of Stokes I flux, which suggests that IRS 5b is phenomenologically similar to other types of flare stars such as RS CVn binaries. The variability timescale of the polarized flux is about a month, and the magnetosphere of IRS 5b must be very large in size. The Class I protostar IRS 7A was detected once in circularly polarized radio emission, even though IRS 7A drives a thermal radio jet. This detection implies that the radio emission from the magnetosphere of a young protostar can escape the absorption by the partially ionized wind at least once in a while. The properties of IRS 7A and IRS 5b suggests that Class I protostars have organized peristellar magnetic fields of a few kilogauss and that the detectability of magnetospheric emission may depend on the evolutionary status of protostar. Also reported is the detection of circularly polarized radio emission toward the variable radio source B5.

preprint2008arXiv

Oxygen isotope anomalies of the Sun and the original environment of the Solar system

We present results from a model of oxygen isotopic anomaly production through selective photodissociation of CO within the collapsing proto-Solar cloud. Our model produces a proto-Sun with a wide range of Delta_17O values depending on the intensity of the ultraviolet radiation field. Dramatically different results from two recent Solar wind oxygen isotope measurements indicate that a variety of compositions remain possible for the solar oxygen isotope composition. However, constrained by other measurements from comets and meteorites, our models imply the birth of the Sun in a stellar cluster with an enhanced radiation field and are therefore consistent with a supernova source for 60Fe in meteorites.

preprint2006arXiv

The Spitzer c2d Survey of Nearby Dense Cores: III: Low Mass Star Formation in a Small Group, L1251B

We present a comprehensive study of a low-mass star-forming region,L1251B, at wavelengths from the near-infrared to the millimeter. L1251B, where only one protostar, IRAS 22376+7455, was known previously, is confirmed to be a small group of protostars based on observations with the Spitzer Space Telescope. The most luminous source of L1251B is located 5" north of the IRAS position. A near-infrared bipolar nebula, which is not associated with the brightest object and is located at the southeast corner of L1251B, has been detected in the IRAC bands. OVRO and SMA interferometric observations indicate that the brightest source and the bipolar nebula source in the IRAC bands are deeply embedded disk sources.Submillimeter continuum observations with single-dish telescopes and the SMA interferometric observations suggest two possible prestellar objects with very high column densities. Outside of the small group, many young stellar object candidates have been detected over a larger region of 12' x 12'. Extended emission to the east of L1251B has been detected at 850 micron; this "east core" may be a site for future star formation since no point source has been detected with IRAC or MIPS. This region is therefore a possible example of low-mass cluster formation, where a small group of pre- and protostellar objects (L1251B) is currently forming, alongside a large starless core (the east core).