Researcher profile

Jeong-Eun Lee

Jeong-Eun Lee contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
31works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

31 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&#34;. 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 &#34;ATOMS,&#34; 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 &#34;ATOMS&#34; 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., &#34;MALT90&#34;, &#34;Orion B&#34;, &#34;EMPIRE&#34;, &#34;ALMA\_IMF&#34;, and &#34;ALMAGAL&#34;, to deepen our understandings of &#34;dense gas&#34; 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 &#34;ATOMS&#34; survey of 146 active Galactic star forming regions. The correlations between $L_{\rm bol}$ and molecular line luminosities $L&#39;_{\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&#39;_{\rm mol}$ relations as their main line emission. The log($L_{\rm bol}$/$L&#39;_{\rm mol}$) for different molecular line tracers have similar distributions. The $L_{\rm bol}$-to-$L&#39;_{\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 &#34;ATOMS&#34; 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.

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 &#34;TUKH122-n,&#34; 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.