Showing 1–50 of 245 results for author: Williams, J P

Disk Evolution Study Through Imaging of Nearby Young Stars (DESTINYS): PDS 111, an old T Tauri star with a young-looking disk

Authors: Annelotte Derkink, Christian Ginski, Paola Pinilla, Nicolas Kurtovic, Lex Kaper, Alex de Koter, Per-Gunnar Valegård, Eric Mamajek, Frank Backs, Myriam Benisty, Til Birnstiel, Gabriele Columba, Carsten Dominik, Antonio Garufi, Michiel Hogerheijde, Rob van Holstein, Jane Huang, François Ménard, Christian Rab, María Claudia Ramírez-Tannus, Álvaro Ribas, Jonathan P. Williams, Alice Zurlo

Abstract: The interplay between T Tauri stars and their circumstellar disks, and how this impacts the onset of planet formation has yet to be established. We studied a seemingly old T Tauri star, PDS 111, and its disk. We analyzed optical, infrared, and sub-millimeter observations obtained with VLT/X-shooter, Mercator/HERMES, TESS, VLT/SPHERE, and ALMA, providing a new view on PDS 111 and its protoplanetary… ▽ More The interplay between T Tauri stars and their circumstellar disks, and how this impacts the onset of planet formation has yet to be established. We studied a seemingly old T Tauri star, PDS 111, and its disk. We analyzed optical, infrared, and sub-millimeter observations obtained with VLT/X-shooter, Mercator/HERMES, TESS, VLT/SPHERE, and ALMA, providing a new view on PDS 111 and its protoplanetary disk. The multi-epoch spectroscopy yields photospheric lines to classify the star, and emission lines to study variability in the hot inner disk and to determine the mass-accretion rate. The SPHERE and ALMA observations are used to characterize the dust distribution of the small and large grains, respectively. PDS 111 is a weak-line T Tauri star with spectral type G2, exhibits strong H$α$ variability and with a low mass-accretion rate of $1-5\times10^{-10}$\,M$_{\odot}$\,yr$^{-1}$. We measured an age of the system of 15.9$^{+1.7}_{-3.7}$ Myr using pre-main sequence tracks. The SPHERE observations show a strongly flaring disk with an asymmetric substructure. The ALMA observations reveal a 30 au cavity in the dust continuum emission with a low contrast asymmetry in the South-West of the disk and a dust disk mass of 45.8\,$M_\oplus$. The $^{12}$CO radial extension is at least three times larger than that of the dust emission. Although the measured age is younger than suggested in literature, PDS 111 still seems relatively old; this provides insight into disk properties at an advanced stage of pre-main sequence evolution. The characteristics of this disk are very similar to its younger counterparts: strongly flaring, an average disk mass, a typical radial extent of the disk gas and dust, and the presence of common substructures. This suggests that disk evolution has not significantly changed the disk properties. These results show similarities with the "Peter Pan disks" around M-dwarfs. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Comments: 23 pages, 23 figures, accepted by A&A (abstract shortened)

Decision theory via model-free generalized fiducial inference

Authors: Jonathan P Williams, Yang Liu

Abstract: Building on the recent development of the model-free generalized fiducial (MFGF) paradigm (Williams, 2023) for predictive inference with finite-sample frequentist validity guarantees, in this paper, we develop an MFGF-based approach to decision theory. Beyond the utility of the new tools we contribute to the field of decision theory, our work establishes a formal connection between decision theori… ▽ More Building on the recent development of the model-free generalized fiducial (MFGF) paradigm (Williams, 2023) for predictive inference with finite-sample frequentist validity guarantees, in this paper, we develop an MFGF-based approach to decision theory. Beyond the utility of the new tools we contribute to the field of decision theory, our work establishes a formal connection between decision theories from the perspectives of fiducial inference, conformal prediction, and imprecise probability theory. In our paper, we establish pointwise and uniform consistency of an {\em MFGF upper risk function} as an approximation to the true risk function via the derivation of nonasymptotic concentration bounds, and our work serves as the foundation for future investigations of the properties of the MFGF upper risk from the perspective of new decision-theoretic, finite-sample validity criterion, as in Martin (2021). △ Less

Submitted 16 May, 2024; originally announced May 2024.

Early Planet Formation in Embedded Disks (eDisk) XV: Influence of Magnetic Field Morphology in Dense Cores on Sizes of Protostellar Disks

Authors: Hsi-Wei Yen, Jonathan P. Williams, Jinshi Sai, Patrick M. Koch, Ilseung Han, Jes K. Jørgensen, Woojin Kwon, Chang Won Lee, Zhi-Yun Li, Leslie W. Looney, Mayank Narang, Nagayoshi Ohashi, Shigehisa Takakuwa, John J. Tobin, Itziar de Gregorio-Monsalvo, Shih-Ping Lai, Jeong-Eun Lee, Kengo Tomida

Abstract: The magnetic field of a molecular cloud core may play a role in the formation of circumstellar disks in the core. We present magnetic field morphologies in protostellar cores of 16 targets in the Atacama Large Millimeter/submillimeter Array large program "Early Planet Formation in Embedded Disks (eDisk)", which resolved their disks with 7 au resolutions. The 0.1-pc scale magnetic field morphologie… ▽ More The magnetic field of a molecular cloud core may play a role in the formation of circumstellar disks in the core. We present magnetic field morphologies in protostellar cores of 16 targets in the Atacama Large Millimeter/submillimeter Array large program "Early Planet Formation in Embedded Disks (eDisk)", which resolved their disks with 7 au resolutions. The 0.1-pc scale magnetic field morphologies were inferred from the James Clerk Maxwell Telescope (JCMT) POL-2 observations. The mean orientations and angular dispersions of the magnetic fields in the dense cores are measured and compared with the radii of the 1.3 mm continuum disks and the dynamically determined protostellar masses from the eDisk program. We observe a significant correlation between the disk radii and the stellar masses. We do not find any statistically significant dependence of the disk radii on the projected misalignment angles between the rotational axes of the disks and the magnetic fields in the dense cores, nor on the angular dispersions of the magnetic fields within these cores. However, when considering the projection effect, we cannot rule out a positive correlation between disk radii and misalignment angles in three-dimensional space. Our results suggest that the morphologies of magnetic fields in dense cores do not play a dominant role in the disk formation process. Instead, the sizes of protostellar disks may be more strongly affected by the amount of mass that has been accreted onto star+disk systems, and possibly other parameters, for example, magnetic field strength, core rotation, and magnetic diffusivity. △ Less

Submitted 14 May, 2024; originally announced May 2024.

Comments: 14 pages, 4 figures, accepted by ApJ

SCExAO/CHARIS Multi-Wavelength, High-Contrast Imaging of the BD+45$^\circ$598 Debris Disk

Authors: Maria Vincent, Kellen Lawson, Thayne Currie, Jonathan P. Williams, Olivier Guyon, Julien Lozi, Vincent Deo, Sébastien Vievard

Abstract: We present a multi-wavelength (1.16$μ$m-2.37$μ$m) view of the debris disk around BD+45$^\circ$598, using the Subaru Coronagraphic Extreme Adaptive Optics system paired with the Coronagraphic High Angular Resolution Imaging Spectrograph. With an assumed age of 23 Myr, this source allows us to study the early evolution of debris disks and search for forming planets. We fit a scattered light model to… ▽ More We present a multi-wavelength (1.16$μ$m-2.37$μ$m) view of the debris disk around BD+45$^\circ$598, using the Subaru Coronagraphic Extreme Adaptive Optics system paired with the Coronagraphic High Angular Resolution Imaging Spectrograph. With an assumed age of 23 Myr, this source allows us to study the early evolution of debris disks and search for forming planets. We fit a scattered light model to our disk using a differential evolution algorithm, and constrain its geometry. We find the disk to have a peak density radius of $R_0 = 109.6$ au, an inclination of $i = 88.1^\circ$, and position angle $PA = 111.0^\circ$. While we do not detect a substellar companion in the disk, our calculated contrast limits indicate sensitivity to planets as small as $\sim 10 M_{\rm Jup}$ at a projected separation of 12 au of the star, and as small as $\sim 4 M_{\rm Jup}$ beyond 38 au. When measuring intensity as a function of wavelength, the disk color constrains the minimum dust grain size within a range of $\sim0.13$ to 1.01$μ$m. △ Less

Submitted 7 May, 2024; originally announced May 2024.

Comments: 22 pages, 12 figures; Accepted to AJ

Discovery of an accretion streamer and a slow wide-angle outflow around FU Orionis

Authors: A. S. Hales, A. Gupta, D. Ruiz-Rodriguez, J. P. Williams, S. Perez, L. Cieza, C. Gonzalez-Ruilova, J. E. Pineda, A. Santamaria-Miranda, J. Tobin, P. Weber, Z. Zhu, A. Zurlo

Abstract: We present ALMA 12-m, 7-m & Total Power (TP) Array observations of the FU Orionis outbursting system, covering spatial scales ranging from 160 to 25,000 au. The high-resolution interferometric data reveals an elongated $^{12}$CO(2-1) feature previously observed at lower resolution in $^{12}$CO(3-2). Kinematic modeling indicates that this feature can be interpreted as an accretion streamer feeding… ▽ More We present ALMA 12-m, 7-m & Total Power (TP) Array observations of the FU Orionis outbursting system, covering spatial scales ranging from 160 to 25,000 au. The high-resolution interferometric data reveals an elongated $^{12}$CO(2-1) feature previously observed at lower resolution in $^{12}$CO(3-2). Kinematic modeling indicates that this feature can be interpreted as an accretion streamer feeding the binary system. The mass infall rate provided by the streamer is significantly lower than the typical stellar accretion rates (even in quiescent states), suggesting that this streamer alone is not massive enough to sustain the enhanced accretion rates characteristic of the outbursting class prototype. The observed streamer may not be directly linked to the current outburst but rather a remnant of a previous, more massive streamer that may have contributed enough to the disk mass to render it unstable and trigger FU Ori's outburst. The new data detects, for the first time, a vast, slow-moving carbon monoxide molecular outflow emerging from this object. To accurately assess the outflow properties (mass, momentum, kinetic energy), we employed $^{13}$CO(2-1) data to correct for optical depth effects. The analysis indicates that the outflow corresponds to swept-up material not associated with the current outburst, similar to slow-molecular outflows observed around other FUor and Class I protostellar objects. △ Less

Submitted 5 May, 2024; originally announced May 2024.

Journal ref: 2024, ApJ, 966, 96H

Large-sample theory for inferential models: a possibilistic Bernstein--von Mises theorem

Authors: Ryan Martin, Jonathan P. Williams

Abstract: The inferential model (IM) framework offers alternatives to the familiar probabilistic (e.g., Bayesian and fiducial) uncertainty quantification in statistical inference. Allowing this uncertainty quantification to be imprecise makes it possible to achieve exact validity and reliability. But is imprecision and exact validity compatible with attainment of the classical notions of statistical efficie… ▽ More The inferential model (IM) framework offers alternatives to the familiar probabilistic (e.g., Bayesian and fiducial) uncertainty quantification in statistical inference. Allowing this uncertainty quantification to be imprecise makes it possible to achieve exact validity and reliability. But is imprecision and exact validity compatible with attainment of the classical notions of statistical efficiency? The present paper offers an affirmative answer to this question via a new possibilistic Bernstein--von Mises theorem that parallels a fundamental result in Bayesian inference. Among other things, our result demonstrates that the IM solution is asymptotically efficient in the sense that its asymptotic credal set is the smallest that contains the Gaussian distribution whose variance agrees with the Cramer--Rao lower bound. △ Less

Submitted 24 April, 2024; originally announced April 2024.

Comments: 8-page conference paper. Comments welcome at https://researchers.one/articles/24.04.00003

Anatomy of the Class I protostar L1489 IRS with NOEMA -- I. Disk, streamers, outflow(s) and bubbles at 3mm

Authors: M. Tanious, R. Le Gal, R. Neri, A. Faure, A. Gupta, C. J. Law, J. Huang, N. Cuello, J. P. Williams, F. Ménard

Abstract: Over the past few years, chemical studies have revealed multiple structures in the vicinity of young stellar objects (YSOs). It has become evident that specific physical conditions are associated with the emission of particular molecular lines, allowing us to use molecular probes of the YSO physics. Consequently, chemical surveys are now necessary to fully constrain the origin of the observed stru… ▽ More Over the past few years, chemical studies have revealed multiple structures in the vicinity of young stellar objects (YSOs). It has become evident that specific physical conditions are associated with the emission of particular molecular lines, allowing us to use molecular probes of the YSO physics. Consequently, chemical surveys are now necessary to fully constrain the origin of the observed structures. Several surveys have been conducted to explore the chemistry of YSOs, focusing on Class 0 and Class II objects. However, our knowledge of intermediate objects, that are Class I objects, remains limited. To bridge the gap and establish the relationship between observed structures and molecular line emission at the Class I evolutionary stage, we investigate the spatial distribution of key molecular gas species in the low-mass Class I protostar L1489 IRS (IRAS 04016+2610), a source part of the ChemYSO survey. We performed a 3mm line survey at high spatial and high spectral resolution using the NOEMA interferometer and the IRAM-30m telescope. We present here the ten brightest lines of our survey, in which we identified a new ~ 3 000 au long streamer in HC3N, C2H, and c-C3H2 emission, likely associated with more localized accretion shocks probed in SO. In addition, two ~ 10 000 au bubbles are seen with the dense molecular tracers HCO+, CS, and HCN around the YSO. Additionally, potential indicators of a second outflow appear in CS and HCN emission, but its nature remains to be confirmed. The late infall identified at large scales may originate from the nearby prestellar core L1489 and is likely responsible for the formation of an external warped disk in this system. The detection of a potential second outflow could be the direct evidence of a binary system. Finally, we hypothesize that the bubbles may result from the magnetic pressure as observed in numerical simulations. △ Less

Submitted 27 March, 2024; originally announced March 2024.

Comments: 16 pages, 12 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 687, A92 (2024)

A Bayesian shrinkage estimator for transfer learning

Authors: Mohamed A. Abba, Jonathan P. Williams, Brian J. Reich

Abstract: Transfer learning (TL) has emerged as a powerful tool to supplement data collected for a target task with data collected for a related source task. The Bayesian framework is natural for TL because information from the source data can be incorporated in the prior distribution for the target data analysis. In this paper, we propose and study Bayesian TL methods for the normal-means problem and multi… ▽ More Transfer learning (TL) has emerged as a powerful tool to supplement data collected for a target task with data collected for a related source task. The Bayesian framework is natural for TL because information from the source data can be incorporated in the prior distribution for the target data analysis. In this paper, we propose and study Bayesian TL methods for the normal-means problem and multiple linear regression. We propose two classes of prior distributions. The first class assumes the difference in the parameters for the source and target tasks is sparse, i.e., many parameters are shared across tasks. The second assumes that none of the parameters are shared across tasks, but the differences are bounded in $\ell_2$-norm. For the sparse case, we propose a Bayes shrinkage estimator with theoretical guarantees under mild assumptions. The proposed methodology is tested on synthetic data and outperforms state-of-the-art TL methods. We then use this method to fine-tune the last layer of a neural network model to predict the molecular gap property in a material science application. We report improved performance compared to classical fine tuning and methods using only the target data. △ Less

Submitted 5 June, 2024; v1 submitted 25 March, 2024; originally announced March 2024.

Direct-imaging Discovery of a Substellar Companion Orbiting the Accelerating Variable Star, HIP 39017

Authors: Taylor L. Tobin, Thayne Currie, Yiting Li, Jeffrey Chilcote, Timothy D. Brandt, Brianna Lacy, Masayuki Kuzuhara, Maria Vincent, Mona El Morsy, Vincent Deo, Jonathan P. Williams, Olivier Guyon, Julien Lozi, Sebastien Vievard, Nour Skaf, Kyohoon Ahn, Tyler Groff, N. Jeremy Kasdin, Taichi Uyama, Motohide Tamura, Aidan Gibbs, Briley L. Lewis, Rachel Bowens-Rubin, Maïssa Salama, Qier An , et al. (1 additional authors not shown)

Abstract: We present the direct-imaging discovery of a substellar companion (a massive planet or low-mass brown dwarf) to the young, $γ$ Doradus-type variable star, HIP 39017 (HD 65526). The companion's SCExAO/CHARIS JHK ($1.1-2.4μ$m) spectrum and Keck/NIRC2 L$^{\prime}$ photometry indicate that it is an L/T transition object. A comparison of the JHK+L$^{\prime}$ spectrum to several atmospheric model grids… ▽ More We present the direct-imaging discovery of a substellar companion (a massive planet or low-mass brown dwarf) to the young, $γ$ Doradus-type variable star, HIP 39017 (HD 65526). The companion's SCExAO/CHARIS JHK ($1.1-2.4μ$m) spectrum and Keck/NIRC2 L$^{\prime}$ photometry indicate that it is an L/T transition object. A comparison of the JHK+L$^{\prime}$ spectrum to several atmospheric model grids finds a significantly better fit to cloudy models than cloudless models. Orbit modeling with relative astrometry and precision stellar astrometry from Hipparcos and Gaia yields a semi-major axis of $23.8^{+8.7}_{-6.1}$ au, a dynamical companion mass of $30^{+31}_{-12}$~M$_J$, and a mass ratio of $\sim$1.9\%, properties most consistent with low-mass brown dwarfs. However, its mass estimated from luminosity models is a lower $\sim$13.8 $M_{\rm J}$ due to an estimated young age ($\lesssim$ 115 Myr); using a weighted posterior distribution informed by conservative mass constraints from luminosity evolutionary models yields a lower dynamical mass of $23.6_{-7.4}^{+9.1}$~M$_J$ and a mass ratio of $\sim$1.4\%. Analysis of the host star's multi-frequency $γ$ Dor-type pulsations, astrometric monitoring of HIP 39017b, and Gaia Data Release 4 astrometry of the star will clarify the system age and better constrain the mass and orbit of the companion. This discovery further reinforces the improved efficiency of targeted direct-imaging campaigns informed by long-baseline, precision stellar astrometry. △ Less

Submitted 15 May, 2024; v1 submitted 6 March, 2024; originally announced March 2024.

Comments: Published in AJ: April 9, 2024

Journal ref: 2024, AJ, 167, 205

Mass assembly in massive star formation: a fragmentation study of ATLASGAL clumps

Authors: Jagadheep D. Pandian, Rwitika Chatterjee, Timea Csengeri, Jonathan P. Williams, Friedrich Wyrowski, Karl M. Menten

Abstract: The mass assembly in star forming regions arises from the hierarchical structure in molecular clouds in tandem with fragmentation at different scales. In this paper, we present a study of the fragmentation of massive clumps covering a range of evolutionary states, selected from the ATLASGAL survey, using the compact configuration of the Submillimeter Array. The observations reveal a wide diversity… ▽ More The mass assembly in star forming regions arises from the hierarchical structure in molecular clouds in tandem with fragmentation at different scales. In this paper, we present a study of the fragmentation of massive clumps covering a range of evolutionary states, selected from the ATLASGAL survey, using the compact configuration of the Submillimeter Array. The observations reveal a wide diversity in the fragmentation properties with about 60% of the sources showing limited to no fragmentation at the 2" scale, or a physical scale of 0.015 - 0.09 pc. We also find several examples where the cores detected with the Submillimeter array are significantly offset from the clump potential suggesting that initial fragmentation does not result in the formation of a large number of Jeans mass fragments. The fraction of the clump mass that is in compact structures is seen to increase with source evolution. We also see a significant correlation between the maximum mass of a fragment and the bolometric luminosity of the parent clump. These suggest that massive star formation proceeds through clump fed core accretion with the initial fragmentation being dependent on the density structure of the clumps and/or magnetic fields. △ Less

Submitted 5 March, 2024; originally announced March 2024.

Comments: Accepted for publication in the Astrophysical Journal

The SPHERE view of the Orion star-forming region

Authors: P. -G. Valegard, C. Ginski, A. Derkink, A. Garufi, C. Dominik, A. Ribas, J. P. Williams, M. Benisty, T. Birnstiel, S. Facchini, G. Columba, M. Hogerheijde, R. G. Van Holstein, J. Huang, M. Kenworthy, C. F. Manara, P. Pinilla, Ch. Rab, R. Sulaiman, A. Zurlo

Abstract: We present SPHERE/IRDIS H-band data for a sample of 23 stars in the Orion Star forming region observed within the DESTINYS (Disk Evolution Study Through Imaging of Nearby Young Stars) program. We use polarization differential imaging in order to detect scattered light from circumstellar dust. From the scattered light observations we characterize the disk orientation, radius and contrast. We analys… ▽ More We present SPHERE/IRDIS H-band data for a sample of 23 stars in the Orion Star forming region observed within the DESTINYS (Disk Evolution Study Through Imaging of Nearby Young Stars) program. We use polarization differential imaging in order to detect scattered light from circumstellar dust. From the scattered light observations we characterize the disk orientation, radius and contrast. We analyse the disks in context of the stellar parameters and the environment of the Orion star-forming region. We use ancillary X-shooter spectroscopic observations to characterize the central stars in the systems. We furthermore use a combination of new and archival ALMA mm-continuum observations to characterize the dust masses present in the circumstellar disks. Within our sample we detect extended circumstellar disks in 10 of 23 systems. Of these, three are exceptionally extended (V351 Ori, V599 Ori and V1012 Ori) and show scattered light asymmetries which may indicate perturbations by embedded planets or (in the case of V599 Ori) by an outer stellar companion. Our high resolution imaging observations are also sensitive to close (sub)stellar companions and we detect 9 such objects in our sample of which 5 were previously unknown. We find in particular a possible sub-stellar companion (either a very low mass star or a high mass brown dwarf) 137 au from the star RY Ori. We find a strong anti-correlation between disk detection and multiplicity, with only 2 of our 10 disk detections located in stellar multiple systems. We also find a correlation between scattered light contrast and the millimetre flux suggesting that disks that have a high dust content are typically bright in near-infrared scattered light. Conversely we do not find significant correlations between scattered light contrast of the disks and the stellar mass or age. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Comments: Accepted by A&A

SMA detection of an extreme millimeter flare from the young class III star HD 283572

Authors: Joshua Bennett Lovell, Garrett K. Keating, David J. Wilner, Sean M. Andrews, Meredith MacGregor, Ramisa Akther Rahman, Ramprasad Rao, Jonathan P. Williams

Abstract: We present evidence of variable 1.3 millimeter emission from the 1-3 Myr, SpT G2-G5 class III YSO, HD~283572. HD~283572 was observed on 8 dates with the Submillimeter Array between 2021 December and 2023 May, a total on-source time of 10.2 hours, probing a range of timescales down to 5.2 seconds. Averaging all data obtained on 2022 Jan 17 shows a 4.4 mJy ($8.8σ$) point source detection with a nega… ▽ More We present evidence of variable 1.3 millimeter emission from the 1-3 Myr, SpT G2-G5 class III YSO, HD~283572. HD~283572 was observed on 8 dates with the Submillimeter Array between 2021 December and 2023 May, a total on-source time of 10.2 hours, probing a range of timescales down to 5.2 seconds. Averaging all data obtained on 2022 Jan 17 shows a 4.4 mJy ($8.8σ$) point source detection with a negative spectral index ($α{=}{-2.7}{\pm}1.2$), with peak emission rising to 13.8 mJy in one 3 minute span, and 25 mJy in one 29.7 second integration ($L_ν=4.7\times10^{17}$ erg s$^{-1}$ Hz$^{-1}$). Combining our data for the other 7 dates shows no detection, with an rms noise of 0.24 mJy beam$^{-1}$. The stochastic millimeter enhancements on time frames of seconds--minutes--hours with negative spectral indices are most plausibly explained by synchrotron or gyro-synchrotron radiation from stellar activity. HD 283572's 1.3 mm light-curve has similarities with variable binaries, suggesting HD 283572's activity may have been triggered by interactions with an as-yet undetected companion. We additionally identify variability of HD 283572 at 10 cm, from VLASS data. This study highlights the challenges of interpreting faint mm emission from evolved YSOs that may host tenuous disks, and suggests that a more detailed temporal analysis of spatially unresolved data is generally warranted. The variability of class III stars may open up new ground for understanding the physics of flares in the context of terrestrial planet formation. △ Less

Submitted 2 February, 2024; originally announced February 2024.

Comments: Accepted for publication in ApJL. 12 pages, inc. 5 figures, 2 tables

Anytime-Valid Generalized Universal Inference on Risk Minimizers

Authors: Neil Dey, Ryan Martin, Jonathan P. Williams

Abstract: A common goal in statistics and machine learning is estimation of unknowns. Point estimates alone are of little value without an accompanying measure of uncertainty, but traditional uncertainty quantification methods, such as confidence sets and p-values, often require strong distributional or structural assumptions that may not be justified in modern problems. The present paper considers a very c… ▽ More A common goal in statistics and machine learning is estimation of unknowns. Point estimates alone are of little value without an accompanying measure of uncertainty, but traditional uncertainty quantification methods, such as confidence sets and p-values, often require strong distributional or structural assumptions that may not be justified in modern problems. The present paper considers a very common case in machine learning, where the quantity of interest is the minimizer of a given risk (expected loss) function. For such cases, we propose a generalized universal procedure for inference on risk minimizers that features a finite-sample, frequentist validity property under mild distributional assumptions. One version of the proposed procedure is shown to be anytime-valid in the sense that it maintains validity properties regardless of the stopping rule used for the data collection process. We show how this anytime-validity property offers protection against certain factors contributing to the replication crisis in science. △ Less

Submitted 31 January, 2024; originally announced February 2024.

Comments: 33 pages, 6 figures

TIPSY: Trajectory of Infalling Particles in Streamers around Young stars. Dynamical analysis of the streamers around S CrA and HL Tau

Authors: Aashish Gupta, Anna Miotello, Jonathan P. Williams, Til Birnstiel, Michael Kuffmeier, Hsi-Wei Yen

Abstract: Context. Elongated trails of infalling gas, often referred to as "streamers," have recently been observed around young stellar objects (YSOs) at different evolutionary stages. This asymmetric infall of material can significantly alter star and planet formation processes, especially in the more evolved YSOs. Aims. In order to ascertain the infalling nature of observed streamer-like structures and… ▽ More Context. Elongated trails of infalling gas, often referred to as "streamers," have recently been observed around young stellar objects (YSOs) at different evolutionary stages. This asymmetric infall of material can significantly alter star and planet formation processes, especially in the more evolved YSOs. Aims. In order to ascertain the infalling nature of observed streamer-like structures and then systematically characterize their dynamics, we developed the code TIPSY (Trajectory of Infalling Particles in Streamers around Young stars). Methods. Using TIPSY, the streamer molecular line emission is first isolated from the disk emission. Then the streamer emission, which is effectively a point cloud in three-dimensional (3D) position-position-velocity space, is simplified to a curve-like representation. The observed streamer curve is then compared to the theoretical trajectories of infalling material. The best-fit trajectories are used to constrain streamer features, such as the specific energy, the specific angular momenta, the infall timescale, and the 3D morphology. Results. We used TIPSY to fit molecular-line ALMA observations of streamers around a Class II binary system, S CrA, and a Class I/II protostar, HL Tau. Our results indicate that both of the streamers are consistent with infalling motion. TIPSY results and mass estimates suggest that S CrA and HL Tau are accreting material at a rate of $\gtrsim27$ M$_{jupiter}$ Myr$^{-1}$ and $\gtrsim5$ M$_{jupiter}$ Myr$^{-1}$, respectively, which can significantly increase the mass budget available to form planets. Conclusions. TIPSY can be used to assess whether the morphology and kinematics of observed streamers are consistent with infalling motion and to characterize their dynamics, which is crucial for quantifying their impact on the protostellar systems. △ Less

Submitted 18 January, 2024; originally announced January 2024.

Comments: Accepted in Astronomy & Astrophysics

Early Planet Formation in Embedded Disks (eDisk) XIV: Flared Dust Distribution and Viscous Accretion Heating of the Disk around R CrA IRS 7B-a

Authors: Shigehisa Takakuwa, Kazuya Saigo, Miyu Kido, Nagayoshi Ohashi, John J. Tobin, Jes K. Jørgensen, Yuri Aikawa, Yusuke Aso, Sacha Gavino, Ilseung Han, Patrick M. Koch, Woojin Kwon, Chang Won Lee, Jeong-Eun Lee, Zhi-Yun Li, Zhe-Yu Daniel Lin, Leslie W. Looney, Shoji Mori, Jinshi Sai, Rajeeb Sharma, Patrick Sheehan, Kengo Tomida, Jonathan P. Williams, Yoshihide Yamato, Hsi-Wei Yen

Abstract: We performed radiative transfer calculations and observing simulations to reproduce the 1.3-mm dust-continuum and C$^{18}$O (2-1) images in the Class I protostar R CrA IRS7B-a, observed with the ALMA Large Program ``Early Planet Formation in Embedded Disks (eDisk)". We found that the dust disk model passively heated by the central protostar cannot reproduce the observed peak brightness temperature… ▽ More We performed radiative transfer calculations and observing simulations to reproduce the 1.3-mm dust-continuum and C$^{18}$O (2-1) images in the Class I protostar R CrA IRS7B-a, observed with the ALMA Large Program ``Early Planet Formation in Embedded Disks (eDisk)". We found that the dust disk model passively heated by the central protostar cannot reproduce the observed peak brightness temperature of the 1.3-mm continuum emission ($\sim$195 K), regardless of the assumptions about the dust opacity. Our calculation suggests that viscous accretion heating in the disk is required to reproduce the observed high brightness temperature. The observed intensity profile of the 1.3-mm dust-continuum emission along the disk minor axis is skewed toward the disk far side. Our modeling reveals that such an asymmetric intensity distribution requires flaring of the dust along the disk's vertical direction with the scale-height following $h/r \sim r^{0.3}$ as function of radius. These results are in sharp contrast to those of Class II disks, which show geometrically flat dust distributions and lower dust temperatures. From our modeling of the C$^{18}$O (2-1) emission, the outermost radius of the gas disk is estimated to be $\sim$80 au, larger than that of the dust disk ($\sim$62 au), to reproduce the observed distribution of the C$^{18}$O (2-1) emission in IRS 7B-a. Our modeling unveils a hot and thick dust disk plus a larger gas disk around one of the eDisk targets, which could be applicable to other protostellar sources in contrast to more evolved sources. △ Less

Submitted 16 January, 2024; originally announced January 2024.

Comments: 26 pages, 13 figures

The Kinematic and Dynamic Properties of HBC 494's Wide-Angle Outflows

Authors: Austen Fourkas, Dary Ruiz-Rodriguez, Lee G. Mundy, Jonathan P. Williams

Abstract: We present Atacama Large Millimeter/sub-millimeter Array (ALMA) Cycle-5 observations of HBC 494, as well as calculations of the kinematic and dynamic variables which represent the object's wide-angle bipolar outflows. HBC 494 is a binary FU Orionis type object located in the Orion A molecular cloud. We take advantage of combining the ALMA main array, Atacama Compact Array (ACA), and Total Power (T… ▽ More We present Atacama Large Millimeter/sub-millimeter Array (ALMA) Cycle-5 observations of HBC 494, as well as calculations of the kinematic and dynamic variables which represent the object's wide-angle bipolar outflows. HBC 494 is a binary FU Orionis type object located in the Orion A molecular cloud. We take advantage of combining the ALMA main array, Atacama Compact Array (ACA), and Total Power (TP) array in order to map HBC 494's outflows and thus, estimate their kinematic parameters with higher accuracy in comparison to prior publications. We use $^{12}$CO, $^{13}$CO, C$^{18}$O and SO observations to describe the object's outflows, envelope, and disc, as well as estimate the mass, momentum, and kinetic energy of the outflows. After correcting for optical opacity near systemic velocities, we estimate a mass of $3.0\times10^{-2}$ M$_{\odot}$ for the southern outflow and $2.8\times10^{-2}$ M$_{\odot}$ for the northern outflow. We report the first detection of a secondary outflow cavity located approximately $15$" north of the central binary system, which could be a remnant of a previous large-scale accretion outburst. Furthermore, we find CO spatial features in HBC 494's outflows corresponding to position angles of $\sim35^{\circ}$ and $\sim145^{\circ}$. This suggests that HBC 494's outflows are most likely a composite of overlapping outflows from two different sources, i.e., HBC 494a and HBC 494b, the two objects in the binary system. △ Less

Submitted 24 January, 2024; v1 submitted 6 January, 2024; originally announced January 2024.

Comments: 26 pages, 24 figures, accepted to MNRAS

Uncertainty quantification in automated valuation models with locally weighted conformal prediction

Authors: Anders Hjort, Gudmund Horn Hermansen, Johan Pensar, Jonathan P. Williams

Abstract: Non-parametric machine learning models, such as random forests and gradient boosted trees, are frequently used to estimate house prices due to their predictive accuracy, but such methods are often limited in their ability to quantify prediction uncertainty. Conformal Prediction (CP) is a model-agnostic framework for constructing confidence sets around machine learning prediction models with minima… ▽ More Non-parametric machine learning models, such as random forests and gradient boosted trees, are frequently used to estimate house prices due to their predictive accuracy, but such methods are often limited in their ability to quantify prediction uncertainty. Conformal Prediction (CP) is a model-agnostic framework for constructing confidence sets around machine learning prediction models with minimal assumptions. However, due to the spatial dependencies observed in house prices, direct application of CP leads to confidence sets that are not calibrated everywhere, i.e., too large of confidence sets in certain geographical regions and too small in others. We survey various approaches to adjust the CP confidence set to account for this and demonstrate their performance on a data set from the housing market in Oslo, Norway. Our findings indicate that calibrating the confidence sets on a \textit{locally weighted} version of the non-conformity scores makes the coverage more consistently calibrated in different geographical regions. We also perform a simulation study on synthetically generated sale prices to empirically explore the performance of CP on housing market data under idealized conditions with known data-generating mechanisms. △ Less

Submitted 11 December, 2023; originally announced December 2023.

Disk Evolution Study Through Imaging of Nearby Young Stars (DESTINYS): HD 34700 A unveils an inner ring

Authors: G. Columba, E. Rigliaco, R. Gratton, D. Mesa, V. D'Orazi, C. Ginski, N. Engler, J. P. Williams, J. Bae, M. Benisty, T. Birnstiel, P. Delorme, C. Dominik, S. Facchini, F. Menard, P. Pinilla, C. Rab, Á. Ribas, V. Squicciarini, R. G. van Holstein, A. Zurlo

Abstract: Context. The study of protoplanetary disks is fundamental to understand their evolution and interaction with the surrounding environment, and to constrain planet formation mechanisms. Aims. We aim at characterising the young binary system HD 34700 A, which shows a wealth of structures. Methods. Taking advantage of the high-contrast imaging instruments SPHERE at the VLT, LMIRCam at the LBT, and… ▽ More Context. The study of protoplanetary disks is fundamental to understand their evolution and interaction with the surrounding environment, and to constrain planet formation mechanisms. Aims. We aim at characterising the young binary system HD 34700 A, which shows a wealth of structures. Methods. Taking advantage of the high-contrast imaging instruments SPHERE at the VLT, LMIRCam at the LBT, and of ALMA observations, we analyse this system at multiple wavelengths. We study the rings and spiral arms morphology and the scattering properties of the dust. We discuss the possible causes of all the observed features. Results. We detect for the first time, in the H$α$ band, a ring extending from $\sim$65 au to ${\sim}$120 au, inside the ring already known from recent studies. These two have different physical and geometrical properties. Based on the scattering properties, the outer ring may consist of grains of typical size $a_{out} > 4 μm$, while the inner ring of smaller grains ($a_{in} <= 0.4 {μm}$). Two extended logarithmic spiral arms stem from opposite sides of the disk. The outer ring appears as a spiral arm itself, with a variable radial distance from the centre and extended substructures. ALMA data confirm the presence of a millimetric dust substructure centred just outside the outer ring, and detect misaligned gas rotation patterns for HD 34700 A and B. Conclusions. The complexity of HD 34700 A, revealed by the variety of observed features, suggests the existence of one or more disk-shaping physical mechanisms. Possible scenarios, compatible with our findings, involve the presence inside the disk of a yet undetected planet of several Jupiter masses and the system interaction with the surroundings by means of gas cloudlet capture or flybys. Further observations with JWST/MIRI or ALMA (gas kinematics) could shed more light on these. △ Less

Submitted 25 October, 2023; originally announced October 2023.

Comments: Accepted for publication on A&A. 14 + 5 pages, 9 + 7 figures (text + appendix)

Journal ref: A&A 681, A19 (2024)

Early Planet Formation in Embedded Disks (eDisk) X: Compact Disks, Extended Infall, and a Fossil Outburst in the Class I Oph IRS43 Binary

Authors: Suchitra Narayanan, Jonathan P. Williams, John J. Tobin, Jes K. Jorgensen, Nagayoshi Ohashi, Zhe-Yu Daniel Lin, Merel L. R. van't Hoff, Zhi-Yun Li, Adele L. Plunkett, Leslie W. Looney, Shigehisa Takakuwa, Hsi-Wei Yen, Yusuke Aso, Christian Flores, Jeong-Eun Lee, Shih-Ping Lai, Woojin Kwon, Itziar de Gregorio-Monsalvo, Rajeeb Sharma, Chang Won Lee

Abstract: We present the first results from the Early Planet Formation in Embedded Disks (eDisk) ALMA Large Program toward Oph IRS43, a binary system of solar mass protostars. The 1.3 mm dust continuum observations resolve a compact disk, ~6au radius, around the northern component and show that the disk around the southern component is even smaller, <~3 au. CO, 13CO, and C18O maps reveal a large cavity in a… ▽ More We present the first results from the Early Planet Formation in Embedded Disks (eDisk) ALMA Large Program toward Oph IRS43, a binary system of solar mass protostars. The 1.3 mm dust continuum observations resolve a compact disk, ~6au radius, around the northern component and show that the disk around the southern component is even smaller, <~3 au. CO, 13CO, and C18O maps reveal a large cavity in a low mass envelope that shows kinematic signatures of rotation and infall extending out to ~ 2000au. An expanding CO bubble centered on the extrapolated location of the source ~130 years ago suggests a recent outburst. Despite the small size of the disks, the overall picture is of a remarkably large and dynamically active region. △ Less

Submitted 23 October, 2023; originally announced October 2023.

Comments: Paper 10 of the ALMA eDisk Large Program. Accepted for publication in ApJ

Early Planet Formation in Embedded Disks (eDisk) XII: Accretion streamers, protoplanetary disk, and outflow in the Class I source Oph IRS63

Authors: Christian Flores, Nagayoshi Ohashi, John J. Tobin, Jes K. Jørgensen, Shigehisa Takakuwa, Zhi-Yun Li, Zhe-Yu Daniel Lin, Merel L. R. van 't Hoff, Adele L. Plunkett, Yoshihide Yamato, Jinshi Sai, Patrick M. Koch, Hsi-Wei Yen, Yuri Aikawa, Yusuke Aso, Itziar de Gregorio-Monsalvo, Miyu Kido, Woojin Kwon, Jeong-Eun Lee, Chang Won Lee, Leslie W. Looney, Alejandro Santamaría-Miranda, Rajeeb Sharma, Travis J. Thieme, Jonathan P. Williams , et al. (3 additional authors not shown)

Abstract: We present ALMA observations of the Class I source Oph IRS63 in the context of the Early Planet Formation in Embedded Disks (eDisk) large program. Our ALMA observations of Oph IRS63 show a myriad of protostellar features, such as a shell-like bipolar outflow (in $^{12}$CO), an extended rotating envelope structure (in $^{13}$CO), a streamer connecting the envelope to the disk (in C$^{18}$O), and se… ▽ More We present ALMA observations of the Class I source Oph IRS63 in the context of the Early Planet Formation in Embedded Disks (eDisk) large program. Our ALMA observations of Oph IRS63 show a myriad of protostellar features, such as a shell-like bipolar outflow (in $^{12}$CO), an extended rotating envelope structure (in $^{13}$CO), a streamer connecting the envelope to the disk (in C$^{18}$O), and several small-scale spiral structures seen towards the edge of the dust continuum (in SO). By analyzing the velocity pattern of $^{13}$CO and C$^{18}$O, we measure a protostellar mass of $\rm M_\star = 0.5 \pm 0.2 $~$\rm M_\odot$ and confirm the presence of a disk rotating at almost Keplerian velocity that extends up to $\sim260$ au. These calculations also show that the gaseous disk is about four times larger than the dust disk, which could indicate dust evolution and radial drift. Furthermore, we model the C$^{18}$O streamer and SO spiral structures as features originating from an infalling rotating structure that continuously feeds the young protostellar disk. We compute an envelope-to-disk mass infall rate of $\sim 10^{-6}$~$\rm M_\odot \, yr^{-1}$ and compare it to the disk-to-star mass accretion rate of $\sim 10^{-8}$~$\rm M_\odot \, yr^{-1}$, from which we infer that the protostellar disk is in a mass build-up phase. At the current mass infall rate, we speculate that soon the disk will become too massive to be gravitationally stable. △ Less

Submitted 23 October, 2023; originally announced October 2023.

Comments: 26 pages and 17 figures

Early Planet Formation in Embedded Disks (eDisk). VIII. A Small Protostellar Disk around the Extremely Low-Mass and Young Class 0 Protostar, IRAS 15398-3359

Authors: Travis J. Thieme, Shih-Ping Lai, Nagayoshi Ohashi, John J. Tobin, Jes K. Jørgensen, Jinshi Sai, Yusuke Aso, Jonathan P. Williams, Yoshihide Yamato, Yuri Aikawa, Itziar de Gregorio-Monsalvo, Ilseung Han, Woojin Kwon, Chang Won Lee, Jeong-Eun Lee, Zhi-Yun Li, Zhe-Yu Daniel Lin, Leslie W. Looney, Suchitra Narayanan, Nguyen Thi Phuong, Adele L. Plunkett, Alejandro Santamaría-Miranda, Rajeeb Sharma, Shigehisa Takakuwa, Hsi-Wei Yen

Abstract: Protostellar disks are a ubiquitous part of the star formation process and the future sites of planet formation. As part of the Early Planet Formation in Embedded Disks (eDisk) large program, we present high-angular resolution dust continuum ($\sim40\,$mas) and molecular line ($\sim150\,$mas) observations of the Class 0 protostar, IRAS 15398-3359. The dust continuum is small, compact, and centrall… ▽ More Protostellar disks are a ubiquitous part of the star formation process and the future sites of planet formation. As part of the Early Planet Formation in Embedded Disks (eDisk) large program, we present high-angular resolution dust continuum ($\sim40\,$mas) and molecular line ($\sim150\,$mas) observations of the Class 0 protostar, IRAS 15398-3359. The dust continuum is small, compact, and centrally peaked, while more extended dust structures are found in the outflow directions. We perform a 2D Gaussian fitting to find the deconvolved size and $2σ$ radius of the dust disk to be $4.5\times2.8\,\mathrm{au}$ and $3.8\,\mathrm{au}$, respectively. We estimate the gas+dust disk mass assuming optically thin continuum emission to be $0.6-1.8\,M_\mathrm{jup}$, indicating a very low-mass disk. The CO isotopologues trace components of the outflows and inner envelope, while SO traces a compact, rotating disk-like component. Using several rotation curve fittings on the PV diagram of the SO emission, the lower limits of the protostellar mass and gas disk radius are $0.022\,M_\odot$ and $31.2\,\mathrm{au}$ from our Modified 2 single power-law fitting. A conservative upper limit of the protostellar mass is inferred to be $0.1\,M_\odot$. The protostellar mass-accretion rate and the specific angular momentum at the protostellar disk edge are found to be between $1.3-6.1\times10^{-6}\,M_\odot\,\mathrm{yr^{-1}}$ and $1.2-3.8\times10^{-4}\,\mathrm{km\,s^{-1}\,pc}$, respectively, with an age estimated between $0.4-7.5\times10^{4}\,$yr. At this young age with no clear substructures in the disk, planet formation would likely not yet have started. This study highlights the importance of high-resolution observations and systematic fitting procedures when deriving dynamical properties of deeply embedded Class 0 protostars. △ Less

Submitted 19 October, 2023; originally announced October 2023.

Comments: 28 pages, 16 figures. Accepted for publication in ApJ as one of the first-look papers of the eDisk ALMA Large Program

Bayes Watch: Bayesian Change-point Detection for Process Monitoring with Fault Detection

Authors: Alexander C. Murph, Curtis B. Storlie, Patrick M. Wilson, Jonathan P. Williams, Jan Hannig

Abstract: When a predictive model is in production, it must be monitored in real-time to ensure that its performance does not suffer due to drift or abrupt changes to data. Ideally, this is done long before learning that the performance of the model itself has dropped by monitoring outcome data. In this paper we consider the problem of monitoring a predictive model that identifies the need for palliative ca… ▽ More When a predictive model is in production, it must be monitored in real-time to ensure that its performance does not suffer due to drift or abrupt changes to data. Ideally, this is done long before learning that the performance of the model itself has dropped by monitoring outcome data. In this paper we consider the problem of monitoring a predictive model that identifies the need for palliative care currently in production at the Mayo Clinic in Rochester, MN. We introduce a framework, called \textit{Bayes Watch}, for detecting change-points in high-dimensional longitudinal data with mixed variable types and missing values and for determining in which variables the change-point occurred. Bayes Watch fits an array of Gaussian Graphical Mixture Models to groupings of homogeneous data in time, called regimes, which are modeled as the observed states of a Markov process with unknown transition probabilities. In doing so, Bayes Watch defines a posterior distribution on a vector of regime assignments, which gives meaningful expressions on the probability of every possible change-point. Bayes Watch also allows for an effective and efficient fault detection system that assesses what features in the data where the most responsible for a given change-point. △ Less

Submitted 4 October, 2023; originally announced October 2023.

Comments: 34 pages without Supplementary Material, 5 figures, 1 table

Testing external photoevaporation in the $σ$-Orionis cluster with spectroscopy and disk mass measurements

Authors: K. Maucó, C. F. Manara, M. Ansdell, G. Bettoni, R. Claes, J. Alcala, A. Miotello, S. Facchini, T. J. Haworth, G. Lodato, J. P. Williams

Abstract: The evolution of protoplanetary disks is regulated by an interplay of several processes, either internal to the system or related to the environment. As most of the stars and planets have formed in massive stellar clusters, studying the effects of UV radiation on disk evolution is of paramount importance. Here we test the impact of external photoevaporation on the evolution of disks in the $σ$ Ori… ▽ More The evolution of protoplanetary disks is regulated by an interplay of several processes, either internal to the system or related to the environment. As most of the stars and planets have formed in massive stellar clusters, studying the effects of UV radiation on disk evolution is of paramount importance. Here we test the impact of external photoevaporation on the evolution of disks in the $σ$ Orionis cluster by conducting the first combined large-scale UV to IR spectroscopic and mm-continuum survey of this region. We study a sample of 50 targets located at increasing distances from the central, OB system $σ$ Ori. We combine new VLT/X-Shooter spectra with new and previously published ALMA measurements of disk dust and gas fluxes and masses. We confirm the previously found decrease of $M_{\rm dust}$ in the inner $\sim$0.5 pc of the cluster. This is particularly evident when considering the disks around the more massive stars ($\ge$ 0.4 $M_{\odot}$), where those located in the inner part ($<$ 0.5 pc) have $M_{\rm dust}$ about an order of magnitude lower than the more distant ones. About half of the sample is located in the region of the $\dot{M}_{\rm acc}$ vs $M_{\rm disk}$ expected by models of external photoevaporation, namely showing shorter disk lifetimes. These are observed for all targets with projected separation from $σ$ Ori $<$ 0.5 pc, proving that the presence of a massive stellar system affects disk evolution. External photoevaporation is a viable mechanism to explain the observed shorter disk lifetimes and lower $M_{\rm dust}$ in the inner $\sim$0.5 pc of the cluster. Follow-up observations of the low stellar mass targets are crucial to confirm the dependence of the external photoevaporation process with stellar host mass. This work confirms that the effects of external photoevaporation are significant down to impinging radiation as low as $\sim 10^{4}$ G$_0$. △ Less

Submitted 11 September, 2023; originally announced September 2023.

Comments: Accepted for publication on Astronomy & Astrophysics. 13 pages, 7 figures + appendix. Abstract abridged to meet arXiv requirements

Journal ref: A&A 679, A82 (2023)

Early Planet Formation in Embedded Disks (eDisk) VI: Kinematic Structures around the Very Low Mass Protostar IRAS 16253-2429

Authors: Yusuke Aso, Woojin Kwon, Nagayoshi Ohashi, Jes K. Jorgensen, John J. Tobin, Yuri Aikawa, Itziar de Gregorio-Monsalvo, Ilseung Han, Miyu Kido, Patrick M. Koch, Shih-Ping Lai, Chang Won Lee, Jeong-Eun Lee, Zhi-Yun Li, Zhe-Yu Daniel Lin, Leslie W. Looney, Suchitra Narayanan, Nguyen Thi Phuong, Jinshi Sai, Kazuya Saigo, Alejandro Santamaria-Miranda, Rajeeb Sharma, Shigehisa Takakuwa, Travis J. Thieme, Kengo Tomida , et al. (2 additional authors not shown)

Abstract: Precise estimates of protostellar masses are crucial to characterize the formation of stars of low masses down to brown-dwarfs (BDs; M* < 0.08 Msun). The most accurate estimation of protostellar mass uses the Keplerian rotation in the circumstellar disk around the protostar. To apply the Keplerian rotation method to a protostar at the low-mass end, we have observed the Class 0 protostar IRAS 16253… ▽ More Precise estimates of protostellar masses are crucial to characterize the formation of stars of low masses down to brown-dwarfs (BDs; M* < 0.08 Msun). The most accurate estimation of protostellar mass uses the Keplerian rotation in the circumstellar disk around the protostar. To apply the Keplerian rotation method to a protostar at the low-mass end, we have observed the Class 0 protostar IRAS 16253-2429 using the Atacama Large Millimeter/submillimeter Array (ALMA) in the 1.3 mm continuum at an angular resolution of 0.07" (10 au), and in the 12CO, C18O, 13CO (J=2-1), and SO (J_N = 6_5-5_4) molecular lines, as part of the ALMA Large Program Early Planet Formation in Embedded Disks (eDisk). The continuum emission traces a non-axisymmetric, disk-like structure perpendicular to the associated 12CO outflow. The position-velocity (PV) diagrams in the C18O and 13CO lines can be interpreted as infalling and rotating motions. In contrast, the PV diagram along the major axis of the disk-like structure in the 12CO line allows us to identify Keplerian rotation. The central stellar mass and the disk radius are estimated to be ~0.12-0.17 Msun and ~13-19 au, respectively. The SO line suggests the existence of an accretion shock at a ring (r~28 au) surrounding the disk and a streamer from the eastern side of the envelope. IRAS 16253-2429 is not a proto-BD but has a central stellar mass close to the BD mass regime, and our results provide a typical picture of such very low-mass protostars. △ Less

Submitted 4 September, 2023; originally announced September 2023.

Comments: 41 pages, 14 figures

Journal ref: ApJ 954 (2023) 101-117

Early Planet Formation in Embedded Disks (eDisk) IX: High-resolution ALMA Observations of the Class 0 Protostar R CrA IRS5N and its surrounding

Authors: Rajeeb Sharma, Jes K. Jørgensen, Sacha Gavino, Nagayoshi Ohashi, John J. Tobin, Zhe-Yu Daniel Lin, Zhi-Yun Li, Shigehisa Takakuwa, Chang Won Lee, Jinshi Sai, Woojin Kwon, Itziar de Gregorio-Monsalvo, Alejandro Santamaría-Miranda, Hsi-Wei Yen, Yuri Aikawa, Yusuke Aso, Shih-Ping Lai, Jeong-Eun Lee, Leslie W. Looney, Nguyen Thi Phuong, Travis J. Thieme, Jonathan P. Williams

Abstract: We present high-resolution, high-sensitivity observations of the Class 0 protostar RCrA IRS5N as part of the Atacama Large Milimeter/submilimeter Array (ALMA) large program Early Planet Formation in Embedded Disks (eDisk). The 1.3 mm continuum emission reveals a flattened continuum structure around IRS5N, consistent with a protostellar disk in the early phases of evolution. The continuum emission… ▽ More We present high-resolution, high-sensitivity observations of the Class 0 protostar RCrA IRS5N as part of the Atacama Large Milimeter/submilimeter Array (ALMA) large program Early Planet Formation in Embedded Disks (eDisk). The 1.3 mm continuum emission reveals a flattened continuum structure around IRS5N, consistent with a protostellar disk in the early phases of evolution. The continuum emission appears smooth and shows no substructures. However, a brightness asymmetry is observed along the minor axis of the disk, suggesting the disk is optically and geometrically thick. We estimate the disk mass to be between 0.007 and 0.02 M$_{\odot}$. Furthermore, molecular emission has been detected from various species, including C$^{18}$O (2$-$1), $^{12}$CO (2$-$1), $^{13}$CO (2$-$1), and H$_2$CO (3$_{0,3}-2_{0,2}$, 3$_{2,1}-2_{2,0}$, and 3$_{2,2}-2_{2,1}$). By conducting a position-velocity analysis of the C$^{18}$O (2$-$1) emission, we find that the disk of IRS5N exhibits characteristics consistent with Keplerian rotation around a central protostar with a mass of approximately 0.3 M$_{\odot}$. Additionally, we observe dust continuum emission from the nearby binary source, IRS5a/b. The emission in $^{12}$CO toward IRS5a/b seems to emanate from IRS5b and flow into IRS5a, suggesting material transport between their mutual orbits. The lack of a detected outflow and large-scale negatives in \tlvco~observed toward IRS5N suggests that much of the flux from IRS5N is being resolved out. Due to this substantial surrounding envelope, the central IRS5N protostar is expected to be significantly more massive in the future. △ Less

Submitted 1 September, 2023; originally announced September 2023.

Comments: 25 pages, 24 figures. Accepted for publication in ApJ

Model-free generalized fiducial inference

Authors: Jonathan P Williams

Abstract: Motivated by the need for the development of safe and reliable methods for uncertainty quantification in machine learning, I propose and develop ideas for a model-free statistical framework for imprecise probabilistic prediction inference. This framework facilitates uncertainty quantification in the form of prediction sets that offer finite sample control of type 1 errors, a property shared with c… ▽ More Motivated by the need for the development of safe and reliable methods for uncertainty quantification in machine learning, I propose and develop ideas for a model-free statistical framework for imprecise probabilistic prediction inference. This framework facilitates uncertainty quantification in the form of prediction sets that offer finite sample control of type 1 errors, a property shared with conformal prediction sets, but this new approach also offers more versatile tools for imprecise probabilistic reasoning. Furthermore, I propose and consider the theoretical and empirical properties of a precise probabilistic approximation to the model-free imprecise framework. Approximating a belief/plausibility measure pair by an [optimal in some sense] probability measure in the credal set is a critical resolution needed for the broader adoption of imprecise probabilistic approaches to inference in statistical and machine learning communities. It is largely undetermined in the statistical and machine learning literatures, more generally, how to properly quantify uncertainty in that there is no generally accepted standard of accountability of stated uncertainties. The research I present in this manuscript is aimed at motivating a framework for statistical inference with reliability and accountability as the guiding principles. △ Less

Submitted 23 July, 2023; originally announced July 2023.

Leveraging arbitrary mobile sensor trajectories with shallow recurrent decoder networks for full-state reconstruction

Authors: Megan R. Ebers, Jan P. Williams, Katherine M. Steele, J. Nathan Kutz

Abstract: Sensing is one of the most fundamental tasks for the monitoring, forecasting and control of complex, spatio-temporal systems. In many applications, a limited number of sensors are mobile and move with the dynamics, with examples including wearable technology, ocean monitoring buoys, and weather balloons. In these dynamic systems (without regions of statistical-independence), the measurement time h… ▽ More Sensing is one of the most fundamental tasks for the monitoring, forecasting and control of complex, spatio-temporal systems. In many applications, a limited number of sensors are mobile and move with the dynamics, with examples including wearable technology, ocean monitoring buoys, and weather balloons. In these dynamic systems (without regions of statistical-independence), the measurement time history encodes a significant amount of information that can be extracted for critical tasks. Most model-free sensing paradigms aim to map current sparse sensor measurements to the high-dimensional state space, ignoring the time-history all together. Using modern deep learning architectures, we show that a sequence-to-vector model, such as an LSTM (long, short-term memory) network, with a decoder network, dynamic trajectory information can be mapped to full state-space estimates. Indeed, we demonstrate that by leveraging mobile sensor trajectories with shallow recurrent decoder networks, we can train the network (i) to accurately reconstruct the full state space using arbitrary dynamical trajectories of the sensors, (ii) the architecture reduces the variance of the mean-square error of the reconstruction error in comparison with immobile sensors, and (iii) the architecture also allows for rapid generalization (parameterization of dynamics) for data outside the training set. Moreover, the path of the sensor can be chosen arbitrarily, provided training data for the spatial trajectory of the sensor is available. The exceptional performance of the network architecture is demonstrated on three applications: turbulent flows, global sea-surface temperature data, and human movement biomechanics. △ Less

Submitted 20 July, 2023; originally announced July 2023.

Comments: 11 pages, 5 figures, 2 tables

Radio multiwavelength analysis of the compact disk CX Tau: Presence of strong free-free variability or anomalous microwave emission

Authors: Pietro Curone, Leonardo Testi, Enrique Macias, Marco Tazzari, Stefano Facchini, Jonathan P. Williams, Cathie J. Clarke, Antonella Natta, Giovanni Rosotti, Claudia Toci, Giuseppe Lodato

Abstract: Protoplanetary disks emit radiation across a broad range of wavelengths, requiring a multiwavelength approach to fully understand their physical mechanisms and how they form planets. Observations at sub-millimeter to centimeter wavelengths can provide insights into the thermal emission from dust, free-free emission from ionized gas, and possible gyro-synchrotron emission from the stellar magnetosp… ▽ More Protoplanetary disks emit radiation across a broad range of wavelengths, requiring a multiwavelength approach to fully understand their physical mechanisms and how they form planets. Observations at sub-millimeter to centimeter wavelengths can provide insights into the thermal emission from dust, free-free emission from ionized gas, and possible gyro-synchrotron emission from the stellar magnetosphere. This work is focused on CX Tau, a ${\sim}0.4\,M_\odot$ star with an extended gas emission and a compact and apparently structureless dust disk, with an average millimeter flux when compared to Class II sources in Taurus. We present observations from the Karl G. Jansky Very Large Array (VLA) across four bands (between 9.0 mm and 6.0 cm) and combine them with archival data from the Atacama Large Millimeter/submillimeter Array (ALMA), the Submillimeter Array (SMA) and the Plateau de Bure Interferometer (PdBI). This multiwavelength approach allows us to separate the dust continuum from other emissions. After isolating the dust thermal emission, we derived an upper limit of the dust disk extent at 1.3 cm which is consistent with theoretical predictions of a radial drift-dominated disk. Centimeter data show a peculiar behavior: deep observations at 6.0 cm did not detect the source, while at 1.3 cm the flux density is anomalously higher than adjacent bands. Intraband spectral indices suggest a dominant contribution from free-free emission, whereas gyro-synchrotron emission is excluded. To explain these observations, we propose a strong variability among the free-free emission with timescales shorter than a month. Another possible interpretation is the presence of anomalous microwave emission from spinning dust grains. △ Less

Submitted 22 August, 2023; v1 submitted 20 July, 2023; originally announced July 2023.

Comments: Accepted for publication in A&A; 13 pages, 13 figures, 1 table

Journal ref: A&A 677, A118 (2023)

Early Planet Formation in Embedded Disks (eDisk) V: Possible Annular Substructure in a Circumstellar Disk in the Ced110 IRS4 System

Authors: Jinshi Sai, Hsi-Wei Yen, Nagayoshi Ohashi, John J. Tobin, Jes K. Jørgensen, Shigehisa Takakuwa, Kazuya Saigo, Yusuke Aso, Zhe-Yu Daniel Lin, Patrick M. Koch, Yuri Aikawa, Christian Flores, Itziar de Gregorio-Monsalvo, Ilseung Han, Miyu Kido, Woojin Kwon, Shih-Ping Lai, Chang Won Lee, Jeong-Eun Lee, Zhi-Yun Li, Leslie W. Looney, Shoji Mori, Nguyen Thi Phuong, Alejandro Santamaría-Miranda, Rajeeb Sharma , et al. (3 additional authors not shown)

Abstract: We have observed the Class 0/I protostellar system Ced110 IRS4 at an angular resolution of $0.05''$ ($\sim$10 au) as a part of the ALMA large program; Early Planet Formation in the Embedded Disks (eDisk). The 1.3 mm dust continuum emission reveals that Ced110 IRS4 is a binary system with a projected separation of $\sim$250 au. The continuum emissions associated with the main source and its compani… ▽ More We have observed the Class 0/I protostellar system Ced110 IRS4 at an angular resolution of $0.05''$ ($\sim$10 au) as a part of the ALMA large program; Early Planet Formation in the Embedded Disks (eDisk). The 1.3 mm dust continuum emission reveals that Ced110 IRS4 is a binary system with a projected separation of $\sim$250 au. The continuum emissions associated with the main source and its companion, named Ced110 IRS4A and IRS4B respectively, exhibit disk-like shapes and likely arise from dust disks around the protostars. The continuum emission of Ced110 IRS4A has a radius of $\sim$91.7 au ($\sim0.485''$), and shows bumps along its major axis with an asymmetry. The bumps can be interpreted as an shallow, ring-like structure at a radius of $\sim$40 au ($\sim0.2''$) in the continuum emission, as demonstrated from two-dimensional intensity distribution models. A rotation curve analysis on the C$^{18}$O and $^{13}$CO $J=2$-1 lines reveals the presence of a Keplerian disk within a radius of 120 au around Ced110 IRS4A, which supports the interpretation that the dust continuum emission arises from a disk. The ring-like structure in the dust continuum emission might indicate a possible, annular substructure in the surface density of the embedded disk, although the possibility that it is an apparent structure due to the optically thick continuum emission cannot be ruled out. △ Less

Submitted 31 August, 2023; v1 submitted 17 July, 2023; originally announced July 2023.

Comments: 32 pages, 23 figures. Accepted for publication in ApJ as one of the first-look papers of the eDisk ALMA Large Program

Early Planet Formation in Embedded Disks (eDisk). VII. Keplerian Disk, Disk Substructure, and Accretion Streamers in the Class 0 Protostar IRAS 16544-1604 in CB 68

Authors: Miyu Kido, Shigehisa Takakuwa, Kazuya Saigo, Nagayoshi Ohashi, John J. Tobin, Jes K, Jørgensen, Yuri Aikawa, Yusuke Aso, Frankie J. Encalada, Christian Flores, Sacha Gavino, Itziar de Gregorio-Monsalvo, Ilseung Han, Shingo Hirano, Patrick M. Koch, Woojin Kwon, Shih-Ping Lai, Chang Won Lee, Jeong-Eun Lee, Zhi-Yun Li, Zhe-Yu Daniel Lin, Leslie W. Looney, Shoji Mori, Suchitra Narayanan , et al. (12 additional authors not shown)

Abstract: We present observations of the Class 0 protostar IRAS 16544-1604 in CB 68 from the ''Early Planet Formation in Embedded Disks (eDisk)'' ALMA Large program. The ALMA observations target continuum and lines at 1.3-mm with an angular resolution of $\sim$5 au. The continuum image reveals a dusty protostellar disk with a radius of $\sim$30 au seen close to edge-on, and asymmetric structures both along… ▽ More We present observations of the Class 0 protostar IRAS 16544-1604 in CB 68 from the ''Early Planet Formation in Embedded Disks (eDisk)'' ALMA Large program. The ALMA observations target continuum and lines at 1.3-mm with an angular resolution of $\sim$5 au. The continuum image reveals a dusty protostellar disk with a radius of $\sim$30 au seen close to edge-on, and asymmetric structures both along the major and minor axes. While the asymmetry along the minor axis can be interpreted as the effect of the dust flaring, the asymmetry along the major axis comes from a real non-axisymmetric structure. The C$^{18}$O image cubes clearly show the gas in the disk that follows a Keplerian rotation pattern around a $\sim$0.14 $M_{\odot}$ central protostar. Furthermore, there are $\sim$1500 au-scale streamer-like features of gas connecting from North-East, North-North-West, and North-West to the disk, as well as the bending outflow as seen in the $^{12}$CO (2-1) emission. At the apparent landing point of NE streamer, there are SO (6$_5$-5$_4$) and SiO (5-4) emission detected. The spatial and velocity structure of NE streamer can be interpreted as a free-falling gas with a conserved specific angular momentum, and the detection of the SO and SiO emission at the tip of the streamer implies presence of accretion shocks. Our eDisk observations have unveiled that the Class 0 protostar in CB 68 has a Keplerian rotating disk with flaring and non-axisymmetric structure associated with accretion streamers and outflows. △ Less

Submitted 27 June, 2023; originally announced June 2023.

Comments: 30 pages, 24 figures, accepted for publication in The Astrophysical Journal as one of the first-look papers of the eDisk ALMA Large Program

Early Planet Formation in Embedded Disks (eDisk). II. Limited Dust Settling and Prominent Snow Surfaces in the Edge-on Class I Disk IRAS 04302+2247

Authors: Zhe-Yu Daniel Lin, Zhi-Yun Li, John J. Tobin, Nagayoshi Ohashi, Jes Kristian Jørgensen, Leslie W. Looney, Yusuke Aso, Shigehisa Takakuwa, Yuri Aikawa, Merel L. R. van 't Hoff, Itziar de Gregorio-Monsalvo, Frankie J. Encalada, Christian Flores, Sacha Gavino, Ilseung Han, Miyu Kido, Patrick M. Koch, Woojin Kwon, Shih-Ping Lai, Chang Won Lee, Jeong-Eun Lee, Nguyen Thi Phuong, Jinshi Sai, Rajeeb Sharma, Patrick Sheehan , et al. (4 additional authors not shown)

Abstract: While dust disks around optically visible, Class II protostars are found to be vertically thin, when and how dust settles to the midplane are unclear. As part of the Atacama Large Millimeter/submillimeter Array (ALMA) large program, Early Planet Formation in Embedded Disks, we analyze the edge-on, embedded, Class I protostar IRAS 04302+2247, also nicknamed the ``Butterfly Star." With a resolution… ▽ More While dust disks around optically visible, Class II protostars are found to be vertically thin, when and how dust settles to the midplane are unclear. As part of the Atacama Large Millimeter/submillimeter Array (ALMA) large program, Early Planet Formation in Embedded Disks, we analyze the edge-on, embedded, Class I protostar IRAS 04302+2247, also nicknamed the ``Butterfly Star." With a resolution of 0.05" (8~au), the 1.3 mm continuum shows an asymmetry along the minor axis which is evidence of an optically thick and geometrically thick disk viewed nearly edge-on. There is no evidence of rings and gaps, which could be due to the lack of radial substructure or the highly inclined and optically thick view. With 0.1" (16~au) resolution, we resolve the 2D snow surfaces, i.e., the boundary region between freeze-out and sublimation, for $^{12}$CO $J$=2--1, $^{13}$CO $J$=2--1, C$^{18}$O $J$=2--1, $H_{2}$CO $J$=$3_{0,3}$--$2_{0,2}$, and SO $J$=$6_{5}$--$5_{4}$, and constrain the CO midplane snow line to $\sim 130$ au. We find Keplerian rotation around a protostar of $1.6 \pm 0.4 M_{\odot}$ using C$^{18}$O. Through forward ray-tracing using RADMC-3D, we find that the dust scale height is $\sim 6$ au at a radius of 100~au from the central star and is comparable to the gas pressure scale height. The results suggest that the dust of this Class~I source has yet to vertically settle significantly. △ Less

Submitted 27 June, 2023; originally announced June 2023.

Comments: 33 pages, 21 figures. Accepted for publication in ApJ as one of the first-look papers of the eDisk ALMA Large Program

Early Planet Formation in Embedded Disks (eDisk) III: A first high-resolution view of sub-mm continuum and molecular line emission toward the Class 0 protostar L1527 IRS

Authors: Merel L. R. van 't Hoff, John J. Tobin, Zhi-Yun Li, Nagayoshi Ohashi, Jes K. Jørgensen, Zhe-Yu Daniel Lin, Yuri Aikawa, Yusuke Aso, Itziar de Gregorio-Monsalvo, Sacha Gavino, Ilseung Han, Patrick M. Koch, Woojin Kwon, Chang Won Lee, Jeong-Eun Lee, Leslie W. Looney, Suchitra Narayanan, Adele Plunkett, Jinshi Sai, Alejandro Santamaría-Miranda, Rajeeb Sharma, Patrick D. Sheehan, Shigehisa Takakuwa, Travis J. Thieme, Jonathan P. Williams , et al. (3 additional authors not shown)

Abstract: Studying the physical and chemical conditions of young embedded disks is crucial to constrain the initial conditions for planet formation. Here, we present Atacama Large Millimeter/submillimeter Array (ALMA) observations of dust continuum at $\sim$0.06" (8 au) resolution and molecular line emission at $\sim$0.17" (24 au) resolution toward the Class 0 protostar L1527 IRS from the Large Program eDis… ▽ More Studying the physical and chemical conditions of young embedded disks is crucial to constrain the initial conditions for planet formation. Here, we present Atacama Large Millimeter/submillimeter Array (ALMA) observations of dust continuum at $\sim$0.06" (8 au) resolution and molecular line emission at $\sim$0.17" (24 au) resolution toward the Class 0 protostar L1527 IRS from the Large Program eDisk (Early Planet Formation in Embedded Disks). The continuum emission is smooth without substructures, but asymmetric along both the major and minor axes of the disk as previously observed. The detected lines of $^{12}$CO, $^{13}$CO, C$^{18}$O, H$_2$CO, c-C$_3$H$_2$, SO, SiO, and DCN trace different components of the protostellar system, with a disk wind potentially visible in $^{12}$CO. The $^{13}$CO brightness temperature and the H$_2$CO line ratio confirm that the disk is too warm for CO freeze out, with the snowline located at $\sim$350 au in the envelope. Both molecules show potential evidence of a temperature increase around the disk-envelope interface. SO seems to originate predominantly in UV-irradiated regions such as the disk surface and the outflow cavity walls rather than at the disk-envelope interface as previously suggested. Finally, the continuum asymmetry along the minor axis is consistent with the inclination derived from the large-scale (100" or 14,000 au) outflow, but opposite to that based on the molecular jet and envelope emission, suggesting a misalignment in the system. Overall, these results highlight the importance of observing multiple molecular species in multiple transitions to characterize the physical and chemical environment of young disks. △ Less

Submitted 27 June, 2023; originally announced June 2023.

Comments: 27 pages, 16 figures, 2 tables, 10 pages appendix with 12 figures. Accepted for publication in ApJ as one of the first-look papers of the eDisk ALMA Large Program

Resolving the Binary Components of the Outbursting Protostar HBC 494 with ALMA

Authors: Pedro Henrique Nogueira, Alice Zurlo, Sebastián Pérez, Camilo González-Ruilova, Lucas A. Cieza, Antonio Hales, Trisha Bhowmik, Dary A. Ruíz-Rodríguez, David A. Principe, Gregory J. Herczeg, Jonathan P. Williams, Jorge Cuadra, Matías Montesinos, Nicolás Cuello, Prachi Chavan, Simon Casassus, Zhaohuan Zhu, Felipe G. Goicovic

Abstract: Episodic accretion is a low-mass pre-main sequence phenomenon characterized by sudden outbursts of enhanced accretion. These objects are classified into two: protostars with elevated levels of accretion that lasts for decades or more, called FUors, and protostars with shorter and repetitive bursts, called EXors. HBC 494 is a FUor object embedded in the Orion Molecular Cloud. Earlier Atacama Large… ▽ More Episodic accretion is a low-mass pre-main sequence phenomenon characterized by sudden outbursts of enhanced accretion. These objects are classified into two: protostars with elevated levels of accretion that lasts for decades or more, called FUors, and protostars with shorter and repetitive bursts, called EXors. HBC 494 is a FUor object embedded in the Orion Molecular Cloud. Earlier Atacama Large (sub-)Millimeter Array (ALMA) continuum observations showed an asymmetry in the disk at 0.''2 resolution. Here, we present follow-up observations at ~0.''03, resolving the system into two components: HBC 494 N (primary) and HBC 494 S (secondary). No circumbinary disk was detected. Both disks are resolved with a projected separation of ~0.''18 (75 au). Their projected dimensions are 84+/-1.8 x 66.9+/-1.5 mas for HBC 494 N and 64.6+/-2.5 x 46.0+/-1.9 mas for HBC 494 S. The disks are almost aligned and with similar inclinations. The observations show that the primary is ~5 times brighter/more massive and ~2 times bigger than the secondary. We notice that the northern component has a similar mass to the FUors, while the southern has to EXors. The HBC 494 disks show individual sizes that are smaller than single eruptive YSOs. In this work, we also report 12CO, 13CO, and C18O molecular line observations. At large scale, the 12CO emission shows bipolar outflows, while the 13CO and C18O maps show a rotating and infalling envelope around the system. At a smaller scale, the 12CO and 13CO moment zero maps show cavities within the continuum disks' area, which may indicate continuum over-subtraction or slow-moving jets and chemical destruction along the line-of-sight. △ Less

Submitted 26 June, 2023; v1 submitted 24 May, 2023; originally announced May 2023.

Comments: Published in MNRAS. Main text (11 pages, 8 figures). Appendix (11 pages, 9 figures). Total (22 pages, 17 figures)

Introduction to Generalized Fiducial Inference

Authors: Alexander C. Murph, Jan Hannig, Jonathan P. Williams

Abstract: Fiducial inference was introduced in the first half of the 20th century by Fisher (1935) as a means to get a posterior-like distribution for a parameter without having to arbitrarily define a prior. While the method originally fell out of favor due to non-exactness issues in multivariate cases, the method has garnered renewed interest in the last decade. This is partly due to the development of ge… ▽ More Fiducial inference was introduced in the first half of the 20th century by Fisher (1935) as a means to get a posterior-like distribution for a parameter without having to arbitrarily define a prior. While the method originally fell out of favor due to non-exactness issues in multivariate cases, the method has garnered renewed interest in the last decade. This is partly due to the development of generalized fiducial inference, which is a fiducial perspective on generalized confidence intervals: a method used to find approximate confidence distributions. In this chapter, we illuminate the usefulness of the fiducial philosophy, introduce the definition of a generalized fiducial distribution, and apply it to interesting, non-trivial inferential examples. △ Less

Submitted 28 February, 2023; originally announced February 2023.

Comments: 20 pages, 4 tables, 11 figures

Valid Inference for Machine Learning Model Parameters

Authors: Neil Dey, Jonathan P. Williams

Abstract: The parameters of a machine learning model are typically learned by minimizing a loss function on a set of training data. However, this can come with the risk of overtraining; in order for the model to generalize well, it is of great importance that we are able to find the optimal parameter for the model on the entire population -- not only on the given training sample. In this paper, we construct… ▽ More The parameters of a machine learning model are typically learned by minimizing a loss function on a set of training data. However, this can come with the risk of overtraining; in order for the model to generalize well, it is of great importance that we are able to find the optimal parameter for the model on the entire population -- not only on the given training sample. In this paper, we construct valid confidence sets for this optimal parameter of a machine learning model, which can be generated using only the training data without any knowledge of the population. We then show that studying the distribution of this confidence set allows us to assign a notion of confidence to arbitrary regions of the parameter space, and we demonstrate that this distribution can be well-approximated using bootstrapping techniques. △ Less

Submitted 9 May, 2024; v1 submitted 21 February, 2023; originally announced February 2023.

Comments: 35 pages, 6 figures

Sensing with shallow recurrent decoder networks

Authors: Jan P. Williams, Olivia Zahn, J. Nathan Kutz

Abstract: Sensing is a universal task in science and engineering. Downstream tasks from sensing include inferring full state estimates of a system (system identification), control decisions, and forecasting. These tasks are exceptionally challenging to achieve with limited sensors, noisy measurements, and corrupt or missing data. We propose a SHallow REcurrent Decoder (SHRED) neural network structure for se… ▽ More Sensing is a universal task in science and engineering. Downstream tasks from sensing include inferring full state estimates of a system (system identification), control decisions, and forecasting. These tasks are exceptionally challenging to achieve with limited sensors, noisy measurements, and corrupt or missing data. We propose a SHallow REcurrent Decoder (SHRED) neural network structure for sensing which incorporates (i) a recurrent neural network (LSTM) to learn a latent representation of the temporal dynamics of the sensors, and (ii) a shallow decoder that learns a mapping between this latent representation and the high-dimensional state space. By explicitly accounting for the time-history, or trajectory, of the sensor measurements, SHRED enables accurate reconstructions with far fewer sensors, outperforms existing techniques when more measurements are available, and is agnostic towards sensor placement. In addition, a compressed representation of the high-dimensional state is directly obtained from sensor measurements, which provides an on-the-fly compression for modeling physical and engineering systems. Forecasting is also achieved from the sensor time-series data alone, producing an efficient paradigm for predicting temporal evolution with an exceptionally limited number of sensors. In the example cases explored, including turbulent flows, complex spatio-temporal dynamics can be characterized with exceedingly limited sensors that can be randomly placed with minimal loss of performance. △ Less

Submitted 5 June, 2024; v1 submitted 27 January, 2023; originally announced January 2023.

Comments: 12 pages, 7 figures

Word Embeddings as Statistical Estimators

Authors: Neil Dey, Matthew Singer, Jonathan P. Williams, Srijan Sengupta

Abstract: Word embeddings are a fundamental tool in natural language processing. Currently, word embedding methods are evaluated on the basis of empirical performance on benchmark data sets, and there is a lack of rigorous understanding of their theoretical properties. This paper studies word embeddings from a statistical theoretical perspective, which is essential for formal inference and uncertainty quant… ▽ More Word embeddings are a fundamental tool in natural language processing. Currently, word embedding methods are evaluated on the basis of empirical performance on benchmark data sets, and there is a lack of rigorous understanding of their theoretical properties. This paper studies word embeddings from a statistical theoretical perspective, which is essential for formal inference and uncertainty quantification. We propose a copula-based statistical model for text data and show that under this model, the now-classical Word2Vec method can be interpreted as a statistical estimation method for estimating the theoretical pointwise mutual information (PMI). Next, by building on the work of Levy and Goldberg (2014), we develop a missing value-based estimator as a statistically tractable and interpretable alternative to the Word2Vec approach. The estimation error of this estimator is comparable to Word2Vec and improves upon the truncation-based method proposed by Levy and Goldberg (2014). The proposed estimator also performs comparably to Word2Vec in a benchmark sentiment analysis task on the IMDb Movie Reviews data set. △ Less

Submitted 17 January, 2023; originally announced January 2023.

Reflections on nebulae around young stars: A systematic search for late-stage infall of material onto Class II disks

Authors: Aashish Gupta, Anna Miotello, Carlo F. Manara, Jonathan P. Williams, Stefano Facchini, Giacomo Beccari, Til Birnstiel, Christian Ginski, Alvaro Hacar, Michael Küffmeier, Leonardo Testi, Lukasz Tychoniec, Hsi-Wei Yen

Abstract: Context. While it is generally assumed that Class II sources evolve largely in isolation from their environment, many still lie close to molecular clouds and may continue to interact with them. This may result in late accretion of material onto the disk that can significantly influence disk processes and planet formation. Aims. In order to systematically study late infall of gas onto disks, we i… ▽ More Context. While it is generally assumed that Class II sources evolve largely in isolation from their environment, many still lie close to molecular clouds and may continue to interact with them. This may result in late accretion of material onto the disk that can significantly influence disk processes and planet formation. Aims. In order to systematically study late infall of gas onto disks, we identify candidate Class II sources in close vicinity to a reflection nebula (RN) that may be undergoing this process. Methods. First we targeted Class II sources with known kilo-au scale gas structures - possibly due to late infall of material - and we searched for RNe in their vicinity in optical and near-infrared images. Second, we compiled a catalogue of Class II sources associated with RNe and looked for the large-scale CO structures in archival ALMA data. Using the catalogues of protostars and RNe, we also estimated the probability of Class II sources interacting with surrounding material. Results. All of the sources with large-scale gas structures also exhibit some reflection nebulosity in their vicinity. Similarly, at least five Class II objects associated with a prominent RNe, and for which adequate ALMA observations are available, were found to have spirals or stream-like structures which may be due to late infall. We report the first detection of these structures around S CrA. Conclusions. Our results suggest that a non-negligible fraction of Class II disks in nearby star-forming regions may be associated with RNe and could therefore be undergoing late accretion of gas. Surveys of RNe and kilo-au scale gas structures around Class II sources will allow us to better understand the frequency and impact of late-infall phenomena. △ Less

Submitted 8 January, 2023; originally announced January 2023.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 670, L8 (2023)

Transfer Learning with Uncertainty Quantification: Random Effect Calibration of Source to Target (RECaST)

Authors: Jimmy Hickey, Jonathan P. Williams, Emily C. Hector

Abstract: Transfer learning uses a data model, trained to make predictions or inferences on data from one population, to make reliable predictions or inferences on data from another population. Most existing transfer learning approaches are based on fine-tuning pre-trained neural network models, and fail to provide crucial uncertainty quantification. We develop a statistical framework for model predictions… ▽ More Transfer learning uses a data model, trained to make predictions or inferences on data from one population, to make reliable predictions or inferences on data from another population. Most existing transfer learning approaches are based on fine-tuning pre-trained neural network models, and fail to provide crucial uncertainty quantification. We develop a statistical framework for model predictions based on transfer learning, called RECaST. The primary mechanism is a Cauchy random effect that recalibrates a source model to a target population; we mathematically and empirically demonstrate the validity of our RECaST approach for transfer learning between linear models, in the sense that prediction sets will achieve their nominal stated coverage, and we numerically illustrate the method's robustness to asymptotic approximations for nonlinear models. Whereas many existing techniques are built on particular source models, RECaST is agnostic to the choice of source model. For example, our RECaST transfer learning approach can be applied to a continuous or discrete data model with linear or logistic regression, deep neural network architectures, etc. Furthermore, RECaST provides uncertainty quantification for predictions, which is mostly absent in the literature. We examine our method's performance in a simulation study and in an application to real hospital data. △ Less

Submitted 29 November, 2022; originally announced November 2022.

Comments: 26 pages, 2 figures, 4 tables

Beyond time-homogeneity for continuous-time multistate Markov models

Authors: Emmett B. Kendall, Jonathan P. Williams, Gudmund H. Hermansen, Frederic Bois, Vo Hong Thanh

Abstract: Multistate Markov models are a canonical parametric approach for data modeling of observed or latent stochastic processes supported on a finite state space. Continuous-time Markov processes describe data that are observed irregularly over time, as is often the case in longitudinal medical data, for example. Assuming that a continuous-time Markov process is time-homogeneous, a closed-form likelihoo… ▽ More Multistate Markov models are a canonical parametric approach for data modeling of observed or latent stochastic processes supported on a finite state space. Continuous-time Markov processes describe data that are observed irregularly over time, as is often the case in longitudinal medical data, for example. Assuming that a continuous-time Markov process is time-homogeneous, a closed-form likelihood function can be derived from the Kolmogorov forward equations -- a system of differential equations with a well-known matrix-exponential solution. Unfortunately, however, the forward equations do not admit an analytical solution for continuous-time, time-inhomogeneous Markov processes, and so researchers and practitioners often make the simplifying assumption that the process is piecewise time-homogeneous. In this paper, we provide intuitions and illustrations of the potential biases for parameter estimation that may ensue in the more realistic scenario that the piecewise-homogeneous assumption is violated, and we advocate for a solution for likelihood computation in a truly time-inhomogeneous fashion. Particular focus is afforded to the context of multistate Markov models that allow for state label misclassifications, which applies more broadly to hidden Markov models (HMMs), and Bayesian computations bypass the necessity for computationally demanding numerical gradient approximations for obtaining maximum likelihood estimates (MLEs). Supplemental materials are available online. △ Less

Submitted 16 July, 2024; v1 submitted 6 November, 2022; originally announced November 2022.

Generalized Fiducial Inference on Differentiable Manifolds

Authors: Alexander C Murph, Jan Hannig, Jonathan P Williams

Abstract: We introduce a novel approach to inference on parameters that take values in a Riemannian manifold embedded in a Euclidean space. Parameter spaces of this form are ubiquitous across many fields, including chemistry, physics, computer graphics, and geology. This new approach uses generalized fiducial inference to obtain a posterior-like distribution on the manifold, without needing to know a parame… ▽ More We introduce a novel approach to inference on parameters that take values in a Riemannian manifold embedded in a Euclidean space. Parameter spaces of this form are ubiquitous across many fields, including chemistry, physics, computer graphics, and geology. This new approach uses generalized fiducial inference to obtain a posterior-like distribution on the manifold, without needing to know a parameterization that maps the constrained space to an unconstrained Euclidean space. The proposed methodology, called the constrained generalized fiducial distribution (CGFD), is obtained by using mathematical tools from Riemannian geometry. A Bernstein-von Mises-type result for the CGFD, which provides intuition for how the desirable asymptotic qualities of the unconstrained generalized fiducial distribution are inherited by the CGFD, is provided. To demonstrate the practical use of the CGFD, we provide three proof-of-concept examples: inference for data from a multivariate normal density with the mean parameters on a sphere, a linear logspline density estimation problem, and a reimagined approach to the AR(1) model, all of which exhibit desirable coverages via simulation. We discuss two Markov chain Monte Carlo algorithms for the exploration of these constrained parameter spaces and adapt them for the CGFD. △ Less

Submitted 8 December, 2022; v1 submitted 30 September, 2022; originally announced September 2022.

Comments: 31 pages, 7 figures

Discovery of a brown dwarf with quasi-spherical mass-loss

Authors: Dary A. Ruíz-Rodríguez, Lucas A. Cieza, Simon Casassus, Victor Almendros-Abad, Paula Jofré, Koraljka Muzic, Karla Peña Ramirez, Grace Batalla-Falcon, Michael M. Dunham, Camilo González-Ruilova, Antonio Hales, Elizabeth Humphreys, Pedro H. Nogueira, Claudia Paladini, John Tobin, Jonathan P. Williams, Alice Zurlo

Abstract: We report the serendipitous discovery of an elliptical shell of CO associated with the faint stellar object SSTc2d J163134.1-24006 as part of the "Ophiuchus Disk Survey Employing ALMA" (ODISEA), a project aiming to study the entire population of protoplanetary disks in the Ophiuchus Molecular Cloud from 230 GHz continuum emission and $^{12}$CO (J=2-1), $^{13}$CO (J=2-1) and C$^{18}$CO (J=2-1) line… ▽ More We report the serendipitous discovery of an elliptical shell of CO associated with the faint stellar object SSTc2d J163134.1-24006 as part of the "Ophiuchus Disk Survey Employing ALMA" (ODISEA), a project aiming to study the entire population of protoplanetary disks in the Ophiuchus Molecular Cloud from 230 GHz continuum emission and $^{12}$CO (J=2-1), $^{13}$CO (J=2-1) and C$^{18}$CO (J=2-1) lines readable in Band-6. Remarkably, we detect a bright $^{12}$CO elliptical shape emission of $\sim$ 3$^{"}$ $\times$ 4$^{"}$ towards SSTc2d J163134.1-24006 without a 230 GHz continuum detection. Based on the observed near-IR spectrum taken with the Very Large Telescope (KMOS), the brightness of the source, its 3-dimensional motion, and Galactic dynamic arguments, we conclude that the source is not a giant star in the distant background ($>$5 - 10 kpc) and is most likely to be a young brown dwarf in the Ophiuchus cloud, at a distance of just $\sim$139 pc. This is the first report of quasi-spherical mass loss in a young brown dwarf. We suggest that the observed shell could be associated with a thermal pulse produced by the fusion of deuterium, which is not yet well understood, but for a sub-stellar object is expected to occur during a short period of time at an age of a few Myr, in agreement with the ages of the objects in the region. Other more exotic scenarios, such as a merger with planetary companions, cannot be ruled out from the current observations. △ Less

Submitted 1 September, 2022; originally announced September 2022.

Comments: 27 pages, 20 Figures. Accepted ApJ

A Chemical Map of the Outbursting V883 Ori system: Vertical and Radial Structures

Authors: D. A. Ruiz-Rodriguez, J. P. Williams, J. H. Kastner, L. Cieza, M. Leemker, D. A. Principe

Abstract: We present the first results of a pilot program to conduct an Atacama Large Millimeter/submillimeter Array (ALMA) Band 6 (211-275 GHz) spectral line study of young stellar objects (YSO) that are undergoing rapid accretion episodes, i.e. FU Ori objects (FUors). Here, we report on molecular emission line observations of the FUor system, V883 Ori. In order to image the FUor object with full coverage… ▽ More We present the first results of a pilot program to conduct an Atacama Large Millimeter/submillimeter Array (ALMA) Band 6 (211-275 GHz) spectral line study of young stellar objects (YSO) that are undergoing rapid accretion episodes, i.e. FU Ori objects (FUors). Here, we report on molecular emission line observations of the FUor system, V883 Ori. In order to image the FUor object with full coverage from ~0.5 arcsec to the map size of ~30 arcsec, i.e. from disc to outflow scales, we combine the ALMA main array (the 12-m array) with the Atacama Compact Array (7-m array) and the total power (TP) array. We detect HCN, HCO$^{+}$, CH$_{3}$OH, SO, DCN, and H$_{2}$CO emissions with most of these lines displaying complex kinematics. From PV diagrams, the detected molecules HCN, HCO$^{+}$, CH$_{3}$OH, DCN, SO, and H$_{2}$CO probe a Keplerian rotating disc in a direction perpendicular to the large-scale outflow detected previously with the $^{12}$CO and $^{13}$CO lines. Additionally, HCN and HCO$^{+}$ reveal kinematic signatures of infall motion. The north outflow is seen in HCO$^{+}$, H$_{2}$CO, and SO emissions. Interestingly, HCO$^{+}$ emission reveals a pronounced inner depression or "hole" with a size comparable to the radial extension estimated for the CH$_{3}$OH and 230 GHz continuum. The inner depression in the integrated HCO$^{+}$ intensity distribution of V883 Ori is most likely the result of optical depth effects, wherein the optically thick nature of the HCO$^{+}$ and continuum emission towards the innermost parts of V883 Ori can result in a continuum subtraction artifact in the final HCO$^{+}$ flux level. △ Less

Submitted 5 July, 2022; originally announced July 2022.

Distribution of solids in the rings of the HD 163296 disk: a multiwavelength study

Authors: G. Guidi, A. Isella, L. Testi, C. J. Chandler, H. B. Liu, H. M. Schmid, G. Rosotti, C. Meng, J. Jennings, J. P. Williams, J. M. Carpenter, I. de Gregorio-Monsalvo, H. Li, S. F. Liu, S. Ortolani, S. P. Quanz, L. Ricci, M. Tazzari

Abstract: In this paper we analyze new observations from ALMA and VLA, at a high angular resolution corresponding to 5 - 8 au, of the protoplanetary disk around HD 163296 to determine the dust spatial distribution and grain properties. We fit the spectral energy distribution as a function of the radius at five wavelengths from 0.9 to 9\,mm, using a simple power law and a physical model based on an analytic… ▽ More In this paper we analyze new observations from ALMA and VLA, at a high angular resolution corresponding to 5 - 8 au, of the protoplanetary disk around HD 163296 to determine the dust spatial distribution and grain properties. We fit the spectral energy distribution as a function of the radius at five wavelengths from 0.9 to 9\,mm, using a simple power law and a physical model based on an analytic description of radiative transfer that includes isothermal scattering. We considered eight dust populations and compared the models' performance using Bayesian evidence. Our analysis shows that the moderately high optical depth ($τ$>1) at $λ\leq$ 1.3 mm in the dust rings artificially lower the millimeter spectral index, which should therefore not be considered as a reliable direct proxy of the dust properties and especially the grain size. We find that the outer disk is composed of small grains on the order of 200 $μ$m with no significant difference between rings at 66 and 100 au and the adjacent gaps, while in the innermost 30 au, larger grains ($\geq$mm) could be present. We show that the assumptions on the dust composition have a strong impact on the derived surface densities and grain size. In particular, increasing the porosity of the grains to 80\% results in a total dust mass about five times higher with respect to grains with 25\% porosity. Finally, we find that the derived opacities as a function of frequency deviate from a simple power law and that grains with a lower porosity seem to better reproduce the observations of HD163296. While we do not find evidence of differential trapping in the rings of HD163296, our overall results are consistent with the postulated presence of giant planets affecting the dust temperature structure and surface density, and possibly originating a second-generation dust population of small grains. △ Less

Submitted 4 July, 2022; originally announced July 2022.

Comments: 27 pages, 28 figures, Accepted for publication in Astronomy and Astrophysics (A&A)

Journal ref: A&A 664, A137 (2022)

A VLA View of the Flared, Asymmetric Disk Around the Class 0 Protostar L1527 IRS

Authors: Patrick D. Sheehan, John J. Tobin, Zhi-Yun Li, Merel L. R. van 't Hoff, Jes K. Jørgensen, Woojin Kwon, Leslie W. Looney, Nagayoshi Ohashi, Shigehisa Takakuwa, Jonathan P. Williams, Yusuke Aso, Sacha Gavino, Itziar de Gregorio-Monsalvo, Ilseung Han, Chang Won Lee, Adele Plunkett, Rajeeb Sharma, Yuri Aikawa, Shih-Ping Lai, Jeong-Eun Lee, Zhe-Yu Daniel Lin, Kazuya Saigo, Kengo Tomida, Hsi-Wei Yen

Abstract: We present high resolution Karl G. Jansky Very Large Array (VLA) observations of the protostar L1527 IRS at 7 mm, 1.3 cm, and 2 cm wavelengths. We detect the edge-on dust disk at all three wavelengths and find that it is asymmetric, with the southern side of the disk brighter than the northern side. We confirm this asymmetry through analytic modeling and also find that the disk is flared at 7 mm.… ▽ More We present high resolution Karl G. Jansky Very Large Array (VLA) observations of the protostar L1527 IRS at 7 mm, 1.3 cm, and 2 cm wavelengths. We detect the edge-on dust disk at all three wavelengths and find that it is asymmetric, with the southern side of the disk brighter than the northern side. We confirm this asymmetry through analytic modeling and also find that the disk is flared at 7 mm. We test the data against models including gap features in the intensity profile, and though we cannot rule such models out, they do not provide a statistically significant improvement in the quality of fit to the data. From these fits, we can however place constraints on allowed properties of any gaps that could be present in the true, underlying intensity profile. The physical nature of the asymmetry is difficult to associate with physical features due to the edge-on nature of the disk, but could be related to spiral arms or asymmetries seen in other imaging of more face-on disks. △ Less

Submitted 27 June, 2022; originally announced June 2022.

Comments: 20 pages, 5 figures, 3 tables; accepted for publication in the Astrophysical Journal

A Statistical Primer on Classical Period-Finding Techniques in Astronomy

Authors: Naomi Giertych, Ahmed Shaban, Pragya Haravu, Jonathan P Williams

Abstract: The aim of our paper is to investigate the properties of the classical phase-dispersion minimization (PDM), analysis of variance (AOV), string-length (SL), and Lomb-Scargle (LS) power statistics from a statistician's perspective. We confirm that when the data are perturbations of a constant function, i.e. under the null hypothesis of no period in the data, a scaled version of the PDM statistic fol… ▽ More The aim of our paper is to investigate the properties of the classical phase-dispersion minimization (PDM), analysis of variance (AOV), string-length (SL), and Lomb-Scargle (LS) power statistics from a statistician's perspective. We confirm that when the data are perturbations of a constant function, i.e. under the null hypothesis of no period in the data, a scaled version of the PDM statistic follows a beta distribution, the AOV statistic follows an F distribution, and the LS power follows a chi-squared distribution with two degrees of freedom. However, the SL statistic does not have a closed-form distribution. We further verify these theoretical distributions through simulations and demonstrate that the extreme values of these statistics (over a range of trial periods), often used for period estimation and determination of the false alarm probability (FAP), follow different distributions than those derived for a single period. We emphasize that multiple-testing considerations are needed to correctly derive FAP bounds. Though, in fact, multiple-testing controls are built into the FAP bound for these extreme-value statistics, e.g. the FAP bound derived specifically for the maximum LS power statistic over a range of trial periods. Additionally, we find that all of these methods are robust to heteroscedastic noise aimed to mimic the degradation or miscalibration of an instrument over time. Finally, we examine the ability of these statistics to detect a non-constant periodic function via simulating data that mimics a well-detached binary system, and we find that the AOV statistic has the most power to detect the correct period, which agrees with what has been observed in practice. △ Less

Submitted 20 March, 2024; v1 submitted 20 May, 2022; originally announced May 2022.

Protostellar and Protoplanetary Disk Masses in the Serpens Region

Authors: Alexa R. Anderson, Jonathan P. Williams, Nienke van der Marel, Charles J. Law, Luca Ricci, John J. Tobin, Simin Tong

Abstract: We present the results from an Atacama Large Millimeter/Submillimeter Array (ALMA) 1.3 mm continuum and $^{12}$CO ($J=2-1$) line survey spread over 10 square degrees in the Serpens star-forming region of 320 young stellar objects, 302 of which are likely members of Serpens (16 Class I, 35 Flat spectrum, 235 Class II, and 16 Class III). From the continuum data, we derive disk dust masses and show t… ▽ More We present the results from an Atacama Large Millimeter/Submillimeter Array (ALMA) 1.3 mm continuum and $^{12}$CO ($J=2-1$) line survey spread over 10 square degrees in the Serpens star-forming region of 320 young stellar objects, 302 of which are likely members of Serpens (16 Class I, 35 Flat spectrum, 235 Class II, and 16 Class III). From the continuum data, we derive disk dust masses and show that they systematically decline from Class I to Flat spectrum to Class II sources. Grouped by stellar evolutionary state, the disk mass distributions are similar to other young ($<3$ Myr) regions, indicating that the large scale environment of a star-forming region does not strongly affect its overall disk dust mass properties. These comparisons between populations reinforce previous conclusions that disks in the Ophiuchus star-forming region have anomalously low masses at all evolutionary stages. Additionally, we find a single deeply embedded protostar that has not been documented elsewhere in the literature and, from the CO line data, 15 protostellar outflows which we catalog here. △ Less

Submitted 19 April, 2022; originally announced April 2022.

Comments: 14 pages, 9 figures, 3 tables, accepted to ApJ

High-resolution ALMA observations of transition disk candidates in Lupus

Authors: Nienke van der Marel, Jonathan P. Williams, Giovanni Picogna, Sierk van Terwisga, Stefano Facchini, Carlo F. Manara, Apostolos Zormpas, Megan Ansdell, .

Abstract: Transition disks with small inner dust cavities are interesting targets for the study of disk clearing mechanisms. Such disks have been identified through a deficit in the infrared part of their SED, but spatially resolved millimeter imaging is required to confirm the presence of an inner dust cavity. We use high-resolution ALMA observations of 30 mas resolution in Band 6 continuum and $^{12}$CO 2… ▽ More Transition disks with small inner dust cavities are interesting targets for the study of disk clearing mechanisms. Such disks have been identified through a deficit in the infrared part of their SED, but spatially resolved millimeter imaging is required to confirm the presence of an inner dust cavity. We use high-resolution ALMA observations of 30 mas resolution in Band 6 continuum and $^{12}$CO 2--1 emission of 10 transition disk candidates in the Lupus star forming region, in order to confirm the presence of inner dust cavities and infer the responsible mechanism. The continuum data are analyzed using visibility modeling and the SEDs are compared with radiative transfer models. Out of the six transition disk candidates selected from their SED, only one disk revealed an inner dust cavity of 4 au in radius. Three of the other disks are highly inclined, which limits the detectability of an inner dust cavity but it is also demonstrated to be the possible cause for the infrared deficit in their SED. The two remaining SED-selected disks are very compact, with dust radii of only $\sim$3 au. From the four candidates selected from low-resolution images, three new transition disks with large inner cavities $>$20 au are identified, bringing the total number of transition disks with large cavities in Lupus to 13. SED-selected transition disks with small cavities are biased towards highly inclined and compact disks, which casts doubt on the use of their occurrence rates in estimating dispersal timescales of photoevaporation. Using newly derived disk dust masses and radii, we re-evaluate the size-luminosity and $M_{\rm dust}-M_{\rm star}$ relations. These relations can be understood if the bright disks are dominated by disks with substructure whereas faint disks are dominated by drift-dominated disks. (Abridged) △ Less

Submitted 18 April, 2022; originally announced April 2022.

Comments: Submitted to A&A. Currently in revision. 28 pages, 13 figures

Disk Evolution Study Through Imaging of Nearby Young Stars (DESTINYS): A Panchromatic View of DO Tau's Complex Kilo-au Environment

Authors: Jane Huang, Christian Ginski, Myriam Benisty, Bin Ren, Alexander J. Bohn, Élodie Choquet, Karin I. Öberg, Álvaro Ribas, Jaehan Bae, Edwin A. Bergin, Til Birnstiel, Yann Boehler, Stefano Facchini, Daniel Harsono, Michiel Hogerheijde, Feng Long, Carlo F. Manara, François Ménard, Paola Pinilla, Christophe Pinte, Christian Rab, Jonathan P. Williams, Alice Zurlo

Abstract: While protoplanetary disks are often treated as isolated systems in planet formation models, observations increasingly suggest that vigorous interactions between Class II disks and their environments are not rare. DO Tau is a T Tauri star that has previously been hypothesized to have undergone a close encounter with the HV Tau system. As part of the DESTINYS ESO Large Programme, we present new VLT… ▽ More While protoplanetary disks are often treated as isolated systems in planet formation models, observations increasingly suggest that vigorous interactions between Class II disks and their environments are not rare. DO Tau is a T Tauri star that has previously been hypothesized to have undergone a close encounter with the HV Tau system. As part of the DESTINYS ESO Large Programme, we present new VLT/SPHERE polarimetric observations of DO Tau and combine them with archival HST scattered light images and ALMA observations of CO isotopologues and CS to map a network of complex structures. The SPHERE and ALMA observations show that the circumstellar disk is connected to arms extending out to several hundred au. HST and ALMA also reveal stream-like structures northeast of DO Tau, some of which are at least several thousand au long. These streams appear not to be gravitationally bound to DO Tau, and comparisons with previous Herschel far-IR observations suggest that the streams are part of a bridge-like structure connecting DO Tau and HV Tau. We also detect a fainter redshifted counterpart to a previously known blueshifted CO outflow. While some of DO Tau's complex structures could be attributed to a recent disk-disk encounter, they might be explained alternatively by interactions with remnant material from the star formation process. These panchromatic observations of DO Tau highlight the need to contextualize the evolution of Class II disks by examining processes occurring over a wide range of size scales. △ Less

Submitted 8 May, 2022; v1 submitted 4 April, 2022; originally announced April 2022.

Comments: 42 pages, 20 figures, accepted by ApJ, reduced data available at https://zenodo.org/record/6408903 (typos from v1 fixed)

Demographics of young stars and their protoplanetary disks: lessons learned on disk evolution and its connection to planet formation

Authors: C. F. Manara, M. Ansdell, G. P. Rosotti, A. M. Hughes, P. J. Armitage, G. Lodato, J. P. Williams

Abstract: Since Protostars and Planets VI (PPVI), our knowledge of the global properties of protoplanetary and debris disks, as well as of young stars, has dramatically improved. At the time of PPVI, mm-observations and optical to near-infrared spectroscopic surveys were largely limited to the Taurus star-forming region, especially of its most massive disk and stellar population. Now, near-complete surveys… ▽ More Since Protostars and Planets VI (PPVI), our knowledge of the global properties of protoplanetary and debris disks, as well as of young stars, has dramatically improved. At the time of PPVI, mm-observations and optical to near-infrared spectroscopic surveys were largely limited to the Taurus star-forming region, especially of its most massive disk and stellar population. Now, near-complete surveys of multiple star-forming regions cover both spectroscopy of young stars and mm interferometry of their protoplanetary disks. This provides an unprecedented statistical sample of stellar masses and mass accretion rates, as well as disk masses and radii, for almost 1000 young stellar objects within 300 pc from us, while also sampling different evolutionary stages, ages, and environments. At the same time, surveys of debris disks are revealing the bulk properties of this class of more evolved objects. This chapter reviews the statistics of these measured global star and disk properties and discusses their constraints on theoretical models describing global disk evolution. Our comparisons of observations to theoretical model predictions extends beyond the traditional viscous evolution framework to include analytical descriptions of magnetic wind effects. Finally, we discuss how recent observational results can provide a framework for models of planet population synthesis and planet formation. △ Less

Submitted 1 March, 2023; v1 submitted 18 March, 2022; originally announced March 2022.

Comments: Review Chapter for Protostars and Planets VII, Editors: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura. Final accepted version after interaction with the referees and community feedback. 21 pages (24 with references), 8 figures. Data table available at http://ppvii.org/chapter/15/