Polarized basis sets of Slater-type orbitals: H to Ne atoms

We present three Slater-type atomic orbital (STO) valence basis (VB) sets for the first and second row atoms, referred to as the VB1, VB2, and VB3 bases. The smallest VB1 basis has the following structure: [3, 1] for the H and He atoms, [5, 1] for Li and Be, and [5, 3, 1] for the B to Ne series. For the VB2 and VB3 bases, both the number of shells and the number of functions per shell are successively increased by one with respect to VB1. With the exception of the H and Li atoms, the exponents for the VB1 bases were obtained by minimizing the sum of the Hartree-Fock (HF) and frozen-core singles and doubles configuration interaction (CISD FC) energies of the respective atoms in their ground state. For H and Li, we minimized the sum of the HF and CISD FC energies of the corresponding diatoms (i.e., of H(2) or Li(2)) plus the ground-state energy of the atom. In the case of the VB2 basis sets, the sum that was minimized also included the energies of the positive and negative ions, and for the VB3 bases, the energies of a few lowest lying excited states of the atom. To account for the core correlations, the VBx (x = 1, 2, and 3) basis sets for the Li to Ne series were enlarged by one function per shell. The exponents of these extended (core-valence, CV) basis sets, referred to, respectively, as the CVBx (x = 1, 2, and 3) bases, were optimized by relying on the same criteria as in the case of the VBx (x = 1, 2, and 3) bases, except that the full CISD rather than CISD FC energies were employed. We show that these polarized STO basis sets provide good HF and CI energies for the ground and excited states of the atoms considered, as well as for the corresponding ions.     

Synthesis, ionisation potentials and electron affinities of hexaazatrinaphthylene derivatives

Several hexaazatrinaphthylene derivatives and a tris(thieno)hexaazatriphenylene derivative have been synthesised by reaction of the appropriate diamines with hexaketocyclohexane. The crystal structure of 2,3,8,9,14,15-hexachloro-5,6,11,12,17,18-hexaazatrinaphthylene has been determined by X-ray diffraction; this reveals a molecular structure in good agreement with that predicted by density functional theory (DFT) calculations and pi-stacking with an average spacing between adjacent molecular planes of 3.18 A. Solid-state ionisation potentials have been measured by using UV photoelectron spectroscopy and fall in the range of 5.99 to 7.76 eV, whereas solid-state electron affinities, measured using inverse photoelectron spectroscopy, vary in the range -2.65 to -4.59 eV. The most easily reduced example is a tris(thieno)hexaazatriphenylene substituted with bis(trifluoromethyl)phenyl groups; DFT calculations suggest that the highly exothermic electron affinity is due both to the replacement of the outermost phenylene rings of hexaazatrinaphthylene with thieno groups and to the presence of electron-withdrawing bis(trifluoromethyl)phenyl groups. The rather exothermic electron affinities, the potential for adopting pi-stacked structures and the low intramolecular reorganisation energies obtained by DFT calculations suggest that some of these molecules may be useful electron-transport materials.     

Transition metals as electron traps. II. Structures,energetics and electron transfer dissociations of ternary Co,Ni and Zn–peptide complexes in the gas phase

Transition metal cations Co²⁺, Ni²⁺ and Zn²⁺ form 1 : 1 : 1 ternary complexes with 2,2′‐bipyridine (bpy) and peptides in aqueous methanol solutions that have been studied for tripeptides GGG and GGL. Electrospray ionization of these solutions produced singly charged [Metal(bpy)(peptide ? H)]⁺ and doubly charged [Metal(bpy)(peptide)]²⁺ ions (Metal = metal ion) that underwent charge reduction by glancing collisions with Cs atoms at 50 and 100 keV collision energies. Electron transfer to [Metal(bpy)(peptide)]²⁺ ions was less than 4.2 eV exoergic and formed abundant fractions of non‐dissociated charge‐reduced intermediates. Charge‐reduced [Metal(bpy)(peptide)]⁺ ions dissociated by the loss of a hydrogen atom, ammonia, water and ligands that depended on the metal ion. The Ni and Co complexes mainly dissociated by the elimination of ammonia, water, and the peptide ligand. The Zn complex dissociated by the elimination of ammonia and bpy. A sequence‐specific fragment was observed only for the Co complex. Electron transfer to [Metal(bpy)(peptide ? H)]⁺ was 0.6–1.6 eV exoergic and formed intermediate radicals that were detected as stable anions after a second electron transfer from Cs. [Metal(bpy)(peptide ? H)] neutrals and their anions dissociated by the loss of bpy and peptide ligands with branching ratios that depended on the metal ion. Optimized structures for several spin states, electron transfer and dissociation energies were addressed by combined density functional theory and Møller–Plesset perturbational calculations to aid interpretation of experimental data. The experimentally observed ligand loss and backbone cleavage in charge‐reduced [Metal(bpy)(peptide)]⁺ complexes correlated with the dissociation energies at the present level of theory. The ligand loss in ⁺CR^? spectra showed overlap of dissociations in charge‐reduced [Metal(bpy)(peptide ? H)] complexes and their anionic counterparts which complicated spectra interpretation and correlation with calculated dissociation energies. Copyright © 2009 John Wiley & Sons, Ltd.     

Density functional study of the structural,electronic, and magnetic properties of neutral and charged rhodium clusters up to tetramer

Model core potential computations were performed for Rh₂, Rh₃, and Rh₄ clusters and their respective cations and anions using the linear combination of Gaussian‐type orbital, nonlocal spin density method. The optimized geometries, electronic and magnetic structures, binding and fragmentation energies, adiabatic ionization potentials, and electron affinities were determined. Results show that the ionization potentials, electron affinities, binding energies, and magnetic moments decrease with the cluster size. For Rh₂ and Rh₃ the most stable structures exhibit ferromagnetic properties, while Rh₄ in its ground state is found to be paramagnetic. The structures of minimum energy for the charged species often differs from the corresponding neutral one. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Theoretical determination of electron affinity and ionization potential of DNA and RNA bases

The ionization potentials and electron affinities of thymine, cytosine, adenine, guanine, and uracil were determined at density functional level using different exchange‐correlation functionals and basis sets. Results showed that the computed ionization potentials are very close to the experimental counterparts. The sign of adiabatic electron affinities of adenine, thymine, and uracil is unaffected by the used level of theory while that for guanine and cytosine depends on both the used potential and basis set. Vertical electron affinities are always negative in agreement with the experimental indications. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 1243–1250, 2000     

The arsenic clusters Asn (n = 1-5) and their anions: structures, thermochemistry, and electron affinities

The molecular structures, electron affinities, and dissociation energies of the As(n)/As(-) (n) (n = 1-5) species have been examined using six density functional theory (DFT) methods. The basis set used in this work is of double-zeta plus polarization quality with additional diffuse s- and p-type functions, denoted DZP++. These methods have been carefully calibrated (Chem Rev 2002, 102, 231) for the prediction of electron affinities. The geometries are fully optimized with each DFT method independently. Three different types of the neutral-anion energy separations reported in this work are the adiabatic electron affinity (EA(ad)), the vertical electron affinity (EA(vert)), and the vertical detachment energy (VDE). The first dissociation energies D(e)(As(n-1)-As) for the neutral As(n) species, as well as those D(e)(As(-) (n-1)-As) and D(e) (As(n-1)-As(-)) for the anionic As(-) (n) species, have also been reported. The most reliable adiabatic electron affinities, obtained at the DZP++ BLYP level of theory, are 0.90 (As), 0.74 (As(2)), 1.30 (As(3)), 0.49 (As(4)), and 3.03 eV (As(5)), respectively. These EA(ad) values for As, As(2), and As(4) are in good agreement with experiment (average absolute error 0.09 eV), but that for As(3) is a bit smaller than the experimental value (1.45 +/- 0.03 eV). The first dissociation energies for the neutral arsenic clusters predicted by the B3LYP method are 3.93 eV (As(2)), 2.04 eV (As(3)), 3.88 eV (As(4)), and 1.49 eV (As(5)). Compared with the available experimental dissociation energies for the neutral clusters, the theoretical predictions are excellent. Two dissociation limits are possible for the arsenic cluster anions. The atomic arsenic results are 3.91 eV (As(-) (2) --> As(-) + As), 2.46 eV (As(-) (3) --> As(-) (2) + As), 3.14 eV (As(-) (4) --> As(-) (3) + As), and 4.01 eV (As(-) (5) --> As(-) (4) + As). For dissociation to neutral arsenic clusters, the predicted dissociation energies are 2.43 eV (As(-) (3) --> As(2) + As(-)), 3.53 eV (As(-) (4) --> As(3) + As(-)), and 3.67 eV (As(-) (5) --> As(4) + As(-)). For the vibrational frequencies of the As(n) series, the BP86 and B3LYP methods produce good results compared with the limited experiments, so the other predictions with these methods should be reliable.     

Theoretical investigation of the structures and properties of fluoromethyl peroxyl radicals

The equilibrium geometries, harmonic vibrational frequencies, charge distributions, spin density distributions, dipole moments, electron affinities (EAs), and C? O bond dissociation energies (BDEs) of HO, CH₃O, CH₂FO, CHF₂O, and CF₃O peroxyl radicals have been calculated using ab initio molecular orbital theory and density functional theory (DFT) at the B3LYP level. The C? H bond dissociation energies of the parent fluoromethanes have been calculated using the same levels of theory. Both the MP2(full) and B3LYP methods, using the 6‐31G(d,p) basis set, are found to be capable of accurately predicting the geometries of peroxyl radicals. Electron correlation accounts for ～25% of the C? H BDE of fluoromethanes and for ～50% of the C? O BDE of the corresponding peroxyl radicals. The B3LYP/6‐31G(d,p) method is found to be comparable to high ab initio levels in predicting C? O BDEs of studied peroxyl radicals and C? H BDEs of the parent alkanes. The progressive fluorine substitution of hydrogen atoms in methyl peroxyl radicals results in shortening of the C? O bond, lengthening of the O? O bond, an increase (decrease) of the spin density on the terminal (inner) oxygen, a decrease in the dipole moments, and an increase in electron affinities. Both C? O BDEs and EAs of peroxyl radicals (RO) correlate well with Taft σ* substituent constants for the R group in peroxyl radicals. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

Geometries,vibrational frequencies,and electron affinities of X2Cl (X=C,Si,Ge) clusters

Ab initio quantum chemical calculations have been performed on X₂Cl^? and X₂Cl (X = C, Si, Ge) clusters. The geometrical structures, vibrational frequencies, electronic properties and dissociation energies are investigated at the Hartree–Fock (HF), Møller–Plesset second‐ and fourth‐order (MP2, MP4), CCSD(T) level with the 6‐311+G(d) basis set. The X₂Cl (X = C, Si, Ge) and X₂Cl^? (X = Si, Ge) take a bent shape obtained at the ground state, while C₂Cl^? has a linear structure. The impact on internal electron transfer between the X₂Cl and the corresponding anional clusters is studied. The three different types of electron affinities (EAs) at the CCSD(T) are reported. The most reliable adiabatic electronic affinities, obtained at the CCSD(T)/cc‐pvqz level of theory, are predicted to be 3.30, 2.62, and 1.98 eV for C₂Cl, Si₂Cl, and Ge₂Cl, respectively. The calculated EAs of C₂Cl and Ge₂Cl are in good agreement with theoretical results reported. The correlation effects and basis sets effects on the geometrical structures and dissociation energies are discussed. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Core ionization energies of atoms and ions calculated using the generalized Sturmian method

The physical event of the umbrella inversion of ammonia has been studied in detail by application of the formalisms of frontier orbital theory, the density functional theory, the localized molecular orbital method, and the energy partitioning analysis. An intuitive structure for the transition state and dynamics of the physical process of structural reorganization prior to inversion have been suggested. The computation starts with the CNDO/2 equilibrium geometry, and thereafter the cycle proceeds through all the conformations of ammonia obtained by varying the ∠HNH angle in steps of 2° from its equilibrium value up to the transition state. The geometry of each conformation is optimized with respect to the length of the N–H bond. The glimpses of the charge density reorganization during the movement of the molecule from equilibrium conformation toward the transition state is computed in terms of dipole moment and the quantum mechanical hybridizations of bond pair and lone pair of N atom through the localized molecular orbitals (LMOs) of all the conformations. Results demonstrate that as the geometry of the molecule begins to evolve through the reorganization of structure, the N–H bond length and the dipole moment begin to decrease, and the trend continues up to the transition state. The dipole moment of the molecule at the suggested transition state is zero. The computed nature of quantum mechanical hybridization of bond pair and lone pair of the N atom as a function of reaction coordinates of the ∠HNH angles reveals that the percentage of s character of the lone pair hybrid decreases and that of the bond pair hybrid forming the σ(N–H) bond increases during the process of geometry reorganization from the equilibrium shape to the transition state. The rationale of the zero dipole moment of the transition state for inversion is not straightforward from its point‐group symmetry, but is self‐evident from its electronic structure drawn in terms of LMOs. The electronic structure of the transition state, which may be drawn in terms of the LMOs, seems to closely reproduce its suggested intuitive structure. The pattern of variation of dipole moment and nature of the changes of the percentage of the s character in the lone pair hybrid creating dipole with the evolution of geometry during the physical process of structural reorganization for the inversion are found to be nicely correlated according to the suggestion of Coulson. The profiles of the increasing strength of the N–H bond and the increasing percentage of s character of the bond pair hybrid of N atom forming this bond as a function of reaction coordinates are also found to be correlated in accordance with the suggestion of Coulson. The profile of global hardness as a function of reaction coordinate seems to demonstrate that the dynamics of the evolution of the molecular structure from equilibrium shape to the transition state following the course of suggested mode of structural reorganization conforms to the principle of maximum hardness (PMH). The profiles of parameters like the energies of highest occupied and lowest unoccupied molecular orbital (HOMO and LUMO), the gap in energy between HOMO and LUMO, the global hardness, the global softness, and chemical potential as a function of reaction coodinates of a continuous structural evolution from equilibrium shape to the transition state mimic the potential energy diagram of the total energy. Both the frontier orbital parameters and the density functional quantities are found to be equally effective and reliable to monitor the process of necessary structural reorganization prior to the inversion of mofecules. An energy partitioning analysis demonstrates that the origin of barrier has no unique single source rather as many as four mutually exclusive but interacting one‐ and two‐center energy terms within the molecule entail the origin and the height of the barrier. From a close analysis of the results, it seems highly probable that the necessary structural reorganization prior to umbrella inversion of ammonia most realistically occurs following the course of normal modes of vibration of the molecule. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 80: 1–26, 2000     

Tautomeric forms of adenine: Vertical ionization energies and Dyson orbitals

For the MP2/6‐311++g(2df,p) optimized geometry of all the 14 adenine tautomers, the first three vertical ionization energies have been calculated using several electron propagator decouplings. The corresponding Dyson orbitals provide detailed insight into the role of structural variations in different adenine tautomers. Changes in the electron binding energies and the corresponding Dyson orbital amplitudes have been correlated with tautomeric proton shifts and changes in conjugation patterns. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Absolute photoionization cross sections of furanic fuels: 2‐ethylfuran, 2‐acetylfuran and furfural

Absolute photoionization cross sections of the molecules 2‐ethylfuran, 2‐acetylfuran and furfural, including partial ionization cross sections for the dissociative ionized fragments, are measured for the first time. These measurements are important because they allow fuel quantification via photoionization mass spectrometry and the development of quantitative kinetic modeling for the complex combustion of potential fuels. The experiments are carried out using synchrotron photoionization mass spectrometry with an orthogonal time‐of‐flight spectrometer used for mass analysis at the Advanced Light Source of Lawrence Berkeley National Laboratory. The CBS‐QB3 calculations of adiabatic ionization energies and appearance energies agree well with the experimental results. Several bond dissociation energies are also derived and presented. Copyright © 2015 John Wiley & Sons, Ltd.     

On the accuracy of correlation energy calculated by the correlation factor method: first- and second-row atoms

Summary The correlation energy of two- and four-isoelectronic series, a representative example for which the existing spin-density functionals fails, is calculated using the Colle and Salvetti method, considering mono- and multideterminantal wave functions. The results are in agreement with experimental data, and show the potentiality of this method when it is applied to wave functions including the most relevant configurational features. Also, results for the ionization energies and electron affinities of first- and second-row atoms are reported.     

Comparative study of metal-porphyrins, -porphyrazines,and -phthalocyanines

A theoretical comparative study of complexes of porphyrin (P), porphyrazine (Pz), and phthalocyanine (Pc) with metal (M) = Fe, Co, Ni, Cu, and Zn has been carried out using a DFT method. The calculations provide a clear elucidation of the ground states for the MP/Pz/Pc molecules and for a series of [MP/Pz/Pc](x-) and [MP/Pz/Pc](y+) ions (x = 1, 2, 3, 4; y = 1, 2). There are significant differences among MP, MPz, and MPc in the electronic structure and other calculated properties. For FeP/Pz and CoP/Pz, the first oxidation occurs at the central metal, while it is the macroring of FePc and CoPc that is the site of oxidation. The smaller coordination cavity results in a stronger ligand field in Pz than in P. However, the benzo annulation produces a surprisingly strong destabilizing effect on the metal-macrocycle bonding. The effects of Cl axial bonding upon the electronic structures of the iron(III) complexes of P, Pz, and Pc were examined, as was the bonding of pyridine (py) to NiP, NiPz, and NiPc. The porphinato core size plays a crucial role in controlling the spin state of Fe(III) in these complexes. FePc(Cl) is predicted to be a pure intermediate-spin system, whereas NiPz(py)(2) and NiPc(py)(2) are metastable in high-spin (S = 1) states. The NiPz/Pc-(py)(2) binding energy curve has only a shallow well that facilitates decomposition of the complex. The NiP-(py)(2) bond energy is small, but the relatively deep well in the binding energy curve ought to make this system stable to decomposition.     

The importance of fully quantitative non-local treatments of exchange energy in Density Functional Theory of atoms,molecules and clusters,characterised by strongly inhomogeneous electron liquids

After a brief summary of some recent work on non-local approximations to the exchange energy and potential in density functional theory, a table is presented for atoms from Be to Ar, with an even number of electrons, having large electron density gradients, which highlights the importance of exchange. Until non-local exchange energy is truly quantitative, refinements of correlation energy may be submerged through ‘noisy’ exchange.     

Correlations of the physicochemical parameters of azaindoles with their reactions

Using the AM1 semiempirical quantum method the enthalpies of formation, ionization energies, electron affinities, energy differences between highest occupied and lowest unoccupied orbitals, atomic charges, bond orders, and dipole moments have been calculated for 4-, 5-, 6-, and 7-azaindoles. A correlation has been built up between the calculated physicochemical parameters and the Hammett para-substituent and inductive constants. The 1H to 7H proton transfer in 7-azaindole has been quantitatively described. __________ Translated from Khimiya Geterotsiklicheskikh Soedinenii, No. 7, pp. 1062–1072, July, 2006.     

SiBr~4及其不同离子态的电子结构研究: 紫外光 电子能谱(PES)和理论研究

本文详细报道了SiBr~4的紫外光电子能谱(PES)及该化合物不同离子态(X^2T~2,A^2T~1,B^2E,C^2T~2等)的电子结构和性质。实验测得对应基态离子态的绝热电离能I~a(X^2T~2←X^1A~1)=10.532eV,X^2T~2离子态的振动频率为(450±30)cm^-^1。结合理论计算对紫外光电子能谱进行了指认和分析,结果表明X^2T~2,A~2T~1两个离子态存在明显的自旋-轨道耦合作用,自旋-轨道耦合导致的分裂分别为：0.27ev和0.53eV。此外,不同计算方法比较显示外层格林函数方法计算得到电离能与实验吻合很好。     

The structure and energetics of triplet [B, C, F, H2]

In the current paper, we discuss our high level quantum chemical results for the structure and energetics of triplet (and hence open-shell) isomers corresponding to the stoichiometry of one boron, carbon, and fluorine apiece, and two hydrogens. While partially bond-ruptured excited ketene- and diazomethane-like H₂C-B-F and the carbene H(F)B-C-H plausibly emerge as the most stable isomers, a variety of novel structural features emerge for the assembled energy minima of at least 16 species. All of these species are compared as well as transition states that connect them. Comparison is made with corresponding forms of the singlet species with this stoichiometry, shown earlier by us to have a rich diversity of structures as well as a large range of energies and relative stabilities.     

Spontaneous ionization of N-alkylphenothiazine molecules adsorbed in channel-type zeolites: effects of alkyl chain length and confinement on electron transfer

The mere mixing of N-alkylphenothiazines with three channel-type acid zeolites with various structures (ferrierite, H-MFI, and mordenite) induces the spontaneous ionization of the heterocyclic molecule in high yield upon adsorption. The diffuse reflectance UV-visible absorption and Raman scattering spectra show that the accessibility of the highly polarizing acid sites is not indispensable to induce the spontaneous ionization process. Due to their particularly low ionization potential values (6.7 eV), the adsorption of the molecules on the external surface or in the inner volume is the key parameter to generate the radical cation. However, the ionization yield and charge stabilization are intimately correlated to the possibility of the zeolites accommodating molecules inside their channels. Moreover, the higher electrostatic field gradient induced by high confinement is required to favor the second ionization and dication formation. The alkyl chain length plays a decisive role by either slowing down the diffusion process or blocking the molecule at the pore entry. Therefore, the efficiency of the ionization process that depends on the number of adsorbed molecules decreases significantly from phenothiazine to the N-alkylphenothiazines. The spectral data demonstrate that deformation of the alkyl group is necessary to allow the diffusion of the molecules into the channels.