Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex systems theory”

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 829 records · Page 46

DyG-DPCD: A Distributed Parallel Community Detection Algorithm for Large-Scale Dynamic Graphs

Dynamic (Temporal) graphs capture the valuable evolution of real-world systems, from the continuously evolving patterns of social interactions and genetic pathways to the dynamic fluctuations of economic forces. Detecting communities for such evolving networks poses unique challenges. Detecting and analyzing the evolution of communities within dynamic graphs unlocks valuable insights into the underlying structural and temporal patterns of real-world systems. However, the sheer volume of modern graph data and the inherent complexity of the temporal dimension pose significant challenges to scalable community detection algorithms. Addressing this gap, our work explores the limited landscape of scalable distributed-memory parallel methods specifically designed for dynamic network community detection. We propose a novel parallel algorithm, DyG-DPCD (Dynamic Graph Distributed Parallel Community Detection), to detect communities in dynamic networks using the Message Passing Interface (MPI) framework. We present a vertex-centric approach, allowing us to detect communities through local optimization. Furthermore, we enhance our baseline algorithm by incorporating three heuristics, which improve the algorithm’s performance significantly while maintaining the quality of the solutions. We demonstrate the efficiency of our algorithm by experimenting on several real-world large-scale networks with hundreds of millions of edges spanning diverse domains. Notably, DyG-DPCD achieves speedups between 25× and 30× for large networks that we experimented on using NERSC compute nodes. In conclusion, our algorithm outperforms the STINGER parallel re-agglomeration algorithm by 30×.

97 MATHEMATICS AND COMPUTING↗

Ripening of Rh Nanoparticle Catalysts in Reverse Water–Gas Shift via a Data-Driven Model Combining Physics, Theory, and Experiment

Degradation via sintering is an ongoing challenge that impedes the broad commercial success of supported metallic nanoparticle catalysts. To mitigate degradation via informed catalyst design and process operations, here we aim to disambiguate the underlying mechanisms of sintering by combining theory and experiment in a quantitative framework. While mechanistic sintering models exist, they only model a single sintering pathway, even though multiple sintering mechanisms can occur simultaneously or dominate at different stages of the process. Data-driven machine learning models have emerged as a means to represent complex processes through data regression. However, machine learning models have very large data needs and lack mechanistic insights due to their black-box encoding. To develop an interpretive model of catalyst degradation via sintering, we constructed a hybrid model combining mechanistic “physics-based” models and data-driven methods to obtain both reliable predictions and mechanistic insights regarding experimentally observed sintering phenomena. Focusing on nanoparticle sintering in the Rh–TiO 2 catalyst for the reverse water–gas shift (RWGS) reaction, the hybrid model couples a mechanistic term for Ostwald ripening with energy values calculated via density functional theory (DFT) with a parametric, data-driven discrepancy function term for unmodeled mechanisms. The hybrid model is trained using Bayesian inference with data collected from small-angle X-ray scattering (SAXS) in situ experiments wherein average nanoparticle diameter versus time was measured at three relevant operating temperatures. The calibrated hybrid model results show that an Ostwald ripening-only model parameterized with fixed DFT energies does not fully capture the time and temperature dependence of the SAXS-observed sintering kinetics, and that an additional functional contribution, or DFT energy calibration, is required to reconcile simulation and experiment. Analysis of the hybrid-model error confirms that the hybrid model outperforms both the purely mechanistic and purely data-driven alternatives in terms of expected predictive accuracy for time-evolving average particle sizes. Furthermore, the results support the hypothesis that the Ostwald ripening mechanism is less important for explaining the sintering phenomena as operating temperature increases under an assumed fixed DFT parameterization. This could be explained in one of two ways: either latent, unmodeled sintering mechanisms dominate at higher temperatures, or the DFT uncertainty increases with temperature. The proposed modeling approach directly links theory to experiments and simulations via a statistical hybrid modeling framework and can be extended to other catalytic systems to improve predictive models and mechanistic understanding.

Bayesian hybrid modeling↗

Learning heterogeneous reaction kinetics from X-ray videos pixel by pixel

Reaction rates at spatially heterogeneous, unstable interfaces are notoriously difficult to quantify, yet are essential in engineering many chemical systems, such as batteries and electrocatalysts. Experimental characterizations of such materials by operando microscopy produce rich image datasets, but data-driven methods to learn physics from these images are still lacking because of the complex coupling of reaction kinetics, surface chemistry and phase separation. Here we show that heterogeneous reaction kinetics can be learned from in situ scanning transmission X-ray microscopy (STXM) images of carbon-coated lithium iron phosphate (LFP) nanoparticles. Combining a large dataset of STXM images with a thermodynamically consistent electrochemical phase-field model, partial differential equation (PDE)-constrained optimization and uncertainty quantification, we extract the free-energy landscape and reaction kinetics and verify their consistency with theoretical models. We also simultaneously learn the spatial heterogeneity of the reaction rate, which closely matches the carbon-coating thickness profiles obtained through Auger electron microscopy (AEM). Across 180,000 image pixels, the mean discrepancy with the learned model is remarkably small (<7%) and comparable with experimental noise. Our results open the possibility of learning nonequilibrium material properties beyond the reach of traditional experimental methods and offer a new non-destructive technique for characterizing and optimizing heterogeneous reactive surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Stabilization of [(N 5 ) 2 BX] 2– and [(N 5 ) 2 B 2 X 2 ] 2– (X = H, F, Cl, Br) by Conjugation and Hyperconjugation Effects

The isolation of nucleophilic boron bases has led to a paradigm shift in boron chemistry. Previous studies of the bis(carbene) borylene complexes revealed that these compounds possess strong donor abilities, and their reaction inertness is due to the large steric hindrance between boron reagents and reactant. In the present study, we have theoretically studied the [(N 5 ) 2 BX] 2– and [(N 5 ) 2 B 2 X 2 ] 2– compounds (X = H, F, Cl, Br). Their electronic structures and properties are discussed by using the NBO, LOL, and ELF methods. We found that both π-conjugation and hyperconjugation effects can effectively stabilize the substituted nucleophilic anionic boron compounds [(N 5 ) 2 BX] 2– and [(N 5 ) 2 B 2 X 2 ] 2– . Substituents, especially X = H, stabilize the boron center through highly delocalized π-bonding, involving the formally “empty” in-plane p orbitals of the boron atom. While the halogen substituents have high electron withdrawal ability, leading to systems being less stable, we suggest the borinium anions [(N 5 ) 2 BH] 2– and [(N 5 ) 2 B 2 H 2 ] 2– as possible synthetic targets of novel environmentally friendly catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Incorporating Coverage-Dependent Reaction Barriers into First-Principles-Based Microkinetic Models: Approaches and Challenges

Mean-field microkinetic models (MKMs) are appealing for their relatively facile construction, computational tractability, and high-throughput catalyst screening capabilities. As such, they will continue to be a valuable tool for materials design in heterogeneous catalysis even as the field aims to describe more complex systems. Numerous prior reports have provided the groundwork for constructing first-principles-based MKMs, including the analysis of strategies for incorporating lateral interactions into thermodynamic parameters (e.g., adsorption energies). Yet, there remains a need for concerted dialogue on methods for calculating and incorporating coverage-dependent kinetic parameters into MKMs. In this Perspective, we assess strategies for doing so, including the corresponding key physical implications and computational challenges. Here, we emphasize that decoupling thermodynamic and kinetic parameters within MKMs can violate thermodynamic consistency and risk unphysical solutions. For some reactions and catalyst materials, scaling relationships can predict coverage-dependent activation energies, but there are several exceptions evident in the literature, indicating that this approach is not universally applicable and that the field could benefit from research aimed at elucidating the limitations. Conducting high-coverage transition state searches is a rigorous but computationally costly strategy, and the effects of various methods for mitigating this cost on resulting energetics have yet to be broadly explored and validated. The goal of this Perspective is to generate discussion on and inspire focused research into the physical relevance of approaches for describing coverage-dependent reaction barriers in MKMs, including the development of computationally tractable methodologies, to advance the applicability of MKMs across diverse reaction chemistries and conditions.

36 MATERIALS SCIENCE↗

Site Mixing for Engineering Magnetic Topological Insulators

The van der Waals compound, MnBi 2 Te 4 , is the first intrinsic magnetic topological insulator, providing a materials platform for exploring exotic quantum phenomena such as the axion insulator state and the quantum anomalous Hall effect. However, intrinsic structural imperfections lead to bulk conductivity, and the roles of magnetic defects are still unknown. With higher concentrations of the same types of magnetic defects, the isostructural compound MnSb 2 Te 4 is a better model system for a systematic investigation of the connections among magnetism, topology, and lattice defects. In this work, the impact of antisite defects on the magnetism and electronic structure is studied in MnSb 2 Te 4 . Mn-Sb site mixing leads to complex magnetic structures and tunes the interlayer magnetic coupling between antiferromagnetic and ferromagnetic. The detailed nonstoichiometry and site mixing of MnSb2Te4 crystals depend on the growth parameters, which can lead to ≈40% of Mn sites occupied by Sb and ≈15% of Sb sites by Mn in as-grown crystals. Single-crystal neutron diffraction and electron microscopy studies show nearly random distribution of the antisite defects. Band structure calculations suggest that the Mn-Sb site mixing favors a ferromagnetic interlayer coupling, consistent with experimental observation, but is detrimental to the band inversion required for a nontrivial topology. Overall, our results suggest a long-range magnetic order of Mn ions sitting on Bi sites in MnBi 2 Te 4 . The effects of site mixing should be considered in all layered heterostructures that consist of alternating magnetic and topological layers, including the entire family of MnTe(Bi 2 Te 3 )n, its Sb analogs, and their solid solution.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A Model for Optimally Allocating Curbside Space Among Competing Uses

The emergence of various new forms of urban mobility services in recent years is leading to new pressures on curbside space. Municipalities, the entities typically responsible for managing the curbside, are in many instances handling these growing pressures by reallocating portions of the curbside away from traditional uses (such as metered and residential parking) in favor of uses such as ridehailing, scooter and bike-share corrals. As yet, however, such actions are being undertaken on an ad-hoc basis, due to the rapidly growing complexity of the curbside and the lack of standard analytical approaches. This lack of analytical capability is due to the traditional focus of transportation network modeling being focused predominantly on the interaction of supply and demand on links and nodes, with limited focus on link edges (the curbside). In this paper we address this research need by proposing a framework for modeling inter-modal competition for curbside space, inspired by the classical Bid-Rent Model of urban land use, intended to support curb managers to move towards maximizing the aspects of economic welfare that relate to curb access. In the bi-level model, choices made by the curbside manager impact travelers’ mode choices, and vice versa. We then present a simple numerical case study to demonstrate the properties of the proposed model, showing its tractability, flexibility, and intuitive sensitivity to systematic variation in inputs. The framework demonstrates the type of adaptive and evolving approach needed to maximize benefits from increasingly dynamic curb management strategies. The paper concludes with a brief discussion of future research needs to advance this line of inquiry.

33 ADVANCED PROPULSION SYSTEMS↗

The presence of Ru metal modifies the behavior of deuterium bound to pyridinic nitrogen in doped graphene-like materials

We used a combination of x-ray photoelectron spectroscopy (XPS) and density functional theory (DFT) to investigate how the presence of a ruthenium (Ru) metal substrate modifies the N-D bond strength of the pyridinic nitrogen species in nitrogen (N)-doped graphene (Gr) on Ru(0001) compared to metal-free highly-oriented pyrolytic graphite (HOPG) substrate. The N-dopants were introduced through low-energy N+/N2+ irradiation. N is embedded in the carbon materials in two predominant configurations: as graphitic N (GN: N substituted in the hexagonal C lattice) and pyridinic N (PN: substitutional N adjacent to a C vacancy). XPS data showed that upon atomic deuterium (D) exposure at 220 K, the PN peak shifted to a higher binding energy by +1.2 eV for HOPG and +1.0 eV for N-doped graphene on Ru(0001), while the GN peak remained unchanged, indicating that the D atoms bound solely to pyridinic N. 85% of PN sites on HOPG can be saturated with D atoms, whereas only ~30% of pyridinic N sites are able to bind D atoms in N-doped Gr/Ru(0001). Our DFT analysis shows that this difference is due to the coordination of PN to Ru atoms, which necessitates bond cleavage of the N-Ru interaction prior to D atom adsorption. D begins desorbing from N-HOPG at ~573 K and is fully desorbed by ~973 K, whereas desorption from N-doped graphene on Ru(0001) begins at ~290 K, with complete desorption observed at ~700 K, indicating that the Ru metal weakens the N-D bond strength. We also studied a high-surface-area, layered, and porous N-doped carbon material, internally labeled NC900, which was synthesized by pyrolysis of glucose and graphitic carbon nitride (g-C3N4) at 1173 K. D exposure caused a +1.1 eV PN peak shift in NC900, with no change in the GN peak. Moreover, the NC900 exhibited the same desorption behavior as HOPG, demonstrating that well-defined model systems can effectively capture the behavior of more complex N-doped carbon materials.

Alupothe Gedara, Buddhika S.↗

A coarse-grained model of clay colloidal aggregation and consolidation with explicit representation of the electrical double layer

The aggregation of clay minerals in liquid water exemplifies colloidal self-assembly in nature. These negatively charged aluminosilicate platelets interact through multiple mechanisms with different sensitivities to particle shape, surface charge, aqueous chemistry, and interparticle distance and exhibit complex aggregation structures. Experiments have difficulty resolving the associated colloidal assemblages at the scale of individual particles. Conversely, all-atom molecular dynamics (MD) simulations provide detailed insight on clay colloidal interaction mechanisms, but they are limited to systems containing a few particles. We develop a new coarse-grained (CG) model capable of representing assemblages of hundreds of clay particles with accuracy approaching that of MD simulations, at a fraction of the computational cost. Our CG model is parameterized based on MD simulations of a pair of smectite clay particles in liquid water. A distinctive feature of our model is that it explicitly represents the electrical double layer (EDL), i.e., the cloud of charge-compensating cations that surrounds the clay particles. Our model captures the simultaneous importance of long-range colloidal interactions (i.e., interactions consistent with simplified analytical models, already included in extant clay CG models) and short-range interactions such as ion correlation and surface and ion hydration effects. The resulting simulations correctly predict, at low solid-water ratios, the existence of ordered arrangements of parallel particles separated by water films with a thickness up to ~10 nm and, at high solid-water ratios, the coexistence of crystalline and osmotic swelling states, in agreement with experimental observations.

54 ENVIRONMENTAL SCIENCES↗

INTEGRATE - Inverse Network Transformations for Efficient Generation of Robust Airfoil and Turbine Enhancements

The INTEGRATE (Inverse Network Transformations for Efficient Generation of Robust Airfoil and Turbine Enhancements) project is developing a new inverse-design capability for the aerodynamic design of wind turbine rotors using invertible neural networks. This AI-based design technology can capture complex non-linear aerodynamic effects while being 100 times faster than design approaches based on computational fluid dynamics. This project enables innovation in wind turbine design by accelerating time to market through higher-accuracy early design iterations to reduce the levelized cost of energy. INVERTIBLE NEURAL NETWORKS Researchers are leveraging a specialized invertible neural network (INN) architecture along with the novel dimension-reduction methods and airfoil/blade shape representations developed by collaborators at the National Institute of Standards and Technology (NIST) learns complex relationships between airfoil or blade shapes and their associated aerodynamic and structural properties. This INN architecture will accelerate designs by providing a cost-effective alternative to current industrial aerodynamic design processes, including: - Blade element momentum (BEM) theory models: limited effectiveness for design of offshore rotors with large, flexible blades where nonlinear aerodynamic effects dominate - Direct design using computational fluid dynamics (CFD): cost-prohibitive - Inverse-design models based on deep neural networks (DNNs): attractive alternative to CFD for 2D design problems, but quickly overwhelmed by the increased number of design variables in 3D problems AUTOMATED COMPUTATIONAL FLUID DYNAMICS FOR TRAINING DATA GENERATION - MERCURY FRAMEWORK The INN is trained on data obtained using the University of Marylands (UMD) Mercury Framework, which has with robust automated mesh generation capabilities and advanced turbulence and transition models validated for wind energy applications. Mercury is a multi-mesh paradigm, heterogeneous CPU-GPU framework. The framework incorporates three flow solvers at UMD, 1) OverTURNS, a structured solver on CPUs, 2) HAMSTR, a line based unstructured solver on CPUs, and 3) GARFIELD, a structured solver on GPUs. The framework is based on Python, that is often used to wrap C or Fortran codes for interoperability with other solvers. Communication between multiple solvers is accomplished with a Topology Independent Overset Grid Assembler (TIOGA). NOVEL AIRFOIL SHAPE REPRESENTATIONS USING GRASSMAN SPACES We developed a novel representation of shapes which decouples affine-style deformations from a rich set of data-driven deformations over a submanifold of the Grassmannian. The Grassmannian representation as an analytic generative model, informed by a database of physically relevant airfoils, offers (i) a rich set of novel 2D airfoil deformations not previously captured in the data , (ii) improved low-dimensional parameter domain for inferential statistics informing design/manufacturing, and (iii) consistent 3D blade representation and perturbation over a sequence of nominal shapes. TECHNOLOGY TRANSFER DEMONSTRATION - COUPLING WITH NREL WISDEM Researchers have integrated the inverse-design tool for 2D airfoils (INN-Airfoil) into WISDEM (Wind Plant Integrated Systems Design and Engineering Model), a multidisciplinary design and optimization framework for assessing the cost of energy, as part of tech-transfer demonstration. The integration of INN-Airfoil into WISDEM allows for the design of airfoils along with the blades that meet the dynamic design constraints on cost of energy, annual energy production, and the capital costs. Through preliminary studies, researchers have shown that the coupled INN-Airfoil + WISDEM approach reduces the cost of energy by around 1% compared to the conventional design approach. This page will serve as a place to easily access all the publications from this work and the repositories for the software developed and released through this pr...

aerodynamics↗

Advances in relaxation and memory effects of magnetic nanoparticles for biomedical applications

Functionalized magnetic nanoparticles are pivotal in magnetic resonance imaging, computed tomography, controlled drug delivery, and hyperthermia treatments due to their exceptional magnetic relaxation and functional properties. The magnetic core composition and structure significantly affects the complex magnetic properties of these nanoparticles necessitating a thorough examination of magnetism fundamentals related to these systems. One important aspect is the ability of magnetic nanoparticles to retain previous magnetic state configurations known as memory effect, primarily governed by domain structure and magnetic anisotropy. Despite its relevance to advanced applications, comprehensive studies on magnetic relaxation and memory effects remain limited. Here, the present review aims to bridge this gap by investigating relaxation mechanisms, synthesis strategies, and applications, fostering further innovation. It investigates the memory effects and their dependence on particle composition and morphology along with key synthesis techniques for large-scale production in industrial adoption. Structured into focused sections on magnetic properties and their influence on biomedical and technological applications, this review provides essential insights into memory effects, magneto-relaxation mechanisms, influencing factors, and both experimental and theoretical methodologies. It also delves into computational modelling and AI-driven design, which are revolutionizing the prediction, discovery, and optimization of materials with tailored properties.

36 MATERIALS SCIENCE↗

Deciphering Reaction Mechanisms of Molecular Proton Reduction Catalysts with Cyclic Voltammetry: Kinetic vs Thermodynamic Control

The kinetics and thermodynamics of elementary reaction steps involved in the catalytic reduction of protons to hydrogen define the reaction landscape for catalysis. The mechanisms can differ in the order of the elementary proton transfer, electron transfer, and bond-forming steps and can be further differentiated by the sites at which protons and electrons localize. Access to fully elucidated mechanistic, kinetic, and thermochemical details of molecular catalysts is crucial to facilitate the development of new catalysts that operate with optimal efficiency, selectivity, and durability. The mechanism by which a catalyst operates, as well as the kinetics and thermodynamics associated with the individual steps, can often be accessed through electroanalytical studies. Here, this Account details the application of cyclic voltammetry to interrogate reaction mechanisms and quantify the kinetics and thermodynamics of elementary reaction steps for a series of molecular catalysts that mediate electrochemical proton reduction. I distinguish the limiting scenarios wherein a catalyst operates under kinetic control vs thermodynamic control, with a focus on detecting how cyclic voltammetry features shift with proton source strength and concentration, as well as scan rate. For systems that operate under kinetic control, catalytic currents are observed at, or slightly positive toward, the formal potential for the redox process that triggers catalysis. Under thermodynamic control, catalytic responses shift as a function of the proton source pKa and effective pH of the solution. After drawing this distinction, we introduce the appropriate voltammetry experiments and accompanying analytical expressions for extracting key metrics from the data. To illustrate analytical strategies to quantify elementary reaction steps of catalysts operating under kinetic control, I describe our studies of proton reduction catalysts Co(dmgBF 2 ) 2 (CH 3 CN) 2 (dmgBF 2 = difluoroboryl-dimethylglyoxime) and [Ni(P 2 Ph N 2 Ph ) 2 ] 2+ (P 2 Ph N 2 Ph = 1,5-phenyl-3,7-phenyl-1,5-diaza-3,7-diphosphacyclooctane). Here, peak shift analysis, foot-of-the-wave analysis, and plateau current analysis are applied to data sets wherein voltammetric response are recorded as a function of catalyst concentration, proton source concentration, proton source strength, and scan rate to quantify rate constants for elementary proton transfer and bond-forming steps in a catalytic cycle. Further, the case study of [Ni(P 2 Ph N 2 Ph ) 2 ] 2+ illustrates how complementary spectroscopic methods can bolster the mechanistic assignment. Collectively, these two studies showcase how detailed mechanistic studies inform on rate-limiting elementary steps in catalysis and other key processes underpinning catalysis. Second, I present analytical strategies to interrogate catalysts operating under thermodynamic control, centered on the case study of [Ni II (P 2 Ph N 2 Bn ) 2 ] 2+ (P 2 Ph N 2 Bn = 1,5-dibenzyl-3,7-diphenyl-1,5-diaza-3,7-diphosphacyclooctane). Here, the application of nonaqueous Pourbaix theory to extract thermodynamic information is introduced, and the construction of a coupled Pourbaix diagram is detailed. This study identifies ligand-based protonation as the key process that places catalysis under thermodynamic control and influences the reaction mechanism. Together, the work detailed in this Account showcases the utility of electroanalytical methods to disentangle complex reaction mechanisms and extract key thermochemical and kinetic parameters for elementary steps of catalysis. Through detailed presentation of the key analytical expressions that underpin these analyses, this Account seeks to facilitate the adoption of cyclic voltammetry by the community to fully extract kinetic, thermochemical, and mechanistic information on electrochemical small-molecule activation.

catalysts↗

Closed-Loop Control of Active Nematic Flows

Stabilizing and shaping autonomous flows of active fluids is a fundamental challenge and a prerequisite for applications. We embed a light-responsive microtubule-based nematic in a proportional-integral control loop that adjusts the applied light intensity in response to real-time measurements of the spatially averaged flow speed. The self-regulating hardware-software-wetware system maintains a target flow speed against external or internal perturbations, including protein aging and aggregation, sample-to-sample variability, and temperature variation. Varying the controller’s gains reveals antagonistic roles between feedback and intrinsic processes, leading to nontrivial dynamics observed in fluctuation spectra. In particular, oscillations emerge from the interplay between the controller, motor binding kinetics, and active hydrodynamic relaxation. Accounting for the underlying binding timescale, our coarse-grained model and nematohydrodynamics simulations corroborate these observations. This work provides insight into the coupled dynamics of controlled active matter, laying the foundation for spatiotemporal patterning of active stress to generate and stabilize new dynamical configurations.

Active nematics↗

Uncovering the mechanism of selective stabilization of high-energy diastereoisomers via inclusion

Abstract Supramolecular systems may be used to stabilize otherwise unstable isomers to find alternative synthetic pathways. It has been reported that cucurbit[8]uril can stabilize trans -I and trans -II Cu $$^{{\textrm{II}}}$$ II cyclam, whereas trans -III is the only non-substituted trans Cu $$^{{\textrm{II}}}$$ II cyclam diastereoisomer found outside of the host molecule experimentally. Quantum chemistry methods can provide valuable insight into the intermolecular interactions involved in these inclusion complexes. All five possible trans diastereoisomers of Cu $$^{{\textrm{II}}}$$ II cyclam were studied within the host molecule to calculate the interaction energy and free energy of association for each complex. The relative free energies of the five free cyclams confirm that trans -I and trans -II are the most energetically accessible diastereoisomers from the initial trans -III starting point. Energy decomposition analysis was used to identify the attractive and repulsive interactions between cyclam and cucurbit[8]uril and showed that trans -II encounters repulsive forces almost three times greater than trans -I, which may explain the 7:3 ratio of trans -I to trans -II within cucurbit[8]uril that occurs experimentally. Optimized complex geometries with trans -III, IV, and V show that the cyclams protrude out of cucurbit[8]uril, whereas trans -I and trans -II become more encapsulated and elongate the host, suggesting that the position of the cyclam is extremely important when forming non-covalent interactions. Our results agree with the experimental findings and provide greater insight into why the most stable isolated cyclam diastereoisomer, trans -III, does not form a complex.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evaluating Implementations of the Immersed Boundary Method in the Weather Research and Forecasting Model

The terrain-following coordinate system used by many atmospheric models can cause numerical instabilities due to discretization errors as resolved terrain slopes increase and the grid becomes highly skewed. The immersed boundary (IB) method, which does not require the grid to conform to the terrain, has been shown to alleviate these errors, and has been used successfully for high-resolution atmospheric simulations over steep terrain, including vertical building surfaces. Since many previous applications of IB methods to atmospheric models have used very fine grid resolution (5 m or less), the present study seeks to evaluate IB method performance over a range of grid resolutions and aspect ratios. Two classes of IB algorithms, velocity reconstruction and shear stress reconstruction, are tested within the common framework of the Weather Research and Forecasting (WRF) Model. Performance is evaluated in two test cases, one with flat terrain and the other with the topography of Askervein Hill, both under neutrally stratified conditions. WRF-IB results are compared to similarity theory, observations, and native WRF results. Despite sensitivity to the location at which the IB intersects the model grid, the velocity reconstruction IB method shows consistent performance when used with a hybrid RANS/LES surface scheme. The shear stress reconstruction IB method is not sensitive to the grid intersection, but is less consistent and near-surface velocity errors can occur at coarse resolutions. This study represents an initial investigation of IB method variability across grid resolutions in WRF. Future work will focus on improving IB method performance at intermediate to coarse resolutions.

54 ENVIRONMENTAL SCIENCES↗

Spin-state gaps and self-interaction-corrected density functional approximations: Octahedral Fe(II) complexes as case study

Accurate prediction of a spin-state energy difference is crucial for understanding the spin crossover phenomena and is very challenging for density functional approximations, especially for local and semi-local approximations due to delocalization errors. Here, we investigate the effect of the self-interaction error removal from the local spin density approximation (LSDA) and Perdew–Burke–Ernzerhof generalized gradient approximation on the spin-state gaps of Fe(II) complexes with various ligands using recently developed locally scaled self-interaction correction (LSIC) by Zope et al. [J. Chem. Phys. 151, 214108 (2019)]. The LSIC method is exact for one-electron density, recovers the uniform electron gas limit of the underlying functional, and approaches the well-known Perdew–Zunger self-interaction correction (PZSIC) as a particular case when the scaling factor is set to unity. Our results, when compared with reference diffusion Monte Carlo results, show that the PZSIC method significantly overestimates spin-state gaps favoring low spin states for all ligands and does not improve upon density functional approximations. The perturbative LSIC-LSDA using PZSIC densities significantly improves the gaps with a mean absolute error of 0.51 eV but slightly overcorrects for the stronger CO ligands. Finally, the quasi-self-consistent LSIC-LSDA, such as coupled-cluster single double and perturbative triple [CCSD(T)], gives a correct sign of spin-state gaps for all ligands with a mean absolute error of 0.56 eV, comparable to that of CCSD(T) (0.49 eV).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Lieb-Robinson Bounds with Exponential-in-Volume Tails

Lieb-Robinson bounds demonstrate the emergence of locality in many-body quantum systems. Intuitively, Lieb-Robinson bounds state that, with local or exponentially decaying interactions, the correlation that can be built up between two sites separated by distance 𝑟 after a time 𝑡 decays as exp (𝑣⁢𝑡 −𝑟), where 𝑣 is the emergent Lieb-Robinson velocity. In many problems, it is important to also capture how much of an operator grows to act on 𝑟 𝑑 sites in 𝑑 spatial dimensions. Perturbation theory and cluster expansion methods suggest that, at short times, these volume-filling operators are suppressed as exp (−𝑟 𝑑 ). We confirm this intuition, showing that, for 𝑟 >𝑣⁢𝑡, the volume-filling operator is suppressed by exp (−(𝑟−𝑣⁢𝑡) 𝑑 /(𝑣⁢𝑡) 𝑑−1 ). This closes a conceptual and practical gap between the cluster expansion and the Lieb-Robinson bound. We then present two very different applications of this new bound. Firstly, we obtain improved bounds on the classical computational resources necessary to simulate many-body dynamics with error tolerance 𝜀 for any finite time 𝑡: as 𝜀 becomes sufficiently small, only 𝜀 −O⁡(𝑡 𝑑−1 ) resources are needed. A protocol that likely saturates this bound is given. Secondly, we prove that disorder operators have volume-law suppression near the “solvable (Ising) point” in quantum phases with spontaneous symmetry breaking, which implies a new diagnostic for distinguishing many-body phases of quantum matter.

computational complexity↗

Connection between the physicochemical characteristics of amorphous carbon thin films and their electrochemical properties

Connecting a material's surface chemistry with its electrocatalytic performance is one of the major questions in analytical electrochemistry. This is especially important in many sensor applications where analytes from complex media need to be measured. Unfortunately, today this connection is still largely missing except perhaps for the most simple ideal model systems. Here we present an approach that can be used to obtain insights about this missing connection and apply it to the case of carbon nanomaterials. In this paper we show that by combining advanced computational techniques augmented by machine learning methods with x-ray absorption spectroscopy (XAS) and electrochemical measurements, it is possible to obtain a deeper understanding of the correlation between local surface chemistry and electrochemical performance. As a test case we show how by computationally assessing the growth of amorphous carbon (a-C) thin films at the atomic level, we can create computational structural motifs that may in turn be used to deconvolute the XAS data from the real samples resulting in local chemical information. Then, by carrying out electrochemical measurements on the same samples from which x-ray spectra were measured and that were further characterized computationally, it is possible to gain insight into the interplay between the local surface chemistry and electrochemical performance. To demonstrate this methodology, we proceed as follows: after assessing the basic electrochemical properties of a-C films, we investigate the effect of short HNO 3 treatment on the sensitivity of these electrodes towards an inner sphere redox probe dopamine to gain knowledge about the influence of altered surface chemistry to observed electrochemical performance. Overall, these results pave the way towards a more general assessment of electrocatalysis in different systems and provide the first steps towards data driven tailoring of electrode surfaces to gain optimal performance in a given application.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗