Search NASA⌕ Search

SEARCH · Search NASA

Results for “(CHEM)”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Rapid water dynamics structures the OH-stretching spectra of solitary water in ionic liquids and dipolar solvents

In a recent study [J. Phys. Chem. B 126, 4584–4598 (2022)], we have used infrared spectroscopy to investigate the solvation and dynamics of solitary water in ionic liquids and dipolar solvents. Complex shapes observed for water OH-stretching bands, common to all high-polarity solvents, were assigned to water in several solvation states. Here, in the present study, classical molecular dynamics simulations of a single water molecule in four ionic liquids and three dipolar solvents were used to test and refine this interpretation. Consistent with past assignments, simulations show solitary water usually donates two hydrogen bonds to distinct solvent molecules. Such symmetrically solvated water produces the primary pair of peaks identified in the OH spectra of water in nearly all solvents. We had further proposed that additional features flanking this main peak are due to asymmetric solvation states, states in which only one OH group makes a hydrogen bond to solvent. Such states were found in significant concentrations in all of the systems simulated. Simulations of the OH stretching spectra using a semiclassical description and the vibrational map developed by Auer and Skinner [J. Chem. Phys. 128, 224511–224512 (2008)] provided semi-quantitative agreement with experiment. Analysis of species-specific spectra confirmed assignment of the additional features in the experimental spectra to asymmetrically solvated water. The simulations also showed that rapid water motions cause a marked motional narrowing compared with the inhomogeneous limit. This narrowing is largely responsible for making the additional features due to minority solvation states manifest in the spectra.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Graph-based quantum response theory and shadow Born–Oppenheimer molecular dynamics

Graph-based linear scaling electronic structure theory for quantum-mechanical molecular dynamics simulations [A. M. N. Niklasson et al., J. Chem. Phys. 144, 234101 (2016)] is adapted to the most recent shadow potential formulations of extended Lagrangian Born–Oppenheimer molecular dynamics, including fractional molecular-orbital occupation numbers [A. M. N. Niklasson, J. Chem. Phys. 152, 104103 (2020) and A. M. N. Niklasson, Eur. Phys. J. B 94, 164 (2021)], which enables stable simulations of sensitive complex chemical systems with unsteady charge solutions. The proposed formulation includes a preconditioned Krylov subspace approximation for the integration of the extended electronic degrees of freedom, which requires quantum response calculations for electronic states with fractional occupation numbers. For the response calculations, we introduce a graph-based canonical quantum perturbation theory that can be performed with the same natural parallelism and linear scaling complexity as the graph-based electronic structure calculations for the unperturbed ground state. Further, the proposed techniques are particularly well-suited for semi-empirical electronic structure theory, and the methods are demonstrated using self-consistent charge density-functional tight-binding theory both for the acceleration of self-consistent field calculations and for quantum-mechanical molecular dynamics simulations. Graph-based techniques combined with the semi-empirical theory enable stable simulations of large, complex chemical systems, including tens-of-thousands of atoms.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Introducing the embedded random phase approximation: H 2 dissociative adsorption on Cu(111) as an exemplar

The random phase approximation (RPA) as a means of treating electron correlation recently has been shown to outperform standard density functional theory (DFT) approximations in a variety of cases. However, the computational cost of the RPA is substantially more than DFT, especially when aiming to study extended surfaces. Properly accounting for sufficient surface ensemble size, Brillouin zone sampling, and vacuum separation of periodic images in standard periodic-planewave-based DFT code raises the cost to achieve converged results. Here, we show that sub-system embedding schemes enable use of the RPA for modeling heterogeneous reactions at reduced computational cost. Further, we explore two different embedded RPA (emb-RPA) approaches, periodic emb-RPA and cluster emb-RPA. We use the (experimentally and theoretically) well-studied H 2 dissociative adsorption on Cu(111) as our exemplar, and first perform full periodic RPA calculations as a benchmark. The full RPA results match well the semi-empirical barrier fit to experimental observables and others derived from high-level computations, e.g., from recent embedded n-electron valence second order perturbation theory [Zhao et al., J. Chem. Theory Comput. 16(11), 7078–7088 (2020)] and quantum Monte Carlo [Doblhoff-Dier et al., J. Chem. Theory Comput. 13(7), 3208–3219 (2017)] simulations. Among the two emb-RPA approaches tested, the cluster emb-RPA accurately reproduces the energy profile (maximum error of 50 meV along the reaction pathway) while reducing the computational cost by approximately two orders of magnitude. We therefore expect that the embedded cluster approach will enable wider RPA implementation in heterogeneous catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Alkali hydroxide (LiOH, NaOH, KOH) in water: Structural and vibrational properties, including neutron scattering results

Structural and vibrational properties of aqueous solutions of alkali hydroxides (LiOH, NaOH, and KOH) are computed using quantum molecular dynamics simulations for solute concentrations ranging between 1 and 10M. Element-resolved partial radial distribution functions, neutron and x-ray structure factors, and angular distribution functions are computed for the three hydroxide solutions as a function of concentration. The vibrational spectra and frequency-dependent conductivity are computed from the Fourier transforms of velocity autocorrelation and current autocorrelation functions. Our results for the structure are validated with the available neutron data for 17M concentration of NaOH in water [Semrouni et al., Phys. Chem. Chem. Phys. 21, 6828 (2019)]. We found that the larger ionic radius [rLi+

Chemistry↗
DOE OSTI2024

A static quantum embedding scheme based on coupled cluster theory

Here, we develop a static quantum embedding scheme that utilizes different levels of approximations to coupled cluster (CC) theory for an active fragment region and its environment. To reduce the computational cost, we solve the local fragment problem using a high-level CC method and address the environment problem with a lower-level Møller–Plesset (MP) perturbative method. This embedding approach inherits many conceptual developments from the hybrid second-order Møller–Plesset (MP2) and CC works by Nooijen [J. Chem. Phys. 111, 10815 (1999)] and Bochevarov and Sherrill [J. Chem. Phys. 122, 234110 (2005)]. We go beyond those works here by primarily targeting a specific localized fragment of a molecule and also introducing an alternative mechanism to relax the environment within this framework. We will call this approach MP-CC. We demonstrate the effectiveness of MP-CC on several potential energy curves and a set of thermochemical reaction energies, using CC with singles and doubles as the fragment solver, and MP2-like treatments of the environment. The results are substantially improved by the inclusion of orbital relaxation in the environment. Using localized bonds as the active fragment, we also report results for N=N bond breaking in azomethane and for the central C–C bond torsion in butadiene. We find that when the fragment Hilbert space size remains fixed (e.g., when determined by an intrinsic atomic orbital approach), the method achieves comparable accuracy with both a small and a large basis set. Additionally, our results indicate that increasing the fragment Hilbert space size systematically enhances the accuracy of observables, approaching the precision of the full CC solver.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

An energy-modified quantum defect method for the analysis of Rydberg spectra: Application to 2-butyne

The high resolution Rydberg absorption spectrum of 2-butyne C 4 H 6 recorded previously at the SOLEIL synchrotron facility has been interpreted using multichannel quantum defect theory (MQDT). The calculations are based on the continuum scattering calculations of Xu et al., J. Chem. Phys. 136, 154303 (2012) and of Jacovella et al., J. Phys. Chem. A 119, 12339 (2015) pertaining to the dipole-allowed excited state symmetries in absorption from the ground state. In contrast to the traditional approach of calculating low-lying electronic states first and then attempting to extend the calculations to ever higher energy, here the analysis proceeds through the extension of these detailed calculations of the electronic continuum scattering down into the discrete region of the spectrum. The continuum reaction matrices and dipole transition moments are adapted to the discrete Rydberg region via the use of an energy-modified formulation of MQDT theory and associated energy dependences of the quantum defects. The analysis reproduces more than 40 Rydberg states from n ≈ 10 down to the 3 d and 4 s levels with an rms error of better than 20 cm −1 . Further, these belong to five Rydberg series with three different molecular symmetries. While the approach predicts many additional series, most of these are calculated and observed to carry only little oscillator strength. The analysis shows that the Rydberg spectrum is dominated by the excitation of an e″ symmetry electron of fδ and gπ type, in line with what previous studies of the above-threshold shape resonance of 2-butyne have shown. The present study is intended to serve as an example showing how first principles continuum calculations may be useful for the interpretation of highly bound discrete states in a range that poses problems for the standard ab initio techniques. The quantitative treatment of the dipole absorption cross sections is deferred to a future paper.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Influence of counterion substitution on the properties of imidazolium-based ionic liquid clusters

Due to their unique physiochemical properties that may be tailored for specific purposes, ionic liquids (ILs) have been investigated for various applications, including chemical separations, catalysis, energy storage, and space propulsion. The different cations and anions comprising ILs may be selected to optimize a range of desired properties, such as thermal stability, ionic conductivity, and volatility, leading to the designation of certain ILs as designer “green” solvents. The effect of counterions on the properties of ILs is of both fundamental scientific interest and technological importance. Herein, we report a systematic experimental and theoretical investigation of the size, charge, stability toward dissociation, and geometric/electronic structure of 1-ethyl-3-methyl imidazolium (EMIM)-based IL clusters containing two different atomic counterions (i.e., bromide [Br − ] and iodide [I − ]). This work extends our studies of EMIM + cations with atomic chloride (Cl − ) and molecular tetrafluoroborate (BF 4 − ) anions reported previously by Baxter et al. [Chem. Mater. 34, 2612 (2022)] and Zhang et al . [J. Phys. Chem. Lett. 11, 6844 (2020)], respectively. Distributions of anionic IL clusters were generated in the gas phase using electrospray ionization and characterized by high mass resolution mass spectrometry, energy-resolved collision-induced dissociation, and negative ion photoelectron spectroscopy experiments. The experimental results reveal anion-dependent trends in the size distribution, relative abundance, ionic charge state, stability toward dissociation, and electron binding energies of the IL clusters. Complementary global optimization theory provides molecular-level insights into the bonding and electronic structure of a selected subset of clusters, including their low energy structures and electrostatic potential maps, and how these fundamental characteristics are influenced by anion substitution. Collectively, our findings demonstrate how the fundamental properties of ILs, which determine their suitability for many applications, may be tuned by substituting counterions. These observations are critical in the sub-nanometer cluster size regime where phenomena do not scale predictably to the bulk phase, and each atom counts toward determining behavior.

cluster↗
DOE OSTI2025

Langevin Dynamics modeling of gas-phase ion-ion recombination (Final Technical Report)

A self-consistent trajectory simulation approach to model MN reactions (Fig. 1) which incorporates the probability of electron transfer as a Monte Carlo operator (Fig. 2) was developed and published as Liu et al. J. Chem. Phys. 159, 114111 (2023). The electron transfer probability p ET estimated using the two-state Landau-Zener (LZ) theory was incorporated into classical trajectory simulations to elicit predictions of MN reaction cross-section σ (vacuum) or rate constant β (finite pressure). Electronic structure calculations with multireference configuration interaction (MRCI) and large correlation consistent basis sets were used to derive inputs to the LZ theory. The key advance of our trajectory simulation approach is the incorporation of electron transfer probability and the inclusion of the effect of ion-neutral interactions on MN using a Langevin representation of the effect of neutral gas on ions. For H + – H - and Li + – H(D) - pairs, our approach quantitatively agrees with measured speed-dependent cross-sections for up to ~10 5 m/s. For the ion pair Ne + – Cl - , our predictions of the MN rate constant at ~1 torr are a factor of ~2 – 3 higher than the experimentally measured value. Similarly, for Xe + – F - in the pressure range of ~20000 – 80000 Pa, our predictions of the MN rate constant are ~20% lower but are in excellent qualitative agreement with experimental data. The paradigm of using trajectory simulations to self-consistently model MN reactions is the basis for inclusion of additional non-classical, and static magnetic and electric field effects. Subsequent work, published as Roy et al. focused on modeling recombination rate constant for three ion pairs (rare gas Ar + cation and halide anions): Ar + – Cl - , Ar + – Br - , Ar + – I - , 2) considering spin-orbit couplings in the electronic structure calculations to obtain high-fidelity estimates of the electron transfer probability and incorporated within the classical trajectory simulations to elicit predictions. In addition to calculations of ion-ion recombination rate constants, a classical trajectory simulation technique (published as Roy et al. J. Chem. Phys. 162(9), 094104 (2023)) that uses quaternions to represent orientation of non-spherical particles (ions or aerosol particles) was developed to simulate the recombination of diatomic or more generally, polyatomic molecules. Finally, several other ion pairs such as Ne + – Cl - , Kr + – Cl - , were explored using the developed semi-classical trajectory simulations to understand various challenges in tackling electronic structure calculations. Using empirical approaches to parameterize the electron transfer radius, trajectory simulations were also used to probe the effect of ion number density on MN rate constant.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Vertical profiles of light absorption and scattering associated with black carbon particle fractions in the springtime Arctic above 79° N

Despite the potential importance of black carbon (BC) for radiative forcing of the Arctic atmosphere, vertically resolved measurements of the particle light scattering coefficient (σ sp ) and light absorption coefficient (σ ap ) in the springtime Arctic atmosphere are infrequent, especially measurements at latitudes at or above 80° N. Here, relationships among vertically distributed aerosol optical properties (σap, σsp and single scattering albedo or SSA), particle microphysics and particle chemistry are examined for a region of the Canadian archipelago between 79.9 and 83.4° N from near the surface to 50 hPa. Airborne data collected during April 2015 are combined with ground-based observations from the observatory at Alert, Nunavut and simulations from the Goddard Earth Observing System (GEOS) model, GEOS-Chem, coupled with the TwO-Moment Aerosol Sectional (TOMAS) model (collectively GEOS-Chem–TOMAS; Kodros et al., 2018) to further our knowledge of the effects of BC on light absorption in the Arctic troposphere. The results are constrained for σsp less than 15 Mm -1 , which represent 98% of the observed σsp, because the single scattering albedo (SSA) has a tendency to be lower at lower σsp, resulting in a larger relative contribution to Arctic warming. At 18.4 m 2 g -1 , the average BC mass absorption coefficient (MAC) from the combined airborne and Alert observations is substantially higher than the two averaged modelled MAC values (13.6 and 9.1 m2 g -1 ) for two different internal mixing assumptions, the latter of which is based on previous observations. The higher observed MAC value may be explained by an underestimation of BC, the presence of small amounts of dust and/or possible differences in BC microphysics and morphologies between the observations and model. In comparing the observations and simulations, we present σap and SSA, as measured, and σ ap 2 and the corresponding SSA to encompass the lower modelled MAC that is more consistent with accepted MAC values. Median values of the measured σap, rBC and the organic component of particles all increase by a factor of 1.8±0.1, going from near-surface to 750 hPa, and values higher than the surface persist to 600 hPa. Modelled BC, organics and σap agree with the near-surface measurements but do not reproduce the higher values observed between 900 and 600 hPa. The differences between modelled and observed optical properties follow the same trend as the differences between the modelled and observed concentrations of the carbonaceous components (black and organic). Model-observation discrepancies may be mostly due to the modelled ejection of biomass burning particles only into the boundary layer at the sources. For the assumption of the observed MAC value, the SSA range between 0.88 and 0.94, which is significantly lower than other recent estimates for the Arctic, in part reflecting the constraint of σsp<15 Mm -1 . The large uncertainties in measuring optical properties and BC, and the large differences between measured and modelled values here and in the literature, argue for improved measurements of BC and light absorption by BC and more vertical profiles of aerosol chemistry, microphysics and other optical properties in the Arctic.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Modeling the smoky troposphere of the southeast Atlantic: a comparison to ORACLES airborne observations from September of 2016

The southeast Atlantic is home to well-defined smoke outflow from Africa coinciding vertically with extensive marine boundary-layer cloud decks, both reaching their climatological maxima in spatial extent around September. A framework is put forth for evaluating the performance of a range of global and regional aerosol models against observations made during the NASA ORACLES (ObseRvations of Aerosols above CLouds and their intEractionS) airborne mission in September 2016. The sparse airborne observations are first aggregated into 2o grid boxes and into three vertical layers: the cloud-topped marine boundary layer (MBL), the layer from cloud top to 3 km, and the 3-6 km layer. Aerosol extensive properties simulated for the entire study region for all September suggest that the 2016 ORACLES observations are reasonably representative of the regional monthly average, with systematic deviations of 30% or less. All six models typically place the bottom of the smoke layer at lower altitudes than do the airborne lidar observations by 300-1400 m, whereas model aerosol top heights are within 0-500 m of the observations. All but one of the models that report carbonaceous aerosol masses underestimate the ratio of particulate extinction to the masses, a proxy for mass extinction efficiency, in 3-6 km. Notable findings on individual models include that WRF-CAM5 predicts the mass of black carbon and organic aerosols with minor (~10% or less) biases. GEOS-5 overestimates the carbonaceous particle masses in the MBL by a factor of 3-6. Extinction coefficients in the free troposphere (FT) and above-cloud aerosol optical depth (ACAOD) are 10-30% lower in WRF-CAM5, 30-50% lower in GEOS-5, 10-40% higher in GEOS-Chem, 10-20% higher in EAM-E3SM except for the practically unbiased 3-6 km extinction, and 20-70% lower in the Unified Model, than the airborne in situ, lidar and sunphotometer measurements. ALADIN-Climate also underestimates the ACAOD, by 30%. GEOS-5 and GEOS-Chem predict carbon monoxide in the MBL with small (10% or less) negative biases, despite their overestimates of carbonaceous aerosol masses. Overall, this study highlights a new approach to utilizing airborne aerosol measurements for model diagnosis.

Shinozuka, Yohei↗
DOE OSTI2020

Air quality impact of the Northern California Camp Fire of November 2018

The Northern California Camp Fire that took place in November 2018 was one of the most damaging environmental events in California history. Here, we analyze ground-based station observations of airborne particulate matter that has a diameter <2.5 µm (PM 2.5 ) across Northern California and conduct numerical simulations of the Camp Fire using the Weather Research and Forecasting model online coupled with chemistry (WRF-Chem). Simulations are evaluated against ground-based observations of PM 2.5 , black carbon, and meteorology, as well as satellite measurements, such as Tropospheric Monitoring Instrument (TROPOMI) aerosol layer height and aerosol index. The Camp Fire led to an increase in Bay Area PM 2.5 to over 50 µg m –3 for nearly 2 weeks, with localized peaks exceeding 300 µg m –3 . Using the Visible Infrared Imaging Radiometer Suite (VIIRS) high-resolution fire detection products, the simulations reproduce the magnitude and evolution of surface PM 2.5 concentrations, especially downwind of the wildfire. The overall spatial patterns of simulated aerosol plumes and their heights are comparable with the latest satellite products from TROPOMI. WRF-Chem sensitivity simulations are carried out to analyze uncertainties that arise from fire emissions, meteorological conditions, feedback of aerosol radiative effects on meteorology, and various physical parameterizations, including the planetary boundary layer model and the plume rise model. Downwind PM 2.5 concentrations are sensitive to both flaming and smoldering emissions over the fire, so the uncertainty in the satellite-derived fire emission products can directly affect the air pollution simulations downwind. Our analysis also shows the importance of land surface and boundary layer parameterization in the fire simulation, which can result in large variations in magnitude and trend of surface PM 2.5 . Inclusion of aerosol radiative feedback moderately improves PM 2.5 simulations, especially over the most polluted days. Results of this study can assist in the development of data assimilation systems as well as air quality forecasting of health exposures and economic impact studies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Understanding and improving model representation of aerosol optical properties for a Chinese haze event measured during KORUS-AQ

KORUS-AQ was an international cooperative air quality field study in South Korea that measured local and remote sources of air pollution affecting the Korean Peninsula during May–June 2016. Some of the largest aerosol mass concentrations were measured during a Chinese haze transport event (24 May). Air quality forecasts using the WRF-Chem model with aerosol optical depth (AOD) data assimilation captured AOD during this pollution episode but overpredicted surface particulate matter concentrations in South Korea, especially PM 2.5 , often by a factor of 2 or larger. Analysis revealed multiple sources of model deficiency related to the calculation of optical properties from aerosol mass that explain these discrepancies. Using in situ observations of aerosol size and composition as inputs to the optical properties calculations showed that using a low-resolution size bin representation (four bins) underestimates the efficiency with which aerosols scatter and absorb light (mass extinction efficiency). Besides using finer-resolution size bins (8–16 bins), it was also necessary to increase the refractive indices and hygroscopicity of select aerosol species within the range of values reported in the literature to achieve better consistency with measured values of the mass extinction efficiency (6.7 m 2 g –1 observed average) and light-scattering enhancement factor (ƒ(RH)) due to aerosol hygroscopic growth (2.2 observed average). Furthermore, an evaluation of the optical properties obtained using modeled aerosol properties revealed the inability of sectional and modal aerosol representations in WRF-Chem to properly reproduce the observed size distribution, with the models displaying a much wider accumulation mode. Other model deficiencies included an underestimate of organic aerosol density (1.0 g cm –3 in the model vs. observed average of 1.5 g cm –3 ) and an overprediction of the fractional contribution of submicron inorganic aerosols other than sulfate, ammonium, nitrate, chloride, and sodium corresponding to mostly dust (17 %–28 % modeled vs. 12 % estimated from observations). These results illustrate the complexity of achieving an accurate model representation of optical properties and provide potential solutions that are relevant to multiple disciplines and applications such as air quality forecasts, health impact assessments, climate projections, solar power forecasts, and aerosol data assimilation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

The mechanisms and seasonal differences of the impact of aerosols on daytime surface urban heat island effect

The urban heat island intensity (UHII) is the temperature difference between urban areas and their rural surroundings. It is commonly attributed to changes in the underlying surface structure caused by urbanization. Air pollution caused by aerosol particles can affect the UHII through changing (1) the surface energy balance by the aerosol radiative effect (ARE) and (2) planetary-boundary-layer (PBL) stability and airflow intensity by modifying thermodynamic structure, which is referred to as the aerosol dynamic effect (ADE). By analyzing satellite data and ground-based observations collected from 2001 to 2010 at 35 cities in China and using the WRF-Chem model, we find that the impact of aerosols on UHII differs considerably: reducing the UHII in summer but increasing the UHII in winter. This seasonal contrast is proposed to be caused by the different strengths of the ARE and ADE between summer and winter. In summer, the ARE on UHII is dominant over the ADE, cooling down surface temperature more strongly in urban areas than in rural areas because of much higher aerosol loading, and offsets the urban heating, therefore weakening UHII. In winter, however, the ADE is more dominant, because aerosols stabilize the PBL more in the polluted condition, weakening the near-surface heat transport over urban areas in both vertical and horizontal directions. This means that the heat accumulated in urban areas is dispersed less effectively, and thus the UHII is enhanced. These findings shed new light on the impact of the interaction between urbanization-induced surface changes and air pollution on urban climate.urban areas and their rural surroundings. It is commonly attributed to changes in the underlying surface structure caused by urbanization. Air pollution caused by aerosol particles can affect the UHII through changing (1) the surface energy balance by the aerosol radiative effect (ARE) and (2) planetary-boundary-layer (PBL) stability and airflow intensity by modifying thermodynamic structure, which is referred to as the aerosol dynamic effect (ADE). By analyzing satellite data and ground-based observations collected from 2001 to 2010 at 35 cities in China and using the WRF-Chem model, we find that the impact of aerosols on UHII differs considerably: reducing the UHII in summer but increasing the UHII in winter. This seasonal contrast is proposed to be caused by the different strengths of the ARE and ADE between summer and winter. In summer, the ARE on UHII is dominant over the ADE, cooling down surface temperature more strongly in urban areas than in rural areas because of much higher aerosol loading, and offsets the urban heating, therefore weakening UHII. In winter, however, the ADE is more dominant, because aerosols stabilize the PBL more in the polluted condition, weakening the near-surface heat transport over urban areas in both vertical and horizontal directions. This means that the heat accumulated in urban areas is dispersed less effectively, and thus the UHII is enhanced. These findings shed new light on the impact of the interaction between urbanization-induced surface changes and air pollution on urban climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Spatial distributions of X CO 2 seasonal cycle amplitude and phase over northern high-latitude regions

Satellite-based observations of atmospheric carbon dioxide (CO 2 ) provide measurements in remote regions, such as the biologically sensitive but undersampled northern high latitudes, and are progressing toward true global data coverage. Recent improvements in satellite retrievals of total column-averaged dry air mole fractions of CO 2 (X CO 2 ) from the NASA Orbiting Carbon Observatory 2 (OCO-2) have allowed for unprecedented data coverage of northern high-latitude regions, while maintaining acceptable accuracy and consistency relative to ground-based observations, and finally providing sufficient data in spring and autumn for analysis of satellite-observed X CO 2 seasonal cycles across a majority of terrestrial northern high-latitude regions. Here, we present an analysis of X CO 2 seasonal cycles calculated from OCO-2 data for temperate, boreal, and tundra regions, subdivided into 5° latitude by 20° longitude zones. We quantify the seasonal cycle amplitudes (SCAs) and the annual half drawdown day (HDD). OCO-2 SCAs are in good agreement with ground-based observations at five high-latitude sites, and OCO-2 SCAs show very close agreement with SCAs calculated for model estimates of X CO 2 from the Copernicus Atmosphere Monitoring Services (CAMS) global inversion-optimized greenhouse gas flux model v19r1 and the CarbonTracker2019 model (CT2019B). Model estimates of X CO 2 from the GEOS-Chem CO 2 simulation version 12.7.2 with underlying biospheric fluxes from CarbonTracker2019 (GC-CT2019) yield SCAs of larger magnitude and spread over a larger range than those from CAMS, CT2019B, or OCO-2; however, GC-CT2019 SCAs still exhibit a very similar spatial distribution across northern high-latitude regions to that from CAMS, CT2019B, and OCO-2. Zones in the Asian boreal forest were found to have exceptionally large SCA and early HDD, and both OCO-2 data and model estimates yield a distinct longitudinal gradient of increasing SCA from west to east across the Eurasian continent. In northern high-latitude regions, spanning latitudes from 47 to 72° N, longitudinal gradients in both SCA and HDD are at least as pronounced as latitudinal gradients, suggesting a role for global atmospheric transport patterns in defining spatial distributions of X CO2 seasonality across these regions. GEOS-Chem surface contact tracers show that the largest X CO 2 SCAs occur in areas with the greatest contact with land surfaces, integrated over 15–30d. The correlation of X CO 2 SCA with these land surface contact tracers is stronger than the correlation of X CO 2 SCA with the SCA of CO 2 fluxes or the total annual CO 2 flux within each 5° latitude by 20° longitude zone. This indicates that accumulation of terrestrial CO 2 flux during atmospheric transport is a major driver of regional variations in X CO 2 SCA.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Aerosols from anthropogenic and biogenic sources and their interactions – modeling aerosol formation, optical properties, and impacts over the central Amazon basin

The Green Ocean Amazon experiment – GoAmazon 2014–2015 – explored the interactions between natural biogenic forest emissions from central Amazonia and urban air pollution from Manaus. Previous GoAmazon 2014–2015 studies showed that nitrogen oxide (NO x = NO + NO 2 ) and sulfur oxide (SO x ) emissions from Manaus strongly interact with biogenic volatile organic compounds (BVOCs), affecting secondary organic aerosol (SOA) formation. In previous studies, ground-based and aircraft measurements provided evidence of SOA formation and strong changes in aerosol composition and properties. Aerosol optical properties also evolve, and their impacts on the Amazonian ecosystem can be significant. As particles age, some processes, such as SOA production, black carbon (BC) deposition, particle growth and the BC lensing effect change the aerosol optical properties, affecting the solar radiation flux at the surface. This study analyzes data and models SOA formation using the Weather Research and Forecasting with Chemistry (WRF-Chem) model to assess the spatial variability in aerosol optical properties as the Manaus plumes interact with the natural atmosphere. The following aerosol optical properties are investigated: single scattering albedo (SSA), asymmetry parameter (g aer ), absorption Ångström exponent (AAE) and scattering Ångström exponent (SAE). These simulations were validated using ground-based measurements at three experimental sites, namely the Amazon Tall Tower Observatory – ATTO (T0a), downtown Manaus (T1), Tiwa Hotel (T2) and Manacapuru (T3), as well as the U.S. Department of Energy (DOE) Gulfstream 1 (G-1) aircraft flights. WRF-Chem simulations were performed over 7 d during March 2014. Results show a mean biogenic SOA (BSOA) mass enrichment of 512 % at the T1 site, 450 % in regions downwind of Manaus, such as the T3 site, and 850 % in areas north of the T3 site in simulations with anthropogenic emissions. The SOA formation is rather fast, with about 80 % of the SOA mass produced in 3–4 h. Comparing the plume from simulations with and without anthropogenic emissions, SSA shows a downwind reduction of approximately 10 %, 11 % and 6 % at the T1, T2 and T3 sites, respectively. Other regions, such as those further downwind of the T3 site, are also affected. The g aer values increased from 0.62 to 0.74 at the T1 site and from 0.67 to 0.72 at the T3 site when anthropogenic emissions are active. During the Manaus plume-aging process, a plume tracking analysis shows an increase in SSA from 0.91 close to Manaus to 0.98 160 km downwind of Manaus as a result of SOA production and BC deposition.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

North China Plain as a hot spot of ozone pollution exacerbated by extreme high temperatures

Abstract. A large population in China has been exposed to both severe ozone (O3) pollution and extreme heat under global warming. Here, the spatiotemporal characteristics of coupled extremes in surface O3 and heat (OPCs) over China are investigated using surface observations, a process-based chemical transport model (GEOS-Chem), and multi-model simulations from Phase 6 of the Coupled Model Intercomparison Project (CMIP6). North China Plain (NCP; 37–41∘ N; 114–120∘ E) is identified as a hot spot of OPCs, where more than half of the O3 pollution days are accompanied by high temperature extremes. OPCs over NCP exceeded 40 d during 2014–2019, exhibiting an increasing trend. Both O3 concentrations and temperatures are elevated during OPCs compared with O3 pollution days occurring individually (OPIs). Therefore, OPCs impose more severe health impacts to humans than OPIs, but the stronger health effects are mainly driven by the higher temperatures. GEOS-Chem simulations further reveal that enhanced chemical production resulting from hot and stable atmospheric conditions under anomalous weather patterns primarily contributes to the exacerbated O3 levels during OPCs. In the future, CMIP6 projections suggest increased occurrences of OPCs over NCP in the middle of this century, but by the end of this century, OPCs may decrease or increase depending on the pollutant emission scenarios. However, for all future scenarios, extreme high temperatures will play an increasingly important role in modulating O3 pollution in a warming climate.

Wang, Pinya (ORCID:0000000170347868)↗
DOE OSTI2022

Global simulations of monoterpene-derived peroxy radical fates and the distributions of highly oxygenated organic molecules (HOMs) and accretion products

Abstract. We evaluate monoterpene-derived peroxy radical (MT-RO2) unimolecular autoxidation and self- and cross-reactions with other RO2 species in the GEOS-Chem global chemical transport model. The formation of associated highly oxygenated organic molecules (HOMs) and accretion products are tracked in competition with other bimolecular reactions. Autoxidation is the dominant fate up to 6–8 km for first-generation MT-RO2, which can undergo unimolecular H shifts. Reaction with NO can be a more common fate for H-shift rate constants < 0.1 s−1 or at altitudes higher than 8 km due to the imposed Arrhenius temperature dependence of unimolecular H shifts. For MT-derived HOM-RO2, generated by multistep autoxidation of first-generation MT-RO2, reaction with other RO2 species is predicted to be the major fate throughout most of the boreal and tropical forest regions, whereas reaction with NO dominates in the temperate and subtropical forests of the Northern Hemisphere. The newly added reactions result in an approximate 4 % global average decrease in HO2 and RO2, mainly due to faster self-/cross-reactions of MT-RO2, but the impact upon HO2, OH, and NOx abundances is only important in the planetary boundary layer (PBL) over portions of tropical forests. Predicted HOM concentrations in MT-rich regions and seasons can exceed total organic aerosol predicted by the standard version of the GEOS-Chem model depending on the parameters used. Comparisons to observations reveal that large uncertainties remain for key reaction parameters and processes, especially with respect to the photochemical lifetime and volatility of HOMs as well as the rates and branching of associated RO2 accretion products. Further observations and laboratory studies related to MT-RO2-derived HOMs and gas-phase RO2 accretion product formation kinetics – especially their atmospheric fate, such as gas–particle partitioning, multiphase chemistry, and net secondary organic aerosol formation – are needed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Contribution of expanded marine sulfur chemistry to the seasonal variability of dimethyl sulfide oxidation products and size-resolved sulfate aerosol

Marine emissions of dimethyl sulfide (DMS) and the subsequent formation of its oxidation products methanesulfonic acid (MSA) and sulfuric acid (H 2 SO 4 ) are well-known natural precursors of atmospheric aerosols, contributing to particle mass and cloud formation over ocean and coastal regions. Despite a long recognized and well-studied role in the marine troposphere, DMS oxidation chemistry remains a work in progress within many current air quality and climate models, with recent advances exploring heterogeneous chemistry and uncovering previously unknown intermediate species. With the identification of additional DMS oxidation path ways and intermediate species that influence the eventual fate of DMS, it is important to understand the impact of these pathways on the overall sulfate aerosol budget and aerosol size distribution. In this work, we update and evaluate the DMS oxidation mechanism of the chemical transport model GEOS-Chem by implementing expanded DMS oxidation pathways in the model. These updates include gas- and aqueous-phase reactions, the formation of the intermediates dimethyl sulfoxide (DMSO) and methanesulfinic acid (MSIA), and cloud loss and aerosol uptake of the recently quantified intermediate hydroperoxymethyl thioformate (HPMTF). We find that this updated mechanism collectively decreases the global mean surface-layer gas-phase sulfur dioxide (SO 2 ) mixing ratio by 40 % and enhances the sulfate aerosol ($SO^{2-}_{4}$) mixing ratio by 17 %. We further perform sensitivity analyses exploring the contribution of cloud loss and aerosol uptake of HPMTF to the overall sulfur budget. Comparing modeled concentrations to available observations, we find improved biases relative to previous studies. To quantify the impacts of these chemistry updates on global particle size distributions and the mass concentration, we use the TwO-Moment Aerosol Sectional (TOMAS) aerosol microphysics module coupled to GEOS-Chem and find that changes in particle formation and growth affect the size distribution of aerosol. With this new DMS-oxidation scheme, the global annual mean surface-layer number concentration of particles with diameters smaller than 80 nm decreases by 16.8 %, with cloud loss processes related to HPMTF being mostly responsible for this reduction. However, the global annual mean number of particles larger than 80 nm (corresponding to particles capable of acting as cloud condensation nuclei, CCN) increases by 3.8 %, suggesting that the new scheme promotes seasonal particle growth to these sizes.

54 ENVIRONMENTAL SCIENCES↗