Ensemble‐Averaged QM/MM Kinetic Isotope Effects for the SN2 Reaction of Cyanide Anions with Chloroethane in DMSO Solution

The existence of solvent fluctuations leads to populations of reactant‐state (RS) and transition‐state (TS) configurations and implies that property calculations must include appropriate averaging over distributions of values for individual configurations. Average kinetic isotope effects 〈KIE〉 for NC^?+EtCl→NCEt+Cl^? in DMSO solution at 30 °C are best obtained as the ratio 〈f_RS〉/〈f_TS〉 of isotopic partition function ratios separately averaged over all RS and TS configurations. In this way the hybrid AM1/OPLS‐AA potential yields 〈KIE〉 values for all six isotopic substitutions (2° α‐²H₂, 2° β‐²H₃, α‐¹¹C/¹⁴C, leaving group ³⁷Cl, and nucleophile ¹³C and ¹⁵N) for this reaction in the correct direction as measured experimentally. These thermally‐averaged calculated KIEs may be compared meaningfully with experiment, and only one of them differs in magnitude from the experimental value by more than one standard deviation from the mean. This success contrasts with previous KIE calculations based upon traditional methods without averaging. The isotopic partition function ratios are best evaluated using all (internal) vibrational and (external) librational frequencies obtained from Hessians determined for subsets of atoms, relaxed to local minima or saddle points, within frozen solvent environments of structures sampled along molecular dynamics trajectories for RS and TS. The current method may perfectly well be implemented with other QM or QM/MM methods, and thus provides a useful tool for investigating KIEs in relation to studies of chemical reaction mechanisms in solution or catalyzed by enzymes.     

Geometric H/D isotope effects and cooperativity of the hydrogen bonds in porphycene.

We investigate the primary, secondary, and vicinal hydrogen/deuterium (H/D) isotope effects on the geometry of the two intramolecular hydrogen bonds in porphycene. Multidimensional potential energy surfaces describing the anharmonic motion in the vicinity of the trans isomer are calculated for the different symmetric (HH/DD) and asymmetric (HD) isotopomers. From the solution of the nuclear Schr?dinger equation the ground-state wavefunction is obtained, which is further used to determine the quantum corrections to the classical equilibrium geometries of the hydrogen bonds and thus the geometric isotope effects. In particular, it is found that the hydrogen bonds are cooperative, that is, both expand simultaneously even in the case of an asymmetric isotopic substitution. The theoretical predictions compare favorably with NMR chemical-shift data.     

Libra: An open‐Source “methodology discovery” library for quantum and classical dynamics simulations

The “methodology discovery” library for quantum and classical dynamics simulations is presented. One of the major foci of the code is on nonadiabatic molecular dynamics simulations with model and atomistic Hamiltonians treated on the same footing. The essential aspects of the methodology, design philosophy, and implementation are discussed. The code capabilities are demonstrated on a number of model and atomistic test cases. It is demonstrated how the library can be used to study methodologies for quantum and classical dynamics, as well as a tool for performing detailed atomistic studies of nonadiabatic processes in molecular systems. The source code and additional information are available on the Web at http://www.acsu.buffalo.edu/~alexeyak/libra/index.html . © 2016 Wiley Periodicals, Inc.     

How Frustrated Lewis Acid/Base Systems Pass through Transition‐State Regions: H2 Cleavage by [tBu3P/B(C6F5)3]

We investigate the transition‐state (TS) region of the potential energy surface (PES) of the reaction tBu₃P+H₂+B(C₆F₅)₃→tBu₃P‐H⁽⁺⁾+^(?)H?B(C₆F₅)₃ and the dynamics of the TS passage at room temperature. Owing to the conformational inertia of the phosphane???borane pocket involving heavy tBu₃P and B(C₆F₅)₃ species and features of the PES E(P???H, B???H | B???P) as a function of P???H, B???H, and B???P distances, a typical reactive scenario for this reaction is a trajectory that is trapped in the TS region for a period of time (about 350 fs on average across all calculated trajectories) in a quasi‐bound state (scattering resonance). The relationship between the timescale of the TS passage and the effective conformational inertia of the phosphane???borane pocket leads to a prediction that isotopically heavier Lewis base/Lewis acid pairs and normal counterparts could give measurably different reaction rates. Herein, the predicted quasi‐bound state could be verified in molecular collision experiments involving femtosecond spectroscopy.     

Proton transfer 200 years after von Grotthuss: insights from ab initio simulations.

In the last decade, ab initio simulations and especially Car-Parrinello molecular dynamics have significantly contributed to the improvement of our understanding of both the physical and chemical properties of water, ice, and hydrogen-bonded systems in general. At the heart of this family of in silico techniques lies the crucial idea of computing the many-body interactions by solving the electronic structure problem "on the fly" as the simulation proceeds, which circumvents the need for pre-parameterized potential models. In particular, the field of proton transfer in hydrogen-bonded networks greatly benefits from these technical advances. Here, several systems of seemingly quite different nature and of increasing complexity, such as Grotthuss diffusion in water, excited-state proton-transfer in solution, phase transitions in ice, and protonated water networks in the membrane protein bacteriorhodopsin, are discussed in the realms of a unifying viewpoint.     

Nuclear quantum effects induce metallization of dense solid molecular hydrogen

We present an accurate computational study of the electronic structure and lattice dynamics of solid molecular hydrogen at high pressure. The band‐gap energies of the , Pc, and structures at pressures of 250, 300, and 350 GPa are calculated using the diffusion quantum Monte Carlo (DMC) method. The atomic configurations are obtained from ab initio path‐integral molecular dynamics (PIMD) simulations at 300 K and 300 GPa to investigate the impact of zero‐point energy and temperature‐induced motion of the protons including anharmonic effects. We find that finite temperature and nuclear quantum effects reduce the band‐gaps substantially, leading to metallization of the and Pc phases via band overlap; the effect on the band‐gap of the structure is less pronounced. Our combined DMC‐PIMD simulations predict that there are no excitonic or quasiparticle energy gaps for the and Pc phases at 300 GPa and 300 K. Our results also indicate a strong correlation between the band‐gap energy and vibron modes. This strong coupling induces a band‐gap reduction of more than 2.46 eV in high‐pressure solid molecular hydrogen. Comparing our DMC‐PIMD with experimental results available, we conclude that none of the structures proposed is a good candidate for phases III and IV of solid hydrogen. © 2017 Wiley Periodicals, Inc.     

The Role of Hydrogen‐Bonding Interactions in the Ultrafast Relaxation Dynamics of the Excited States of 3‐ and 4‐Aminofluoren‐9‐ones

The dynamics of the excited states of 3‐ and 4‐aminofluoren‐9‐ones (3AF and 4AF, respectively) are investigated in different kinds of solvents by using a subpicosecond time‐resolved absorption spectroscopic technique. They undergo hydrogen‐bonding interaction with protic solvents in both the ground and excited states. However, this interaction is more significant in the lowest excited singlet (S₁) state because of its substantial intramolecular charge‐transfer character. Significant differences in the spectroscopic characteristics and temporal dynamics of the S₁ states of 3AF and 4AF in aprotic and protic solvents reveal that the intermolecular hydrogen‐bonding interaction between the S₁ state and protic solvents plays an important role in its relaxation process. Perfect linear correlation between the relaxation times of the S₁ state and the longitudinal relaxation times (τ_L) of alcoholic solvents confirms the prediction regarding the solvation process via hydrogen‐bond reorganization. In the case of weakly interacting systems, the relaxation process can be well described by a dipolar solvation‐like process involving rotation of the OH groups of the alcoholic solvents, whereas in solvents having a strong hydrogen‐bond‐donating ability, for example, methanol and trifluoroethanol, it involves the conversion of the non‐hydrogen‐bonded form to the hydrogen‐bonded complex of the S₁ state. Efficient radiationless deactivation of the S₁ state of the aminofluorenones by protic solvents is successfully explained by the energy‐gap law, by using the energy of the fully solvated S₁ state determined from the time‐resolved spectroscopic data.     

Proton and Deuteron Transfer Reactions in Molecular Nanoclusters

Proton and deuteron transfer rates and mechanisms are studied in polar molecular nanoclusters. The cluster environment strongly influences the reaction rate and the nature of these changes is studied as a function of the cluster size. The stabilities of the covalent reactant and polar product states change with cluster size and this effect alters both the equilibrium properties and transfer rate. The proton and deuteron are light quantum particles and the quantum character of the rate process is reflected in the magnitude of the kinetic isotope effect. Our mixed quantum–classical rate simulations indicate that the magnitude of the isotope effect decreases as the cluster size increases. More generally, our study shows how quantum effects combined with structural nanosolvation effects can lead to changes in reaction rates and mechanisms which should be applicable to many quantum charge transfer reactions in molecular nanoclusters.     

Ultrafast Relaxation Dynamics of the Excited States of 1‐Amino‐ and 1‐(N,N‐Dimethylamino)‐fluoren‐9‐ones

The dynamics of the excited states of 1‐aminofluoren‐9‐one (1AF) and 1‐(N,N‐dimethylamino)‐fluoren‐9‐one (1DMAF) are investigated by using steady‐state absorption and fluorescence as well as subpicosecond time‐resolved absorption spectroscopic techniques. Following photoexcitation of 1AF, which exists in the intramolecular hydrogen‐bonded form in aprotic solvents, the excited‐state intramolecular proton‐transfer reaction is the only relaxation process observed in the excited singlet (S₁) state. However, in protic solvents, the intramolecular hydrogen bond is disrupted in the excited state and an intermolecular hydrogen bond is formed with the solvent leading to reorganization of the hydrogen‐bond network structure of the solvent. The latter takes place in the timescale of the process of solvation dynamics. In the case of 1DMAF, the main relaxation pathway for the locally excited singlet, S₁(LE), or S₁(ICT) state is the configurational relaxation, via nearly barrierless twisting of the dimethylamino group to form the twisted intramolecular charge‐transfer, S₁(TICT), state. A crossing between the excited‐state and ground‐state potential energy curves is responsible for the fast, radiationless deactivation and nonemissive character of the S₁(TICT) state in polar solvents, both aprotic and protic. However, in viscous but strong hydrogen‐bond‐donating solvents, such as ethylene glycol and glycerol, crossing between the potential energy surfaces for the ground electronic state and the hydrogen‐bonded complex formed between the S₁(TICT) state and the solvent is possibly avoided and the hydrogen‐bonded complex is weakly emissive.     

Nuclear quantum effects at aqueous metal interfaces captured by molecular dynamics simulations

In this review, we summarize the recent development in modeling nuclear quantum effects at aqueous metal interfaces. First, we review the nuclear quantum effects on the water-metal interface at ultrahigh vacuum. Then, we illustrate the nuclear quantum effects at the potential of zero charge conditions. At last, we give some outlook for the perspective work in modeling the nuclear quantum effects at electrochemical interfaces and some practical simulation strategies.     

Ab initio static and molecular dynamics study of 4-styrylpyridine.

We report an in‐depth theoretical study of 4‐styrylpyridine in its singlet S₀ ground state. The geometries and the relative stabilities of the trans and cis isomers were investigated within density functional theory (DFT) as well as within Hartree–Fock (HF), second‐order Møller–Plesset (MP2), and coupled cluster (CC) theories. The DFT calculations were performed using the B3LYP and PBE functionals, with basis sets of different qualities, and gave results that are very consistent with each other. The molecular structure is thus predicted to be planar at the energy minimum, which is associated with the trans conformation, and to become markedly twisted at the minimum of higher energy, which is associated with the cis conformation. The results of the calculations performed with the post‐HF methods approach those obtained with the DFT methods, provided that the level of treatment of the electronic correlation is high enough and that sufficiently flexible basis sets are used. Calculations carried out within DFT also allowed the determination of the geometry and the energy of the molecule at the biradicaloid transition state associated with the thermal cis?trans isomerization and at the transition states associated with the enantiomerization of the cis isomer and with the rotations of the pyridinyl and phenyl groups in the trans and cis isomers. Car–Parrinello molecular dynamics simulations were also performed at 50, 150, and 300 K using the PBE functional. The studies allowed us to evidence the highly flexible nature of the molecule in both conformations. In particular, the trans isomer was found to exist mainly in a nonplanar form at finite temperatures, while the rotation of the pyridinyl ring in the cis isomer was incidentally observed to take place within ≈1 ps during the simulation carried out at 150 K on this isomer.     

Finite element interpolation for combined classical/quantum mechanical molecular dynamics simulations

A method is presented to interpolate the potential energy function for a part of a system consisting of a few degrees of freedom, such as a molecule in solution. The method is based on a modified finite element (FE) interpolation scheme. The aim is to save computer time when expensive methods such as quantum-chemical calculations are used to determine the potential energy function. The expensive calculations are only carried out if the molecule explores new unknown regions of the conformation space. If the molecule resides in regions previously explored, a cheap interpolation is performed instead of an expensive calculation, using known neighboring points. We report the interpolation techniques for the energies and the forces of the molecule, the handling of the FE mesh, and an application to a simple test example in molecular dynamics (MD) simulations. Good performance of the method was obtained (especially for MD simulations with a preceding Monte Carlo mesh generation) without losing accuracy. © 1997 John Wiley & Sons, Inc. J Comput Chem 18 : 1484–1495, 1997