Insights into the Electronic Dynamics in Chemical Reactions

How to search for and reveal a hidden intermediate? The ELF topological description of non-synchronicity in double proton transfer reactions under oriented external electric field

CONTEXT

The nature of double intermolecular proton transfer was studied with the ELF topological approach in two model dimers (the formic acid homodimer and the 1,2,3-triazole-guanidine heterodimer) under an oriented external electric field. It has been shown that each of the two dimers can have either a one-step (one transition state structure) or two-step (two transition state structures) reaction path, depending on the intensity and orientation of the external electric field. The presence of a singularly broad shoulder (plateau in the case of homodimer and plateau-like for heterodimer) around the formal transition state structure results from the strong asynchronicity of the reaction. A careful ELF topological analysis of the nature of protons, hydride (localized) or roaming (delocalized) proton, along the reaction path allowed us to unambiguously classify the one-step mechanisms governing the double-proton transfer reactions into three distinct classes: (1) concerted-synchronous, when two events (roaming proton regions) completely overlap, (2) concerted-asynchronous, when two events (roaming proton regions) partially overlap, and (3) two-stage one-step non-concerted, when two roaming proton regions are separated by a "hidden intermediate region". All the structures belonging to this separatrix region are of the zwitterion form.

METHODS

Geometry optimization of the stationary points on the potential energy surface was performed using density functional theory-wB97XD functional-in combination with the 6-311+ +G(2d, 2p) basis set for all the atoms. All first-principles calculations were performed using the Gaussian 09 quantum chemical packages. We also used the electron localization function (ELF) to reveal the nature of the proton along the reaction path: a bound proton (hydride) becomes a roaming proton (carrying a tiny negative charge ≈ 0.3 e) exchanging with two adjacent atoms via two attractors (topological critical points with (3, - 3) signature). The ELF analyses were performed using the TopMod package.

Ab initio and semi-empirical Molecular Dynamics simulations of chemical reactions in isolated molecules and in clusters

Recent progress in “on-the-fly” trajectory simulations of molecular reactions, using different electronic structure methods is discussed, with analysis of the insights that such calculations can provide and of the strengths and limitations of the algorithms available.

Theory of divalent main group H2 activation: electronics and quasiclassical trajectories

Density functional theory (DFT), absolutely localized molecular orbital (ALMO) analysis, and quasiclassical trajectories (QCTs) were used to study the structure, barrier heights, thermodynamics, electronic properties, and dynamics of dihydrogen (H2) activation by singlet divalent main group compounds (ER2; E = C, Si, Ge). ALMO energy and charge decomposition calculations reveal that in the transition state CR2 acts as an ambiphile toward H2 because of equal forward-bonding and back-bonding orbital stabilization while SiR2 and GeR2 act as nucleophiles with dominant orbital energy stabilization arising from ER2 to H2 donation. Frontier molecular orbital (FMO) energy gaps do not provide a reasonable estimate of energy stabilization gained between the ER2 and H2 in the transition state or an accurate description of the nucleophilic versus electrophilic character because of electron repulsion and orbital overlap influences that are neglected. In CR2 transition states, forward-bonding and back-bonding are maximized in the nonleast motion geometry. In contrast, SiR2/GeR2 transition states have side-on geometries to avoid electron-electron repulsion. Electron repulsion, rather than orbital interactions, also determines the relative barrier heights of CR2 versus SiR2/GeR2 reactions. Examination of barrier heights and reaction energies shows a clear kinetic-thermodynamic relationship for ER2 activation of H2. A computational survey of R groups on ER2 divalent atom centers was performed to explore the possibility for H2 activation to occur with a low barrier and thermodynamically reversible. QCTs show that dihydrogen approach and reaction with CR2 may involve geometries significantly different than the static transition-state structure. In contrast, trajectories for dihydrogen addition to SiR2 involve geometries close to the side-on approach suggested by the static transition-state structure. QCTs also demonstrate that addition of H2 to CR2 and SiR2 is dynamically concerted with the average time gap of bond formation between E-H bonds of approximately 11 and 21 fs, respectively.

Valence-bond description of chemical reactions on Born-Oppenheimer molecular dynamics trajectories

The nature of chemical bonds on dynamic paths was investigated using the complete active space valence-bond (CASVB) method and the Born-Oppenheimer dynamics. To extract the chemical bond picture during reactions, a scheme to collect contributions from several VB (resonance) structures into a small numbers of indices was introduced. In this scheme, a tree diagram for the VB structures is constructed with the numbers of the ionic bonds treated as generation. A pair of VB structures is related to each other if one VB structure is transferred into the other by changing a covalent bond to an ionic bond. The former and latter VB structures are named parent and child structures, respectively. The weights of the bond pictures are computed as the sum of the CASVB occupation numbers running from the top generation to the bottom along the descent of the VB structures. Thus, a number of CASVB occupation numbers are collected into a small number of indices, and a clear bond picture may be obtained from the CASVB wave function. The scheme was applied to the hydrogen exchange reaction H(2)+F-->H+HF and the Diels-Alder reaction C(5)H(6)(cyclopentadiene)+CH(2)=CH(2)(ethylene)-->C(7)H(10)(norbornene). In both the reactions, the scheme gave a clear picture for the Born-Oppenheimer dynamics trajectories. The reconstruction of the bonds during reactions was well described by following the temporal changes in weight.

When is a molecule properly solvated by a continuum model or in a cluster ansatz? A first-principles simulation of alanine hydration

In order to test the validity of the cluster ansatz approach as well as of the continuum model approach and to learn about the solvation shell, we carried out first-principles molecular dynamics simulations of the alanine hydration. Our calculations contained one alanine molecule dissolved in 60 water molecules. Dipole moments of individual molecules were derived by means of maximally localized Wannier functions. We observed an average dipole moment of about 16.0 D for alanine and of about 3.3 D for water. In particular, the average water dipole moment in proximity of alanine's COO(-) group decayed continously with increasing distance, while, surprisingly, close to the CH3 and NH3+ group, the dipole moment first rose before its value dropped. In a cluster ansatz approach, we considered snapshots of alanine surrounded by different water molecule shells. The dipole moments from the cluster approaches utilizing both maximally localized Wannier functions as well as natural population analysis served to approximate the dipole moments of the total trajectory. Sufficient convergence of the cluster ansatz approach is found for either of the two solvent shells around the polar groups and one solvent shell around the apolar groups or two solvent shells around the polar groups surrounded by a dieletric continuum.

Multicoordinate Driven Method for Approximating Enzymatic Reaction Paths: Automatic Definition of the Reaction Coordinate Using a Subset of Chemical Coordinates

We present a generalization of the reaction coordinate driven method to find reaction paths and transition states for complicated chemical processes, especially enzymatic reactions. The method is based on the definition of a subset of chemical coordinates; it is simple, robust, and suitable to calculate one or more alternative pathways, intermediate minima, and transition-state geometries. Though the results are approximate and the computational cost is relatively high, the method works for large systems, where others often fail. It also works when a certain reaction path competes with others having a lower energy barrier. Accordingly, the procedure is appropriate to test hypothetical reaction mechanisms for complicated systems and provides good initial guesses for more accurate methods. We present tests on a number of simple reactions and on several complicated chemical transformations and compare the results with those obtained by other methods. Calculation of the reaction path for the enzymatic hydrolysis of the substrate by dUTPase for an active-site model with 85 atoms, including several loosely bound water molecules, indicates that the method is feasible for the study of enzyme mechanisms.

The Microscopic Switching Mechanism of a [2]Catenane

We apply numerical simulations at an all-atom level to investigate the switching mechanism of a [2]catenane, a prototype of a molecular machine. This system is able to switch reversibly between two different stable states, upon external stimuli, with a time scale ranging from microseconds up to milliseconds, well over the typical domain of molecular dynamics (MD) computer simulations. However, combining a strategy recently developed for investigating rare events with ordinary MD, we are able to unravel the microscopic mechanism of the conformational rearrangements involved in the switching process, including dynamical effects. Along the path that connects the product and reactant state, we find several intermediate states characterized by pi-pi stacking interactions and hydrogen bonds. Moreover, counterions interact strongly with the system in a correlated way, in agreement with recent static calculations performed on [2]rotaxanes.