An official website of the United States government
Abstract Modeling circumstellar disks in the traditional sense carries the assumption that the dust density distribution can be accurately described with a fixed parametric form. Furthermore, commonly used algorithms for subtracting the stellar point-spread function (PSF) distort the true morphology of the faint underlying disk structure, especially dusty features that are located at small angular separations. These phenomena often lead to significant residuals with parametric disk models and make it difficult to measure the full realizable range of the scattering function of the dust. We address these challenges withffortissimo, a novel, pixel-based freeform forward modeling pipeline designed to characterize extended objects in images reduced using the Karhunen–Loéve Image Projection algorithm. We built this pipeline within the framework ofJax, which is a machine learning library in Python that enables efficient optimization through automatic differentiation (“autodiff”) and graphics processing unit–accelerated array computations. Using visible light images of the disk around HR 4796A taken by the “extreme” Magellan Adaptive Optics instrument (MagAO-X), we show that our data-driven freeform models excel at fitting a complex dust distribution and can infer the dust scattering properties even through PSF subtraction artifacts. Additionally, we demonstrate the potential for retrieving spatial dust features beyond the diffraction limit of the telescope. We note that there are remaining challenges to address before precision photometry using these freeform models is advised. These include a better background, wind-driven halo, and speckle characterization, as preventing the freeform models from learning these noise artifacts is currently difficult.
more »
« less
This content will become publicly available on April 16, 2027
A Multiband Study of the HR 4796A Disk in the Optical Using MagAO-X
We present total intensity images of the debris disk around HR 4796A from observations spanning 2023 to 2025 with the Magellan extreme adaptive optics instrument (MagAO-X). We detected the disk at high signal-to-noise ratios at g′(527 nm), r′(615 nm), i′(762 nm), and z′(909 nm). Additionally, we present images collected using the “star-hopping” technique that show the entirety of the disk, including the dramatic forward-scattering at the minor axis. We subjected our images to a battery of modeling techniques to constrain the geometry and photometry of the disk. Leveraging our clear detections of the disk’s minor axis, we modeled the scattering phase function (SPF) using a basis of the Legendre polynomials. To mitigate self-subtraction artifacts in our angular differential imaging, we implemented a forward-modeling pipeline that generates a pixel-based freeform disk forward model leading to a deconvolved image of the disk. Our best-fit disk models reveal (1) highly forward-scattering SPFs with a minimum at the ∼65∘ scattering angle, (2) a faint halo of dust just exterior to the spine of the disk that is not well-described by a broken power-law density profile, (3) a red spectral slope for the dust, and finally (4) a compact, clump-like feature in the freeform disk models. Our empirically measured SPFs suggest that the scattering is dominated by large, highly absorptive grains. However, we emphasize the need for testing advanced irregular grain models using our SPFs to learn more about the physical and chemical properties of this complex system.
more »
« less
- Award ID(s):
- 2307612
- PAR ID:
- 10687929
- Author(s) / Creator(s):
- ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more »
- Publisher / Repository:
- American Astronomical Society
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 1001
- Issue:
- 2
- ISSN:
- 0004-637X
- Page Range / eLocation ID:
- 189
- Subject(s) / Keyword(s):
- Debris disks Adaptive optics High contrast techniques
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The HH 111 protostellar disk has recently been found to host a pair of spiral arms. Here we report the dust polarization results in the disk as well as the inner envelope around it, obtained with the Atacama Large Millimeter/submillimeter Array in continuum atλ∼ 870μm and ∼0.″05 resolution. In the inner envelope, polarization is detected with a polarization degree of ∼6% and an orientation almost everywhere parallel to the minor axis of the disk and thus likely to be due to the dust grains magnetically aligned mainly by toroidal fields. In the disk, the polarization orientation is roughly azimuthal on the far side and becomes parallel to the minor axis on the near side, with a polarization gap in between on the far side near the central protostar. The disk polarization degree is ∼2%. The polarized intensity is higher on the near side than the far side, showing a near–far side asymmetry. More importantly, the polarized intensity and thus polarization degree are lower in the spiral arms but higher in between the arms, showing an anticorrelation of the polarized intensity with the spiral arms. Our modeling results indicate that this anticorrelation is useful for constraining the polarization mechanism and is consistent with the dust self-scattering by the grains that have grown to a size of ∼150μm. The interarms are sandwiched and illuminated by two brighter spiral arms and thus have higher polarized intensity. Our dust self-scattering model can also reproduce the observed polarization orientation parallel to the minor axis on the near side and the observed azimuthal polarization orientation at the two disk edges in the major axis. Further modeling work is needed to study how to reproduce the observed near–far side asymmetry in the polarized intensity and the observed azimuthal polarization orientation on the far side.more » « less
-
Context. T Tauri stars are low-mass young stars whose disks provide the setting for planet formation. Despite this, their structure is poorly understood. We present new infrared interferometric observations of the SU Aurigae circumstellar environment that offer resolution that is three times higher and a better baseline position angle coverage than previous observations. Aims. We aim to investigate the characteristics of the circumstellar material around SU Aur, constrain the disk geometry, composition and inner dust rim structure. Methods. The CHARA array offers unique opportunities for long baseline observations, with baselines up to 331 m. Using the CLIMB three-telescope combiner in the K -band allows us to measure visibilities as well as closure phase. We undertook image reconstruction for model-independent analysis, and fitted geometric models such as Gaussian and ring distributions. Additionally, the fitting of radiative transfer models constrain the physical parameters of the disk. For the first time, a dusty disk wind is introduced to the radiative transfer code TORUS to model protoplanetary disks. Our implementation is motivated by theoretical models of dusty disk winds, where magnetic field lines drive dust above the disk plane close to the sublimation zone. Results. Image reconstruction reveals an inclined disk with slight asymmetry along its minor-axis, likely due to inclination effects obscuring the inner disk rim through absorption of incident star light on the near-side and thermal re-emission and scattering of the far-side. Geometric modelling of a skewed ring finds the inner rim at 0.17 ± 0.02 au with an inclination of 50.9 ± 1.0° and minor axis position angle 60.8 ± 1.2°. Radiative transfer modelling shows a flared disk with an inner radius at 0.18 au which implies a grain size of 0.4 μ m assuming astronomical silicates and a scale height of 15.0 at 100 au. Among the tested radiative transfer models, only the dusty disk wind successfully accounts for the K -band excess by introducing dust above the mid-plane.more » « less
-
Abstract We performed radiative transfer calculations and observing simulations to reproduce the 1.3 mm dust-continuum and C18O (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 a dust disk model passively heated by the central protostar cannot reproduce the observed peak brightness temperature of the 1.3 mm continuum emission (∼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 far side of the disk. Our modeling reveals that this asymmetric intensity distribution requires flaring of the dust along the disk vertical direction with the scale height followingh/r∼r0.3as a 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 C18O (2–1) emission, the outermost radius of the gas disk is estimated to be ∼80 au, which is larger than that of the dust disk (∼62 au), to reproduce the observed distribution of the C18O (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.more » « less
-
Abstract We present 870 μ m Atacama Large Millimeter/submillimeter Array polarization observations of thermal dust emission from the iconic, edge-on debris disk β Pic. While the spatially resolved map does not exhibit detectable polarized dust emission, we detect polarization at the ∼3 σ level when averaging the emission across the entire disk. The corresponding polarization fraction is P frac = 0.51% ± 0.19%. The polarization position angle χ is aligned with the minor axis of the disk, as expected from models of dust grains aligned via radiative alignment torques (RAT) with respect to a toroidal magnetic field ( B -RAT) or with respect to the anisotropy in the radiation field ( k -RAT). When averaging the polarized emission across the outer versus inner thirds of the disk, we find that the polarization arises primarily from the SW third. We perform synthetic observations assuming grain alignment via both k -RAT and B -RAT. Both models produce polarization fractions close to our observed value when the emission is averaged across the entire disk. When we average the models in the inner versus outer thirds of the disk, we find that k -RAT is the likely mechanism producing the polarized emission in β Pic. A comparison of timescales relevant to grain alignment also yields the same conclusion. For dust grains with realistic aspect ratios (i.e., s > 1.1), our models imply low grain-alignment efficiencies.more » « less
