Edge-on Disk Research Articles

JWST Imaging of Edge-on Protoplanetary Disks. III. Drastic Morphological Transformation Across the Mid-infrared in Oph163131

We present JWST broadband images of the highly inclined protoplanetary disk SSTc2d J163131.2-242627 (Oph163131) from 2.0 to 21 μm. The images show a remarkable evolution in disk structure with wavelength, quite different from previous JWST observations of other edge-on disks. At 2.0 and 4.4 μm, Oph163131 shows two scattering surfaces separated by a dark lane, typical of highly inclined disks. Starting at 7.7 μm, however, (1) the two linear nebulosities flanking the dark lane disappear; (2) the brighter nebula tracing the disk upper surface transitions into a compact central source distinctly larger than the JWST point-spread function and whose intrinsic size increases with wavelength; and (3) patches of extended emission appear at low latitudes, and at surprisingly large radii nearly twice that of the scattered light seen with Hubble Space Telescope and NIRCam, and of the gas. We interpret the compact central source as thermal emission from the star and the inner disk that is not seen directly, but which instead is able to progressively propagate to greater distances at longer wavelengths. The lack of sharp-edged structures in the extended patchy emission argues against the presence of shocks and suggests photoexcitation or stochastic heating of material smoothly flowing away from the star along the disk surface. Finally, the dark lane thickness decreases significantly between 0.6 and 4.4 μm, which indicates that the surface layers of Oph163131 lack grains larger than 1 μm.

EDIG-CHANGES

Context. The kinematic information of the extraplanar diffuse ionized gas (eDIG) around galaxies provides clues to the origin of the gas. Aims. The eDIG-CHANGES project studies the physical and kinematic properties of the eDIG around the CHANG-ES sample of nearby edge-on disk galaxies. Methods. We use a novel multi-slit narrow-band spectroscopy technique to obtain the spatial distribution of the spectral properties of the ionized gas around NGC 891, which is often regarded as an analogue of the Milky Way. We developed specific data reduction procedures for the multi-slit narrow-band spectroscopy data taken with the MDM 2.4 m telescope. The data presented in this paper cover the Hα and [N II]λλ6548, 6583 Å emission lines. Results. The eDIG traced by the Hα and [N II] lines shows an obvious asymmetric morphology, being brighter in the northeastern part of the galactic disk and extending a few kiloparsecs above and below the disk. Global variations in the [N II]/Hα line ratio suggest additional heating mechanisms for the eDIG at large heights beyond photoionization. We also construct position-velocity (PV) diagrams of the eDIG based on our optical multi-slit spectroscopy data and compare them to similar PV diagrams constructed with the H I data. The dynamics of the two gas phases are generally consistent with each other. Modelling the rotation curves at different heights from the galactic mid-plane suggests a vertical negative gradient in turnover radius and maximum rotation velocity, with magnitudes of approximately 3 kpc kpc−1 and 22 − 25 km s−1 kpc−1, respectively. Conclusions. Our measured vertical gradients of the rotation curve parameters suggest significant differential rotation of the ionized gas in the halo, often referred to as the lagging eDIG. Systematic study of the lagging eDIG, using the multi-slit narrow-band spectroscopy technique developed in our eDIG-CHANGES project, will help us to better understand the dynamics of the ionized gas in the halo.

Recovering chemical bimodalities in observed edge-on stellar disks: Insights from AURIGA simulations

The well-known bimodal distribution of Milky Way disk stars in the [α/Fe]–metallicity plane is often used to define thick and thin disks. In external edge-on galaxies, there have been attempts to identify this type of bimodality using integral field spectroscopy (IFS) data. However, for unresolved stellar populations, observations only contain integrated information, making these studies challenging. We assessed the ability to recover chemical bimodalities in IFS observations of edge-on galaxies, using 24 Milky Way-mass galaxies from the AURIGA zoom-in cosmological simulations. We first analyzed the distribution of single stellar particles in the [Mg/Fe]–[Fe/H] plane, finding that bimodality is frequent but not ubiquitous and often unclear. Then we produced mock IFS [Mg/Fe] and [Fe/H] maps of galaxies seen edge on, and considered integrated stellar-population properties (projected and spatially binned). We investigated how the distribution of stars in the [Mg/Fe]–[Fe/H] plane is affected by edge-on projection and spatial binning. Bimodality is preserved, while distributions change their shapes. Naturally, broad distributions of individual star particles are narrowed into smaller [Mg/Fe] and [Fe/H] ranges for spatial bins. We observe continuous distributions from high [Mg/Fe] and low [Fe/H], to lower [Mg/Fe] values and higher [Fe/H]. Despite being continuous, these distributions are bimodal in most cases. The overlap in [Fe/H] is small, and different [Mg/Fe] components show up as peaks instead of sequences (even when the latter are present for individual particles). The larger the spatial bins, the narrower the [Mg/Fe]–[Fe/H] distribution. This narrowing helps amplify the density of different [Mg/Fe] peaks, often leading to a clearer bimodality in mock IFS observations than for original star particles. We also assessed the correspondence of chemical bimodalities with the distinction between geometric thick and thin disks. Their individual particles have different distributions, but mostly overlap in [Mg/Fe] and [Fe/H]. However, integrated properties of geometric thick and thin disks in mock maps do mostly segregate into different regions of the [Mg/Fe]–[Fe/H] plane. In bimodal distributions, they correspond to the two distinct peaks. Our results show that this approach can be used for bimodality studies in future IFS observations of edge-on external galaxies.

CLOUDY modeling suggests a diversity of ionization mechanisms for diffuse extraplanar gas

Context. The ionization of diffuse gas located far above the energetic midplane OB stars poses a challenge to the commonly accepted notion that radiation from OB stars is the primary ionization source for gas in galaxies. Aims. We investigated the sources of ionizing radiation, specifically leaking midplane H II regions and/or in situ hot low-mass evolved stars (HOLMES), in extraplanar diffuse ionized gas (eDIG) in a sample of eight nearby (17−52 Mpc) edge-on disk galaxies observed with the Multi Unit Spectroscopic Explorer (MUSE). Methods. We constructed a model for the photoionization of eDIG clouds and the propagation of ionizing radiation through the eDIG using subsequent runs of CLOUDY photoionization code. Our model includes radiation originating both from midplane OB stars and in situ evolved stars and its dilution and processing as it propagates in the eDIG. Results. We fit the model to the data using the vertical line ratio profiles of our sample galaxies, and find that while the ionization by in situ evolved stars is insignificant for most of the galaxies in our sample, it may be able to explain the enhanced high-ionization lines in the eDIG of the green valley galaxy ESO 544−27. Conclusions. Our results show that while leaking radiation from midplane H II regions is the primary ionization source for eDIG, in situ evolved stars can play a significant part in ionizing extraplanar gas in galaxies with low star forming rates.

An Edge-on Regular Disk Galaxy at z = 5.289

While rotation-supported gas disks are known to exist as early as at z ≈ 7, it is still a general belief that stellar disks form late in the Universe. This picture is now being challenged by the observations from the James Webb Space Telescope (JWST), which have revealed a large number of disk-like galaxies that could be at z > 3, with some being candidates at z > 7. As an early formation of stellar disks will greatly impact our theory of galaxy formation and evolution, it is important to determine when such systems first emerged. Here we present D-CEERS-RUBIES-z5289 at z = 5.289 ± 0.001, the second confirmed stellar disk at z > 5, discovered using the archival JWST NIRCam imaging and NIRSpec spectroscopic data. This galaxy has a highly regular edge-on disk morphology, extends to ∼6.2 kpc along its major axis, and has an effective radius of ∼1.3–1.4 kpc. Such a large stellar disk is yet to be produced in numerical simulations. By analyzing its 10-band spectral energy distribution using four different tools, we find that it has a high stellar mass of 109.5–10.0 M ⊙. Its age is in the range of 330–510 Myr, and it has a mild star formation rate of 10–30 M ⊙ yr−1. While the current spectroscopic data do not allow the derivation of its rotation curve, the width of its Hα line from the partial slit coverage on one side of the disk reaches ∼345 km s−1, which suggests that it could have a significant contribution from rotation.

Illuminating the Incidence of Extraplanar Dust Using Ultraviolet Reflection Nebulae with GALEX

Circumgalactic dust grains trace the circulation of mass and metals between star-forming regions and gaseous galactic halos, giving insight into feedback and tidal stripping processes. We perform a search for ultraviolet (UV) reflection nebulae produced by extraplanar dust around 551 nearby (D < 100 Mpc), edge-on disk galaxies using archival near-ultraviolet and far-ultraviolet images from the Galaxy Evolution Explorer (GALEX), accounting for the point-spread function (FWHM = 4″–5″). We detect extraplanar emission ubiquitously in stacks of galaxies binned by morphology and star formation rate, with scale heights of h h = 1–2.3 kpc and ≈10% of the total (reddened) flux in the galaxy found beyond the B-band isophotal level of μ B = 25 mag arcsec−2. This emission is detected in 7% of the individual galaxies, and an additional one-third have at least 5% of their total flux found beyond μ B = 25 mag arcsec−2 in a disk component. The extraplanar luminosities and colors are consistent with reflection nebulae rather than stellar halos and indicate that, on average, disk galaxies have an extraplanar dust mass of 5%–15% of that in their interstellar medium. This suggests that recycled material composes at least a third of the inner circumgalactic medium (R < 10 kpc) in ∼L* galaxies.

CHANG-ES. XXXII. Spatially Resolved Thermal–Nonthermal Separation from Radio Data Alone—New Probes into NGC 3044 and NGC 5775

We have carried out spatially resolved thermal–nonthermal separation on two edge-on galaxies, NGC 3044 and NGC 5775, using only radio data. Narrowband imaging within a frequency band that is almost contiguous from 1.25 to 7.02 GHz (L band, S band, and C band) has allowed us to fit spectra and construct thermal, nonthermal, and nonthermal spectral index maps. This method does not require any ancillary Hα and IR data or rely on dust corrections that are challenging in edge-on galaxies. For NGC 3044, at 15″ resolution, we find a median thermal fraction of ∼13% with an estimated uncertainty in this fraction of ∼50% at 4.13 GHz. This compares well with the Hα mixture method results. We uncovered evidence for a vertical outflow feature reaching at least z ∼ 3.5 kpc in projection above the plane, reminiscent of M82's starburst wind. For the higher star formation rate galaxy, NGC 5775, at 12″ resolution, we find a median thermal fraction of 44% at 4.13 GHz with an estimated error on this fraction of 17%. Both galaxies show a change of slope (flattening) in L band. These results suggest that a radio-only method for separating thermal from nonthermal emission is not only feasible, but able to reveal new features that might otherwise be obscured in edge-on disks.

Resolving Twin Jets and Twin Disks with JWST and ALMA: The Young WL 20 Multiple System

We report the discovery of jets emanating from pre-main-sequence objects exclusively at mid-infrared wavelengths, enabled by the superb sensitivity of JWST’s Mid-Infrared Instrument Medium-Resolution Spectrometer. These jets are observed only in lines of [Ni ii], [Fe ii], [Ar ii], and [Ne ii]. The H2 emission, imaged in eight distinct transitions, has a completely different morphology, exhibiting a wide-angled, biconical shape, symmetrically distributed about the jet axes. Synergistic high-resolution Atacama Large Millimeter/submillimeter Array (ALMA) observations resolve a pair of side-by-side edge-on accretion disks lying at the origin of the twin mid-infrared jets. Assuming coevality of the components of the young multiple system under investigation, the system age is at least (2–2.5) × 106 yr, despite the discrepantly younger age inferred from the spectral energy distribution of the combined edge-on disk sources. The later system evolutionary stage is corroborated by ALMA observations of CO(2−1), 13CO(2−1), and C18O(2−1), which show no traces of molecular outflows or remnant cavity walls. Consequently, the observed H2 structures must have their origins in wide-angled disk winds, in the absence of any ambient, swept-up gas. In the context of recent studies of protostars, we propose an outflow evolutionary scenario in which the molecular gas component dominates in the youngest sources, whereas the fast, ionized jets dominate in the oldest sources, as is the case for the twin jets discovered in the WL 20 system.

ϵ Sagittarii: An Extreme Rapid Rotator with a Decretion Disk

We report high-precision, multiwavelength linear-polarization observations of the bright B9 (or A0) star ϵ Sagittarii. The polarization shows the distinctive wavelength dependence expected for a rapidly rotating star. Analysis of the polarization data reveals an angular rotation rate ω (=Ω/Ωcrit) of 0.995 or greater, the highest yet measured for a star in our Galaxy. An additional wavelength-independent polarization component is attributed to electron scattering in a low-density, edge-on gas disk that also produces the narrow absorption components seen in the spectrum. Several properties of the star (polarization due to a disk, occasional weak Hα emission, and multiple periodicities seen in space photometry) resemble those of Be stars, but the level of activity in all cases is much lower than that of typical Be stars. The stellar properties are inconsistent with single-rotating-star evolutionary tracks, indicating that it is most likely a product of binary interaction. The star is an excellent candidate for observation by interferometry, optical spectropolarimetry to detect the Öhman effect, and ultraviolet polarimetry, any of which would allow its extreme rotation to be tested and its stellar properties to be refined.

The characterisation of water ice in debris discs: implications for JWST scattered light observations

Abstract Water ice plays a crucial role throughout the different stages of planetary evolution and is abundant in the Universe. However, its presence and nature in debris discs of exoplanetary systems are not yet strongly established observationally. In this study, we quantify and discuss the impact of ice parameters such as volume fraction ${\mathcal {F}}_{\rm ice}$, blow-out grain size, size distribution, and its phase on the observational appearance of debris discs, considering the diverse nature of these systems around stellar spectral types ranging from A to M. Our findings reveal that the prominent ice features at approximately 2.7 and 3.3 μm depend on both the water ice fraction ${\mathcal {F}}_{\rm ice}$ and the scattering angle, with backscattering geometries yielding the most prominent signatures. When the phase function is considered and data are not background limited, strong forward and backward scattering (near edge-on discs) are expected to yield the strongest detections in images/spectra for A or F-type stars, while scattering angle matters less for later type stars. The Fresnel peak at 3.1 μm serves as a viable discriminant for the transitional phase (crystalline/amorphous), while simultaneously constraining the water ice temperature. For JWST imaging, we find that the F356W and F444W filter combination is most effective for constraining the grain size distribution, while the F356W and F277W filter combination provides better constraints on the ice fraction ${\mathcal {F}}_{\rm ice}$ in debris discs. However, degeneracy between the grain size distribution and ice fraction when using photometric flux ratios means that obtaining robust constraints will likely require more than two filters, or spectroscopic data.

The β Pictoris b Hill sphere transit campaign

The β Pictoris system is the closest known stellar system with directly detected gas giant planets, an edge-on circumstellar disc, and evidence of falling sublimating bodies and transiting exocomets. The inner planet, β Pictoris c, has also been indirectly detected with radial velocity (RV) measurements. The star is a known δ Scuti pulsator, and the long-term stability of these pulsations opens up the possibility of indirectly detecting the gas giant planets through time delays of the pulsations due to a varying light travel time. We search for phase shifts in the δ Scuti pulsations consistent with the known planets β Pictoris b and c and carry out an analysis of the stellar pulsations of β Pictoris over a multi-year timescale. We used photometric data collected by the BRITE-Constellation, bRing, ASTEP, and TESS to derive a list of the strongest and most significant δ Scuti pulsations. We carried out an analysis with the open-source python package maelstrom to study the stability of the pulsation modes of β Pictoris in order to determine the long-term trends in the observed pulsations. We did not detect the expected signal for β Pictoris b or β Pictoris c. The expected time delay is 6 s for β Pictoris c and 24 s for β Pictoris b. With simulations, we determined that the photometric noise in all the combined data sets cannot reach the sensitivity needed to detect the expected timing drifts. An analysis of the pulsational modes of β Pictoris using maelstrom showed that the modes themselves drift on the timescale of a year, fundamentally limiting our ability to detect exoplanets around β Pictoris via pulsation timing.

Modelling predicts a molecule-rich disc around the AGB star L2 Puppis

ABSTRACT The nearby oxygen-rich AGB star L$_2$ Pup hosts a well-studied nearly edge-on disc. To date, discs around AGB stars have not been chemically studied in detail. By combining a parametrization commonly used for protoplanetary discs and archival ALMA observations, we retrieved an updated density and temperature structure of this disc. This physical model was then used as input to the first chemical model of an AGB disc. The model shows that the physical structure of the disc has a large impact on its chemistry, with certain species showing large changes in column density relative to a radial outflow, indicating that chemistry could be used as a tracer of discs that cannot be directly imaged. Despite its oxygen-rich nature, the daughter species formed within the disc are surprisingly carbon-rich. Two chemical regimes can be distinguished: cosmic-ray induced chemistry in the midplane and photochemistry induced by the interstellar radiation field in the outer regions. Certain complex organic molecules are formed in the midplane. This occurs via gas-phase chemistry only, as the disc is too warm for dust-gas chemistry. The photochemistry in the outer regions leads to the efficient formation of (long) carbon-chains. The predictions of the model allow us to tentatively put the disc’s age $\lesssim 10^5$ yr. Additional observations are necessary to better constrain the physical structure of L$_2$ Pup’s disc and are essential to test the predictions made by the chemical model. Our exploratory work paves the way for a more general study of the chemistry of AGB discs.

EDIG-CHANGES. II. Project Design and Initial Results on NGC 3556

The extraplanar diffuse ionized gas (eDIG) represents ionized gases traced by optical/UV lines beyond the stellar extent of galaxies. We herein introduce a novel multislit narrow-band spectroscopy method to conduct spatially resolved spectroscopy of the eDIG around a sample of nearby edge-on disk galaxies (eDIG-CHANGES). In this paper, we introduce the project design and major scientific goals, as well as a pilot study of NGC 3556 (M108). The eDIG is detected to a vertical extent of a few kiloparsecs above the disk, comparable to the X-ray and radio images. We do not see significant vertical variation of the [N ii]/Hα line ratio. A rough examination of the pressure balance between different circumgalactic medium phases indicates the magnetic field is in a rough pressure balance with the X-ray emitting hot gas and may play an important role in the global motion of both the eDIG and the hot gas in the lower halo. At the location of an Hubble Space Telescope/Cosmic Origins Spectrograph observed UV bright background active galactic nucleus ∼29 kpc from the center of NGC 3556, the magnetic pressure is much lower than that of the hot gas and the ionized gas traced by UV absorption lines, although the extrapolation of the pressure profiles may cause some biases in this comparison. By comparing the position–velocity diagrams of the optical and CO lines, we also find the dynamics of the two gas phases are consistent with each other, with no evidence of a global inflow/outflow and a maximum rotation velocity of ∼150 km s−1.

Dracula’s Chivito: Discovery of a Large Edge-on Protoplanetary Disk with Pan-STARRS

We report the serendipitous discovery of a large edge-on protoplanetary disk in Pan-STARRS (PS1) images. PS1 has five broadband filters designated as g P1, r P1, i P1, z P1, and y P1 with mean wavelengths 4866, 6215, 7545, 8679, and 9633 Å, respectively. The disk’s apparent size in the PS1 images is ≈11″, making this one of the largest known disks on the sky. It is likely a young system, still surrounded by the envelope, which is very faint but still visible in the PS1 images in the northern part (alternatively this structure could be filaments from the disk itself). We use the PS1 magnitudes and other available photometric data to construct the spectral energy distribution of the disk. An optical spectrum indicates that the obscured star is hot, most likely of type late A. We adopt a distance of 300 pc for this object based on Gaia Data Release 3 extinctions. We model the system using the HOCHUNK3D radiative transfer software and find that the system is consistent with a hot star of effective temperature 8000 K surrounded by a disk of size 1650 au and mass 0.2 M ⊙ at inclination 82°.

High-resolution Pan-STARRS and SMA Observations of IRAS 23077+6707: A Giant Edge-on Protoplanetary Disk

We present resolved images of IRAS 23077+6707 (“Dracula’s Chivito”) in 1.3 mm/225 GHz thermal dust and CO gas emission with the Submillimeter Array (SMA) and optical (0.5–0.8 μm) scattered light with the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS). The Pan-STARRS data show a bipolar distribution of optically scattering dust that is characteristic for disks observed at high inclinations. Its scattered light emission spans ∼14″, with two highly asymmetric filaments extending along the upper bounds of each nebula by ∼9″. The SMA data measure 1.3 mm continuum dust as well as 12CO, 13CO, and C18O J = 2 − 1 line emission over 12″–14″ extents, with the gas presenting the typical morphology of a disk in Keplerian rotation, in both position–velocity space and in each CO line spectrum. IRAS 23077+6707 has no reported distance estimate, but if it is located in the Cepheus star-forming region (180–800 pc), it would have a radius spanning thousands of astronomical units. Taken together, we infer IRAS 23077+6707 to be a giant and gas-rich edge-on protoplanetary disk, which to our knowledge is the largest in extent so far discovered.

Length and Velocity Scales in Protoplanetary Disk Turbulence

In the theory of protoplanetary disk turbulence, a widely adopted ansatz, or assumption, is that the turnover frequency of the largest turbulent eddy, Ω L , is the local Keplerian frequency Ω K . In terms of the standard dimensionless Shakura–Sunyaev α parameter that quantifies turbulent viscosity or diffusivity, this assumption leads to characteristic length and velocity scales given respectively by αH and αc , in which H and c are the local gas scale height and sound speed. However, this assumption is not applicable in cases when turbulence is forced numerically or driven by some natural processes such as vertical shear instability. Here, we explore the more general case where Ω L ≥ Ω K and show that, under these conditions, the characteristic length and velocity scales are respectively α/R′H and αR′c , where R′≡ΩL/ΩK is twice the Rossby number. It follows that α=α˜/R′ , where α˜c is the root-mean-square average of the turbulent velocities. Properly allowing for this effect naturally explains the reduced particle scale heights produced in shearing box simulations of particles in forced turbulence, and it may help with interpreting recent edge-on disk observations; more general implications for observations are also presented. For R′>1 , the effective particle Stokes numbers are increased, which has implications for particle collision dynamics and growth, as well as for planetesimal formation.

JWST MIRI MRS Images of Disk Winds, Water, and CO in an Edge-on Protoplanetary Disk

We present JWST MIRI MRS observations of the edge-on protoplanetary disk around the young subsolar-mass star Tau 042021, acquired as part of the Cycle 1 GO program “Mapping Inclined Disk Astrochemical Signatures.” These data resolve the mid-IR spatial distributions of H2, revealing X-shaped emission extending to ∼200 au above the disk midplane with a semiopening angle of 35° ± 5°. We do not velocity-resolve the gas in the spectral images, but the measured semiopening angle of the H2 is consistent with a magnetohydrodynamic wind origin. A collimated, bipolar jet is seen in forbidden emission lines from [Ne ii], [Ne iii], [Ni ii], [Fe ii], [Ar ii], and [S iii]. Extended H2O and CO emission lines are also detected, reaching diameters of ∼90 and 190 au, respectively. Hot molecular emission is not expected at such radii, and we interpret its extended spatial distribution as scattering of inner disk molecular emission by dust grains in the outer disk surface. H i recombination lines, characteristic of inner disk accretion shocks, are similarly extended and are likely also scattered light from the innermost star–disk interface. Finally, we detect extended polycyclic aromatic hydrocarbon (PAH) emission at 11.3 μm cospatial with the scattered-light continuum, making this the first low-mass T Tauri star around which extended PAHs have been confirmed, to our knowledge. MIRI MRS line images of edge-on disks provide an unprecedented window into the outflow, accretion, and scattering processes within protoplanetary disks, allowing us to constrain the disk lifetimes and accretion and mass-loss mechanisms.

Deprojection and stellar dynamical modelling of boxy/peanut bars in edge-on discs

ABSTRACT We present a new method to infer the 3D luminosity distributions of edge-on barred galaxies with boxy-peanut/X (BP/X) shaped structures from their 2D surface brightness distributions. Our method relies on forward modelling of newly introduced parametric 3D density distributions for the BP/X bar, disc and other components using an existing image fitting software package (imfit). We validate our method using an N-body simulation of a barred disc galaxy with a moderately strong BP/X shape. For fixed orientation angles, the derived 3D BP/X-shaped density distribution is shown to yield a gravitational potential that is accurate to at least 5 per cent and forces that are accurate to at least 15 per cent, with average errors being $\sim 1.5~{{\ \rm per\ cent}}$ for both. When additional quantities of interest, such as the orientation of the bar to the line of sight, its pattern speed, and the stellar mass-to-light ratio are unknown they can be recovered to high accuracy by providing the parametric density distribution to the Schwarzschild modelling code FORSTAND. We also explore the ability of our models to recover the mass of the central supermassive black hole. This method is the first to be able to accurately recover both the orientation of the bar to the line of sight and its pattern speed when the disc is perfectly edge-on.

JWST Imaging of Edge-on Protoplanetary Disks. I. Fully Vertically Mixed 10 μm Grains in the Outer Regions of a 1000 au Disk

Scattered light imaging of protoplanetary disks provides key insights on the geometry and dust properties in the disk surface. Here, we present James Webb Space Telescope (JWST) 2–21 μm images of a 1000 au radius edge-on protoplanetary disk surrounding an 0.4 M ⊙ young star in Taurus, Two Micron All Sky Survey (2MASS) J04202144 + 2813491. These observations represent the longest wavelengths at which a protoplanetary disk is spatially resolved in scattered light. We combine these observations with Hubble Space Telescope optical images and Atacama Large Millimeter/submillimeter Array continuum and CO mapping. We find that the changes in the scattered light disk morphology are remarkably small across a factor of 30 in wavelength, indicating that dust in the disk surface layers is characterized by an almost gray opacity law. Using radiative transfer models, we conclude that grains up to ≳10 μm in size are fully coupled to the gas in this system, whereas grains ≳100 μm are strongly settled toward the midplane. Further analyses of these observations, and similar ones of other edge-on disks, will provide strong empirical constraints on disk dynamics and evolution and grain growth models. In addition, the 7.7 and 12. μm JWST images reveal an X-shaped feature located above the warm molecular layer traced by CO line emission. The highest elevations at which this feature is detectable roughly match the maximal extent of the disk in visible wavelength scattered light as well as of an unusual kinematic signature in CO. We propose that these phenomena could be related to a disk wind entraining small dust grains.

JWST Imaging of Edge-on Protoplanetary Disks. II. Appearance of Edge-on Disks with a Tilted Inner Region: Case Study of IRAS04302+2247

We present James Webb Space Telescope imaging from 2 to 21 μm of the edge-on protoplanetary disk around the embedded young star IRAS04302+2247. The structure of the source shows two reflection nebulae separated by a dark lane. The source extent is dominated by the extended filamentary envelope at ∼4.4 μm and shorter wavelengths, transitioning at 7.7 μm and longer wavelengths to more compact lobes of scattered light from the disk itself. The dark lane thickness does not vary significantly with wavelength, which we interpret as an indication for intermediate-sized (∼10 μm) grains in the upper layers of the disk. Intriguingly, we find that the brightest nebula of IRAS40302 switches side between 12.8 and 21 μm. We explore the effect of a tilted inner region on the general appearance of edge-on disks. We find that radiative transfer models of a disk including a tilted inner region can reproduce an inversion in the brightest nebula. In addition, for specific orientations, the model predicts strong lateral asymmetries, which can occur for more than half possible viewing azimuths. A large number of edge-on protoplanetary disks observed in scattered light show such lateral asymmetries (15/20), which suggests that a large fraction of protoplanetary disks might host a tilted inner region. Stellar spots may also induce lateral asymmetries, which are expected to vary over a significantly shorter timescale. Variability studies of edge-on disks would allow us to test the dominant scenario for the origin of these asymmetries.