Search

Atmospheric escape of close-in exoplanets can be driven by high energy radiation from the host star. The planetary outflows interacting with the stellar wind may generate observable transit signals that depend on the strength of the stellar wind. We perform detailed radiation-hydrodynamics simulations of the atmospheric escape of hot Jupiters with including the wind from the host star in a self-consistent, dynamically coupled manner. We show that the planetary outflow is shaped by the balance between its thermal pressure and the ram pressure of the stellar wind. We use the simulation outputs to calculate the Lyman-$\alpha$ and H$\alpha$ transit signatures. Strong winds can confine the outflow and decrease the Lyman-$\alpha$ transit depth. Contrastingly, the wind effect on H$\alpha$ is weak because of the small contribution from the uppermost atmosphere of the planet. Observing both of the lines is important to understand the effect of the UV radiation and wind from the host. The atmospheric mass-loss rate is approximately independent of the strength of the wind. We also discuss the effect of the coronal mass ejections on the signatures. We argue that around M dwarfs the effect can be significant in every transit., Comment: 5 pages, 4 figures, accepted for publication in MNRAS

Transit observations have revealed the existence of atmospheric escape in several hot Jupiters. High energy photons from the host star heat the upper atmosphere and drive the hydrodynamic escape. The escaping atmosphere can interact with the stellar wind from the host star. We run radiation hydrodynamics simulations with non-equilibrium chemistry to investigate the wind effects on the escape and the transit signature. Our simulations follow the planetary outflow driven by the photoionization heating and the wind interaction in a dynamically coupled, self-consistent manner. We show that the planetary mass-loss rate is almost independent of the wind strength, which however affects the Ly-α transit depth considerably. But the Hα transit depth is almost independent of the wind strength because it is largely caused by the lower hot layer. We argue that observations of both lines can solve the degeneracy between the EUV flux from the host and the wind strength.

Photoelectric effect of dust grains by UV radiation is an important process for disk heating, but as a disk evolves, the amount of dust grains decreases. Photoeaporation is a disk dispersal process, which is caused by high-energy radiation. We perform a set of photoevaporation simulations solving hydrodynamics with radiative transfer and non-equilibrium chemistry in a self-consistent way. We run a series of simulations with varying the dust-to-gas mass ratio in a range . We show that H2 pumping and X-ray heating mainly contribute to the disk heating in case of and photoelectric effect mainly heats the gas in cases. The mass-loss profile changes significantly with respect to the main heating process. The outer disk is more efficiently dispersed when photoelectric effect is the main heating source.

Imai et al., The chemical diversity of low-mass protostellar sources has so far been recognized, and environmental effects are invoked as its origin. In this context, observations of isolated protostellar sources without the influence of nearby objects are of particular importance. Here, we report the chemical and physical structures of the low-mass Class 0 protostellar source IRAS 16544−1604 in the Bok globule CB 68, based on 1.3 mm Atacama Large Millimeter/submillimeter Array observations at a spatial resolution of ∼70 au that were conducted as part of the large program FAUST. Three interstellar saturated complex organic molecules (iCOMs), CH3OH, HCOOCH3, and CH3OCH3, are detected toward the protostar. The rotation temperature and the emitting region size for CH3OH are derived to be 131 ± 11 K and ∼10 au, respectively. The detection of iCOMs in close proximity to the protostar indicates that CB 68 harbors a hot corino. The kinematic structure of the C18O, CH3OH, and OCS lines is explained by an infalling–rotating envelope model, and the protostellar mass and the radius of the centrifugal barrier are estimated to be 0.08–0.30 M⊙ and, This paper makes use of the following ALMA data: ADS/JAO.ALMA#2016.1.01376.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), NSC and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. The authors acknowledge financial support from JSPS and MAEE under the Japan–France Integrated Action Program (SAKURA: 25765VC). D.J. is supported by NRC Canada and by an NSERC Discovery Grant. E.B. and G.B. acknowledge funding from the European Research Council under the European Union's Horizon 2020 research and innovation program, for the project "The Dawn of Organic Chemistry," grant agreement No. 741002. G.B. acknowledges support from the PID2020-117710GB-I00 grant funded by MCIN/AEI/10.13039/501100011033 and from the Unit of Excellence María de Maeztu 2020–2023 award to the Institute of Cosmos Sciences (CEX2019-00918-M). I.J.-S. has received partial support from the Spanish State Research Agency (AEI) through project Nos. PID2019-105552RB-C41 and MDM-2017-0737 Unidad de Excelencia "Maria de Maeztu" - Centro de Astrobiologia (CSIC-INTA). This study is supported by grants-in-aid from the Ministry of Education, Culture, Sports, Science, and Technology of Japan (18H05222, 18J11010, 19H05069, 19K14753, and 21K13954).

Ring structures are observed by (sub-)millimeter dust continuum emission in various circumstellar disks from early stages of Class 0 and I to late stage of Class II young stellar objects (YSOs). In this paper, we study one of the possible scenarios of such ring formation in early stage, which is coagulation of dust aggregates. The dust grains grow in an inside-out manner because the growth timescale is roughly proportional to the orbital period. The boundary of the dust evolution can be regarded as the growth front, where the growth time is comparable to the disk age. With radiative transfer calculations based on the dust coagulation model, we find that the growth front can be observed as a ring structure because dust surface density is sharply changed at this position. Furthermore, we confirm that the observed ring positions in the YSOs with an age of $\lesssim1$ Myr are consistent with the growth front. The growth front could be important to create the ring structure in particular for early stage of the disk evolution such as Class 0 and I sources., 15 pages, 12 figures, 1 table, accepted for publication in ApJ

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

We study the early growth of massive seed black holes (BHs) via accretion in protogalactic nuclei where the stellar bulge component is assembled, performing axisymmetric two-dimensional radiation hydrodynamical simulations. We find that when a seed BH with $M_\bullet \sim 10^5~M_\odot$ is embedded in dense metal-poor gas ($Z=0.01~Z_\odot$) with a density of $\gtrsim 100~{\rm cm}^{-3}$ and bulge stars with a total mass of $M_\star \gtrsim 100~M_\bullet$, a massive gaseous disk feeds the BH efficiently at rates of $\gtrsim 0.3-1~M_\odot~{\rm yr}^{-1}$ and the BH mass increases nearly tenfold within $\sim 2$ Myr. This rapid accretion phase lasts until a good fraction of the gas bounded within the bulge accretes onto the BH, although the feeding rate is regulated owing to strong outflows driven by ionizing radiation emitted from the accreting BH. The transient growing mode can be triggered for seed BHs formed in massive dark-matter halos with masses of $\gtrsim 10^9~M_\odot$ at $z\sim 15-20$ (the virial temperature is $T_{\rm vir}\simeq 10^5~{\rm K}$). The host halos are heavier and rarer than those of typical first galaxies, but are more likely to end up in quasar hosts by $z\simeq 6$. This mechanism naturally yields a mass ratio of $M_\bullet/M_\star >0.01$ higher than the value seen in the local universe and the existence of such overmassive BHs provides us a unique opportunity of detecting highly accreting seed BHs at $z\sim 15$ with AB magnitude of $m_{\rm AB} \sim26 - 29$ mag at $2~μ{\rm m}$ (rest-frame 10 eV) by the upcoming observations by the James Webb Space Telescope and Nancy Grace Roman Space Telescope., 28 pages, 11 figures, accepted for publication to ApJ. Movies of the simulation are available at https://www.youtube.com/playlist?list=PLqDrKvxmu0l4fyJ74u5ZYyQ8PmwuNZA32

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.

Debris disks are classically considered to be gas-less systems, but recent (sub)millimeter observations have detected tens of those with rich gas content. The origin of the gas component remains unclear; namely, it can be protoplanetary remnants and/or secondary products deriving from large bodies. In order to be protoplanetary in origin, the gas component of the parental protoplanetary disk is required to survive for $\gtrsim10{\,\rm Myr}$. However, previous models predict $\lesssim 10{\,\rm Myr}$ lifetimes because of efficient photoevaporation at the late stage of disk evolution. In the present study, we investigate photoevaporation of gas-rich, optically-thin disks around intermediate-mass stars at a late stage of the disk evolution. The evolved system is modeled as those where radiation force is sufficiently strong to continuously blow out small grains ($\lesssim 4 {\,\rm \mu m}$), which are an essential component for driving photoevaporation via photoelectric heating induced by stellar far-ultraviolet (FUV). We find that the grain depletion reduces photoelectric heating, so that FUV photoevaporation is not excited. Extreme-ultraviolet (EUV) photoevaporation is dominant and yields a mass-loss rate of $2$--$5\times10^{-10}(\Phi_{\rm EUV}/10^{41}{\,\rm s}^{-1})^{1/2}\,M_\odot\,{\rm yr}^{-1}$, where $\Phi_{\rm EUV}$ is the EUV emission rate. The estimated lifetimes of the gas component are $\sim 50 (M_{\rm disk}/10^{-2}\,M_\odot)(\Phi_{\rm EUV}/10^{41}\,{\rm s}^{-1})^{1/2}\,{\rm Myr}$ and depend on the ``initial'' disk mass at the point small grains have been depleted in the system. With an order estimation, we show that the gas component can survive for a much longer time around A-type stars than lower-mass stars. This trend is consistent with the higher frequency of gas-rich debris disks around A-type stars, implying the possibility of the gas component being protoplanetary remnants., Comment: Accepted for publication in ApJ

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

We report a study of the low-mass Class 0 multiple system VLA 1623AB in the Ophiuchus star-forming region, using H13CO+ (J = 3–2), CS (J = 5–4), and CCH (N = 3–2) lines as part of the ALMA Large Program FAUST. The analysis of the velocity fields revealed the rotation motion in the envelope and the velocity gradients in the outflows (about 2000 au down to 50 au). We further investigated the rotation of the circumbinary VLA 1623A disk, as well as the VLA 1623B disk. We found that the minor axis of the circumbinary disk of VLA 1623A is misaligned by about 12° with respect to the large-scale outflow and the rotation axis of the envelope. In contrast, the minor axis of the circumbinary disk is parallel to the large-scale magnetic field according to previous dust polarization observations, suggesting that the misalignment may be caused by the different directions of the envelope rotation and the magnetic field. If the velocity gradient of the outflow is caused by rotation, the outflow has a constant angular momentum and the launching radius is estimated to be 5–16 au, although it cannot be ruled out that the velocity gradient is driven by entrainments of the two high-velocity outflows. Furthermore, we detected for the first time a velocity gradient associated with rotation toward the VLA 16293B disk. The velocity gradient is opposite to the one from the large-scale envelope, outflow, and circumbinary disk. The origin of its opposite gradient is also discussed.