On the Gas‐Phase Reactivity of Complexed OH+ with Halogenated Alkanes

OH⁺ is an extraordinarily strong oxidant. Complexed forms (L? OH⁺), such as H₂OOH⁺, H₃NOH⁺, or iron–porphyrin‐OH⁺ are the anticipated oxidants in many chemical reactions. While these molecules are typically not stable in solution, their isolation can be achieved in the gas phase. We report a systematic survey of the influence on L on the reactivity of L? OH⁺ towards alkanes and halogenated alkanes, showing the tremendous influence of L on the reactivity of L? OH⁺. With the help of with quantum chemical calculations, detailed mechanistic insights on these very general reactions are gained. The gas‐phase pseudo‐first‐order reaction rates of H₂OOH⁺, H₃NOH⁺, and protonated 4‐picoline‐N‐oxide towards isobutane and different halogenated alkanes C_nH_2n+1Cl (n=1–4), HCF₃, CF₄, and CF₂Cl₂ have been determined by means of Fourier transform ion cyclotron resonance meaurements. Reaction rates for H₂OOH⁺ are generally fast (7.2×10^?10–3.0×10^?9 cm³ mol^?1 s^?1) and only in the cases HCF₃ and CF₄ no reactivity is observed. In contrast to this H₃NOH⁺ only reacts with tC₄H₉Cl (k_obs=9.2×10^?10), while 4‐CH₃‐C₅H₄N‐OH⁺ is completely unreactive. While H₂OOH⁺ oxidizes alkanes by an initial hydride abstraction upon formation of a carbocation, it reacts with halogenated alkanes at the chlorine atom. Two mechanistic scenarios, namely oxidation at the halogen atom or proton transfer are found. Accurate proton affinities for HOOH, NH₂OH, a series of alkanes C_nH_2n+2 (n=1–4), and halogenated alkanes C_nH_2n+1Cl (n=1–4), HCF₃, CF₄, and CF₂Cl₂, were calculated by using the G3 method and are in excellent agreement with experimental values, where available. The G3 enthalpies of reaction are also consistent with the observed products. The tendency for oxidation of alkanes by hydride abstraction is expressed in terms of G3 hydride affinities of the corresponding cationic products C_nH_2n+1⁺ (n=1–4) and C_nH_2nCl⁺ (n=1–4). The hypersurface for the reaction of H₂OOH⁺ with CH₃Cl and C₂H₅Cl was calculated at the B3 LYP, MP2, and G3_m* level, underlining the three mechanistic scenarios in which the reaction is either induced by oxidation at the hydrogen or the halogen atom, or by proton transfer.     

The effect of basis set superposition error on the convergence of interaction energies

 For the intermolecular interaction energies of ion-water clusters [OH⁻(H₂O) _n (n=1,2), F⁻(H₂O), Cl⁻(H₂O), H₃O⁺(H₂O) _n (n=1,2), and NH₄ ⁺(H₂O) _n (n=1,2)] calculated with correlation-consistent basis sets at MP2, MP4, QCISD(T), and CCSD(T) levels, the basis set superposition error is nearly zero in the complete basis set (CBS) limit. That is, the counterpoise-uncorrected intermolecular interaction energies are nearly equal to the counterpoise-corrected intermolecular interaction energies in the CBS limit. When the basis set is smaller, the counterpoise-uncorrected intermolecular interaction energies are more reliable than the counterpoise-corrected intermolecular interaction energies. The counterpoise-uncorrected intermolecular interaction energies evaluated using the MP2/aug-cc-pVDZ level is reliable. Received: 14 March 2001 / Accepted: 25 April 2001 / Published online: 9 August 2001     

Theoretical studies on structures and electronic spectra of linear carbon chains C2nH+ (n = 1−5)

The density functional theory (DFT) and the complete active space self‐consistent‐field (CASSCF) method have been used for full geometry optimization of carbon chains C_2nH⁺ (n = 1–5) in their ground states and selected excited states, respectively. Calculations show that C_2nH⁺ (n = 1–5) have stable linear structures with the ground state of X³Π for C₂H⁺ or X³Σ^? for other species. The excited‐state properties of C_2nH⁺ have been investigated by the multiconfigurational second‐order perturbation theory (CASPT2), and predicted vertical excitation energies show good agreement with the available experimental values. On the basis of our calculations, the unsolved observed bands in previous experiments have been interpreted. CASSCF/CASPT2 calculations also have been used to explore the vertical emission energy of selected low‐lying states in C_2nH⁺ (n = 1–5). Present results indicate that the predicted vertical excitation and emission energies of C_2nH⁺ have similar size dependences, and they gradually decrease as the chain size increases. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

A B3LYP study on counterpoise-corrected geometry optimizations for hydrated complexes of [K(H2O)n] and [Na(H2O)n]

We report the basis set dependencies and the basis set superposition errors for the hydrated complexes of K⁺ and Na⁺ ions in relation to the recent studies of the KcsA potassium channel. The basis set superposition errors are estimated by the geometry optimizations at the counterpoise-corrected B3LYP level. The counterpoise optimizations alter the hydration distances by about 0.02–0.03 Å. The enthalpies and free energies for K⁺ + n(H₂O) → [K(H₂O)_n]⁺ and Na⁺ + n(H₂O) → [Na(H₂O)_n]⁺ (n = 1–6) are compared between the theoretical and experimental values. The results show that the addition of diffuse functions to K, Na, and O species are effective. However, it is also found that the counterpoise corrections using diffuse functions work so as to underestimate the free energies for the complexes with increasing the hydration number. The stabilization energies in aqueous solution are larger for a Na⁺ ion than for a K⁺ ion, suggesting the contributions of their dehydration processes to the ion selectivity of the KcsA potassium channel. The changes in coordination distance between the isolated [K(H₂O)₈]⁺ and the [K(H₂O)₈]⁺ in the KcsA potassium channel indicate the importance of hydrogen bondings between the first hydration shell and the outer hydration shells.     

Monohydration selectivity mechanism of alkali cations in the potassium channel of excitable biomembranes

Nonempirical quantum chemical method Hartree–Fock–Roothan LCAO SCF MO in a two-exponent Dunning basis with the use of an extended set of Gaussian functions by Huzinaga–Dunning with consideration of electron correlation according to the Meller–Plesset theory of excitations of the second order was used to study monohydrates of Li⁺, Na⁺, K⁺, and HCOO^? ions. The indicated basis was supplemented with polarization functions of d-type on the O atom and of p-type on the hydrogen atom as well as with diffusion functions of p-type on the oxygen atom. It has been found that binding energies of the water molecule with Li⁺, Na⁺ appeared to be higher and with K⁺ lower than with HCOO^? · H₂O. Potential curve shapes of K⁺ + H₂O and HCOO^? + H₂O reactions are shown to be similar. The molecular mechanism of K⁺ channel selectivity of an excitable membrane is explained on the basis of the obtained calculations.     

Structures,Energies, and Fundamental Vibrations of the Sulfonium Ions H3Sn+ (n = 1–4)

Using ab initio MO calculations at the MP2/6‐311G(2df,2pd) level of theory the most stable structures of the following seven ions were determined: H₃S⁺ (C_3v), H₂S–SH⁺ (C_s), H₂S–S–SH⁺ (C₁), HS–S(H)–SH⁺ (C₁), H₂S–S–S–SH⁺ (C₁), HS–S(H)–S–SH⁺ (C₁) and S(SH)₃⁺ (C₃). In the case of the isomeric H₃S₃⁺ cations the species protonated at the terminal sulfur atom is most stable while in the case of the H₃S₄⁺ ions the protonation at the β sulfur atom is energetically preferred. However, the energy differences between isomeric cations are rather small. At the same level of theory the wavenumbers of the harmonic fundamental vibrations were calculated and compared to the available experimental data leading to a support for the existing assignments in certain cases but in some cases to revisions. The reaction enthalpies and Gibbs free energies of the proton transfer reactions H₂S_n + H₂S_n+1 → H₃S_n⁺ + HS_n+1^– were calculated by the G2 method. For n = 1–3 the enthalpies are found in the range 639–731 kJ mol^–1.     

Formation of H3O+ from alcohols and ethers induced by intense laser fields

The processes of H₃O⁺ production from alcohols (ethanol, 2‐propanol, 1‐propanol, 2‐butanol) and ethers (diethyl ether and ethyl methyl ether), and their deuterium‐substituted species, by intense laser fields (800 nm, 100 fs, ～1 × 10¹⁴ W/cm) were investigated through time‐of‐flight (TOF) mass spectrometry. H₃O⁺ formation was observed for all these compounds except for ethyl methyl ether. From the analysis of TOF signals of H_(3?n)D_nO⁺ (n = 0, 1, 2, and 3) that have expanding tails with increasing flight time, it has been confirmed that the reaction proceeds through metastable dissociation from the intermediate species C₂H_(5?m)D_mO⁺(m = 0–5). The common shape of the H_(3?n)D_nO⁺ signal profiles contains two major distributions in the time constant, i.e., fast and slow components of <50 ns and ～500 ns, respectively. The H_(3?n)D_nO⁺ branching ratio is interpreted to be the result of complete scrambling of four hydrogen atoms at the C? C site in C₂H₄‐OH⁺, and partial exchange (18–38%) of a hydrogen atom in the OH group with four other hydrogen atoms within 1 ns prior to H_(3?n)D_nO⁺ production. Ab initio calculations for the isomers and transition states of C₂H₅O⁺ were also performed, and the observed H_(3?n)D_nO⁺ production mechanism has been discussed. In addition, a stable isomer having a complex structure and two isomerization pathways were discovered to contribute to the H₃O⁺ formation process. Copyright © 2010 John Wiley & Sons, Ltd.     

A density functional theory study of the Cu+·(CO)n (n = 1–3) complexes

Density functional theory calculations, with an effective core potential for the copper ion, and large polarized basis set functions have been used to construct the potential energy surface of the Cu⁺·(CO)_n (n = 1–3) complexes. A linear configuration is obtained for the global minimum of the Cu⁺·CO and Cu⁺·(CO)₂ complexes with a bond dissociation energy (BDE) of 35.9 and 40.0 kcal mol^-1, respectively. For the Cu⁺·(CO)₃ complex, a trigonal planar geometry is obtained for the global minimum with a BDE of 16.5 kcal mol^?1. C-coordinated copper ion complexes exhibit stronger binding energy than O-coordinated complexes as a result of C_lp → 4s σ-donation. The computed sequential BDEs of Cu⁺·(CO)_n (n = 1–4) complexes agree well with experimental findings, in which the electrostatic energy and σ-donation play an important role in the observed trend.     

Synthesis and Crystal Structures of Alkali and Alkaline Earth Metal Complexes of 3,5-Dinitropyrazolyl-4-carboxylate

The synthesis and crystal structures of 3,5-dinitro-1H-pyrazolyl-4-carboxylic acid (H₂dnpzc) and its four complexes with Ca²⁺, Ba²⁺, Na⁺ and K⁺ are reported in this paper. Ca(dnpzc) · 5H₂O exhibits a 1D polymeric structure, whereas Ba(dnpzc) · 4H₂O possesses a 2D structure. The structure of Na₂(dnpzc) · 4H₂O consists of 2D layers of [Na(dnpzc)]^–_n and 1D chains of [Na(H₂O)₃]⁺_n. K₂(dnpzc) · H₂O has a true 3D structure. It was observed that the doubly deprotonated ligand (dnpzc^2–) can act as a versatile bridge to form polymeric structures by varying combinations of its 8 potential donor atoms (two carboxy O atoms, two pyrazolyl N atoms and four nitro O atoms). Particularly in the structure of K₂(dnpzc) · H₂O, all the 8 donor atoms of dnpzc^2– take part in the coordination and as many as 10 potassium atoms are connected by one ligand.     

Structure and energy of the positively ionized water clusters

Geometry optimization of small (H₂O)_n⁺ clusters (n ≤ 4) at the UHF/4–31 + + G^** level indicates that the cations consist of two fragments: the OH radical and the H_2n−1 O⁺_n−1 ion. The latter can be considered as a thermodynamically stable combination of a distorted H₃O⁺ ion and (n−2) H₂O molecules. The H bond between the fragments becomes weaker with increasing cluster size. Extrapolation of the adiabatic ionization potentials calculated for the (H₂O)_n oligomers (n ≤ 4) at the MP2 level to an infinite cluster size provides the value of approximately 8.7 eV, which can be presumably necessary for the ionization of liquid water in a vacuum. © 1997 John Wiley & Sons, Inc.     

Electron impact ionization of ethylene–methanol heteroclusters: Stable configurations and mechanisms in intracluster ion–molecule reactions

Reactions that proceed within mixed ethylene–methanol cluster ions were studied using an electron impact time-of-flight mass spectrometer. The ion abundance ratio, [(C₂H₄)_n(CH₃OH)_mH⁺]/[(C₂H₄)_n(CH₃OH)_m⁺], shows a propensity to increase as the ethylene/methanol mixing ratio increases, indicating that the proton is preferentially bound to a methanol molecule in the heterocluster ions. The results from isotope-labelling experiments indicate that the effective formation of a protonated heterocluster is responsible for ethylene molecules in the clusters. The observed (C₂H₄)_n(CH₃OH)_m⁺ and (C₂H₄)_n(CH₃OH)_m–1CH₃O⁺ ions are interpreted as a consequence of the ion–neutral complex and intracluster ion–molecule reaction, respectively. Experimental evidence for the stable configurations of heterocluster species is found from the distinct abundance distributions of these ions and also from the observation of fragment peaks in the mass spectra. Investigations on the relative cluster ion distribution under various conditions suggest that (C₂H₄)_n(CH₃OH)_mH⁺ ions with n + m ≤ 3 have particularly stable structures. The result is understood on the basis of ion–molecule condensation reactions, leading to the formation of fragment ions, $ {\rm CH}_2=\!=\mathop {\rm O}\limits^ + {\rm CH}_3 $ and (CH₃OH)H₃O⁺, and the effective stabilization by a polar molecule. The reaction energies of proposed mechanisms are presented for (C₂H₄)_n(CH₃OH)_mH⁺(n + m ≤ 3) using semi-empirical molecular orbital calculations.     

Unimolecular and hydrolysis channels for the detachment of water from microsolvated alkaline earth dication (Mg2+, Ca2+, Sr2+, Ba2+) clusters

We examine theoretically the three channels that are associated with the detachment of a single water molecule from the aqueous clusters of the alkaline earth dications, [M(H₂O) _n ]²⁺, M = Mg, Ca, Sr, Ba, n ≤ 6. These are the unimolecular water loss (M²⁺(H₂O) _n?1 + H₂O) and the two hydrolysis channels resulting the loss of hydronium ([MOH(H₂O) _n?2]⁺ + H₃O⁺) and Zundel ([MOH(H₂O) _n?3]⁺ + H₃O⁺(H₂O)) cations. Minimum energy paths (MEPs) corresponding to those three channels were constructed at the Møller–Plesset second order perturbation (MP2) level of theory with basis sets of double- and triple-ζ quality. We furthermore investigated the water and hydronium loss channels from the mono-hydroxide water clusters with up to four water molecules, [MOH(H₂O) _n ]⁺, 1 ≤ n ≤ 4. Our results indicate the preference of the hydronium loss and possibly the Zundel-cation loss channels for the smallest size clusters, whereas the unimolecular water loss channel is preferred for the larger ones as well as the mono-hydroxide clusters. Although the charge separation (hydronium and Zundel-cation loss) channels produce more stable products when compared to the ones for the unimolecular water loss, they also require the surmounting of high-energy barriers, a fact that makes the experimental observation of fragments related to these hydrolysis channels difficult.     

Electron Capture by the Thiyl Radical and Disulfide Bond: Ligand Effects on the Reduction Potential

The effect of non‐polar and polar ligands and of monovalent cations on the one‐electron reduction potential of the thiyl radical and the disulfide bond was evaluated. The reduction potentials E° for the CH₃S^.–n L/CH₃S^?–n L and CH₃SSCH₃–L/CH₃SSCH₃^.?–L redox couples were calculated at the B3LYP, M06‐2X and MP2 levels of theory, with n=1, 2 and L=CH₄, C₂H₄, H₂O, CH₃OH, NH₃, CH₃COOH, CH₃CONH₂, NH₄⁺, Na⁺, K⁺ and Li⁺. Non‐polar ligands decrease the E° value of the thiyl radical and disulfide bond, while neutral polar ligands favour electron uptake. Charged polar ligands and cations favour electron capture by the thiyl radical while disfavouring electron uptake by the disulfide bond. Thus, the same type of ligand can have a different effect on E° depending on the redox couple. Therefore, properties of an isolated ligand cannot uniquely determine E°. The ligand effects on E° are discussed in terms of the vertical electron affinity and reorganization energy, as well as molecular orbital theory. For a given redox couple, the ligand type influences the nature of the anion formed upon electron capture and the corresponding reorganization process towards the reduced geometry.     

Amine-Poly(Hydrogen Fluoride) Solid Complexes: New Studies of Formation,Crystal Structures,and Hn−1Fn− Ion Diversity

Triethylamine tris(hydrogen fluoride), a versatile fluorinating agent , solidifies at −27 °C forming pairs of hydrogen-bonded Et₃NH⁺ and H₂F₃⁻ ions. Six further low-melting ionic amine–HF adducts have also been identified and their structures determined. The H_n−1F_n⁻ ions with n here taking values between 3 and 7 are classified along with those observed in other crystals.     

Cationic Complexes of Hydrogen with Helium

High‐resolution mass spectra of helium nanodroplets doped with hydrogen or deuterium reveal that copious amounts of helium can be bound to H⁺, H₂⁺, H₃⁺, and larger hydrogen‐cluster ions. All conceivable He_nH_x⁺ stoichiometries are identified if their mass is below the limit of ≈120 u set by the resolution of the spectrometer. Anomalies in the ion yields of He_nH_x⁺ for x=1, 2, or 3, and n≤30 reveal particularly stable cluster ions. Our results for He_nH₁⁺ are consistent with conclusions drawn from previous experimental and theoretical studies which were limited to smaller cluster ions. The He_nH₃⁺ series exhibits a pronounced anomaly at n=12 which was outside the reliable range of earlier experiments. Contrary to findings reported for other diatomic dopant molecules, the monomer ion (i.e. H₂⁺) retains helium with much greater efficiency than hydrogen‐cluster ions.     

A comparison of the unimolecular decomposition pathways of some C7 and C8 ions in the gas phase

[CnH_2n?3]⁺ and [C_nH_2n?4]⁺·(n = 7, 8) ions have been generated in the mass spectrometer from C_nH_2n?3 Br (n = 7, 8) precursors and from two steroids. The relative abundances of competing ‘metastable transitionss’ indicate (partial) isomerization to a common structure (or mixture of structures) prior to decomposition in most examples of all four types of ions. In contrast, [C₈H₁₀O]⁺· and [C₈H₁₂O]⁺· ions, generated from different sources as molecular ions and by fragmentation of steroids, do not decompose through common-intermediates.     

Theoretical study of silicon–sulfur clusters (SiS2)n− (n=1–6)

The possible geometrical structures and relative stability of silicon–sulfur clusters (SiS₂) (n=1–6) are explored by means of density functional theory (DFT) quantum chemical calculations. We also compare DFT with second‐order Møller–Plesset (MP2) and Hartree–Fock (HF) methods. The effects of polarization functions, diffuse functions, and electron correlation are included in MP2 and B3LYP quantum chemical calculations, and B3LYP is effective in larger cluster structure optimization, so we can conclude that the DFT approach is useful in establishing trends. The electronic structures and vibrational spectra of the most stable geometrical structures of (SiS₂)_n⁻ are analyzed by B3LYP. As a result, the regularity of the (SiS₂)_n⁻ cluster growing is obtained, and the calculation may predict the formation mechanism of the (SiS₂)_n⁻ cluster. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 280–290, 2001     

The application of cholesky decomposition in valence bond calculation

The Cholesky decomposition (CD) technique, used to approximate the two‐electron repulsion integrals (ERIs), is applied to the valence bond self‐consistent field (VBSCF) method. Test calculations on ethylene, C_2nH_2n+2, and C_2nH_4n?2 molecules (n = 1–7) show that the performance of the VBSCF method is much improved using the CD technique, and thus, the integral transformation from basis functions to VB orbitals is no longer the bottleneck in VBSCF calculations. The errors of the CD‐based ERIs and of the total energy are controlled by the CD threshold, for which a value of 10^?6 ensures to control the total energy error within 10^?6 Hartree. © 2016 Wiley Periodicals, Inc.     

Analytical TDHF second derivatives of dynamic electronic polarizability with respect to nuclear coordinates. Application to the dynamic ZPVA correction

Analytical calculation of second derivatives of the dynamic electronic polarizability with respect to vibrational normal coordinates has been implemented at the TDHF level of approximation. Important simplifications are obtained by using the 2n+1 rule and interchange relations. In application to small molecules the first‐order ZPVA correction is found to have greater frequency dispersion than the electronic polarizability, but the former remains small in comparison. For the linear polyenes, C_2nH_2n+2 (n=1–3), the percentage ZPVA correction decreases with chain length for all ?ω<0.10 a.u. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1920–1932, 2001