Showing 1–15 of 15 results for author: Seo, Y M

Retrievals of Protoplanetary Disk Parameters using Thermochemical Models: I. Disk Gas Mass from Hydrogen Deuteride Spectroscopy

Authors: Young Min Seo, Karen Willacy, Geoffrey Bryden, Dariusz C. Lis, Paul F. Goldsmith, Klaus M. Pontoppidan, Wing-Fai Thi

Abstract: We discuss statistical relationships between the mass of protoplanetary disks and the hydrogen deuteride (HD) line emission and the dust spectral energy distribution (SED) determined using 3000 ProDiMo disk models. The models have 15 free parameters describing disk physical properties, the central star, and the local radiation field. The sampling of physical parameters is done using a Monte Carlo… ▽ More We discuss statistical relationships between the mass of protoplanetary disks and the hydrogen deuteride (HD) line emission and the dust spectral energy distribution (SED) determined using 3000 ProDiMo disk models. The models have 15 free parameters describing disk physical properties, the central star, and the local radiation field. The sampling of physical parameters is done using a Monte Carlo approach to evaluate the probability density functions of observables as a function of physical parameters. We find that the HD fractional abundance is almost constant even though the UV flux varies by several orders of magnitude. Probing the statistical relation between the physical quantities and the HD flux, we find that low-mass (optically thin) disks display a tight correlation between the average disk gas temperature and HD line flux, while massive disks show no such correlation. We demonstrate that the central star luminosity, disk size, dust size distribution, and HD flux may be used to determine the disk gas mass to within a factor of three. We also find that the far-IR and sub-mm/mm SEDs and the HD flux may serve as strong constraints for determining the disk gas mass to within a factor of two. If the HD lines are fully spectrally resolved ($R\gtrsim 1.5\times10^6, Δv=0.2~\rm km\,s^{-1}$), the 56 $μ$m and 112 $μ$m HD line profiles alone may constrain the disk gas mass to within a factor of two. △ Less

Submitted 24 April, 2024; v1 submitted 24 October, 2023; originally announced October 2023.

Comments: 31 pages, 16 figures, ApJ accepted, preprint

Statistical Prediction of [CII] Observations by Constructing Probability Density Functions using SOFIA, Herschel, and Spitzer Observations

Authors: Young Min Seo, Karen Willacy, Umaa Rebbapragada

Abstract: We present a statistical algorithm for predicting the [CII] emission from Herschel and Spitzer continuum images using probability density functions between the [CII] emission and continuum emission. The [CII] emission at 158 $μ$m is a critical tracer in studying the life cycle of interstellar medium and galaxy evolution. Unfortunately, its frequency is in the far infrared (FIR), which is opaque th… ▽ More We present a statistical algorithm for predicting the [CII] emission from Herschel and Spitzer continuum images using probability density functions between the [CII] emission and continuum emission. The [CII] emission at 158 $μ$m is a critical tracer in studying the life cycle of interstellar medium and galaxy evolution. Unfortunately, its frequency is in the far infrared (FIR), which is opaque through the troposphere and cannot be observed from the ground except for highly red-shifted sources (z $\gtrsim$ 2). Typically [CII] observations of closer regions have been carried out using suborbital or space observatories. Given the high cost of these facilities and limited time availability, it is important to have highly efficient observations/operations in terms of maximizing science returns. This requires accurate prediction of the strength of emission lines and, therefore, the time required for their observation. However, [CII] emission has been hard to predict due to a lack of strong correlations with other observables. Here we adopt a new approach to making accurate predictions of [CII] emission by relating this emission simultaneously to several tracers of dust emission in the same region. This is done using a statistical methodology utilizing probability density functions (PDFs) among [CII] emission and Spitzer IRAC and Herschel PACS/SPIRE images. Our test result toward a star-forming region, RCW 120, demonstrates that our methodology delivers high-quality predictions with less than 30\% uncertainties over 80\% of the entire observation area, which is more than sufficient to test observation feasibility and maximize science return. The {\it pickle} dump files storing the PDFs and trained neural network module are accessible upon request and will support future far-infrared missions, for example, GUSTO and FIR Probe. △ Less

Submitted 21 August, 2023; originally announced August 2023.

Comments: 7 figures

Journal ref: Journal of Astronomical Instrumentation, 2023, 10.1142/S2251171723500071

Alignment of dense molecular core morphology and velocity gradients with ambient magnetic fields

Authors: A. Pandhi, R. K. Friesen, L. Fissel, J. E. Pineda, P. Caselli, M. C-Y. Chen, J. Di Francesco, A. Ginsburg, H. Kirk, P. C. Myers, S. S. R. Offner, A. Punanova, F. Quan, E. Redaelli, E. Rosolowsky, S. Scibelli, Y. M. Seo, Y. Shirley

Abstract: Studies of dense core morphologies and their orientations with respect to gas flows and the local magnetic field have been limited to only a small sample of cores with spectroscopic data. Leveraging the Green Bank Ammonia Survey alongside existing sub-millimeter continuum observations and Planck dust polarization, we produce a cross-matched catalogue of 399 dense cores with estimates of core morph… ▽ More Studies of dense core morphologies and their orientations with respect to gas flows and the local magnetic field have been limited to only a small sample of cores with spectroscopic data. Leveraging the Green Bank Ammonia Survey alongside existing sub-millimeter continuum observations and Planck dust polarization, we produce a cross-matched catalogue of 399 dense cores with estimates of core morphology, size, mass, specific angular momentum, and magnetic field orientation. Of the 399 cores, 329 exhibit 2D $\mathrm{v}_\mathrm{LSR}$ maps that are well fit with a linear gradient, consistent with rotation projected on the sky. We find a best-fit specific angular momentum and core size relationship of $J/M \propto R^{1.82 \pm 0.10}$, suggesting that core velocity gradients originate from a combination of solid body rotation and turbulent motions. Most cores have no preferred orientation between the axis of core elongation, velocity gradient direction, and the ambient magnetic field orientation, favouring a triaxial and weakly magnetized origin. We find, however, strong evidence for a preferred anti-alignment between the core elongation axis and magnetic field for protostellar cores, revealing a change in orientation from starless and prestellar populations that may result from gravitational contraction in a magnetically-regulated (but not dominant) environment. We also find marginal evidence for anti-alignment between the core velocity gradient and magnetic field orientation in the L1228 and L1251 regions of Cepheus, suggesting a preferred orientation with respect to magnetic fields may be more prevalent in regions with locally ordered fields. △ Less

Submitted 24 July, 2023; originally announced July 2023.

Comments: 33 pages, 28 figures, accepted to MNRAS

Probing Polarization and the Role of Magnetic Fields in Cloud Destruction in the Keyhole Nebula

Authors: Young Min Seo, C. Darren Dowell, Paul F. Goldsmith, Jorge L. Pineda, Liton Majumdar

Abstract: We present polarimetric observations of the Keyhole Nebula in the Carina Nebula Complex carried out using the Stratospheric Observatory for Infrared Astronomy. The Keyhole Nebula located to the west of $η$ Carinae is believed to be disturbed by the stellar winds from the star. We observed the Keyhole Nebula at 89 $μ$m wavelength with the HAWC+ instrument. The observations cover the entire Keyhole… ▽ More We present polarimetric observations of the Keyhole Nebula in the Carina Nebula Complex carried out using the Stratospheric Observatory for Infrared Astronomy. The Keyhole Nebula located to the west of $η$ Carinae is believed to be disturbed by the stellar winds from the star. We observed the Keyhole Nebula at 89 $μ$m wavelength with the HAWC+ instrument. The observations cover the entire Keyhole Nebula spanning 8$'$ by 5$'$ with central position RA = 10:44:43 and Dec = -59:38:04. The typical uncertainty of polarization measurement is less than 0.5\% in the region with intensity above 5,500 MJy sr$^{-1}$. The polarization has a mean of 2.4\% with a standard deviation of 1.6\% in the region above this intensity, similar to values in other high--mass star--forming regions. The magnetic field orientation in the bar--shaped structure is similar to the large--scale magnetic field orientation. On the other hand, the magnetic field direction in the loop is not aligned with the large--scale magnetic fields but has tight alignment with the loop itself. Analysis of the magnetic field angles and the gas turbulence suggests that the field strength is $\sim$70 $μ$G in the loop. A simple comparison of the magnetic field tension to the ram pressure of $η$ Carinae's stellar wind suggests that the magnetic fields in the Keyhole Nebula are not strong enough to maintain the current structure against the impact of the stellar wind, and that the role of the magnetic field in resisting stellar feedback in the Keyhole Nebula is limited. △ Less

Submitted 21 June, 2021; originally announced June 2021.

Comments: ApJ accepted. This is the version of the article before peer review or editing, as submitted by an author to ApJ. IOP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it

Applications of Machine Learning Algorithms In Processing Terahertz Spectroscopic Data

Authors: Young Min Seo, Paul F. Goldsmith, Volker Tolls, Russell Shipman, Craig Kulesa, William Peters, Christopher Walker, Gary Melnick

Abstract: We present the data reduction software and the distribution of Level 1 and Level 2 products of the Stratospheric Terahertz Observatory 2 (STO2). STO2, a balloon-borne Terahertz telescope, surveyed star-forming regions and the Galactic plane and produced approximately 300,000 spectra. The data are largely similar to spectra typically produced by single-dish radio telescopes. However, a fraction of… ▽ More We present the data reduction software and the distribution of Level 1 and Level 2 products of the Stratospheric Terahertz Observatory 2 (STO2). STO2, a balloon-borne Terahertz telescope, surveyed star-forming regions and the Galactic plane and produced approximately 300,000 spectra. The data are largely similar to spectra typically produced by single-dish radio telescopes. However, a fraction of the data contained rapidly varying fringe/baseline features and drift noise, which could not be adequately corrected using conventional data reduction software. To process the entire science data of the STO2 mission, we have adopted a new method to find proper off-source spectra to reduce large-amplitude fringes and new algorithms including Asymmetric Least Square (ALS), Independent Component Analysis (ICA), and Density-based spatial clustering of applications with noise (DBSCAN). The STO2 data reduction software efficiently reduced the amplitude of fringes from a few hundred to 10 K and resulted in baselines of amplitude down to a few K. The Level 1 products typically have the noise of a few K in [CII] spectra and ~1 K in [NII] spectra. Using a regridding algorithm, we made spectral maps of star-forming regions and the Galactic plane survey using an algorithm employing a Bessel-Gaussian kernel. Level 1 and 2 products are available to the astronomical community through the STO2 data server and the DataVerse. The software is also accessible to the public through Github. The detailed addresses are given in Section 4 of the paper on data distribution. △ Less

Submitted 2 September, 2020; originally announced September 2020.

Comments: Accepted to Journal of Astronomical Instrumentation

Probing ISM Structure in Trumpler 14 & Carina I Using The Stratospheric Terahertz Observatory 2

Authors: Young Min Seo, Paul F. Goldsmith, Chris Walker, David J. Hollenbach, Mark G. Wolfire, Craig Kulesa, Volker Tolls, Pietro N. Bernasconi, Umit Kavak, Floris F. S. van der Tak, Russ Shipman, Jian Rong Gao, Alexander Tielens, Michael G. Burton, Harold Yorke, Erick Young, William L. Peters, Abram Young, Christopher Groppi, Kristina Davis, Jorge L. Pineda, William D. Langer, Jonathan H. Kawamura, Antony Stark, Gary Melnick , et al. (4 additional authors not shown)

Abstract: We present observations of the Trumpler 14/Carina I region carried out using the Stratospheric Terahertz Observatory 2 (STO2). The Trumpler 14/Carina I region is in the west part of the Carina Nebula Complex, which is one of the most extreme star-forming regions in the Milky Way. We observed Trumpler 14/Carina I in the 158 $μ$m transition of [C\,{\sc ii}] with a spatial resolution of 48$''$ and a… ▽ More We present observations of the Trumpler 14/Carina I region carried out using the Stratospheric Terahertz Observatory 2 (STO2). The Trumpler 14/Carina I region is in the west part of the Carina Nebula Complex, which is one of the most extreme star-forming regions in the Milky Way. We observed Trumpler 14/Carina I in the 158 $μ$m transition of [C\,{\sc ii}] with a spatial resolution of 48$''$ and a velocity resolution of 0.17 km s$^{-1}$. The observations cover a 0.25$^\circ$ by 0.28$^\circ$ area with central position {\it l} = 297.34$^\circ$, {\it b} = -0.60$^\circ$. The kinematics show that bright [C\,{\sc ii}] structures are spatially and spectrally correlated with the surfaces of CO clouds, tracing the photodissociation region and ionization front of each molecular cloud. Along 7 lines of sight that traverse Tr 14 into the dark ridge to the southwest, we find that the [C\,{\sc ii}] luminosity from the HII region is 3.7 times that from the PDR. In same los we find in the PDRs an average ratio of 1:4.1:5.6 for the mass in atomic gas:dark-CO gas: molecular gas traced by CO. Comparing multiple gas tracers including HI 21cm, [C\,{\sc ii}], CO, and radio recombination lines, we find that the HII regions of the Carina Nebula Complex are well-described as HII regions with one-side freely expanding towards us, consistent with the champagne model of ionized gas evolution. The dispersal of the GMC in this region is dominated by EUV photoevaporation; the dispersal timescale is 20-30 Myr. △ Less

Submitted 22 May, 2019; v1 submitted 22 March, 2019; originally announced March 2019.

Comments: ApJ accepted

An Ammonia Spectral Map of the L1495-B218 Filaments in the Taurus Molecular Cloud: II CCS & HC$_7$N Chemistry and Three Modes of Star Formation in the Filaments

Authors: Young Min Seo, Liton Majumdar, Paul F. Goldsmith, Yancy L. Shirley, Karen Willacy, Derek Ward-Thompson, Rachel Friesen, David Frayer, Sarah E. Church, Dongwoo Chung, Kieran Cleary, Nichol Cunningham, Kiruthika Devaraj, Dennis Egan, Todd Gaier, Rohit Gawande, Joshua O. Gundersen, Andrew I. Harris, Pekka Kangaslahti, Anthony C. S. Readhead, Lorene Samoska, Matthew Sieth, Michael Stennes, Patricia Voll, Steve White

Abstract: We present deep CCS and HC$_7$N observations of the L1495-B218 filaments in the Taurus molecular cloud obtained using the K-band focal plane array on the 100m Green Bank Telescope. We observed the L1495-B218 filaments in CCS $J_N$ = 2$_1$$-$1$_0$ and HC$_7$N $J$ = 21$-$20 with a spectral resolution of 0.038 km s$^{-1}$ and an angular resolution of 31$''$. We observed strong CCS emission in both ev… ▽ More We present deep CCS and HC$_7$N observations of the L1495-B218 filaments in the Taurus molecular cloud obtained using the K-band focal plane array on the 100m Green Bank Telescope. We observed the L1495-B218 filaments in CCS $J_N$ = 2$_1$$-$1$_0$ and HC$_7$N $J$ = 21$-$20 with a spectral resolution of 0.038 km s$^{-1}$ and an angular resolution of 31$''$. We observed strong CCS emission in both evolved and young regions and weak emission in two evolved regions. HC$_7$N emission is observed only in L1495A-N and L1521D. We find that CCS and HC$_7$N intensity peaks do not coincide with NH$_3$ or dust continuum intensity peaks. We also find that the fractional abundance of CCS does not show a clear correlation with the dynamical evolutionary stage of dense cores. Our findings and chemical modeling indicate that the fractional abundances of CCS and HC$_7$N are sensitive to the initial gas-phase C/O ratio, and they are good tracers of young condensed gas only when the initial C/O is close to solar value. Kinematic analysis using multiple lines including NH$_3$, HC$_7$N, CCS, CO, HCN, \& HCO$^+$ suggests that there may be three different star formation modes in the L1495-B218 filaments. At the hub of the filaments, L1495A/B7N has formed a stellar cluster with large-scale inward flows (fast mode), while L1521D, a core embedded in a filament, is slowly contracting due to its self-gravity (slow mode). There is also one isolated core that appears to be marginally stable and may undergo quasi-static evolution (isolated mode). △ Less

Submitted 14 December, 2018; originally announced December 2018.

Comments: ApJ accepted

Droplets I: Pressure-Dominated Sub-0.1 pc Coherent Structures in L1688 and B18

Authors: Hope How-Huan Chen, Jaime E. Pineda, Alyssa A. Goodman, Andreas Burkert, Stella S. R. Offner, Rachel K. Friesen, Philip C. Myers, Felipe Alves, Hector G. Arce, Paola Caselli, Ana Chacon-Tanarro, Michael Chun-Yuan Chen, James Di Francesco, Adam Ginsburg, Jared Keown, Helen Kirk, Peter G. Martin, Christopher Matzner, Anna Punanova, Elena Redaelli, Erik Rosolowsky, Samantha Scibelli, Young Min Seo, Yancy Shirley, Ayushi Singh

Abstract: We present the observation and analysis of newly discovered coherent structures in the L1688 region of Ophiuchus and the B18 region of Taurus. Using data from the Green Bank Ammonia Survey (GAS), we identify regions of high density and near-constant, almost-thermal, velocity dispersion. Eighteen coherent structures are revealed, twelve in L1688 and six in B18, each of which shows a sharp "transiti… ▽ More We present the observation and analysis of newly discovered coherent structures in the L1688 region of Ophiuchus and the B18 region of Taurus. Using data from the Green Bank Ammonia Survey (GAS), we identify regions of high density and near-constant, almost-thermal, velocity dispersion. Eighteen coherent structures are revealed, twelve in L1688 and six in B18, each of which shows a sharp "transition to coherence" in velocity dispersion around its periphery. The identification of these structures provides a chance to study the coherent structures in molecular clouds statistically. The identified coherent structures have a typical radius of 0.04 pc and a typical mass of 0.4 Msun, generally smaller than previously known coherent cores identified by Goodman et al. (1998), Caselli et al. (2002), and Pineda et al. (2010). We call these structures "droplets." We find that unlike previously known coherent cores, these structures are not virially bound by self-gravity and are instead predominantly confined by ambient pressure. The droplets have density profiles shallower than a critical Bonnor-Ebert sphere, and they have a velocity (VLSR) distribution consistent with the dense gas motions traced by NH3 emission. These results point to a potential formation mechanism through pressure compression and turbulent processes in the dense gas. We present a comparison with a magnetohydrodynamic simulation of a star-forming region, and we speculate on the relationship of droplets with larger, gravitationally bound coherent cores, as well as on the role that droplets and other coherent structures play in the star formation process. △ Less

Submitted 15 May, 2019; v1 submitted 26 September, 2018; originally announced September 2018.

Comments: Accepted by ApJ in April, 2019

Journal ref: 2019ApJ...877...93C

The Green Bank Ammonia Survey: Unveiling the Dynamics of the Barnard 59 star-forming Clump

Authors: E. Redaelli, F. O. Alves, P. Caselli, J. E. Pineda, R. K. Friesen, A. Chacón-Tanarro, C. D. Matzner, A. Ginsburg, E. Rosolowsky, J. Keown, S. S. R. Offner, J. Di Francesco, H. Kirk, P. C. Myers, A. Hacar, A. Cimatti, H. H. Chen, M. C. Chen, Y. M. Seo, K. I. Lee

Abstract: Understanding the early stages of star formation is a research field of ongoing development, both theoretically and observationally. In this context, molecular data have been continuously providing observational constraints on the gas dynamics at different excitation conditions and depths in the sources. We have investigated the Barnard 59 core, the only active site of star formation in the Pipe N… ▽ More Understanding the early stages of star formation is a research field of ongoing development, both theoretically and observationally. In this context, molecular data have been continuously providing observational constraints on the gas dynamics at different excitation conditions and depths in the sources. We have investigated the Barnard 59 core, the only active site of star formation in the Pipe Nebula, to achieve a comprehensive view of the kinematic properties of the source. These information were derived by simultaneously fitting ammonia inversion transition lines (1,1) and (2,2). Our analysis unveils the imprint of protostellar feedback, such as increasing line widths, temperature and turbulent motions in our molecular data. Combined with complementary observations of dust thermal emission, we estimate that the core is gravitationally bound following a virial analysis. If the core is not contracting, another source of internal pressure, most likely the magnetic field, is supporting it against gravitational collapse and limits its star formation efficiency. △ Less

Submitted 6 November, 2017; originally announced November 2017.

Comments: 18 pages, 18 figures

The Green Bank Ammonia Survey: Observations of Hierarchical Dense Gas Structures in Cepheus-L1251

Authors: Jared Keown, James Di Francesco, Helen Kirk, Rachel K. Friesen, Jaime E. Pineda, Erik Rosolowsky, Adam Ginsburg, Stella S. R. Offner, Paola Caselli, Felipe Alves, Ana Chacón-Tanarro, Anna Punanova, Elena Redaelli, Young Min Seo, Christopher D. Matzner, Michael Chun-Yuan Chen, Alyssa A. Goodman, How-Huan Chen, Yancy Shirley, Ayushi Singh, Hector G. Arce, Peter Martin, Philip C. Myers

Abstract: We use Green Bank Ammonia Survey observations of NH$_3$ (1,1) and (2,2) emission with 32'' FWHM resolution from a ~ 10 pc$^{2}$ portion of the Cepheus-L1251 molecular cloud to identify hierarchical dense gas structures. Our dendrogram analysis of the NH$_3$ data results in 22 top-level structures, which reside within 13 lower-level, parent structures. The structures are compact (0.01 pc… ▽ More We use Green Bank Ammonia Survey observations of NH$_3$ (1,1) and (2,2) emission with 32'' FWHM resolution from a ~ 10 pc$^{2}$ portion of the Cepheus-L1251 molecular cloud to identify hierarchical dense gas structures. Our dendrogram analysis of the NH$_3$ data results in 22 top-level structures, which reside within 13 lower-level, parent structures. The structures are compact (0.01 pc $\lesssim R_{eff} \lesssim$ 0.1 pc) and are spatially correlated with the highest H$_2$ column density portions of the cloud. We also compare the ammonia data to a catalog of dense cores identified by higher-resolution (18.2'' FWHM) Herschel Space Observatory observations of dust continuum emission from Cepheus-L1251. Maps of kinetic gas temperature, velocity dispersion, and NH$_3$ column density, derived from detailed modeling of the NH$_3$ data, are used to investigate the stability and chemistry of the ammonia-identified and Herschel-identified structures. We show that the dust and dense gas in the structures have similar temperatures, with median $T_{dust}$ and $T_K$ measurements of 11.7 $\pm$ 1.1 K and 10.3 $\pm$ 2.0 K, respectively. Based on a virial analysis, we find that the ammonia-identified structures are gravitationally dominated, yet may be in or near a state of virial equilibrium. Meanwhile, the majority of the Herschel-identified dense cores appear to be not bound by their own gravity and instead confined by external pressure. CCS $(2_0-1_0)$ and HC$_5$N $(9-8)$ emission from the region reveal broader line widths and centroid velocity offsets when compared to the NH$_3$ (1,1) emission in some cases, likely due to these carbon-based molecules tracing the turbulent outer layers of the dense cores. △ Less

Submitted 12 October, 2017; originally announced October 2017.

Comments: Accepted for publication in ApJ

The Green Bank Ammonia Survey: Dense Cores Under Pressure in Orion A

Authors: Helen Kirk, Rachel K. Friesen, Jaime E. Pineda, Erik Rosolowsky, Stella S. R. Offner, Christopher D. Matzner, Philip C. Myers, James Di Francesco, Paola Caselli, Felipe O. Alves, Ana Chacón-Tanarro, How-Huan Chen, Michael Chun-Yuan Chen, Jared Keown, Anna Punanova, Young Min Seo, Yancy Shirley, Adam Ginsburg, Christine Hall, Ayushi Singh, Héctor G. Arce, Alyssa A. Goodman, Peter Martin, Elena Redaelli

Abstract: We use gas temperature and velocity dispersion data from the Green Bank Ammonia Survey and core masses and sizes from the James Clerk Maxwell Telescope Gould Belt Survey to estimate the virial states of dense cores within the Orion A molecular cloud. Surprisingly, we find that almost none of the dense cores are sufficiently massive to be bound when considering only the balance between self-gravity… ▽ More We use gas temperature and velocity dispersion data from the Green Bank Ammonia Survey and core masses and sizes from the James Clerk Maxwell Telescope Gould Belt Survey to estimate the virial states of dense cores within the Orion A molecular cloud. Surprisingly, we find that almost none of the dense cores are sufficiently massive to be bound when considering only the balance between self-gravity and the thermal and non-thermal motions present in the dense gas. Including the additional pressure binding imposed by the weight of the ambient molecular cloud material and additional smaller pressure terms, however, suggests that most of the dense cores are pressure confined. △ Less

Submitted 17 August, 2017; originally announced August 2017.

Comments: Accepted for publication in ApJ Aug 10, 2017

The Role of Non-ionizing Radiation Pressure in Star Formation: The Stability of Cores and Filaments

Authors: Young Min Seo, Andrew N. Youdin

Abstract: Stars form when filaments and dense cores in molecular clouds fragment and collapse due to self-gravity. In the most basic analyses of gravitational stability, the competition between self-gravity and thermal pressure sets the critical (i.e. maximum stable) mass of spheres and the critical line density of cylinders. Previous work has considered additional support from magnetic fields and turbulenc… ▽ More Stars form when filaments and dense cores in molecular clouds fragment and collapse due to self-gravity. In the most basic analyses of gravitational stability, the competition between self-gravity and thermal pressure sets the critical (i.e. maximum stable) mass of spheres and the critical line density of cylinders. Previous work has considered additional support from magnetic fields and turbulence. Here, we consider the effects of non-ionizing radiation, specifically the inward radiation pressure force that acts on dense structures embedded in an isotropic radiation field. Using hydrostatic, isothermal models, we find that irradiation lowers the critical mass and line density for gravitational collapse, and can thus act as a trigger for star formation. For structures with moderate central densities, $\sim10^3$ cm$^{-3}$, the interstellar radiation field in the Solar vicinity has an order unity effect on stability thresholds. For more evolved objects with higher central densities, a significant lowering of stability thresholds requires stronger irradiation, as can be found closer to the Galactic center or near stellar associations. Even when strong sources of ionizing radiation are absent or extincted, our study shows that interstellar irradiation can significantly influence the star formation process. △ Less

Submitted 7 June, 2016; originally announced June 2016.

Comments: 13pages, 10 figures, accepted to MNRAS

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

Authors: Young Min Seo, Yancy L. Shirley, Paul Goldsmith, Derek Ward-Thompson, Jason M. Kirk, Markus Schmalzl, Jeong-Eun Lee, Rachel Friesen, Glen Langston, Joe Masters, Robert W. Garwood

Abstract: 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… ▽ More 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. △ Less

Submitted 17 March, 2015; v1 submitted 17 March, 2015; originally announced March 2015.

Comments: ApJ Accepted

On the internal dynamics of starless cores: stability of starless cores with internal motions and collapse dynamics

Authors: Young Min Seo, Seung Soo Hong, Yancy L. Shirley

Abstract: In order to understand the collapse dynamics of observed low-mass starless cores, we revise the conventional stability condition of hydrostatic Bonnor-Ebert spheres to take internal motions into account. Because observed starless cores resemble Bonnor-Ebert density structures, the stability and dynamics of the starless cores are frequently analyzed by comparing to the conventional stability condit… ▽ More In order to understand the collapse dynamics of observed low-mass starless cores, we revise the conventional stability condition of hydrostatic Bonnor-Ebert spheres to take internal motions into account. Because observed starless cores resemble Bonnor-Ebert density structures, the stability and dynamics of the starless cores are frequently analyzed by comparing to the conventional stability condition of a hydrostatic Bonnor-Ebert sphere. However, starless cores are not hydrostatic but have observed internal motions. In this study, we take gaseous spheres with a homologous internal velocity field and derive stability conditions of the spheres utilizing a virial analysis. We propose two limiting models of spontaneous gravitational collapse: the collapse of critical Bonnor-Ebert spheres and uniform density spheres. The collapse of these two limiting models are intended to provide the lower and the upper limits, respectively, of the infall speeds for a given density structure. The results of our study suggest that the stability condition sensitively depends on internal motions. A homologous inward motion with a transonic speed can reduce the critical size compared to the static Bonnor-Ebert sphere by more than a factor of two. As an application of the two limiting models of spontaneous gravitational collapse, we compare the density structures and infall speeds of the observed starless cores L63, L1544, L1689B, and L694-2 to the two limiting models. L1689B and L694-2 seem to have been perturbed to result in faster infall motions than for spontaneous gravitational collapse. △ Less

Submitted 10 April, 2013; originally announced April 2013.

Gravitational Instability of Rotating, Pressure-Confined, Polytropic Gas Disks With Vertical Stratification

Authors: Jeong-Gyu Kim, Woong-Tae Kim, Young Min Seo, Seung Soo Hong

Abstract: We investigate gravitational instability (GI) of rotating, vertically-stratified, pressure-confined, polytropic gas disks using linear stability analysis as well as analytic approximations. The disks are initially in vertical hydrostatic equilibrium and bounded by a constant external pressure. We find that GI of a pressure-confined disk is in general a mixed mode of the conventional Jeans and dist… ▽ More We investigate gravitational instability (GI) of rotating, vertically-stratified, pressure-confined, polytropic gas disks using linear stability analysis as well as analytic approximations. The disks are initially in vertical hydrostatic equilibrium and bounded by a constant external pressure. We find that GI of a pressure-confined disk is in general a mixed mode of the conventional Jeans and distortional instabilities, and is thus an unstable version of acoustic-surface-gravity waves. The Jeans mode dominates in weakly confined disks or disks with rigid boundaries. When the disk has free boundaries and is strongly pressure-confined, on the other hand, the mixed GI is dominated by the distortional mode that is surface-gravity waves driven unstable under own gravity and thus incompressible. We demonstrate that the Jeans mode is gravity-modified acoustic waves rather than inertial waves and that inertial waves are almost unaffected by self-gravity. We derive an analytic expression for the effective sound speed c_eff of acoustic-surface-gravity waves. We also find expressions for the gravity reduction factors relative to a razor-thin counterpart, appropriate for the Jeans and distortional modes. The usual razor-thin dispersion relation after correcting for c_eff and the reduction factors closely matches the numerical results obtained by solving a full set of linearized equations. The effective sound speed generalizes the Toomre stability parameter of the Jeans mode to allow for the mixed GI of vertically-stratified, pressure-confined disks. △ Less

Submitted 23 October, 2012; originally announced October 2012.

Comments: 18 pages, 13 figures; Accepted for publication in the ApJ