Part I. Approximate Hartree-Fock wavefunctions one-electron properties and electronic structure of the water molecule. Part II. Perturbation-variational calculation of the nuclear spin-spin isotope coupling constant in HD

Part I

Several approximate Hartree-Fock SCF wavefunctions for the ground electronic state of the water molecule have been obtained using an increasing number of multicenter s, p, and d Slater-type atomic orbitals as basis sets. The predicted charge distribution has been extensively tested at each stage by calculating the electric dipole moment, molecular quadrupole moment, diamagnetic shielding, Hellmann-Feynman forces, and electric field gradients at both the hydrogen and the oxygen nuclei. It was found that a carefully optimized minimal basis set suffices to describe the electronic charge distribution adequately except in the vicinity of the oxygen nucleus. Our calculations indicate, for example, that the correct prediction of the field gradient at this nucleus requires a more flexible linear combination of p-orbitals centered on this nucleus than that in the minimal basis set. Theoretical values for the molecular octopole moment components are also reported.

Part II

The perturbation-variational theory of R. M. Pitzer for nuclear spin-spin coupling constants is applied to the HD molecule. The zero-order molecular orbital is described in terms of a single 1s Slater-type basis function centered on each nucleus. The first-order molecular orbital is expressed in terms of these two functions plus one singular basis function each of the types e^-r/r and e^-r ln r centered on one of the nuclei. The new kinds of molecular integrals were evaluated to high accuracy using numerical and analytical means. The value of the HD spin-spin coupling constant calculated with this near-minimal set of basis functions is J_HD = +96.6 cps. This represents an improvement over the previous calculated value of +120 cps obtained without using the logarithmic basis function but is still considerably off in magnitude compared with the experimental measurement of J_HD = +43 0 ± 0.5 cps.

Triel Bonds, pi-Hole-pi-Electrons Interactions in Complexes of Boron and Aluminium Trihalides and Trihydrides with Acetylene and Ethylene

MP2/aug-cc-pVTZ calculations were performed on complexes of aluminium and boron trihydrides and trihalides with acetylene and ethylene. These complexes are linked through triel bonds where the triel center (B or Al) is characterized by the Lewis acid properties through its -hole region while -electrons of C2H2 or C2H4 molecule play the role of the Lewis base. Some of these interactions possess characteristics of covalent bonds, i.e., the Al--electrons links as well as the interaction in the BH3-C2H2 complex. The triel--electrons interactions are classified sometimes as the 3c-2e bonds. In the case of boron trihydrides, these interactions are often the preliminary stages of the hydroboration reaction. The Quantum Theory of Atoms in Molecules as well as the Natural Bond Orbitals approach are applied here to characterize the -hole--electrons interactions.

Computational study on the addition of HCN to methanimine catalyzed by formamidine and formamide

In this paper, we have carried out a theoretical study on the addition of HCN to methanimine with formamidine or formamide using second-order Moller-Plesset perturbation (MP2) method with 6-31 + G(d,p) basis sets. At MP2 level. a high-energy, intermediate has been located for each pathway. The addition of HCN to methanimine with formamidine has the lowest free energy barrier according to the calculations at MP2 level. (C) 2004 Elsevier B.V. All rights reserved.

A DFT study on isomorphously substituted MCM-22 zeolite

The density functional theory has been used to study the isomorphously substituted MCM-22 zeolite for the first time. The effect of the basis sets on the calculation results is discussed in details. Data of several index properties for characterizing the relative acidity of T-MCM-22 (T = B, Al, Ga, and Fe), including proton affinity, bond length and bond angle, OH stretching frequency, and charge on the acidic proton, show that the acidity of T-MCM-22 increases in the sequence of B-MCM-22 < Fe-MCM-22 < Ga-MCM-22 < Al-MCM-22. After making a correction, the calculated OH stretching frequencies for Al-MCM-22 and Fe-MCM-22 show a reasonable agreement with the experimental data. On the basis of an equilibrium structure of the B-MCM-22 zeolite, the effect of the B element in the synthesis of the Ti-MCM-22 is also discussed. The adding of the B element during the synthesis of the Ti-MCM-22 can decrease greatly the Ti substitution energy because of the forming of a structure quite similar to the terminal silanol group. The results can provide some constructively information for zeolite synthesis.

Theoretical study on analytical potential function and spectroscopic parameters for CaF molecule

The equilibrium properties and potential energy curves of the ground electronic state of CaF have been calculated using the Brueckner Doubles calculation with a triples contribution added [BD(T)] and the gradient-corrected density functional theory with three-parameter exact exchange mixing (B3LY-P) method, with 6-311 + G*,6-311 + G(2df,2pd) and 6-311 + G(3df,3pd) basis sets. All the computational PECs are fitted to analytical potential energy functions using Murrell-Sorbie, Huxley and Tang-Toennies potentials. Based on this, the spectroscopic parameters are calculated, and then compared with some other theoretical and experimental data. (C) 2004 Elsevier B.V. All rights reserved.

Density functional studies of diatomic LaO to LuO

Bond distances, vibrational frequencies, electron affinities, ionization potentials, dissociation energies and dipole moments of the title molecules in neutral, positively and negatively charged ions were studied by use of density functional method. Ground electronic state was assigned for each molecule. The bonding patterns were analyzed and compared with both the available data and across the series. It was found that besides ionic component, covalent bonds are formed between the metal s, d and f orbitals and oxygen p orbitals. Contrary to the well known lanthanide contraction, the bond distance is not regular from LaO to LuO for both neutral and charged molecules. An obvious population at 5d orbital was observed through the lanthanide series. 4f electrons also participate the chemical bonding for CeO to NdO and TbO to TmO. For EuO, GdO, YbO and LuO, 4f electrons tend to be localized. The spin multiplicity is regular for neutral and charged molecules. The spin multiplicity of the charged molecules can be obtained by -1 (or +1 for TbO+, DyO+, YbO- and YbO+) compared with the corresponding neutral molecules.

Structural stability and electronic state of transition metal trimers

Ground state geometries were searched for transition metal trimers Sc-3, Y-3, La-3, Lu-3, Ti-3, Zr-3, and Hf-3 by density functional methods. For all the studied trimers, our calculation indicates that the ground state geometries are either equilateral triangle (Zr-3 and Hf-3) or near equilateral triangle (Ti-3, Sc-3, Y-3, La-3, and Lu-3). For rare earth trimers Sc-3, Y-3, La-3, and Lu-3, isosceles triangle (near equilateral triangle) at quartet state is the ground state. Isosceles triangle at doublet state is the competitive candidate for the ground state. For Zr-3 and Hf-3, equilateral triangle at singlet state is the most stable. For Ti-3, isosceles triangle (near equilateral triangle) at quintet state gives the ground state. For Sc-3, Zr-3, and Hf-3, where experimental results are available, the predicted geometries are in agreement with experiment in which the ground state is equilateral triangle (Zr-3) or fluxional (Sc-3 and Hf-3). For Y-3, the calculated geometry is in agreement with experimental observation and previous theoretical study that Y-3 is a bent molecule for the ground state.

Ground state of lutetium dimer by density functional methods

Bond distances, vibrational frequencies and dissociation energies for the ground state of Lu-2 were studied by density functional methods B3LYP, B3PW91, BLYP, BHLYP, BP86, B3P86, MPW1PW91, PBE1PBE and SVWN with CEP-121G and SDD basis sets. Singlet state is predicted to be the most stable. CEP-121G has a better overall performance than SDD. At CEP-121G basis set, all density functional methods used in this study perform well in reproducing the spectroscopic constants.

Electron affinities and ionization potentials of 4d and 5d transition metal atoms by CCSD(T), MP2 and density functional theory

The electron affinities and ionization potentials of 4d and 5d transition metal atoms were studied by CCSD(T), MP2 and density functional methods. The calculated results indicate that density functional method B3LYP has the best overall performance in predicting both electron affinity and ionization potential. SVWN gives largest IP and EA for 4d and 5d atoms. For the two basis sets used in this study, LANL2DZ and SDD, the performance of B3LYP/SDD combination is better than B3LYP/LANL2DZ, in particular for electron affinity calculation.

Investigation of numerical atomic orbitals for first-principles calculations of the electronic and transport properties of silicon nanowire structures

This thesis is focused on the application of numerical atomic basis sets in studies of the structural, electronic and transport properties of silicon nanowire structures from first-principles within the framework of Density Functional Theory. First we critically examine the applied methodology and then offer predictions regarding the transport properties and realisation of silicon nanowire devices. The performance of numerical atomic orbitals is benchmarked against calculations performed with plane waves basis sets. After establishing the convergence of total energy and electronic structure calculations with increasing basis size we have shown that their quality greatly improves with the optimisation of the contraction for a fixed basis size. The double zeta polarised basis offers a reasonable approximation to study structural and electronic properties and transferability exists between various nanowire structures. This is most important to reduce the computational cost. The impact of basis sets on transport properties in silicon nanowires with oxygen and dopant impurities have also been studied. It is found that whilst transmission features quantitatively converge with increasing contraction there is a weaker dependence on basis set for the mean free path; the double zeta polarised basis offers a good compromise whereas the single zeta basis set yields qualitatively reasonable results. Studying the transport properties of nanowire-based transistor setups with p+-n-p+ and p+-i-p+ doping profiles it is shown that charge self-consistency affects the I-V characteristics more significantly than the basis set choice. It is predicted that such ultrascaled (3 nm length) transistors would show degraded performance due to relatively high source-drain tunnelling currents. Finally, it is shown the hole mobility of Si nanowires nominally doped with boron decreases monotonically with decreasing width at fixed doping density and increasing dopant concentration. Significant mobility variations are identified which can explain experimental observations.

Accurate localized resolution of identity approach for linear-scaling hybrid density functionals and for many-body perturbation theory

© 2015 IOP Publishing Ltd and Deutsche Physikalische Gesellschaft.A key component in calculations of exchange and correlation energies is the Coulomb operator, which requires the evaluation of two-electron integrals. For localized basis sets, these four-center integrals are most efficiently evaluated with the resolution of identity (RI) technique, which expands basis-function products in an auxiliary basis. In this work we show the practical applicability of a localized RI-variant ('RI-LVL'), which expands products of basis functions only in the subset of those auxiliary basis functions which are located at the same atoms as the basis functions. We demonstrate the accuracy of RI-LVL for Hartree-Fock calculations, for the PBE0 hybrid density functional, as well as for RPA and MP2 perturbation theory. Molecular test sets used include the S22 set of weakly interacting molecules, the G3 test set, as well as the G2-1 and BH76 test sets, and heavy elements including titanium dioxide, copper and gold clusters. Our RI-LVL implementation paves the way for linear-scaling RI-based hybrid functional calculations for large systems and for all-electron many-body perturbation theory with significantly reduced computational and memory cost.

12C2H2-Ar van der Waals complex

New theoretical and experimental results on the acetylene-Ar van der Waals complex are presented and the literature is reviewed. New ab initio calculations at the MP2 level were performed using large basis sets with diffuse functions and taking into account the basis set superposition error. It was found that the structure of acetylene is not significantly altered by the complexation and that its vibrational frequencies are only slightly lowered. Finally, it was observed that the calculated properties of the complex (structure, vibrational spectrum, bond dissociation energy) are not sensitive to the structure imposed on acetylene. Experimentally, acetylene-Ar was produced in a supersonic expansion under experimental conditions corresponding to 9 K rotational temperature. Thanks to the performances of CW-CRDS detection, the Ka = 0 ← 1, 1 ← 0, and 2 ← 1 sub- bands of the v1 + v3 band could be recorded and resolved and most of their lines assigned. Upper-state rotational constants were fitted, however not including the upper Ka = 2 state, which shows K-doubling the opposite of the expected. The Lorentzian width of most line profiles sets the mean lifetime to some 7.5 ns. Local perturbations affecting line positions and/or line widths are demonstrated. Additional series of lines tentatively attributed to acetylene-Ar are discussed.© 2009 American Chemical Society.

Accurate ground-state potential energy surfaces of the C2H2–Kr and C2H2–Xe van der Waals complexes

Accurate ab initio intermolecular potential energy surfaces (IPES) have been obtained for the first time for the ground electronic state of the C 2H2-Kr and C2H2-Xe van der Waals complexes. Extensive tests, including complete basis set and all-electron scalar relativistic results, support their calculation at the CCSD(T) level of theory, using small-core relativistic pseudopotentials for the rare-gas atoms and aug-cc-pVQZ basis sets extended with a set of 3s3p2d1f1g mid-bond functions. All results are corrected for the basis set superposition error. The importance of the scalar relativistic and rare-gas outer-core (n.1)d correlation effects is investigated. The calculated IPES, adjusted to analytical functions, are characterized by global minima corresponding to skew T-shaped geometries, in which the Jacobi vector positioning the rare-gas atom with respect to the center of mass of the C2H2 moiety corresponds to distances of 4.064 and 4.229Å, and angles of 65.22° and 68.67° for C 2H2-Kr and C2H2-Xe, respectively. The interaction energy of both complexes is estimated to be -151.88 (1.817 kJ mol-1) and -182.76 cm-1 (2.186 kJ mol-1), respectively. The evolution of the topology of the IPES as a function of the rare-gas atom, from He to Xe, is also discussed. © 2012 Taylor and Francis.

Equilibrium structure of the hydrogen bonded dimer H2O ⋯ HF

The equilibrium structure of the hydrogen bonded complex H2O HF has been calculated ab initio using the CCSD(T) method with basis sets up to sextuple- quality with diffuse functions and taking into account the basis set superposition error correction. The calculations carried out confirm the importance of diffuse functions and of counterpoise correction to obtain an accurate geometry. The most important point is that the basis set convergence is extremely slow and, for this reason an accurate ab initio structure requires a very large basis set. Nevertheless, the ab initio structure is significantly different from the experimental r0 and rm structures. Analysis of the basis set convergence and of the approximations used for the determination of the experimental structures indicates that the ab initio structure is expected to be more reliable.

Many-body theory of positron-atom interactions

A many-body theory approach is developed for the problem of positron-atom scattering and annihilation. Strong electron- positron correlations are included nonperturbatively through the calculation of the electron-positron vertex function. It corresponds to the sum of an infinite series of ladder diagrams, and describes the physical effect of virtual positronium formation. The vertex function is used to calculate the positron-atom correlation potential and nonlocal corrections to the electron-positron annihilation vertex. Numerically, we make use of B-spline basis sets, which ensures rapid convergence of the sums over intermediate states. We have also devised an extrapolation procedure that allows one to achieve convergence with respect to the number of intermediate- state orbital angular momenta included in the calculations. As a test, the present formalism is applied to positron scattering and annihilation on hydrogen, where it is exact. Our results agree with those of accurate variational calculations. We also examine in detail the properties of the large correlation corrections to the annihilation vertex.