Ab initio studies of negative ion-molecule(s) clusters present in the atmosphere. II. OH−(H2O)n for n = 0, 2, 3, 4

Ab initio calculations are performed with 6–31G basis set to study the geometry and binding of the H₃O, H₅O, H₇O, and H₉O complexes. The H₃O complex is also investigated with the 6–31 G* basis set and MP 2 (Moller–Plesset perturbation theory of second order).     

The structures,thermochemistry, and electron affinities of hydrogenated silicon clusters Si6Hn/Si6H (n = 3–14)

The hydrogenated silicon clusters structures, electron affinities, and dissociation energies of the Si₆H_n/Si₆H (n = 3?14) species have been systematically investigated by means of three density functional theory (DFT) methods. The basis set used in this work is of double‐ζ plus polarization quality with additional diffuse s‐ and p‐type functions, denoted DZP++. The geometries are fully optimized with each DFT method independently. Three different types of energy separations presented in this work are the adiabatic electron affinity (EA_ad), the vertical electron affinity (EA_vert), and the vertical detachment energy (VDE). The first Si? H dissociation energies D_e (Si₆H_n→ Si₆H_n?1+H) for the neutral Si₆H_n and D_e (Si₆H→Si₆H+H) for the anionic Si₆H species have also been reported. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Aromaticity of ionic structures: Investigation and application of NICS value and 4n + 2 rule

At DFT/B3LYP/6‐31G** theoretical level, C₆H and C (n = 0, ?2, and +2), C₆H and C (n = 0, ±2, ±4, and ±6), C₆H (n = 0–6), as well as C₆H₆‐A and C₆‐A (A = Be, B, N, O, Mg, Al, Si, S, and Fe) structures were investigated. Comparing NICS values of C₆H and C (n = 0, ?2, and +2), we discovered that C₆H, C₆H were antiaromatic, and C₆H₆, C₆, C, C had aromaticity with negative NICS values. According to research of C₆H and C (n = 0, ±2, ±4, ±6), C₆H (n = 0–6), we sustained that their σ and π orbit were different and the locations of electrons were difficult to confirm in ionic structures. Thus, neither 4n + 2 rule nor NICS values can precisely estimate the aromaticity of ionic structures. Besides, through WBI (NBO) research of C₆H₆‐A and C₆‐A (A = Be, B, N, O, Mg, Al, Si, S, and Fe) structures, we found that C₆H₆ was easy to accept electrons, contrarily, C₆ was prone to bestowing electrons. Moreover, C₆H₆ took the symmetrical carbon atoms form feeble interaction or bond, and C₆ used all carbon atoms to impact with other atom. C₆H₆ generated two contrapuntal single bonds with oxygen, sulfur, and nitrogen atoms, whereas C₆ molecule formed double bond with oxygen and nitrogen atoms, two conjoint single bonds with sulfur atom. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Structure and dissociation energy of weakly bound H (n = 5−8) complexes

The geometrical parameters, vibrational frequencies, and dissociation energies for H (n = 5–8) clusters have been investigated using high level ab initio quantum mechanical techniques with large basis sets. The highest level of theory employed in this study is TZ2P CCSD(T). The C₁ structure of H is predicted to be a global minimum, while the C_s structure of H is calculated to be a transition state. Harmonic vibrational frequencies are also determined at the DZP and TZ2P CCSD levels of theory. The dissociation energies, D_e, for H (n = 5–8) have been predicted using energy differences at each optimized geometry, and zero‐point vibrational energies (ZPVEs) are considered to compare with experimental values. The dissociation energies (D_o) have been predicted to be 1.69, 1.65, 1.65, and 1.46 kcal · mol for H, H, H (C₁ symmetry) and H, respectively, at the TZ2P CCSD(T) level of theory. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Mechanistic investigation of vibrational fine structure in e‐H2 scattering using local complex potential‐based time dependent wave packet approach

The structural features of vibrational excitation cross‐sections in resonant e‐H₂ scattering have been investigated using a time dependent wave packet approach and a local complex potential to describe the ²Σ H anion. An analysis of the partial contributions to the vibrational excitation cross‐sections reveals that all features of the excitation profile result from simple interference between bound vibrational levels of H₂ and H. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

The adsorption of CO2, H2CO3, HCO3− and CO32− on Cu2O (111) surface: First‐principles study

The adsorption of CO₂, and its derivatives, H₂CO₃, HCO, and CO, on Cu₂O (111) surface has been investigated by first‐principles calculations based on the density functional theory at B3LYP hybrid functional level. The Cu₂O (111) surface has been modeled using an embedded cluster method,in which the quantum clusters plus some ab initio ion model potentials were inserted in an array of point charges. On the surface, H₂CO₃ was dissociated into an H⁺ and an HCO ion. Among the CO₂ species, HCO was the only activated species on the surface. The results suggest that the reduction of CO₂ on Cu₂O (111) surface can start from the form of HCO. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Reaction mechanism of platinum dimer cation with ammonia based on the relativistic density functional study

The gas‐phase reactions between Pt and NH₃ have been investigated using the relativistic density functional approach (ZORA‐PW91/TZ2P). The quartet and doublet potential energy surfaces of Pt + NH₃ have been explored. The minimum energy reaction path proceeds through the following steps: Pt(⁴Σ_u) + NH₃ → q‐1 → d‐2 → d‐3 → d‐4 → d‐Pt₂NH⁺ + H₂. In the whole reaction pathway, the step of d‐2 → d‐3 is the rate‐determining step with a energy barrier of 36.1 kcal/mol, and exoergicity of the whole reaction is 12.0 kcal/mol. When Pt₂NH⁺ reacts with NH₃ again, there are two rival reaction paths in the doublet state. One is degradation of NH and another is loss of H₂. In the case of degradation of NH, the activation energy is only 3.4 kcal/mol, and the overall reaction is exothermic by 8.9 kcal/mol. Thus, this reaction is favored both thermodynamically and kinetically. However, in the case of loss of H₂, the rate‐determining step's energy barrier is 64.3 kcal/mol and the overall reaction is endothermic by 8.5 kcal/mol, so it is difficult to take place. Predicted relative energies and barriers along the suggested reaction paths are in reasonable agreement with experimental observations. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Theoretical study of one‐electron bonds in a series of high‐spin lithium–beryllium–hydrogen clusters: “Valence shell single‐electron repulsion” rule and electron localization function analysis

A series of high‐spin clusters containing Li, H, and Be in which the valence shell molecular orbitals (MOs) are occupied by a single electron has been characterized using ab initio and density functional theory (DFT) calculations. A first type (⁵Li₂, ⁿ⁺¹LiH_n+ (n = 2–5), ⁸Li₂H) possesses only one electron pair in the lowest MO, with bond energies of ～3 kcal/mol. In a second type, all the MOs are singly occupied, which results in highly excited species that nevertheless constitute a marked minimum on their potential energy surface (PES). Thus, it is possible to design a larger panel of structures (⁸LiBe, ⁷Li₂, ⁸Li, ⁴LiH⁺, ⁶BeH, ⁿ⁺³LiH (n = 3, 4), ⁿ⁺²LiH (n = 4–6), ⁸Li₂H, ⁹Li₂H, ²²Li₃Be₃ and ²²Li₆H), single‐electron equivalent to doublet “classical” molecules ranging from CO to C₆H₆. The geometrical structure is studied in relation to the valence shell single‐electron repulsion (VSEPR) theory and the electron localization function (ELF) is analyzed, revealing a striking similarity with the corresponding structure having paired electrons. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

D∞h B2(BS)2−/2− and Td B(BS)4−: Boron sulfide clusters containing BB multiple bonds and B− tetrahedral centers

A density functional theory investigation on the geometrical and electronic properties of B₄S (B₂(BS)) and B₅S (B(BS)) clusters has been performed in this work. Both the doublet B₂(BS) ([S?B? BB? B?S]^?) (D_∞h, ²Π_u) and the singlet B₂(BS) ([S?B? B?B? B?S]^2?) (D_∞h, ¹Σ) proved to have perfect linear ground‐state structures containing a multiply bonded BB core (BB or B?B) terminated with two BS groups, while T_d B(BS) turned out to possess a perfect B^? tetrahedral center directly corrected to four BS groups, similar to the corresponding boron hydride molecules of D_∞h B₂H, D_∞h B₂H, and T_d BH, respectively. B₄S₂ and B₅S₄ neutrals, however, appeared to be much different: they favor a planar fan‐shaped C_2v B₄S₂ (a di‐S‐bridged B₄ rhombus) and a planar kite‐like C_2v B₅S₄ (a di‐S‐bridged B₃ triangle bonded to two BS groups), respectively. One‐electron detachment energies and symmetrical stretching vibrational frequencies are calculated for D_∞h B₂(BS) and T_d B(BS) monoanions to facilitate their future characterizations. Neutral salts of B₂(BS)₂Li₂ with an elusive B?B triple bond and B(BS)₄Li containing a tetrahedral B^? center are predicted possible to be targeted in experiments. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Leading‐order behavior of the correlation energy in the uniform electron gas

We show that, in the high‐density limit, restricted Møller‐Plesset (RMP) perturbation theory yields E = π^?2(1 ? ln 2) ln r_s + O(r) for the correlation energy per electron in the uniform electron gas, where r_s is the Seitz radius. This contradicts an earlier derivation which yielded E = O(ln|ln r_s|). The reason for the discrepancy is explained. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Theoretical study on the gas‐phase reactivity of halogenated alkylperoxyl radicals toward alkenes

The electrophilic additions of hydroperoxyl (HO ${}_{\text{2}}^{\text{⋅}}$), alkylperoxyl (RO ${}_{\text{2}}^{\text{⋅}}$), and halogenated alkylperoxyl radicals to ethylene were studied using the AM1 and PM3 semiempirical MO methods at the SCF/UHF level. Reactantlike transition states were predicted for the title additions. The AM1 activation enthalpies (ΔH ${}_{\text{f}}^{\text{*}}$) were found to be increased in the order HO ${}_{\text{2}}^{\text{⋅}}$<CH₃O ${}_{\text{2}}^{\text{⋅}}$<C₂H₅O ${}_{\text{2}}^{\text{⋅}}$<i‐C₃H₇O ${}_{\text{2}}^{\text{⋅}}$. The reactivity of an alkylperoxyl radical toward ethylene was found to be increased as the degree of halogen substitution on the alkyl group increased. A good correlation was established between ΔH ${}_{\text{f}}^{\text{*}}$ and the Taft polar substituent constants, σ*. The Evans–Polanyi correlation between ΔH ${}_{\text{f}}^{\text{*}}$ and ΔH ${}_{\text{r}}^{\text{°}}$ was justified and the validity of the Hammond postulate was indicated. The calculated results were compared with the available experimental findings. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 71: 273–283, 1999     

Ab initio calculation of electronic states of the ions H4+ and H5+

Electronic state calculations for the ions H₄⁺ (with symmetries D ₄ and C _2v) and H (with symmetries D ₅ and D _2d) are made using the valence-bond method. All the configurations obtained from the given set of 1s-functions of Slater type are taken into account. Space functions are used throughout the computation (“spin-free quantum chemistry”). Preliminary quasidiagonalization of the secular equation is implemented by the construction of the multiplet eigenfunctions ^2S+1Γ^(α) from the initial variational functions. The results of the calculations are as follows: the ion H is unstable, the ion H is stable with equilibrium nuclear conformation of symmetry D _2d and with the energy of dissociation into H and H₂ near 4 eV.     

The formation of highly excited H in the reaction H2+(v) + H2 → H + H

A quasiclassical trajectory surface hopping method has been used to study H(v) + H₂ → H + H for v = 0, 3, 7, 10, 13, and 17 with an emphasis on determining the H internal energy and angular momentum distributions for high v. For v = 13 and 17, significant cross sections are found for producing H at energies above its dissociation energy. An average metastable H lifetime of 11.5 ps for v = 13 and 4.7 ps for v = 17 is found, but there is also a much longer lived component to the lifetime distributions that is more important for v = 13 than for v = 17. Some of the longer lived metastables correspond to high angular momentum orbiting states of H, but other sources of metastability are also present.     

The role of d orbitals in tetrahedral hybrid orbitals of oxy-ions

The hybrid orbitals of tetrahedral oxy-ions containing some d character have been calculated by maximum overlap method. The d characters of hybrid orbitals increase in the order of SiO, PO, SO, ClO, and decrease in order of GeO, AsO, SeO, BrO. The bond strengths are also obtained for these ions. The hybrid Orbital of VO, CrO, and MnO are of the type d³s as the result of calculation.     

Dipole,quadrupole, octupole,and dipole–octupole polarizabilities at real and imaginary frequencies for H,HE, and H2 and the dispersion-energy coefficients for interactions between them

We have calculated certain dynamic polarizabilities (for both real and imaginary frequencies) for H, He, and H₂ and the dispersion-energy coefficients for long-range interactions between them. We have done so in a sum-over-states formalism with explicitly electron-correlated wave functions to describe the states. To be precise, we have determined the dipole (α₁), quadrupole (α₂), and octupole (α₃) polarizabilities of H and He for real frequencies (ω) in a range between zero and the first electronic-transition frequency and for imaginary frequencies (iω) on a 32-point Gauss-Legendre grid running from zero to ?ω = 20 E_h, and for H₂, we have found the dipole (α), quadrupole (C), and dipole–octupole (E) polarizability tensors for the same real and imaginary frequencies. The dispersion-energy coefficients, obtained by combining the sum-over-states for-malism for the polarizabilities with analytic integration over ω, gave values of C₆, C₈, and C₁₀ for the atom–atom systems; C, C, C, C, and C for the atom–diatom systems; and C₆, C and C for the H₂? H₂ system. Nearly all the results are considered to be more reliable than those hitherto published and some have been obtained for the first time, e.g., C(iω), E(ω), and E(iω) for H₂ and C, C, and C for the H? H₂ system. © 1993 John Wiley & Sons, Inc.     

Variation of the polarizability of the hydrogen molecule ion and the hydrogen molecule with internuclear separation

The variation of the polarizability of H and H₂ with internuclear separation R = 1.6 – R = 2.4 a.u. for H and R = 1.0 – R = 2.0 a.u. for H₂ is determined using a variational method suggested by Das and Bersohn. From these data, values of 〈α〉_0,J for which nuclear motion due to zero point vibration and centrifugal stretching is taken into account, are calculated at 300°K. The relative percent increases of the motion averaged values compared to the equilibrium values are as follows: 10.50% for H and 6.52% for H₂.     

Structure and stability of B6, B,and B clusters

We investigated various isomers of B₆, B, and B clusters with ab initio [Hartree–Fock (HF), MP2)] and density functional theory (DFT) methods. Ten B₆ isomers, 6 B isomers, and 6 B isomers are determined to be local minima on their potential energy hypersurfaces by the HF, B3LYP, B3PW91, and MP2 methods. Fourteen of these structures are first reported. The most stable neutral B₆ cluster is the capped pentagonal pyramid (C_5v), in agreement with the results reported previously. Hexagon B (C_2h) isomer and fan‐shaped B (C_2v) isomer are found to be the most stable on the cationic and anionic energy hypersurfaces, respectively. Natural bond orbital analysis suggests that there are three‐centered bonds in the most stable B₆ neutral and ionic clusters. The multicentered bonds are responsible for the special stability of the lowest‐energy isomer. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem 94: 269–278, 2003     

Role of angular electron pair correlation in stabilizing C

A pseudo‐potential that was successfully employed in an earlier study by the Compton group is used to describe the binding of a single electron to a C₆₀ molecule to form C. Then, the interaction of a second electron with the C anion is treated in two manners. First, as performed in the earlier Compton study, a mean‐field (i.e., Hartree–Fock) approach is used to estimate the C‐to‐C energy difference for the singlet state of the dianion and, much as in the earlier study, this dianion is predicted to be unstable by ∼0.4 eV. Second, for this same singlet state, a configuration interaction wave function is employed that allows for the angular correlation of the two excess electrons, allowing them to avoid one another by moving on opposite sides of the C₆₀ skeleton. The energy of the dianion is lowered by 0.3 eV when angular correlation is included, suggesting that the singlet dianion is unstable with respect to electron loss by only ∼0.1 eV. A Coulomb barrier (>1 eV high) and angular momentum barriers then combine to trap electrons of singlet C from detaching, thus producing the very long observed lifetimes. In addition, the energy of the lowest triplet state of C is also discussed. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Convergence accelerator approach for the evaluation of some three‐electron integrals containing explicit rij factors

A simplified analysis is presented for the evaluation of the three‐electron one‐center integrals of the form ∫rrrrrred r ₁d r ₂d r ₃, for the cases i, j, k, ≥−2, l=−2, m≥−1, n≥−1. These integrals arise in the calculation of lower bounds for energy levels and certain relativistic corrections to the energy when Hylleraas‐type basis sets are employed. Convergence accelerator techniques are employed to obtain a reasonable number of digits of precision, without excessive CPU requirements. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 72: 93–99, 1999     

Negative ion fragmentations of deprotonated peptides containing post‐translational modifications. An unusual cyclisation/rearrangement involving phosphotyrosine; a joint experimental and theoretical study

The characteristic fragmentations of a pTyr group in the negative ion electrospray mass spectrum of the [M–H]^? anion of a peptide or protein involve the formation of PO (m/z 79) and the corresponding [(M‐H)^?–HPO₃]^? species. In some tetrapeptides where pTyr is the third residue, these characteristic anion fragmentations are accompanied by ions corresponding to H₂PO and [(M‐H)^?–H₃PO₄]^? (these are fragmentations normally indicating the presence of pSer or pThr). These product ions are formed by rearrangement processes which involve initial nucleophilic attack of a C‐terminal ‐CO [or ‐C(?NH)O^?] group at the phosphorus of the Tyr side chain [an S_N2(P) reaction]. The rearrangement reactions have been studied by ab initio calculations at the HF/6‐31+G(d)//AM1 level of theory. The study suggests the possibility of two processes following the initial S_N2(P) reaction. In the rearrangement (involving a C‐terminal carboxylate anion) with the lower energy reaction profile, the formation of the H₂PO and [(M‐H)^?–H₃PO₄]^? anions is endothermic by 180 and 318 kJ mol^?1, respectively, with a maximum barrier (to a transition state) of 229 kJ mol^?1. The energy required to form H₂PO by this rearrangement process is (i) more than that necessary to effect the characteristic formation of PO from pTyr, but (ii) comparable with that required to effect the characteristic α, β and γ backbone cleavages of peptide negative ions. Copyright © 2009 John Wiley & Sons, Ltd.