Academia Sinica Logo

Understanding orbital evolution and planet formation through the lenses of N-body simulations is crucial for chaotic dynamical systems. On the other hand, atmospheric and climate modeling can provide us with unique insights into atmospheric evolution, habitability, and the observation prospects of rocky extrasolar planets. In this talk, I will present results from two separate studies combining techniques of N-body simulations in conjunction with atmospheric and climate modeling. In the first study, we used a three-dimensional N-Rigid-Body integrator and an intermediately-complex general circulation model to simulate the evolving climates of TRAPPIST-1 e and f with different orbital and spin evolution pathways. We find that sporadic libration and rotation induced by planetary interactions, such as that due to mean motion resonances (MMRs) in compact planetary systems may destabilize attendant exoplanets away from synchronized states (or 1:1 spin-orbit ratio). Planet f perturbed by MMR effects with sporadic spin-variations are colder and dryer compared to their synchronized counterparts due to the zonal drift of the substellar point away from open ocean basins of their initial eyeball states. On the other hand, the differences between perturbed and synchronized planet e are minor due to higher instellation, warmer surfaces, and reduced climate hysteresis. This is the first study to incorporate the time-dependent outcomes of direct gravitational N-Rigid-Body simulations into 3D climate modeling of extrasolar planets and our results show that planets at the outer edge of the habitable zones in compact multi-planet systems are vulnerable to rapid global glaciations. In the absence of external mechanisms such as orbital forcing or tidal heating, these planets could be trapped in permanent snowball states. In a second study, we examine the building blocks of pre-planetary materials using the outputs of an N-body planet formation paradigm and a Python-based volatile accretion model. For Earth-like protoplanets, our results suggest that volatile elemental ratios (e.g., C/N, C/H) are driven by the complex interplay between delivery, atmospheric ablation, and mantle degassing.