We overview various efforts within the DART Investigation Team’s Ejecta Working Group to predict the characteristics, quantity, dynamical behavior, and observability of DART impact ejecta. We discuss various methodologies for simulation of the impact/cratering process with their advantages and drawbacks in relation to initializing ejecta for subsequent dynamical propagation through and away from the Didymos system. We discuss the most relevant forces acting on ejecta once decoupled from Dimorphos’s surface and highlight various software packages we have developed and used to dynamically simulate ejecta under the action of those forces. With some additional software packages, we explore the influence of additional perturbing effects, such as interparticle collisions within true N-body codes and nonspherical and rotating particles’ interplay with solar radiation pressure. We find that early-timescale and close-proximity ejecta evolution is highly sensitive to some of these effects (e.g., collisions) while relatively insensitive to other factors. We present a methodology for turning the time-evolving size- and spatially discretized number density field output from ejecta simulations into synthetic images for multiple platforms/cameras over wide-ranging vantage points and timescales. We present such simulated images and apply preliminary analyses to them for nominal and off-nominal cases bracketing realistic total mass of ejecta and ejecta cumulative size–frequency distribution slope. Our analyses foreshadow the information content we may be able to extract from the actual images taken during and after the DART encounter by both LICIACube and Earth-vicinity telescopes.

We have monitored the Didymos–Dimorphos binary asteroid in spectropolarimetric mode in the optical range before and after the DART impact. The ultimate goal was to obtain constraints on the characteristics of the ejected dust for modeling purposes. Before impact, Didymos exhibited a linear polarization rapidly increasing with phase angle, reaching a level of ∼5% in the blue and ∼4.5% in the red. The shape of the polarization spectrum was anticorrelated with that of its reflectance spectrum, which appeared typical of an S-class asteroid. After impact, the level of polarization dropped by about 1 percentage point (pp) in the blue band and about 0.5 pp in the red band, then continued to linearly increase with phase angle, with a slope similar to that measured prior to impact. The polarization spectra, once normalized by their values at an arbitrary wavelength, show very little or no change over the course of all observations before and after impact. The lack of any remarkable change in the shape of the polarization spectrum after impact suggests that the way in which polarization varies with wavelength depends on the composition of the scattering material, rather than on its structure, be this a surface or a debris cloud.

We present Hubble Space Telescope Near-Infrared Camera and Multi-Object Spectrometer coronagraphic imaging polarimetry of the TW Hydrae protoplanetary disk. These observations simultaneously measure the total and polarized intensity, allowing direct measurement of the polarization fraction across the disk. In accord with the self-shadowing hypothesis recently proposed by Debes et al., we find that the total and polarized intensity of the disk exhibits strong azimuthal asymmetries at projected distances consistent with the previously reported bright and dark ring-shaped structures (∼45-99 au). The sinusoidal-like variations possess a maximum brightness at position angles near ∼268-300 • and are up to ∼28% stronger in total intensity. Furthermore, significant radial and azimuthal variations are also detected in the polarization fraction of the disk. In particular, we find that regions of lower polarization fraction are associated with annuli of increased surface brightness, suggesting that the relative proportion of multiple-to-single scattering is greater along the ring and gap structures. Moreover, we find strong (∼20%) azimuthal variation in the polarization fraction along the shadowed region of the disk. Further investigation reveals that the azimuthal variation is not the result of disk flaring effects, but instead from a decrease in the relative contribution of multiple-to-single scattering within the shadowed region. Employing a two-layer scattering surface, we hypothesize that the diminished contribution in multiple scattering may result from shadowing by an inclined inner disk, which prevents direct stellar light from reaching the optically thick underlying surface component.

This paper is dedicated to a long-standing problem of the shape of the negative branch of polarization (NBP) for Jupiter's moon Europa, determination of which is crucial for the characterization of the icy regolith on this satellite and similar objects, as well as for further progress in understanding light scattering by particulate surfaces. To establish the shape of Europa's NBP, in 2018–2021 we accomplished high-precision disk-integrated polarimetry of Europa in the UBVR I bands using the identical two-channel photoelectric polarimeters mounted on the 2.6 m Shajn reflector of the Crimean Astrophysical Observatory and the 2 m telescope of the Peak Terskol Observatory. We found that the polarization dependence on the phase angle in each filter is an asymmetric curve with a sharp polarization minimum P min ≈ − 0.3 % at phase angle α min ≤ 0 .° 4 , after which the polarization degree gradually increases to positive values, passing the inversion angle at α inv ≈ 6° − 7°. Within the error limits, the parameters P min, α min , and α inv of the NPB are independent of the wavelength in the visible spectrum. The polarization curve clearly demonstrates the so-called polarization opposition effect (POE). Our analysis of the previous and new polarimetric observations of Europa allows us to conclude that the POE is caused by coherent backscattering of sunlight on microscopic icy grains covering Europa’s surface. Computer modeling with the numerical radiative transfer coherent backscattering method demonstrates that the best fit to the polarimetric observations and geometric albedo of Europa is provided by a regolith layer of elementary single-scattering albedo ∼0.985 and extinction mean free path length 2π l/λ eff ≈ 150, λ eff representing the effective wavelength in the UBVR I spectral bands.

Context. Comets are conglomerates of ice and dust particles, the latter of which encode information on changes in the radiative and thermal environments. Dust displays distinctive scattered and thermal radiation in the visible and mid-infrared (MIR) wavelengths, respectively, based on its inherent characteristics. Aims. We aim to identify a possible correlation between the properties of scattered and thermal radiation from dust and the principal dust characteristics responsible for this relationship, and therefrom gain insights into comet evolution. Methods. We use the NASA/PDS archival polarimetric data on cometary dust in the red (0.62−0.73 μm) and K (2.00−2.39 μm) domains to leverage the relative excess of the polarisation degree of a comet to the average trend at the given phase angle (Pexcess) as a metric of the dust’s scattered light characteristics. The flux excess of silicate emissions to the continuum around 10 μm (FSi/Fcont) is adopted from previous studies as a metric of the dust’s MIR feature. Results. The two observables – Pexcess and FSi/Fcont – show a positive correlation when Pexcess is measured in the K domain (Spearman’s rank correlation coefficient ρ = 0.71−0.19+0.10). No significant correlation was identified in the red domain (ρ = 0.13−0.15+0.16). The gas-rich comets have systematically weaker FSi/Fcont than the dust-rich ones, and yet both groups retain the same overall tendency with different slope values. Conclusions. The observed positive correlation between the two metrics indicates that composition is a peripheral factor in characterising the dust’s polarimetric and silicate emission properties. The systematic difference in FSi/Fcont for gas-rich versus dust-rich comets would instead correspond to the difference in their dust size distribution. Hence, our results suggest that the current MIR spectral models of cometary dust, which search for a minimum χ2 fit by considering various dust properties simultaneously, should prioritise the dust size and porosity over the composition. With light scattering being sensitive to different size scales in two wavebands, we expect the K-domain polarimetry to be sensitive to the properties of dust aggregates, such as size and porosity, which might have been influenced by evolutionary processes. On the other hand, the red-domain polarimetry reflects the characteristics of sub-micrometre constituents in the aggregate.