DFT Virtual Screening Identifies Rhodium-Amidinate Complexes as Potential Homogeneous Catalysts for Methane-to-Methanol Oxidation
Explore the source record for details and available documents.
catalysis (homogeneous), catalysis (heterogeneous)↗
SEARCH · Search NASA
Results for “explicit correlation”
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.
Instability and dynamics of two nonlinearly coupled intense laser beams in a quantum plasma
We consider nonlinear interactions between two relativistically strong laser beams and a quantum plasma composed of degenerate electron fluids and immobile ions. The collective behavior of degenerate electrons is modeled by quantum hydrodynamic equations composed of the electron continuity, quantum electron momentum (QEM) equation, as well as the Poisson and Maxwell equations. The QEM equation accounts the quantum statistical electron pressure, the quantum electron recoil due to electron tunneling through the quantum Bohm potential, electron-exchange, and electron-correlation effects caused by electron spin, and relativistic ponderomotive forces (RPFs) of two circularly polarized electromagnetic (CPEM) beams. The dynamics of the latter are governed by nonlinear wave equations that include nonlinear currents arising from the relativistic electron mass increase in the CPEM wave fields, as well as from the beating of the electron quiver velocity and electron density variations reinforced by the RPFs of the two CPEM waves. Furthermore, nonlinear electron density variations associated with the driven (by the RPFs) quantum electron plasma oscillations obey a coupled nonlinear Schroedinger and Poisson equations. The nonlinearly coupled equations for our purposes are then used to obtain a general dispersion relation (GDR) for studying the parametric instabilities and the localization of CPEM wave packets in a quantum plasma. Numerical analyses of the GDR reveal that the growth rate of a fastest growing parametrically unstable mode is in agreement with the result that has been deduced from numerical simulations of the governing nonlinear equations. Explicit numerical results for two-dimensional (2D) localized CPEM wave packets at nanoscales are also presented. Possible applications of our investigation to intense laser-solid density compressed plasma experiments are highlighted.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Reducing the Cost of Energy Differences in Variational Monte Carlo with Spotlight Sampling
Here, we investigate an approximate sampling scheme that can significantly reduce the cost scaling of variational Monte Carlo when it is employed to predict the energy differences associated with local chemical changes. Inspired by side-chaining and embedding methods, this spotlight sampling approach adopts an approximate fragmented Hamiltonian and correlated sampling to reduce cost scaling to the point that it is essentially linear with system size, with the potential to go sublinear if certain conditions are met. In tests on bond stretching energies in alcohols, hydrogen dimer chains, and molecules with various degrees of π-system delocalization, we observe the anticipated linear scaling and an explicit cost crossover with standard variational Monte Carlo.
Bumann, Sonja [University of California, Berkeley,↗
On the non-linear scale of cosmological perturbation theory
We discuss the convergence of cosmological perturbation theory. We prove that the polynomial enhancement of the non-linear corrections expected from the effects of soft modes is absent in equal-time correlators like the power or bispectrum. We first show this at leading order by resumming the most important corrections of soft modes to an arbitrary skeleton of hard fluctuations. We derive the same result in the eikonal approximation, which also allows us to show the absence of enhancement at any order. We complement the proof by an explicit calculation of the power spectrum at two-loop order, and by further numerical checks at higher orders. Using these insights, we argue that the modification of the power spectrum from soft modes corresponds at most to logarithmic corrections at any order in perturbation theory. Finally, we discuss the asymptotic behavior in the large and small momentum regimes and identify the expansion parameter pertinent to non-linear corrections.
79 ASTRONOMY AND ASTROPHYSICS↗
Hamming distance and the onset of quantum criticality
Simulating models for quantum correlated matter unveils the inherent limitations of deterministic classical computations. In particular, in the case of quantum Monte Carlo methods, this is manifested by the emergence of negative weight configurations in the sampling, that is, the sign problem (SP). There have been several recent calculations which exploit the SP to locate underlying critical behavior. Here, utilizing a metric that quantifies phase-space ergodicity in such sampling, the Hamming distance, we suggest a significant advance on these ideas to extract the location of quantum critical points in various fermionic models, in spite of the presence of a severe SP. Combined with other methods, exact diagonalization in our case, it elucidates both the nature of the different phases as well as their location, as we demonstrate explicitly for the honeycomb and triangular Hubbard models, in both their U(1) and SU(2) forms. Furthermore, our approach exemplifies a possible path allowing the exploration of the phase diagram of a variety of fermionic quantum models hitherto considered to be impractical via quantum Monte Carlo simulations.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Mean field model parameterization to recover coverage-dependent kinetics
Lateral interactions between adsorbates introduce coverage dependence into adsorption energies and activation barriers of surface reactions. Lattice-based kinetic Monte Carlo (kMC) simulations can capture these interactions quantitatively but are laborious to parameterize and solve. Mean field models are more tractable, but protocols to construct and parameterize them are unclear. Here we explore the ability of a coverage-aware mean-field model to map to a lattice-kMC model of a generic two-step reaction network, including quasi-equilibrated adsorption and rate-limiting dissociation steps. We derive expressions for mean-field and coverage-dependent adsorption energies and dissociation barriers and parameterize against lattice-kMC predictions. We show that the parameterized mean-field rates correlate with ground truth lattice-kMC results across a wide range of reaction conditions and identify regions where the mean field fails. The mean field model similarly captures kMC-derived rate-order, Arrhenius and Sabatier plots at a greatly reduced computational cost. Further, the results provide guidance for parameterizing mean-field models, benchmarked against explicit lattice-based approaches for incorporating the influence of coverage effects.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Communication: A new class of non-empirical explicit density functionals on the third rung of Jacob’s ladder
We construct an orbital-free non-empirical meta-generalized gradient approximation (GGA) functional, which depends explicitly on density through the density overlap regions indicator [P. de Silva and C. Corminboeuf, J. Chem. Theory Comput. 10, 3745 (2014)]. The functional does not depend on either the kinetic energy density or the density Laplacian; therefore, it opens a new class of meta-GGA functionals. By construction, our meta-GGA yields exact exchange and correlation energy for the hydrogen atom and recovers the second order gradient expansion for exchange in the slowly varying limit. We show that for molecular systems, overall performance is better than non-empirical GGAs. For atomization energies, performance is on par with revTPSS, without any dependence on Kohn-Sham orbitals.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
SODAs: sparse optimization for the discovery of differential and algebraic equations
Differential-algebraic equations (DAEs) integrate ordinary differential equations (ODEs) with algebraic constraints, providing a fundamental framework for developing models of dynamical systems characterized by time-scale separation, conservation laws and physical constraints. While sparse optimization has revolutionized model development by allowing data-driven discovery of parsimonious models from a library of possible equations, existing approaches for dynamical systems assume DAEs can be reduced to ODEs by eliminating variables before model discovery. This assumption limits the applicability of such methods for DAE systems with unknown constraints and time scales. We introduce sparse optimization for differential-algebraic systems (SODAs), a data-driven method for the identification of DAEs in their explicit form. By discovering the algebraic and dynamic components sequentially without prior identification of the algebraic variables, this approach leads to a sequence of convex optimization problems. It has the advantage of discovering interpretable models that preserve the structure of the underlying physical system. To this end, SODAs improves since SODAs is singular numerical stability when handling high correlations between library terms, caused by near-perfect algebraic relationships, by iteratively refining the conditioning of the candidate library. We demonstrate the performance of our method on biological, mechanical and electrical systems, showcasing its robustness to noise in both simulated time series and real-time experimental data.
DAE↗
Stochastic approach to correlations beyond the mean field with the Skyrme interaction
Large-scale calculation based on the multi-configuration Skyrme density functional theory is performed for the light N=Z even-even nucleus, {sup 12}C. Stochastic procedures and the imaginary-time evolution are utilized to prepare many Slater determinants. Each state is projected on eigenstates of parity and angular momentum. Then, performing the configuration mixing calculation with the Skyrme Hamiltonian, we obtain low-lying energy-eigenstates and their explicit wave functions. The generated wave functions are completely free from any assumption and symmetry restriction. Excitation spectra and transition probabilities are well reproduced, not only for the ground-state band, but for negative-parity excited states and the Hoyle state.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Shock compression behavior of stainless steel 316L octet-truss lattice structures
Lattice structures offer desirable mechanical properties for applications of energy absorption and impact mitigation but limited research has been carried out on their shock compression behavior. In this work, the shock compression behavior of stainless steel 316L (SS316L) octet-truss lattice structures was investigated through experimental techniques and numerical simulations. Plate impact experiments with high-speed imaging were conducted at impact velocities of 270 – 390 m/s on lattice specimens with 5x5x10 unit cell geometries additively manufactured (AM) using direct metal laser sintering. High-speed imaging together with digital image correlation was used to extract full-field measurements and define a two-wave structure consisting of an elastic wave and planar compaction (shock) wave which propagated along the impact direction. A linear shock velocity versus particle velocity relation was found to approximate the measurements with a unit slope and a linear fit constant equal to the crushing speed. Furthermore, the shock velocity versus particle velocity relation, full-field measurements, and elastic limit together with the Eulerian form of the Rankine-Hugoniot jump conditions were used to find relations for the stress and internal energy behind the shock. Stress behind the shock increased with relative density and particle velocity, and specific internal energy converged to a single curve similar to that of bulk AM SS316L. Explicit finite element analysis using the Johnson-Cook constitutive model demonstrated similar shock behavior observed in experiments and a linear shock velocity versus particle velocity relation and corresponding Hugoniot calculations were found to be in agreement with experimental results. Numerical simulations confirmed negligible effects of exterior versus interior measurements and further validated the application of one-dimensional shock theory.
42 ENGINEERING↗
A modal wave-packet model for the multi-mode Richtmyer–Meshkov instability
In this work, a model for multimode perturbations subject to the Richtmyer–Meshkov (RM) instability is presented and compared with simulations and experiments for conditions relevant to inertial confinement fusion. The model utilizes the single mode response to the RM impulse whereby its amplitude h(k, t) first grows with an initial velocity V 0 ∝ kh(k, 0) that eventually decays in time as 1/kV 0 t. Both the growth and saturation stages are subject to nonlinearities since they depend explicitly on the initial amplitude. However, rather than using the individual mode amplitude h(k, t), nonlinearity is taken to occur when the root-mean-square amplitude h rms (k, t) of a wave-packet within wavenumbers k ± δk becomes comparable to 1/k. This is done because nearby sidebands can act in unison for an auto-correlation distance 1/δk beyond nonlinearity as observed in the beam-plasma instability. Thus, the nonlinear saturation amplitude for each mode is reduced from the usual 1/k by a phase space factor that depends on the physical dimensionality, as in the Haan model for the Rayleigh–Taylor instability. In addition, for RM, the average value of kh rms for the initial spectrum is used to calculate a nonlinear factor F NL that reduces V 0 , as observed for single modes. For broadband perturbations, the model describes self-similar growth ∝t θ as successively longer wavelength modes reach saturation. The growing and saturated modes must be discerned because only the former promote θ and are enhanced by reshock and spherical convergence. All of these flows are described here by the model in good agreement with simulations and experiments.
42 ENGINEERING↗
Modification of the Jet Energy-Energy Correlator in Cold Nuclear Matter
We compute medium corrections to the energy-energy correlator (EEC) for jets in electron-nucleus collisions at leading order in the QCD coupling and the interaction of the jet with the medium. We derive an analytical expression for the modification of the EEC as a function of the opening angle and show that the modification is strongest at large angles within the jet cone. We obtain explicit results for the dependence of the modification on the jet energy, the scattering power of cold nuclear matter, and the path length within the medium. We extend our calculations to gluon jets in proton-nucleus collisions and compare our results with recent preliminary data for proton-lead collisions at the LHC. Here, we also discuss the role of comovers on the EEC in p + Pb collisions.
Hard scattering↗
Stochastic generation of explicit pore structures by thresholding Gaussian random fields
We provide a description and computational investigation of an efficient method to stochastically generate realistic pore structures. Smolarkiewicz and Winter introduced this specific method in pores resolving simulation of Darcy flows (Smolarkiewicz and Winter, 2010 [1]) without giving a complete formal description or analysis of the method, or indicating how to control the parameterization of the ensemble. We address both issues in this paper. The method consists of two steps. First, a realization of a correlated Gaussian field, or topography, is produced by convolving a prescribed kernel with an initial field of independent, identically distributed random variables. The intrinsic length scales of the kernel determine the correlation structure of the topography. Next, a sample pore space is generated by applying a level threshold to the Gaussian field realization: points are assigned to the void phase or the solid phase depending on whether the topography over them is above or below the threshold. Hence, the topology and geometry of the pore space depend on the form of the kernel and the level threshold. Manipulating these two user prescribed quantities allows good control of pore space observables, in particular the Minkowski functionals. Extensions of the method to generate media with multiple pore structures and preferential flow directions are also discussed. To demonstrate its usefulness, the method is used to generate a pore space with physical and hydrological properties similar to a sample of Berea sandstone. -- Graphical abstract: -- Highlights: •An efficient method to stochastically generate realistic pore structures is provided. •Samples are generated by applying a level threshold to a Gaussian field realization. •Two user prescribed quantities determine the topology and geometry of the pore space. •Multiple pore structures and preferential flow directions can be produced. •A pore space based on Berea sandstone is generated.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Spatial and temporal agreement in climate model simulations of the Interdecadal Pacific Oscillation
Accelerated warming and hiatus periods in the long-term rise of Global Mean Surface Temperature (GMST) have, in recent decades, been associated with the Interdecadal Pacific Oscillation (IPO). Critically, decadal climate prediction relies on the skill of state-of-the-art climate models to reliably represent these low-frequency climate variations. We undertake a systematic evaluation of the simulation of the IPO in the suite of Coupled Model Intercomparison Project 5 (CMIP5) models. We track the IPO in pre-industrial (control) and all-forcings (historical) experiments using the IPO tripole index (TPI). The TPI is explicitly aligned with the observed spatial pattern of the IPO, and circumvents assumptions about the nature of global warming. We find that many models underestimate the ratio of decadal-to-total variance in sea surface temperatures (SSTs). However, the basin-wide spatial pattern of positive and negative phases of the IPO are simulated reasonably well, with spatial pattern correlation coefficients between observations and models spanning the range 0.4–0.8. Deficiencies are mainly in the extratropical Pacific. Models that better capture the spatial pattern of the IPO also tend to more realistically simulate the ratio of decadal to total variance. Of the 13% of model centuries that have a fractional bias in the decadal-to-total TPI variance of 0.2 or less, 84% also have a spatial pattern correlation coefficient with the observed pattern exceeding 0.5. This result is highly consistent across both IPO positive and negative phases. This is evidence that the IPO is related to one or more inherent dynamical mechanisms of the climate system.
54 ENVIRONMENTAL SCIENCES↗
Six-point functions and collisions in the black hole interior
In the eternal AdS black hole geometry, we consider two signals sent from the boundaries into the black hole interior shared between the two asymptotic regions. We compute three different out-of-time-order six-point functions to quantify various properties of the collision of these signals behind the horizons: (i) We diagnose the strength of the collision by probing the two-signal state on a late time slice with boundary operators. (ii) We quantify two-sided operator growth, which provides a dual description of the signals meeting in the black hole interior, in terms of the quantum butterfly effect and quantum circuits. (iii) We consider an explicit coupling between the left and right CFTs to make the wormhole traversable and extract information about the collision product from behind the horizon. At a technical level, our results rely on the method of eikonal resummation to obtain the relevant gravitational contributions to Lorentzian six-point functions at all orders in the GN-expansion. We observe that such correlation functions display an intriguing factorization property. We corroborate these results with geodesic computations of six-point functions in two- and three-dimensional gravity.
2D Gravity↗
Kβ X-ray Emission Spectroscopy of Cu(I)-Lytic Polysaccharide Monooxygenase: Direct Observation of the Frontier Molecular Orbital for H 2 O 2 Activation
Lytic polysaccharide monooxygenases (LPMOs) catalyze the degradation of recalcitrant carbohydrate polysaccharide substrates. These enzymes are characterized by a mononuclear Cu(I) active site with a three-coordinate T-shaped “His-brace” configuration including the N-terminal histidine and its amine group as ligands. Here, this study explicitly investigates the electronic structure of the d 10 Cu(I) active site in a LPMO using Kβ X-ray emission spectroscopy (XES). The lack of inversion symmetry in the His-brace site enables the 3d/p mixing required for intensity in the Kβ valence-to-core (VtC) XES spectrum of Cu(I)-LPMO. These Kβ XES data are correlated to density functional theory (DFT) calculations to define the bonding, and in particular, the frontier molecular orbital (FMO) of the Cu(I) site. These experimentally validated DFT calculations are used to evaluate the reaction coordinate for homolytic cleavage of the H 2 O 2 O–O bond and understand the contribution of this FMO to the low barrier of this reaction and how the geometric and electronic structure of the Cu(I)-LPMO site is activated for rapid reactivity with H 2 O 2 .
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
A practical computational method for the anisotropic redshift-space three-point correlation function
Here, we present an algorithm enabling computation of the anisotropic redshift-space galaxy threepoint correlation function ( 3 PCF) scaling as N 2 , with N the number of galaxies. Our previous work showed how to compute the isotropic 3 PCF with this scaling by expanding the radially binned density field around each galaxy in the survey into spherical harmonics and combining these coefficients to form multipole moments. The N 2 scaling occurred because this approach never explicitly required the relative angle between a galaxy pair about the primary galaxy. Here, we generalize this work, demonstrating that in the presence of azimuthally symmetric anisotropy produced by redshift-space distortions (RSD), the 3 PCF can be described by two triangle side lengths, two independent total angular momenta, and a spin. This basis for the anisotropic 3 PCF allows its computation with negligible additional work over the isotropic 3 PCF. We also present the covariance matrix of the anisotropic 3 PCF measured in this basis. Our algorithm tracks the full 5 D redshift-space 3 PCF, uses an accurate line of sight to each triplet, is exact in angle, and easily handles edge correction. It will enable use of the anisotropic large-scale 3 PCF as a probe of RSD in current and upcoming large-scale redshift surveys.
79 ASTRONOMY AND ASTROPHYSICS↗
LA100 Equity Strategies. Chapter 12: Distribution Grid Upgrades for Equitable Resilience and Solar, Storage, and Electric Vehicle Access
The LA100 Equity Strategies project integrates community guidance with robust research, modeling, and analysis to identify strategy options that can increase equitable outcomes in Los Angeles' clean energy transition. As Los Angeles transitions toward clean energy, existing distribution grid infrastructure will need to be updated and expanded to support reliable service during routine operations, enable interconnection with distributed energy resources and electrified loads, and provide access to energy-related services during disasters. This chapter focuses on equity in distribution grid upgrades, reliability, and resilience in Los Angeles. Specifically, NREL performed grid upgrade and resilience analyses using a detailed model of the distribution grid and income-differentiated household load profiles, electric vehicle (EV) adoption patterns, distributed solar adoption, and grid reliability to explore two key questions to inform how the City of Los Angeles can ensure a resilient and reliable distribution grid for all communities during the clean energy transition: Where can distribution system upgrades can be prioritized to enable equitable access to, and adoption of, clean energy technologies and how can Los Angeles provide equitable, resilient access to electricity-related services (e.g., health care, food) during disaster events like earthquakes and flooding? The electric distribution system is the "last mile" of the grid, linking the multistate bulk power system with customers; new loads, including EVs; and distributed energy resources, such as customer and community solar and storage. This analysis focuses on the 4.8-kilovolt (kV) system, including service transformers that represent the utility-side of the grid connection for most residential customers. Chapter 17 looks at the customer-side of the grid connection with a focus on electric panel upgrade needs. The transition toward clean energy can put additional stress on the distribution system from distributed energy resources and electrification - especially EVs and increased use of electricity for heating, cooling, cooking, and hot water. This stress, measured here as the number of equipment overloads and voltage violations, correlates strongly to grid reliability and therefore is used as a proxy for understanding additional upgrades needed and to help ensure equitable access to electrification and distributed energy resources. NREL also conducted community resilience analysis to examine customer-level access to both electricity and a larger range of services, such as hospitals and grocery stores during a disaster. This analysis explicitly considers equity to understand differences in current resilience and resilience strategies to effectively improve critical services access for all Angelenos. Research was guided by input from the community engagement process, and associated equity strategies are presented in alignment with that guidance.
14 SOLAR ENERGY↗