This paper presents a detailed hydrostatic model of the upper atmosphere of HD 189733b, with the goal of constraining its temperature, particle densities, and radiation field over the pressure range 10 −4 −10 µbar, where the observed Hα transmission spectrum is produced. The atomic hydrogen level population is computed including both collisional and radiative transition rates. The Lyα resonant scattering is computed using a Monte-Carlo simulation. The model transmission spectra are in broad agreement with the data. Excitation of the H(2 ) population is mainly by Lyα radiative excitation due to the large Lyα intensity. The density of H(2 ) is nearly flat over two decades in pressure, and is optically thick to Hα. Additional models computed for a range of the stellar Lyman continuum (LyC) flux suggest that the variability in Hα transit depth may be due to the variability in the stellar LyC. Since metal lines provide the dominant cooling of this part of the atmosphere, the atmosphere structure is sensitive to the density of species such as Mg and Na which may themselves be constrained by observations. Since the Hα and Na D lines have comparable absorption depths, we argue that the center of the Na D lines are also formed in the atomic layer where the Hα line is formed.

We investigate giant molecular cloud collisions and their ability to induce gravitational instability and thus star formation. This mechanism may be a major driver of star formation activity in galactic disks. We carry out a series of 3D, magnetohydrodynamics (MHD), adaptive mesh refinement simulations to study how cloud collisions trigger formation of dense filaments and clumps. Heating and cooling functions are implemented based on photodissociation region models that span the atomic-to-molecular transition and can return detailed diagnostic information. The clouds are initialized with supersonic turbulence and a range of magnetic field strengths and orientations. Collisions at various velocities and impact parameters are investigated. Comparing and contrasting colliding and non-colliding cases, we characterize morphologies of dense gas, magnetic field structure, cloud kinematic signatures, and cloud dynamics. We present key observational diagnostics of cloud collisions, especially: relative orientations between magnetic fields and density structures, like filaments; 13 CO( J=2-1), 13 CO( J=3-2), and 12 CO( J=8-7) integrated intensity maps and spectra; and cloud virial parameters. We compare these results to observed Galactic clouds.

We study giant molecular cloud (GMC) collisions and their ability to trigger star cluster formation. We further develop our three-dimensional magnetized, turbulent, colliding GMC simulations by implementing star formation subgrid models. Two such models are explored: (1) "Density-Regulated," i.e., fixed efficiency per free-fall time above a set density threshold and (2) "Magnetically Regulated," i.e., fixed efficiency per free-fall time in regions that are magnetically supercritical. Variations of parameters associated with these models are also explored. In the non-colliding simulations, the overall level of star formation is sensitive to model parameter choices that relate to effective density thresholds. In the GMC collision simulations, the final star formation rates and efficiencies are relatively independent of these parameters. Between the non-colliding and colliding cases, we compare the morphologies of the resulting star clusters, properties of star-forming gas, time evolution of the star formation rate (SFR), spatial clustering of the stars, and resulting kinematics of the stars in comparison to the natal gas. We find that typical collisions, by creating larger amounts of dense gas, trigger earlier and enhanced star formation, resulting in 10 times higher SFRs and efficiencies. The star clusters formed from GMC collisions show greater spatial substructure and more disturbed kinematics.

Gas giant exoplanets orbiting at close distances to the parent star are subjected to large radiation and stellar wind fluxes. In this paper, hydrodynamic simulations of the planetary upper atmosphere and its interaction with the stellar wind are carried out to understand the possible flow regimes and how they affect the Lyα transmission spectrum. Following Tremblin and Chiang, charge exchange reactions are included to explore the role of energetic atoms as compared to thermal particles. In order to understand the role of the tail as compared to the leading edge of the planetary gas, the simulations were carried out under axisymmetry, and photoionization and stellar wind electron impact ionization reactions were included to limit the extent of the neutrals away from the planet. By varying the planetary gas temperature, two regimes are found. At high temperature, a supersonic planetary wind is found, which is turned around by the stellar wind and forms a tail behind the planet. At lower temperatures, the planetary wind is shut off when the stellar wind penetrates inside where the sonic point would have been. In this regime mass is lost by viscous interaction at the boundary between planetary and stellar wind gases. Absorption by cold hydrogen atoms is large near the planetary surface, and decreases away from the planet as expected. The hot hydrogen absorption is in an annulus and typically dominated by the tail, at large impact parameter, rather than by the thin leading edge of the mixing layer near the substellar point.

It has been suggested that hot stellar wind gas in a bow shock around an exoplanet is sufficiently opaque to absorb stellar photons and give rise to an observable transit depth at optical and UV wavelengths. In the first part of this paper, we use the CLOUDY plasma simulation code to model the absorption from X-ray to radio wavelengths by 1-D slabs of gas in coronal equilibrium with varying densities (10 4 − 10 8 cm −3 ) and temperatures (2000 − 10 6 K) illuminated by a solar spectrum. For slabs at coronal temperatures (10 6 K) and densities even orders of magnitude larger than expected for the compressed stellar wind (10 4 − 10 5 cm −3 ), we find optical depths orders of magnitude too small (> 3 × 10 −7 ) to explain the ∼ 3% UV transit depths seen with Hubble. Using this result and our modeling of slabs with lower temperatures (2000 − 10 4 K), the conclusion is that the UV transits of WASP-12b and HD 189733b are likely due to atoms originating in the planet, as the stellar wind is too highly ionized. A corollary of this result is that transport of neutral atoms from the denser planetary atmosphere outward must be a primary consideration when constructing physical models. In the second part of this paper, additional calculations using CLOUDY are carried out to model a slab of planetary gas in radiative and thermal equilibrium with the stellar radiation field. Promising sources of opacity from the X-ray to radio wavelengths are discussed, some of which are not yet observed.