Planetary Architectures Research Articles

Direct imaging of exoplanets: Legacy and prospects

Understanding how giant and terrestrial planets form and evolve, what is their internal structure and that of their atmosphere, represents one of the major challenges of modern astronomy, which is directly connected to the ultimate search for life at the horizon 2030–2050. However, several astrophysical (understanding of the formation and physics of giant and terrestrial exoplanets), biological (identification of the best biomarkers) and technological (technical innovations for the new generations of telescopes and instruments) obstacles must be overcome. From the astrophysical point of view, it is indeed crucial to understand the mechanisms of formation and evolution of giant planets, including planet and disk interactions, which will completely sculpt the planetary architectures and thus dominate the formation of terrestrial planets, especially in regions around the host star capable of supporting life. It is also important to develop dedicated instrumentation and techniques to study in their totality the population of giant and terrestrial planets, but also to reveal in the near future the first biological markers of life in the atmospheres of terrestrial planets. In that perspective, direct imaging from ground-based observatories or in space is playing a central role in concert with other observing techniques. In this paper, I will briefly review the genesis of this observing technique, the main instrumental innovation and challenges, stellar targets and surveys, to then present the main results obtained so far about the physics and the mechanisms of formation and evolution of young giant planets and planetary system architectures. I will then present the exciting perspectives offered by the upcoming generation of planet imagers about to come online, particularly on the future extremely large telescopes. On the timescale of a human Life, we may well be witnessing the first discovery of an exoplanet and the first detection of indices of life in the atmosphere of a nearby exo-Earth!

Friends Not Foes: Strong Correlation between Inner Super-Earths and Outer Gas Giants

The connection between outer gas giants and inner super-Earths reflects their formation and evolutionary histories. Past work exploring this link has suggested a tentative positive correlation between these two populations, but these studies have been limited by small sample sizes and in some cases sample biases. Here we take a new look at this connection with a sample of 184 super-Earth systems with publicly available radial velocity data and resolved outer gas giants. We calculate the frequency of outer gas giants (GG) in super-Earth (SE) systems, dividing our sample into metal-rich ([Fe/H] > 0) and metal-poor ([Fe/H]≤0) hosts. We find P(GG∣SE, [Fe/H]>0) = 28.0−4.6+4.9% and P(GG∣SE, [Fe/H]≤0) = 4.5−1.9+2.6% . Comparing these conditional occurrence rates to field giant occurrence rates from Rosenthal et al. (2021), we show that there is a distinct positive correlation between inner super-Earths and outer gas giants for metal-rich host stars at the 2.7σ level, but this correlation disappears for metal-poor systems. We further find that, around metal-rich stars, the GG/SE correlation enhances slightly for systems with giants that are more distant (beyond 3 au), more eccentric (e > 0.2), and/or in multi-gas giant systems. Such trends disappear around metal-poor stars with the exception of systems of multiple giants in which we observe a tentative anti-correlation. Our findings highlight the critical role metallicity (disk solid budget) plays in shaping the overall planetary architecture.

White dwarf pollution: one star or two?

ABSTRACT The accretion of tidally disrupted planetary bodies is the current consensus model for the presence of photospheric metals commonly detected in white dwarfs. While most dynamical studies have considered a single star and associated planetary instabilities, several investigations have instead considered the influence of widely bound stellar companions as potential drivers of white dwarf pollution. This study examines the prevalence of wide binaries among polluted white dwarfs using Gaia Data Release 3 astrometry, where three samples are investigated: 71 DAZ stars with metals detected in the ultraviolet using Hubble Space Telescope, and two groups of DZ stars identified via Sloan Digital Sky Survey spectroscopy, comprised of 116 warmer and 101 cooler sources. Each sample was searched for spatially resolved, comoving companions, and compared to the same analysis of thousands of field white dwarfs within overlapping regions of the Gaia Hertzsprung–Russell diagram. The wide binary fraction of the DAZ sample is $10.6_{-3.2}^{+3.9}$ per cent, and within $1\sigma$ of the corresponding field. However, the search yields wide binary fractions of less than 1.8 per cent for the two independent DZ star catalogues, which are each distinct from their fields by more than $3\sigma$. Both sets of results support that pollution in white dwarfs is not the result of stellar companions, and the delivery of metals to white dwarf surfaces is caused by major planets. The discrepancy between the DAZ and DZ star wide binary fractions cannot be caused by white dwarf spectral evolution, suggesting these two populations may have distinct planetary architectures.

Dynamical Architectures of S-type Transiting Planets in Binaries. I. Target Selection Using Hipparcos and Gaia Proper Motion Anomalies* * Also referred to as Hipparcos and Gaia astrometric accelerations.

The effect of stellar multiplicity on planetary architecture and orbital dynamics provides an important context for exoplanet demographics. We present a volume-limited catalog of up to 300 pc of 66 stars hosting planets and planet candidates from Kepler, K2, and TESS with significant Hipparcos-Gaia proper motion anomalies, which indicates the presence of companions. We assess the reliability of each transiting planet candidate using ground-based follow-up observations, and find that the TESS Objects of Interest (TOIs) with significant proper anomalies show nearly four times more false positives due to eclipsing binaries compared to TOIs with marginal proper anomalies. In addition, we find tentative evidence that orbital periods of planets orbiting TOIs with significant proper anomalies are shorter than those orbiting TOIs without significant proper anomalies, consistent with the scenario that stellar companions can truncate planet-forming disks. Furthermore, TOIs with significant proper anomalies exhibit lower Gaia differential velocities in comparison to field stars with significant proper anomalies, suggesting that planets are more likely to form in binary systems with low-mass substellar companions or stellar companions at wider separation. Finally, we characterize the three-dimensional architecture of LTT 1445 ABC using radial velocities, absolute astrometry from Gaia and Hipparcos, and relative astrometry from imaging. Our analysis reveals that LTT 1445 is a nearly flat system, with a mutual inclination of ∼2.°88 between the orbit of BC around A and that of C around B. This coplanarity may explain why multiple planets around LTT 1445 A survive in the dynamically hostile environments of this system.

Earths Are Not Super-Earths, Saturns Are Not Jupiters: Imprints of Pressure-bump Planet Formation on Planetary Architectures

In protoplanetary disks, sufficiently massive planets excite pressure bumps, which can then be preferred locations for forming new planet cores. We discuss how this loop may affect the architecture of multiplanet systems and compare our predictions with observations. Our main prediction is that low-mass planets and giant planets can each be divided into two subpopulations with different levels of mass uniformity. Low-mass planets that can and cannot reach the pebble isolation mass (the minimum mass required to produce a pressure bump) develop into intra-system similarity “super-Earths” and more diverse “Earths,” respectively. Gas giants that do and do not accrete envelopes quickly develop into similar “Jupiters” and more diverse “Saturns,” respectively. Super-Earths prefer to form long chains via repeated pressure-bump planet formation, while Jupiter formation is usually terminated at pairs or triplets due to dynamical instability. These predictions are broadly consistent with observations. In particular, we discover a previously overlooked mass uniformity dichotomy among the observed populations of both low-mass planets (Earths versus super-Earths) and gas giants (Saturns versus Jupiters). For low-mass planets, planets well below the pebble isolation mass (≲3 M ⊕ or ≲1.5 R ⊕ for Sun-like stars) show significantly higher intra-system pairwise mass differences than planets around the pebble isolation mass. For gas giants, the period ratios of intra-system pairs show a bimodal distribution, which can be interpreted as two subpopulations with different levels of mass uniformity. These findings suggest that pressure-bump planet formation could be an important ingredient in shaping planetary architectures.

Planetary companions orbiting the M dwarfs GJ 724 and GJ 3988

We report the discovery of two exoplanets around the M dwarfs GJ 724 and GJ 3988 using the radial velocity (RV) method. We obtained a total of 153 3.5 m Calar Alto/CARMENES spectra for both targets and measured their RVs and activity indicators. We also added archival ESO/HARPS data for GJ 724 and infrared RV measurements from Subaru/IRD for GJ 3988. We searched for periodic and stable signals to subsequently construct Keplerian models, considering different numbers of planets, and we selected the best models based on their Bayesian evidence. Gaussian process (GP) regression was included in some models to account for activity signals. For both systems, the best model corresponds to one single planet. The minimum masses are 10.75−0.87+0.96 and 3.69−0.41+0.42 Earth-masses for GJ 724 b and GJ 3988 b, respectively. Both planets have short periods (P < 10 d) and, therefore, they orbit their star closely (a < 0.05 au). GJ 724 b has an eccentric orbit (e = 0.577−0.052+0.055), whereas the orbit of GJ 3988 b is circular. The high eccentricity of GJ 724 b makes it the most eccentric single exoplanet (to this date) around an M dwarf. Thus, we suggest a further analysis to understand its configuration in the context of planetary formation and architecture. In contrast, GJ 3988 b is an example of a common type of planet around mid-M dwarfs.

Planets Across Space and Time (PAST). IV. The Occurrence and Architecture of Kepler Planetary Systems as a Function of Kinematic Age Revealed by the LAMOST–Gaia–Kepler Sample

One of the fundamental questions in astronomy is how planetary systems form and evolve. Measuring the planetary occurrence and architecture as a function of time directly addresses this question. In the fourth paper of the Planets Across Space and Time series, we investigate the occurrence and architecture of Kepler planetary systems as a function of kinematic age by using the LAMOST–Gaia–Kepler sample. To isolate the age effect, other stellar properties (e.g., metallicity) have been controlled. We find the following results. (1) The fraction of stars with Kepler-like planets (F Kep) is about 50% for all stars; no significant trend is found between F Kep and age. (2) The average planet multiplicity () exhibits a decreasing trend (∼2σ significance) with age. It decreases from ∼ 3 for stars younger than 1 Gyr to ∼ 1.8 for stars of about 8 Gyr. (3) The number of planets per star (η = F Kep× ) also shows a decreasing trend (∼2σ–3σ significance). It decreases from η ∼ 1.6–1.7 for young stars to η ∼ 1.0 for old stars. (4) The mutual orbital inclination of the planets (σ i,k ) increases from to as the stars age from 0.5 to 8 Gyr with a best fit of . Interestingly, the solar system also fits such a trend. The fact that F Kep remains relatively constant at approximately ∼ 50% across different ages suggests the robustness of planet formation throughout the history of the Galaxy. The age dependence of and σ i,k demonstrates that the planetary architecture is evolving, and planetary systems generally become dynamically hotter with fewer planets as they age.

The GAPS programme at TNG

Context. Observational data from space- and ground-based campaigns have revealed that the 10-30 Ma old V1298Tau star hosts a compact and massive system of four planets. Mass estimates are available for the two outer giant planets and point to unexpectedly high densities for their young ages. Aims. We investigate the formation of these two outermost giant planets, V1298 Tau b and e, and the present dynamical state of V1298 Tau’s global architecture in order to shed light on the history of this young and peculiar extrasolar system. Methods. We performed detailed N-body simulations to explore the link between the densities of V1298 Tau b and e and their migration and accretion of planetesimals within the native circumstellar disk. We combined N-body simulations and the normalized angular momentum deficit (NAMD) analysis of the architecture to characterize V1298 Tau’s dynamical state and connect it to the formation history of the system. We searched for outer planetary companions to constrain V1298 Tau’s planetary architecture and the extension of its primordial circumstellar disk. Results. The high densities of V1298 Tau b and e suggest they formed at quite a distance from their host star, likely beyond the CO2 snowline. The higher nominal density of V1298 Tau e suggests it formed farther out than V1298 Tau b. The current architecture of V1298 Tau is not characterized by resonant chains. Planet-planet scattering with an outer giant planet is the most likely cause for the lack of a resonant chain between V1298 Tau’s planets, but currently our search for outer companions using SPHERE and Gaia observations can exclude only the presence of planets more massive than 2 MJ. Conclusions. The most plausible scenario for V1298 Tau’s formation is that the system formed by convergent migration and resonant trapping of planets born in a compact and plausibly massive disk. In the wake of their migration, V1298 Tau b and e would have left a dynamically excited protoplanetary disk, naturally creating the conditions for the later breaking of the resonant chain by planet-planet scattering.

Can comets deliver prebiotic molecules to rocky exoplanets?

In this work, we consider the potential of cometary impacts to deliver complex organic molecules and the prebiotic building blocks required for life to rocky exoplanets. Numerical experiments have demonstrated that for these molecules to survive, impacts at very low velocities are required. This work shows that for comets scattered from beyond the snow-line into the habitable zone, the minimum impact velocity is always lower for planets orbiting Solar-type stars than M-dwarfs. Using both an analytical model and numerical N-body simulations, we show that the lowest velocity impacts occur onto planets in tightly packed planetary systems around high-mass (i.e. Solar-mass) stars, enabling the intact delivery of complex organic molecules. Impacts onto planets around low-mass stars are found to be very sensitive to the planetary architecture, with the survival of complex prebiotic molecules potentially impossible in loosely packed systems. Rocky planets around M-dwarfs also suffer significantly more high velocity impacts, potentially posing unique challenges for life on these planets. In the scenario that cometary delivery is important for the origins of life, this study predicts the presence of biosignatures will be correlated with (i) decreasing planetary mass (i.e. escape velocity), (ii) increasing stellar-mass and (iii) decreasing planetary separation (i.e. exoplanets in tightly-packed systems).

Refining the Stellar Parameters of τ Ceti: a Pole-on Solar Analog

To accurately characterize the planets a star may be hosting, stellar parameters must first be well determined. τ Ceti is a nearby solar analog and often a target for exoplanet searches. Uncertainties in the observed rotational velocities have made constraining τ Ceti’s inclination difficult. For planet candidates from radial velocity (RV) observations, this leads to substantial uncertainties in the planetary masses, as only the minimum mass () can be constrained with RV. In this paper, we used new long-baseline optical interferometric data from the CHARA Array with the MIRC-X beam combiner and extreme precision spectroscopic data from the Lowell Discovery Telescope with EXPRES to improve constraints on the stellar parameters of τ Ceti. Additional archival data were obtained from a Tennessee State University Automatic Photometric Telescope and the Mount Wilson Observatory HK project. These new and archival data sets led to improved stellar parameter determinations, including a limb-darkened angular diameter of 2.019 ± 0.012 mas and rotation period of 46 ± 4 days. By combining parameters from our data sets, we obtained an estimate for the stellar inclination of 7° ± 7°. This nearly pole-on orientation has implications for the previously reported exoplanets. An analysis of the system dynamics suggests that the planetary architecture described by Feng et al. may not retain long-term stability for low orbital inclinations. Additionally, the inclination of τ Ceti reveals a misalignment between the inclinations of the stellar rotation axis and the previously measured debris disk rotation axis (i disk = 35° ± 10°).

Chasing extreme planetary architectures

Context. Planet(s) in binaries are unique architectures for testing predictions of planetary formation and evolution theories in very hostile environments. Their presence in such a highly perturbed region poses a clear challenge from the early phase of planetesimals accretion to the dynamical evolution and stability through a very long period of time (several Gyrs in some case). Aims. The combination of radial velocity, speckle interferometry, high-contrast imaging and high-precision astrometry with interfer-ometry, offers a unique and unprecedented set of observable to push the exploration of the extreme planetary system HD 196885. By dissecting the physical and orbital properties of each component, we aim at shedding light on its global architecture and stability. Methods. We used the IRDIS dual-band imager of SPHERE at VLT, and the speckle interferometric camera HRCAM of SOAR, to acquire high-angular resolution images of HD 196885 AB between 2015 and 2020. Radial velocity observations started in 1983 with CORAVEL at OHP have been extended with a continuous monitoring with CORALIE at La Silla, and ELODIE and SOPHIE at OHP over almost 40 yr extending the radial velocity measurements HD 196885 A and resolving both the binary companion and the inner giant planet HD 196885 Ab. Finally, we took advantage of the exquisite astrometric precision of the dual-field mode of VLTI/GRAVITY (down to 30 µas) to monitor the relative position of HD 196885 A and B to search for the 3.6 yr astrometric wobble of the circumpri-mary planet Ab imprinted on the binary separation. Results. Our observations enable to accurately constrain the orbital properties of the binary HD 196885 AB, seen on an inclined and retrograde orbit (iAB = 120.43 deg) with a semi-major axis of 19.78 au, and an eccentricity of 0.417. The GRAVITY measurements confirm for the first time the nature of the inner planet HD 196885 Ab by rejecting all families of pole-on solutions in the stellar or brown dwarf masses. The most favored island of solutions is associated with a Jupiter-like planet (MAb = 3.39 MJup), with moderate eccentricity (eAaAb = 0.44), and inclination close to 143.04 deg. This results points toward a significant mutual inclination (Ф = 24.36 deg) between the orbital planes (relative to the star) of the binary companion B and the planet Ab. Our dynamical simulations indicate that the system is dynamically stable over time. Eccentricity and mutual inclination variations could be expected for moderate von Zipele Kozai Lidov cycles that may affect the inner planet.

Planet formation via pebble accretion in externally photoevaporating discs

ABSTRACT We demonstrate that planet formation via pebble accretion is sensitive to external photoevaporation of the outer disc. In pebble accretion, planets grow by accreting from a flux of solids (pebbles) that radially drift inwards from the pebble production front. If external photoevaporation truncates the outer disc fast enough, it can shorten the time before the pebble production front reaches the disc outer edge, cutting off the supply of pebble flux for accretion, hence limiting the pebble mass reservoir for planet growth. Conversely, cloud shielding can protect the disc from strong external photoevaporation and preserve the pebble reservoir. Because grain growth and drift can occur quickly, shielding even on a short time-scale (<1 Myr) can have a non-linear impact on the properties of planets growing by pebble accretion. For example, a $10^{-3}\, \mathrm{ M}_{\oplus }$ planetary seed at 25 au stays at 25 au with a lunar mass if the disc is immediately irradiated by a 103 G0 field, but grows and migrates to be approximately Earth-like in both mass and orbital radius if the disc is shielded for just 1 Myr. In NGC 2024, external photoevaporation is thought to happen to discs that are <0.5 Myr old, which coupled with the results here suggests that the exact planetary parameters can be very sensitive to the star-forming environment. Universal shielding for time-scales of at least ${\sim} 1.5\,$ Myr would be required to completely nullify the environmental impact on planetary architectures.

The Dynamical Consequences of a Super-Earth in the Solar System

Placing the architecture of the solar system within the broader context of planetary architectures is one of the primary topics of interest within planetary science. Exoplanet discoveries have revealed a large range of system architectures, many of which differ substantially from the solar system’s model. One particular feature of exoplanet demographics is the relative prevalence of super-Earth planets, for which the solar system lacks a suitable analog, presenting a challenge to modeling their interiors and atmospheres. Here we present the results of a large suite of dynamical simulations that insert a hypothetical planet in the mass range 1–10 M ⊕ within the semimajor axis range 2–4 au, between the orbits of Mars and Jupiter. We show that, although the system dynamics remain largely unaffected when the additional planet is placed near 3 au, Mercury experiences substantial instability when the additional planet lies in the range 3.1–4.0 au, and perturbations to the Martian orbit primarily result when the additional planet lies in the range 2.0–2.7 au. We further show that, although Jupiter and Saturn experience relatively small orbital perturbations, the angular momentum transferred to the ice giants can result in their ejection from the system at key resonance locations of the additional planet. We discuss the implications of these results for the architecture of the inner and outer solar system planets, and for exoplanetary systems.

Finding order in planetary architectures

Finding order in planetary architectures

Formation of rocky super-earths from a narrow ring of planetesimals

The formation of super-Earths, the most abundant planets in the Galaxy, remains elusive. These planets have masses that typically exceed that of the Earth by a factor of a few, appear to be predominantly rocky, although often surrounded by H/He atmospheres, and frequently occur in multiples. Moreover, planets that encircle the same star tend to have similar masses and radii, whereas those belonging to different systems exhibit remarkable overall diversity. Here we advance a theoretical picture for rocky planet formation that satisfies the aforementioned constraints: building upon recent work, which has demonstrated that planetesimals can form rapidly at discrete locations in the disk, we propose that super-Earths originate inside rings of silicate-rich planetesimals at approximately ~1 au. Within the context of this picture, we show that planets grow primarily through pairwise collisions among rocky planetesimals until they achieve terminal masses that are regulated by isolation and orbital migration. We quantify our model with numerical simulations and demonstrate that our synthetic planetary systems bear a close resemblance to compact, multi-resonant progenitors of the observed population of short-period extrasolar planets. A planetary origin model that forms exoplanets from a narrow ring of silicate material at a stellocentric distance of 1 au is able to explain the physical properties of super-Earths and reproduce the ‘peas in a pod’ pattern of uniformity within planetary architecture.

Prediction of an Earth-sized Planet Formed in the Habitable Zone of the SPECULOOS-2 System

Transiting planets in the nearby TRAPPIST-1 system provide rare examples of habitable zone (HZ), Earth-sized planets that can be characterized via transmission spectroscopy. However, these present-day HZ planets likely formed interior to HZ and probably underwent an evolution very different from that of Earth. We present the integrative analysis of the planetary architecture of the recently discovered planetary system SPECULOOS-2. Our analysis answers the question: If there are additional exoplanets in the SPECULOOS-2 system, what are their orbital and physical properties? We predict an Earth-sized planet in the habitable zone (P ∼ 14.5–18.2 days). In contrast to TRAPPIST-1e, this predicted planet most likely completed its formation inside the habitable zone. If confirmed, this planet will offer an Earth-sized, habitable zone planet that is likely to have an evolutionary path more similar to Earth than those in the TRAPPIST-1 system.

System Architecture and Planetary Obliquity: Implications for Long-term Habitability

In the search for life beyond our solar system, attention should be focused on those planets that have the potential to maintain habitable conditions over the prolonged periods of time needed for the emergence and expansion of life as we know it. The observable planetary architecture is one of the determinants for long-term habitability as it controls the orbital evolution and ultimately the stellar fluxes received by the planet. With an ensemble of n-body simulations and obliquity models of hypothetical planetary systems, we demonstrate that the amplitude and period of the eccentricity, obliquity, and precession cycles of an Earth-like planet are sensitive to the orbital characteristics of a giant companion planet. A series of transient, ocean-coupled climate simulations show how these characteristics of astronomical cycles are decisive for the evolving surface conditions and long-term fractional habitability relative to the modern Earth. The habitability of Earth-like planets increases with the eccentricity of a Jupiter-like companion, provided that the mean obliquity is sufficiently low to maintain temperate temperatures over large parts of its surface throughout the orbital year. A giant companion closer in results in shorter eccentricity cycles of an Earth-like planet but longer, high-amplitude, obliquity cycles. The period and amplitude of obliquity cycles can be estimated to first order from the orbital pathways calculated by the n-body simulations. In the majority of simulations, the obliquity amplitude relates directly to the orbital inclination whereas the period of the obliquity cycle is a function of the nodal precession and the proximity of the giant companion.

The Hot Neptune WASP-166 b with ESPRESSO – I. Refining the planetary architecture and stellar variability

ABSTRACT In this paper, we present high-resolution spectroscopic transit observations from ESPRESSO of the super-Neptune WASP-166 b. In addition to spectroscopic ESPRESSO data, we analyse photometric data from TESS of six WASP-166 b transits along with simultaneous NGTS observations of the ESPRESSO runs. These observations were used to fit for the planetary parameters as well as assessing the level of stellar activity (e.g. spot crossings, flares) present during the ESPRESSO observations. We utilize the reloaded Rossiter McLaughlin (RRM) technique to spatially resolve the stellar surface, characterizing the centre-to-limb convection-induced variations, and to refine the star–planet obliquity. We find WASP-166 b has a projected obliquity of $\lambda = -15.52^{+2.85}_{-2.76}\, ^{\circ }$ and vsin (i) = 4.97 ± 0.09 km s−1 which is consistent with the literature. We were able to characterize centre-to-limb convective variations as a result of granulation on the surface of the star on the order of a few km s−1 for the first time. We modelled the centre-to-limb convective variations using a linear, quadratic, and cubic model with the cubic being preferred. In addition, by modelling the differential rotation and centre-to-limb convective variations simultaneously, we were able to retrieve a potential antisolar differential rotational shear (α ∼ −0.5) and stellar inclination (i* either 42.03$^{+9.13}_{-9.60}\, ^{\circ }$ or 133.64$^{+8.42}_{-7.98}\, ^{\circ }$ if the star is pointing towards or away from us). Finally, we investigate how the shape of the cross-correlation functions change as a function of limb angle and compare our results to magnetohydrodynamic simulations.

Trends in Silicates in the β Pictoris Disk

While β Pic is known to host silicates in ring-like structures, whether the properties of these silicate dust vary with stellocentric distance remains an open question. We re-analyze the β Pictoris debris disk spectrum from the Spitzer Infrared Spectrograph (IRS) and a new Infrared Telescope Facility Spectrograph and Imager spectrum to investigate trends in Fe/Mg ratio, shape, and crystallinity in grains as a function of wavelength, a proxy for stellocentric distance. By analyzing a re-calibrated and re-extracted spectrum, we identify a new 18 μm forsterite emission feature and recover a 23 μm forsterite emission feature with a substantially larger line-to-continuum ratio than previously reported. We find that these prominent spectral features are primarily produced by small submicron-sized grains, which are continuously generated and replenished from planetesimal collisions in the disk and can elucidate their parent bodies’ composition. We discover three trends about these small grains: as stellocentric distance increases, (1) small silicate grains become more crystalline (less amorphous), (2) they become more irregular in shape, and (3) for crystalline silicate grains, the Fe/Mg ratio decreases. Applying these trends to β Pic’s planetary architecture, we find that the dust population exterior to the orbits of β Pic b and c differs substantially in crystallinity and shape. We also find a tentative 3–5 μm dust excess due to spatially unresolved hot dust emission close to the star. From our findings, we infer that the surfaces of large planetesimals are more Fe-rich and collisionally processed closer to the star but more Fe-poor and primordial farther from the star.

An Integrative Analysis of the Rich Planetary System of the Nearby Star e Eridani: Ideal Targets for Exoplanet Imaging and Biosignature Searches

e Eridani, the fifth-closest Sun-like star, hosts at least three planets and could possibly harbor more. However, the veracity of the planet candidates in the system and its full planetary architecture remain unknown. Here we analyze the planetary architecture of e Eridani via DYNAMITE, a method providing an integrative assessment of the system architecture (and possibly yet-undetected planets) by combining statistical, exoplanet-population-level knowledge with incomplete but specific information available on the system. DYNAMITE predicts the most likely location of an additional planet in the system based on the Kepler population demographic information from more than 2000 planets. Additionally, we analyze the dynamical stability of e Eridani system via N-body simulations. Our DYNAMITE and dynamical stability analyses provide support for planet candidates g, c, and f, and also predict one additional planet candidate with an orbital period between 549–733 days, in the habitable zone of the system. We find that planet candidate f, if it exists, would also lie in the habitable zone. Our dynamical stability analysis also shows that the e Eridani planetary eccentricities, as reported, do not allow for a stable system, suggesting that they are lower. We introduce a new statistical approach for estimating the equilibrium and surface temperatures of exoplanets, based on a prior from the planetary albedo distribution. e Eridani is a rich planetary system with a possibility of containing two potentially habitable planets, and its vicinity to our solar system makes it an important target for future imaging studies and biosignature searches.