Search

We have observed the very low-mass Class 0 protostar IRAS 15398−3359 at scales ranging from 50 to 1800 au, as part of the Atacama Large Millimeter/Submillimeter Array Large Program FAUST. We uncover a linear feature, visible in H2CO, SO, and C18O line emission, which extends from the source in a direction almost perpendicular to the known active outflow. Molecular line emission from H2CO, SO, SiO, and CH3OH further reveals an arc-like structure connected to the outer end of the linear feature and separated from the protostar, IRAS 15398−3359, by 1200 au. The arc-like structure is blueshifted with respect to the systemic velocity. A velocity gradient of 1.2 km s−1 over 1200 au along the linear feature seen in the H2CO emission connects the protostar and the arc-like structure kinematically. SO, SiO, and CH3OH are known to trace shocks, and we interpret the arc-like structure as a relic shock region produced by an outflow previously launched by IRAS 15398−3359. The velocity gradient along the linear structure can be explained as relic outflow motion. The origins of the newly observed arc-like structure and extended linear feature are discussed in relation to turbulent motions within the protostellar core and episodic accretion events during the earliest stage of protostellar evolution. The authors thank the anonymous reviewer for invaluable comments. They are also grateful to Hauyu Baobab Liu for useful discussions. This paper makes use of the following ALMA data set: ADS/JAO.ALMA#2018.1.01205.L (PI: Satoshi Yamamoto). ALMA is a partnership of the ESO (representing its member states), the NSF (USA) and NINS (Japan), together with the NRC (Canada) and the NSC and ASIAA (Taiwan), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by the ESO, the AUI/NRAO, and the NAOJ. The authors thank to the ALMA staff for their excellent support. This work is supported by the projects PRIN-INAF-MAIN-STREAM 2017 Protoplanetary disks seen through the eyes of new-generation instruments, ARIEL and the astrochemical link between circumstellar disks and planets, and the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programmes: The Dawn of Organic Chemistry (DOC), grant agreement No. 741002, and Astro-Chemistry Origins (ACO), grant No. 811312. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. I.J.-S. has received partial support from the Spanish FEDER (project number ESP2017-86582-C4-1-R) and the State Research Agency (AEI; project number PID2019-105552RB-C41). D.J. is supported by NRC Canada and by an NSERC Discovery Grant. This project is also supported by a Grant-in-Aid from Japan Society for the Promotion of Science (KAKENHI: Nos. 18H05222, 19H05069, 19K14753). Yuki Okoda thanks the Advanced Leading Graduate Course for Photon Science (ALPS) and Japan Society for the Promotion of Science (JSPS) for financial support. Peerreview

The density distribution of the inter-galactic medium is an uncertain but highly important issue in the study of cosmic reionization. It is expected that there are abundant gas clouds hosted by low-mass "minihalos" in the early universe, which act as photon sinks until photoevaporated by the emerging ultra-violet background (UVB) radiation. We perform a suite of radiation hydrodynamics simulations to study the photoevaporation of minihalos. Our simulations follow hydrodynamics, non-equilibrium chemistry, and the associated cooling processes in a self-consistent manner. We conduct a parametric study by considering a wide range of gas metallicity ($0\,Z_\odot \leq Z \leq 10^{-3}\,Z_\odot$), halo mass ($10^3 M_\odot \leq M \leq 10^8 M_\odot$), UVB intensity ($0.01 \leq J_{21} \leq 1$), and turn-on redshift of ionizing sources ($10\leq z_{\rm IN} \leq 20$). We show that small halos are evaporated in a few tens million years, whereas larger mass halos survive for ten times longer. We show that the gas mass evolution of a minihalo can be characterized by a scaling parameter that is given by a combination of the halo mass, background radiation intensity, and redshift. Efficient radiative cooling in metal-enriched halos induces fast condensation of the gas to form a dense, self-shielded core. The cold, dense core can become gravitationally unstable in halos with high metallicities. Early metal enrichment may allow star formation in minihalos during cosmic reionization.

We analyze multi-frequency, high-resolution continuum data obtained by ALMA and JVLA to study detailed structure of the dust distribution in the infant disk of a Class~0/I source, L1527 IRS. We find three clumps aligning in the north-south direction in the $7 {\rm \,mm}$ radio continuum image. The three clumps remain even after subtracting free-free contamination, which is estimated from the $1.3{\rm \,cm}$ continuum observations. The northern and southern clumps are located at a distance of $\sim 15{\rm \,au}$ from the central clump and are likely optically thick at $7{\rm \,mm}$ wavelength. The clumps have similar integrated intensities. The symmetric physical properties could be realized when a dust ring or spiral arms around the central protostar is projected to the plane of the sky. We demonstrates for the first time that such substructure may form even in the disk-forming stage, where the surrounding materials actively accrete toward a disk-protostar system.

Context. In recent years hydrodynamical (HD) models have become important to describe the gas kinematics in protoplanetary disks, especially in combination with models of photoevaporation and/or magnetically driven winds. Our aim is to investigate how vertical shear instability (VSI) could influence the thermally driven winds on the surface of protoplanetary disks. Aims. In this first part of the project, we focus on diagnosing the conditions of the VSI at the highest numerical resolution ever recorded, and suggest at what resolution per scale height we obtain convergence. At the same time, we want to investigate the vertical extent of VSI activity. Finally, we determine the regions where extreme UV (EUV), far-UV (FUV), and X-ray photons are dominant in the disk. Methods. We perform global HD simulations using the PLUTO code. We adopt a global isothermal accretion disk setup, 2.5D (2 dimensions, 3 components) which covers a radial domain from 0.5 to 5.0 and an approximately full meridional extension. Our simulation runs cover a resolution from 12 to 203 cells per scale height. Results. We determine 50 cells per scale height to be the lower limit to resolve the VSI. For higher resolutions, ≥50 cells per scale height, we observe the convergence for the saturation level of the kinetic energy. We are also able to identify the growth of the “body” modes, with higher growth rate for higher resolution. Full energy saturation and a turbulent steady state is reached after 70 local orbits. We determine the location of the EUV heated region defined by Σr = 1019 cm−2 to be at HR ~ 9.7 and the FUV–X-ray heated boundary layer defined by Σr = 1022 cm−2 to be at HR ~ 6.2, making it necessary to introduce a hot atmosphere. For the first time we report the presence of small-scale vortices in the r − Z plane between the characteristic layers of large-scale vertical velocity motions. Such vortices could lead to dust concentration, promoting grain growth. Our results highlight the importance of combining photoevaporation processes in the future high-resolution studies of turbulence and accretion processes in disks.

Recent multi-wavelength observations suggest that inner parts of protoplanetary disks (PPDs) have shorter lifetimes for heavier host stars. Since PPDs around high-mass stars are irradiated by strong ultra-violet radiation, photoevaporation may provide an explanation for the observed trend. We perform radiation hydrodynamics simulations of photoevaporation of PPDs for a wide range of host star mass of $M_* =0.5$-$7.0 M_{\odot}$. We derive disk mass-loss rate $\dot{M}$, which has strong stellar dependence as $\dot{M} \approx 7.30\times10^{-9}(M_{*}/M_{\odot})^{2}M_{\odot}\rm{yr}^{-1}$. The absolute value of $\dot{M}$ scales with the adopted far-ultraviolet and X-ray luminosities. We derive the surface mass-loss rates and provide polynomial function fits to them. We also develop a semi-analytic model that well reproduces the derived mass-loss rates. The estimated inner disk lifetime decreases as the host star mass increases, in agreement with the observational trend. We thus argue that photoevaporation is a major physical mechanism for PPD dispersal for a wide range of the stellar mass and can account for the observed stellar mass dependence of the inner disk lifetime., Comment: 15 pages, 12figures, ApJ in press

The super-Eddington accretion onto intermediate seed BHs is a potential formation mode of supermassive black holes exceeding $10^9~M_\odot$ in the early universe. We here investigate how such rapid accretion may occur with finite amounts of heavy elements contained in the gas and dust. In our 1D radiation-hydrodynamics simulations, the radiative transfer is solved for both the direct UV lights emitted by an accretion disk and the diffuse IR lights thermally emitted by dust grains. Our results show that the radiative force by the IR lights causes a strong feedback to regulate the mass accretion. The resulting mean accretion rate is lower with the higher metallicity, and there is the critical metallicity $Z \sim 10^{-2}~Z_\odot$, above which the super-Eddington accretion is prevented by the radiation pressure of the IR lights. With this taken into account, we examine if the dusty super-Eddington accretion occurs in young galaxies using a simple model. We show that a sufficient number of galaxies at $z \gtrsim 10$ can be such potential sites if BHs accrete the cold dense gas with $T \sim 10^2$ K, approximately the thermal equilibrium value at $Z = 10^{-2}~Z_\odot$. We argue that the efficiency of the BH growth via the rapid accretion depends on the metallicity, and that the metallicity slightly lower than $10^{-2}~Z_\odot$ provides a chance for the most efficient growth., Comment: 15 pages, 9 figures. Accepted for publication in MNRAS. Eq. (14) and (21) were corrected

We perform a suite of 3D radiation hydrodynamics simulations of photoevaporation of molecular gas clumps illuminated by external massive stars. We study the fate of solar-mass clumps and derive their lifetimes with varying the gas metallicity over a range of $10^{-3} \,Z_\odot \leq Z \leq Z_\odot $. Our simulations incorporate radiation transfer of far ultraviolet (FUV) and extreme ultraviolet (EUV) photons, and follow atomic/molecular line cooling and dust-gas collisional cooling. Nonequilibrium chemistry is coupled with the radiative transfer and hydrodynamics in a self-consistent manner. We show that radiation-driven shocks compress gas clumps to have a volume that is set by the pressure-equilibrium with the hot ambient gas. Radiative cooling enables metal-rich clumps to condense and to have small surface areas, where photoevaporative flows are launched. For our fiducial set-up with an O-type star at a distance of 0.1 parsec, the resulting photoevaporation rate is as small as $\sim 10^{-5} M_{\odot}/{\rm yr}$ for metal-rich clumps, but is larger for metal-poor clumps that have larger surface areas. The clumps are continuously accelerated away from the radiation source by the so-called rocket effect, and can travel over $\sim$1 parsec within the lifetime. We also study photoevaporation of clumps in a photo-dissociation region. Photoelectric heating is inefficient for metal-poor clumps that contain a smaller amount of grains, and thus they survive for over $10^5$ years. We conclude that the gas metallicity strongly affects the clump lifetime and thus determines the strength of feedback from massive stars in star-forming regions.

We perform a suite of radiation hydrodynamics simulations of photoevaporating disks with varying the metallicity in a wide range of $10^{-3} \, Z_\odot \leq Z \leq 10^{0.5} \, Z_\odot $. We follow the disk evolution for over $\sim 5000$ years by solving hydrodynamics, radiative transfer, and non-equilibrium chemistry. Our chemistry model is updated from the first paper of this series by adding X-ray ionization and heating. We study the metallicity dependence of the disk photoevaporation rate and examine the importance of X-ray radiation. In the fiducial case with solar metallicity, including the X-ray effects does not significantly increase the photoevaporation rate when compared to the case with ultra-violet (UV) radiation only. At sub-solar metallicities in the range of $Z \gtrsim 10^{-1.5} \, Z_\odot $, the photoevaporation rate increases as metallicity decreases owing to the reduced opacity of the disk medium. The result is consistent with the observational trend that disk lifetimes are shorter in low metallicity environments. Contrastingly, the photoevaporation rate decreases at even lower metallicities of $Z \lesssim 10^{-1.5} \, Z_\odot $, because dust-gas collisional cooling remains efficient compared to far UV photoelectric heating whose efficiency depends on metallicity. The net cooling in the interior of the disk suppresses the photoevaporation. However, adding X-ray radiation significantly increases the photoevaporation rate, especially at $Z \sim 10^{-2}\, Z_\odot$. Although the X-ray radiation itself does not drive strong photoevaporative flows, X-rays penetrate deep into the neutral region in the disk, increase the ionization degree there, and reduce positive charges of grains. Consequently, the effect of photoelectric heating by far UV radiation is strengthened by the X-rays and enhances the disk photoevaporation., Comment: 13 pages, 8 figures

Protoplanetary disks are thought to have lifetimes of several million years in the solar neighborhood, but recent observations suggest that the disk lifetimes are shorter in a low metallicity environment. We perform a suite of radiation hydrodynamics simulations of photoevaporation of protoplanetary disks to study the disk structure and its long-term evolution of $\sim 10000$ years, and the metallicity dependence of mass-loss rate. Our simulations follow hydrodynamics, extreme and far ultra-violet radiative transfer, and non-equilibrium chemistry in a self-consistent manner. Dust grain temperatures are also calculated consistently by solving the radiative transfer of the stellar irradiation and grain (re-)emission. We vary the disk gas metallicity over a wide range of $10^{-4}~ Z_\odot \leq Z \leq 10 ~Z_\odot$. The photoevaporation rate is lower with higher metallicity in the range of $10^{-1} \,Z_\odot \lesssim Z \lesssim 10 \,Z_\odot$, because dust shielding effectively prevents far-ultra violet (FUV) photons from penetrating into and heating the dense regions of the disk. The photoevaporation rate sharply declines at even lower metallicities in $10^{-2} \,Z_\odot \lesssim Z \lesssim 10^{-1}\,Z_\odot$, because FUV photoelectric heating becomes less effective than dust-gas collisional cooling. The temperature in the neutral region decreases, and photoevaporative flows are excited only in an outer region of the disk. At $10^{-4}\,Z_\odot \leq Z \lesssim 10^{-2}\,Z_\odot$, HI photoionization heating acts as a dominant gas heating process and drives photoevaporative flows with roughly a constant rate. The typical disk lifetime is shorter at $Z=0.3~Z_\odot$ than at $Z = Z_\odot$, being consistent with recent observations of the extreme outer galaxy., Comment: 25 pages, 17 figures, to appear in ApJ