Search NASA⌕ Search

SEARCH · Search NASA

Results for “Stellar feedback”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Nonlinear evolution of protostellar disks and light modulations in young stellar objects

An evolutionary model of dynamical processes in protostellar disks is described and illustrated with graphs of typical results. The effective transport mechanisms are discussed, including thermal convection, nonaxisymmetric gravitational instabilities in the outer regions of disks, and wave propagation. Consideration is then given to the stages of dynamical evolution, FU Ori outburst phenomena, unsteady accretion-disk flows, and nonlinear feedback as a mechanism to modulate mass transfer. The simulations show that mass redistribution is determined by angular-momentum transfer, which in turn is regulated by the effective viscosity generated by convectively driven turbulence. Significant mass transfer occurs as a result of mixing of infalling material with disk gas and is affected by the tidal torque associated with the growth of nonaxisymmetric disturbances in the outer disk. The time scale for disk evolution is found to be about 1 Myr.

Lin, D. N. C.↗
DOE OSTI2022

Star formation and AGN feedback in the local Universe: Combining LOFAR and MaNGA

The effect of active galactic nuclei (AGN) on their host galaxies – in particular their levels of star formation – remains one of the key outstanding questions of galaxy evolution. Successful cosmological models of galaxy evolution require a fraction of energy released by an AGN to be redistributed into the interstellar medium to reproduce the observed stellar mass and luminosity function and to prevent the formation of over-massive galaxies. Observations have confirmed that the radio-AGN population is energetically capable of heating and redistributing gas at all phases, however, direct evidence of AGN enhancing or quenching star formation remains rare. With modern, deep radio surveys and large integral field spectroscopy (IFS) surveys, we can detect fainter synchrotron emission from AGN jets and accurately probe the star-forming properties of galaxies, respectively. Here, in this paper, we combine data from the LOw Frequency ARray Two-meter Sky Survey (LoTSS) with data from one of the largest optical IFS surveys, Mapping Nearby Galaxies at Apache Point Observatory (MaNGA), to probe the star-forming properties of 307 local (z < 0.15) galaxies that host radio-detected AGN (RDAGN). We compare our results to a robust control sample of non-active galaxies that each match the stellar mass, redshift, visual morphology, and inclination of a RDAGN host. We find that RDAGN and control galaxies have broad star-formation rate (SFR) distributions, typically lie below the SFMS, and have negative stellar light-weighted age gradients. These results indicate that AGN selected based on their current activity are not responsible for suppressing their host galaxies’ star formation. Rather, our results support the maintenance mode role that radio AGN are expected to have in the local Universe.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Overview of Wendelstein 7-X high-performance operation

The Wendelstein 7-X (W7-X) stellarator has completed two consecutive experimental campaigns OP 2.2 (Sep.-Dec. 2024) and OP 2.3 (Feb.-May 2025) under a new operational strategy enabling more than one year of uninterrupted device availability. This approach, supported by exceptionally high subsystem reliability, allowed sustained high-efficiency plasma operations with up to 80–100 discharges per day across a broad range of magnetic configurations. Several key technical upgrades-most notably the first operation of a 1.5 MW class steady-state gyrotron, a new steady-state pellet injector, and advanced real-time feedback control systems significantly enhanced heating, fueling, and plasma control capabilities. Together, these improvements enabled major advances in long-pulse performance, high-β operation, and confinement optimization. Long-pulse discharges achieved 1.8 GJ of injected energy under fully detached divertor conditions, while reduced-field scenarios facilitated record volume-averaged β values approaching 3%. High-performance plasmas with centrally peaked density profiles, created via neutral beam injection (NBI) or sustained pellet fueling, demonstrated strongly reduced turbulent transport and stellarator-record fusion triple products. Complementary studies of power exhaust and divertor heat loads revealed the role of scrape-off-layer drift physics in shaping strike-line patterns under attached conditions. Together, the results from OP 2.2 and OP 2.3 significantly expand the operational space of W7-X and strengthen its role as a leading platform for steady-state stellarator research and reactor-relevant plasma scenarios.

Wendelstein 7-X↗
DOE OSTI2024

Ne viii in the Warm-hot Circumgalactic Medium of FIRE Simulations and in Observations

Abstract The properties of warm-hot gas around ∼ L * galaxies can be studied with absorption lines from highly ionized metals. We predict Ne viii column densities from cosmological zoom-in simulations of halos with masses in ∼10 12 and ∼10 13 M ☉ from the Feedback in Realistic Environments (FIRE) project. Ne viii traces the volume-filling, virial-temperature gas in ∼10 12 M ☉ halos. In ∼10 13 M ☉ halos the Ne viii gas is clumpier, and biased toward the cooler part of the warm-hot phase. We compare the simulations to observations from the COS Absorption Survey of Baryon Harbors (or CASBaH) and COS Ultraviolet Baryon Survey (or CUBS). We show that when inferring halo masses from stellar masses to compare simulated and observed halos, it is important to account for the scatter in the stellar-mass–halo-mass relation, especially at M ⋆ ≳ 10 10.5 M ☉ . Median Ne viii columns in the fiducial FIRE-2 model are about as high as observed upper limits allow, while the simulations analyzed do not reproduce the highest observed columns. This suggests that the median Ne viii profiles predicted by the simulations are consistent with observations, but that the simulations may underpredict the scatter. We find similar agreement with analytical models that assume a product of the halo gas fraction and metallicity (relative to solar) ∼0.1, indicating that observations are consistent with plausible circumgalactic medium temperatures, metallicities, and gas masses. Variants of the FIRE simulations with a modified supernova feedback model and/or active galactic nuclei feedback included (as well as some other cosmological simulations from the literature) more systematically underpredict Ne viii columns. The circumgalactic Ne viii observations therefore provide valuable constraints on simulations that otherwise predict realistic galaxy properties.

Wijers, Nastasha A. (ORCID:0000000163747185)↗
DOE OSTI2025

A New Framework for Interstellar Medium Emission Line Models: Connecting Multiscale Simulations across Cosmological Volumes

The James Webb Space Telescope (JWST) and Atacama Large Millimeter/submillimeter Array have detected emission lines from the ionized interstellar medium (ISM) in some of the first galaxies at z ≳ 6. These measurements present an opportunity to better understand galaxy assembly histories and may allow important tests of state-of-the-art galaxy formation simulations. It is challenging, however, to model these lines in their proper cosmological context. In order to meet this challenge, we introduce a novel subgrid line emission modeling framework. The framework uses the high-z zoom-in simulation suite from the Feedback in Realistic Environments (FIRE) collaboration. The line emission signals from H II regions within each simulated FIRE galaxy are modeled using the semianalytic HIIL INES code. A machine learning approach is then used to determine the conditional probability distribution for the line luminosity to stellar-mass ratio from the H II regions around each simulated stellar particle. This conditional probability distribution can then be applied to predict the line luminosities around stellar particles in lower-resolution, yet larger volume cosmological simulations. As an example, we apply this approach to the IllustrisTNG simulations at z = 6. The resulting predictions for the [O II ], [O III ], and Balmer line luminosities as a function of star formation rate agree well with current observations. Our predictions differ, however, from related works in the literature, which lack detailed subgrid ISM models. This highlights the importance of our multiscale simulation modeling framework. Finally, we provide forecasts for future line luminosity function measurements from the JWST and quantify the cosmic variance in such surveys.

(ISM:) H II regions↗
NASA NTRS2017

Constraints on Climate and Habitability for Earth-like Exoplanets Determined from a General Circulation Model

Conventional definitions of habitability require abundant liquid surface water to exist continuously over geologic timescales. Water in each of its thermodynamic phases interacts with solar and thermal radiation and is the cause for strong climatic feedbacks. Thus, assessments of the habitable zone require models to include a complete treatment of the hydrological cycle over geologic time. Here, we use the Community Atmosphere Model from the National Center for Atmospheric Research to study the evolution of climate for an Earth-like planet at constant CO2, under a wide range of stellar fluxes from F-, G-, and K-dwarf main sequence stars. Around each star we find four stable climate states defined by mutually exclusive global mean surface temperatures (Ts); snowball (Ts ≤ 235 K), waterbelt (235 K ≤ Ts ≤ 250 K), temperate (275 K ≤ Ts ≤ 315 K), and moist greenhouse (Ts ≥ 330 K). Each is separated by abrupt climatic transitions. Waterbelt, temperate, and cooler moist greenhouse climates can maintain open-ocean against both sea ice albedo and hydrogen escape processes respectively, and thus constitute habitable worlds. We consider the warmest possible habitable planet as having Ts ∼ 355 K, at which point diffusion limited water-loss could remove an Earth ocean in ∼1 Gyr. Without long timescale regulation of non-condensable greenhouse species at Earth-like temperatures and pressures, such as CO2, habitability can be maintained for an upper limit of ∼2.2, ∼2.4, and ∼4.7 Gyr around F-, G-, and K-dwarf stars respectively, due to main sequence brightening.

Eric T Wolf↗
NASA NTRS2017

Constraints on Climate and Habitability for Earth-like Exoplanets Determined from a General Circulation Model

Conventional definitions of habitability require abundant liquid surface water to exist continuously over geologic timescales. Water in each of its thermodynamic phases interacts with solar and thermal radiation and is the cause for strong climatic feedbacks. Thus, assessments of the habitable zone require models to include a complete treatment of the hydrological cycle over geologic time. Here, we use the Community Atmosphere Model from the National Center for Atmospheric Research to study the evolution of climate for an Earth-like planet at constant CO2, under a wide range of stellar fluxes from F-, G-, and K-dwarf main sequence stars. Around each star we find four stable climate states defined by mutually exclusive global mean surface temperatures (T(sub s)); snowball (T(sub s) ≼ 235 K), waterbelt (235 K ≼ T(sub s) ≼ 250 K), temperate (275 K ≼ T(sub s) ≼ 315 K), and moist greenhouse (T(sub s) ≽ 330 K). Each is separated by abrupt climatic transitions. Waterbelt, temperate, and cooler moist greenhouse climates can maintain open-ocean against both sea ice albedo and hydrogen escape processes respectively, and thus constitute habitable worlds. We consider the warmest possible habitable planet as having T(sub s) ∼ 355 K, at which point diffusion limited water-loss could remove an Earth ocean in ∼1 Gyr. Without long timescale regulation of non-condensable greenhouse species at Earth-like temperatures and pressures, such as CO2, habitability can be maintained for an upper limit of ∼2.2, ∼2.4, and ∼4.7 Gyr around F-, G-, and K-dwarf stars respectively, due to main sequence brightening.

Astrobiology↗
DOE OSTI2023

H i -rich but low star formation galaxies in MaNGA: physical properties and comparison to control samples

ABSTRACT Gas-rich galaxies are typically star forming. We make use of H i-MaNGA, a programme of H i follow-up for the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey of the Sloan Digital Sky Surveys, to construct a sample of unusual neutral hydrogen (H i, 21 cm)-rich galaxies that have low star formation rates (SFRs), using infrared colour from the Wide-field Infrared Survey Explorer as a proxy for specific SFR. Out of a set of 1575 MaNGA galaxies with H i-MaNGA detections, we find that 83 (5 per cent) meet our selection criteria to be H i rich with low SFR. We construct two stellar mass-matched control samples: H i-rich galaxies with typical SFR (high SF control) and H i-poor galaxies with low SFR (low H i control). We investigate the properties of each of these samples, comparing physical parameters such as ionization state maps, stellar and ionized gas velocity and dispersion, environment measures, metallicity, and morphology to search for the reasons why these unusual H i-rich galaxies are not forming stars. We find evidence for recent external accretion of gas in some galaxies (via high counter-rotating fractions), along with some evidence for active galactic nucleus (AGN) feedback (from a high central low-ionization emission-line region and/or red geyser fraction), and bar quenching (via an enhanced strong bar fraction). Some galaxies in the sample are consistent with simply having their H i in a high angular momentum, large-radius, low-density disc. We conclude that no single physical process can explain all H i-rich, low-SFR galaxies.

Astronomy & Astrophysics↗
DOE OSTI2026

The DREAMS Project: Disentangling the Impact of Halo-to-halo Variance and Baryonic Feedback on Milky Way Satellite Galaxies

We analyze the properties of satellite galaxies around 1024 Milky Way–mass hosts from the DREAMS Project, simulated within a ΛCDM cosmology. Utilizing the TNG galaxy-formation model, the DREAMS simulations incorporate both baryonic physics and cosmological uncertainties for a large sample of galaxies with diverse environments and formation histories. We investigate the relative impact of the physical uncertainty from the galaxy-formation model on predicted satellite properties using four metrics: the satellite stellar mass function, radial distribution, inner slope of the dark matter density profile, and stellar half-light radius. We compare these predictions to observations from the SAGA survey and the DREAMS N-body simulations and find that uncertainties from baryonic physics modeling are subdominant to the scatter arising from halo-to-halo variance. Where baryonic modeling does affect satellites, the supernova wind energy has the largest effect on the satellite properties that we investigate. Specifically, increased supernova wind energy suppresses the stellar mass of satellites and results in more extended stellar half-light radii. The adopted wind speed has only a minor impact, and other astrophysical and cosmological parameters show no measurable effect. Our findings highlight the robustness of satellite properties against uncertainties in baryonic physics modeling.

Rose, Jonah C. [Princeton University, NJ (United S↗
DOE OSTI2021

Comparing galaxy formation in the L-GALAXIES semi-analytical model and the IllustrisTNG simulations

ABSTRACT We perform a comparison, object by object and statistically, between the Munich semi-analytical model, L-GALAXIES, and the IllustrisTNG hydrodynamical simulations. By running L-GALAXIES on the IllustrisTNG dark matter-only merger trees, we identify the same galaxies in the two models. This allows us to compare the stellar mass, star formation rate, and gas content of galaxies, as well as the baryonic content of subhaloes and haloes in the two models. We find that both the stellar mass functions and the stellar masses of individual galaxies agree to better than ${\sim} 0.2\,$dex. On the other hand, specific star formation rates and gas contents can differ more substantially. At z = 0, the transition between low-mass star-forming galaxies and high-mass quenched galaxies occurs at a stellar mass scale ${\sim} 0.5\,$dex lower in IllustrisTNG than that in L-GALAXIES. IllustrisTNG also produces substantially more quenched galaxies at higher redshifts. Both models predict a halo baryon fraction close to the cosmic value for clusters, but IllustrisTNG predicts lower baryon fractions in group environments. These differences are primarily due to differences in modelling feedback from stars and supermassive black holes. The gas content and star formation rates of galaxies in and around clusters and groups differ substantially, with IllustrisTNG satellites less star forming and less gas rich. We show that environmental processes such as ram-pressure stripping are stronger and operate to larger distances and for a broader host mass range in IllustrisTNG. We suggest that the treatment of galaxy evolution in the semi-analytic model needs to be improved by prescriptions that capture local environmental effects more accurately.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2006

Instrument Pointing Control System for the Stellar Interferometry Mission - Planet Quest

This paper describes the high precision Instrument Pointing Control System (PCS) for the Stellar Interferometry Mission (SIM) - Planet Quest. The PCS system provides front-end pointing, compensation for spacecraft motion, and feedforward stabilization, which are needed for proper interference. Optical interferometric measurements require very precise pointing (0.03 as, 1-(sigma) radial) for maximizing the interference pattern visibility. This requirement is achieved by fine pointing control of articulating pointing mirrors with feedback from angle tracking cameras. The overall pointing system design concept is presentcd. Functional requirements and an acquisition concept are given. Guide and Science pointing control loops are discussed. Simulation analyses demonstrate the feasibility of the design.

control↗
NASA NTRS2020

The Warm Gas in the MW: A Kinematical Model

We develop a kinematical model for the Milky Way Si IV-bearing gas to determine its density distribution and kinematics. This model is constrained by a column density line-shape sample extracted from the Hubble Space Telescope/Cosmic Origins Spectrograph archival data, which contains 186 active galactic nucleus sight lines. We find that the Si IV ion density distribution is dominated by an extended disk along the z-direction (above or below the midplane), i.e., n(z) = n0 exp(-(z z/0 ) 0.82 , where z0 is the scale height of - 6.3+1.5 -1.6 kpc (northern hemisphere) and 3.6+0.9 -1.0 kpc (southern hemisphere). The density distribution of the disk in the radial direction shows a sharp edge at 15–20 kpc given by, n(rXY) = n0 exp(-(rXY/r0) 3.36 , where r0≈12.5±0.6 kpc. The difference of density distributions over rXY and z directions indicates that the warm gas traced by Si IV is mainly associated with disk processes (e.g., feedback or cycling gas) rather than accretion. We estimate the mass of the warm gas (within 50 kpc) is log(M(50kpc) M(sun)) ≈ 8.1 (assuming Z≈0.5 Z(sun)), and a 3σ upper limit of log(M(250kpc) M(sun)) ≈ 9.1 (excluding the Magellanic system). Kinematically, the warm gas disk is nearly co-rotating with the stellar disk at vrot = 215 ± 3 km s-1, which lags the midplane rotation by about 8 km s-1 kpc-1 (within 5 kpc). Meanwhile, we note that the warm gas in the northern hemisphere has significant accretion with vacc of 69±7 km s-1 at 10 kpc (an accretion rate of -0.60+0.13- 0.11 M(sun) yr-1), while in the southern hemisphere, there is no measurable accretion, with an upper limit of 0.4M(sun) yr−1.

Zhijie Qu↗
NASA NTRS2012

Evidence for Ultra-Fast Outflows in Radio-Quiet AGNs: III - Location and Energetics

Using the results of a previous X-ray photo-ionization modelling of blue-shifted Fe K absorption lines on a sample of 42 local radio-quiet AGNs observed with XMM-Newton, in this letter we estimate the location and energetics of the associated ultrafast outflows (UFOs). Due to significant uncertainties, we are essentially able to place only lower/upper limits. On average, their location is in the interval approx.0.0003-0.03pc (approx.10(exp 2)-10(exp 4)tau(sub s) from the central black hole, consistent with what is expected for accretion disk winds/outflows. The mass outflow rates are constrained between approx.0.01- 1 Stellar Mass/y, corresponding to approx. or >5-10% of the accretion rates. The average lower-upper limits on the mechanical power are logE(sub K) approx. or = 42.6-44.6 erg/s. However, the minimum possible value of the ratio between the mechanical power and bolometric luminosity is constrained to be comparable or higher than the minimum required by simulations of feedback induced by winds/outflows. Therefore, this work demonstrates that UFOs are indeed capable to provide a significant contribution to the AGN r.osmological feedback, in agreement with theoretical expectations and the recent observation of interactions between AGN outflows and the interstellar medium in several Seyferts galaxies .

Tombesi, F.↗
DOE OSTI2023

The State of the Molecular Gas in Post-starburst Galaxies

The molecular gas in galaxies traces both the fuel for star formation and the processes that can enhance or suppress star formation. Observations of the molecular gas state can thus point to when and why galaxies stop forming stars. In this study, we present Atacama Large Millimeter/submillimeter Array observations of the molecular gas in galaxies evolving through the post-starburst phase. These galaxies have low current star formation rates (SFRs), regardless of the SFR tracer used, with recent starbursts ending within the last 600 Myr. We present CO (3–2) observations for three post-starburst galaxies, and dense gas HCN/HCO+/HNC (1–0) observations for six (four new) post-starburst galaxies. The post-starbursts have low excitation traced by the CO spectral-line energy distribution up to CO (3–2), more similar to early-type than starburst galaxies. The low excitation indicates that lower density rather than high temperatures may suppress star formation during the post-starburst phase. One galaxy displays a blueshifted outflow traced by CO (3–2). MaNGA observations show that the ionized gas velocity is disturbed relative to the stellar velocity field, with a blueshifted component aligned with the molecular gas outflow, suggestive of a multiphase outflow. Low ratios of HCO + /CO, indicating low fractions of dense molecular gas relative to the total molecular gas, are seen throughout post-starburst phase, except for the youngest post-starburst galaxy considered here. These observations indicate that the impact of any feedback or quenching processes may be limited to low excitation and weak outflows in the cold molecular gas during the post-starburst phase.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1993

The evolution of the intergalactic medium and the origin of the galaxy luminosity function

The coupling of the Press and Schechter prescription with the CDM scenario and the Hoyle-Rees-Ostriker cooling criterion leads to a galaxy formation scenario in which galaxies are overproduced by a large factor. Although star formation might be suppressed in the smaller halos, a large amount of energy per galactic mass is needed to account for the present number density of galaxies. The evolution of the intergalactic medium (IGM) provides a simple criterion to prevent galaxy formation without requiring feedback, since halos with small virial temperatures are not able to retain the infalling hot gas of the IGM. If the ionizing background has decreased since z is approximately 1 - 2, then this criterion explains the slope of the luminosity function at the faint end. In addition, this scenario predicts two populations of dwarf galaxies, well differentiated in age, gas content, stellar populations, and clustering properties, which can be identified with dE and dIm galaxies.

Valls-Gabaud, David↗
NASA NTRS2013

Low Masses and High Redshifts: The Evolution of the Mass-Metallicity Relation

We present the first robust measurement of the high redshift mass-metallicity (MZ) relation at 10(exp 8) < M/Stellar Mass < or approx. 10(exp 10), obtained by stacking spectra of 83 emission-line galaxies with secure redshifts between 1.3 < or approx. z < or approx. 2.3. For these redshifts, infrared grism spectroscopy with the Hubble Space Telescope Wide Field Camera 3 is sensitive to the R23 metallicity diagnostic: ([O II] (lambda)(lambda)3726, 3729 + [OIII] (lambda)(lambda)4959, 5007)/H(beta). Using spectra stacked in four mass quartiles, we find a MZ relation that declines significantly with decreasing mass, extending from 12+log(O/H) = 8.8 at M = 10(exp 9.8) Stellar Mass to 12+log(O/H)= 8.2 at M = 10(exp 8.2) Stellar Mass. After correcting for systematic offsets between metallicity indicators, we compare our MZ relation to measurements from the stacked spectra of galaxies with M > or approx. 10(exp 9.5) Stellar Mass and z approx. 2.3. Within the statistical uncertainties, our MZ relation agrees with the z approx. 2.3 result, particularly since our somewhat higher metallicities (by around 0.1 dex) are qualitatively consistent with the lower mean redshift (z = 1.76) of our sample. For the masses probed by our data, the MZ relation shows a steep slope which is suggestive of feedback from energy-driven winds, and a cosmological downsizing evolution where high mass galaxies reach the local MZ relation at earlier times. In addition, we show that our sample falls on an extrapolation of the star-forming main sequence (the SFR-M∗ relation) at this redshift. This result indicates that grism emission-line selected samples do not have preferentially high star formation rates (SFRs). Finally, we report no evidence for evolution of the mass-metallicity-SFR plane; our stack-averaged measurements show excellent agreement with the local relation.

Hgih Redshifts↗
DOE OSTI2024

Constraining neutrino-DM interactions with Milky Way dwarf spheroidals and supernova neutrinos

We constrain the neutrino-dark matter cross section using properties of the dark matter density profiles of Milky Way dwarf spheroidal galaxies. The constraint arises from core-collapse supernova neutrinos scattering on dark matter as a form of energy injection, allowing the transformation of the dark matter density profile from a cusped profile to a flatter profile. We assume a standard cosmology of dark energy and cold, collisionless, and non-self-interacting dark matter. By requiring that the dark matter cores do not lose too much mass or overshoot constraints from stellar kinematics, we place an upper limit on the cross section of σ ν − DM ( E ν = 15 MeV , m χ ≲ 130 GeV ) ≈ 3.4 × 10 − 23 cm 2 and σ ν − DM ( E ν = 15 MeV , m χ ≳ 130 GeV ) ≈ 3.2 × 10 − 27 ( m χ 1 GeV ) 2 cm 2 , which is stronger than previous bounds for these energies. Consideration of baryonic feedback or host galaxy effects on the dark matter profile can strengthen this constraint. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
NASA NTRS2019

Enhanced Habitability on High Obliquity Bodies near the Outer Edge of the Habitable Zone of Sun-Like Stars

High obliquity planets represent potentially extreme limits of terrestrial climate, as they exhibit large seasonality, a reversed annual-mean pole-to-equator gradient of stellar heating, and novel cryospheres. A suite of 3D global climate model simulations is performed for low and high obliquity planets with various stellar fluxes, CO2 concentrations, and initial conditions to explore the propensity for high obliquity climates to undergo global glaciation. We also simulate planets with thick CO2 or H2 atmospheres, such as those expected to develop near or beyond the outer edge of the habitable zone. We show that high obliquity planets are hotter than their low obliquity counterparts due to ice-albedo feedbacks for cold climates, and water vapor in warm climates. We suggest that the water vapor greenhouse trapping is greater on high obliquity bodies for a given global-mean temperature due to the different dynamical regimes that occur between the two states. While equatorial ice belts are stable at high obliquity in some climate regimes, it is substantially harder to achieve global glaciation than for a low obliquity planet. Temperate polar conditions can be present at high obliquity at forcings for which low obliquity planets would be in a hard snowball state. Furthermore, open ocean can persist even in the winter hemisphere and when global-mean temperatures are well below freezing. However, the influence of obliquity diminishes for dense atmospheres, in agreement with calculations from 1D energy balance models.

Planet and satellite atmospheres↗