Ab initio molecular orbital calculations

A series of non-empirical calculations on furan, pyrrole and 1,2,5-oxadiazole are reported in which the effect of polarisation functions added to the minimal 7s 3p basis on each atom is studied. The effect on these planar molecules is largely through the rather than the-system. A comparison with the results of work with scaled functions is reported. Both series are shown to lead to much improved agreement with the electron spectroscopy energy levels. The effect on the dipole moments of these changes in basis is more variable but, with the exception of furan, the agreement with experiment is improved in the present method.

---

Zusammenfassung Für die Moleküle Furan, Pyrrol und 1,2,5-Oxadiazol wurde eine Reihe von nichtempirischen Rechnungen durchgeführt, in denen der Einfluß von zusätzlichen Polarisationsfunktionen zur minimalen 7s 3p-Basis an jedem Atom untersucht wird. Die Ergebnisse werden mehr durch die Art der Beschreibung des Systems der-Elektronen als durch diejenige der-Elektronen beeinflußt. Ein Vergleich mit den Ergebnissen bei Verwendung skalierter Funktionen wird durchgeführt. Beide Reihen von Ergebnissen zeigen eine verbesserte Übereinstimmung zu den Energiemeßwerten der Elektronenspektroskopie. Die Änderungen des berechneten Dipolmoments bei derartigen Basisvariationen sind größer als bei früheren Methoden. Die Übereinstimmung mit dem Experiment wird, mit Ausnahme von Furan, jedoch verbessert.

     

Structures and relative stability of medium-sized silicon clusters. IV. Motif-based low-lying clusters Si21-Si30

Structures and relative stability of four families of low-lying silicon clusters in the size range of Sin(n=21-30) are studied, wherein two families of the clusters show prolate structures while the third one shows near-spherical structures. The prolate clusters in the first family can be assembled by connecting two small-sized magic clusters Sin (n=6, 7, 9, or 10) via a fused-puckered-hexagonal-ring Si9 unit (a fragment of bulk diamond silicon), while those in the second family can be constructed on the basis of a structural motif consisting of a puckered-hexagonal-ring Si6 unit (also a fragment of bulk diamond silicon) and a small-sized magic cluster Sin (n=6, 7, 9, or 10). For Si21-Si29, the predicted lowest-energy clusters (except Si27) exhibit prolate structures. For clusters larger than Si25, the third family of near-spherical clusters becomes energetically competitive. These near-spherical clusters all exhibit endohedral cagedlike structures, and the cages are mostly homologue to the carbon-fullerene cages which consist of pentagons and hexagons exclusively. In addition, for Si26-Si30, we construct a new (fourth) family of low-lying clusters which have "Y-shaped" three-arm structures, where each arm is a small-sized magic cluster (Si6, Si7, or Si10). Density-functional calculation with the B3LYP functional shows that this new family of clusters is also energetically competitive, compared to the two prolate and one near-spherical low-lying families.     

Structures and relative stability of medium-sized silicon clusters. V. Low-lying endohedral fullerenelike clusters Si31-Si40 and Si45

We have performed unconstrained search for low-lying structures of medium-sized silicon clusters Si(31)-Si(40) and Si(45), by means of the minimum-hopping global optimization method coupled with a density-functional based tight-binding model of silicon. Subsequent geometric optimization by using density-functional theory with the PBE, BLYP, and B3LYP functionals was carried out to determine the relative stability of various candidate low-lying silicon clusters obtained from the unconstrained search. The low-lying characteristics of these clusters can be affirmed by comparing the binding energies per atom of these clusters with previously determined lowest-energy clusters(Si(n)) in the size range of 21相似文献   

Ab initio molecular orbital calculations for butadiene and its ions

Ab initio calculations are reported for cis and trans butadiene and some of their ions. The calculations are compared with semi-empirical results, and used to predict coupling constants.     

Ab initio molecular orbital calculations on the pyrazine-lithium ion pair

The results of ab initio molecular orbital calculations are presented for the pyrazine^?/Li⁺ ion pair. It is shown that the lowest energy conformation is where the lithium occupies a position in the plane of the pyrazine ring, along the C ₂ axis passing through the nitrogen atoms.     

Ab initio calculations of structure and stability of small boron nitride clusters

Clusters of boron nitride B_xN_x (x = 1–4, 12, 15, 30) were investigated by the Hartree-Fock and density functional methods using the 6-31G* basis. It was found that linear, cyclic, and shell structures are stable against minor deformations of the B_xN_x cluster. Inclusion of electron correlation in calculation markedly changes the electron density distribution and the structure of the clusters.     

Ab initio molecular orbital theory study of GaAs clusters: The geometry

We have studied the structure and geometry of neutral and charged atomic clusters consisting of Ga and As atoms via ab initio Hartree–Fock (HF) and second‐order Møller–Plesset methods. The Ga_mAs_n cluster with m≠n composition prefers a nontetrahedral geometry in the charge neutral (q=0) state. These clusters tend to be stable in tetrahedral geometry when appropriately charged. The Ga_mAs_n cluster with m=n composition (1:1 ratio of Ga to As atoms) tends to be stable in a tetrahedral geometry in the charge neutral (q=0) state. With increasing size of the cluster, the geometry of Ga_nAs_n cluster approaches the zinc‐belende‐type crystalline structure. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 77: 563–573, 2000     

Ab initio molecular orbital calculations on beryllium and magnesium atom reactions with water

Ab initio molecular orbital calculations are reported for beryllium and magnesium atom complexes with water, Be. OH₂ and Mg OH₂, and for species resulting from insertion of the metal atom(s) into the water molecule, HBeOH, HMgOH, HBeOBeH, and HMgOMgH     

Ab initio molecular orbital calculations on the water-carbon dioxide system: Molecular complexes

An ab initio LCAO MO SCF calculation has been carried out on three different complexes between water and carbon dioxide. It has been found that a planar complex in which the carbon atom is bound to the oxygen atom in the water is the most stable, with a calculated energy of formation of ?25.7 kJ/mol.     

Ab initio molecular orbital calculations on some selected boron(I) hydrides

Ab initio LCAO-MO-SCF Hartree-Fock-Roothaan calculations have been carried out for a series of closed-shell boron hydrides (B³⁺, B⁺ BH, BH₂⁻,BH₂⁺, BH₃ and BH₄⁻). Koopmans' theorem vertical ionization potentials for the core and highest occupied molecular orbitals are presented. Proton and hydride affinities of some of these molecules have been calculated along with the energies of reaction between various pairs of these boron hydrides.     

Ab initio molecular orbital calculations on the associated complexes of lithium cyanide with ammonia

Ab initio molecular orbital (MO ) calculations with the 3-21G and 6-31G basis sets are carried out on a series of complexes of NH₃ with Li⁺, C?N^?, LiCN, and its isomer LiNC. The BSSE -corrected interaction energies, geometrical parameters, internal force constants, and harmonic vibrational frequencies are evaluated for 15 species. Complexes with trifurcated (C_3v) structures are calculated to be saddle points on the potential energy surfaces and have one imaginary frequency each. Calculated energies, geometrical parameters, internal force constants, and harmonic vibrational frequencies of the various species considered are discussed in terms of the nature of association of LiCN with ammonia. The vibrational frequencies of the relevant complexed species are compared with the experimental frequencies reported earlier for solutions of lithium cyanide in liquid ammonia. © 1995 John Wiley & Sons, Inc.     

Structures and stability of medium-sized silicon clusters. III. Reexamination of motif transition in growth pattern from Si15 to Si20

It has been established from experiments that stable medium-sized ionic clusters Si15-Si20 are prolate in shape. Density-functional theories (DFTs) also predict that nearly all low-lying neutral clusters in this size range are prolate in shape. Moreover, most of them are built onto two generic structural motifs, either the tricapped-trigonal-prism (TTP) Si9 motif or the six/six Si6Si6 (sixfold-puckered hexagonal ring Si6 plus six-atom tetragonal bipyramid Si6) motif. However, it appears that the exact location of the TTP-to-six/six motif transition is dependent on the functional (e.g., PBE or BLYP) used in the DFT calculations. Here, we present total-energy calculations for two series of clusters (one series containing six/six motif and the other containing the TTP motif) in the size range of Si16-Si20. The calculations were based on all-electron DFT methods with a medium [6-311G (2d)] and a large (cc-pVTZ) basis sets, as well as coupled-cluster single and double substitutions (including triple excitations) [CCSD(T)] method with a modest (cc-pVDZ) basis set. In the DFT calculations, two popular hybrid density functionals, the B3LYP and PBE1PBE, were selected. It is found that the B3LYP total-energy calculations slightly favor the six/six motif, whereas the PBE1PBE calculations slightly favor the TTP motif. The CCSD(T) total-energy calculations, however, show that isomers based on the six/six motif are energetically slightly favorable in the size range of Si16-Si20. Hence, the TTP-to-six/six motif transition is more likely to occur at Si16.     

Ab initio molecular dynamics simulation of liquid Al88Si12 alloys

First-principles molecular dynamics simulations are carried out to study the structures, dynamics, and electronic properties of liquid Al88Si12 in the temperature ranging from 898 to 1298 K. The temperature dependence of static structure factors, pair correlation functions, and electronic density-of-states are investigated. The structural properties obtained from the simulations are in good agreement with the x-ray diffraction experimental results.     

Ab initio molecular orbital studies of polyethylene deformation

A molecular LCAO Hartree-Fock procedure was used to calculate total energies of axially stretched normal paraffins containing up to nine carbon atoms. The results are used to model the mechanical properties of polyethylene.     

Ab initio molecular orbital study of 1-fluorosilatrane

The molecular geometry of 1-fluorosilatrane was optimized fully by restricted Hartree–Fock (HF) calculations using the 3-21G, 3-21G(d) and 6-31G(d) basis sets, with the aim of locating the positions of the local minima on the energy hypersurface. The optimized geometries were compared with available experimental (X-ray and ED) and semiempirical data. The ab initio calculations using polarized basis sets are in good agreement with those of previously reported semiempirical calculations, giving a slightly longer equilibrium Si? N distance (～ 256 pm) in the case of the endo minimum. However, the exo minimum predicted by the semiempirical methods is not supported. There was no experimental evidence for the existence of this exo minimum, and the present ab initio calculations suggest that it is highly unstable. There is considerable disagreement among the experimental results in the C? N and C? C bond lengths in various silatranes, their differences being as large as 13 pm. The present calculations predict that these differences may appear because the silatrane skeleton is flexible with low-energy, large-amplitude internal motions which introduce considerable uncertainties into the position of ring carbon atoms. © 1994 by John Wiley & Sons, Inc.     

Ab initio molecular orbital calculations for 3,6-dihydro-1,2-dithiin and 3,6-dihydro-1,2-dioxin

Optimized geometries and total energies for the conformers of 3,6-dihydro-1,2-dithiin ( 2 ) and 3,6-dihydro-1,2-dioxin ( 3 ) were calculated at several ab initio MO levels: RHF/3-21G(*), RHF/6-31G*, MP2/6-31G*, and MP2/6-31G*/ /RHF/3-21G(*). For the dioxin, in addition to the above levels the corresponding nonextended basis sets ab initio methods were also carried out. The dithiin results are compared with those of simple disulfanes, HSSH and (CH₃)₂S₂, whose optimized geometries agree closely with the observed structures, which is the gauche (C₂ symmetry). For the disulfanes, the gauche geometries from RHF/3-21G(*) are in good agreement with the observed structure while the RHF/3-21G results best fit the dioxin. Pertinent structural data at the RHF/3-21G(*) for the half-chair (C₂) dithiin are: bond lengths, ? SS? , ? CS? , ? CC?, and ? C?C? , 2.050, 1.817, 1.515, and 1.317 Å, respectively; bond angles, CSS, ?CCS, and C?CS, 98.0, 114.2, and 127.8°, respectively; CSSC dihedral angle of 63.2°; and twist angle of 36.5°. The total energy for half-chair dithiin at MP2/6-31G*//RHF/3-21G(*) is less than the planar (C_2v) and the half-boat (C_s) structures by 69.67 and 29.05 kJ/mol, respectively. The calculated structural data (vs. observed) at RHF/3-21G for the half-chair dioxin are: bond lengths, ? OO? , ? CO? , ? CC?, and C?C, 1.464 (1.463), 1.454, 1.509, and 1.313 Å (1.338 Å), respectively; bond angles, COO, ?CCO, and C?CO, 105.0, 109.8 (110.3), and 120.7° (119.9°), respectively; COOC dihedral angle of 79.7° (80 ± 2°); and twist angle of 39.0 (38.3°). The total energy for half-chair dioxin at MP2/6-31G//RHF/3-21G is less than the planar and the half-boat structures by 70.35 and 42.85 kJ/mol, respectively. The total energies calculated at the extended basis sets (*) ab initio levels for the C₂ symmetry dioxin are much lower than those of the nonextended basis sets. © John Wiley & Sons, Inc.     

Ab initio and indo molecular orbital calculations of the geometries and electronic structures of substituted sulfines

Molecular orbital calculations are presented for the parent sulfine and some mono-and dihalogensubstituted sulfines, using ab initio and INDO methods. A partial geometry optimization was performed for nine different sulfines. Charge distributions, potential surfaces and dipole moments were calculated from the wave functions of the optimized geometries. Cis-trans interconversion barriers and electronic spectra are also given. The atomic charges of the S and O atoms are insensitive to substitutions at carbon, and substituents greatly influence the potential in the environment of the molecule. The implications of the results for the chemical behavior of sulfine derivates are briefly discussed and compared with experimental data.     

Ab initio molecular orbital calculations of the first two adiabatic ionizations of the water dimer

Ab initio molecular orbital methods have been used to study the first two adiabatic ionizations of the water dimer. Both ionizations cause considerable geometrical rearrangements in the water dimer structure. Good agreement is found with experiment on the difference between the vertical and adiabatic energies for the first ionization.