Search

General circulation models (GCMs) provide context for interpreting multiwavelength, multiphase data of the atmospheres of tidally locked exoplanets. In the current study, the non-hydrostatic THOR GCM is coupled with the HELIOS radiative transfer solver for the first time, supported by an equilibrium chemistry solver (FastChem), opacity calculator (HELIOS-K), and Mie scattering code (LX-MIE). To accurately treat the scattering of radiation by medium-sized to large aerosols/condensates, improved two-stream radiative transfer is implemented within a GCM for the first time. Multiple scattering is implemented using a Thomas algorithm formulation of the two-stream flux solutions, which decreases the computational time by about 2 orders of magnitude compared to the iterative method used in past versions of HELIOS. As a case study, we present four GCMs of the hot Jupiter WASP-43b, where we compare the temperature, velocity, entropy, and streamfunction, as well as the synthetic spectra and phase curves, of runs using regular versus improved two-stream radiative transfer and isothermal versus non-isothermal layers. While the global climate is qualitatively robust, the synthetic spectra and phase curves are sensitive to these details. A THOR + HELIOS WASP-43b GCM (horizontal resolution of about 4 deg on the sphere and with 40 radial points) with multiwavelength radiative transfer (30 k-table bins) running for 3000 Earth days (864 000 time-steps) takes about 19–26 d to complete depending on the type of GPU. [ABSTRACT FROM AUTHOR]

Accurate line lists are crucial for correctly modelling a variety of astrophysical phenomena, including stellar photospheres and the atmospheres of extrasolar planets. This paper presents a new line database Toto for the main isotopologues of titanium oxide (TiO): |$^{46}\text{Ti}^{16}\text{O}$|⁠ , |$^{47}\text{Ti}^{16}\text{O}$|⁠ , |$^{48}\text{Ti}^{16}\text{O}$|⁠ , |$^{49}\text{Ti}^{16}\text{O}$|⁠ , and |$^{50}\text{Ti}^{16}\text{O}$|⁠. The |$^{48}\text{Ti}^{16}\text{O}$| line list contains transitions with wave-numbers up to 30 000 cm−1, i.e. longwards of 0.33 μm. The Toto line list includes all dipole-allowed transitions between 13 low-lying electronic states (X 3Δ, a1Δ, d 1Σ+, E 3Π, A 3Φ, B3Π, C 3Δ, b 1Π, c1Φ, f1Δ, e 1Σ+). Ab initio potential energy curves (PECs) are computed at the icMRCI level and combined with spin–orbit and other coupling curves. These PECs and couplings are iteratively refined to match known empirical energy levels. Accurate line intensities are generated using ab initio dipole moment curves. The Toto line lists are appropriate for temperatures below 5000 K and contain 30 million transitions for |$^{48}\text{Ti}^{16}\text{O}$|⁠ ; it is made available in electronic form via the CDS data centre and via www.exomol.com. Tests of the line lists show greatly improved agreement with observed spectra for objects such as M-dwarfs GJ876 and GL581. [ABSTRACT FROM AUTHOR]

Current observational data of exoplanets are providing increasing detail of their 3D atmospheric structures. As characterization efforts expand in scope, the need to develop consistent 3D radiative-transfer methods becomes more pertinent as the complex atmospheric properties of exoplanets are required to be modelled together consistently. We aim to compare the transmission and emission spectra results of a 3D Monte Carlo radiative transfer (MCRT) model to contemporary radiative-transfer suites. We perform several benchmarking tests of an MCRT code, Cloudy Monte Carlo Radiative Transfer (cmcrt), to transmission and emission spectra model output. We add flexibility to the model through the use of k -distribution tables as input opacities. We present a hybrid MCRT and ray tracing methodology for the calculation of transmission spectra with a multiple scattering component. cmcrt compares well to the transmission spectra benchmarks at the 10s of ppm level. Emission spectra benchmarks are consistent to within 10 per cent of the 1D models. We suggest that differences in the benchmark results are likely caused by geometric effects between plane-parallel and spherical models. In a practical application, we post-process a cloudy 3D HD 189733b GCM model and compare to available observational data. Our results suggest the core methodology and algorithms of cmcrt produce consistent results to contemporary radiative transfer suites. 3D MCRT methods are highly suitable for detailed post-processing of cloudy and non-cloudy 1D and 3D exoplanet atmosphere simulations in instances where atmospheric inhomogeneities, significant limb effects/geometry or multiple scattering components are important considerations. [ABSTRACT FROM AUTHOR]

Terrestrial planets are thought to be the result of a vast number of gravitational interactions and collisions between smaller bodies. We use numerical simulations to show that practically identical initial conditions result in a wide array of final planetary configurations. This is a result of the chaotic evolution of trajectories which are highly sensitive to minuscule displacements. We determine that differences between systems evolved from virtually identical initial conditions can be larger than the differences between systems evolved from very different initial conditions. This implies that individual simulations lack predictive power. For example, there is not a reproducible mapping between the initial and final surface density profiles. However, some key global properties can still be extracted if the statistical spread across many simulations is considered. Based on these spreads, we explore the collisional growth and orbital properties of terrestrial planets, which assemble from different initial conditions (we vary the initial planetesimal distribution, planetesimal masses, and giant planet orbits.). Confirming past work, we find that the resulting planetary systems are sculpted by sweeping secular resonances. Configurations with giant planets on eccentric orbits produce fewer and more massive terrestrial planets on tighter orbits than those with giants on circular orbits. This is further enhanced if the initial mass distribution is biased to the inner regions. In all cases, the outer edge of the system is set by the final location of the ν6 resonance and we find that the mass distribution peaks at the ν5 resonance. Using existing observations, we find that extrasolar systems follow similar trends. Although differences between our numerical modelling and exoplanetary systems remain, we suggest that CoRoT-7, HD 20003 and HD 20781 may host undetected giant planets. [ABSTRACT FROM AUTHOR]