A new form of electric-field-dependent basis functions for the calculation of electric polarizabilities and dipole moments

An ab-initio molecular orbital theory of electrical polarization is presented in which the molecular orbitals are written as linear combinations of atomic functions which depend explicitly on the strength of a uniform external electric field. The wavefunctions in the presence of such a field are determined using self-consistent field perturbation theory. It is shown that the use of field-dependent atomic functions provides an efficient technique for the calculation of electric polarizability tensors. Polarizability tensors and electric-dipole moments calculated using both a minimal and a split-valence-shell basis set are compared with experimental results. Both polarizability-tensor components and dipole moments are seriously underestimated at the minimal bases-set level. The split-valence basis approach yields substantially better results; the calculated values at this level are in reasonable agreement with the corresponding experimental values. The experimental ordering of isotropic polarizabilities for a set of small molecules is duplicated quite closely by both the minimal and the split-valence-shell calculations.     

Direct inversion of the iterative subspace with contracted planewave basis functions

Ways to reduce the computational cost of periodic electronic structure calculations by using basis functions corresponding to linear combinations of planewaves have been examined recently. These contracted planewave (CPW) basis functions correspond to Fourier series representations of atom‐centered basis functions, and thus provide access to some beneficial properties of planewave (PW) and localized basis functions. This study reports the development and assessment of a direct inversion of the iterative subspace (DIIS) method that employs unique properties of CPW basis functions to efficiently converge electronic wavefunctions. This method relies on access to a PW‐based representation of the electronic structure to provide a means of efficiently evaluating matrix–vector products involving the application of the Fock matrix to the occupied molecular orbitals. These matrix–vector products are transformed into a form permitting the use of direct diagonalization techniques and DIIS methods typically employed with atom‐centered basis sets. The abilities of this method are assessed through periodic Hartree–Fock calculations of a range of molecules and solid‐state systems. The results show that the method reported in this study is approximately five times faster than CPW‐based calculations in which the entire Fock matrix is calculated. This method is also found to be weakly dependent upon the size of the basis set, thus permitting the use of larger CPW basis sets to increase variational flexibility with a minor impact on computational performance. © 2018 Wiley Periodicals, Inc.     

Are atoms sic to molecular electronic wavefunctions? II. Analysis of fors orbitals

Electronic wavefunctions that describe molecules in the full optimized reaction space (FORS) are multiconfigurational wavefunctions which are invariant under non-singular linear transformations of the occupied molecular orbitals. They offer therefore a considerably wider scope for orbital interpretations than the single-configuration Hartree-Fock approximation. For example they can be analyzed in terms of natural MOs and in terms of localized MOs. The latter turn out to be remarkably atomic in character and a new localization procedure can be formulated which yields atom-adapted molecular orbitals. These have the character of minimal-basis-set AOs that are optimally adapted to the molecular environment and furnish an unambigious atomic population analysis. On the other hand, chemically adapted molecular orbitals can be defined by an appropriate compromise between natural orbitals and localized orbitals. The freedom to use, as configuration-generating molecular orbitals, atom-adapted FORS MOs as well as chemically adapted FORS MOs makes FORS wavefunctions particularly suitable for chemical interpretations. The ensuing analysis establishes the minimal basis set (in molecule-adapted form) as a theoretically sound concept for the understanding of accurate molecular wavefunctions. An illustrative example is discussed.     

Structure of X-ray photoelectron, emission, and conversion spectra and molecular orbitals of uranium compounds

The fine structure of X-ray photoelectron spectra of uranium compounds in the range of electron binding energies from 0 to ∼50 eV is largely determined by the electrons of the outer and inner valence molecular orbitals arising from the valence atomic shells, including the U6p and Lns low-energy occupied atomic shells. This result is in agreement with the data of the electronic structure calculations of these compounds and confirmed by the nuclear electron (conversion) and X-ray emission spectroscopic investigations. It is shown that the fine structure of X-ray photoelectron spectra associated with the electrons of inner valence molecular orbitals makes it possible to judge the participation of the electrons of the occupied atomic shells in chemical bonding, the structure of the nearest environments of the atom, and the bond lengths in the compounds. The overall contribution of the electrons of these molecular orbitals to the absolute value of binding energy may prove to be comparable to the contribution of the electrons of the outer valence molecular orbitals to atomic bonding. This is a new and important fact in chemistry. Translated fromZhurnal Strukturnoi Khimii, Vol. 39, No. 6, pp. 1037–1046, November–December, 1998.     

Full configuration interaction calculation of the low lying valence and Rydberg states of BeH

The all-electron full configuration interaction (FCI) vertical excitation energies for some low lying valence and Rydberg excited states of BeH are presented in this article. A basis set of valence atomic natural orbitals has been augmented with a series of Rydberg orbitals that have been generated as centered onto the Be atom. The resulting basis set can be described as 4s2p1d/2s1p (Be/H) + 4s4p3d. It allows to calculate Rydberg states up to n= {3,4,5} of the s, p, and d series of Rydberg states. The FCI vertical ionization potential for the same basis set and geometry amounts to 8.298 eV. Other properties such as FCI electric dipole and quadrupole moments and FCI transition dipole and quadrupole moments have also been calculated. The results provide a set of benchmark values for energies, wave functions, properties, and transition properties for the five electron BeH molecule. Most of the states have large multiconfigurational character in spite of their essentially single excited nature and a number of them present an important Rydberg-valence mixing that is achieved through the mixed nature of the particle MO of the single excitations.     

Molecular polarizabilities and induced dipole moments in molecular mechanics

Molecular polarizabilities may be divided into either atomic contributions or bond contributions. The common way to estimate molecular polarizabilities is to assign atomic or bond parameters for each atom or bond type to fit experimental or quantum mechanical results. In this study we have taken a different approach. A general formula based on MM3 force constants and bond lengths was used to compute bond polarizabilities and molecular polarizabilities. New parameters for polarizabilities are not required. A fair agreement between experimental and computed molecular polarizabilities was obtained, with a RMS deviation of 0.82 ų (11.7%) and signed average error of 0.01 ų for a broad selection of 57 molecules studied. Two methods, the many‐body interaction and the pair‐interaction approaches, have been used to study induced dipole moments using the bond polarizabilities estimated from the new formula. The pair‐interaction approximation, which involves much less computation than the many‐body interaction approach, gives a satisfactory representation of induced dipole interaction. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 813–825, 2000     

A Combined Experimental and Theoretical Electron Density Study of Intra‐ and Intermolecular Interactions in Thiourea S,S‐Dioxide

The thiourea S,S‐dioxide molecule is recognized as a zwitterion with a high dipole moment and an unusually long C? S bond. The molecule has a most interesting set of intermolecular interactions in the crystalline state—a relatively strong O???H? N hydrogen bond and very weak intermolecular C???S and N???O interactions. The molecule has C_s symmetry, and each oxygen atom is hydrogen‐bonded to two hydrogen atoms with O???H? N distances of 2.837 and 2.826 Å and angles of 176.61 and 158.38°. The electron density distribution is obtained both from Xray diffraction data at 110 K and from a periodic density functional theory (DFT) calculation. Bond characterization is made in terms of the analysis of topological properties. The covalent characters of the C? N, N? H, C? S, and S? O bonds are apparent, and the agreement on the topological properties between experiment and theory is adequate. The features of the Laplacian distributions, bond paths, and atomic domains are comparable. In a systematic approach, DFT calculations are performed based on a monomer, a dimer, a heptamer, and a crystal to see the effect on the electron density distribution due to the intermolecular interactions. The dipole moment of the molecule is enhanced in the solid state. The typical values of ρ_b and H_b of the hydrogen bonds and weak intermolecular C???S and N???O interactions are given. All the interactions are verified by the location of the bond critical point and its associated topological properties. The isovalue surface of Laplacian charge density and the detailed atomic graph around each atomic site reveal the shape of the valence‐shell charge concentration and provide a reasonable interpretation of the bonding of each atom.     

电负性研究(Ⅱ)—氢的电负性

氢的电负性值是氢元素性质的重要参数，１９３２年Pauling犤１～３犦定量确定氢的相对电负性值等于２．１，１９６１年Allred犤４，５犦用更准确的实验数据对Paul－ing电负性标度进行了修正，氢的电负性值被确定为２．２，目前这两个数值都在采用。元素的电负性值是与元素的性质紧密相关的，一个合适的电负性标度应该至少反映所有重要元素的电负性值，氢的化合物比任何其它元素都多，理应有一个基本的准确电负性值，然而一些电负性标度中却缺乏这样的数据。在Murphy等四人犤６犦最近发表的论文中，对Pauling电负性标度又进行了深入考查与…     

Molecule‐adapted basis sets optimized with a quantum Monte Carlo method

The Monte Carlo simulated annealing method is adapted to optimize correlated Gaussian‐type functions in nonrelativistic molecular environments. Starting from an atom‐centered atomic Gaussian basis set, the uncontracted functions are reoptimized in the molecular environments corresponding to the H₂O, CN^?, N₂, CO, BF, NO⁺, CO₂, and CS systems. These new molecular adapted basis sets are used to calculate total energies, harmonic vibrational frequencies, and equilibrium geometries at a correlated level of theory. The present methodology is a simple and effective way to improve molecular correlated wave functions, without the need to enlarge the molecular basis set. Additionally, this methodology can be used to generate hierarchical sequences of molecular basis sets with increasing size, which are relevant to establish complete basis set limits. © 2014 Wiley Periodicals, Inc.     

Molecule intrinsic minimal basis sets. I. Exact resolution of ab initio optimized molecular orbitals in terms of deformed atomic minimal-basis orbitals

A method is presented for expressing the occupied self-consistent-field (SCF) orbitals of a molecule exactly in terms of chemically deformed atomic minimal-basis-set orbitals that deviate as little as possible from free-atom SCF minimal-basis orbitals. The molecular orbitals referred to are the exact SCF orbitals, the free-atom orbitals referred to are the exact atomic SCF orbitals, and the formulation of the deformed "quasiatomic minimal-basis-sets" is independent of the calculational atomic orbital basis used. The resulting resolution of molecular orbitals in terms of quasiatomic minimal basis set orbitals is therefore intrinsic to the exact molecular wave functions. The deformations are analyzed in terms of interatomic contributions. The Mulliken population analysis is formulated in terms of the quasiatomic minimal-basis orbitals. In the virtual SCF orbital space the method leads to a quantitative ab initio formulation of the qualitative model of virtual valence orbitals, which are useful for calculating electron correlation and the interpretation of reactions. The method is applicable to Kohn-Sham density functional theory orbitals and is easily generalized to valence MCSCF orbitals.     

Orthogonalized atomic orbitals and the interpretation of valence bond wavefunctions

The use of orthogonalized atomic orbitals in valence bond type wavefunctions is critically examined. By analyzing a test case in detail, it is shown that the simple interpretation of such wavefunctions in terms of weighted chemical formulae can be misleading.     

Local Hartree–Fock orbitals using a three‐level optimization strategy for the energy

Using the three‐level energy optimization procedure combined with a refined version of the least‐change strategy for the orbitals—where an explicit localization is performed at the valence basis level—it is shown how to more efficiently determine a set of local Hartree–Fock orbitals. Further, a core–valence separation of the least‐change occupied orbital space is introduced. Numerical results comparing valence basis localized orbitals and canonical molecular orbitals as starting guesses for the full basis localization are presented. The results show that the localization of the occupied orbitals may be performed at a small computational cost if valence basis localized orbitals are used as a starting guess. For the unoccupied space, about half the number of iterations are required if valence localized orbitals are used as a starting guess compared to a canonical set of unoccupied Hartree–Fock orbitals. Different local minima may be obtained when different starting guesses are used. However, the different minima all correspond to orbitals with approximately the same locality. © 2013 Wiley Periodicals, Inc.     

Ab initio assessment of the electronic structure of 5α-reduced progestins

Progesterone (P) yields to 5α-reduced progestins, namely 5α-pregnanedione (DHP), tetrahydroprogesterone (THP), and allopregnanolone (ALLO-P). The geometries and electronic structure of these steroids were assessed by ab initio calculations using the 6-31G* basis set. The parameters measured were bond distances, valence angles, and dihedral angles. Likewise, the following were calculated: total energy; frontier orbitals, i.e., highest occupied molecular orbital (HOMO); lowest unoccupied molecular orbital (LUMO); dipole moment; atomic charges; and electrostatic potentials. The frontier orbitals of P were located at the π-double bond. However, the HOMO of the 5α-progestins was extended into the molecule, while the LUMO was confined at the C20 carbonyl group. The atomic charges, electronic density surfaces and electrostatic potentials showed patterns according to the stereochemical arrangement of the C3 and C20 carbonyl and hydroxyl functional groups. Interestingly, P and THP showed the larger dipole moment and high electronic density at the A-ring because the double bond and the 3α-hydroxy group, respectively. The present results might explain to some extent the metabolism of the studied progestins. Similarly, some physicochemical properties, such as dipole moments and electrostatic potentials, seem related with important biological actions such as uterine contractility and control of gonadotropin secretion. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 67: 329–338, 1998     

Interatomic and intraatomic analysis of the density matrix: applications to chlorobenzenes

Bond orders and valence indices have been evaluated employing Mayer’s definitions with orthogonalized atomic orbitals (OAO) obtained from L?wdin orthogonalization over an STO-3G basis set in anab initio formalism. It has been observed that the eigenvalues of the submatrices associated with bond order orbitals. natural hybrid orbitals and natural bond orbitals also reproduce the same values of the bond orders and the valence indices which in turn are quite close to the classical values. Bond orders obtained by a similarity transformation of theab initio density matrix differ appreciably in numerical magnitude.     

Application of Löwdin's canonical orthogonalization method to the Slater‐type orbital configuration‐interaction basis set

We apply Löwdin's canonical orthogonalization method to investigate the linearly dependent problem arising from the variational calculation of atomic systems using Slater‐type orbital configuration‐interaction (STO‐CI) basis functions. With a specific arithmetic precision used in numerical computations, the nonorthogonal STO‐CI basis is easily linearly dependent when the number of basis functions is sufficiently large. We show that Löwdin's canonical orthogonalization method can successfully overcome such problem and simultaneously reduce the dimension of basis set. This is illustrated first through an S‐wave model He atom, and then the real two‐electron atoms in both the ground and excited states. In all of these calculations, the variational bound state energies of the two‐electron systems are obtained in reasonably high accuracy using over‐redundant STO‐CI bases, however, without using extended high‐precision technique. © 2015 Wiley Periodicals, Inc.     

The Fock Matrix Analysis for Atomic Orbitals in Molecular Orbitals I. A New Look on the Covalent Bond in H2?

The character of the molecular orbitals can be better accounted for in terms of molecular adapted atomic orbitals and the Fock matrix expanded in these atomic orbital sets. A clean‐cut and unique criterion for the diradicals and the covalent bonds can be given for the molecular orbitals in both restricted and unrestricted Hartree‐Fock wavefunctions. Instead of the picture that overlap charge migrates into the bonding region, the new analysis displays another picture that the charge densities for the electrons with α and β spins give rise to two opposite spin density shifts. If the α one shifts from atom A toward atom B then it is vice versa for the β one. The spin density shifts proceed until the bonding molecular orbitals form.     

Geometry-dependent atomic charges: Methodology and application to alkanes,aldehydes, ketones,and amides

A general methodology for deriving geometry-dependent atomic charges is presented. The main ingredient of the method is a model that describes the molecular dipole moment in terms of geometry-dependent point charges. The parameters of the model are determined from ab initio calculations of molecular dipole moments and their Cartesian derivatives at various molecular geometries. Transferability of the parameters is built into the model by fitting ab initio calculations for various molecules simultaneously. The results show that charge flux along the bonds is a major contributing factor to the geometry dependence of the atomic charges, with additional contributions from fluxes along valence angles and adjacent bonds. Torsion flux is found to be smaller in magnitude than the bond and valence angle fluxes but is not always unimportant. A set of electrostatic parameters is presented for alkanes, aldehydes, ketones, and amides. Transferability of these parameters for a host of molecules is established to within 3 ?5% error in the predicted dipole moments. A possible extension of the method to include atomic dipoles is outlined. With the inclusion of such atomic dipoles and with the set of transferable point charges and charge flux parameters, it is demonstrated that molecular electrostatic potentials as well as electrostatic forces on nuclei can be reproduced much better than is possible with other models (such as potential derived charges). © 1995 by John Wiley & Sons, Inc.