Resonant structures in exozodiacal clouds created by exo-Earths in the habitable zone of late-type stars

TL;DR

The paper extends resonant-dust modeling to a broad spectral range (F4–M4) by incorporating spectral-type–dependent stellar wind drag alongside PR drag, focusing on Earth-like planets in the habitable zone. Using MERCURY6-based simulations, it demonstrates that resonant rings form across all spectral types, but stellar wind drag dominates for old M-dwarfs, reducing ring contrast and optical depth relative to wind-insensitive cases. Across spectral types, ring contrasts and multi-size optical depths generally rise for lower-mass stars, while asymmetric MIR flux peaks at ~10 μm are strongest for K-type stars due to a balance between ring strength and emitting area. The work provides semi-analytical fits to estimate resonant-ring strength and highlights implications for MIR interferometry missions like LIFE, emphasizing the need to account for spectral-type wind effects in exozodi models and target prioritization.

Abstract

Earth-like exoplanets can create resonant structures in exozodiacal dust through mean motion resonances (MMRs). These structures not only suggest the presence of such planets, but also act as potential noise sources in future mid-infrared (MIR) nulling interferometry observations. We aim to investigate how resonant structures in exozodiacal dust vary across stellar spectral types (F4--M4), and to evaluate how stellar wind drag affects their morphology and brightness in mature planetary systems. We conducted numerical simulations of dust dynamics, extending earlier studies by including spectral type variation in stellar wind drag in addition to Poynting-Robertson (PR) drag. Our models represented systems of a few Gyr hosting an Earth-like exoplanet in the habitable zone (HZ). We produced spatially resolved maps of optical depth and thermal emission for different stellar spectral types. Our simulations showed that resonant ring structures were formed for all stellar spectral types considered. In particular, we found that stellar wind drag played a critical role in shaping dust dynamics around old M-type stars, where it could dominate over PR drag by a factor of approximately 44. This reduced the contrast of resonant rings relative to the background disk, compared to cases without spectral type variation in stellar wind. Across different spectral types, the optical depth contrast of the resonant ring increased for lower-mass stars, assuming a fixed background level. Asymmetric thermal emission distributions were derived across all spectral types, which peaked for K-type stars. Our findings highlight the importance of incorporating both resonant dynamics and stellar wind effects when modeling exozodiacal dust around stars of different spectral types.

Figures (17)

• Figure 1: (Left) Example of face-on surface density distribution of dust in a planet's co-rotating frame, with a K4-type star located at the origin. The white dot indicates an Earth-like planet placed at the inner boundary of the CHZ. A total of 100 dust particles with a size of $50~\mathrm{\mu m}$ are used. The resonant ring structure is visible around the planet's position. The color scale indicates surface density in arbitrary units, with lighter colors corresponding to higher densities. (Right) Azimuthally averaged surface density of dust in arbitrary normalized units, based on the distribution in the left panel. The background disk appears as a nearly constant plateau, while the resonant ring structure shows a density enhancement.
• Figure 2: $\langle Q_\text{PR} \rangle$ values for grain sizes $s \sim 0.1$--$300~\mathrm{\mu m}$ (left) and the corresponding $\beta_\text{PR}$ values (right; analogous to Figure 7 of RW2020). Solid lines show numerical values for different stellar spectral types used in this study, while the dashed lines in the right panel indicate the analytical $\beta_\text{PR}$ values using $\langle Q_\text{PR} \rangle$ = 1 (Eq. \ref{['eq:beta_pr']}). Black dots mark the discrete grain sizes adopted in the main study. The black dotted line in the right panel denotes $\beta_\text{PR} = 0.5$, and sizes above this threshold, corresponding to grains blown out of the system for circular orbits, are marked with '$\times$'. The size $\sim 50 ~ \mathrm{\mu m}$ used in the pilot study is safely higher than the blowout size and meets the assumption of $\langle Q_\text{PR} \rangle \sim 1$.
• Figure 3: $\psi$ values given by Eqs. \ref{['eq: psi']} and \ref{['eq:massloss']} for main-sequence stars of spectral types F4 to M4 at an age of $\sim$ 4 Gyr (2--3 Gyr for F-type stars). The four spectral types used in this study are marked in blue, green, yellow, and red, respectively. The dashed line shows $\psi = 1$, where the stellar wind drag equals the PR drag.
• Figure 4: Contrast values from the pilot study, plotted against the $a_\text{p}^{1/2}/\beta_\text{PR}$ term for Model I (left) and Model II (right). Data corresponding to F4, G4, K4, and M4 are shown in different colors: blue, green, yellow, and red, respectively. Planetary masses of 0.5, 1.0, and 2.0 $M_\oplus$ are represented by different shapes, with solid lines indicating the fits from SK2008 for each $M_\text{p}$, increasing from bottom to top. For each color and symbol, the three points from left to right correspond to planets located at the inner, middle, and outer edges of the habitable zone. We note that these results are shown for $s = 50 ~ \mathrm{\mu m}$; the trends do not apply to grains smaller than $\sim 10 ~ \mathrm{\mu m}$ around F4, $\sim 3 ~ \mathrm{\mu m}$ around G4, and $\sim 1 ~ \mathrm{\mu m}$ around K4 and M4-type stars, which are blown out or fail to form resonant rings (see Section \ref{['sec:res_main_C']}).
• Figure 5: Contrast values from the main study plotted against grain sizes. Cases where $s < s_\text{BO,eff}$ are shown as $\times$ marks, which overlap with other data points up to $1~\mathrm{\mu m}$. Dashed lines in the upper panel indicate the size-combined contrast values of optical depth ($\langle C_\tau \rangle$; see Section \ref{['sec:res_tau_F']}) for each spectral type. Colors and symbols follow those used in Fig. \ref{['fig:contrast_pilot']}, but only cases with $M_\text{p} = 1~ M_\oplus$ and $a_\text{p} = \text{HZ}_\text{mid}$ are considered here, resulting in 8 points for each spectral type. The contrast for $300~\mathrm{\mu m}$ dust around the K-type star is only slightly lower than that for M-type stars, nearly overlapping and thus hidden from view.