Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density 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.

At least 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Thermodynamics–Structure–Dynamics Correlations and Nonuniversal Effects in the Elastically Collective Activated Hopping Theory of Glass-Forming Liquids

We employ the microscopic Elastically Collective Nonlinear Langevin Equation (ECNLE) theory of activated dynamics in combination with crystal-avoiding simulations to study four inter-related questions for metastable monodisperse hard sphere fluids. The first is how significantly improved integral equation theory structural input (Modified-Verlet (MV) closure) changes the dynamical predictions of ECNLE theory. The main consequence is a modest enhancement of the importance of the collective elastic barrier relative to its local cage contribution, which increases the alpha relaxation time and fragility relative to prior results based on the Percus–Yevick closure. Second, ECNLE-MV theory predictions for the alpha time and self-diffusion constant in the metastable regime are quantitatively compared to our new simulations. The small adjustment of a numerical prefactor that enters the collective elastic barrier leads to quantitative agreement over three decades. Third, using the more accurate MV structural input, ECNLE theory is shown to predict thermodynamics–structure–dynamics “correlations” based on various long and short wavelength scalar properties all related to static two-point collective density fluctuations. The logarithm of the alpha relaxation time scales as a power law with these scalar metrics with an exponent that is significantly lower in the less dense noncooperative activated regime compared to the very dense highly cooperative regime. However, the discovered correlation of activated relaxation with a thermodynamic property (dimensionless compressibility) is not causal in ECNLE theory, but rather reflects a strong connection between the local structural quantities that quantify kinetic constraints in the theory with the amplitude of long wavelength density fluctuations. Fourth, the consequences of chemically specific nonuniversalities associated with the onset condition and relative importance of collective elasticity are studied. In conclusion, the predicted thermodynamics–structure–dynamics correlations are found to be robust, albeit with nontrivial shifts of the onset condition.

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

Evidence of Ba-substitution induced spin-canting in the magnetic Weyl semimetal EuCd 2 As 2

Recently EuCd 2 As 2 was predicted to be a magnetic Weyl semimetal with a lone pair of Weyl nodes generated by A-type antiferromagnetism and protected by a rotational symmetry. However, it was soon discovered that the actual magnetic structure broke the rotational symmetry and internal pressure was later suggested as a route to stabilize the desired magnetic state. In this work we test this prediction by synthesizing a series of Eu 1-x BaxCd 2 As 2 single crystals and studying their structural, magnetic, and transport properties via both experimental techniques and first-principles calculations. We find that small concentrations of Ba (~3%–10%) lead to a small out-of-plane canting of the Eu moment. However, for higher concentrations this effect is suppressed and a nearly in-plane model is recovered. Studying the transport properties we find that all compositions show evidence of an anomalous Hall effect dominated by the intrinsic mechanism as well as large negative magnetoresistances in the longitudinal channel. A nonmonotonic evolution of the transport properties is seen across the series which correlates to the proposed canting suggesting canting may enhance the topological effects. Careful density functional theory calculations using an all-electron approach revise prior predictions finding a purely ferromagnetic ground state with in-plane moments for both the EuCd 2 As 2 and Eu 0.5 Ba 0.5 Cd 2 As 2 compounds, corroborating our experimental findings. This work suggests that Ba substitution can tune the magnetic properties in unexpected ways which correlate to changes in measures of topological properties, encouraging future work to locate the ideal Ba concentration for Eu moment canting.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

False metals, real insulators, and degenerate gapped metals

This paper deals with a significant family of compounds predicted by simplistic electronic structure theory to be metals but are, in fact, insulators. This false metallic state has been traditionally attributed in the literature to reflect the absence of proper treatment of electron-electron correlation (“Mott insulators”) whereas, in fact, even mean-field like density functional theory describes the insulating phase correctly if the restrictions posed on the simplistic theory are avoided. Such unwarranted restrictions included different forms of disallowing symmetry breaking described in this article. As the science and technology of conductors have transitioned from studying simple elemental metals such as Al or Cu to compound conductors such as binary or ternary oxides and pnictides, a special class of degenerate but gapped metals has been noticed. Their presumed electronic configurations show the Fermi level inside the conduction band or valence band, yet there is an “internal band gap” between the principal band edges. The significance of this electronic configuration is that it might be unstable toward the formation of states inside the internal band gap when the formation of such states costs less energy than the energy gained by transferring carriers from the conduction band to these lower energy acceptor states, changing the original (false) metal to an insulator. The analogous process also exists for degenerate but gapped metals with the Fermi level inside the valence band, where the energy gain is defined by transfer of electrons from the donor level to the unoccupied part of the valence band. We focus here on the fact that numerous electronic structure methodologies have overlooked some physical factors that could stabilize the insulating alternative, predicting instead false metals that do not really exist (note that this is in general not a physical phase transition, but a correction of a previous error in theory that led to a false prediction of a metal). Such errors include: (i) ignoring spin symmetry breaking, such as disallowing magnetic spin ordering in CuBi2O4 or disallowing the formation of polymorphous spin networks in paramagnetic LaTiO 3 and YTiO 3 ; (ii) ignoring structural symmetry breaking, e.g., not enabling energy-lowering bond disproportionation (Li-doped TiO 2 , SrBiO 3 , and rare-earth nickelates), or not exploring pseudo-Jahn–Teller-like distortions in LaMnO 3 , or disallowing spontaneous formation of ordered vacancy compounds in Ba 4 As 3 and Ag 3 Al 22 O 34 ; and (iii) ignoring spin–orbit coupling forcing false metallic states in CaIrO 3 and Sr 2 IrO 4 . The distinction between false metals vs real insulators is important because (a) predicting theoretically that a given compound is metal even though it is found to be an insulator often creates the temptation to invoke high order novel physical effects (such as correlation in d-electron Mott insulators) to explain what was in effect caused by a more mundane artifact in a lower-level mean-field band theory, (b) recent prediction of exotic physical effects such as topological semimetals were unfortunately based on the above compounds that were misconstrued by theory to be metal, but are now recognized to be stable insulators not hosting exotic effects, and (c) practical technological applications based on stable degenerate but gapped metals such as transparent conductors or electrides for catalysis must rely on the systematically correct and reliable theoretical classification of metals vs insulators.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Triangular lattice Hubbard model physics at intermediate temperatures

Moiré systems offer an exciting playground to study the many-body effects of strongly correlated electrons in regimes that are not easily accessible in conventional material settings. Motivated by a recent experiment on WSe 2 /WS 2 moiré bilayers [Y. Tang et al., Nature (London) 579, 353 (2020)], which realized a triangular superlattice with a small hopping (of approximately 10 Kelvin), with tunable density of holes, we explore the Hubbard model on the triangular lattice for intermediate temperatures t ≲ T < U. Employing finite-temperature Lanczos calculations and closely following the fitting protocols used in the experiment, we recover the observed trends in the reported Curie-Weiss temperature Θ with filling, using the reported interaction strength U/t=20. Here, we focus on the large increase of |Θ| on decreasing the density below half filling and the sign change of Θ at higher fillings, which signals the onset of ferromagnetism. The increase in |Θ| is also seen in the t–J model (the low-energy limit of the Hubbard model) in the intermediate temperature range, which we clarify is opposite to the trend in its high temperature limit. Differences between the low, intermediate, and high temperature regimes are discussed. Our numerical calculations also capture the crossover between short-range antiferromagnetic to ferromagnetic order in the intermediate temperature regime, a result broadly consistent with the experimental findings. We find that this behavior is a finite-temperature remnant of the underlying zero-temperature phase transition, which we explore with ground-state density matrix renormalization group calculations. We provide evidence of ferromagnetism characterized by weak (but robust) correlations that explain the small Θ seen in the experiment.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Thermoelectricity and electronic correlation enhancement in FeS by light Se doping

We report thermoelectric studies of FeS 1–x Se x (x = 0, 0.06) superconducting single crystals that feature high irreversibility fields and critical current density J c comparable to materials with much higher superconducting critical temperatures (T c 's). The ratio of T c to the Fermi temperature T F is very small, indicating weak electronic correlations. With a slight selenium substitution on sulfur site in FeS both T c /T F and the effective mass m* rise considerably, implying increase in electronic correlation of the bulk conducting states. The first-principle calculations show rise of the density of states at the Fermi level in FeS 0.94 Se 0.06 when compared to FeS, which is related not only to Fe but also to chalcogen-derived electronic states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Chiral effective field theory calculation of neutrino reactions in warm neutron-rich matter

Neutrino scattering and absorption rates of relevance to supernovae and neutron star mergers are obtained from nuclear matter dynamical structure functions that encode many-body effects from nuclear mean fields and correlations. In this study, we employ nuclear interactions from chiral effective field theory to calculate the density, spin, isospin, and spin-isospin response functions of warm β-equilibrium nuclear matter. We include corrections to the single-particle energies in the mean field approximation as well as vertex corrections resummed in the random phase approximation (RPA), including both direct and exchange diagrams. We find that correlations included through the RPA redistribute the strength of the response to higher energy for neutrino absorption and lower energy for antineutrino absorption. This tends to suppress the absorption rate of electron neutrinos across all relevant energy scales. In contrast, the inclusion of RPA correlations enhances the electron antineutrino absorption rate at low energy and suppresses the rate at high energy. These effects are especially important at high-density and in the vicinity of the neutrino decoupling region. Implications for heavy element nucleosynthesis, electromagnetic signatures of compact object mergers, supernova dynamics, and neutrino detection from galactic supernovae are discussed briefly.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Cross-correlation of SPT-3G D1 CMB lensing and DES Y3 galaxy lensing

Measurements of the weak lensing of galaxies and of the cosmic microwave background (CMB) provide direct probes of the cosmic matter density field, but the two observables are sensitive to different spatial scales, redshift ranges, and survey systematics. Their cross-correlation thus enables consistency checks of the theoretical model and of potential systematics in either dataset. We present measurements of the cross-correlation between CMB lensing and cosmic shear over $\sim$1,300 deg$^2$ of the sky using the SPT-3G D1 CMB lensing maps and the Dark Energy Survey Year 3 (DES Y3) shear catalogs. For the first time, we measure this cross-correlation at high significance ($\sim 14σ$) when using a polarization-only CMB lensing reconstruction that is expected to be robust against biases induced by extragalactic foregrounds. We test a variety of other CMB lensing estimators that include temperature information and exhibit different tradeoffs between foreground biases and noise, as well as a shear sample that consists of blue, star-forming galaxies and has been shown to be less impacted by galaxy intrinsic alignments. Assuming $Λ$CDM and marginalizing over uncertainties in intrinsic alignments, baryonic feedback, and various nuisance parameters, we obtain a constraint on the amplitude of matter clustering $S_8 \equiv σ_8 \sqrt{Ω_m / 0.3} = 0.833^{+0.047}_{-0.061}$, consistent with both the primary CMB results from Planck and shear-only results from DES Y3. By combining our measurement with Planck, we find mild constraints on the astrophysical processes that impact the cross-correlation. We obtain a constraint on the intrinsic alignment amplitude of the DES sample that is competitive with that from shear-only analyses, and we find a lower limit on the strength of baryonic feedback.

Ouellette, A. [Illinois U., Urbana] (ORCID:0000000↗
DOE OSTI2020

TurboRVB: A many-body toolkit for ab initio electronic simulations by quantum Monte Carlo

TurboRVB is a computational package for ab initio Quantum Monte Carlo (QMC) simulations of both molecular and bulk electronic systems. The code implements two types of well established QMC algorithms: Variational Monte Carlo (VMC) and diffusion Monte Carlo in its robust and efficient lattice regularized variant. A key feature of the code is the possibility of using strongly correlated many-body wave functions (WFs), capable of describing several materials with very high accuracy, even when standard mean-field approaches [e.g., density functional theory (DFT)] fail. The electronic WF is obtained by applying a Jastrow factor, which takes into account dynamical correlations, to the most general mean-field ground state, written either as an antisymmetrized geminal power with spin-singlet pairing or as a Pfaffian, including both singlet and triplet correlations. This WF can be viewed as an efficient implementation of the so-called resonating valence bond (RVB) Ansatz, first proposed by Pauling and Anderson in quantum chemistry [L. Pauling, The Nature of the Chemical Bond (Cornell University Press, 1960)] and condensed matter physics [P.W. Anderson, Mat. Res. Bull 8, 153 (1973)], respectively. The RVB Ansatz implemented in TurboRVB has a large variational freedom, including the Jastrow correlated Slater determinant as its simplest, but nontrivial case. Moreover, it has the remarkable advantage of remaining with an affordable computational cost, proportional to the one spent for the evaluation of a single Slater determinant. Therefore, its application to large systems is computationally feasible. The WF is expanded in a localized basis set. Several basis set functions are implemented, such as Gaussian, Slater, and mixed types, with no restriction on the choice of their contraction. The code implements the adjoint algorithmic differentiation that enables a very efficient evaluation of energy derivatives, comprising the ionic forces. Thus, one can perform structural optimizations and molecular dynamics in the canonical NVT ensemble at the VMC level. For the electronic part, a full WF optimization (Jastrow and antisymmetric parts together) is made possible, thanks to state-of-the-art stochastic algorithms for energy minimization. In the optimization procedure, the first guess can be obtained at the mean-field level by a built-in DFT driver. The code was efficiently parallelized by using a hybrid MPI-OpenMP protocol, which is also an ideal environment for exploiting the computational power of modern Graphics Processing Unit accelerators.

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

Experimental study on the role of the target electron temperature as a key parameter linking recycling to plasma performance in JET-ILW*

Changes in global and edge plasma parameters (H98(y,2), dimensionless collisionality ν*, core density peaking, separatrix density ne,sep) with variations in the D2 fueling rate and divertor configuration are unified into a single trend when mapped to $\langle$Te,ot$\rangle$, the spatially averaged spectroscopically derived outer target electron temperature. Dedicated JET with the ITER-like wall (JET-ILW) experiments in combination with an extended JET-ILW database of unseeded low-triangularity H-mode plasmas spanning a wide range of D2 fueling rates, Ip, Bt and heating power have demonstrated the importance of $\langle$Te,ot$\rangle$ as a key physics parameter linking the recycling particle source and detachment with plasma performance. The remarkably robust H98(y,2) trend with $\langle$Te,ot$\rangle$ is connected to a strong inverse correlation between $\langle$Te,ot$\rangle$, ne,sep and ν*, thus directly linking changes in the divertor recycling moderated by $\langle$Te,ot$\rangle$ with the previously established relationship between ν*, core density peaking and core pressure resulting in a degradation in core plasma performance with decreasing $\langle$Te,ot$\rangle$ (increasing ν*). Additionally, a strong inverse correlation between the separatrix to pedestal density ratio, ne,sep/ne,ped, and $\langle$Te,ot$\rangle$ is also established, with the rise in ne,sep/ne,ped saturating at $\langle$Te,ot$\rangle$ > 10 eV. A strong reduction in H98(y,2) is observed as $\langle$Te,ot$\rangle$ is driven from 30 to 10 eV via additional D2 gas fueling, while the divertor remains attached. Consequently, the pronounced performance degradation in attached divertor conditions has implications for impurity seeding radiative divertor scenarios, in which H98(y,2) is already low (~0.7) before impurities are injected into the plasma since moderate gas fueling rates are required to promote high divertor neutral pressure. A favorable pedestal pressure, pe,ped, dependence on Ip has also been observed, with an overall increase in pe,ped at Ip = 3.4 MA as $\langle$Te,ot$\rangle$ is driven down from attached to high-recycling divertor conditions. In contrast, pe,ped is reduced with decreasing $\langle$Te,ot$\rangle$ in the lower Ip branches. Further work is needed to (i) clarify the potential role of edge opacity on the observed favorable pedestal pressure Ip scaling; as well as to (ii) project the global and edge plasma performance trends with $\langle$Te,ot$\rangle$ to reactor-scale devices to improve predictive capability of the coupling between recycling and confined plasma fueling in what are foreseen to be more opaque edge plasma conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Aspects of the normal state resistivity of cuprate superconductors

Planar normal state resistivity data taken from three families of cuprate superconductors are compared here with theoretical calculations from the recent extremely correlated Fermi liquid theory (ECFL). The two hole doped cuprate materials LSCO and BSLCO and the electron doped material LCCO have yielded rich data sets at several densities δ and temperatures T, thereby enabling a systematic comparison with theory. The recent ECFL resistivity calculations for the highly correlated t-t'-J model by us give the resistivity for a wide set of model parameters. After using X-ray diffraction and angle resolved photoemission data to fix parameters appearing in the theoretical resistivity, only one parameter, the magnitude of the hopping t, remains undetermined. For each data set, the slope of the experimental resistivity at a single temperature-density point is sufficient to determine t, and hence the resistivity on absolute scale at all remaining densities and temperatures. This procedure is demonstrated to give a fair account of the entire data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Molecular basis for higher affinity of SARS-CoV-2 spike RBD for human ACE2 receptor

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has caused substantially more infections, deaths, and economic disruptions than the 2002-2003 SARSCoV. The key to understanding SARS-CoV-2's higher infectivity lies partly in its host receptor recognition mechanism. Experiments show that the human angiotensin converting enzyme 2 (ACE2) protein, which serves as the primary receptor for both CoVs, binds to the receptor binding domain (RBD) of CoV-2's spike protein stronger than SARS-CoV's spike RBD. The molecular basis for this difference in binding affinity, however, remains unexplained from X-ray structures. To go beyond insights gained from X-ray structures and investigate the role of thermal fluctuations in structure, we employ all-atom molecular dynamics simulations. Microseconds-long simulations reveal that while CoV and CoV-2 spike-ACE2 interfaces have similar conformational binding modes, CoV-2 spike interacts with ACE2 via a larger combinatorics of polar contacts, and on average, makes 45% more polar contacts. Correlation analysis and thermodynamic calculations indicate that these differences in the density and dynamics of polar contacts arise from differences in spatial arrangements of interfacial residues, and dynamical coupling between interfacial and non-interfacial residues. Furthermore, these results recommend that ongoing efforts to design spike-ACE2 peptide blockers will benefit from incorporating dynamical information as well as allosteric coupling effects.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

A direct numerical simulation of Jet A flame kernel quenching

The safe operation of aeronautical engines requires an understanding of flame ignition, propagation and extinction. In this study, direct numerical simulations are performed using a 29 species reduced chemical mechanism for jet fuel surrogate Jet A to understand the flame quenching process. Here, initially laminar spherical flames of varying sizes and equivalence ratios are subject to an identical periodic domain of decaying and isotropic high intensity turbulence with a turbulent Reynolds number of 2400. All cases become quenched, except for the larger kernel with lower Karlovitz number. An analysis of the flame structure shows broadened preheat zone, flame shortening on the product side, differential species diffusion and partial fuel pyrolysis in the fresh mixture. Two extinction mechanisms are identified arising from flame shortening and high flame stretch. Flame shortening occurs due to turbulence-chemistry interactions that resemble the flame–flame interaction in a laminar counterflow reactant-to-reactant configuration, which contorts and breaks up the ignition kernel. Flame stretch is a local effect that attenuates the heat release rate and causes the flame to retreat towards the product mixtures, similar to what has been observed for reactant-to-product laminar counterflow flames. Chemical explosive mode analysis was also performed to quantify the flame structure and local combustion mode. The diffusion–reaction balance in pinched-off flame islands favors extinction of these smaller structures, while auto-ignition modes are observed within the flame kernel after fresh mixture is engulfed and preheated in the product kernel. Statistics of the density-weighted displacement speed conditional on local combustion mode indicates strong correlation between the local extinction mode and negative displacement speed. The local balance between diffusion and reaction ultimately determines the propensity for local extinction in both laminar and turbulent flames, the extent of which has an impact on global flame propagation.

42 ENGINEERING↗
DOE OSTI2021

Unified descriptor for enhanced critical heat flux during pool boiling of hemi-wicking surfaces

Boiling heat transfer is dictated by interfacial phenomena at the three-phase contact line where vapor bubbles form on the surface. Structured surfaces have shown significant enhancement in critical heat flux (CHF) during pool boiling by tailoring interfacial phenomena. This CHF enhancement has been primarily explained by two structural effects: roughness, which extends the contact line length at the bubble base, and wickability, the ability to imbibe liquid through surface structures by capillary pumping. Here, in this work, we show that CHF enhancement on structured surfaces cannot be described by roughness or wickability alone. This result was confirmed using systematically designed micropillar surfaces with controlled roughness and wickability. Further, we performed a scaling analysis and derived a unified descriptor, which represents the combined effects of thin film density and volumetric wicking rate. This unified descriptor shows a reasonable correlation with CHF values with our experiments and literature data. This work provides important insights in understanding the role of surface structures on CHF enhancement, thereby providing guidelines for the systematic design of surface structures for enhanced pool boiling heat transfer.

42 ENGINEERING↗
DOE OSTI2025

Understanding the impact of the gas diffusion layer structure on catalyst utilization in the PEM water electrolyzer

A multiphysics half-cell model of a polymer electrolyte membrane water electrolyzer (PEMWE) was developed to probe impacts of the detailed 3-dimensional pore structure of the gas diffusion layer (GDL) on performance characteristics. We show that pores in the titanium GDL mesh led to significant underutilization of the catalyst layer (CL), with only 45% of the catalyst effectively utilized. This contradicts the assumption of uniform electron flow across the CL, as shown in graphical abstract (a), as near-zero current was observed near GDL pore regions and the current distribution in CL was influenced by GDL structure, as shown in graphical abstract (b). Instead, oxygen generation was primarily concentrated under the solid titanium regions, diffusing out around the pore walls. High current density peaks were also noted at the GDL-catalyst contact, correlating with degradation hotspots that were directly observed in companion experiments, as shown in graphical abstract (c). Collectively, these findings point to the critical importance of the heterogeneous GDL porous architecture not only for PEMWE efficiency but also for uneven degradation of the CL.

COMSOL↗
DOE OSTI2025

Processing-dependent chemical ordering in Cu 3 Au characterized via non-destructive Bragg coherent diffraction imaging

Of current importance for alloy design is controlling chemical ordering through processing routes to optimize an alloy's mechanical properties for a desired application. However, characterization of chemical ordering remains an ongoing challenge, particularly when nondestructive characterization is needed. Here, in this study, Bragg coherent diffraction imaging is used to reconstruct morphology and lattice displacement in model Cu 3 Au nanocrystals that have undergone different heat treatments to produce variation in chemical ordering. The magnitudes and distributions of the scattering amplitudes (proportional to electron density) and lattice strains within these crystals are then analyzed to correlate them to the expected amount of chemical ordering present. Nanocrystals with increased amounts of ordering are found to generally have less extreme strains present and reduced strain distribution widths. In addition, statistical correlations are found between the spatial arrangement of scattering amplitude and lattice strains.

Warren, Nathaniel [Pennsylvania State Univ., Unive↗
DOE OSTI2020

Isolating the Contributions of Specific Network Sites to the Diffuse Vibrational Spectrum of Interfacial Water with Isotopomer-Selective Spectroscopy of Cold Clusters

Decoding the structural information contained in the interfacial vibrational spectrum of water requires understanding how the spectral signatures of individual water molecules respond to their local hydrogen bonding environments. In this study, we isolated the contributions for the five classes of sites that differ according to the number of donor (D) and acceptor (A) hydrogen bonds that characterize each site. These patterns were measured by exploiting the unique properties of the water cluster cage structures formed in the gas phase upon hydration of a series of cations M+·(H2O)n (M = Li, Na, Cs, NH4, CH3NH3, H3O, and n = 5, 20-22). This selection of ions was chosen to systematically express the A, AD, AAD, ADD, and AADD hydrogen bonding motifs. The spectral signatures of each site were measured using two-color, IR-IR isotopomer-selective photofragmentation vibrational spectroscopy of the cryogenically cooled, mass selected cluster ions in which a single intact H2O is introduced without isotopic scrambling, an important advantage afforded by the cluster regime. The resulting patterns provide an unprecedented picture of the intrinsic line shapes and spectral complexities associated with excitation of the individual OH groups, as well as the correlation between the frequencies of the two OH groups on the same water molecule, as a function of network site. The properties of the surrounding water network that govern this frequency map are evaluated by dissecting electronic structure calculations that explore how changes in the nearby network structures, both within and beyond the first hydration shell, affect the local frequency of an OH oscillator. The qualitative trends are recovered with a simple model that correlates the OH frequency with the network-modulated local electron density in the center of the OH bond.

Yang, Nan↗
DOE OSTI2021

Prediction of Transition-State Scaling Relationships and Universal Transition-State Vibrational and Entropic Correlations for Dehydrogenations

Linear scaling relationships (LSRs) and Brønsted–Evans–Polanyi (BEP) or transition-state scaling (TSS) relations aid with the prediction of electronic energies. However, temperature effects and pre-exponentials are often taken as constants across metal surfaces or a homologous series. Vibrational scaling relationships (VSRs) offer a way to determine such parameters. Transition-state VSRs (TSVSRs) between local minima and transition states of AH X (A = C, N, O) surface diffusions correlate with BEP relations and broaden to thermochemical property scaling. Using density functional theory, we extend TSVSRs to AHX dehydrogenation reactions on transition-metal surfaces, relating vibrational modes of local minima to transition states. We first predict the slopes of the TSS relations by incorporating bond angles using the Slater–Koster structure factors and hybridization through crystal orbital overlap population analysis and energy overlap integrals between adsorbates and metal surfaces. Additionally, we uncover universal thermochemical property scaling, enabling the estimation of entropies and temperature corrections to enthalpies across a homologous series. Here, we demonstrate both significant vibrational corrections in reactions with low intrinsic electronic barriers and considerable variation in the pre-exponential of a simple dehydrogenation reaction across metals and AH X adsorbates.

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

Role of Dynamic Polarization Interactions in the Electrical Double Layer at Calcite (104) Interfaces with Aqueous Solutions

Reactions at mineral interfaces with aqueous solutions control many geochemical and biogeochemical processes in the Earth’s critical zone. At the molecular level, insights into important properties such as the structure of the electrical double layer (EDL) at specific mineral interfaces continue to improve due to the increasing fidelity of laboratory instrumentation and computational approaches. However, molecular simulation approaches suffer from limited reach into relevant scales of time, length, and system complexity. To span this gap, a novel hybrid approach that couples first principle plane-wave density functional theory (DFT) with classical DFT (cDFT) is demonstrated and applied to calcite (104) interfaces with various electrolytes. In this approach, a region of interest described using DFT interacts with the surrounding medium described using cDFT to arrive at a self-consistent ground state. Benchmarking against experimental observations and entirely first principle DFT simulations demonstrates that this hybrid model efficiently encompasses the key short-range and collective interactions in the EDL. Simulations of calcite (104)/solution interfaces reveal the key static and dynamic polarization interactions that give rise to structuring of ions and water. Ion hydration interactions have the strongest effect on the depth of the first minimum in the density distribution of counterions at the surface, and the position and width of the first density peak is largely determined by the strength of ion-correlation forces. Finer details of ion distributions are controlled by mutual polarization of the calcite surface and interfacial electrolyte. Finally, this new ability to efficiently and rigorously predict EDL structure at mineral surfaces in contact with complex solutions paves the way to accurately modeling sorption, nucleation, dissolution, and growth in realistic systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗