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     

Evaluation of two‐center,three‐ and four‐electron integrals over Slater‐type orbitals in elliptical coordinates

An algorithm for evaluation of two‐center, three‐electron integrals with the correlation factors of the type rr and rrr as well as four‐electron integrals with the correlation factors rrr and rrr in the Slater basis is presented. This problem has been solved here in elliptical coordinates, using the generalized and modified form of the Neumann expansion of the interelectronic distance function r for k ≥ ?1. Some numerical results are also included. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

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     

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     

Correlated n1,3S states for coulomb three‐body systems

n^1,3S (n = 1 ? 4) states for atomic three‐body systems are studied with the Angular Correlated Configuration Interaction method. A recently proposed angularly correlated basis set is used to construct, simultaneously and with a single diagonalization, ground and excited states wave functions which: (i) satisfy exactly Kato cusp conditions at the two‐body coalescence points; (ii) involve only linear parameters; (iii) show a fast convergency rate for the energy; and (iv) form an orthogonal set. The efficiency of the method is illustrated by the study a variety of three‐body atomic systems [m m m] with two negatively charged light particles, with diverse masses m and m, and a heavy positively charged nucleus m. The calculated ground 1¹S and excited n^1,3S (n = 2 ? 4) state energies are compared with those given in the literature, when available. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

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     

New implementation of molecular double point‐group symmetry in four‐component relativistic Gaussian‐type spinors

A new, practical implementation of double‐group symmetry to relativistic Gaussian spinors is presented for four‐component relativistic molecular calculations. We show that the systematic adaptability to irreducible representations under arbitrary point‐group symmetry, as well as Kramers (time‐reversal) symmetry, is inherent in the present basis spinors, which possess the analytic structure of Dirac atomic spinors. The implementation of double‐group symmetry entails significant computational efficiencies in the relativistic second‐order Møller–Plesset perturbation calculation on Au₂ and the density functional theory (DFT) calculation with the B3LYP functional on octahedral UF₆, in which the highest symmetries used are, respectively, C and D. The four‐component B3LYP equilibrium geometry of UF₆ is reported. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Maximum carbonyl‐coordination number of scandium Computational study of Sc(CO)n (n = 1–7), Sc(CO)7− and Sc(CO)63−

Understanding the maximum bonding ability is very important with the potential both to design new compounds and to broaden chemists' imagination. While the coordination ability of the late transition metals has been richly understood, that of scandium is very poor. In this work, a detailed computational study on the equilibrium geometries, stability and vibrational frequencies of a series of Sc(CO)_n (n = 1–7), Sc(CO) and Sc(CO) is reported using density functional theory functionals and the coupled cluster (single‐point) method with 6‐311+G(3df) basis set. It was shown that the obtained sequential and average CO binding energies of Sc(CO)_n (n = 4–7), Sc(CO) and Sc(CO) are comparable to those of the experimentally known species, i.e., smaller Sc‐carbonyls (n ≤3) and the analog Ti(CO)₇⁺. Thus, the studied high scandium carbonyls could all be experimentally accessible. In addition, the studied Sc(CO)_n generally favor the low‐spin ground state (doublet) structures except ScCO and Sc(CO)₃ that are in the quartet states. The previously uncertain spectrum bands were assigned to Sc(CO)₄ and Sc(CO)₅ in this work. In all, the appreciable stability suggested that the last 18‐electron first‐row transition metal carbonyls, that is, Sc(CO) and Sc(CO), could be accessible in experiment. © 2013 Wiley Periodicals, Inc.     

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     

Probing the ultrafast photoelectron spectra of Br2 molecule

The time‐dependent‐wave‐packet method is applied to study the ionization of Br₂ molecule with four ionization processes. The ground state absorption makes the photoelectron to be left in the three final ionic states: Br (X²∑), Br (A²∏_u), and Br (B²∑), and each population of these ionic states is related with the laser intensities. The information of the dissociation can be got by analyzing the photoelectron features of the transient wave packet, which also suggests that an ionization process occurs during the dissociation, and the Br atoms that mainly resulted from the dissociation of Br₂ (C¹∏_u) are ionized at later time delays as the dissociation is nearly complete. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Structure and stability of high‐spin Aun(n = 2–8) clusters

The structures and relative stability of the maximum‐spin ⁿ⁺¹Au_n and ⁿAu (n = 2–8) clusters have been determined by density‐functional theory. The structure optimizations and vibrational frequency analysis are performed with the gradient‐corrections of Perdew along with his 1981 local correlation functional, combined with SBKJC effective core potential, augmented in the valence basis set by a set of f functions. We predicted the existence of a number of previously unknown isomers. The energetic and electronic properties of the small high‐spin gold clusters are strongly dependent on sizes. The high‐spin clusters tend to holding three‐dimensional geometry rather than planar form preferred in low‐spin situations. In whole high‐spin Au_n (n = 2–8) neutral and cationic species, ⁵Au₄, ²Au, and ⁴Au are predicted to be of high stability, which can be explained by valence bond theory. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Spectroscopic parameter and molecular constant investigations for low‐lying electronic states of P ion

The potential energy curves (PECs) of three low‐lying electronic states of P ion, X²Π_u, A²Σ, and B²Σ, have been studied using the full valence complete active space self‐consistent field method followed by the highly accurate valence internally contracted multireference configuration interaction (MRCI) approach and MRCI with Davidson correction (+Q). The correlation‐consistent basis sets, aug‐cc‐pV5Z and aug‐cc‐pV6Z, are used and the total energies are extrapolated to the complete basis set limit. Using these PECs obtained with the MRCI+Q/56‐extrapolation, the spectroscopic parameters for these electronic states are determined and compared in detail with experimental data and those of previous studies reported in the literature. The comparison shows that excellent agreement exists between the present results and the available experiments. The first 40 vibrational states for the three electronic states are also computed when the rotational quantum number J equals zero. For each vibrational state, the vibrational level G(υ), inertial rotation constant B_υ, and centrifugal distortion constant D_υ are determined when J = 0, which are in good accord with the available measurements. © 2012 Wiley Periodicals, Inc.     

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     

Theoretical study of the gas‐phase SN2 reactions of X− with CH3OY (X,Y = Cl,Br, I)

The gas‐phase nucleophilic substitution reactions at saturated oxygen X^? + CH₃OY (X, Y = Cl, Br, I) have been investigated at the level of CCSD(T)/6‐311+G(2df,p)//B3LYP/6‐311+G(2df,p). The calculated results indicate that X^? preferably attacks oxygen atom of CH₃OY via a S_N2 pathway. The central barriers and overall barriers are respectively in good agreement with both the predictions of Marcus equation and its modification, respectively. Central barrier heights (ΔH and ΔH) correlate well with the charges (Q) of the leaving groups (Y), Wiberg bond orders (BO) and the elongation of the bonds (O? Y and O? X) in the transition structures. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

An Ab initio theoretical investigation on the geometrical and electronic structures of BAun−/0 (n = 1–4) clusters

An ab initio theoretical investigation on the geometrical and electronic structures and photoelectron spectroscopies (PES) of BAu_n^?/0 (n = 1–4) auroboranes has been performed in this work. Density functional theory and coupled cluster method (CCSD(T)) calculations indicate that BAu (n = 1–4) clusters with n‐Au terminals possess similar geometrical structures and bonding patterns with the corresponding boron hydrides BH. The PES spectra of BAu (n = 1–4) anions have been simulated computationally to facilitate their future experimental characterizations. In this series, the T_d BAu anion appears to be unique and particularly interesting: it possesses a perfect tetrahedral geometry and has the highest vertical electron detachment energy (VDE = 3.69 eV), largest HOMO‐LUMO gap (ΔE_gap = 3.0 eV), and the highest first excitation energy (E_ex = 2.18 eV). The possibility to use the tetrahedral BAu unit as the building block of Li⁺[BAu₄]^? ion‐pair and other [BAu₄]^?‐containing inorganic solids is discussed. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

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     

Field‐assisted dissociative ionization of CH2I2 induced by femtosecond laser field

The field‐assisted dissociative ionization of CH₂I irradiated by a 60‐fs 800‐nm laser with different laser intensities (1–4 × 10¹⁴ W/cm²) is studied both experimentally and theoretically. The different fragmentation patterns are observed in the experiment with a time‐of‐flight mass spectrometer. In the theoretical aspect, the Gaussian 03 program is applied to calculate the potential energies of CH₂I as functions of the C? I and C? H bond distances and I? C? I bond angle under external field with different intensities. The calculations explain our experimental observations well. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

An insight into optical and EPR properties of AgCl and AgF complexes through MS–Xα and SCCEH calculations

MS-Xα and SCCEH calculations on the Ag²⁺ complexes AgF and AgCl (displaying an elongated D_4h symmetry) have been carried out for a better understanding of their experimental optical and EPR properties. As salient features, the present work supports that the unpaired electron in AgCl spends a little more time on ligands than on Ag²⁺, in agreement with the previous analysis of EPR and optical data for KCl:Ag²⁺. Furthermore, the five experimental optical transitions observed in that case are reasonably assigned. The first transition (observed at 12,500 cm^?1) is assigned to a jump involving the 5a_1g orbital built mainly (∽70%) from 3p orbitals of axial ligands, a fact that reflects the distinct level scheme for AgCl when compared to that for more ionic complexes. Calculations on AgF and AgF performed as a function of the equatorial Ag²⁺ –F^? distance led to a reasonable understanding of experimental gyromagnetic and superhyperfine tensors displayed by Ag²⁺ in fluorides. The different relative decrease undergone by g‖– go (8%) and g ? – go (28%) on passing from CsCdF₃:Ag²⁺ to RbCdF₃:Ag²⁺ is shown to be consistent with the formation of AgF and AgF complexes, respectively, related to the different substitutional position of Ag²⁺ in such lattices. The decrement of about 8.5% experienced by both g‖ – go and g? – go values on going from CsCdF₃:Ag²⁺ to NaF:Ag²⁺ is pointed out to reflect the different electrostatic potential (exerted by the rest of the lattice upon the complex) seen by AgF embedded in NaCl or perovskite-type lattices. © 1994 John Wiley & Sons, Inc.