We present a new set of solar metallicity atmosphere and evolutionary models for very cool brown dwarfs and self-luminous giant exoplanets, which we term ATMO 2020. Atmosphere models are generated with our state-of-the-art 1D radiative-convective equilibrium code ATMO, and are used as surface boundary conditions to calculate the interior structure and evolution of $0.001-0.075\,\mathrm{M_{\odot}}$ objects. Our models include several key improvements to the input physics used in previous models available in the literature. Most notably, the use of a new H-He equation of state including ab initio quantum molecular dynamics calculations has raised the mass by $\sim1-2\%$ at the stellar-substellar boundary and has altered the cooling tracks around the hydrogen and deuterium burning minimum masses. A second key improvement concerns updated molecular opacities in our atmosphere model ATMO, which now contains significantly more line transitions required to accurately capture the opacity in these hot atmospheres. This leads to warmer atmospheric temperature structures, further changing the cooling curves and predicted emission spectra of substellar objects. We present significant improvement for the treatment of the collisionally broadened potassium resonance doublet, and highlight the importance of these lines in shaping the red-optical and near-infrared spectrum of brown dwarfs. We generate three different grids of model simulations, one using equilibrium chemistry and two using non-equilibrium chemistry due to vertical mixing, all three computed self-consistently with the pressure-temperature structure of the atmosphere. We show the impact of vertical mixing on emission spectra and in colour-magnitude diagrams, highlighting how the $3.5-5.5\,\mathrm{\mu m}$ flux window can be used to calibrate vertical mixing in cool T-Y spectral type objects.

Atmospheric compositions can provide powerful diagnostics of formation and migration histories of planetary systems. We investigate constraints on atmospheric abundances of H$_2$O, Na, and K, in a sample of transiting exoplanets using latest transmission spectra and new H$_2$ broadened opacities of Na and K. Our sample of 19 exoplanets spans from cool mini-Neptunes to hot Jupiters, with equilibrium temperatures between $\sim$300 and 2700 K. Using homogeneous Bayesian retrievals we report atmospheric abundances of Na, K, and H$_2$O, and their detection significances, confirming 6 planets with strong Na detections, 6 with K, and 14 with H$_2$O. We find a mass-metallicity trend of increasing H$_2$O abundances with decreasing mass, spanning generally substellar values for gas giants and stellar/superstellar for Neptunes and mini-Neptunes. However, the overall trend in H$_2$O abundances, from mini-Neptunes to hot Jupiters, is significantly lower than the mass-metallicity relation for carbon in the solar system giant planets and similar predictions for exoplanets. On the other hand, the Na and K abundances for the gas giants are stellar or superstellar, consistent with each other, and generally consistent with the solar system metallicity trend. The H$_2$O abundances in hot gas giants are likely due to low oxygen abundances relative to other elements rather than low overall metallicities, and provide new constraints on their formation mechanisms. The differing trends in the abundances of species argue against the use of chemical equilibrium models with metallicity as one free parameter in atmospheric retrievals, as different elements can be differently enhanced.

The opacity of alkali atoms, most importantly of Na and K, plays a crucial role in the atmospheres of brown dwarfs and exoplanets. We present a comprehensive study of NaH2 collisional profiles at temperatures from 500 to 3000 K, the temperatures prevailing in the atmosphere of brown dwarfs and Jupiter-mass planets.The relevant H2 perturber densities reach several 10^19 cm^-3 in hot Teff > 1500 K Jupiter-mass planets and can exceed 10^20 cm^-3 for more massive or cooler objects. Accurate pressure-broadened profiles that are valid at high densities of H2 should be incorporated into spectral models. Unified profiles of sodium perturbed by molecular hydrogen were calculated in the semi-classical approach using up-to-date molecular data.New NaH2 collisional profiles and their effects on the synthetic spectra of brown dwarfs and hot Jupiters computed with petitCODE are presented.

T-type dwarfs present a broad and shallow absorption feature centred around 6950 A in the blue wing of the K doublet at 0.77 micron which resembles in depth and shape the satellite absorption predicted by detailed collisional broadening profiles. In our previous work, the predicted line satellite position was however somewhat too blue compared to the observed feature. In this paper we investigate whether new calculations of the energy surfaces of the potentials in the K-H_2 system, including spin-orbit coupling, result in a closer coincidence of the satellite with the observed position. We also investigate the extent to which CaH absorption bands contribute to the feature. We present model atmospheres and synthetic spectra, including gravitational settling for an improved description of depth-dependent abundances of refractory elements, and based on new K-H_2 line profiles using improved interaction potentials. By comparison with a high signal-to-noise optical spectrum of the T1 dwarf epsilon Indi Ba, we find that these new models do reproduce the observed feature, while CaH does not contribute for the atmospheric parameters considered. We also find that CaH is settled out so deep into the atmosphere that even turbulent vertical mixing would appear insufficient to bring significant amounts of CaH to the photosphere in dwarfs later than ~L5. We conclude that previous identification of the feature at this location in T and late L dwarf spectra with CaH was erroneous, as expected on physical grounds: calcium condenses onto grains in early L dwarfs and thus should have settled out of the photosphere in cooler brown dwarfs. This finding revokes one observational verification for the cloud-clearing theory: a gradual clearing of the cloud cover in early T dwarfs.

The competition between the formation and destruction of coronene clusters under interstellar conditions is investigated theoretically. The unimolecular nucleation of neutral clusters is simulated with an atomic model combining an explicit classical force field and a quantum tight-binding approach. Evaporation rates are calculated in the framework of the phase space theory and are inserted in an infrared emission model and compared with the growth rate constants. It is found that, in interstellar conditions, most collisions lead to cluster growth. The time evolution of small clusters (containing up to 312 carbon atoms) was specifically investigated under the physical conditions of the northern photodissociation region of NGC 7023. These clusters are found to be thermally photoevaporated much faster than they are reformed, thus providing an interpretation for the lowest limit of the interstellar cluster size distribution inferred from observations. The effects of ionizing the clusters and density heterogeneities are also considered. Based on our results, the possibility that PAH clusters could be formed in PDRs is critically discussed.