Engineering Papers⌕ Search

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.

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

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Absence of Long-Range Magnetic Ordering in a Trirutile High-Entropy Oxide (Mn 0.2 Fe 0.2 Co 0.2 Ni 0.2 Cu 0.2 )Ta 1.92 O 6–δ

Functionalities of solid-state materials are usually considered to be dependent on their crystal structures. The limited structural types observed in the emerging high-entropy oxides put constraints on the exploration of their physical properties and potential applications. Herein, we synthesized the first high-entropy oxide in a trirutile structure, (Mn 0.2 Fe 0.2 Co 0.2 Ni 0.2 Cu 0.2 )Ta 1.92 O 6–δ , and investigated its magnetism. The phase purity and high-entropy nature were confirmed by powder Xray diffraction and energy-dispersive spectroscopy, respectively. X-ray photoelectron spectroscopy indicated divalent Mn, Co, Ni, and Cu along with trivalent Fe. Magnetic property measurements showed antiferromagnetic coupling and potential short-range magnetic ordering below ~4 K. The temperature-dependent heat capacity data measured under zero and high magnetic fields confirmed the lack of long-range magnetic ordering and a possible low-temperature phonon excitation. The discovery of the first trirutile high-entropy oxide opens a new pathway for studying the relationship between the highly disordered atomic arrangement and magnetic interaction. Furthermore, it provides a new direction for exploring the functionalities of highentropy oxides.

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

Analysis of phase stability and chemical segregation in the Mo-V alloys using a generalized embedded atom method potential

A new interatomic potential for the Mo-V system is introduced to facilitate the study of phase stability and mechanical properties at lower temperatures. This potential is based on a generalization of the embedded atom method and includes contributions from embedding energy, explicit two- and three-body interactions and nonlocal many-body interaction terms. The parameters of the potential are optimized by using data from ab initio density functional theory (DFT) calculations. The potential is rigorously validated across a range of physical properties, such as elastic constants, equation of states, phonon dispersion curves, point defect properties and melting temperatures for different compositions. Even though our potential is trained on a small dataset, its accuracy is comparable to available machine learning potentials for Mo and V. Furthermore, our results show that an ordered B2 phase is stable at low temperatures in alloys containing 50% V, but the solid solution phase is stable above 800 K. However, such long-range ordering is not observed in V-rich or Mo-rich alloys. In addition, our results show that V segregates to dislocation cores and grain boundaries.

36 MATERIALS SCIENCE↗
DOE OSTI2024

N and OH-Immobilized Cu 3 Clusters In Situ Reconstructed from Single-Metal Sites for Efficient CO 2 Electromethanation in Bicontinuous Mesochannels

Cu-based catalysts hold promise for electrifying CO 2 to produce methane, an extensively used fuel. However, the activity and selectivity remain insufficient due to the lack of catalyst design principles to steer complex CO 2 reduction pathways. Herein, we develop a concept to design carbon-supported Cu catalysts by regulating Cu active sites’ atomic-scale structures and engineering the carbon support’s mesoscale architecture. This aims to provide a favorable local reaction microenvironment for a selective CO 2 reduction pathway to methane. In situ X-ray absorption and Raman spectroscopy analyses reveal the dynamic reconstruction of nitrogen and hydroxyl-immobilized Cu 3 (N,OH-Cu 3 ) clusters derived from atomically dispersed Cu–N 3 sites under realistic CO 2 reduction conditions. The N,OH-Cu 3 sites possess moderate *CO adsorption affinity and a low barrier for *CO hydrogenation, enabling intrinsically selective CO 2 -to-CH 4 reduction compared to the C–C coupling with a high energy barrier. Importantly, a block copolymer-derived carbon fiber support with interconnected mesopores is constructed. The unique long-range mesochannels offer an H 2 O-deficient microenvironment and prolong the transport path for the CO intermediate, which could suppress the hydrogen evolution reaction and favor deep CO 2 reduction toward methane formation. Thus, the newly developed catalyst consisting of in situ constructed N,OH-Cu 3 active sites embedded into bicontinuous carbon mesochannels achieved an unprecedented Faradaic efficiency of 74.2% for the CO 2 reduction to methane at an industry-level current density of 300 mA cm –2 . In conclusion, this work explores effective concepts for steering desirable reaction pathways in complex interfacial catalytic systems via modulating active site structures at the atomic level and engineering pore architectures of supports on the mesoscale to create favorable microenvironments.

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

Establishing coherent momentum-space electronic states in locally ordered materials

Rich momentum-dependent electronic structure naturally arises in solids with long-range crystalline symmetry. Reliable and scalable quantum technologies rely on materials that are either not perfect crystals or non-crystalline, breaking translational symmetry. This poses the fundamental questions of whether coherent momentum-dependent electronic states can arise without long-range order, and how they can be characterized. Here we investigate Bi 2 Se 3 , which exists in crystalline, nanocrystalline, and amorphous forms, allowing direct comparisons between varying degrees of spatial ordering. Through angle-resolved photoemission spectroscopy, we show for the first time momentum-dependent band structure with Fermi surface repetitions in an amorphous solid. The experimental data is complemented by a model that accurately reproduces the vertical, dispersive features as well as the replication at higher momenta in the amorphous form. These results reveal that well-defined real-space length scales are sufficient to produce dispersive band structures, and that photoemission can expose the imprint of these length scales on the electronic structure.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Multi-Scale Integrated Monitoring System for Enhancing Methane Emission Detection, Quantification & Prediction

This report details the progress and findings of a comprehensive study on reviewing existing solutions, identifying technology gaps, and formulating an “all-in-one” integrated strategy for developing the next-generation multiscale methane monitoring and modeling platform, conducted under grant number DE-FE0032292. Co-led by Dr. David Ebert, Dr. Binbin Weng, and Dr. Chenghao Wang at the University of Oklahoma, the project’s goal was to develop an integrated approach for building this engineering platform to detect, quantify, and mitigate methane emissions across various temporal scale, spatial scales, and sectors. The planning grant study began with an extensive review of various methane sensing and monitoring technologies and systems, surveying over 100 technology providers globally. This review revealed the prevalence of optical methods over chemical methods in commercially available sensors, with Non-Dispersive Infrared (NDIR), Tunable Diode Laser Absorption Spectroscopy (TDLAS), and Optical Gas Imaging (OGI) cameras being the most prevalent options. A trend towards more advanced optical techniques was observed, driven by increased regulatory focus and technological advancements. The technical evaluation of these sensing technologies provided crucial insights into their capabilities and limitations. The study examined emerging technologies such as Differential Absorption LiDAR (DIAL), which show promise for high-precision and long-range detection. The team then investigated the features and application bandwidth of various sensing platforms, including handheld, fixed/stationary, mobile, aerials, and spaceborne monitors. Pilot field studies were conducted to assess the capabilities of solutions for different emission scenarios. Field work with sensor deployments was conducted at three distinct site types: an oil & gas industry site, a cattle ranching operation, and a waste processing facility. The team also conducted a thorough review of methane flux inverse modeling approaches, focused on physically based methods. These approaches were categorized into simple, intermediate, and advanced methods. A realtime WRF-GHG (Weather Research and Forecasting-Greenhouse Gas) modeling system was developed and applied, incorporating multiple data sources to guide field experiments and inform methane plume detection. The project identified and analyzed numerous categories of methane data sources, including satellite measurements, ground-based sensors, and inventory databases. Key platforms examined include EDGAR, EPA GHGI, NASA TROPOMI, Carbon Mapper, and Climate TRACE, among others. The team proposed an architecture for a comprehensive methane monitoring platform. This system incorporates multi-source data acquisition, advanced data processing and assimilation, interactive visualization tools, and analytical capabilities for emissions forecasting and scenario analysis. The proposed platform aims to provide a user-friendly interface catering to various stakeholders, from researchers to policymakers. The architecture includes sophisticated data ingestion methods, a centralized data warehouse, and advanced analytical tools for data fusion and interpretation. To ensure the relevance and effectiveness of the proposed system, a comprehensive survey was conducted to gather stakeholder input on system requirements. Key findings include a strong need for integrating various data types and formats, a preference for real-time data updates and advanced visualization tools, and a demand for user-friendly interfaces catering to different expertise levels.

03 NATURAL GAS↗
DOE OSTI2023

Softening of a flat phonon mode in the kagome ScV 6 Sn 6

Geometrically frustrated kagome lattices are raising as novel platforms to engineer correlated topological electron flat bands that are prominent to electronic instabilities. Here, we demonstrate a phonon softening at the k z = π plane in ScV 6 Sn 6 . The low energy longitudinal phonon collapses at ~98 K and q = $\tfrac{1}{3}$ $\tfrac{1}{3}$ $\tfrac{1}{2}$ due to the electron-phonon interaction, without the emergence of long-range charge order which sets in at a different propagation vector q CDW = $\tfrac{1}{3}$ $\tfrac{1}{3}$ $\tfrac{1}{3}$. Theoretical calculations corroborate the experimental finding to indicate that the leading instability is located at $\tfrac{1}{3}$ $\tfrac{1}{3}$ $\tfrac{1}{2}$ of a rather flat mode. We relate the phonon renormalization to the orbital-resolved susceptibility of the trigonal Sn atoms and explain the approximately flat phonon dispersion. Our data report the first example of the collapse of a kagome bosonic mode and promote the 166 compounds of kagomes as primary candidates to explore correlated flat phonon-topological flat electron physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Optical filters made from random metasurfaces using Bayesian optimization

Abstract We theoretically investigate the ability to design optical filters from a single material and a single layer of randomly dispersed resonant dielectric particles, defining a random metasurface. Using a Bayesian and generalized Mie inverse-design approach, we design particle radii distributions that give rise to longpass, shortpass, bandpass, and bandstop spectral bands in the infrared. The optical response is shown to be directly related to electric and magnetic multipole scattering of the constituent particles and their near field coupling. We discuss the effect of the particle size distribution and particle–particle coupling interactions on filter design in random systems lacking long-range order.

Materials Science↗
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

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↗
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↗
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 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↗