Envelope-embedded protostellar disks represent the earliest stage of protoplanetary disk evolution, but their masses and sizes are difficult to measure because disk emission is entangled with the envelope. We analyze 26 protostellar disk systems in Taurus using ALMA Band 7 (345 GHz; ~0.3'') and VLA Ka-band (33 GHz; ~0.2'') continuum observations, together with molecular-line data to constrain disk gas masses. At 345 GHz, the median flux density, dust mass, and radius are 71 mJy, 5.5 M_Mearth, and 28 AU, with 68% ranges of 54-107 mJy, 3.9-9.4 M_Mearth, and 25-39 AU. At 33 GHz, the corresponding medians are 0.43 mJy, 39 M_Mearth, and 32 AU, with ranges of 0.41-0.80 mJy, 34-52 M_Mearth, and 29-33 AU. Taurus Class I disks are fainter and less massive than those in Orion, comparable to Perseus Class I disks but fainter than Perseus Class 0 disks, and brighter and more massive than those in Ophiuchus. Within Taurus, Class 0/I disks are brighter than Class II disks at both frequencies, although their inferred dust masses are comparable at 345 GHz and slightly higher at 33 GHz. Radiative-transfer modeling of CO isotopologue emission yields a median gas mass of 6.7 x 10^-4 M_Msun. The resulting CO-inferred gas-to-dust ratios span a broad range, with a mean of 147 +/- 75, a median of 26, and a 16th-84th percentile range of 8-147. This distribution overlaps the Taurus Class II population at the low end and ISM-like or higher values, including the AGE-PRO Ophiuchus Class 0/I population, at the high end.

We present new ALMA 0.9 mm and VLA 9 mm observations in the Taurus Molecular Cloud (TMC) of 25 protostellar systems, containing 40 protostars, observed at 0.3" ( 20 au) resolution. Within separations of 18-10,000 au, the ALMA/VLA-observed Taurus sample has a multiplicity fraction (MF), defined as the fraction of systems with at least one companion, of 0.50 +/- 0.07, and a companion fraction (CF), defined as the average number of companions per system, of 0.58 +/- 0.20. To build a more complete census of protostellar multiplicity in this region, we supplement the observed sample with 24 protostars (12 protostellar systems and 5 additional companions associated with systems we observed) previously identified through archival infrared or ALMA observations. Together, these 64 individual protostars (37 systems) define our Taurus+ sample, for which we measure higher values of 0.53 +/- 0.06 and 0.72 +/- 0.19 for the MF and CF, respectively. These multiplicity statistics in the TMC are notably higher than those reported in the more clustered star-forming regions of Orion and Perseus at the 3-4 sigma level, suggesting that Taurus may preserve a larger fraction of primordial multiples. The separation distributions in our samples show populations of both close and wide multiples, but a deficit at intermediate separations of 200-300 au. This pattern may suggest two distinct formation pathways: close binaries (<200 au) arising primarily from disk fragmentation, and wide multiples (>1000 au) from core fragmentation.

Kernel phase interferometry (KPI) is a post-processing technique that treats a conventional telescope as an interferometer by accurately modeling a telescope pupil as an array of virtual subapertures. KPI provides angular resolution within the diffraction limit by eliminating instrumental phase errors to first order. It has been successfully demonstrated to boost angular resolution on both space- and ground-based observatories, and is especially useful for enhancing space telescopes, as their diameters are smaller than the largest ground-based facilities. Here we present the first demonstration of KPI on JWST/MIRI data at 7.7 microns, 10 microns, and 15 microns. We generate contrast curves for 16 white dwarfs from the MIRI Exoplanets Orbiting White dwarfs (MEOW) Survey, finding significantly deeper contrast at small angular separations compared to traditional imaging with JWST/MIRI, down to within λ/D. Additionally, we use our KPI setup to successfully recover four known companions orbiting white dwarfs and brown dwarfs. This analysis shows that at these wavelengths KPI can uniquely access the orbital parameter space where inward-migrating post-main-sequence giant exoplanets are now thought to exist. We discuss the prospects for applying KPI to a larger sample of white dwarfs observed with JWST, increasing the volume of directly imaged close-in post-main-sequence exoplanets.

The HD 34282 transition disk is a favorable location for protoplanet searches given its complex substructures that may be driven by planet formation. Early millimeter continuum studies revealed a cavity with a stellar-centric radius of 0 .'' 24 and a ring that extends to 1 .'' 15, which implied the presence of a 50 MJ companion with a separation of similar to 0 .'' 1. We observed HD 34282 with Keck/NIRC2 in L ' using nonredundant masking interferometry in order to simultaneously search for protoplanets at small separations (down to similar to 0.5 lambda/D) and characterize the inner substructures of the transition disk. We fit geometric models to the data and find that a morphology with two disk components separated by a gap best describes the data. The inner component is comprised of scattering dust with a near-circular brightness distribution on the sky. The outer component can be modeled as a ring with a cavity that extends to similar to 0 .'' 26 in radius, which is about 11 degrees more inclined than measurements from past J-band observations of material near this region. We do not detect the 50 MJ companion but place an upper limit of similar to 40 MJ on companions at greater than or similar to lambda/D (0 .'' 08). We find that we are sensitive to a planet mass-times-accretion rate of similar to 1 x 10-6 MJ2 yr-1, or an accretion rate that could form a Jupiter analog in greater than or similar to 1 Myr.

Protoplanetary disks are the best places for observing planetary embryos. Direct imaging coupled with interferometric techniques, such as non-redundant masking (NRM), can help us better understand gas giant accretion timescales and dynamical interactions by resolving protoplanetary disks that exhibit evidence of planet formation. By using NRM we can achieve angular resolution down to and within the diffraction limit, and image planet formation on solar system scales (down to similar to 3-7 AU for K and L band, respectively) given the distances to most young stars (similar to 150 pc). We present progress on a NRM imaging survey designed to search for protoplanets embedded in protoplanetary disks. The goals of this survey are to detect and characterize protoplanets at solar system scales in a significant sample of protoplanetary disks and to characterize disk structure and dynamical interactions. From this survey, we can place constraints on the underlying protoplanet population and timescales under which giant gas planets form at spatial separations down to similar to 3-7 AU.

SCALES (Slicer Combined with Array of Lenslets for Exoplanet Spectroscopy) is a 2 - 5 micron high-contrast lenslet-based integral field spectrograph (IFS) designed to characterize exoplanets and their atmospheres. The SCALES medium-spectral-resolution mode uses a lenslet subarray with a 0.34 x 0.36 arcsecond field of view which allows for exoplanet characterization at increased spectral resolution. We explore the sensitivity limitations of this mode by simulating planet detections in the presence of realistic noise sources. We use the SCALES simulator scalessim to generate high-fidelity mock observations of planets that include speckle noise from their host stars, as well as other atmospheric and instrumental noise effects. We employ both angular and reference differential imaging as methods of disentangling speckle noise from the injected planet signals. These simulations allow us to assess the feasibility of speckle deconvolution for SCALES medium resolution data, and to test whether one approach outperforms another based on planet angular separations and contrasts.

We present Karl G. Jansky Very Large Array observations covering the NGC 1977 region at 3.0, 6.4, and 15.0 GHz. We search for compact radio sources and detect continuum emission from 34 NGC 1977 cluster members and 37 background objects. Of the 34 radio-detected cluster members, 3 are associated with known proplyds in NGC 1977, 22 are associated with additional young stellar objects in NGC 1977, and 9 are newly-identified cluster members. We examine the radio spectral energy distributions, circular polarization, and variability of the detected NGC 1977 sources, and identify 10 new candidate proplyds whose radio fluxes are dominated by optically thin free-free emission. We use measurements of free-free emission to calculate the mass-loss rates of known proplyds and new candidate proplyds in NGC 1977, and find values $\sim10^{-9}-10^{-8}$ M$_{\odot}$ yr$^{-1}$, which are lower than the mass-loss rates measured towards proplyds in the Orion Nebula Cluster, but consistent with the mass-loss rates predicted by external photoevaporation models for spatially-extended disks that are irradiated by the typical external UV fields encountered in NGC 1977. Finally, we show that photoevaporative disk winds in NGC 1977 may be illuminated by internal or external sources of ionization, depending on their positions within the cluster. This study provides new constraints on disk properties in a clustered star-forming region with a weaker UV environment than the Orion Nebula Cluster, but a stronger UV environment than low-mass star-forming regions like Taurus. Such intermediate UV environments represent the typical conditions of Galactic star and planet formation.

We present a direct imaging study of V892 Tau, a young Herbig Ae/Be star with a close-in stellar companion and circumbinary disk. Our observations consist of images acquired via Keck II/NIRC2 with nonredundant masking and the pyramid wavefront sensor at K ′ band (2.12 μ m) and L ′ band (3.78 μ m). Sensitivity to low-mass accreting companions and cool disk material is high at L ′ band, while complimentary observations at K ′ band probe hotter material with higher angular resolution. These multiwavelength, multiepoch data allow us to differentiate the secondary stellar emission from disk emission and deeply probe the structure of the circumbinary disk at small angular separations. We constrain architectural properties of the system by fitting geometric disk and companion models to the K ′ - and L ′ -band data. From these models, we constrain the astrometric and photometric properties of the stellar binary and update the orbit, placing the tightest estimates to date on the V892 Tau orbital parameters. We also constrain the geometric structure of the circumbinary disk, and resolve a circumprimary disk for the first time.

The Slicer Combined with an Array of Lenslets for Exoplanet Spectroscopy (SCALES) instrument is a lenslet-based integral field spectrograph that will operate at 2 to 5 microns, imaging and characterizing colder (and thus older) planets than current high-contrast instruments. Its spatial resolution for distant science targets and/or close-in disks and companions could be improved via interferometric techniques such as sparse aperture masking. We introduce a nascent Python package, NRM-artist, that we use to design several SCALES masks to be non-redundant and to have uniform coverage in Fourier space. We generate high-fidelity mock SCALES data using the scalessim package for SCALES’ low spectral resolution modes across its 2 to 5 micron bandpass. We include realistic noise from astrophysical and instrument sources, including Keck adaptive optics and Poisson noise. We inject planet and disk signals into the mock datasets and subsequently recover them to test the performance of SCALES sparse aperture masking and to determine the sensitivity of various mask designs to different science signals.

We present the highest-angular-resolution infrared monitoring of LkCa 15, a young solar analog hosting a transition disk. This system has been the subject of a number of direct-imaging studies from the millimeter through the optical, which have revealed multiple protoplanetary disk rings as well as three orbiting protoplanet candidates detected in infrared continuum emission (one of which was simultaneously seen at H α ). We use high-angular-resolution infrared imaging from 2014 to 2020 to systematically monitor these infrared signals and determine their physical origin. We find that three self-luminous protoplanets cannot explain the positional evolution of the infrared sources since the longer time baseline images lack the coherent orbital motion that would be expected for companions. However, the data still strongly prefer a time-variable morphology that cannot be reproduced by static scattered-light disk models. The multiepoch observations suggest the presence of complex and dynamic substructures moving through the forward-scattering side of the disk at ∼20 au or quickly varying shadowing by closer-in material. We explore whether the previous H α detection of one candidate would be inconsistent with this scenario and in the process develop an analytical signal-to-noise penalty for H α excesses detected near forward-scattered light. Under these new noise considerations, the H α detection is not strongly inconsistent with forward scattering, making the dynamic LkCa 15 disk a natural explanation for both the infrared and H α data.

Abstract As gas giant planets and brown dwarfs radiate away the residual heat from their formation, they cool through a spectral type transition from L to T, which encompasses the dissipation of cloud opacity and the appearance of strong methane absorption. While there are hundreds of known T-type brown dwarfs, the fi rst generation of directly-imaged exoplanets were all L-type. Recently, Kuzuhara et al. (2013) announced the discovery of GJ 504 b, the fi rst T dwarf exoplanet. GJ 504 b provides a unique opportunity to study the atmosphere of a new type of exoplanet with a ∼ 500 K temperature that bridges the gap between the fi rst directly imaged planets ( ∼ 1000 K) and our own Solar System’s Jupiter ( ∼ 130 K). We observed GJ 504 b in three narrow L-band fi lters (3.71, 3.88, and 4.00 µm), spanning the red end of the broad methane fundamental absorption feature (3.3 µm) as part of the LEECH exoplanet imaging survey. By compar-ing our new photometry and literature photometry to a grid of custom model atmospheres, we were able to fi t GJ 504 b’s unusual spectral energy distribution for the fi rst time. We fi nd that GJ 504 b is well-fi t by models with the follow-ing parameters: Teff =544±10 K, g/s2, [M/H]=0.60±0.12, cloud opacity parameter of fsed = 2 − 5, R=0.96±0.07 RJup, and log(L)=-6.13±0.03 L ⊙ , im-plying a hot start mass of 3-30 Mjup for a conservative age range of 0.1-6.5 Gyr. Of particular interest, our model fi ts suggest that GJ 504 b has a super-stellar metallicity. Since planet formation can create objects with non-stellar metal-licities, while binary star formation cannot, this result suggests that GJ 504 b formed like a planet, not like a binary companion.

The Orion Nebula Cluster (ONC) hosts protoplanetary disks experiencing external photoevaporation by the cluster’s intense UV field. These “proplyds” are comprised of a disk surrounded by an ionization front. We present ALMA Band 3 (3.1 mm) continuum observations of 12 proplyds. Thermal emission from the dust disks and free–free emission from the ionization fronts are both detected, and the high-resolution (0.″057) of the observations allows us to spatially isolate these two components. The morphology is unique compared to images at shorter (sub)millimeter wavelengths, which only detect the disks, and images at longer centimeter wavelengths, which only detect the ionization fronts. The disks are small ( r d = 6.4–38 au), likely due to truncation by ongoing photoevaporation. They have low spectral indices ( α ≲ 2.1) measured between Bands 7 and 3, suggesting the dust emission is optically thick. They harbor tens of Earth masses of dust as computed from the millimeter flux using the standard method although their true masses may be larger due to the high optical depth. We derive their photoevaporative mass-loss rates in two ways: first, by invoking ionization equilibrium and second, by using the brightness of the free–free emission to compute the density of the outflow. We find decent agreement between these measurements and M ̇ = 0.6–18.4 × 10 −7 M ⊙ yr −1 . The photoevaporation timescales are generally shorter than the ∼1 Myr age of the ONC, underscoring the known “proplyd lifetime problem.” Disk masses that are underestimated due to being optically thick remains one explanation to ease this discrepancy.

The Slicer Combined with Array of Lenslets for Exoplanet Spectroscopy (SCALES) is a 2 μm to 5 μm, high-contrast Integral Field Spectrograph (IFS) currently being built for Keck Observatory. With both low (R ≲ 250) and medium (R approximately 3500 to 7000) spectral resolution IFS modes, SCALES will detect and characterize significantly colder exoplanets than those accessible with near-infrared (approximately 1 μm to 2 μm) high-contrast spectrographs. This will lead to new progress in exoplanet atmospheric studies, including detailed characterization of benchmark systems that will advance the state of the art of atmospheric modeling. SCALES’ unique modes, while designed specifically for direct exoplanet characterization, will enable a broader range of novel (exo)planetary observations as well as galactic and extragalactic studies. Here we present the science cases that drive the design of SCALES. We describe an end-to-end instrument simulator that we use to track requirements and show simulations of expected science yields for each driving science case. We conclude with a discussion of preparations for early science when the instrument sees first light in approximately 2025.

The stellar cluster environment is expected to play a central role in the evolution of circumstellar disks. We use thermochemical modeling to constrain the dust and gas masses, disk sizes, UV and X-ray radiation fields, viewing geometries, and central stellar masses of 20 class II disks in the Orion Nebula Cluster (ONC). We fit a large grid of disk models to 350 GHz continuum, CO J = 3 − 2, and HCO ^+ J = 4 − 3 Atacama Large Millimeter/submillimeter Array observations of each target, and we introduce a procedure for modeling interferometric observations of gas disks detected in absorption against a bright molecular cloud background. We find that the ONC disks are massive and compact, with typical radii <100 au, gas masses ≥10 ^−3 M _⊙ , and gas-to-dust ratios ≥100. The interstellar‐medium‐like gas-to-dust ratios derived from our modeling suggest that compact, externally irradiated disks in the ONC are less prone to gas-phase CO depletion than the massive and extended gas disks that are commonly found in nearby low-mass star-forming regions. The presence of massive gas disks indicates that external photoevaporation may have only recently begun operating in the ONC; though it remains unclear whether other cluster members are older and more evaporated than the ones in our sample. Finally, we compare our dynamically derived stellar masses with the stellar masses predicted from evolutionary models and find excellent agreement. Our study has significantly increased the number of dynamical mass measurements in the mass range ≤0.5 M _⊙ , demonstrating that the ONC is an ideal region for obtaining large samples of dynamical mass measurements toward low-mass M-dwarfs.

Multiwavelength high-resolution imaging of protoplanetary disks has revealed the presence of multiple, varied substructures in their dust and gas components, which might be signposts of young, forming planetary systems. AB Aurigae bears an emblematic (pre)transitional disk showing spiral structures observed in the inner cavity of the disk in both the submillimeter (Atacama Large Millimeter/submillimeter Array (ALMA); 1.3 mm, 12CO) and near-infrared (Spectro-polarimetric High-contrast Exoplanet Research; 1.5–2.5 μm) wavelengths, which have been claimed to arise from dynamical interactions with a massive companion. In this work, we present new deep K s (2.16 μm) and L′ (3.7 μm) band images of AB Aurigae obtained with the L/M-band Infrared Camera on the Large Binocular Telescope, aimed for the detection of both planetary companions and extended disk structures. No point source is recovered, in particular at the outer regions of the disk, where a putative candidate (ρ = 0.″681, PA = 7.°6) had been previously claimed. The nature of a second innermost planet candidate (ρ = 0.″16, PA = 203.°9) cannot be investigated by the new data. We are able to derive 5σ detection limits in both magnitude and mass for the system, going from 14 M Jup at 0.″3 (49 au) down to 3–4 M Jup at 0.″6 (98 au) and beyond, based on the ATMO 2020 evolutionary models. We detect the inner spiral structures (<0.″5) resolved in both CO and polarimetric H-band observations. We also recover the ring structure of the system at larger separation (0.″5–0.″7) showing a clear southeast/northwest asymmetry. This structure, observed for the first time at L′ band, remains interior to the dust cavity seen at ALMA, suggesting an efficient dust trapping mechanism at play in the disk.

Betelgeuse has experienced a sudden shift in its brightness and dimmed mysteriously. This is likely caused by a hot blob of plasma ejected from Betelgeuse and then cooled to obscuring dust. If true, it is a remarkable opportunity to directly witness the formation of dust around a red supergiant star. Today's optical telescope facilities are not optimized for monitoring the Betelgeuse surface, so in this work, we propose a low-cost optical interferometer. The facility will consist of 12 x 4 inch optical telescopes mounted to the surface of a large radio dish for model-independent aperture synthesis imaging; polarization-maintaining single-mode fibers will carry the coherent beams from the individual optical telescopes to an all-in-one beam combiner. A fast steering mirror assisted fiber injection system guides the flux into fibers. A metrology system senses vibration-induced piston errors in optical fibers, and these errors are corrected using fast-steering delay lines. We will present the design.

We present Atacama Large Millimeter Array CO(3$-$2) and HCO$^+$(4$-$3) observations covering the central $1\rlap{.}'5$$\times$$1\rlap{.}'5$ region of the Orion Nebula Cluster (ONC). The unprecedented level of sensitivity ($\sim$0.1 mJy beam$^{-1}$) and angular resolution ($\sim$$0\rlap{.}''09 \approx 35$ AU) of these line observations enable us to search for gas-disk detections towards the known positions of submillimeter-detected dust disks in this region. We detect 23 disks in gas: 17 in CO(3$-$2), 17 in HCO$^+$(4$-$3), and 11 in both lines. Depending on where the sources are located in the ONC, we see the line detections in emission, in absorption against the warm background, or in both emission and absorption. We spectrally resolve the gas with $0.5$ km s$^{-1}$ channels, and find that the kinematics of most sources are consistent with Keplerian rotation. We measure the distribution of gas-disk sizes and find typical radii of $\sim$50-200 AU. As such, gas disks in the ONC are compact in comparison with the gas disks seen in low-density star-forming regions. Gas sizes are universally larger than the dust sizes. However, the gas and dust sizes are not strongly correlated. We find a positive correlation between gas size and distance from the massive star $\theta^1$ Ori C, indicating that disks in the ONC are influenced by photoionization. Finally, we use the observed kinematics of the detected gas lines to model Keplerian rotation and infer the masses of the central pre-main-sequence stars. Our dynamically-derived stellar masses are not consistent with the spectroscopically-derived masses, and we discuss possible reasons for this discrepancy.

Measuring the masses of protoplanetary disks is crucial for understanding their planet-forming potential. Typically, dust masses are derived from (sub-)millimeter flux density measurements plus assumptions for the opacity, temperature, and optical depth of the dust. Here we use radiative transfer models to quantify the validity of these assumptions with the aim of improving the accuracy of disk dust mass measurements. We first carry out a controlled exploration of disk parameter space. We find that the disk temperature is a strong function of disk size, while the optical depth depends on both disk size and dust mass. The millimeter-wavelength spectral index can be significantly shallower than the naive expectation due to a combination of optical depth and deviations from the Rayleigh-Jeans regime. We fit radiative transfer models to the spectral energy distributions (SEDs) of 132 disks in the Taurus-Auriga region using a Markov chain Monte Carlo approach. We used all available data to produce the most complete SEDs used in any extant modeling study. We perform the fitting twice: first with unconstrained disk sizes and again imposing the disk size-brightness relation inferred for sources in Taurus. This constraint generally forces the disks to be smaller, warmer, and more optically thick. From both sets of fits, we find disks to be similar to 1-5 times more massive than when derived using (sub-) millimeter measurements and common assumptions. With the uncertainties derived from our model fitting, the previously measured dust mass-stellar mass correlation is present in our study but only significant at the 2 sigma level.