Abstract. Secondary organic aerosol derived from isopreneepoxydiols (IEPOX-SOA) is thought to contribute the dominant fraction oftotal isoprene SOA, but the current volatility-based lumped SOAparameterizations are not appropriate to represent the reactive uptake ofIEPOX onto acidified aerosols. A full explicit modeling of this chemistryis however computationally expensive owing to the many species and reactionstracked, which makes it difficult to include it in chemistry–climate modelsfor long-term studies. Here we present three simplified parameterizations(version 1.0) for IEPOX-SOA simulation, based on an approximateanalytical/fitting solution of the IEPOX-SOA yield and formation timescale.The yield and timescale can then be directly calculated using the globalmodel fields of oxidants, NO, aerosol pH and other key properties, and drydeposition rates. The advantage of the proposed parameterizations is thatthey do not require the simulation of the intermediates while retaining thekey physicochemical dependencies. We have implemented the newparameterizations into the GEOS-Chem v11-02-rc chemical transport model,which has two empirical treatments for isoprene SOA (the volatility-basis-set, VBS, approach and a fixed 3 % yield parameterization), and comparedall of them to the case with detailed fully explicit chemistry. The bestparameterization (PAR3) captures the global tropospheric burden of IEPOX-SOAand its spatiotemporal distribution (R2=0.94) vs. thosesimulated by the full chemistry, while being more computationally efficient(∼5 times faster), and accurately captures the response tochanges in NOx and SO2 emissions. On the other hand, the constant3 % yield that is now the default in GEOS-Chem deviates strongly (R2=0.66), as does the VBS (R2=0.47, 49 % underestimation), withneither parameterization capturing the response to emission changes. Withthe advent of new mass spectrometry instrumentation, many detailed SOAmechanisms are being developed, which will challenge global and especiallyclimate models with their computational cost. The methods developed in thisstudy can be applied to other SOA pathways, which can allow includingaccurate SOA simulations in climate and global modeling studies in thefuture.
more »
« less
Active Site Engineering via Optimizing the Heterogeneous Support Structure for Single-Atom Catalysis
Supported single-atom catalysts show a large range of activities and selectivities that depend on the local environment of the catalytic sites. A theory-based optimization strategy is presented that is based on a density functional theory determination of the transition states and intermediates for a low-dimensional coordinate representation of the heterogeneity of the active sites. The approach is applied to a vanadium catalyst on an amorphous SiO2 support that involves a large kinetic network described using a full chemistry model. Without assuming a priori scaling relations or mechanism reduction, the optimal state of heterogeneity is found to lie at atomic configurations where the activation energies for two distinct key chemical processes are equal. It is found a posteriori that the behavior of the system is consistent with linear free energy scaling relations in the randomness parameters. The energetic span theory proves quite useful in reducing the full chemistry model to a small number of key reactions. The use of a nonlinear optimization algorithm in combination with energetic span theory provides significant simplification in treating disordered systems.
more »
« less
- Award ID(s):
- 1664555
- PAR ID:
- 10481306
- Publisher / Repository:
- American Chemical Society
- Date Published:
- Journal Name:
- The Journal of Physical Chemistry C
- Volume:
- 127
- Issue:
- 34
- ISSN:
- 1932-7447
- Page Range / eLocation ID:
- 16901 to 16913
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract In the design of stellarators, energetic particle confinement is a critical point of concern which remains challenging to study from a numerical point of view. Standard Monte Carlo (MC) analyses are highly expensive because a large number of particle trajectories need to be integrated over long time scales, and small time steps must be taken to accurately capture the features of the wide variety of trajectories. Even when they are based on guiding center trajectories, as opposed to full-orbit trajectories, these standard MC studies are too expensive to be included in most stellarator optimization codes. We present the first multifidelity Monte Carlo (MFMC) scheme for accelerating the estimation of energetic particle confinement in stellarators. Our approach relies on a two-level hierarchy, in which a guiding center model serves as the high-fidelity model, and a data-driven linear interpolant is leveraged as the low-fidelity surrogate model. We apply MFMC to the study of energetic particle confinement in a four-period quasi-helically symmetric stellarator, assessing various metrics of confinement. Stemming from the very high computational efficiency of our surrogate model as well as its sufficient correlation to the high-fidelity model, we obtain speedups of up to 10 with MFMC compared to standard MC.more » « less
-
Abstract Whether or not coherent magnetospheric whistler waves play important roles in the pitch‐angle scattering of energetic particles is a crucial question in magnetospheric physics. The interaction of a thermal distribution of energetic particles with coherent whistler waves is thus investigated. The distribution is prescribed by the Maxwell‐Jüttner distribution, which is a relativistic generalization of the Maxwell‐Boltzmann distribution. Coherent whistler waves are modeled by circularly polarized waves propagating parallel to the background magnetic field. It is shown that for parameters relevant to magnetospheric chorus, a significant fraction (1–5%) of the energetic particle population undergoes drastic, nondiffusive pitch‐angle scattering by coherent chorus. The scaling of this fraction with the wave amplitude may also explain the association of relativistic microbursts to large‐amplitude chorus. A much improved condition for large pitch‐angle scattering is presented that is related to, but may or may not include the exact resonance condition depending on the particle's initial conditions. The theory reveals a critical mechanism not contained in the widely used second‐order trapping theory.more » « less
-
The scaling relations between city attributes and population are emergent and ubiquitous aspects of urban growth. Quantifying these relations and understanding their theoretical foundation, however, is difficult due to the challenge of defining city boundaries and a lack of historical data to study city dynamics over time and space. To address this issue, we analyze scaling between city infrastructure and population across 857 metropolitan areas in the conterminous United States over an unprecedented 115 years (1900–2015) using dasymetrically refined historical population estimates, historical urban road network models, and multi-temporal settlement data to define dynamic city boundaries. We demonstrate that urban scaling exponents closely match theoretical models over a century. Despite some close quantitative agreement with theory, the empirical scaling relations unexpectedly vary across regions. Our analysis of scaling coefficients, meanwhile, reveals that contemporary cities use more developed land and kilometers of road than cities of similar population in 1900, which has serious implications for urban development and impacts on the local environment. Overall, our results provide a new way to study urban systems based on novel, geohistorical data.more » « less
-
Second-order Møller–Plesset perturbation theory (MP2) provides a valuable alternative to density functional theory for modeling problems in organic and biological chemistry. However, MP2 suffers from known limitations in the description of van der Waals (London) dispersion interactions and reaction thermochemistry. Here, a spin-component-scaled, dispersion-corrected MP2 model (SCS-MP2D) is proposed that addresses these weaknesses. The dispersion correction, which is based on Grimme's D3 formalism, replaces the uncoupled Hartree–Fock dispersion inherent in MP2 with a more robust coupled Kohn–Sham treatment. The spin-component scaling of the residual MP2 correlation energy then reduces the remaining errors in the model. This two-part correction strategy solves the problem found in earlier spin-component-scaled MP2 models where completely different spin-scaling parameters were needed for describing reaction energies versus intermolecular interactions. Results on 18 benchmark data sets and two challenging potential energy curves demonstrate that SCS-MP2D considerably improves upon the accuracy of MP2 for intermolecular interactions, conformational energies, and reaction energies. Its accuracy and computational cost are competitive with state-of-the-art density functionals such as DSD-BLYP-D3(BJ), revDSD-PBEP86-D3(BJ), ωB97X-V, and ωB97M-V for systems with ∼100 atoms.more » « less
An official website of the United States government

