Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Delocalization”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Understanding of three different polyvinylpyrrolidone (PVP) based battery binders blends on graphene surfaces from first principles via DFT simulations

Binders play a crucial role in binding electrodes, current collectors, and conductive agents during battery manufacturing. The selection of the binders and their blends considerably affect various properties such as mechanical durability, adhesion performance, ionic/electronic conductivities, and solid electrolyte interface stability for lithium-ion batteries. In this work, polyvinylidene di-fluoride (PVDF), polyacrylic acid (PAA), polyvinylpyrrolidone (PVP), and lithiated polyacrylic acid (Li-PAA) were selected as the model binders, and their blends were studied on graphene surfaces. Further, among them, PVDF is hardly soluble in water unlike the other binders, so in practice, the blend of PVP:PVDF is not thought to be a proper combination. However, the PVP:PVDF binder has also been included because the purpose of this work is to calculate several electronic properties such as binding energies, intermolecular interaction energies, bond critical points, electron density, and the Laplacian of electron density at critical bonding regions of all these blends at the molecular level via density functional theory (DFT) simulations in order to evaluate and compare how is the interaction strength and bonding type of three different lithium-ion batteries (LIBs) binder blends on graphene surfaces. Most stable binder pairs and their binding mechanisms on graphene surfaces were studied. The delocalization of lithium in Li-PAA was studied with the presence of water.

25 ENERGY STORAGE↗
DOE OSTI2023

E-beam-enhanced solid-state mechanical amorphization of α-quartz: Reduced deformation barrier via localized excess electrons as network modifiers

Under hydrostatic pressure, α-quartz (α-SiO 2 ) undergoes solid-state mechanical amorphization wherein the interpenetration of [SiO 4 ] 2- tetrahedra occurs and the material loses crystallinity. This phase transformation requires a high hydrostatic pressure of 14GPa because the repulsive forces resulting from the ionic nature of the Si–O bonds prevent the severe distortion of the atomic configuration. Herein, we experimentally and computationally demonstrate that e-beam irradiation changes the nature of the interatomic bonds in α-quartz and enhances the solid-state mechanical amorphization at nanoscale. Specifically, during in situ uniaxial compression, a larger permanent deformation occurs in α-quartz submicron pillars compressed during e-beam irradiation than in those without e-beam irradiation. Microstructural analysis reveals that the large permanent deformation under e-beam irradiation originates from the enhanced mechanical amorphization of α-quartz and the subsequent viscoplastic deformation of the amorphized region. Further, atomic-scale simulations suggest that the delocalized excess electrons introduced by e-beam irradiation move to highly distorted atomic configurations and alleviate the repulsive force, thus reducing the barrier to the solid-state mechanical amorphization. In conclusion, these findings deepen our understanding of electron–matter interactions and can be extended to new glass forming and processing technologies at nano- and microscale.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Laser powder bed fusion of oxide dispersion-strengthened IN718 alloys: A complementary study on microstructure and mechanical properties

In this study, two new grades of oxide dispersion strengthened (ODS) Inconel 718 (IN718) alloys were designed by the thermochemical CALPHAD method and produced by laser powder bed fusion (LPBF) technique. Alloys designated as IN718-YF and IN718-YFH, that consist Y 2 O 3 –FeO and Y 2 O 3 –FeO–Hf, respectively, were fabricated with >99.9 % densification using optimized process parameters. CALPHAD calculations were highly consistent with experimental findings, highlighting the formation of Al-containing Y–Ti–O and Y–Hf–O nano-oxides in both alloy types. Texture analyzes revealed no significant texture development in as-built (AB) or heat-treated (HT) alloys. Heat treatment was applied at 1050 °C for 1 h to enhance nano-oxide density. Further, the nano-oxide number density remained similar in IN718-YF while it decreased in IN718-YFH alloy as a result of carbide formation after the heat treatment. Besides, formation of secondary γ' particles was observed in the IN718-YFH/HT alloy. Even though the yield strengths of IN718-YF and IN718-YFH alloys in both AB and HT conditions were similar, the ductility of IN718-YFH was ~50 % less in almost all conditions compared to the ductility of IN718-YF. This has been shown to be as a result of irregular shaped micron-sized Y-Hf-O oxides, martensite formation in AB condition, increased amount of carbides and existence of secondary γ' particles in HT condition in IN718-YFH. High density of stacking faults (SF) forming at the interface of the nano-oxides have been detected in IN718-YF alloys. Besides dislocation/nanoparticle interactions, SFs which are responsible for the delocalization of the deformation improve the ductility of IN718-YF alloys. Overall, high temperature mechanical tests exhibit that both alloys have higher strength with improved ductility compared to the standard IN718 alloys, indicating the contribution of the nano-oxides.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Bay capping via acetylene addition to polycyclic aromatic hydrocarbons: Mechanism and kinetics

The bay-capping mechanism on PAH armchair edges and the kinetics of acetylene addition to 6–6–5 and 5–6–5 bays have been explored by ab initio/RRKM-ME calculations. The bays on the edges were modeled by C 21 H 11 and C 20 H 9 radicals produced by H abstractions from 7H-benzo[c]cyclopenta[e]pyrene and dicyclopenta[cf]pyrene. The C 20 H 9 + C 2 H 2 reaction is shown to have a low entrance barrier and to rapidly form the capped product, indaceno[2,1,8,7-cdefg]pyrene, along with ethynyl substituted dicyclopenta[cf]pyrene at temperatures above 1400 K. The reactivity of C 21 H 11 is shown to be governed by the location of the unpaired electron; the π radical R1 formed by H abstraction from the CH 2 group in 7H-benzo[c]cyclopenta[e]pyrene reacts with C 2 H 2 very slowly owing to a high entrance barrier, with the bay-capping rate constant approaching 10 -16 cm 3 molecule -1 s -1 only at temperatures above 2000K. Further, this result reaffirms that the growth of π aryl radicals via acetylene addition is inefficient and reflects the generally low reactivity of such radicals where the spin density is highly delocalized over the entire polyaromatic system. Alternatively, the σ C 21 H 11 radical R2 produced by H abstraction from the five-membered ring at the bay rapidly reacts with C 2 H 2 forming the bay-capped product, with the rate constant on the order of 10 -12 cm 3 molecule -1 s -1 at T≥1500K. Rate constants for the capping reactions at the 6–6–5 and 5–6–5 bays are compared with those at the 6–0–6, 6–6–6, and 6–5–6 bays. The site-specific bay-capping rate constants have been utilized in kMC simulations of the PAH growth and the results showed measurable differences when the 6–6–5 and 5–6–5 bay-capping reactions are taken into account, including an increase of the growth rate and the formation of closed-shell PAH and a rise of the number of embedded five-membered rings accompanied with a slight decrease of their overall amount.

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

Tuning the mechanical properties of nanoglass-metallic glass composites with brick and mortar designs

We use molecular dynamics simulations to characterize the tensile deformation and failure of metallic glass (MG)-nanoglass (NG) composites designed with a brick and mortar architecture consisting of a 3 nm-grain-size NG matrix and a MG second phase. Results show that the overall strength of the matrix is significantly improved by the second phase without sacrificing the large NG matrix ductility. The mechanical synergy of the two phases is optimized by arranging the MG bricks in a staggered way, which effectively delocalizes the plastic deformation, hindering the buildup of local stress hot spots and the generation of critical shear bands.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Local structural ordering determines the mechanical damage tolerance of amorphous grain boundary complexions

Amorphous grain boundary complexions act as toughening features within a microstructure because they can absorb dislocations more efficiently than traditional grain boundaries. This toughening effect should be a strong function of the local internal structure of the complexion, which has recently been shown to be determined by grain boundary crystallography. To test this hypothesis, molecular dynamics are used here to simulate dislocation absorption and damage nucleation for complexions with different distributions of structural short-range order. Here, the complexion with a more disordered structure away from the dislocation absorption site is actually found to better resist crack nucleation, as damage tolerance requires delocalized deformation and the operation of shear-transformation zones through the complexion thickness. The more damage tolerant complexion accommodates plastic strain efficiently within the entire complexion, providing the key mechanistic insight that local patterning and asymmetry of structural short-range order controls the toughening effect of amorphous complexions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Detecting Magnon-phonon coupling in yttrium iron garnet with variable temperature STEM-EELS

Magnons, quanta of spin wave excitations in magnetically ordered materials, have been identified as candidates for several potentially transformative technologies in recent years. Macroscopic techniques, such as neutron scattering or Raman spectroscopy, can be used to identify and analyze magnons, but provide relatively delocalized information about the sample. Understanding how the bonding and local structure of a material interacts with, and influences, the magnon population in a material is a crucial step toward the ability to produce any real-world application utilizing magnons. By leveraging the combined spatial resolution of scanning transmission electron microscopy (STEM) and the energy resolution of monochromated electron energy-loss spectroscopy (EELS) nanoscale analysis of magnons can be performed. While the weak interaction of magnons with the electron beam makes magnon EELS challenging on reasonable timescales, magnon-phonon coupling can be leveraged to understand magnons through their effect on the more easily measured phonons. Here, we examine yttrium iron garnet (YIG) flakes, and demonstrate non-linear, temperature-dependent shifts in the phonon frequencies, consistent with previously described magnon-phonon coupling effects. The ability to measure the temperature-dependence of vibrational frequencies with high precision in individual nanoscale flakes, demonstrates the ability to study magnon-phonon coupling in the STEM with unprecedented spatial resolution.

Reifsnyder, Alexander [ORNL]↗
DOE OSTI2024

Development of (NO)Fe(N 2 S 2 ) as a Metallodithiolate Spin Probe Ligand: A Case Study Approach

The ubiquity of sulfur–metal connections in nature inspires the design of bi- and multimetallic systems in synthetic inorganic chemistry. Common motifs for biocatalysts developed in evolutionary biology include the placement of metals in close proximity with flexible sulfur bridges as well as the presence of π-acidic/delocalizing ligands. This Account will delve into the development of a (NO)Fe(N 2 S 2 ) metallodithiolate ligand that harnesses these principles. The Fe(NO) unit is the centroid of a N 2 S 2 donor field, which as a whole is capable of serving as a redox-active, bidentate S-donor ligand. Its paramagnetism as well as the ν(NO) vibrational monitor can be exploited in the development of new classes of heterobimetallic complexes. We offer four examples in which the unpaired electron on the {Fe(NO)} 7 unit is spin-paired with adjacent paramagnets in proximal and distal positions. First, the exceptional stability of the (NO)Fe(N 2 S 2 )-Fe(NO) 2 platform, which permits its isolation and structural characterization at three distinct redox levels, is linked to the charge delocalization occurring on both the Fe(NO) and the Fe(NO) 2 supports. This accommodates the formation of a rare nonheme {Fe(NO)} 8 triplet state, with a linear configuration. A subsequent FeNi complex, featuring redox-active ligands on both metals (NO on iron and dithiolene on nickel), displayed unexpected physical properties. Our research showed good reversibility in two redox processes, allowing isolation in reduced and oxidized forms. Various spectroscopic and crystallographic analyses confirmed these states, and Mössbauer data supported the redox change at the iron site upon reduction. Oxidation of the complex produced a dimeric dication, revealing an intriguing magnetic behavior. The monomer appears as a spin-coupled diradical between {Fe(NO)} 7 and the nickel dithiolene monoradical, while dimerization couples the latter radical units via a Ni 2 S 2 rhomb. Magnetic data (SQUID) on the dimer dication found a singlet ground state with a thermally accessible triplet state that is responsible for magnetism. A theoretical model built on an H 4 chain explains this unexpected ferromagnetic low-energy triplet state arising from the antiferromagnetic coupling of a four-radical molecular conglomerate. For comparison, two (NO)Fe(N 2 S 2 ) were connected through diamagnetic group 10 cations producing diradical trimetallic complexes. Antiferromagnetic coupling is observed between {Fe(NO)} 7 units, with exchange coupling constants (J) of -3, -23, and -124 cm –1 for Ni II , Pd II , and Pt II , respectively. This trend is explained by the enhanced covalency and polarizability of sulfur-dense metallodithiolate ligands. A central paramagnetic trans-Cr(NO)(MeCN) receiver unit core results in a cissoid structural topology, influenced by the stereoactivity of the lone pair(s) on the sulfur donors. This {Cr(NO)} 5 radical bridge, unlike all previous cases, finds the coupling between the distal Fe(NO) radicals to be ferromagnetic (J = 24 cm –1 ). The stability and predictability of this S = 1/2 moiety and the steric/electronic properties of the bridging thiolate sulfurs suggest it to be a likely candidate for the development of novel molecular (magnetic) compounds and possibly materials. The role of synthetic inorganic chemistry in designing synthons that permit connections of the (NO)Fe(N 2 S 2 ) metalloligand is highlighted as well as the properties of the heterobi- and polymetallic complexes derived therefrom.

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

Insights into Copper Sulfide Formation from Cu and S K edge XAS and DFT studies

An understanding of the fundamentals of the reaction between CuO with trace amounts of H 2 S to form CuS products is critical for the optimal utilization of this process in sulfur removal applications. Unfortunately, CuS is a complex material, featuring various Cu 2-x S compounds (with 0 ≤ x ≤ 1), distorted crystal phases, and varying electronic structures and coordination environments of Cu and S ions. In this work, we combine ex situ and in situ X-ray absorption spectroscopy (XAS) at S and Cu K edges, fixed bed sorption experiments, DFT simulations, and other characterization techniques to speciate the CuS products formed at different temperatures (298–383 K) and from CuO sorbents with different crystallite sizes (2.8–40 nm). The results of our analysis identify the formation of a distorted CuS layer at the surface of CuO crystals with disulfide groups with shorter Cu–S bonds and higher delocalization of the positive charge of the Cu center into (S 1– )2. This distorted CuS layer dominates the XAS signal at lower temperatures (298–323 K) and at the initial stages of sulfidation at higher temperatures (353 and 383 K) where conversion is low (<40%). First-principles atomistic simulations confirm the thermodynamic favorability of the formation of surface (S 1– )2 on both CuO (111) and ($\bar{1}11$) surfaces, providing further support for our experimental observations. Furthermore, these simulations reveal that the presence of disulfide bonds stabilized surface hydroxyl groups, leading to lower Gibbs Free Energies of their surface migration.

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

Methane Generation from CO 2 with a Molecular Rhenium Catalyst

The atomic-level tunability of molecular structures is a compelling reason to develop homogeneous catalysts for challenging reactions such as the electrochemical reduction of carbon dioxide to valuable C 1 –C n products. Of particular interest is methane, the largest component of natural gas. Herein, we report a series of three isomeric rhenium tricarbonyl complexes coordinated by the asymmetric diimine ligands 2-(isoquinolin-1-yl)-4,5-dihydrooxazole ( quin-1-oxa ), 2-(quinolin-2-yl)-4,5-dihydrooxazole ( quin-2-oxa ), and 2-(isoquinolin-3-yl)-4,5-dihydrooxazole ( quin-3-oxa ) that catalyze the reduction of CO 2 to carbon monoxide and methane, albeit the latter with a low efficiency. To our knowledge, these complexes are the first examples of rhenium(I) catalysts capable of converting carbon dioxide into methane. Re(quin-1-oxa)(CO) 3 Cl ( 1 ), Re(quin-2-oxa)(CO) 3 Cl ( 2 ), and Re(quin-3-oxa)(CO) 3 Cl ( 3 ) were characterized and studied using a variety of electrochemical and spectroscopic techniques. In bulk electrolysis experiments, the three complexes reduce CO 2 to CO and CH 4 . When the controlled-potential electrolysis experiments are performed at -2.5 V (vs Fc +/0 ) and in the presence of the Brønsted acid 2,2,2-trifluoroethanol, methane is produced with turnover numbers that range from 1.3 to 1.8. Isotope labeling experiments using 13 CO 2 atmosphere produce 13 CH 4 ( m / z = 17) confirming that methane originates from CO 2 reduction. Theoretical calculations are performed to investigate the mechanistic aspects of the 8e – /8H + reduction of CO 2 to CH 4 . Overall, a ligand-assisted pathway is proposed to be an efficient pathway in the formation of CH 4 . Delocalization of the electron density on the (iso)quinoline moiety upon reduction stabilizes the key carbonyl intermediate leading to additional reactivity of this ligand. These results should aid the development of more robust catalytic systems that produce CH 4 from CO 2 .

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

Origin of Ferromagnetic Exchange Coupling in Donor–Acceptor Biradical Analogues of Charge-Separated Excited States

A new donor–acceptor biradical complex, Tp Cum,Me Zn(SQ-VD) (Tp Cum,Me Zn + = zinc(II) hydro-tris(3-cumenyl-5-methylpyrazolyl)borate complex cation; SQ = orthosemiquinone; VD = oxoverdazyl), which is a ground-state analogue of a charge-separated excited state, has been synthesized and structurally characterized. Here, the magnetic exchange interaction between the S = 1/2 SQ and the S = 1/2 VD within the SQ-VD biradical ligand is observed to be ferromagnetic, with $J_{SQ-VD}$ = +77 cm –1 (H = -2$J_{SQ-VD}\bar{S}_{SQ}·\bar{S}_{VD}$) determined from an analysis of the variable-temperature magnetic susceptibility data. The pairwise biradical exchange interaction in Tp Cum,Me Zn(SQ-VD) can be compared with that of the related donor–acceptor biradical complex Tp Cum,Me Zn(SQ-NN) (NN = nitronyl nitroxide, S = 1/2), where $J_{SQ-NN}$ ≅ +550 cm –1 . This represents a dramatic reduction in the biradical exchange by a factor of ~7, despite the isolobal nature of the VD and NN acceptor radical SOMOs. Computations assessing the magnitude of the exchange were performed using a broken-symmetry density functional theory (DFT) approach. These computations are in good agreement with those computed at the CASSCF NEVPT2 level, which also reveals an S = 1 triplet ground state as observed in the magnetic susceptibility measurements. A combination of electronic absorption spectroscopy and CASSCF computations has been used to elucidate the electronic origin of the large difference in the magnitude of the biradical exchange coupling between Tp Cum,Me Zn(SQ-VD) and Tp Cum,Me Zn(SQ-NN). A Valence Bond Configuration Interaction (VBCI) model was previously employed to highlight the importance of mixing an SQ SOMO → NNLUMO charge transfer configuration into the electronic ground state to facilitate the stabilization of the high-spin triplet (S = 1) ground state in Tp Cum,Me Zn(SQ-NN). Here, CASSCF computations confirm the importance of mixing the pendant radical (e.g., VD, NN) LUMO (VD LUMO and NN LUMO ) with the SOMO of the SQ radical (SQ SOMO ) for stabilizing the triplet, in addition to spin polarization and charge transfer contributions to the exchange. An important electronic structure difference between Tp Cum,Me Zn(SQ-VD) and Tp Cum,Me Zn(SQ-NN), which leads to their different exchange couplings, is the reduced admixture of excited states that promote ferromagnetic exchange into the Tp Cum,Me Zn(SQ-VD) ground state, and the intrinsically weaker mixing between the VD LUMO and the SQ SOMO compared to that observed for Tp Cum,Me Zn(SQ-NN), where this orbital mixing is significant. The results of this comparative study contribute to a greater understanding of biradical exchange interactions, which are important to our understanding of excited-state singlet–triplet energy gaps, electron delocalization, and the generation of electron spin polarization in both the ground and excited states of (bpy)Pt(CAT-radical) complexes.

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

Stabilization of [(N 5 ) 2 BX] 2– and [(N 5 ) 2 B 2 X 2 ] 2– (X = H, F, Cl, Br) by Conjugation and Hyperconjugation Effects

The isolation of nucleophilic boron bases has led to a paradigm shift in boron chemistry. Previous studies of the bis(carbene) borylene complexes revealed that these compounds possess strong donor abilities, and their reaction inertness is due to the large steric hindrance between boron reagents and reactant. In the present study, we have theoretically studied the [(N 5 ) 2 BX] 2– and [(N 5 ) 2 B 2 X 2 ] 2– compounds (X = H, F, Cl, Br). Their electronic structures and properties are discussed by using the NBO, LOL, and ELF methods. We found that both π-conjugation and hyperconjugation effects can effectively stabilize the substituted nucleophilic anionic boron compounds [(N 5 ) 2 BX] 2– and [(N 5 ) 2 B 2 X 2 ] 2– . Substituents, especially X = H, stabilize the boron center through highly delocalized π-bonding, involving the formally “empty” in-plane p orbitals of the boron atom. While the halogen substituents have high electron withdrawal ability, leading to systems being less stable, we suggest the borinium anions [(N 5 ) 2 BH] 2– and [(N 5 ) 2 B 2 H 2 ] 2– as possible synthetic targets of novel environmentally friendly catalysts.

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

Third-Order Møller–Plesset Theory Made More Useful? The Role of Density Functional Theory Orbitals

The practical utility of Møller–Plesset (MP) perturbation theory is severely constrained by the use of Hartree–Fock (HF) orbitals. It has recently been shown that the use of regularized orbital-optimized MP2 orbitals and scaling of MP3 energy could lead to a significant reduction in MP3 error. In this work, we examine whether density functional theory (DFT)-optimized orbitals can be similarly employed to improve the performance of MP theory at both the MP2 and MP3 levels. We find that the use of DFT orbitals leads to significantly improved performance for prediction of thermochemistry, barrier heights, noncovalent interactions, and dipole moments relative to the standard HF-based MP theory. Indeed, MP3 (with or without scaling) with DFT orbitals is found to surpass the accuracy of coupled-cluster singles and doubles (CCSD) for several data sets. We also found that the results are not particularly functional sensitive in most cases (although range-separated hybrid functionals with low delocalization error perform the best). As such, MP3 based on DFT orbitals thus appears to be an efficient, noniterative O(N 6 ) scaling wave-function approach for single-reference electronic structure computations. Scaled MP2 with DFT orbitals is also found to be quite accurate in many cases, although modern double hybrid functionals are likely to be considerably more accurate.

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

Comparative Density Functional Theory Study of Magnetic Exchange Couplings in Dinuclear Transition-Metal Complexes

Multicenter transition-metal complexes (MCTMs) with magnetically interacting ions have been proposed as components for information-processing devices and storage units. For any practical application of MCTMs as magnetic units, it is crucial to characterize their magnetic behavior, and in particular, the isotropic magnetic exchange coupling, J, between its magnetic centers. Due to the large size of typical MCTMs, density functional theory is the only practical electronic structure method for evaluating the J coupling. Here, we assess the accuracy of different density functional approximations for predicting the magnetic couplings of eight dinuclear transition-metal complexes, including five dimanganese, two dicopper, and one divanadium with known reliable experimental J couplings spanning from ferromagnetic to strong antiferromagnetic. The density functionals considered include global hybrid functionals which mix semilocal density functional approximations and exact exchange with a fixed admixing parameter, six local hybrid functionals where the admixing parameters are extended to be spatially dependent, the SCAN and r 2 SCAN meta-generalized gradient approximations (GGAs), and two widely used GGAs. We found that global hybrids tested in this work have a tendency to over-correct the error in magnetic coupling parameters from the Perdew–Burke–Ernzerhof (PBE) GGA as seen for manganese complexes. The performance of local hybrid density functionals shows no improvement in terms of bias and is scattered without a clear trend, suggesting that more efforts are needed for the extension from global to local hybrid density functionals for this particular property. The SCAN and r 2 SCAN meta-GGAs are found to perform as well as benchmark global hybrids on most tested complexes. We further analyze the charge density redistribution of meta-GGAs as well as global and local hybrid density functionals with respect to that of PBE, in connection to the self-interaction error or delocalization error.

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

Nuclear–Electronic Orbital QM/MM Approach: Geometry Optimizations and Molecular Dynamics

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods allow simulations of chemical reactions in atomistic solvent and heterogeneous environments such as proteins. Herein, the nuclear–electronic orbital (NEO) QM/MM approach is introduced to enable the quantization of specified nuclei, typically protons, in the QM region using a method such as NEO-density functional theory (NEO-DFT). This approach includes proton delocalization, polarization, anharmonicity, and zero-point energy in geometry optimizations and dynamics. Expressions for the energies and analytical gradients associated with the NEO-QM/MM method, as well as the previously developed polarizable continuum model (NEO-PCM), are provided. Geometry optimizations of small organic molecules hydrogen bonded to water in either dielectric continuum solvent or explicit atomistic solvent illustrate that aqueous solvation can strengthen hydrogen-bonding interactions for the systems studied, as indicated by shorter intermolecular distances at the hydrogen-bond interface. We then performed a real-time direct dynamics simulation of a phenol molecule in explicit water using the NEO-QM/MM method. Furthermore, these developments and initial examples provide the foundation for future studies of nuclear–electronic quantum dynamics in complex chemical and biological environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2021

Charge Transfer and Spin Dynamics in a Zinc Porphyrin Donor Covalently Linked to One or Two Naphthalenediimide Acceptors

We report quantum coherence effects on charge transfer and spin dynamics in a system having two degenerate electron acceptors are studied using a zinc 5,10,15-tri(n-pentyI)-20-phenylporphyrin (ZnP) electron donor covalently linked to either one or two naphthalene-1,8:4,5-bis(dicarboximide) (NDI) electron acceptors using an anthracene (An) spacer, ZnP-An-NDI (1) and ZnP-An-NDI 2 (2), respectively. Following photoexcitation of 1 and 2 in toluene at 295 K, femtosecond transient absorption spectroscopy shows that the electron transfer (ET) rate constant for 2 is about three times larger than that of 1, which can be accounted for by the statistical nature of incoherent ET as well as the electron couplings for the charge separation reactions. In contrast, the rate constant for charge recombination (CR) of 1 is about 25% faster than that of 2. Using femtosecond transient infrared spectroscopy and theoretical analysis, we find that the electron on NDI 2 ∙- in 2 localizes onto one of the two NDIs prior to CR, thus precluding electronically coherent CR from NDI 2 ∙- . Conversely, CR in both 1 and 2 is spin coherent as indicated by the observation of a resonance in the 3* ZnP yield following CR as a function of applied magnetic field, giving spin-spin exchange interaction energies of 2J = 210 and 236 mT, respectively, where the line width of the resonance for 2 is greater than 1. These data show that while CR is a spin-coherent process, incoherent hopping of the electron between the two NDIs in 2, consistent with the lack of delocalization noted above, results in greater spin decoherence in 2 relative to 1.

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

Electronic Structure of Metallophthalocyanines, MPc (M = Fe, Co, Ni, Cu, Zn, Mg) and Fluorinated MPc

In this work, we compute the electronic structure and optical excitation energies of metal-free and transition-metal phthalocyanines (H2Pc and MPc for M = Fe, Co, Ni, Cu, Zn, Mg) using density functional theory with optimally tuned range-separated hybrid functionals (OT-RSH). We show that the OT-RSH approach provides photoemission spectra in quantitative agreement with experiments as well as optical band gaps within 10% of their experimental values, capturing the interplay of localized d-states and delocalized pi-pi* states for these organometallic compounds. We examine the tunability of MPcs and H 2 Pc through fluorination, resulting in quasi-rigid shifts of the molecular orbital energies by up to 0.7 eV. Our comprehensive data set provides a new computational benchmark for gas-phase phthalocyanines, significantly improving upon other density-functional-theory-based approaches.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗