SEARCH · Engineering Papers

Results for “long-range dispersion”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Resolution of zigzag magnetic correlations in Na-deficient Na x IrO 3 without long-range ordering

The materials search for Kitaev quantum spin liquids led to the discovery of many honeycomb lattice materials. Much attention has been paid to materials without magnetic order down to the lowest temperatures. The newly synthesized Na-deficient Na x IrO 3 has been found to bear no sign of long-range magnetic order above 1 K from physical property measurements. In this paper, we report momentum-and energy-resolved excitation spectra in Na-deficient Na x IrO 3 measured using a resonant inelastic x-ray scattering spectrometer. Orbital excitation spectra show that the octahedral and trigonal crystal field splittings are larger in Na x IrO 3 than in Na 2 IrO 3 . On the other hand, the low-energy spectrum at low temperature shows a wave-vector dispersion and a spectral weight distribution that are similar to those of Na 2 IrO 3 , revealing that the two-dimensional zigzag magnetic correlations in Na x IrO 3 are similar to those in Na 2 IrO 3 in terms of the ordered magnetic moment direction and three dynamically fluctuating zigzag orders. The azimuth angle dependence of the low-energy spectrum corroborates these results. The two-dimensional zigzag magnetic correlations rapidly weaken until 50 K. At high temperatures, the spectral weight distribution of the low-energy excitation resembles that of the pure Kitaev model, indicating that the Kitaev interaction dominates the dynamic magnetic response at high temperature. Here we suggest that the larger crystal field and distortion and the weakened longer-range Heisenberg exchange interactions due to the Na deficiency in Na x IrO 3 contribute to bring Na x IrO 3 away from the zigzag long-range magnetic order phase.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Theoretical study of the crystal and electronic properties of α - RuI 3

The material α-RuCl 3 , with a two-dimensional Ru honeycomb sublattice, has attracted considerable attention because it may be a realization of the Kitaev quantum spin liquid. Recently, a new honeycomb material, α-RuCl 3 , was prepared under moderately high pressure, and it is stable under ambient conditions. However, different from α-RuCl 3 , α-RuI 3 was reported to be a paramagnetic metal without long-range magnetic order down to 0.35 K. Here, the structural and electronic properties of the quasi-two-dimensional α-RuI 3 are theoretically studied. First, based on first-principles density functional theory calculations, the ABC stacking honeycomb-layer R$\bar{3}$ (No. 148) structure is found to be the most likely stacking order for α-RuI 3 along the c axis. Furthermore, both R$\bar{3}$ and P$\bar{3}$1c are dynamically stable because no imaginary frequency modes were obtained in the phononic dispersion spectrum without Hubbard U. Moreover, the different physical behavior of α-RuI 3 compared to α-RuCl 3 can be understood naturally. The strong hybridization between Ru 4d and I 5p orbitals decreases the “effective” atomic Hubbard repulsion, leading the electrons of RuI3 to be less localized than in RuCl 3 . As a consequence, the effective electronic correlation is reduced from Cl to I, leading to the metallic nature of α-RuI 3 . Based on the DFT+U (Ueff=2 eV) plus spin-orbital coupling, we obtained a spin-orbit Mott insulating behavior for α-RuCl 3 and, with the same procedure, a metallic behavior for α-RuI 3 , in good agreement with experimental results. Furthermore, when introducing large (unrealistic) U eff =6 eV, the spin-orbit Mott gap opens in α-RuI 3 as well, supporting the physical picture we are proposing. Our results provide guidance to experimentalists and theorists working on two-dimensional transition metal tri-iodide layered materials.

2-dimensional systems↗

DOE OSTI2021

Simplified tuning of long-range corrected density functionals for use in symmetry-adapted perturbation theory

Long considered a failure, second-order symmetry-adapted perturbation theory (SAPT) based on Kohn–Sham orbitals, or SAPT0(KS), can be resurrected for semiquantitative purposes using long-range corrected density functionals whose asymptotic behavior is adjusted separately for each monomer. Here, as in other contexts, correct asymptotic behavior can be enforced via “optimal tuning” based on the ionization energy theorem of density functional theory, but the tuning procedure is tedious, expensive for large systems, and comes with a troubling dependence on system size. Here, we show that essentially identical results are obtained using a fast, convenient, and automated tuning procedure based on the size of the exchange hole. In conjunction with “extended” (X)SAPT methods that improve the description of dispersion, this procedure achieves benchmark-quality interaction energies, along with the usual SAPT energy decomposition, without the hassle of system-specific tuning.

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA NTRS2022

Investigating Particle-Particle Electrostatic Effects on Charged Lunar Dust Transport via Discrete Element Modeling

NASA surface exploration missions have always seen negative effects of dust including the Apollo missions. The astronaut-witnessed unusual behavior of the dust particles that surround the vehicle after engine cutoff has the potential to have more of an influence on surface systems dust loading than the high velocity lunar rocket plume ejecta in the landing process. The levitation and transport of the fine components of regolith on lunar surface has been linked to electrostatic effects and electric field, but so far there is no accurate model considering the inter-particle electrostatic interactions, especially when the particles are charged by rocket plume or other mechanical interactions due to exploration activities. This study is proposed to investigate the dynamics of charged lunar regolith with a discrete element method (DEM) approach focusing on the inter-particle interactions and contact charge transfer. The grain dynamics is coupled with mechanical and electrical particle interactions, and both short- and long-range interactions between spherical particles are incorporated. A tribo-charging model based on instantaneous collisions between particles is adopted and validated by comparing the simulation results to existing experimental data. Sensitivity analysis is conducted to quantify the effects of initial charge, tribo-charging, and E-field on transport of lunar dust based on JSC-1 simulants with a radius of 50 lm. DEM simulations are also conducted in a near realistic lunar environment with the estimations of initial conditions that shows the difference in position and velocity distributions between charged particles and uncharged particles. The results indicate that the charged dust particles have higher dispersion of position and velocity by several orders of magnitude due to electrostatic effects. This provides a potential explanation for the phenomena of the approximately 30 s dust lofting following Apollo Lunar Module landing.

Lunar dust↗

DOE OSTI2021

Spatiotemporal spread of perturbations in power-law models at low temperatures: Exact results for classical out-of-time-order correlators

Here, we present exact results for the classical version of the out-of-time-order commutator (OTOC) for a family of power-law models consisting of N particles in one dimension and confined by an external harmonic potential. These particles are interacting via power-law interaction of the form ∝ $Σ_{i,j=(i≠j)}^{N}$ |x i -x j | -k for every k>1 where x i is the position of the i th particle. We present numerical results for the OTOC for finite N at low temperatures and short enough times so that the system is well approximated by the linearized dynamics around the many-body ground state. In the large-N limit, we compute the ground-state dispersion relation in the absence of external harmonic potential exactly and use it to arrive at analytical results for OTOC. We find excellent agreement between our analytical results and the numerics. We further obtain analytical results in the limit where only linear and leading nonlinear (in momentum) terms in the dispersion relation are included. The resulting OTOC is in agreement with numerics in the vicinity of the edge of the “light cone.” We find remarkably distinct features in OTOC below and above k=3 in terms of going from non-Airy behavior (1 < k < 3) to an Airy universality class (k>3). We present certain additional rich features for the case k=2 that stem from the underlying integrability of the Calogero-Moser model. We present a field theory approach that also assists in understanding certain aspects of OTOC such as the sound speed. Our findings are a step forward towards a more general understanding of the spatiotemporal spread of perturbations in long-range interacting systems.

1-dimensional systems↗

NASA NTRS2017

Basalt Weathering in a Cold and Icy Climate: Three Sisters, Oregon as an Analog for Early Mars

There is abundant evidence for liquid water on early Mars, but the debate remains whether early Mars was warm and wet or cold and icy with punctuated periods of melting. To further investigate the hypothesis of a cold and icy early Mars, we collected rocks and sediments from the Collier and Diller glacial valleys in the Three Sisters volcanic complex in Oregon. We analyzed rocks and sediments with X-ray diffraction (XRD), scanning and transmission electron microscopies with energy dispersive spectroscopy (SEM, TEM, EDS), and visible, short-wave infrared (VSWIR) and thermal-IR (TIR) spectroscopies to characterize chemical weathering and sediment transport through the valleys. Here, we focus on the composition and mineralogy of the weathering products and how they compare to those identified on the martian surface. Phyllosilicates (smectite), zeolites, and poorly crystalline phases were discovered in pro- and supra-glacial sediments, whereas Si-rich regelation films were found on hand samples and boulders in the proglacial valleys. Most phyllosilicates and zeolites are likely detrital, originating from hydrothermally altered units on North Sister. TEM-EDS analyses of the <2 um size fraction of glacial flour samples demonstrate a variety of poorly crystalline (i.e., no long-range crystallographic order) phases: iron oxides, devitrified volcanic glass, and Fe-Si-Al phases. The CheMin XRD on the Curiosity rover in Gale crater has identified significant amounts of X-ray amorphous materials in all samples measured to date. The amorphous component is likely a combination of silicates, iron oxides, and sulfates. Although we have not yet observed amorphous sulfate in the samples from Three Sisters, the variety of poorly crystalline weathering products found at this site is consistent with the variable composition of the X-ray amorphous component identified by CheMin. We suggest that these amorphous phases on Mars could have formed in a similarly cold and icy environment.

Rampe, E. B.↗

DOE OSTI2020

The effect of local chemical ordering on Frank-Read source activation in a refractory multi-principal element alloy

In this work, we investigate the operation of Frank-Read (FR) sources in a refractory multi-principal element alloy (MPEA). Simulations of discrete dislocation motion in MPEAs is enabled by the development of a phase field dislocation dynamics model that treats the atomic-scale fluctuations in lattice energies across the glide plane, present due to local ordering in the chemical composition within the nominally random MPEA atomic structure. We consider, through simulation, a range of length scales over which ordering occurs, varying from short-range lengths, a few times dislocation core width, to long-range lengths, an order of magnitude longer than the core. Characteristic of this body-centered cubic MPEA, the simulations also include screw/edge character dependence in glide resistance, as informed by atomic scale simulation. The critical stresses to activate the source for the same source size are found to be statistically distributed, as a direct consequence of the underlying variation in lattice energy. Furthermore, analysis of critical state for activating edge and screw FR sources in the MPEA reveals that FR source operation occurs via a two-step mechanism, involving athermal kink-pair formation, unlike the conventional FR source operation in a material with no composition fluctuations. This mechanism lowers the average critical stress required to activate the FR source and causes the statistical dispersion in critical stress to depend on the range of composition ordering. More importantly, it leads to a more severe dependence of source strength on FR source length than predicted by line tension alone.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Exfoliation of Cu-Containing Poly(triazine imide): From Three-Dimensional to Two-Dimensional Particle Morphology

Controlling the morphological parameters of extended covalent organic frameworks remains challenging and represents an important yet often elusive metric of consideration. Typically, carbon nitride materials possess local ordering but remain largely amorphous in terms of their long-range order and orientation. This study probes the synthesis of a crystalline carbon nitride, poly(triazine imide) lithium bromide which possesses an atomically-precise extended structure, and demonstrates its exfoliation into a two-dimensional hexagonal sheet-like morphology. Furthermore, a previously unreported carbon nitride material, poly(triazine imide) copper bromide, or PTI-CuBr, was developed through an additional flux-assisted cation-exchange process and is shown to retain its internal Cu cations during solvothermal exfoliation. Characterization by dynamic light scattering and high-angle annular dark-field scanning electron microscopy reveals the morphological changes and captures the high aspect ratio of the thin carbon nitride sheets with <10 nm thickness while maintaining hundreds of nm in width. Additional characterization by energy-dispersive spectroscopy and X-ray photoelectron spectroscopy confirms that the Cu:Br:N molar ratio was maintained within the extended layers throughout the exfoliation process. This top-down synthesis approach differs from typical methods that isolate thin sheets for subsequent metal−cation coordination and illustrates the importance of maintaining oxygen-free conditions to minimize copper clustering. Thus, this new approach is demonstrated to provide a consistent and more homogeneous occupancy of the PTI pore spaces throughout the carbon nitride framework.

Exfoliation↗

DOE OSTI2023

Following the Higgs mode across the BCS-BEC crossover in two dimensions

Although substantial effort has been dedicated to analyzing the Higgs (amplitude) mode in superconducting systems, there are relatively few studies of the Higgs peak's dispersion and width, quantities which are observable in spectroscopic measurements. These properties can be obtained from the location of the pole of the amplitude (Higgs) susceptibility in the lower half plane of complex frequency. We analyze the behavior of the Higgs mode in a two-dimensional neutral fermionic superfluid at T = 0 throughout the crossover from Bardeen-Cooper-Schrieffer (BCS) to Bose-Einstein condensation (BEC) behavior. Here, this occurs when the dressed chemical potential μ changes sign from positive to negative. For μ > 0, we find a pole in the Higgs susceptibility in the lower half plane of frequency and demonstrate that it leads to a well-defined peak in spectroscopic probes. For μ < 0, the pole still exists but is “hidden,” not giving rise to a peak in spectroscopic probes. Extending this analysis to a charged superconductor, we find that the Higgs mode is unaffected by the long-range Coulomb interaction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Hybrid Nanocrystals of Small Molecules and Chemically Disordered Polymers

Organic crystals formed by small molecules can be highly functional but are often brittle or insoluble structures with limited possibilities for use or processing from a liquid phase. A possible solution is the nanoscale integration of polymers into organic crystals without sacrificing long-range order and therefore function. This enables the organic crystals to benefit from the advantageous mechanical and chemical properties of the polymeric component. We report here on a strategy in which small molecules co-crystallize with side chains of chemically disordered polymers to create hybrid nanostructures containing a highly ordered lattice. Synchrotron X-ray scattering, absorption spectroscopy, and coarse-grained molecular dynamics simulations reveal that the polymer backbones form an “exo-crystalline” layer of disordered chains that wrap around the nanostructures becoming a handle for interesting properties. The morphology of this “hybrid bonding polymer” nanostructure is dictated by the competition between the polymers’ entropy and the enthalpy of the lattice allowing for control over the aspect ratio of the nanocrystal by changing the degree of polymer integration. Further, we observed that nanostructures with an exo-crystalline layer of polymer exhibit enhanced fracture strength, self-healing capacity, and dispersion in water, which benefits their use as light-harvesting assemblies in photocatalysis. Guided by computation, future work could also explore these hybrid nanostructures as building blocks for functional materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Nature of polariton transport in a Fabry-Perot cavity

Fabry-Perot microcavities can strongly enhance interactions between light and molecules, leading to the formation of hybrid light-matter states known as polaritons. Polaritons possess much smaller effective masses and much larger group velocities when the molecules are resonant with cavity modes that have finite (nonzero) in-plane wave vectors, giving rise to the possibilities of long-range and ultrafast ballistic transport. In this paper, we present the results of numerical simulations of the ultrafast ballistic transport phenomenon in real space and time during and after initialization with a short, spatially localized pulse. Here, we address this problem with two approaches: the standard transfer-matrix method (TMM) for planar structures and a second simulation based on the numerical solution of the Maxwell-Bloch equations that can be used for general configurations (complex metasurfaces, for instance) and excitation modes. The agreement between the TMM and the full numerical calculation with the Maxwell-Bloch equations when applied to multilayer planar configurations provides proof of the validity of both approaches for the present analyses. Overall, we find that the transport of the molecular excitons inside the cavity synchronize with the evolution of the enhanced electromagnetic field inside the cavity. Moreover, the synchronized transport rate is in good agreement with the group velocities predicted from a calculated dispersion relation across a wide range of frequencies. Finally, we relate the group velocity to the Hopfield coefficient obtained from quantum modeling and suggest that the dependence of light-matter coupling on the in-plane wave vector can be an important but overlooked factor for understanding the transport behavior of polaritons. These simulations provide an intuitive tool for understanding the collective motion of light and excitons and helps us to better understand how experimental observations of polaritons should be interpreted.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1999

Microgravity Experiments to Evaluate Electrostatic Forces in Controlling Cohesion and Adhesion of Granular Materials

The bulk behavior of dispersed, fluidized, or undispersed stationary granular systems cannot be fully understood in terms of adhesive/cohesive properties without understanding the role of electrostatic forces acting at the level of the grains themselves. When grains adhere to a surface, or come in contact with one another in a stationary bulk mass, it is difficult to measure the forces acting on the grains, and the forces themselves that induced the cohesion and adhesion are changed. Even if a single grain were to be scrutinized in the laboratory, it might be difficult, perhaps impossible, to define the distribution and character of surface charging and the three-dimensional relationship that charges (electrons, holes) have to one another. The hypothesis that we propose to test in microgravity (for dielectric materials) is that adhesion and cohesion of granular matter are mediated primarily by dipole forces that do not require the presence of a net charge; in fact, nominally electrically neutral materials should express adhesive and cohesive behavior when the neutrality results from a balance of positive and negative charge carriers. Moreover, the use of net charge alone as a measure of the electrical nature of grain-to-grain relationships within a granular mass may be misleading. We believe that the dipole forces arise from the presence of randomly-distributed positive and negative fixed charge carriers on grains that give rise to a resultant dipole moment. These dipole forces have long-range attraction. Random charges are created whenever there is triboelectrical activity of a granular mass, that is, whenever the grains experience contact/separation sequences or friction.

Marshall, J.↗

DOE OSTI2022

Particle phase-state variability in the North Atlantic free troposphere during summertime is determined by atmospheric transport patterns and sources

Abstract. Free tropospheric aerosol particles have important but poorly constrained climate effects due to transformations of their physicochemical properties during long-range transport. In this study, we investigate the chemical composition and provide an overview of the phase states of individual particles that have undergone long-range transport over the North Atlantic Ocean in June and July 2014, 2015, and 2017 to the Observatory of Mount Pico (OMP) in the Azores. The OMP is an ideal site for studying long-range-transported free tropospheric particles because local emissions have a negligible influence and contributions from the boundary layer are rare. We used the FLEXible PARTicle Lagrangian particle dispersion model (FLEXPART) to determine the origins and transport trajectories of sampled air masses and found that most of them originated from North America and recirculated over the North Atlantic Ocean. The FLEXPART analysis showed that the sampled air masses were highly aged (average plume age >10 d). Size-resolved chemical compositions of individual particles were probed using computer-controlled scanning electron microscopy with an energy-dispersive X-ray spectrometer (CCSEM-EDX) and scanning transmission X-ray microscopy with near-edge X-ray absorption fine structure spectroscopy (STXM-NEXAFS). CCSEM-EDX results showed that the most abundant particle types were carbonaceous (∼ 29.9 % to 82.0 %), sea salt (∼ 0.3 % to 31.6 %), and sea salt with sulfate (∼ 2.4 % to 31.5 %). We used a tilted stage interfaced within an environmental scanning electron microscope (ESEM) to determine the phase states of individual submicron particles. We found that most particles (∼ 47 % to 99 %) were in the liquid state at the time of collection due to inorganic inclusions. Moreover, we also observed substantial fractions of solid and semisolid particles (∼ 0 % to 30 % and ∼ 1 % to 42 %, respectively) during different transport patterns and events, reflecting the particles' phase-state variability for different atmospheric transport events and sources. Combining phase state measurements with FLEXPART CO tracer analysis, we found that wildfire-influenced plumes can result in particles with a wide range of viscosities after long-range transport in the free troposphere. We also used temperature and RH values extracted from the Global Forecast System (GFS) along the FLEXPART-simulated path to predict the phase state of the particles during transport and found that neglecting internal mixing with inorganics would lead to an overestimation of the viscosity of free tropospheric particles. Our findings warrant future investigation aiming at the quantitative assessment of the influence of internal mixing on the phase states of the individual particles. This study also provides insights into the chemical composition and phase state of free tropospheric particles, which can help models to reduce uncertainties about the effects of ambient aerosol particles on climate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Integrated Quantum-Classical Protocol for the Realistic Description of Solvated Multinuclear Mixed-Valence Transition-Metal Complexes and Their Solvatochromic Properties

Linear cyanide-bridged polymetallic complexes, which undergo photoinduced metal-to-metal charge transfer, represent prototypical systems for studying long-range electron-transfer reactions and understanding the role played by specific solute–solvent interactions in modulating the excited-state dynamics. Here, to tackle this problem, while achieving a statistically meaningful description of the solvent and of its relaxation, one needs a computational approach capable of handling large polynuclear transition-metal complexes, both in their ground and excited states, as well as the ability to follow their dynamics in several environments up to nanosecond time scales. Here, we present a mixed quantum classical approach, which combines large-scale molecular dynamics (MD) simulations based on an accurate quantum mechanically derived force field (QMD-FF) and self-consistent QMD polarized point charges, with IR and UV–vis spectral calculations to model the solvation dynamics and optical properties of a cyano-bridged trinuclear mixed-valence compound (trans-[(NC) 5 Fe III (μ-CN)Ru II (pyridine) 4 (μ-NC)Fe III (CN) 5 ] 4– ). We demonstrate the reliability of the QMD-FF/MD approach in sampling the solute conformational space and capturing the local solute–solvent interactions by comparing the results with higher-level quantum mechanics/molecular mechanics (QM/MM) MD reference data. The IR spectra calculated along the classical MD trajectories in different solvents correctly predict the red shift of the CN stretching band in the aprotic medium (acetonitrile) and the subtle differences measured in water and methanol, respectively. By explicitly including the solvent molecules around the cyanide ligands and calculating the thermal averaged absorption spectra using time-dependent density functional theory calculations within the Tamm–Dancoff approximation, the experimental solvatochromic shift is quantitatively reproduced going from water to methanol, while it is overestimated for acetonitrile. This discrepancy can likely be traced back to the lack of important dispersion interactions between the solvent cyano groups and the pyridine substituents in our micro solvation model. The proposed protocol is applied to the ground state in water, methanol, and acetonitrile and can be flexibly generalized to study excited-state nonequilibrium solvation dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Phonon gas model for thermal conductivity of dense, strongly interacting liquids

Developing predictive thermal property models for liquids based on microscopic principles has been elusive. The difficulty is that liquids have gas-like and solid-like attributes that are at odds when considering the frameworks of microscopic models: Models for gases are simple due to randomness and low density, whereas models for crystalline solids rely on symmetry and long-range order for easier calculation. The short-range order in liquids does, however, provide structure to neighboring molecules similar to amorphous solids, and there have been recent advances indicating that collective vibrational modes store heat in liquids. Models combining Debye approximations from solid-state physics and Frenkel’s theory of liquids can accurately predict the heat capacity of liquids. Phonon-like dispersions in liquids have also been widely observed in neutron scattering experiments. These developments motivate us to propose a model where high-frequency vibrational modes, which travel at the speed of sound and have a mean free path on the order of the average intermolecular distance, conduct heat in liquids. We use this liquid phonon gas model to calculate the thermal conductivity of liquids with varying intermolecular interaction energies from strongest to weakest—Coulomb, hydrogen-bonding, Keesom, and London dispersion energy. Generally, the model is more accurate as the intermolecular interaction energy and density of liquids increase. The calculated thermal conductivity of Coulombic-bound molten sodium nitrate and hydrogen-bonded water is within 1.46% and 2.98% of the experimentally measured values, respectively, across their entire temperature ranges. Further modal analysis of the velocity and the mean free path of collective vibrations could establish the liquid phonon gas model as an accurate model for weakly interacting liquids as well.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Electronic structure theory on modeling short-range noncovalent interactions between amino acids

While short-range noncovalent interactions (NCIs) are proving to be of importance in many chemical and biological systems, these atypical bindings happen within the so-called van der Waals envelope and pose an enormous challenge for current computational methods. We introduce SNCIAA, a database of 723 benchmark interaction energies of short-range noncovalent interactions between neutral/charged amino acids originated from protein x-ray crystal structures at the “gold standard” coupled-cluster with singles, doubles, and perturbative triples/complete basis set [CCSD(T)/CBS] level of theory with a mean absolute binding uncertainty less than 0.1 kcal/mol. Subsequently, a systematic assessment of commonly used computational methods, such as the second-order Møller−Plesset theory (MP2), density functional theory (DFT), symmetry-adapted perturbation theory (SAPT), composite electronic-structure methods, semiempirical approaches, and the physical-based potentials with machine learning (IPML) on SNCIAA is carried out. It is shown that the inclusion of dispersion corrections is essential even though these dimers are dominated by electrostatics, such as hydrogen bondings and salt bridges. Overall, MP2, ωB97M-V, and B3LYP+D4 turned out to be the most reliable methods for the description of short-range NCIs even in strongly attractive/repulsive complexes. SAPT is also recommended in describing short-range NCIs only if the δMP2 correction has been included. The good performance of IPML for dimers at close-equilibrium and long-range conditions is not transferable to the short-range. We expect that SNCIAA will assist the development/improvement/validation of computational methods, such as DFT, force-fields, and ML models, in describing NCIs across entire potential energy surfaces (short-, intermediate-, and long-range NCIs) on the same footing.

Chemistry↗

DOE OSTI2023

Accurate prediction of global-density-dependent range-separation parameters based on machine learning

In this work, we develop an accurate and efficient XGBoost machine learning model for predicting the global-density-dependent range-separation parameter, ωGDD, for long-range corrected functional (LRC)-ωPBE. This ωGDDML model has been built using a wide range of systems (11 466 complexes, ten different elements, and up to 139 heavy atoms) with fingerprints for the local atomic environment and histograms of distances for the long-range atomic correlation for mapping the quantum mechanical range-separation values. The promising performance on the testing set with 7046 complexes shows a mean absolute error of 0.001 117 a0−1 and only five systems (0.07%) with an absolute error larger than 0.01 a0−1, which indicates the good transferability of our ωGDDML model. In addition, the only required input to obtain ωGDDML is the Cartesian coordinates without electronic structure calculations, thereby enabling rapid predictions. LRC-ωPBE(ωGDDML) is used to predict polarizabilities for a series of oligomers, where polarizabilities are sensitive to the asymptotic density decay and are crucial in a variety of applications, including the calculations of dispersion corrections and refractive index, and surpasses the performance of all other popular density functionals except for the non-tuned LRC-ωPBE. Finally, LRC-ωPBE (ωGDDML) combined with (extended) symmetry-adapted perturbation theory is used in calculating noncovalent interactions to further show that the traditional ab initio system-specific tuning procedure can be bypassed. The present study not only provides an accurate and efficient way to determine the range-separation parameter for LRC-ωPBE but also shows the synergistic benefits of fusing the power of physically inspired density functional LRC-ωPBE and the data-driven ωGDDML model.

Chemistry↗

DOE OSTI2025

Wildfire towers drive firebrand lofting: insights from coupled fire-atmosphere model simulations

Wildfire behavior is shaped by complex fire dynamics, with firebrands playing a critical role in spot fire ignition and fire spread. While previous studies have explored firebrand generation and transport, the specific role of towers and troughs from wildland fires in the lofting of firebrands remains unquantified. This study addresses that gap by using physics-based coupled fire-atmosphere model simulations to examine how wildfire towers (updrafts) and troughs (downdrafts) influence firebrand lofting. Our results show that the majority of firebrands (78.85%) are lofted from towers, where strong updrafts drive long-range transport. In contrast, only 21.15% of firebrands are lofted within troughs, where downdrafts cause most firebrands to fall near the fireline. We also find that firebrand size significantly influences lofting behavior, with smaller particles (1 mm radius) exhibiting the strongest correlation with updraft intensity. These findings highlight the dominant role of wildfire towers in promoting long-distance firebrand dispersal—an essential factor in rapid wildfire growth and wildland-urban interface (WUI) fire risks. By quantifying the relationship between firebrand lofting and fire-induced atmospheric features, this study provides critical insights to improve spot fire modeling, support mitigation planning, and enhance firefighter and WUI community safety in spot fire-prone regions.

54 ENVIRONMENTAL SCIENCES↗