Nuclear quantum and H/D isotope effects on three-centered bonding diborane: Path integral molecular dynamics simulations

Nuclear quantum and H/D isotope effects of bridging and terminal hydrogen atoms of diborane (B₂H₆) molecules were systematically studied by classical ab initio molecular dynamics (CLMD) and ab initio path integral molecular dynamics (PIMD) simulations with BHandHLYP/6-31++G** level of theory at room temperature (298.15 K). Calculated results clearly show that H/D isotope effect appears in the distribution of hydrogen (deuterium) of B₂H₆ (B₂D₆). Geometry of B₂H₆ also plays a significant role in the nuclear quantum effect proved by PIMD simulations, but slightly deviated from its equilibrium structure when simulated via CLMD simulation. The bond lengths between boron atoms R (B1 … B2) and the bridging hydrogen atoms R_HH (H_B1 … H_B2) of the B₂H₆ molecule obtained from PIMD simulations are slightly longer than those of the deuterated form of the diborane (B₂D₆) molecule. The principal component analysis (PCA) was also employed to distinguish the important modes of bridging hydrogen as related to the nuclear quantum and H/D isotope effects. The highest level of contribution obtained from PCA of PIMD simulations is bending, while various mixed vibrations with less contribution were also found. Therefore, the nuclear quantum and H/D isotope effects need to be taken into account for a better understanding of diborane geometry.     

Multi-component molecular orbital study of isotope effects on lithium hydride molecules with the configuration interaction scheme

The multi-component molecular orbital method, which can take account of the quantum effect of the electrons and nuclei, is applied to the calculation of lithium hydride isotope species with the configuration interaction (CI) scheme. The optimum basis set functions for quantum nuclei are proposed by the fully variational procedure under single electronic–single nuclear excitation CI level. The average internuclear distances and dipole moments for isotopic lithium hydride molecules calculated with small basis functions are reasonable agreement with the corresponding experimental values.     

Nonadiabatic relaxation dynamics of water anion cluster and its isotope effects by ring‐polymer molecular dynamics simulation

We have applied a recently developed hybrid quantum ring‐polymer molecular dynamics method to the nonadiabatic p → s relaxation dynamics in water anion clusters to understand the isotope effects observed in previous experiments. The average relaxation times for (H₂O)₅₀^? and (D₂O)₅₀^? were calculated at 120 and 207 fs, respectively, and are comparable to the experimental results. Therefore, we conclude that nuclear quantum effects play an essential role in understanding the observed isotope effects for water anion cluster nonadiabatic dynamics. The nonadiabatic relaxation mechanisms are also discussed in detail. © 2014 Wiley Periodicals, Inc.     

Solvent isotope effects on the solubility of argon and pair interactions in aqueous solution of hexamethylenetetramine

The solubilities of gaseous argon in solutions of hexamethylenetetramine (hmta) in H₂O and D₂O were measured at different concentrations, at five temperatures in the range of 283–318 K, and at a partial gas pressure of 101325 Pa by microvolumetry with an accuracy of ≤0.3%. The standard Setchenov coefficients for argon solutions and the thermodynamic parameters of Ar-hmta pairwise interactions were calculated from the data on the solubility. Interactions between hydrated Ar and hmta molecules are characterized by mutual repulsions. The structural states of water in the hydration sphere of hmta and in the bulk of solution are substantially different. Due to the difference in the mechanisms of hydration of hmta and Ar, the addition of the first component leads to a decrease in the solubility of the second component. Published inIzvestiya Akademii Nauk. Seriya Khimicheskaya, No 8, pp. 1389–1394, August, 2000.     

Subchronic toxicological effects of aquatic nitrobenzene on Medaka and Chinese rare minnow

Based on the simulative exposure experiments of nitrobenzene to two small experimental fishes, Medaka (Oryzias latipes) and Chinese Rare Minnow (Gobiocypris rarus), it was found that nitrobenzene could decrease in a linear way in the static aquatic system. Both fishes could accumulate dose-related levels of nitrobenzene quickly and eliminate the compound rapidly when they were transferred to clean water. The Chinese rare minnow showed more sensitivity to the acute toxicity exposure of nitrobenzene than Medaka. Typical molecular biomarkers in oxidative defense system including superoxide dismu- tase (SOD), catalase (CAT) showed the chemical induced alterations. The changes of acetylcholi- nesterase (AChE) activities indicated that nitrobenzene might affect the normal neural function. Nitro- benzene exposure could also lead to obvious damaging effects on the target organs, such as gill and liver by hitopathological studies.     

The any particle molecular orbital approach: A short review of the theory and applications

The any particle molecular orbital (APMO) approach extends regular electronic structure methods to study atomic and molecular systems in which electrons and other particles are treated simultaneously as quantum waves. A number of electronic structure methodologies have been extended under the APMO framework and applied to investigate nuclear quantum effects including isotope effects and nuclear delocalization and to calculate proton binding energies and affinities. In addition, APMO methodologies have been employed to analyze physical and chemical properties of atomic and molecular systems containing exotic subatomic particles.     

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.     

Ab initio path‐integral calculations of kinetic and equilibrium isotope effects on base‐catalyzed RNA transphosphorylation models

Detailed understandings of the reaction mechanisms of RNA catalysis in various environments can have profound importance for many applications, ranging from the design of new biotechnologies to the unraveling of the evolutionary origin of life. An integral step in the nucleolytic RNA catalysis is self‐cleavage of RNA strands by 2′‐O‐transphosphorylation. Key to elucidating a reaction mechanism is determining the molecular structure and bonding characteristics of transition state. A direct and powerful probe of transition state is measuring isotope effects on biochemical reactions, particularly if we can reproduce isotope effect values from quantum calculations. This article significantly extends the scope of our previous joint experimental and theoretical work in examining isotope effects on enzymatic and nonenzymatic 2′‐O‐transphosphorylation reaction models that mimic reactions catalyzed by RNA enzymes (ribozymes), and protein enzymes such as ribonuclease A (RNase A). Native reactions are studied, as well as reactions with thio substitutions representing chemical modifications often used in experiments to probe mechanism. Here, we report and compare results from eight levels of electronic‐structure calculations for constructing the potential energy surfaces in kinetic and equilibrium isotope effects (KIE and EIE) computations, including a “gold‐standard” coupled‐cluster level of theory [CCSD(T)]. In addition to the widely used Bigeleisen equation for estimating KIE and EIE values, internuclear anharmonicity and quantum tunneling effects were also computed using our recently developed ab initio path‐integral method, that is, automated integration‐free path‐integral method. The results of this work establish an important set of benchmarks that serve to guide calculations of KIE and EIE for RNA catalysis. © 2014 Wiley Periodicals, Inc.     

Hydrogen bonding versus van der Waals interactions: competitive influence of noncovalent interactions on 2D self-assembly at the liquid-solid interface

The structures of the self-assembled monolayers of various 4-alkoxybenzoic acids physisorbed at the liquid-solid interface were established by employing scanning tunnelling microscopy (STM). This study has been essentially undertaken to explore the competitive influence of van der Waals and hydrogen-bonding interactions on the process of two-dimensional self-assembly. These acid derivatives form hydrogen-bonded dimers as expected; however, the dimers organise themselves in the form of relatively complex lamellae. The characteristic feature of these lamellae is the presence of regular discommensurations or kinks along the lamella propagation direction. The formation of kinked lamellae is discussed in light of the registry mechanism of the alkyl chains with the underlying graphite substrate. The location of the kinks along a lamella depends on the number (odd or even) of carbon atoms in the alkyl chain. This result indicates that concerted van der Waals interactions of the alkyl chain units introduce the odd/even chain-length effect on the surface-assembled supramolecular patterns. The odd/even effects are retained even upon complexation with a hydrogen-bond acceptor. However, as the solvent is changed from 1-phenyloctane to 1-octanoic acid, the kinked lamellae as well as the odd/even effects disappear. This solvent-induced convergence of supramolecular patterns is attained by means of co-crystallisation of octanoic acid molecules in the 2D crystal lattice, which is evident from high-resolution STM images. The solvent co-adsorption phenomenon is discussed in terms of competing van der Waals and hydrogen-bonding interactions.     

Nuclear Dynamics and Electronic Effects of Hydrogen on Solid Surfaces

Hydrogen is involved in a variety of chemical processes on surfaces. While hydrogen exhibits vibrational and rotational dynamics in its adsorption state, it in some cases undergoes diffusion into the substrate as well as on the surface, and participates in chemical reactions. Furthermore, hydrogen exchanges an electron with surfaces having a significant effect on the surface electronic structure. In this personal account, we review our recent studies on surface nuclear dynamics of hydrogen, hydrogen transport across surfaces, catalytic hydrogenation/isotope exchange reactions, and charge transfer between the surface and hydrogen by using a depth‐resolved technique of nuclear reaction analysis and a quantum‐state‐selective detection of resonance enhanced multiphoton ionization in combination with surface science techniques. As a future prospect, we refer to ultraslow μ spin rotation spectroscopy for a direct probe of the hydrogen charge state at surfaces.     

The scaled-particle theory and solvent isotope effects on the thermodynamic properties of inert solutes in water and methanol

The scaled-particle theory has been applied to the calculation of the thermodynamic changes associated with the formation of a cavity in several isotopic varieties of liquid water and methanol. From these results, the thermodynamic functions for the transfer of a cavity (or a hard-sphere solute) have been computed for the following solvent pairs: H₂OD₂O, H₂OH₂ ¹⁸O, H₂ ¹⁸OD₂ ¹⁸O, D₂OD₂ ¹⁸O, CH₃OHCH₃OD. For the last two of these solvents, density measurements required for the calculations were carried out as a function of temperature. The calculated deuterium solvent isotope effect on the heats and entropies of hard-sphere solutes in water is much greater than the¹⁸O isotope effect; the former also exhibits a more pronounced temperature dependence. The transfer functions computed for hard-sphere solutes are compared to experimental data on the transfer of various solutes from H₂O to D₂O and from CH₃OH to CH₃OD. In most of the cases examined, the cavity effect accounts for a large part of the transfer quantities measured for rare gases, hydrocarbons, and solutes containing a significant hydrocarbon substituent.     

Hydrogen bonding and dissociation effects on the gas phase proton transfer reactions of ozone

Recently, the proton affinity (PA) of ozone was experimentally determined by Cacace and Speranza [Science (1994) 265: 208] using a bracketing technique that involved the proton transfer (PT) reactions: O₃H⁺+B⇒O₃+BH⁺; for different Br?nsted bases B. These authors showed that the simple collision model is not adequate to describe PT. We now present a theoretical model reflecting this bracketing procedure by explicitly introducing H-bonding complexing, dissociation and PT contributions, to discuss the kinetic model that assumes that PT occurs through one elementary step. The methods used include semiempirical density functional theory and ab initio Hartree-Fock methods. The procedure is gauged by using estimated PA of ozone obtained from deprotonation reactions including the Br?nsted bases BNH₃, H₂O, HOCl, SO₂, CH₃F and Kr. The PA-obtained range was from 145.3 to 160.3 kcal/mol, in fair agreement with the experimental value of 148.0±3 kcal/mol. The model seems to be fairly independent of the reference bases used to evaluate the PA. H-bonding effects appear to be a determining factor to explain collision efficiencies. Received: 5 August 1997 / Accepted: 25 September 1997     

惰性气体原子与HF分子相互作用的分子形貌动态研究

研究了惰性气体原子(Rg=He,Ne,Ar和Kr)分别向HF分子两端原子接近过程中分子形貌的动态变化,使用MELD精密从头算方法中的CISD/6-311++G(2d,2p)方法,计算了FH-Rg和HF-Rg2个二聚体的单电子作用势、分子平面的边界轮廓以及相应界面上的电子密度.研究表明,惰性气体原子Rg分别向HF的两端接近时,双方的外形在接近的方向上均有伸手欢迎的表示,分子的前沿电子密度也相应的发生变化,对惰性键提供了一种新的认识.     

Electrophilic fluorination of aromatic compounds with NF type reagents: kinetic isotope effects and mechanism

H/D Isotope effects in fluorination of aromatic compounds with NF type reagents have been studied to reveal the reaction mechanism. The results obtained are consistent with a polar S_EAr mechanism. Small deuterium isotope effects (k_H/k_D = 0.86-0.99) show that decomposition of a Wheland-type intermediate is not rate determining. The first example of a 1,2-hydrogen shift accompanying electrophilic fluorination of arenes has been observed in the fluorination of 1,3,5-trideuterobenzene.     

The mechanism of ethylene polymerization by nickel salicylaldiminato catalysts - Agostic interactions and their kinetic isotope effects

The ethylene polymerization reaction of a neutral nickel catalyst was studied by DFT calculations at the Becke3LYP/6-31G(d) level of theory. As in related cases a β-agostic bond stabilizes the nickel alkyl ground states. Transition states for the insertion of the olefin show a distinct α-agostic interaction, which has not been observed for late metal polymerization catalysts before. An ethylene-alkyl complex was identified as the resting state of the reaction. The overall barrier height of the reaction amounts to 17.54 kcal/mol, which slightly increases to 17.60 kcal/mol for the polymerization of deuterated ethylene. Therefore, a small positive kinetic isotope effect (k_H/k_D = 1.09) can be calculated, which is caused by the α-agostic interaction in the transition state. A comparison to other late metal based polymerization systems reveals that the ethylene coordination step of highly active catalysts is significantly lower in energy compared to catalysts which are only moderately active.