Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dimer”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Systematic Evaluation of Counterpoise Correction in Density Functional Theory

A widespread belief persists that the Boys–Bernardi function counterpoise (CP) procedure “overcorrects” supramolecular interaction energies for the effects of basis-set superposition error. To the extent that this is true for correlated wave function methods, it is usually an artifact of low-quality basis sets. The question has not been considered systematically in the context of density functional theory, however, where basis-set convergence is generally less problematic. We present a systematic assessment of the CP procedure for a representative set of functionals and basis sets, considering both benchmark data sets of small dimers and larger supramolecular complexes. The latter include layered composite polymers with ~150 atoms and ligand–protein models with ~300 atoms. Provided that CP correction is used, we find that intermolecular interaction energies of nearly complete-basis quality can be obtained using only double-ζ basis sets. Furthermore, this is less expensive as compared to triple-ζ basis sets without CP correction. CP-corrected interaction energies are less sensitive to the presence of diffuse basis functions as compared to uncorrected energies, which is important because diffuse functions are expensive and often numerically problematic for large systems. Our results upend the conventional wisdom that CP “overcorrects” for basis-set incompleteness. In small basis sets, CP correction is mandatory in order to demonstrate that the results do not rest on error cancellation.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Protein–Ligand Binding Free-Energy Calculations with ARROW–A Purely First-Principles Parameterized Polarizable Force Field

Protein–ligand binding free-energy calculations using molecular dynamics (MD) simulations have emerged as a powerful tool for in silico drug design. Here, we present results obtained with the ARROW force field (FF)–a multipolar polarizable and physics-based model with all parameters fitted entirely to high-level ab initio quantum mechanical (QM) calculations. ARROW has already proven its ability to determine solvation free energy of arbitrary neutral compounds with unprecedented accuracy. The ARROW FF parameterization is now extended to include coverage of all amino acids including charged groups, allowing molecular simulations of a series of protein–ligand systems and prediction of their relative binding free energies. We ensure adequate sampling by applying a novel technique that is based on coupling the Hamiltonian Replica exchange (HREX) with a conformation reservoir generated via potential softening and nonequilibrium MD. ARROW provides predictions with near chemical accuracy (mean absolute error of ~0.5 kcal/mol) for two of the three protein systems studied here (MCL1 and Thrombin). The third protein system (CDK2) reveals the difficulty in accurately describing dimer interaction energies involving polar and charged species. Overall, for all of the three protein systems studied here, ARROW FF predicts relative binding free energies of ligands with a similar accuracy level as leading nonpolarizable force fields.

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

Stochastic Real-Time Second-Order Green’s Function Theory for Neutral Excitations in Molecules and Nanostructures

We present a real-time second-order Green’s function (GF) method for computing excited states in molecules and nanostructures, with a computational scaling of O(N e 3 ), where N e is the number of electrons. The cubic scaling is achieved by adopting the stochastic resolution of the identity to decouple the 4-index electron repulsion integrals. To improve the time propagation and the spectral resolution, we adopt the dynamic mode decomposition technique and assess the accuracy and efficiency of the combined approach for a chain of hydrogen dimer molecules of different lengths. We find that the stochastic implementation accurately reproduces the deterministic results for the electronic dynamics and excitation energies. Furthermore, we provide a detailed analysis of the statistical errors, bias, and long-time extrapolation. Overall, the approach offers an efficient route to investigate excited states in extended systems with open or closed boundary conditions.

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

How Does HF-DFT Achieve Chemical Accuracy for Water Clusters?

Bolstered by recent calculations of exact functional-driven errors (FEs) and density-driven errors (DEs) of semilocal density functionals in the water dimer binding energy , we investigate approximate FEs and DEs in neutral water clusters containing up to 20 monomers, charged water clusters, and alkali- and halide-water clusters. Our proxy for the exact density is r 2 SCAN 50, a 50% global hybrid of exact exchange with r 2 SCAN, which may be less correct than r 2 SCAN for the compact water monomer but importantly more correct for long-range electron transfers in the noncompact water clusters. We show that SCAN makes substantially larger FEs for neutral water clusters than r 2 SCAN, while both make essentially the same DEs. Unlike the case for barrier heights, these FEs are small in a relative sense and become large in an absolute sense only due to an increase in cluster size. SCAN@HF, short for SCAN evaluated on the Hartree–Fock (HF) density, produces a cancellation of errors that makes it chemically accurate for predicting the absolute binding energies of water clusters. Likewise, adding a long-range dispersion correction to r 2 SCAN@HF, as in the composite method HF-r 2 SCAN-DC4, makes its FE more negative than in r 2 SCAN@HF, permitting a near-perfect cancellation of FE and DE. r 2 SCAN by itself (and even more so, r 2 SCAN evaluated on the r 2 SCAN 50 density), is almost perfect for the energy differences between water hexamers, and thus probably also for liquid water away from the boiling point. Thus, the accuracy of composite methods like SCAN@HF and HF-r 2 SCAN-DC4 is not due to the HF density being closer to the exact density, but to a compensation of errors from its greater degree of localization. We also give an argument for the approximate reliability of this unconventional error cancellation for diverse molecular properties. Lastly, we confirm this unconventional error cancellation for the SCAN description of the water trimer via Kohn–Sham inversion of the CCSD(T) density.

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

Polynomial Scaling Localized Active Space Unitary Selective Coupled Cluster Singles and Doubles

We present a polynomial-scaling algorithm for the localized active space unitary selective coupled cluster singles and doubles (LAS-USCCSD) method. In this approach, cluster excitations are selected based on a threshold ϵ determined by the absolute gradients of the LAS-UCCSD energy with respect to cluster amplitudes. Using the generalized Wick’s theorem for multireference wave functions, we derive the gradient expression as a polynomial function of one-, two-, and three-body reduced density matrices and 1- and 2-electron integrals, valid for any multireference wave function. The resulting gradient implementation exhibits a memory scaling of 𝒪(N 6 ), with N spin orbitals in the combined active space of all fragments. The variational quantum eigensolver is used to optimize the selected cluster excitations on a quantum simulator. Furthermore, by plotting the energy error, defined as the difference between the LAS-USCCSD and corresponding CASCI energies, against the inverse cluster amplitude selection threshold (ϵ –1 ) for polyene chains containing 2 to 5 π-bond units, we establish a relationship between the energy error and the threshold. To further validate the accuracy of LAS-USCCSD, we computed the cis–trans isomerization energy of stilbene (a 20-qubit system) and the magnetic coupling constant of the tris-hydroxo-bridged chromium dimer [Cr 2 (OH) 3 (NH 3 ) 6 ] 3+ (evaluated as both 12- and 20-qubit systems) using the Qiskit-Qulacs simulator. Assessing such examples is important to determine the practical feasibility of quantum simulations for chemically realistic systems. Toward this goal, with the LAS-USCCSD algorithm we estimated the quantum resources required for simulating an active space of (30e,22o) in [Cr 2 (OH) 3 (NH 3 ) 6 ] 3+ , a size that remains beyond the reach of current quantum simulators for accurate treatment.

Algorithms↗
DOE OSTI2025

Reducing the Cost of Energy Differences in Variational Monte Carlo with Spotlight Sampling

Here, we investigate an approximate sampling scheme that can significantly reduce the cost scaling of variational Monte Carlo when it is employed to predict the energy differences associated with local chemical changes. Inspired by side-chaining and embedding methods, this spotlight sampling approach adopts an approximate fragmented Hamiltonian and correlated sampling to reduce cost scaling to the point that it is essentially linear with system size, with the potential to go sublinear if certain conditions are met. In tests on bond stretching energies in alcohols, hydrogen dimer chains, and molecules with various degrees of π-system delocalization, we observe the anticipated linear scaling and an explicit cost crossover with standard variational Monte Carlo.

Bumann, Sonja [University of California, Berkeley,↗
DOE OSTI2025

Overset-Grid Method with Smooth Orbital Partitioning for Molecular Scattering Calculations

To solve molecular photoionization and electron scattering problems, we use an overset-grid representation of electronic continuum functions, which has an extended central spherical grid that overlaps small spherical grids (subgrids) centered on each atom of a polyatomic molecule. Here, in this work, we present an improved algorithm that smoothly partitions the total wave function between the central grid and the atomic subgrids. The smooth partitioning allows one to use approximately one-fourth the number of partial waves on the central grid compared to our previous implementation with switching functions. The resulting numerical method for treating electron scattering and photoionization of polyatomic molecules combines the accuracy and flexibility of pure numerical grid representations with the rapid convergence of hybrid combinations of atom-centered basis-set expansions and grid methods. The overset-grid representation is implemented using the complex Kohn variational principle for scattering and photoionization amplitudes. The faster convergence with respect to the number of central grid partial waves is demonstrated and accuracy is verified by comparisons with the previous implementation and with far more computationally demanding single-center numerical expansions in electron-molecule scattering and photoionization calculations on the neon dimer (Ne 2 ) system, carbon tetrafluoride (CF 4 ) molecule, and the pyridine (C 5 H 5 N) molecule in the static-exchange approximation.

Molecules↗
DOE OSTI2025

A Universal Augmentation Framework for Long-Range Electrostatics in Machine Learning Interatomic Potentials

Most current machine learning interatomic potentials (MLIPs) rely on short-range approximations, without explicit treatment of long-range electrostatics. To address this, we recently developed the Latent Ewald Summation (LES) method, which infers electrostatic interactions, polarization, and Born effective charges (BECs), just by learning from energy and force training data. Here, in this study, we present LES as a standalone library, compatible with any short-range MLIP, and demonstrate its integration with methods such as MACE, NequIP, Allegro, CACE, CHGNet, and UMA. We benchmark LES-enhanced models on distinct systems, including bulk water, polar dipeptides, and gold dimer adsorption on defective substrates, and show that LES not only captures correct electrostatics but also improves accuracy. Additionally, we scale LES to large and chemically diverse data by training MACELES-OFF on the SPICE set containing molecules and clusters, making a universal MLIP with electrostatics for organic systems, including biomolecules. MACELES-OFF is more accurate than its short-range counterpart (MACE-OFF) trained on the same data set, predicts dipoles and BECs reliably, and has better descriptions of bulk liquids. By enabling efficient long-range electrostatics without directly training on electrical properties, LES paves the way for electrostatic foundation MLIPs.

Kim, Dongjin [University of California, Berkeley, ↗
DOE OSTI2026

Local Pair Natural Orbital-Based Coupled-Cluster Theory through Full Quadruples (DLPNO–CCSDTQ)

In this work, we implement a local pair natural orbitalbased coupled-cluster method through the full treatment of quadruple excitations (CCSDTQ). The domain-based local pair natural orbital (DLPNO) approach, which has successfully been applied to lower levels of coupled-cluster theory, is utilized in our algorithm, and thus our algorithm is called DLPNO-CCSDTQ. For simplicity in the working equations and in the implementation, we t 1 -dress the twoelectron integrals as well as Fock matrix elements. Our method can recover CCSDTQ-CCSDT and CCSDTQ-CCSDT(Q) energy differences on the order of 0.01−0.05 kcal mol −1 , even at a loose quadruples natural orbital (QNO) occupation number cutoff of 3.33 × 10 −6 . To highlight the capabilities of our code and its potential future applications, we showcase computations that would be intractable with canonical CCSDTQ, such as the benzene dimer, (H 2 O) 17 , and adamantane. With sufficient computing resources, computations up to 15 heavy atoms (40 atoms overall) may be feasible for fully bonded 3D systems.

Cluster chemistry↗
DOE OSTI2026

Linear-Scaling Asymmetric Triples Correction through the Solution of the DLPNO–CCSD Lambda Equations: DLPNO–CCSD(T) Λ

In this research, we derive equations for solving for the stationary points of the DLPNO–CCSD Lagrangian, in the t 1 -transformed formalism introduced earlier and as currently implemented in the P SI 4 quantum chemistry software package. These lambda equations in the local pair natural orbital basis allow for the evaluation of CCSD(T) Λ energetics with linear-scaling computational effort, also known as the asymmetric triples correction. This DLPNO–CCSD(T) Λ method allows for accurate triples contributions to be computed for larger molecules, especially in cases that CCSD(T) is known to be insufficient, such as with multireference systems and bond-breaking systems. We showcase the accuracy of our code on reaction energies, barrier heights, and noncovalent interaction energies. Also showcased are the capabilities of our code by evaluating DLPNO–CCSD(T) Λ energetics on large noncovalent dimers up to 112 atoms, as well as a rhodium catalyst complex containing 66 atoms.

Cluster chemistry↗
DOE OSTI2020

Data-Driven Many-Body Models for Molecular Fluids: CO 2 /H 2 O Mixtures as a Case Study

In this study, we extend the scope of the many-body TTM-nrg and MB-nrg potential energy functions (PEFs), originally introduced for halide ion–water and alkali-metal ion–water interactions, to the modeling of carbon dioxide (CO2) and water (H2O) mixtures as prototypical examples of molecular fluids. Both TTM-nrg and MB-nrg PEFs are derived entirely from electronic structure data obtained at the coupled cluster level of theory and are, by construction, compatible with MB-pol, a many-body PEF that has been shown to accurately reproduce the properties of water. Although both TTM-nrg and MB-nrg PEFs adopt the same functional forms for describing permanent electrostatics, polarization, and dispersion, they differ in the representation of short-range contributions, with the TTM-nrg PEFs relying on conventional Born–Mayer expressions and the MB-nrg PEFs employing multidimensional permutationally invariant polynomials. By providing a physically correct description of many-body effects at both short and long ranges, the MB-nrg PEFs are shown to quantitatively represent the global potential energy surfaces of the CO 2 –CO 2 and CO 2 –H 2 O dimers and the energetics of small clusters, as well as to correctly reproduce various properties in both gas and liquid phases. Furthermore, building upon previous studies of aqueous systems, our analysis provides further evidence for the accuracy and efficiency of the MB-nrg framework in representing molecular interactions in fluid mixtures at different temperature and pressure conditions.

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

Competition between Hückel’s Rule and Jahn–Teller Distortion in Small Carbon Rings: A Quantum Monte Carlo Study

Structural properties and energetics of carbon rings are studied with the diffusion Monte Carlo (DMC) method. Our DMC-based geometry optimization reveals that both polyynic C 4n and cumulenic C 4n+2 rings exhibit bond length alternations for n >= 3, which is understood to be due to Jahn-Teller distortions. The bond length alternation even in a cumulenic (4n + 2) carbon ring was experimentally observed in a recently synthesized C 18 molecule. From a comparison of the DMC cohesive energies of C 4n with those of C 4n + 2 , we present a comprehensive picture of the competition between Huckel's rule and Jahn-Teller distortion in small carbon rings; the former is more dominant than the latter for n < 5 where C 4n + 2 rings are more stable than C 4n , while C 4n rings are as stable as C 4n + 2 for n < 5 where dimerization effects due to Jahn-Teller distortion are more important.

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

Electron-Beam-Induced Molecular Plasmon Excitation and Energy Transfer in Silver Molecular Nanowires

In this work, we investigate (1) electron-beam-induced plasmon absorption spectra of Ag molecular nanowire dimers and (2) electron-beam-induced energy transfer between two nanowires. We employ linear-response time-dependent density functional theory (TDDFT) and real-time TDDFT methods to simulate the electron-beam-induced plasmonic excitations, dynamics, and corresponding electron energy loss spectrum for small models of a single molecular nanowire with four Ag atoms and for two Ag nanowires. An array of different relative orientations of nanowires and of different initial excitation conditions resulting from applying an electron beam at different positions with respect to the Ag nanowires is investigated. The results demonstrate (1) an electron beam can induce plasmonic excitations from the molecular Ag nanowire ground state to the excited states that are both optically allowed and forbidden, (2) a tunability for selective excitations that can be controlled by the position of a focused electron beam, and (3) kinetic and dynamic behaviors of time-dependent electron-beam-induced energy transfer between two Ag molecular nanowires depend on the position of the beam source and nanowire separation distance, providing insights into the spatial dependences of plasmonic couplings in nanowire arrays.

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

Bond Dissociation Energies in Heavy Element Chalcogen and Halogen Small Molecules

Thermodynamic properties including bond dissociation energies, heats of formation, and gas phase acidities for the hydrides and dimers of chalcogens and halogens: H 2 Y, HX, Y 2 and X 2 for Y = Se, Te, At and X = Br, I, and At have been predicted using the Feller-Peterson-Dixon composite correlated molecular orbital theory approach. A full 4-component CCSD(T) approach was used to calculate the spin-orbit effects on thermodynamic properties, except for Se 2 where the AoC-DHF value was used due to strong multireference effects in Se 2 for the SO calculations. The calculated results show that the At 2 bond dissociation energy is quite small, 19.5 kcal/mol, with much of the low bond energy due to spin orbit effects. H 2 Po is not predicted to be stable to dehydrogenation to Po + H 2 in terms of the free energy at 298 K. In the gas phase, HAt is predicted to be a stronger acid than H 2 SO 4 . Here, the current results provide insights into potential difficulties in the actual experimental observation of such species for heavy elements.

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

Interaction–Deletion: A Composite Energy Method for the Optimization of Molecular Systems Selectively Removing Specific Nonbonded Interactions

The complex interactions between different portions of a large molecule can be challenging to analyze through traditional electronic structure calculations. Moreover, standard methods cannot easily quantify the physical consequences of individual pairwise interactions inside a molecule. By creating a set of molecular fragments, we propose a composite energy method to explore changes in a molecule caused by removing selected nonbonded interactions between different molecular portions. Energies and forces are easily obtained with this composite approach, allowing geometry optimizations that lead to chemically meaningful structures that describe how the omitted interactions contribute to changes in the local geometrical minima. We illustrate the application of our new hybrid scheme by computing the influence of intramolecular hydrogen-bonding interactions in two small molecules: 1,6-(tG + G + TG + G + g – )-hexanediol and a cyclic analogue, cis-1,4-cyclohexanediol. The resulting structural and energetic changes are interpreted to yield key physical insights and quantify concepts such as “preparation energy” or “reorganization energy”. We demonstrate that the composite method can be extended to larger molecular systems by showing its application on a Si(100) surface model containing interactions between dissociated ammonia molecules on adjacent surface dimers. The scheme’s efficacy is also tested by applying it to systems having multiple intramolecular interactions, viz., 310-polyglycine and H+GPGG. Furthermore, the cooperative nature of intramolecular hydrogen bonds is explored by using interaction–deletion in 2-nitrobenzene-1,3-diol.

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

Pulse Radiolysis Investigation of Radicals Derived from Water-Soluble Cyanine Dyes: Implications for Super-resolution Microscopy

Light-induced blinking, an inherent feature of many forms of super-resolution microscopy, has been linked to transient reduction of the fluorescent cyanine dye used as an imaging agent. There is, however, only scant literature information related to one-electron reduced cyanine dyes, especially in an aqueous environment. Here, we examine a small series of cyanine dyes, possessing disparate π-conjugation lengths, under selective reducing or oxidizing conditions. The experiment allows recording of both differential absorption spectra and decay kinetics of the resultant one-electron reduced or oxidized transient species in water. Relative to the ground state, absorption transitions for the various radicals are weak and somewhat broadened but do allow correlation with the π-conjugation length. In all cases, absorption maxima lie to the blue of the main ground-state transition. Under anaerobic conditions, the transient species decay on the microsecond to millisecond time scale, with the mean lifetime depending on molecular structure, radiation dose, and dye concentration. The experimental absorption spectra recorded for the one-electron reduced radicals and the presumed dimer cation radical compare well to spectra obtained from time-dependent density functional theory calculations. Furthermore, the results allow conclusions to be drawn regarding the plausibility of the reduced species being responsible for light-induced blinking in direct stochastic optical reconstruction microscopy.

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

A Nonorthogonal Configuration Interaction Approach to Singlet Fission in Perylenediimide Compounds

Perylenediimide molecules constitute a family of chromophores that undergo singlet fission, a process in which an excited singlet state converts into lower energy triplets on two neighboring molecules, potentially increasing the efficiency of organic solar cells. Here, the nonorthogonal configuration interaction method is applied to study the effect of the different crystal packing of various perylenediimide derivatives on the relative energies of the singlet and triplet states, the intermolecular electronic couplings, and the relative rates for singlet fission. The analysis of the wave functions and electronic couplings reveals that charge transfer states play an important role in the singlet fission mechanism. Dimer conformations where the PDI molecules are at large displacements along the long axis and short on the short axis are posed as the most favorable for singlet fission. The role of the substituent at the imide group has been inspected concluding that, although it has no effect in the energies, for some conformations it significantly influences the electronic couplings, and therefore, replacing this substituent with hydrogen may introduce artifacts in the computational modeling of the PDI molecules.

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

Excited-State Charge Transfer in Push–Pull Platinum(II) π-Extended Porphyrins Fused with Pentacenequinone

Platinum(II) π-extended porphyrins fused with pentacenequinone and dihydropentacene have been successfully synthesized. These porphyrins were investigated using various techniques including absorption, steady-state, and time-resolved phosphorescence spectroscopy and differential pulse voltammetry. UV–vis absorption spectra of pentacenequinone-fused porphyrins (SW-Pt1 and SW-Pt2) showed unusually broad and nontypical absorption patterns. Phosphorescence spectra of SW-Pt1, SW-Pt2, and SW-Pt3 displayed similar emissions in the 704–706 nm region indicating electronic transitions of similar origin; however, the triplet lifetimes were found to be quenched in the case of both SW-Pt1 and SW-Pt2, suggesting the occurrence of excited-state events. Facile reductions were obtained for both the pentacene-quinone-fused monomer (SW-Pt2) and dimer (SW-Pt1) and were identified to be located at the pentacenequinone components. Here, the observed orbital segregations for SW-Pt2 and SW-Pt1 from DFT calculations supported the possibility of charge transfer in these push–pull systems. Interestingly, the established energy level diagram revealed that the charge transfer from the triplet excited Pt porphyrin is thermodynamically an uphill process. Systematic studies involving both femtosecond and nanosecond transient absorption techniques revealed that the singlet excited Pt porphyrins undergo an intermediate charge transfer state prior to populating the triplet state, providing a plausible explanation for phosphorescence quenching. The lifetime of the intermediate charge transfer states was found to be 25.9 and 5.68 ps, respectively, for SW-Pt1 and SW-Pt2.

Charge transfer↗