Disk Atmosphere Research Articles

Effect of Time-varying X-Ray Emission from Stellar Flares on the Ionization of Protoplanetary Disks

X-rays have significant impacts on cold, weakly ionized protoplanetary disks by increasing the ionization rate and driving chemical reactions. Stellar flares are explosions that emit intense X-rays and are the unique source of hard X-rays with an energy of ≳10 keV in the protoplanetary disk systems. Hard X-rays should be carefully taken into account in models as they can reach the disk midplane as a result of scattering in the disk atmospheres. However, previous models are insufficient to predict the hard X-ray spectra because of simplifications in flare models. We develop a model of X-ray spectra of stellar flares based on observations and flare theories. The flare temperature and nonthermal electron emissions are modeled as functions of flare energy, which allows us to better predict the hard X-ray photon flux than before. Using our X-ray model, we conduct radiative transfer calculations to investigate the impact of flare hard X-rays on disk ionization, with a particular focus on the protoplanetary disk around a T Tauri star. We demonstrate that for a flare with an energy of 1035 erg, X-ray photons with ≳5 keV increase the ionization rates more than galactic cosmic rays down to z ≈ 0.1R. The contribution of flare X-rays to the ionization at the midplane depends on the disk parameters such as disk mass and dust settling. We also find that the 10 yr averaged X-rays from multiple flares could certainly contribute to the ionization. These results emphasize the importance of stellar flares on the disk evolution.

XRISM Spectroscopy of the Fe Kα Emission Line in the Seyfert Active Galactic Nucleus NGC 4151 Reveals the Disk, Broad-line Region, and Torus

We present an analysis of the first two XRISM/Resolve spectra of the well-known Seyfert-1.5 active galactic nucleus (AGN) in NGC 4151, obtained in 2023 December. Our work focuses on the nature of the narrow Fe K α emission line at 6.4 keV, the strongest and most common X-ray line observed in AGN. The total line is found to consist of three components. Even the narrowest component of the line is resolved with evident Fe K α,1 (6.404 keV) and K α,2 (6.391 keV) contributions in a 2:1 flux ratio, fully consistent with neutral gas with negligible bulk velocity. Subject to the limitations of our models, the narrowest and intermediate-width components are consistent with emission from optically thin gas, suggesting that they arise in a disk atmosphere and/or wind. Modeling the three line components in terms of Keplerian broadening, they are readily associated with (1) the inner wall of the “torus,” (2) the innermost optical “broad-line region” (or “X-ray BLR”), and (3) a region with a radius of r ≃ 100 GM/c 2 that may signal a warp in the accretion disk. Viable alternative explanations of the broadest component include a fast-wind component and/or scattering; however, we find evidence of variability in the narrow Fe K α line complex on timescales consistent with small radii. The best-fit models are statistically superior to simple Voigt functions, but when fit with Voigt profiles the time-averaged lines are consistent with a projected velocity broadening of FWHM =1600−200+400kms−1 . Overall, the resolution and sensitivity of XRISM show that the narrow Fe K line in AGN is an effective probe of all key parts of the accretion flow, as it is currently understood. We discuss the implications of these findings for our understanding of AGN accretion, future studies with XRISM, and X-ray-based black hole mass measurements.

An Expanding Accretion Disk and a Warm Disk Wind as Seen in the Spectral Evolution of HBC 722

We present a comprehensive analysis of the post-outburst evolution of the FU Ori object HBC 722 in optical/near-infrared (NIR) photometry and spectroscopy. Using a modified viscous accretion disk model, we fit the outburst epoch spectral energy distribution to determine the physical parameters of the disk, including Ṁacc=10−4.0M⊙ yr−1, R inner = 3.65 R ⊙, i = 79°, and a maximum disk temperature of Tmax=5700 K. We then use a decade of optical/NIR spectra to demonstrate a changing accretion rate drives the visible-range photometric variation, while the NIR shows the outer radius of the active accretion disk expands outward as the outburst progresses. We also identify the major components of the disk system: a plane-parallel disk atmosphere in Keplerian rotation and a two-part warm disk wind that is collimated near the star and wide-angle at larger radii. The wind is traced by classic wind lines, and appears as a narrow, low-velocity, deep absorption component in several atomic lines spanning the visible spectrum and in the CO 2.29 μm band. We compare the wind lines to those computed from wind models for other FU Ori systems and rapidly accreting young stellar disks and find a 4000–6000 K wind can explain the observed line profiles. Fitting the progenitor spectrum, we find M * = 0.2 M ⊙ and Ṁprogenitor=7.8×10−8M⊙yr−1 . Finally, we discuss HBC 722 relative to V960 Mon, another FU Ori object we have previously studied in detail.

PRODIGE – envelope to disk with NOEMA

Context. Complex organic molecules (COMs) have been found toward low-mass protostars, but the origins of the COM emission are still unclear. It can be associated with, for example, hot corinos, outflows, and/or accretion shock and disk atmospheres. Aims. We aim to disentangle the origin of the COM emission toward the chemically rich protobinary system SVS13A using six O-bearing COMs. Methods. We conducted NOrthern Extended Millimeter Array observations toward SVS13A as part of the PROtostars & DIsks: Global Evolution (PRODIGE) program. Our previous DCN observations reveal a possible infalling streamer, which may affect the chemistry of the central protobinary by inducing accretion outbursts and/or shocked gas. We further analyzed six O-bearing COMs: CH3OH, aGg’- (CH2OH)2, C2H5OH, CH2(OH)CHO, CH3CHO, and CH3OCHO. Although the COM emission is not spatially resolved, we constrained the source sizes to ≲0.3–0.4 arcsec (90–120 au) by conducting uv-domain Gaussian fitting. Interestingly, the high-spectral-resolution data reveal complex line profiles with multiple peaks; although the line emission is likely dominated by the secondary, VLA4A, at VLSR = 7.36 km s−1, the numbers of peaks (~2–5), the velocities, and the linewidths of these six O-bearing COMs are different. The local thermodynamic equilibrium (LTE) fitting unveils differences in excitation temperatures and emitting areas among these COMs. We further conducted multiple-velocity-component LTE fitting to decompose the line emission into different kinematic components. As a result, the emission of these COMs is decomposed into up to six velocity components from the LTE modeling. The physical conditions (temperature, column density, and source size) of these components from each COM are obtained, and Markov chain Monte Carlo sampling was performed to test the fitting results. Results. We find a variety in excitation temperatures (100–500 K) and source sizes (D ~ 10–70 au) from these kinematic components from different COMs. The emission of each COM can trace several components, and different COMs most likely trace different regions. Conclusions. Given this complex structure, we suggest that the central region is inhomogeneous and unlikely to be heated by only protostellar radiation. We conclude that accretion shocks induced by the large-scale infalling streamer likely exist and contribute to the complexity of the COM emission. This underlines the importance of high-spectral-resolution data when analyzing COM emission in protostars and deriving relative COM abundances.

X-Ray and Radio Monitoring of the Neutron Star Low-mass X-Ray Binary 1A 1744-361: Quasiperiodic Oscillations, Transient Ejections, and a Disk Atmosphere

We report on X-ray (NICER/NuSTAR/MAXI/Swift) and radio (MeerKAT) timing and spectroscopic analysis from a 3 month monitoring campaign in 2022 of a high-intensity outburst of the dipping neutron star low-mass X-ray binary 1A 1744−361. The 0.5–6.8 keV NICER X-ray hardness–intensity and color–color diagrams of the observations throughout the outburst suggest that 1A 1744−361 spent most of its outburst in an atoll-state, but we show that the source exhibited Z-state-like properties at the peak of the outburst, similar to a small sample of other atoll-state sources. A timing analysis with NICER data revealed several instances of an ≈8 Hz quasiperiodic oscillation (QPO; fractional rms amplitudes of ∼5%) around the peak of the outburst, the first from this source, which we connect to the normal branch QPOs seen in the Z-state. Our observations of 1A 1744−361 are fully consistent with the idea of the mass accretion rate being the main distinguishing parameter between atoll- and Z-states. Radio monitoring data by MeerKAT suggests that the source was at its radio-brightest during the outburst peak, and that the source transitioned from the “island” spectral state to the “banana” state within ∼3 days of the outburst onset, launching transient jet ejecta. The observations present the strongest evidence for radio flaring, including jet ejecta, during the island-to-banana spectral state transition at low accretion rates (atoll-state). The source also exhibited Fe xxv, Fe xxvi Kα, and Kβ X-ray absorption lines, whose origins likely lie in an accretion disk atmosphere.

Mapping the Vertical Gas Structure of the Planet-hosting PDS 70 Disk

PDS 70 hosts two massive, still-accreting planets and the inclined orientation of its protoplanetary disk presents a unique opportunity to directly probe the vertical gas structure of a planet-hosting disk. Here, we use high-spatial-resolution (≈0.″1; 10 au) observations in a set of CO isotopologue lines and HCO+ J = 4−3 to map the full 2D (r, z) disk structure from the disk atmosphere, as traced by 12CO, to closer to the midplane, as probed by less abundant isotopologues and HCO+. In the PDS 70 disk, 12CO traces a height of z/r ≈ 0.3, 13CO is found at z/r ≈ 0.1, and C18O originates at, or near, the midplane. The HCO+ surface arises from z/r ≈ 0.2 and is one of the few non-CO emission surfaces constrained with high-fidelity in disks to date. In the 12CO J = 3−2 line, we resolve a vertical dip and steep rise in height at the cavity wall, making PDS 70 the first transition disk where this effect is directly seen in line-emitting heights. In the outer disk, the CO emission heights of PDS 70 appear typical for its stellar mass and disk size and are not substantially altered by the two inner embedded planets. By combining CO isotopologue and HCO+ lines, we derive the 2D gas temperature structure and estimate a midplane CO snowline of ≈ 56–85 au. This implies that both PDS 70b and 70c are located interior to the CO snowline and are likely accreting gas with a high C/O ratio of ≈ 1.0, which provides context for future planetary atmospheric measurements from, e.g., JWST, and for properly modeling their formation histories.

JOYS: MIRI/MRS spectroscopy of gas-phase molecules from the high-mass star-forming region IRAS 23385+6053

Context. Space-based mid-infrared (IR) spectroscopy is a powerful tool for the characterization of important star formation tracers of warm gas which are unobservable from the ground. The previous mid-IR spectra of bright high-mass protostars with the Infrared Space Observatory (ISO) in the hot-core phase typically show strong absorption features from molecules such as CO2, C2H2, and HCN. However, little is known about their fainter counterparts at earlier stages. Aims. We aim to characterize the gas-phase molecular features in James Webb Space Telescope Mid-Infrared Instrument Medium Resolution Spectrometer (MIRI/MRS) spectra of the young and clustered high-mass star-forming region IRAS 23385+6053. Methods. Spectra were extracted from several locations in the MIRI/MRS field of view, targeting two mid-IR sources tracing embedded massive protostars as well as three H2 bright outflow knots at distances of >8000 au from the multiple. Molecular features in the spectra were fit with local thermodynamic equilibrium (LTE) slab models, with their caveats discussed in detail. Results. Rich molecular spectra with emission from CO, H2, HD, H2O, C2H2, HCN, CO2, and OH are detected towards the two mid-IR sources. However, only CO and OH are seen towards the brightest H2 knot positions, suggesting that the majority of the observed species are associated with disks or hot core regions rather than outflows or shocks. The LTE model fits to 12CO2, C2H2, HCN emission suggest warm 120–200 K emission arising from a disk surface around one or both protostars. The abundances of CO2 and C2H2 of ~10−7 are consistent with previous observations of high-mass protostars. Weak ~500 K H2O emission at ~6–7 µm is detected towards one mid-IR source, whereas 250–1050 K H2O absorption is found in the other. The H2O absorption may occur in the disk atmosphere due to strong accretion-heating of the midplane, or in a disk wind viewed at an ideal angle for absorption. CO emission may originate in the hot inner disk or outflow shocks, but NIRSpec data covering the 4.6 µm band head are required to determine the physical conditions of the CO gas, as the high temperatures seen in the MIRI data may be due to optical depth. OH emission is detected towards both mid-IR source positions and one of the shocks, and is likely excited by water photodissociation or chemical formation pumping in a highly non-LTE manner. Conclusions. The observed molecular spectra are consistent with disks having already formed around two protostars in the young IRAS 23385+6054 system. Molecular features mostly appear in emission from a variety of species, in contrast to the more evolved hot core phase protostars which typically show only absorption; however, further observations of young high-mass protostars are needed to disentangle geometry and viewing angle effects from evolution.

High-energy neutrinos from merging stellar-mass black holes in active galactic nuclei accretion disc

ABSTRACT A population of binary stellar-mass black hole (BBH) mergers are believed to occur embedded in the accretion disc of active galactic nuclei (AGNs). In this Letter, we demonstrate that the jets from these BBH mergers can propagate collimatedly within the disc atmosphere along with a forward shock and a reverse shock forming at the jet head. Efficient proton acceleration by these shocks is usually expected before the breakout, leading to the production of TeV−PeV neutrinos through interactions between these protons and electron-radiating photons via photon–meson production. AGN BBH mergers occurring in the outer regions of the disc are more likely to produce more powerful neutrino bursts. Taking the host AGN properties of the potential GW190521 electromagnetic (EM) counterpart as an example, one expects ≳1 neutrino events detectable by IceCube if the jet is on-axis and the radial location of the merger is R ≳ 105Rg, where Rg is the gravitational radius of the supermassive BH. Neutrino bursts from AGN BBH mergers could be detected by IceCube following the observation of gravitational waves (GWs), serving as precursor signals before the detection of EM breakout signals. AGN BBH mergers are potential target sources for future joint GW, neutrino, and EM multi-messenger observations.

First detection of X-ray polarization in thermal state of LMC X-3: spectro-polarimetric study with IXPE

ABSTRACT We report a comprehensive spectro-polarimetric study of the black hole binary LMC X-3 using simultaneous Imaging X-ray Polarimetry Explorer (IXPE), NICER, and NuSTAR observations in 0.5–20 keV energy band. The broad-band energy spectrum (0.5–20 keV) with NICER and NuSTAR is well described by the disc emission of temperature ∼1.1 keV and a weak Comptonizing tail beyond ∼10 keV. This evidently suggests a disc-dominated spectral state of the source with disc contribution of ∼96 per cent. The lack of variability (rms ∼ 0.5 per cent) in the power spectrum further corroborates the high/soft nature of the source. A significant polarization degree (PD) of 3.04 ± 0.40 per cent (>7σ) at a polarization angle of −44.24° ± 3.77° (>7σ) is found in 2–8 keV energy range of IXPE. In addition, PD is seen to increase with energy up to ∼4.35 ± 0.98 per cent (>3σ) in 4–8 keV band. Further, we attempt to constrain the source spin (a*) using broad-band spectral modelling that indicates a weakly rotating ‘hole’ in LMC X-3 with $a_{*}=0.273_{-0.012}^{+0.011}\,\text{to}\,0.295_{-0.021}^{+0.008}$ (90 per cent confidence). Based on the spectro-polarimetric results, we infer that the polarization in LMC X-3 is resulted possibly due to the combined effects of the direct and/or reflected emissions from a partially ionized disc atmosphere. Finally, we discuss the relevance of our findings.

A resolved rotating disk wind from a young T Tauri star in the Bok globule CB 26

Context. The disk-outflow connection plays a key role in extracting excess angular momentum from a forming protostar. Although indications of jet rotation have been reported for a few objects, observational constraints of outflow rotation are still very scarce. We have previously reported the discovery of a small collimated molecular outflow from the edge-on T Tauri star-disk system in the Bok globule CB 26 that shows a peculiar velocity pattern, reminiscent of an outflow that corotates with the Keplerian disk. However, we could not ultimately exclude possible alternative explanations for the origin of the observed velocity field. Aims. We report new, high angular resolution millimeter-interferometric observations of CB 26 with the aim of revealing the morphology and kinematics of the outflow at the disk-outflow interface to unambiguously discriminate between the possible alternative explanations for the observed peculiar velocity pattern. Methods. The IRAM PdBI array and the 30 m telescope were used to observe HCO+(1–0) and H13CO+(1–0) at 3.3 mm and 12CO(2–1) at 1.3 mm in three configurations plus zerospacing, resulting in spectral line maps with angular resolutions of 3.″5 and 0.″5, respectively. The SMA was used to observe the HCO+(3–2) line at 1.1 mm with an angular resolution of 1.″35. Additional earlier observations of 13CO(1-0) at 2.7mm with an angular resolution of 1.″0, obtained with OVRO, are also used for the analysis. Using a physical model of the disk, which was derived from the dust continuum emission, we employed chemo-dynamical modeling combined with line radiative transfer calculations to constrain kinematic parameters of the system and to construct a model of the CO emission from the disk that allowed us to separate the emission of the disk from that of the outflow. Results. Our observations confirm the disk-wind nature of the rotating molecular outflow from CB 26 - YSO 1. The new high-resolution data reveal an X-shaped morphology of the CO emission close to the disk, and vertical streaks extending from the disk surface with a small half-opening angle of ≈7°, which can be traced out to vertical heights of ≈500 au. We interpret this emission as the combination of the disk atmosphere and a well-collimated disk wind, of which we mainly see the outer walls of the outflow cone. The decomposition of this emission into a contribution from the disk atmosphere and the disk wind allowed us to trace the disk wind down to vertical heights of ≈40 au, where it is launched from the surface of the flared disk at radii of RL ≈ 20–45 au. The disk wind is rotating with the same orientation and speed as the Keplerian disk and the velocity structure of the cone walls along the flow is consistent with angular momentum conservation. The observed CO outflow has a total gas mass of ≈ 10−3 M⊙, a dynamical age of τdyn ≈ 740 yr, and a total momentum flux of ṖCO ≈ 1.0 × 10−5 M⊙ km s−1 yr−1, which is nearly three orders of magnitude larger than the maximum thrust that can be provided by the luminosity of the central star. Conclusions. We conclude that photoevaporation cannot be the main driving mechanism for this outflow, but it must be predominantly a magnetohydrodynamic disk wind. It is thus far the best-resolved rotating disk wind observed to be launched from a circumstellar disk in Keplerian rotation around a low-mass young stellar object (YSO), albeit also the one with the largest launch radius. It confirms the observed trend that disk winds from Class I YSOs with transitional disks have much larger launch radii than jets ejected from Class 0 protostars.

Modeling the Multiwavelength Evolution of the V960 Mon System

Abstract We study the evolution of the FU Ori object V960 Mon since its outburst, using available multiwavelength photometric time series over 8 yr, complemented by several epochs of moderate-dispersion spectrophotometry. We find that the source fading can be well-described by a decrease in the temperature of the inner disk, which results from a combination of decreasing accretion rate and increasing inner disk radius. We model the system with a disk atmosphere model that produces the observed variations in multiband photometry (this paper) and high-resolution spectral lines (a companion paper).

Iceline variations driven by protoplanetary disc gaps

ABSTRACT The composition of forming planets is strongly affected by the protoplanetary disc’s thermal structure. This thermal structure is predominantly set by dust radiative transfer and viscous (accretional) heating and can be impacted by gaps – regions of low dust and gas density that can occur when planets form. The effect of variations in dust surface density on disc temperature has been poorly understood to date. In this work, we use the radiative transfer code MCMax to model the 2D dust thermal structure with individual gaps corresponding to planets with masses of 0.1 MJ –5 MJ and orbital radii of 3, 5, and 10 au. Low dust opacity in the gap allows radiation to penetrate deeper and warm the mid-plane by up to 16 K, but only for gaps located in the region of the disc where stellar irradiation is the dominant source of heating. In viscously heated regions, the mid-plane of the gap is relatively cooler by up to 100 K. Outside of the gap, broad radial oscillations in heating and cooling are present due to disc flaring. These thermal features affect local dust–gas segregation of volatile elements (H2O, CH4, CO2, and CO). We find that icelines experience dramatic shifts relative to gapless models: up to 6.5 au (or 71 per cent) closer to the star and 4.3 au (or 100 per cent) closer to the mid-plane. While quantitative predictions of iceline deviations will require more sophisticated models, which include transport, sublimation/condensation kinetics, and gas–dust thermal decoupling in the disc atmosphere, our results suggest that planet-induced iceline variations represent a potential feedback from the planet on to the composition of material it is accreting.

Line-driven winds from variable accretion discs

ABSTRACTWe use numerical hydrodynamics simulations to study line-driven winds launched from an accreting α-disc. Building on previous work where the driving radiation field is static, we compute a time-dependent radiation flux from the local, variable accretion rate of the disc. We find that prior to the establishment of a steady state in the disc, variations of $\sim 15{{\ \rm per\ cent}}$ in disc luminosity correlate with variations of ∼ 2–3 in the mass flux of the wind. After a steady state is reached, when luminosity variations drop to $\sim 3{{\ \rm per\ cent}}$, these correlations vanish as the variability in the mass flux is dominated by the intrinsic variability of the winds. This is especially evident in lower luminosity runs where intrinsic variability is higher due to a greater prevalence of failed winds. The changing mass flux occurs primarily due to the formation of clumps and voids near the disc atmosphere that propagate out into the low velocity part of the flow, a process that can be influenced by local variations in disc intensity. By computing the normalized standard deviation of the mass outflow, we show that the impact of luminosity variations on mass outflow is more visible at higher luminosity. However, the absolute change in mass outflow due to luminosity increases is larger for lower luminosity models due to the luminosity-mass flux scaling relation becoming steeper. We further discuss implications for CVs and AGN and observational prospects.

UV-driven chemistry as a signpost of late-stage planet formation

The chemical reservoir within protoplanetary disks has a direct impact on planetary compositions and the potential for life. A long-lived carbon- and nitrogen-rich chemistry at cold temperatures (≤ 50 K) is observed within cold and evolved planet-forming disks. This is evidenced by bright emission from small organic radicals in 1–10 Myr aged systems that would otherwise have frozen out onto grains within 1 Myr. We explain how the chemistry of a planet-forming disk evolves from a cosmic-ray/X-ray-dominated regime to a ultraviolet-dominated chemical equilibrium. This, in turn, will bring about a temporal transition in the chemical reservoir from which planets will accrete. This photochemical dominated gas phase chemistry develops as dust evolves via growth, settling and drift, and the small grain population is depleted from the disk atmosphere. A higher gas-to-dust mass ratio allows for deeper penetration of ultraviolet photons is coupled with a carbon-rich gas (C/O > 1) to form carbon-bearing radicals and ions. This further results in gas phase formation of organic molecules, which then would be accreted by any actively forming planets present in the evolved disk.

An Electron-scattering Time Delay in Black Hole Accretion Disks

Universal to black hole X-ray binaries, the high-frequency soft lag gets longer during the hard-to-intermediate state transition, evolving from ≲1 to ∼10 ms. The soft lag production mechanism is thermal disk reprocessing of nonthermal coronal irradiation. X-ray reverberation models account for the light-travel time delay external to the disk, but assume instantaneous reprocessing of the irradiation inside the electron-scattering-dominated disk atmosphere. We model this neglected scattering time delay as a random walk within an α-disk atmosphere, with approximate opacities. To explain soft lag trends, we consider a limiting case of the scattering time delay that we dub the thermalization time delay, t th; this is the time for irradiation to scatter its way down to the effective photosphere, where it gets thermalized, and then scatter its way back out. We demonstrate that t th plausibly evolves from being inconsequential for low mass accretion rates characteristic of the hard state, to rivaling or exceeding the light-travel time delay for characteristic of the intermediate state. However, our crude model confines t th to a narrow annulus near peak accretion power dissipation, so cannot yet explain in detail the anomalously long-duration soft lags associated with larger disk radii. We call for time-dependent models with accurate opacities to assess the potential relevance of a scattering delay.

Stratified Distribution of Organic Molecules at the Planet-formation Scale in the HH 212 Disk Atmosphere

Formamide (NH2CHO) is considered an important prebiotic molecule because of its potential to form peptide bonds. It was recently detected in the atmosphere of the HH 212 protostellar disk on the solar system scale where planets will form. Here we have mapped it and its potential parent molecules HNCO and H2CO, along with other molecules CH3OH and CH3CHO, in the disk atmosphere, studying its formation mechanism. Interestingly, we find a stratified distribution of these molecules, with the outer emission radius increasing from ∼24 au for NH2CHO and HNCO, to 36 au for CH3CHO, to 40 au for CH3OH, and then to 48 au for H2CO. More importantly, we find that the increasing order of the outer emission radius of NH2CHO, CH3OH, and H2CO is consistent with the decreasing order of their binding energies, supporting that they are thermally desorbed from the ice mantle on dust grains. We also find that HNCO, which has much lower binding energy than NH2CHO, has almost the same spatial distribution, kinematics, and temperature as NH2CHO, and is thus more likely a daughter species of desorbed NH2CHO. On the other hand, we find that H2CO has a more extended spatial distribution with different kinematics from NH2CHO, thus questioning whether it can be the gas-phase parent molecule of NH2CHO.

Probing Magnetic Fields in Protoplanetary Disk Atmospheres through Polarized Near-IR Light Scattered by Aligned Grains

Magnetic fields play essential roles in protoplanetary disks. Magnetic fields in the disk atmosphere are of particular interest, as they are connected to the wind-launching mechanism. In this work, we study the polarization of the light scattered off of magnetically aligned grains in the disk atmosphere, focusing on the deviation of the polarization orientation from the canonical azimuthal direction, which may be detectable in near-IR polarimetry with instruments such as VLT/SPHERE. We show with a simple disk model that the polarization can even be oriented along the radial (rather than azimuthal) direction, especially in highly inclined disks with toroidally dominated magnetic fields. This polarization reversal is caused by the anisotropy in the polarizability of aligned grains and is thus a telltale sign of such grains. We show that the near-IR light is scattered mostly by μm-sized grains or smaller at the τ = 1 surface and such grains can be magnetically aligned if they contain superparamagnetic inclusions. For comparison with observations, we generate synthetic maps of the ratios of U ϕ /I and Q ϕ /I, which can be used to infer the existence of (magnetically) aligned grains through a negative Q ϕ (polarization reversal) and/or a significant level of U ϕ /I. We show that two features observed in the existing data, an asymmetric distribution of U ϕ with respect to the disk minor axis and a spatial distribution of U ϕ that is predominantly positive or negative, are incompatible with scattering by spherical grains in an axisymmetric disk. They provide indirect evidence for scattering by aligned nonspherical grains.

Water Shielding in the Terrestrial Planet-forming Zone: Implication for Inner Disk Organics

The chemical composition of the inner region of protoplanetary disks can trace the composition of planetary-building material. The exact elemental composition of the inner disk has not yet been measured and tensions between models and observations still exist. Recent advancements have shown UV shielding to be able to increase the emission of organics. Here, we expand on these models and investigate how UV shielding may impact chemical composition in the inner 5 au. In this work, we use the model from Bosman et al. and expand it with a larger chemical network. We focus on the chemical abundances in the upper disk atmosphere where the effects of water UV shielding are most prominent and molecular lines originate. We find rich carbon and nitrogen chemistry with enhanced abundances of C2H2, CH4, HCN, CH3CN, and NH3 by >3 orders of magnitude. This is caused by the self-shielding of H2O, which locks oxygen in water. This subsequently results in a suppression of oxygen-containing species like CO and CO2. The increase in C2H2 seen in the model with the inclusion of water UV shielding allows us to explain the observed C2H2 abundance without resorting to elevated C/O ratios as water UV shielding induced an effectively oxygen-poor environment in oxygen-rich gas. Thus, water UV shielding is important for reproducing the observed abundances of hydrocarbons and nitriles. From our model result, species like CH4, NH3, and NO are expected to be observable with the James Webb Space Telescope (JWST).

CO Line Emission Surfaces and Vertical Structure in Midinclination Protoplanetary Disks

High spatial resolution CO observations of midinclination (≈30°–75°) protoplanetary disks offer an opportunity to study the vertical distribution of CO emission and temperature. The asymmetry of line emission relative to the disk major axis allows for a direct mapping of the emission height above the midplane, and for optically thick, spatially resolved emission in LTE, the intensity is a measure of the local gas temperature. Our analysis of Atacama Large Millimeter/submillimeter Array archival data yields CO emission surfaces, dynamically constrained stellar host masses, and disk atmosphere gas temperatures for the disks around the following: HD 142666, MY Lup, V4046 Sgr, HD 100546, GW Lup, WaOph 6, DoAr 25, Sz 91, CI Tau, and DM Tau. These sources span a wide range in stellar masses (0.50–2.10 M ⊙), ages (∼0.3–23 Myr), and CO gas radial emission extents (≈200–1000 au). This sample nearly triples the number of disks with mapped emission surfaces and confirms the wide diversity in line emitting heights (z/r ≈ 0.1 to ≳0.5) hinted at in previous studies. We compute the radial and vertical CO gas temperature distributions for each disk. A few disks show local temperature dips or enhancements, some of which correspond to dust substructures or the proposed locations of embedded planets. Several emission surfaces also show vertical substructures, which all align with rings and gaps in the millimeter dust. Combining our sample with literature sources, we find that CO line emitting heights weakly decline with stellar mass and gas temperature, which, despite large scatter, is consistent with simple scaling relations. We also observe a correlation between CO emission height and disk size, which is due to the flared structure of disks. Overall, CO emission surfaces trace ≈2–5× gas pressure scale heights (Hg) and could potentially be calibrated as empirical tracers of Hg.

Water UV-shielding in the Terrestrial Planet-forming Zone: Implications from Water Emission

Mid-infrared spectroscopy is one of the few ways to observe the composition of the terrestrial planet-forming zone, the inner few astronomical units, of protoplanetary disks. The species currently detected in the disk atmosphere, for example, CO, CO2, H2O, and C2H2, are theoretically enough to constrain the C/O ratio on the disk surface. However, thermochemical models have difficulties in reproducing the full array of detected species in the mid-infrared simultaneously. In an effort to get closer to the observed spectra, we have included water UV-shielding as well as more efficient chemical heating into the thermochemical code Dust and Lines. We find that both are required to match the observed emission spectrum. Efficient chemical heating, in addition to traditional heating from UV photons, is necessary to elevate the temperature of the water-emitting layer to match the observed excitation temperature of water. We find that water UV-shielding stops UV photons from reaching deep into the disk, cooling down the lower layers with a higher column. These two effects create a hot emitting layer of water with a column of 1–10 × 1018 cm−2. This is only 1%–10% of the water column above the dust τ = 1 surface at mid-infrared wavelengths in the models and represents <1% of the total water column.