Flexible 3D models to explore the vast diversity of terrestrial planets and interpret observational data are still in their early stages. In this work, we present OASIS: a novel and flexible 3D virtual planet laboratory. With OASIS we envision a platform that couples self-consistently seven individual modules representing the main physical and chemical processes that shape planetary environments. Additionally, OASIS is capable of producing simulated spectra from different instruments and observational techniques. In this work we focus on the benchmark test of coupling four of the physical modules: fluid dynamics, radiation, turbulence and surface/soil. To test the OASIS platform, we produced 3D simulations of the Venus climate and its atmospheric circulation and study how the modeled atmosphere changes with various cloud covers, atmospheric heat capacity, and surface friction. 3D simulations of Venus are challenging because they require long integration times with a computationally expensive radiative transfer code. By comparing OASIS results with observational data, we verify that the new model is able to successfully simulate Venus. With simulated spectra produced directly from the 3D simulations, we explore the capabilities of future missions, like LUVOIR, to observe Venus analogs located at a distance of 10 pc. With OASIS, we have taken the first steps to build a sophisticated and very flexible platform capable of studying the environment of terrestrial planets, which will be an essential tool to characterize observed terrestrial planets and plan future observations.

THOR is the first open-source general circulation model (GCM) developed from scratch to study the atmospheres and climates of exoplanets, free from Earth- or Solar System-centric tunings. It solves the general non-hydrostatic Euler equations (instead of the primitive equations) on a sphere using the icosahedral grid. In the current study, we report major upgrades to THOR, building upon the work of Mendon\c{c}a et al. (2016). First, while the Horizontally Explicit Vertically Implicit (HEVI) integration scheme is the same as that described in Mendon\c{c}a et al. (2016), we provide a clearer description of the scheme and improved its implementation in the code. The differences in implementation between the hydrostatic shallow (HSS), quasi-hydrostatic deep (QHD) and non-hydrostatic deep (NHD) treatments are fully detailed. Second, standard physics modules are added: two-stream, double-gray radiative transfer and dry convective adjustment. Third, THOR is tested on additional benchmarks: tidally-locked Earth, deep hot Jupiter, acoustic wave, and gravity wave. Fourth, we report that differences between the hydrostatic and non-hydrostatic simulations are negligible in the Earth case, but pronounced in the hot Jupiter case. Finally, the effects of the so-called "sponge layer", a form of drag implemented in most GCMs to provide numerical stability, are examined. Overall, these upgrades have improved the flexibility, user-friendliness, and stability of THOR.

The atmospheric circulation in the upper atmosphere of hot Jupiter planets is strongly influenced by the incoming stellar radiation. In this work we explore the results from a 3D atmospheric model and revisit the main processes driving the circulation in hot Jupiter planets. We use the angular momentum transport as a diagnostic and carry out a Fourier analysis to identify the atmospheric waves involved. We find that the coupling between the angular momentum transported horizontally by the semi-diurnal tide and the mean circulation is the mechanism responsible for producing the strong jet at low latitudes. Our simulations indicate the possible formation of atmospheric indirect cells at low latitudes. The formation of these cells is induced by the presence of the semi-diurnal tide that is driven by the stellar irradiation. The tropical circulation has an important impact transporting heat and momentum from the upper towards the lower atmosphere. One of the consequences of this heat and momentum transport is a global increase of the temperature. We show that the initial conditions do not affect the output of the reference simulation. However, when the period of rotation of the planet was increased ($P_{rot} > 5$ Earth days), vertical transport by stationary waves became stronger, transient waves became non-negligible, and Coriolis influence less dominant, which allowed a steady state with a strong retrograde jet to be stable. We found that at least two statically steady state solutions exist for the same planet parameters.

Dust is ubiquitous in the Universe and its influence on the observed Electromagnetic (EM) radiation needs to be correctly addressed. In recent years it became clear that scattering of EM radiation from interstellar dust grains could change the local properties of the observed Cosmic Microwave Background (CMB) radiation. Here we consider the relevant processes of emission and scattering of EM radiation from spinning dust particles, and discuss their possible influence on the CMB. In particular, we show that scattered radiation can establish a correlation between different spectral components of galactic dipolar emission. This could explain the observed correlation between the CMB and the 100-micron thermal emission form interstellar dust. Another important property of CMB is related with its polarisation anisotropies, and the observation of a cosmological B-mode. We show that scattering of CMB radiation from dust grains in the presence of a static magnetic field could indeed create a B-mode spectral component, which is very similar to that due to primordial gravitational waves. This can be described by a kind of Cotton-Mutton effect on the CMB radiation.

The multi-messenger observation of GW170817 enabled the first historic measure of the Hubble-Lema\^itre constant via a standard siren, so called in analogy to standard candles that enabled the measure of the luminosity distance versus redshift relationship at small redshift. In the next decades third generation observatories are expected to detect hundreds to thousand gravitational wave events from compact binary coalescences with potentially a joint electromagnetic counterpart. In the present work we show how future standard siren detections can be used within the framework of Bayesian model selection to discriminate between cosmological models differing by the parametrization of the late time acceleration. In particular we show that the standard LambdaCDM model can be confirmed or ruled out in case the uncertainty in the gravitational determination of the luminosity distance is reduced with respect to current expectations, e.g. by combining detections from multiple detectors.

We present new methodological features and physical ingredients included in the 1D radiative transfer code HELIOS, improving the hemispheric two-stream formalism. We conduct a thorough intercomparison survey with several established forward models, including COOLTLUSTY, PHOENIX, and find satisfactory consistency with their results. Then, we explore the impact of (i) different groups of opacity sources, (ii) a stellar path length adjustment, and (iii) a scattering correction on self-consistently calculated atmospheric temperatures and planetary emission spectra. First, we observe that temperature-pressure (T-P) profiles are very sensitive to the opacities included, with metal oxides, hydrides, the alkali atoms (and ionized hydrogen) playing an important role for the absorption of shortwave radiation (in very hot surroundings). Moreover, if these species are sufficiently abundant, they are likely to induce non-monotonic T-P profiles. Second, without the stellar path length adjustment, the incoming stellar flux is significantly underestimated for zenith angles above 80{\deg}, which somewhat affects the upper atmospheric temperatures and the planetary emission. Third, the scattering correction improves the accuracy of the computation of the reflected stellar light by ~10%. We use HELIOS to calculate a grid of cloud-free atmospheres in radiative-convective equilibrium for self-luminous planets for a range of effective temperatures, surface gravities, metallicities, and C/O ratios, to be used by planetary evolution studies. Furthermore, we calculate dayside temperatures and secondary eclipse spectra for a sample of exoplanets for varying chemistry and heat redistribution. These results may be used to make predictions on the feasibility of atmospheric characterizations with future observations.

Spectral features in the observed spectra of exoplanets depend on the composition of their atmospheres. A good knowledge of the main atmospheric processes that drive the chemical distribution is therefore essential to interpret exoplanetary spectra. An atmosphere reaches chemical equilibrium if the rates of the forward and backward chemical reactions converge to the same value. However, there are atmospheric processes, such as atmospheric transport, that destabilize this equilibrium. In this work we study the changes in composition driven by a 3D wind field in WASP-43b using our Global Circulation Model, THOR. Our model uses validated temperature- and pressure-dependent chemical timescales that allow us to explore the disequilibrium chemistry of CO, CO$_2$, H$_2$O and CH$_4$. In WASP-43b the formation of the equatorial jet has an important impact in the chemical distribution of the different species across the atmosphere. At low latitudes the chemistry is longitudinally quenched, except for CO$_2$ at solar abundances. The polar vortexes have a distinct chemical distribution since these are regions with lower temperature and atmospheric mixing. Vertical and latitudinal mixing have a secondary impact in the chemical transport. We determine graphically the effect of disequilibrium on observed emission spectra. Our results do not show any significant differences in the emission spectra between the equilibrium and disequilibrium solutions for C/O = 0.5. However, if C/O is increased to 2.0, differences in the spectra due to the disequilibrium chemistry of CH$_4$ become non-negligible. In some spectral ranges the emission spectra can have more than 15$\%$ departures from the equilibrium solution.

Motivated by the work of Cooper & Showman, we revisit the chemical relaxation method, which seeks to enhance the computational efficiency of chemical-kinetics calculations by replacing the chemical network with a handful of independent source/sink terms. Chemical relaxation solves the evolution of the system and can treat disequilibrium chemistry, as the source/sink terms are driven towards chemical equilibrium on a prescribed chemical timescale, but it has surprisingly never been validated. First, we generalize the treatment by forgoing the use of a single chemical timescale, instead developing a pathway analysis tool that allows us to identify the rate-limiting reaction as a function of temperature and pressure. For the interconversion between methane and carbon monoxide and between ammonia, and molecular nitrogen, we identify the key rate-limiting reactions for conditions relevant to currently characterizable exo-atmospheres (500-3000 K, 0.1 mbar to 1 kbar). Second, we extend chemical relaxation to include carbon dioxide and water. Third, we examine the role of metallicity and carbon-to-oxygen ratio in chemical relaxation. Fourth, we apply our pathway analysis tool to diagnose the differences between our chemical network and that of Moses and Venot. Finally, we validate the chemical relaxation method against full chemical kinetics calculations in one dimension. For WASP-18b-, HD 189733b- and GJ 1214-b-like atmospheres, we show that chemical relaxation is mostly accurate to within an order of magnitude, a factor of 2 and $\sim 10\%$, respectively. The level of accuracy attained allows for the chemical relaxation method to be included in three-dimensional general circulation models.

Recently acquired Hubble and Spitzer phase curves of the short-period hot Jupiter WASP-43b make it an ideal target for confronting theory with data. On the observational front, we re-analyze the 3.6 and 4.5 $\mu$m Spitzer phase curves and demonstrate that our improved analysis better removes residual red noise due to intra-pixel sensitivity, which leads to greater fluxes emanating from the nightside of WASP-43b, thus reducing the tension between theory and data. On the theoretical front, we construct cloudfree and cloudy atmospheres of WASP-43b using our Global Circulation Model (GCM), THOR, which solves the non-hydrostatic Euler equations (compared to GCMs that typically solve the hydrostatic primitive equations). The cloudfree atmosphere produces a reasonable fit to the dayside emission spectrum. The multi-phase emission spectra constrain the cloud deck to be confined to the nightside and have a finite cloud-top pressure. The multi-wavelength phase curves are naturally consistent with our cloudy atmospheres, except for the 4.5 $\mu$m phase curve, which requires the presence of enhanced carbon dioxide in the atmosphere of WASP-43b. Multi-phase emission spectra at higher spectral resolution, as may be obtained using the James Webb Space Telescope, and a reflected-light phase curve at visible wavelengths would further constrain the properties of clouds in WASP-43b.

We present an open-source retrieval code named HELIOS-Retrieval (hereafter HELIOS-R), designed to obtain chemical abundances and temperature-pressure profiles from inverting the measured spectra of exoplanetary atmospheres. In the current implementation, we use an exact solution of the radiative transfer equation, in the pure absorption limit, in our forward model, which allows us to analytically integrate over all of the outgoing rays (instead of performing Gaussian quadrature). Two chemistry models are considered: unconstrained chemistry (where the mixing ratios are treated as free parameters) and equilibrium chemistry (enforced via analytical formulae, where only the elemental abundances are free parameters). The nested sampling algorithm allows us to formally implement Occam's Razor based on a comparison of the Bayesian evidence between models. We perform a retrieval analysis on the measured spectra of the HR 8799b, c, d and e directly imaged exoplanets. Chemical equilibrium is disfavored by the Bayesian evidence for HR 8799b, c and d. We find supersolar C/O, C/H and O/H values for the outer HR 8799b and c exoplanets, while the inner HR 8799d and e exoplanets have substellar C/O, substellar C/H and superstellar O/H values. If these retrieved properties are representative of the bulk compositions of the exoplanets, then they are inconsistent with formation via gravitational instability (without late-time accretion) and consistent with a core accretion scenario in which late-time accretion of ices occurred differently for the inner and outer exoplanets. For HR 8799e, we find that spectroscopy in the K band is crucial for constraining C/O and C/H. HELIOS-R is publicly available as part of the Exoclimes Simulation Platform (ESP; www.exoclime.org).

The work investigates the coupling between ion-acoustic waves and neutrino flavor oscillations in a non-relativistic electron-ion plasma under the influence of a mixed neutrino beam. Neutrino oscillations are mediated by the flavor polarization vector dynamics in a material medium. The linear dispersion relation around homogeneous static equilibria is developed. When resonant with the ion-acoustic mode, the neutrino flavor oscillations can transfer energy to the plasma exciting a new fast unstable mode in extreme astrophysical scenarios. The growth rate and the unstable wavelengths are determined in typical type II supernovae parameters. The predictions can be useful for a new indirect probe on neutrino oscillations in nature.

We present the open-source radiative transfer code named HELIOS, which is constructed for studying exoplanetary atmospheres. In its initial version, the model atmospheres of HELIOS are one-dimensional and plane-parallel, and the equation of radiative transfer is solved in the two-stream approximation with non-isotropic scattering. A small set of the main infrared absorbers is employed, computed with the opacity calculator HELIOS-K and combined using a correlated-$k$ approximation. The molecular abundances originate from validated analytical formulae for equilibrium chemistry. We compare HELIOS with the work of Miller-Ricci & Fortney using a model of GJ 1214b, and perform several tests, where we find: model atmospheres with single-temperature layers struggle to converge to radiative equilibrium; $k$-distribution tables constructed with $\gtrsim 0.01$ cm$^{-1}$ resolution in the opacity function ($ \lesssim 10^3$ points per wavenumber bin) may result in errors $\gtrsim 1$-10 % in the synthetic spectra; and a diffusivity factor of 2 approximates well the exact radiative transfer solution in the limit of pure absorption. We construct "null-hypothesis" models (chemical equilibrium, radiative equilibrium and solar element abundances) for 6 hot Jupiters. We find that the dayside emission spectra of HD 189733b and WASP-43b are consistent with the null hypothesis, while it consistently under-predicts the observed fluxes of WASP-8b, WASP-12b, WASP-14b and WASP-33b. We demonstrate that our results are somewhat insensitive to the choice of stellar models (blackbody, Kurucz or PHOENIX) and metallicity, but are strongly affected by higher carbon-to-oxygen ratios. The code is publicly available as part of the Exoclimes Simulation Platform (ESP; exoclime.net).

The atmospheric circulation in Venus is well known to exhibit strong super-rotation. However, the atmospheric mechanisms responsible for the formation of this super-rotation are still not fully understood. In this work, we developed a new Venus general circulation model to study the most likely mechanisms driving the atmosphere to the current observed circulation. Our model includes a new radiative transfer, convection and suitably adapted boundary layer schemes and a dynamical core that takes into account the dependence of the heat capacity at constant pressure with temperature. The new Venus model is able to simulate a super-rotation phenomenon in the cloud region quantitatively similar to the one observed. The mechanisms maintaining the strong winds in the cloud region were found in the model results to be a combination of zonal mean circulation, thermal tides and transient waves. In this process, the semi-diurnal tide excited in the upper clouds has a key contribution in transporting axial angular momentum mainly from the upper atmosphere towards the cloud region. The magnitude of the super-rotation in the cloud region is sensitive to various radiative parameters such as the amount of solar radiative energy absorbed by the surface, which controls the static stability near the surface. In this work, we also discuss the main difficulties in representing the flow below the cloud base in Venus atmospheric models. Our new radiative scheme is more suitable for 3D Venus climate models than those used in previous work due to its easy adaptability to different atmospheric conditions. This flexibility of the model was crucial to explore the uncertainties in the lower atmospheric conditions and may also be used in the future to explore, for example, dynamical-radiative-microphysical feedbacks.

We have designed and developed, from scratch, a global circulation model named THOR that solves the three-dimensional non-hydrostatic Euler equations. Our general approach lifts the commonly used assumptions of a shallow atmosphere and hydrostatic equilibrium. We solve the "pole problem" (where converging meridians on a sphere lead to increasingly smaller time steps near the poles) by implementing an icosahedral grid. Irregularities in the grid, which lead to grid imprinting, are smoothed using the "spring dynamics" technique. We validate our implementation of spring dynamics by examining calculations of the divergence and gradient of test functions. To prevent the computational time step from being bottlenecked by having to resolve sound waves, we implement a split-explicit method together with a horizontally explicit and vertically implicit integration. We validate our global circulation model by reproducing the Earth and also the hot Jupiter-like benchmark tests. THOR was designed to run on Graphics Processing Units (GPUs), which allows for physics modules (radiative transfer, clouds, chemistry) to be added in the future, and is part of the open-source Exoclimes Simulation Platform (ESP; www.exoclime.org).

Turbulent radiation flow is ubiquitous in many physical systems where light-matter interaction becomes relevant. Photon bubbling, in particular, has been identified as the main source of turbulent radiation transport in many astrophysical objects, such as stars and accretion disks. This mechanism takes place when radiation trapping in optically dense media becomes unstable, leading to the energy dissipation from the larger to the smaller bubbles. Here, we report on the observation of photon bubble turbulence in cold atomic gases in the presence of multiple scattering of light. The instability is theoretically explained by a fluid description for the atom density coupled to a diffusive transport equation for the photons, which is known to be accurate in the multiple scattering regime investigated here. We determine the power spectrum of the atom density fluctuations, which displays an unusual $\sim k^{-4}$ scaling, and entails a complex underlying turbulent dynamics resulting from the formation of dynamical bubble-like structures. We derive a power spectrum from the theoretical photon bubble model which, to a high level of accuracy, explains the observations. The experimental results reported here, along with the theoretical model we developed may shed light on the analogue photon bubble instabilities in astrophysical scenarios.

We present a comprehensive analytical study of radiative transfer using the method of moments and include the effects of non-isotropic scattering in the coherent limit. Within this unified formalism, we derive the governing equations and solutions describing two-stream radiative transfer (which approximates the passage of radiation as a pair of outgoing and incoming fluxes), flux-limited diffusion (which describes radiative transfer in the deep interior) and solutions for the temperature-pressure profiles. Generally, the problem is mathematically under-determined unless a set of closures (Eddington coefficients) is specified. We demonstrate that the hemispheric (or hemi-isotropic) closure naturally derives from the radiative transfer equation if energy conservation is obeyed, while the Eddington closure produces spurious enhancements of both reflected light and thermal emission. We concoct recipes for implementing two-stream radiative transfer in stand-alone numerical calculations and general circulation models. We use our two-stream solutions to construct toy models of the runaway greenhouse effect. We present a new solution for temperature-pressure profiles with a non-constant optical opacity and elucidate the effects of non-isotropic scattering in the optical and infrared. We derive generalized expressions for the spherical and Bond albedos and the photon deposition depth. We demonstrate that the value of the optical depth corresponding to the photosphere is not always 2/3 (Milne's solution) and depends on a combination of stellar irradiation, internal heat and the properties of scattering both in optical and infrared. Finally, we derive generalized expressions for the total, net, outgoing and incoming fluxes in the convective regime.

We analyse the Anderson-Higgs mechanism of photon mass acquisition in a plasma and study the contribution to the mass from the orbital angular momentum acquired by a beam of photons when it crosses a spatially structured charge distribution. To this end we apply Proca-Maxwell equations in a static plasma with a particular spatial distribution of free charges, notably a plasma vortex, that is able to impose orbital angular momentum (OAM) onto light. In addition to the mass acquisition of the conventional Anderson-Higgs mechanism, we find that the photon acquires an additional mass from the OAM and that this mass reduces the Proca photon mass.

Electromagnetic waves with an azimuthal phase shift are known to have a well defined orbital angular momentum. Different methods that allow for the detection of the angular momentum are proposed. For some, we discuss the required experimental setup and explore the range of applicability.

It is shown that an electron-neutrino beam, propagating in a background plasma, can be decomposed into orbital momentum (OAM) states, similar to the OAM photon states. Coupling between different OAM neutrino states, in the presence of a plasma vortex, is considered. We show that plasma vorticity can be transfered to the neutrino beam, which is relevant to the understanding of the neutrino sources in astrophysics. Observation of neutrino OAM states could eventually become possible.