RESUMO
The enzyme ribonucleotide reductase (RNR), which catalyzes the reduction of ribonucleotides to deoxynucleotides, is vital for DNA synthesis, replication, and repair in all living organisms. Its mechanism requires long-range radical translocation over â¼32 Å through two protein subunits and the intervening aqueous interface. Herein, a kinetic model is designed to describe reversible radical transfer in Escherichia coli RNR. This model is based on experimentally studied photoRNR systems that allow the photochemical injection of a radical at a specific tyrosine residue, Y356, using a photosensitizer. The radical then transfers across the interface to another tyrosine residue, Y731, and continues until it reaches a cysteine residue, C439, which is primed for catalysis. This kinetic model includes radical injection, an off-pathway sink, radical transfer between pairs of residues along the pathway, and the conformational flipping motion of Y731 at the interface. Most of the input rate constants for this kinetic model are obtained from previous experimental measurements and quantum mechanical/molecular mechanical free-energy simulations. Ranges for the rate constants corresponding to radical transfer across the interface are determined by fitting to the experimentally measured Y356 radical decay times in photoRNR systems. This kinetic model illuminates the time evolution of radical transport along the tyrosine and cysteine residues following radical injection. Further analysis identifies the individual rate constants that may be tuned to alter the timescale and probability of the injected radical reaching C439. The insights gained from this kinetic model are relevant to biochemical understanding and protein-engineering efforts with potential pharmacological implications.
Assuntos
Cisteína , Proteínas de Escherichia coli , Escherichia coli , Ribonucleotídeo Redutases , Cisteína/química , Escherichia coli/enzimologia , Proteínas de Escherichia coli/química , Modelos Químicos , Simulação de Dinâmica Molecular , Ribonucleotídeo Redutases/química , Termodinâmica , Tirosina/químicaRESUMO
Redox-active residues, such as tyrosine and tryptophan, play important roles in a wide range of biological processes. The α3Y de novo protein, which is composed of three α helices and a tyrosine residue Y32, provides a platform for investigating the redox properties of tyrosine in a well-defined protein environment. Herein, the proton-coupled electron transfer (PCET) reaction that occurs upon oxidation of tyrosine in this model protein by a ruthenium photosensitizer is studied by using a vibronically nonadiabatic PCET theory that includes hydrogen tunneling and excited vibronic states. The input quantities to the analytical nonadiabatic rate constant expression, such as the diabatic proton potential energy curves and associated proton vibrational wave functions, reorganization energy, and proton donor-acceptor distribution functions, are obtained from density functional theory calculations on model systems and molecular dynamics simulations of the solvated α3Y protein. Two possible proton acceptors, namely, water or a glutamate residue in the protein scaffold, are explored. The PCET rate constant is greater when glutamate is the proton acceptor, mainly due to the more favorable driving force and shorter equilibrium proton donor-acceptor distance, although contributions from excited vibronic states mitigate these effects. Nevertheless, water could be the dominant proton acceptor if its equilibrium constant associated with hydrogen bond formation is significantly greater than that for glutamate. Although these calculations do not definitively identify the proton acceptor for this PCET reaction, they elucidate the conditions under which each proton acceptor can be favored. These insights have implications for tyrosine-based PCET in a wide variety of biochemical processes.
Assuntos
Simulação de Dinâmica Molecular , Oxirredução , Prótons , Tirosina , Tirosina/química , Tirosina/metabolismo , Transporte de Elétrons , Teoria da Densidade FuncionalRESUMO
Immobilized cobalt phthalocyanine (CoPc) is a highly promising architecture for the six-proton, six-electron reduction of CO2 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 CO2 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 CO2 reduction intermediate is accompanied by electron abstraction from the graphitic surface to the adsorbate without thermodynamically stable intermediates. 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 CO2 reduction chemistry.
RESUMO
The proton-coupled electron transfer (PCET) mechanism for the reaction Mox-OH + e- + H+ â Mred-OH2 was determined through the kinetic resolution of the independent electron transfer (ET) and proton transfer (PT) steps. The reaction of interest was triggered by visible light excitation of [RuII(tpy)(bpy')H2O]2+, RuII-OH2, where tpy is 2,2':6',2â³-terpyridine and bpy' is 4,4'-diaminopropylsilatrane-2,2'-bipyridine, anchored to In2O3: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 at pH = 5-10 revealed a stepwise electron transfer-proton transfer (ET-PT) mechanism, while kinetic measurements made below pKa(RuIII-OH/OH2) = 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 multichannel kinetic model was used to construct a PCET zone diagram and supported the assignment of an ET-PT mechanism at pH = 5-10. Ultimately, this study represents a unique example among Mox-OH/Mred-OH2 reactivity where the protonation and oxidation states of the intermediate were kinetically and spectrally resolved to firmly establish the PCET mechanism.
RESUMO
Colloidal gold nanoparticles (AuNPs) have myriad scientific and technological applications, but their fundamental redox chemistry is underexplored. Reported here are titration studies of oxidation and reduction reactions of aqueous AuNP colloids, which show that the AuNPs bind substantial hydrogen (electrons + protons) under mild conditions. The 5 nm AuNPs are reduced to a similar extent with reductants from borohydrides to H2 and are reoxidized back essentially to their original state by oxidants, including O2. The reactions were monitored via surface plasmon resonance (SPR) optical absorption, which was shown to be much more sensitive to surface H than to changes in solution conditions. Reductions with H2 occurred without pH changes, demonstrating that hydrogenation forms surface H rather than releasing H+. Computational studies suggested that an SPR blueshift was expected for H atom addition, while just electron addition likely would have caused a redshift. Titrations consistently showed a maximum redox change of the 5 nm NPs, independent of the reagent, corresponding to 9% of the total gold or â¼30% hydrogen surface coverage (â¼370 H per AuNP). Larger AuNPs showed smaller maximum fractional surface coverages. We conclude that H binds to the edge, corner, and defect sites of the AuNPs, which explains the stoichiometric limitation and the size effect. The finding of substantial and stable hydrogen on the AuNP surface under mild reducing conditions has potential implications for various applications of AuNPs in reducing environments, from catalysis to biomedicine. This finding contrasts with the behavior of bulk gold and with the typical electron-focused perspective in this field.
RESUMO
Chemists have long been fascinated by chirality, water, and interfaces, making tremendous progress in each research area. However, the chemistry emerging from the interplay of chirality, water, and interfaces has been difficult to study due to technical challenges, creating a barrier to elucidating biological functions at interfaces. Most biopolymers (proteins, DNA, and RNA) fold into macroscopic chiral structures to perform biological functions. Their folding requires water, but water behaves differently at interfaces where the bulk water hydrogen-bonding network terminates. A question arises as to how water molecules rearrange to minimize free energy at interfaces while stabilizing the macroscopic folding of biopolymers to support biological function. This question is central to solving many research challenges, including the molecular origin of biological homochirality, folding and insertion of proteins into cell membranes, and the design of heterogeneous biocatalysts. Researchers can resolve these challenges if they have the theoretical tools to accurately predict molecular behaviors of water and biopolymers at various interfaces. However, developing such tools requires validation by the experimental data. These experimental data are scarce because few physical methods can simultaneously distinguish chiral folding of the biopolymers, separate signals of interfaces from the overwhelming background of bulk solvent, and differentiate water in hydration shells of the polymers from water elsewhere.We recently illustrated these very capacities of chirality-sensitive vibrational sum frequency generation spectroscopy (chiral SFG). While chiral SFG theory dictates that the method is surface-specific under the condition of electronic nonresonance, we show the method can distinguish chiral folding of proteins and DNA and probe water structures in the first hydration shell of proteins at interfaces. Using amide I signals, we observe protein folding into ß-sheets without background signals from α-helices and disordered structures at interfaces, thereby demonstrating the effect of 2D crowding on protein folding. Also, chiral SFG signals of C-H stretches are silent from single-stranded DNA, but prominent for canonical antiparallel duplexes as well as noncanonical parallel duplexes at interfaces, allowing for sensing DNA secondary structures and hybridization. In establishing chiral SFG for detecting protein hydration structures, we observe an H218O isotopic shift that reveals water contribution to the chiral SFG spectra. Additionally, the phase of the O-H stretching bands flips when the protein chirality is switched from L to D. These experimental results agree with our simulated chiral SFG spectra of water hydrating the ß-sheet protein at the vacuum-water interface. The simulations further reveal that over 90% of the total chiral SFG signal comes from water in the first hydration shell. We conclude that the chiral SFG signals originate from achiral water molecules that assemble around the protein into a chiral supramolecular structure with chirality transferred from the protein. As water O-H stretches can reveal hydrogen-bonding interactions, chiral SFG shows promise in probing the structures and dynamics of water-biopolymer interactions at interfaces. Altogether, our work has created an experimental and computational framework for chiral SFG to elucidate biological functions at interfaces, setting the stage for probing the intricate chemical interplay of chirality, water, and interfaces.
Assuntos
Proteínas , Água , Proteínas/química , Análise Espectral/métodos , Dobramento de Proteína , HidrogênioRESUMO
Proton-coupled electron transfer (PCET) plays an essential role in a wide range of electrocatalytic processes. A vast array of theoretical and computational methods have been developed to study electrochemical PCET. These methods can be used to calculate redox potentials and pKa values for molecular electrocatalysts, proton-coupled redox potentials and bond dissociation free energies for PCET at metal and semiconductor interfaces, and reorganization energies associated with electrochemical PCET. Periodic density functional theory can also be used to compute PCET activation energies and perform molecular dynamics simulations of electrochemical interfaces. Various approaches for maintaining a constant electrode potential in electronic structure calculations and modeling complex interactions in the electric double layer (EDL) have been developed. Theoretical formulations for both homogeneous and heterogeneous electrochemical PCET spanning the adiabatic, nonadiabatic, and solvent-controlled regimes have been developed and provide analytical expressions for the rate constants and current densities as functions of applied potential. The quantum mechanical treatment of the proton and inclusion of excited vibronic states have been shown to be critical for describing experimental data, such as Tafel slopes and potential-dependent kinetic isotope effects. The calculated rate constants can be used as input to microkinetic models and voltammogram simulations to elucidate complex electrocatalytic processes.
Assuntos
Elétrons , Prótons , Transporte de Elétrons , Cinética , Simulação de Dinâmica MolecularRESUMO
Hydrogen energy-based electrochemical energy conversion technologies offer the promise of enabling a transition of the global energy landscape from fossil fuels to renewable energy. Here, we present a comprehensive review of the fundamentals of electrocatalysis in alkaline media and applications in alkaline-based energy technologies, particularly alkaline fuel cells and water electrolyzers. Anion exchange (alkaline) membrane fuel cells (AEMFCs) enable the use of nonprecious electrocatalysts for the sluggish oxygen reduction reaction (ORR), relative to proton exchange membrane fuel cells (PEMFCs), which require Pt-based electrocatalysts. However, the hydrogen oxidation reaction (HOR) kinetics is significantly slower in alkaline media than in acidic media. Understanding these phenomena requires applying theoretical and experimental methods to unravel molecular-level thermodynamics and kinetics of hydrogen and oxygen electrocatalysis and, particularly, the proton-coupled electron transfer (PCET) process that takes place in a proton-deficient alkaline media. Extensive electrochemical and spectroscopic studies, on single-crystal Pt and metal oxides, have contributed to the development of activity descriptors, as well as the identification of the nature of active sites, and the rate-determining steps of the HOR and ORR. Among these, the structure and reactivity of interfacial water serve as key potential and pH-dependent kinetic factors that are helping elucidate the origins of the HOR and ORR activity differences in acids and bases. Additionally, deliberately modulating and controlling catalyst-support interactions have provided valuable insights for enhancing catalyst accessibility and durability during operation. The design and synthesis of highly conductive and durable alkaline membranes/ionomers have enabled AEMFCs to reach initial performance metrics equal to or higher than those of PEMFCs. We emphasize the importance of using membrane electrode assemblies (MEAs) to integrate the often separately pursued/optimized electrocatalyst/support and membranes/ionomer components. Operando/in situ methods, at multiscales, and ab initio simulations provide a mechanistic understanding of electron, ion, and mass transport at catalyst/ionomer/membrane interfaces and the necessary guidance to achieve fuel cell operation in air over thousands of hours. We hope that this Review will serve as a roadmap for advancing the scientific understanding of the fundamental factors governing electrochemical energy conversion in alkaline media with the ultimate goal of achieving ultralow Pt or precious-metal-free high-performance and durable alkaline fuel cells and related technologies.
Assuntos
Fontes de Energia Elétrica , Prótons , Hidrogênio/química , Oxigênio/química , ÁguaRESUMO
In the recently discovered proton-coupled energy transfer (PCEnT) mechanism, the transfer of electronic excitation energy between donor and acceptor chromophores is coupled to a proton transfer reaction. Herein, we develop a general theory for PCEnT and derive an analytical expression for the nonadiabatic PCEnT rate constant. This theory treats the transferring hydrogen nucleus quantum mechanically and describes the PCEnT process in terms of nonadiabatic transitions between reactant and product electron-proton vibronic states. The rate constant is expressed as a summation over these vibronic states, and the contribution of each pair of vibronic states depends on the square of the vibronic coupling as well as the spectral convolution integral, which can be viewed as a generalization of the Förster-type spectral overlap integral for vibronic rather than electronic states. The convolution integral also accounts for the common vibrational modes shared by the donor and acceptor chromophores for intramolecular PCEnT. We apply this theory to model systems to investigate the key features of PCEnT processes. The excited vibronic states can contribute significantly to the total PCEnT rate constant, and the common modes can either slow down or speed up the process. Because the pairs of vibronic states that contribute the most to the PCEnT rate constant may correspond to spectroscopically dark states, PCEnT could occur even when there is no apparent overlap between the donor emission and acceptor absorption spectra. This theory will assist in the interpretation of experimental data and will guide the design of additional PCEnT systems.
RESUMO
Chirality-selective vibrational sum frequency generation (chiral SFG) spectroscopy has emerged as a powerful technique for the study of biomolecular hydration water due to its sensitivity to the induced chirality of the first hydration shell. Thus far, water O-H vibrational bands in phase-resolved heterodyne chiral SFG spectra have been fit using one Lorentzian function per vibrational band, and the resulting fit has been used to infer the underlying frequency distribution. Here, we show that this approach may not correctly reveal the structure and dynamics of hydration water. Our analysis illustrates that the chiral SFG responses of symmetric and asymmetric O-H stretch modes of water have opposite phase and equal magnitude and are separated in energy by intramolecular vibrational coupling and a heterogeneous environment. The sum of the symmetric and asymmetric responses implies that an O-H stretch in a heterodyne chiral SFG spectrum should appear as two peaks with opposite phase and equal amplitude. Using pairs of Lorentzian functions to fit water O-H stretch vibrational bands, we improve spectral fitting of previously acquired experimental spectra of model ß-sheet proteins and reduce the number of free parameters. The fitting allows us to estimate the vibrational frequency distribution and thus reveals the molecular interactions of water in hydration shells of biomolecules directly from chiral SFG spectra.
RESUMO
Vibrational strong coupling (VSC) provides a novel means to modify chemical reactions and energy transfer pathways. To efficiently model chemical dynamics under VSC in the collective regime, herein a hybrid quantum mechanical/molecular mechanical (QM/MM) cavity molecular dynamics (CavMD) scheme is developed and applied to an experimentally studied chemical system. This approach can achieve linear scaling with respect to the number of molecules for a dilute solution under VSC by assuming that each QM solute molecule is surrounded by an independent MM solvent bath. Application of this approach to a dilute solution of Fe(CO)5 in n-dodecane under VSC demonstrates polariton dephasing to the dark modes and polariton-enhanced molecular nonlinear absorption. These simulations predict that strongly exciting the lower polariton may provide an energy transfer pathway that selectively excites the equatorial CO vibrations rather than the axial CO vibrations. Moreover, these simulations also directly probe the cavity effect on the dynamics of the Fe(CO)5 Berry pseudorotation reaction for comparison to recent two-dimensional infrared spectroscopy experiments. This theoretical approach is applicable to a wide range of other polaritonic systems and provides a tool for exploring the use of VSC for selective infrared photochemistry.
Assuntos
Simulação de Dinâmica Molecular , Vibração , Espectrofotometria Infravermelho/métodos , Fenômenos Químicos , Transferência de EnergiaRESUMO
Leveraging localized surface plasmon resonances of metal nanoparticles to trigger chemical reactions is a promising approach for heterogeneous catalysis. First-principles modeling of such processes is challenging due to the large number of electrons and electronic excited states as well as the significance of nuclear quantum effects when hydrogen is involved. Herein, the nonadiabatic nuclear-electronic quantum dynamics of plasmon-induced H2 photodissociation near an Al13- cluster is simulated with real-time nuclear-electronic orbital time-dependent density functional theory (RT-NEO-TDDFT). This approach propagates the nonequilibrium quantum dynamics of both electrons and protons. The plasmonic oscillations are shown to inject hot electrons into the antibonding orbital of H2, thereby inducing H2 dissociation. The quantum mechanical treatment of the hydrogen nuclei leads to faster H2 photodissociation and slightly larger isotope effects. Analysis of the nonequilibrium electronic density suggests that these findings stem from enhanced excited-state electronic coupling between the plasmonic mode and the H2 antibonding orbital due to proton delocalization or zero-point energy effects. Given the low computational overhead for including nuclear quantum effects with the RT-NEO-TDDFT approach, this work paves the way for simulating nonadiabatic nuclear-electronic quantum dynamics in other plasmonic systems.
RESUMO
Ribonucleotide reductase (RNR) regulates DNA synthesis and repair in all organisms. The mechanism of Escherichia coli RNR requires radical transfer over a proton-coupled electron transfer (PCET) pathway spanning â¼32 Å across two protein subunits. A key step along this pathway is the interfacial PCET reaction between Y356 in the ß subunit and Y731 in the α subunit. Herein, this PCET reaction between two tyrosines across an aqueous interface is explored with classical molecular dynamics and quantum mechanical/molecular mechanical (QM/MM) free energy simulations. The simulations suggest that the water-mediated mechanism involving double proton transfer through an intervening water molecule is thermodynamically and kinetically unfavorable. The direct PCET mechanism between Y356 and Y731 becomes feasible when Y731 is flipped toward the interface and is predicted to be approximately isoergic with a relatively low free energy barrier. This direct mechanism is facilitated by the hydrogen bonding of water to both Y356 and Y731. These simulations provide fundamental insights into radical transfer across aqueous interfaces.
Assuntos
Ribonucleotídeo Redutases , Tirosina , Tirosina/metabolismo , Prótons , Ribonucleotídeo Redutases/metabolismo , Elétrons , Modelos Moleculares , Transporte de Elétrons , Escherichia coli/metabolismo , ÁguaRESUMO
Interfacial electric fields play a critical role in electrocatalysis and are often characterized by using vibrational probes attached to an electrode surface. Understanding the physical principles dictating the impact of the applied electrode potential on the vibrational probe frequency is important. Herein, a comparative study is performed for two molecular probes attached to a gold electrode. Both probes contain a nitrile (CN) group, but 4-mercaptobenzonitrile (4-MBN) exhibits continuous conjugation from the electrode through the nitrile group, whereas this conjugation is interrupted for 2-(4-mercaptophenyl)acetonitrile (4-MPCN). Periodic density functional theory calculations predict that the CN vibrational frequency shift of the 4-MBN system is dominated by induction, which is a through-bond polarization effect, leading to a strong potential dependence that does not depend significantly on the orientation of the CN bond relative to the surface. In contrast, the CN vibrational frequency shift of the 4-MPCN system is influenced less by induction and more by through-space electric field effects, leading to a weaker potential dependence and a greater orientation dependence. These theoretical predictions were confirmed by surface-enhanced Raman spectroscopy experiments. Balancing through-bond and through-space electrostatic effects may assist in the fundamental understanding and design of electrocatalytic systems.
RESUMO
Understanding breaking and formation of Lewis bonds at an electrified interface is relevant to a large range of phenomena, including electrocatalysis and electroadsorption. The complexities of interfacial environments and associated reactions often impede a systematic understanding of this type of bond at interfaces. To address this challenge, we report the creation of a main group classic Lewis acid-base adduct on an electrode surface and its behavior under varying electrode potentials. The Lewis base is a self-assembled monolayer of mercaptopyridine and the Lewis acid is BF3, forming a Lewis bond between nitrogen and boron. The bond is stable at positive potentials but cleaves at potentials more negative of approximately -0.3 V vs Ag/AgCl without an associated current. We also show that if the Lewis acid BF3 is supplied from a reservoir of Li+BF4- electrolyte, the cleavage is completely reversible. We propose that the N-B Lewis bond is affected both by the field-induced intramolecular polarization (electroinduction) and by the ionic structures and ionic equilibria near the electrode. Our results indicate that the second effect is responsible for the Lewis bond cleavage at negative potentials. This work is relevant to understanding the fundamentals of electrocatalytic and electroadsorption processes.
RESUMO
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. This work highlights the importance of considering multiple competing channels simultaneously for PCET processes.
RESUMO
The controlled generation of nitric oxide (NO) from endogenous sources, such as S-nitrosoglutathione (GSNO), has significant implications for biomedical implants due to the vasodilatory and other beneficial properties of NO. The water-stable metal-organic framework (MOF) Cu-1,3,5-tris[1H-1,2,3-triazol-5-yl]benzene has been shown to catalyze the production of NO and glutathione disulfide (GSSG) from GSNO in aqueous solution as well as in blood. Previous experimental work provided kinetic data for the catalysis of the 2GSNO â 2NO + GSSG reaction, leading to various proposed mechanisms. Herein, this catalytic process is examined using density functional theory. Minimal functional models of the Cu-MOF cluster and glutathione moieties are established, and three distinct catalytic mechanisms are explored. The most thermodynamically favorable mechanism studied is consistent with prior experimental findings. This mechanism involves coordination of GSNO to copper via sulfur rather than nitrogen and requires a reductive elimination that produces a Cu(I) intermediate, implicating a redox-active copper site. The experimentally observed inhibition of reactivity at high pH values is explained in terms of deprotonation of a triazole linker, which decreases the structural stability of the Cu(I) intermediate. These fundamental mechanistic insights may be generally applicable to other MOF catalysts for NO generation.
Assuntos
Estruturas Metalorgânicas , Óxido Nítrico , Óxido Nítrico/química , S-Nitrosoglutationa , Cobre/farmacologia , Dissulfeto de Glutationa , Glutationa/química , CatáliseRESUMO
The coupled quantum dynamics of electrons and protons is ubiquitous in many dynamical processes involving light-matter interaction, such as solar energy conversion in chemical systems and photosynthesis. A first-principles description of such nuclear-electronic quantum dynamics requires not only the time-dependent treatment of nonequilibrium electron dynamics but also that of quantum protons. Quantum mechanical correlation between electrons and protons adds further complexity to such coupled dynamics. Here we extend real-time nuclear-electronic orbital time-dependent density functional theory (RT-NEO-TDDFT) to periodic systems and perform first-principles simulations of coupled quantum dynamics of electrons and protons in complex heterogeneous systems. The process studied is an electronically excited-state intramolecular proton transfer of o-hydroxybenzaldehyde in water and at a silicon (111) semiconductor-molecule interface. These simulations illustrate how environments such as hydrogen-bonding water molecules and an extended material surface impact the dynamical process on the atomistic level. Depending on how the molecule is chemisorbed on the surface, excited-state electron transfer from the molecule to the semiconductor surface can inhibit ultrafast proton transfer within the molecule. This Letter elucidates how heterogeneous environments influence the balance between the quantum mechanical proton transfer and excited electron dynamics. The periodic RT-NEO-TDDFT approach is applicable to a wide range of other photoinduced heterogeneous processes.
RESUMO
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. These calculations highlight the relationship among solvation formalism, nuclear polarization, and energetics.
RESUMO
Cryogenic ion vibrational spectra of D2-tagged cyanobenzoate (CBA) derivatives are obtained and analyzed to characterize the intrinsic spectroscopic responses of the -CO2- headgroup to its location on the ring in both the isolated anions and the cationic complexes with divalent metal ions, M2+ (M = Mg, Ca, Sr). The benzonitrile functionality establishes the different ring isomers (para, meta, ortho) according to the location of the carboxylate and provides an additional reporter on the molecular response to the proximal charge center. The aromatic carboxylates display shifts slightly smaller than those observed for a related aliphatic system upon metal ion complexation. Although the CBA anions display very similar band patterns for all three ring positions, upon complexation with metal ions, the ortho isomer yields dramatically different spectral responses in both the -CO2- moiety and the CN group. This behavior is traced to the emergence of a tridentate binding motif unique to the ortho isomer in which the metal ions bind to both the oxygen atoms of the carboxylate group and the N atom of the cyano group. In that configuration, the -CO2- moiety is oriented perpendicular to the phenyl ring, and the CN stretching fundamental is both strong and red-shifted relative to its behavior in the isolated neutral. The behaviors of the metal-bound ortho complexes occur in contrast to the usual blue shifts associated with "Lewis" type binding of metal ions end-on to -CN. The origins of these spectroscopic features are analyzed with the aid of electronic structure calculations, which also explore differences expected for complexation of monovalent cations to the ortho carboxylate. The resulting insights have implications for understanding the balance between electrostatic and steric interactions at metal binding sites in chemical and biological systems.