Engineering Papers⌕ Search

Engineering topics

Hammes-Schiffer, Sharon

Publications and source records attributed to Hammes-Schiffer, Sharon.

18 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Lagrangian formulation of nuclear–electronic orbital Ehrenfest dynamics with real-time TDDFT for extended periodic systems

Here, we present a Lagrangian-based implementation of Ehrenfest dynamics with nuclear–electronic orbital (NEO) theory and real-time time-dependent density functional theory for extended periodic systems. In addition to a quantum dynamical treatment of electrons and selected protons, this approach allows for the classical movement of all other nuclei to be taken into account in simulations of condensed matter systems. Furthermore, we introduce a Lagrangian formulation for the traveling proton basis approach and propose new schemes to enhance its application for extended periodic systems. Validation and proof-of-principle applications are performed on electronically excited proton transfer in the o-hydroxybenzaldehyde molecule with explicit solvating water molecules. These simulations demonstrate the importance of solvation dynamics and a quantum treatment of transferring protons. This work broadens the applicability of the NEO Ehrenfest dynamics approach for studying complex heterogeneous systems in the condensed phase.

Calculus of variations↗
DOE OSTI2024

The role of an intramolecular hydrogen bond in the redox properties of carboxylic acid naphthoquinones

A bioinspired naphthoquinone model of the quinones in photosynthetic reaction centers but bearing an intramolecular hydrogen-bonded carboxylic acid has been synthesized and characterized electrochemically, spectroscopically, and computationally to provide mechanistic insight into the role of proton-coupled electron transfer (PCET) of quinone reduction in photosynthesis. The reduction potential of this construct is 370 mV more positive than the unsubstituted naphthoquinone. In addition to the reversible cyclic voltammetry, infrared spectroelectrochemistry confirms that the naphthoquinone/naphthoquinone radical anion couple is fully reversible. Calculated redox potentials agree with the experimental trends arising from the intramolecular hydrogen bond. Molecular electrostatic potentials illustrate the reversible proton transfer driving forces, and analysis of the computed vibrational spectra supports the possibility of a combination of electron transfer and PCET processes. The significance of PCET, reversibility, and redox potential management relevant to the design of artificial photosynthetic assemblies involving PCET processes is discussed.

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

Proton-Coupled Electron Transfer Mechanisms for CO 2 Reduction to Methanol Catalyzed by Surface-Immobilized Cobalt Phthalocyanine

Immobilized cobalt phthalocyanine (CoPc) is a highly promising architecture for the six-proton, six-electron reduction of CO 2 to methanol. This electroreduction process relies on proton-coupled electron transfer (PCET) reactions that can occur by sequential or concerted mechanisms. Immobilization on a conductive support such as carbon nanotubes or graphitic flakes can fundamentally alter the PCET mechanisms. We use density functional theory (DFT) calculations of CoPc adsorbed on an explicit graphitic surface model to investigate intermediates in the electroreduction of CO 2 to methanol. Our calculations show that the alignment of the CoPc and graphitic electronic states influences the reductive chemistry. These calculations also distinguish between charging the graphitic surface and reducing the CoPc and adsorbed intermediates as electrons are added to the system. This analysis allows us to identify the chemical transformations that are likely to be concerted PCET, defined for these systems as the mechanism in which protonation of a CO 2 reduction intermediate is accompanied by electron abstraction from the graphitic surface to the adsorbate without thermodynamically stable intermediates. Furthermore, this work establishes a mechanistic pathway for methanol production that is consistent with experimental observations and provides fundamental insight into how immobilization of the CoPc impacts its CO 2 reduction chemistry.

Adsorption↗
DOE OSTI2024

Roadmap on methods and software for electronic structure based simulations in chemistry and materials

This Roadmap article provides a succinct, comprehensive overview of the state of electronic structure methods and software for molecular and materials simulations. Seventeen distinct sections collect insights by 51 leading scientists in the field. Each contribution addresses the status of a particular area, as well as current challenges and anticipated future advances, with a particular eye towards software related aspects and providing key references for further reading. Foundational sections cover density functional theory and its implementation in real-world simulation frameworks, Green's function based many-body perturbation theory, wave-function based and stochastic electronic structure approaches, relativistic effects and semiempirical electronic structure theory approaches. Subsequent sections cover nuclear quantum effects, real-time propagation of the electronic structure, challenges for computational spectroscopy simulations, and exploration of complex potential energy surfaces. The final sections summarize practical aspects, including computational workflows for complex simulation tasks, the impact of current and future high-performance computing architectures, software engineering practices, education and training to maintain and broaden the community, as well as the status of and needs for electronic structure based modeling from the vantage point of industry environments. Overall, the field of electronic structure software and method development continues to unlock immense opportunities for future scientific discovery, based on the growing ability of computations to reveal complex phenomena, processes and properties that are determined by the make-up of matter at the atomic scale, with high precision.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A Many-Body Perspective of Nuclear Quantum Effects in Aqueous Clusters

Nuclear quantum effects play an important role in the structure and thermodynamics of aqueous systems. Here, by performing a many-body expansion with nuclear-electronic orbital (NEO) theory, we show that proton quantization can give rise to significant energetic contributions for many-body interactions spanning several molecules in single-point energy calculations of water clusters. Although zero-point motion produces a large increase in energy at the one-body level, nuclear quantum effects serve to stabilize higher-order molecular interactions. These results are significant because they demonstrate that nuclear quantum effects play a nontrivial role in many-body interactions of aqueous systems. Our approach also provides a pathway for incorporating nuclear quantum effects into water potential energy surfaces. The NEO approach is advantageous for many-body expansion analyses because it includes nuclear quantum effects directly in the energies.

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

Direct Evidence for a Sequential Electron Transfer–Proton Transfer Mechanism in the PCET Reduction of a Metal Hydroxide Catalyst

Here, the proton-coupled electron transfer (PCET) mechanism for the reaction M ox -OH + e − + H + → M red -OH 2 was de-termined through kinetic resolution of the independent electron transfer (ET) and proton transfer (PT) steps. The reaction of interest was triggered by visible light excitation of [Ru II (tpy)(bpy′)H 2 O] 2+ , Ru II -OH 2 , where tpy is 2,2′:6′,2″-terpyridine and bpy′ is 4,4′-diaminopropylsilatrane-2,2′-bipyridine, anchored to In 2 O 3 :Sn (ITO) thin films in aqueous solutions. Interfacial kinetics for the PCET reduction reaction were quantified by nanosecond transient absorption spectroscopy as a function of solution pH and applied potential. Data acquired from pH = 5-10 revealed a stepwise electron-transfer proton-transfer (ET-PT) mechanism, while kinetic measurements made below the pK a (Ru III -OH/OH 2 ) = 1.3 were used to study the analogous interfacial reaction where electron transfer was the only mechanistic step. Analysis of this data with a recently reported multi-channel kinetic model was used to construct a PCET zone diagram and supported the assignment of an ET-PT mechanism from pH = 5-10. Ultimately, this study represents a unique example amongst M ox -OH/M red -OH 2 reactivity where the protona-tion and oxidation state of the intermediate was kinetically and spectrally resolved to firmly establish the PCET mechanism.

Charge transfer↗
DOE OSTI2023

Assessing Implicit and Explicit Polarizable Solvation Models for Nuclear–Electronic Orbital Systems: Quantum Proton Polarization and Solvation Energetics

Accurate simulations of many chemical processes require the inclusion of both nuclear quantum effects and a solvent environment. The nuclear–electronic orbital (NEO) approach, which treats electrons and select nuclei quantum mechanically on the same level, combined with a polarizable continuum model (PCM) for the solvent environment, addresses this challenge in a computationally practical manner. In this work, the NEO-PCM approach is extended beyond the IEF-PCM (integral equation formalism PCM) and C-PCM (conductor PCM) approaches to the SS(V)PE (surface and simulation of volume polarization for electrostatics) and ddCOSMO (domain decomposed conductor-like screening model) approaches. IEF-PCM, SS(V)PE, C-PCM, and ddCOSMO all exhibit similar solvation energies as well as comparable nuclear polarization within the NEO framework. The calculations show that the nuclear density does not leak out of the molecular cavity because it is much more localized than the electronic density. Finally, the polarization of quantized protons is analyzed in both continuum solvent and explicit solvent environments described by the polarizable MB-pol model, illustrating the impact of specific hydrogen-bonding interactions captured only by explicit solvation. Furthermore, these calculations highlight the relationship among solvation formalism, nuclear polarization, and energetics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Nuclear–Electronic Orbital Quantum Mechanical/Molecular Mechanical Real-Time Dynamics

Simulating the nuclear–electronic quantum dynamics of large-scale molecular systems in the condensed phase is key for studying biologically and chemically important processes such as proton transfer and proton-coupled electron transfer reactions. Herein, the real-time nuclear–electronic orbital time-dependent density functional theory (RT-NEO-TDDFT) approach is combined with a hybrid quantum mechanical/molecular mechanical (QM/MM) strategy to enable the accurate description of coupled nuclear–electronic quantum dynamics in the presence of heterogeneous environments such as solvent or proteins. The densities of the electrons and quantum protons are propagated in real time, while the other nuclei are propagated classically on the instantaneous electron–proton vibronic surface. This approach is applied to phenol bound to lysozyme, intramolecular proton transfer in malonaldehyde, and nonequilibrium excited-state intramolecular proton transfer in o-hydroxybenzaldehyde. Furthermore, these examples illustrate that the RT-NEO-TDDFT framework, coupled with an atomistic representation of the environment, allows the simulation of condensed-phase systems that exhibit significant nuclear quantum effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Multicomponent Cholesky Decomposition: Application to Nuclear–Electronic Orbital Theory

The Cholesky decomposition technique is commonly used to reduce the memory requirement for storing two-particle repulsion integrals in quantum chemistry calculations that use atomic orbital bases. However, when quantum methods use multicomponent bases, such as nuclear–electronic orbitals, additional challenges are introduced due to asymmetric two-particle integrals. This work proposes several multicomponent Cholesky decomposition methods for calculations using nuclear–electronic orbital density functional theory. To analyze the errors in different Cholesky decomposition components, benchmark calculations using water clusters are carried out. The largest benchmark calculation is a water cluster (H 2 O) 27 where all 54 protons are treated quantum mechanically. Furthermore, this study provides energetic and complexity analyses to demonstrate the accuracy and performance of the proposed multicomponent Cholesky decomposition method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

General Kinetic Model for pH Dependence of Proton-Coupled Electron Transfer: Application to an Electrochemical Water Oxidation System

The pH dependence of proton-coupled electron transfer (PCET) reactions, which are critical to many chemical and biological processes, is a powerful probe for elucidating their fundamental mechanisms. Herein, a general, multichannel kinetic model is introduced to describe the pH dependence of both homogeneous and electrochemical PCET reactions. According to this model, a weak pH dependence can arise from the competition among multiple sequential and concerted PCET channels involving different forms of the redox species, such as protonated and deprotonated forms, as well as different proton donors and acceptors. The contribution of each channel is influenced by the relative populations of the reactant species, which often depend strongly on pH, leading to complex pH dependence of PCET apparent rate constants. This model is used to explain the origins of the experimentally observed weak pH dependence of the electrochemical PCET apparent rate constant for a ruthenium-based water oxidation catalyst attached to a tin-doped In2O3 (ITO) surface. The weak pH dependence is found to arise from the intrinsic differences in the rate constants of participating channels and the dependence of their relative contributions on pH. This model predicts that the apparent maximum rate constant will become pH-independent at higher pH, which is confirmed by experimental measurements. Our analysis also suggests that the dominant channels are electron transfer at lower pH and sequential PCET via electron transfer followed by fast proton transfer at higher pH. Furthermore, this work highlights the importance of considering multiple competing channels simultaneously for PCET processes.

Catalysts↗
DOE OSTI2023

Synthesis and Surface Attachment of Molecular Re(I) Hydride Species with Silatrane Functionalized Bipyridyl Ligands

Three molecular Re hydrides of the form ( R bpy)Re(CO) 3 H with 2,2′-bipyridine (bpy) ligands containing silatrane functional groups for surface attachment on metal oxide surfaces were synthesized. IR spectroscopy and cyclic voltammetry (CV) demonstrated that the complexes containing the silatrane functional groups have electronic properties similar to those of a control compound, which did not contain functional groups for attachment. Additionally, in a similar fashion to the control compound, the silatrane containing Re hydrides are electrocatalysts for the reduction of CO 2 to CO in solution. The silatrane containing complexes were immobilized on a thin layer of TiO 2 on Si, and the resulting composites were characterized using X-ray photoelectron and IR spectroscopy as well as cyclic voltammetry in the dark and under illumination. Control experiments indicated that the hydride complexes are not stable on the surface and degrade to species which contain a bpy ligand, three CO ligands, and an unknown ligand in the sixth site. Similarly, when one of the silatrane containing Re hydride complexes was immobilized on Si nanoparticles with a thin layer of SiO 2 or silica nanoparticles, the hydride ligand was lost. Density functional theory calculations were used to corroborate the observed behavior of hydride species on a surface. Altogether, this work demonstrates the difficulties associated with attaching well-defined molecular hydride complexes to metal oxide surfaces.

Anions↗
DOE OSTI2023

Proton-regulated alcohol oxidation for high-capacity ketone-based flow battery anolyte

Redox flow batteries have a unique architecture that potentially enables cost-effective long-duration energy storage to address the intermittency introduced by increased renewable integration for the decarbonization of the electric power sector. Targeted molecular engineering has demonstrated electrochemical reversibility in natively redox-inactive ketone molecules in aqueous electrolytes. Yet, the kinetics of fluorenone-based flow batteries continue to be limited by slow alcohol oxidation. We show how strategically designed proton regulators can accelerate alcohol oxidation and thus enhance battery kinetics. Fluorenone-based flow batteries with the organic additive ß-cyclodextrin demonstrate enhanced rate capability, high capacity, and long cycling. This study opens a new avenue to improve the kinetics of aqueous organic flow batteries by modulating the reaction pathway with a homogeneous catalyst.

25 ENERGY STORAGE↗
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 OSTI2023

Solvent Induced Proton Polarization within the Nuclear–Electronic Orbital Framework

To explicitly account for nuclear quantum effects and solvent environments in simulations of chemical processes, the nuclear-electronic orbital approach is coupled with a polarizable continuum model (PCM). This NEO–PCM approach is used to explore the influence of solvation on nuclear polarization through applications to a water dimer and a set of protonated water tetramers. Nuclear polarization in these species is analyzed in terms of changes in proton density and oxygen–hydrogen bond length. Solvation is shown to enhance nuclear polarization with increasing dielectric constant. For the water dimer, the internal, hydrogen-bonded proton is shown to polarize more than the external, free proton. Moreover, proton quantization leads to greater solvent polarization through their mutual polarization. Furthermore, these calculations highlight the complex interplay among electronic, nuclear, and solvent polarization in chemical systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗