Regioselective Oxidative Ring Opening of Cyclopropyl Silyl Ethers: A Quantum Chemical Study

In contrast to the structurally and configurationally stable alkyl‐ or aryl‐substituted cyclopropyl radical cations, cyclopropyl silyl ethers undergo spontaneous ring opening upon oxidation whereby the endocyclic C? C(O‐TMS) bond is cleaved with remarkable selectivity. DFT calculations on 1‐trimethylsilyloxybicyclo[4.1.0]heptane show that this selectivity arises from the topology of the potential surface of the corresponding radical cation which is initially generated in a very steep region of the potential surface from where the steepest descent leads to cleavage of the endocyclic rather than the lateral C? C(OTMS) bond. Cleavage of the lateral bond leads to interesting conformational changes which are explored in detail.     

Ionic liquid salt bridge based on tributyl(2-methoxyethyl)phosphonium bis(pentafluoroethanesulfonyl)amide for stable liquid junction potentials in highly diluted aqueous electrolyte solutions

A moderately hydrophobic ionic liquid, tributyl(2-methoxyethyl)phosphonium bis(pentafluoroethanesulfonyl)amide ([TBMOEP⁺][C₂C₂N⁻]), shows a very stable liquid junction potential upon contact with an aqueous solution whose ionic strength is as low as 1 μmol dm⁻³. The stability with the maximum excursion of the potential within ±0.5 mV for 30 min is very promising for accurate determination of pH and other single ion activities potentiometrically.     

Theoretical study of mechanism of cycloaddition reaction between dimethylmethylene carbene and formaldehyde

The cycloaddition mechanism of forming a polycyclic compound between singlet dimethylmethylene carbene(R1) and formaldehyde(R2) has been investigated with MP2/6‐31G* method, including geometry optimization and vibrational analysis for the involved stationary points on the potential energy surface. The energies of the different conformations are calculated with CCSD(T)//MP2/6‐31G* method. From the potential energy profile, it can be predicted that the dominant reaction pathway of the cycloadditional reaction between singlet dimethylmethylene carbene and formaldehyde consists of two steps: (1) the two reactants(R1, R2) firstly form an energy‐enricheded intermediate (INT1a) through a barrier‐free exothermic reaction of ΔE = 11.3 kJ/mol. (2) Intermediate (INT1a) then isomerizes to a three‐membered product (P1) via a transition state (TS1a) with an energy barrier of 20.0 kJ/mol. The dominant reaction has an excellent selectivity and differs considerably from its competitive reactions in reaction rate. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Theoretical infrared line shapes of H‐bonds within the strong anharmonic coupling theory and Fermi resonances effects

In this article, we extend a previous work toward presenting a theoretical study of the effects of Fermi resonances and the fundamental anharmonic coupling parameter α between the high‐frequency mode and the H‐bond bridge. The model incorporates (i) both intrinsic anharmonicities of the fast mode (double well potential) and the H‐bond Bridge (Morse potential), (ii) strong anharmonic coupling theory, (iii) Fermi resonances by the aid of an anharmonic coupling between the fast mode and one or several harmonic bending modes, (iv) quadratic modulation of both the angular frequency and the equilibrium position of the X? …Y stretching mode on the intermonomer ? H… motions, and (v) the quantum direct (fast and bending modes) and indirect dampings (slow mode). The IR spectral density is obtained by Fourier transform of the autocorrelation function of the transition dipole moment operator of the X? H bond. The numerical calculation shows that Fermi resonances generate very complicated profiles with multisubstructure and also provide a direct evidence of Fermi resonances which were predicted to be a major feature of H‐bonds. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Theoretical prediction on HBeN− and HNBe− anions using multiconfiguration second‐order perturbation theory

The HBeN^? and HNBe^? anions have been investigated for the first time using the CASSCF, CASPT2, and DFT/B3LYP methods with the contracted atomic natural orbital (ANO) and cc‐pVTZ basis sets. The geometries of all stationary points along the potential energy surfaces were optimized at the CASSCF/ANO and B3LYP/cc‐pVTZ levels. The ground and the first excited states of HBeN^? are predicted to be X²Π and A²Σ⁺ states, respectively. It was predicted that the ground state of HNBe^? is X²Σ⁺ state. The A²Π state of HNBe^? has unique imaginary frequency. A bend local minimum M1 was found along the 1²A″ potential energy surface and the A²Π state of HNBe^? should be the transition state of the isomerization reactions for M1 ? M1. The CASPT2/ANO potential energy curves of isomerization reactions were calculated as a function of HBeN bond angle. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Structural and electronic study of neutral,positive, and negative small rhodium clusters [Rhn,Rh, and Rh]

We have carried out a systematic study for the determination of the structure and the fundamental state of neutral and ionic small rhodium clusters [Rh_n, Rh, Rh (n = 5–10)] using ab initio Hartree‐Fock methods with a LANL2DZ basis set. A range of spin multiplicities is studied for each cluster. We present the bond lengths and angles and geometric configuration adopted for the clusters in its minimum energy conformation showing the differences when the clusters have different number of unpaired electrons. In addition, we report the vertical ionization potential and the adiabatic potential calculated by the Koopman's theorem. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem 110:2541–2547, 2010     

Potentiometric determination of coextraction constants of potassium salts in ion-selective electrodes utilizing a nitrobenzene liquid membrane phase

A theoretical treatment of potentiometric data is applied to calculate coextraction constants (K_IA) for three potassium salts from water into a liquid nitrobenzene phase. The experiment involves treating nitrobenzene as a membrane and contacting it with two aqueous solutions of different ion activities. In the presence of either a cation or anion exchanger, the ratio of activities of ions in the two aqueous phases gives rise to a potential difference across the membrane that depends upon the nature and charge of the counter ion of the ion-exchanger in excess. Here, the cation exchanger was chosen to be potassium tetrakis(4-chlorophenyl)borate (KTpClPB) and the anion exchanger was tetradodecylammonium chloride (TDDACl). TDDACl was incrementally added to the nitrobenzene phase containing a fixed concentration of KTpClPB, and the corresponding emf was recorded as a function of concentration of TDDACl. The membrane changes from one with cation exchanger properties (excess KTpClPB) to one with anion exchanger properties (excess TDDACl). The potential difference and shape of the titration curve can be predicted by theory based on the phase boundary potential model. Log(K_IA) values calculated for KCl, KNO₃ and KClO₄ in nitrobenzene were found as: −10.53 (± 0.09), −8.16 (± 0.05) and −5.63 (± 0.03) respectively, in accordance with the Hofmeister series of lipophilicity, and similar to those observed in PVC membranes containing other plasticizers. The method presented here offers the advantage over other methods to calculate K_IA, in that it is relatively experimentally simple without compromising the accuracy of the calculated coextraction constants. The ability to titrate directly into the liquid membrane phase affords a higher precision compared to the preparation of a series of PVC/plasticizer membranes with different compositions.     

New angle‐dependent potential energy function for backbone–backbone hydrogen bond in protein–protein interactions

Backbone–backbone hydrogen bonds (BBHBs) are one of the most abundant interactions at the interface of protein–protein complex. Here, we propose an angle‐dependent potential energy function for BBHB based on density functional theory (DFT) calculations and the operation of a genetic algorithm to find the optimal parameters in the potential energy function. The angular part of the energy funtion is assumed to be the product of the power series of sine and cosine functions with respect to the two angles associated with BBHB. Two radial functions are taken into account in this study: Morse and Leonard‐Jones 12‐10 potential functions. Of these two functions under consideration, the former is found to be more accurate than the latter in terms of predicting the binding energies obtained from DFT calculations. The new HB potential function also compares well with the knowledge‐based potential derived by applying Boltzmann statistics for a variety of protein–protein complexes in protein data bank. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

An ab initio potential energy surface and vibrational energy levels of ZnH2

A three‐dimensional potential energy surface of the electronic ground state of ZnH₂ (${X}^1\sum _g^ +$ ) molecule is constructed from more than 7500 ab initio points calculated at the internally contracted multireference configuration interaction with the Davidson correction (icMRCI+Q) level employing large basis sets. The calculated relative energies of various dissociation reactions are in good agreement with the previous theoretical/experimental values. Low‐lying vibrational energy levels of ZnH₂, ZnD₂, and HZnD are calculated on the three‐dimensional potential energy surface using the Lanczos algorithm, and found to be in good agreement with the available experimental band origins and the previous theoretical values. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Accurate free energy calculation along optimized paths

The path‐based methods of free energy calculation, such as thermodynamic integration and free energy perturbation, are simple in theory, but difficult in practice because in most cases smooth paths do not exist, especially for large molecules. In this article, we present a novel method to build the transition path of a peptide. We use harmonic potentials to restrain its nonhydrogen atom dihedrals in the initial state and set the equilibrium angles of the potentials as those in the final state. Through a series of steps of geometrical optimization, we can construct a smooth and short path from the initial state to the final state. This path can be used to calculate free energy difference. To validate this method, we apply it to a small 10‐ALA peptide and find that the calculated free energy changes in helix‐helix and helix‐hairpin transitions are both self‐convergent and cross‐convergent. We also calculate the free energy differences between different stable states of β‐hairpin trpzip2, and the results show that this method is more efficient than the conventional molecular dynamics method in accurate free energy calculation. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010