The quest for best suited references for configuration interaction singles calculations of core excited states

Near edge X‐ray absorption fine structure (NEXAFS) simulations based on the conventional configuration interaction singles (CIS) lead to excitation energies, which are systematically blue shifted. Using a (restricted) open shell core hole reference instead of the Hartree Fock (HF) ground state orbitals improves (Decleva et al., Chem. Phys., 1992, 168, 51) excitation energies and the shape of the spectra significantly. In this work, we systematically vary the underlying SCF approaches, that is, based on HF or density functional theory, to identify best suited reference orbitals using a series of small test molecules. We compare the energies of the K edges and NEXAFS spectra to experimental data. The main improvement compared to conventional CIS, that is, using HF ground state orbitals, is due to the electrostatic influence of the core hole. Different SCF approaches, density functionals, or the use of fractional occupations lead only to comparably small changes. Furthermore, to account for bigger systems, we adapt the core‐valence separation for our approach. We demonstrate that the good quality of the spectrum is not influenced by this approximation when used together with the non‐separated ground state wave function. Simultaneously, the computational demands are reduced remarkably. © 2016 Wiley Periodicals, Inc.     

Effect of modified virtual orbitals on local Coulomb correlation holes in FCN

Local Coulomb correlation hole distribution functions may be used to assess the extent to which electron correlation effects are present in large scale SCF + CI wave functions. From a set of modified virtual orbitals, ordered according to their interaction with the SCF configuration, we have constructed a limited SCF + CI wave function with improved convergence characteristics with respect to that formed from the canonical virtual orbital set. These wave functions, of the same size yet with different energies, have been used to examine the range and depth of local Coulomb correlation holes in FCN. In all cases, the depth of the local Coulomb hole is no more than 10% or so of that of the corresponding Fermi hole, and the range Fermi correlation is generally less than that of Fermi correlation. This is particularly marked in the high density regions around the nuclei. The significance of our results is discussed in relation to a recent proposal for the incorporation of Coulomb correlation into the local exchange method.     

Electron correlation effects in SCF+CI wave functions—Fermi and coulomb correlation holes in HCN

The Coulomb correlation hole distribution function has been computed with respect to various reference centers in the HCN molecule, using standard SCF +CI type wave functions. The extent to which statistical correlation between unlike-spin electrons is introduced into an SCF wave function through the inclusion of configuration interaction has been assessed by an examination of the range and depth of such holes, and compared with the behavior of analogous Fermi distribution functions. Our results show that the range of Fermi correlation is consistently longer than that of the corresponding Coulomb correlation.     

Density matrix approach to orbital relaxation dynamics in ionization

Dynamical response of electrons to a hole generated during ionization is formulated in time domain with the density matrix equations in the time‐dependent unrestricted Hartree–Fock approximation. Time evolutions of orbital energies and electron‐density distributions are computed for K‐shell and M‐shell ionizations of a Na atom by taking into account nonlinear coupling of density matrices beyond linear response. When the hole is generated so slowly that the adiabatic theorem is satisfied, the simulation eventually converges to the state of a fully relaxed Na⁺ ion. A rapid generation of a K‐shell hole (within about 1 fs) leads to a breakdown of the adiabatic theorem, triggering a collective oscillation of the electrons with the period of sub‐femtoseconds. The shake‐up effect associated with strong orbital relaxation in inner‐shell ionization is manifested as a mixing of occupied and unoccupied states in the density matrix.     

Relaxed core projector-augmented-wave method

We extend the full-potential projector-augmented-wave method beyond the frozen core approximation, i.e., include the self-consistent optimization of the core charge density, in such a manner that the valence wave functions remain orthogonal to the core. The method consists of an on-the-fly repseudization of the all-electron problem, solving for the self-consistent core charge density within a spherical approximation. The key ideas in our procedure are to keep the projector functions fixed throughout the electronic minimization and to derive the new pseudopartial waves from these original projector functions, at each step of the electronic minimization procedure. Results of relaxed core calculations for atomic interconfigurational energies, structural energy differences between bulk phases of Fe, atomization energies of a subset of Pople's G2-1 set, and the Rh 3d surface core level shifts for the (log3 x log3)-Rh(111) surface at 1/3 CO coverage are presented.     

Analysis of molecular orbital wave functions in terms of valence bond functions for molecular fragments. I. Theory

A method of expansion of molecular orbital wave functions into valence bond (VB ) functions is extended to molecular fragments. The wave function is projected onto a basis of mixed determinants, involving molecular orbitals as well as fragment atomic orbitals, and is further expressed as a linear combination of VB functions, characteristic of structural formulas of the fragment but whose remaining bonds are frozen. Structural weights for the fragment are deduced from this expression. Delocalized molecular orbitals are used as a startpoint, as they are after an ordinary SCF calculation. Wave functions of medium-sized molecules may be analyzed with reasonable storage requirements in a computer.     

The interaction of alkali metal cations with oxygen-containing ligands

Ab initio SCF calculations on the interaction of Li⁺ cation with H₂O and H₂CO using two basis sets are presented. Partitioning of SCF energies of interaction into Coulomb-, exchange- and delocalization energies has been performed. Coulomb- and delocalization energies are compared with classical electrostatic and polarization energies. A detailed analysis of the calculated wave functions demonstrates that in the complexes investigated here, charge transfer is of minor importance only. Polarization of the molecules in the strong inhomogeneous field of the cation leads to complicated electron density rearrangements which can be interpreted most easily in terms of polarization of individual localized MO's.     

Physical nature of the chemical bond II. Valence atomic orbital and energy partitioning studies of linear nitriles

The valence atomic orbitals (VAO 's) of several linear nitriles are determined using non-empirical SCF –LCAO –MO wave functions expanded in a minimal (CN^?, HCN, FCN, C₂N₂), double-zeta (CN^?, HCN), or double-zeta + polarization (HCN) basis of Slater atomic orbitals (AO 's). The molecular energy of each system (except the double-zeta + polarization HCN system) is partitioned according to the procedure of Ruedenberg to obtain numerical values of nitrile C and N atomic and C?N bond components of the energy. In addition, the nitrile results are compared with minimal AO basis results obtained previously by other authors for homonuclear diatomics, diatomic hydrides and H₂O. The numerical data are used to test the internal self-consistency of the various definitions entering the partitioning method, i.e. whether or not analogous quantities assume similar values in chemically similar situations. The analysis of nitrile SCF –MO wave functions in terms of the set of VAO 's characteristic of the system under consideration is shown to be a promising approach to the problem of extracting useful information from the wave functions. In general, numerical results for the nitrile systems studied are fairly consistent with the concepts on which the partitioning method is based: promotion, quasi-classical interaction, sharing penetration, sharing interference and charge transfer. However, the VAO expansions for several energy components need to be investigated further and possibly revised.     

Satellite bands in the X-Ray Photoelectron spectrum of formaldehyde

The energies and intensities of satellite bands in the X-ray photoelectron spectrum of formaldehyde have been calculated using molecular wave functions obtained from a combination of multi-configuration self-consistent field theory and configuration-mixing formalism. The intensities of monopole allowed shake-up processes have been obtained from the sudden approximation formula whose derivation by an alternative route is detailed. It is found that the satellite spectra of carbon and oxygen hole states are predicted to differ substantially; giving insight into the atomic orbital character of the molecular orbitals concerned. The significance of the appearance of satellite bands on a quantitative analytical use of X-ray photoelectron spectroscopy is pointed out. The possibilities and condition for the appearance of monopole forbidden valence shell shape-up processes is discussed.     

Conformational preferences of XONO2 systems (X = H,F, Cl,CH3) from ab initio techniques

Ab initio SCF and MP2 methods have been used to calculate the geometries and relative energies of both planar and nonplanar rotamers of several nitrates with a DZP basis set. The planar arrangement of atoms is found to be the lower energy configuration in all cases. The interconversion between rotamers is analyzed by partitioning the contributions to the total SCF energy in a variety of ways.     

Frozen local hole approximation

The frozen local hole approximation (FLHA) is an adiabatic approximation which is aimed to simplify the correlation calculations of valence and conduction bands of solids and polymers or, more generally, of the ionization potentials and electron affinities of any large system. Within this approximation correlated local hole states (CLHSs) are explicitly generated by correlating local Hartree-Fock (HF) hole states, i.e., (N-1)-particle determinants in which the electron has been removed from a local occupied orbital. The hole orbital and its occupancy are kept frozen during these correlation calculations, implying a rather stringent configuration selection. Effective Hamilton matrix elements are then evaluated with the above CLHSs; diagonalization finally yields the desired correlation corrections for the cationic hole states. We compare and analyze the results of the FLHA with the results of a full multireference configuration interaction with single and double excitations calculation for two prototype model systems, (H2)n ladders and H-(Be)n-H chains. Excellent numerical agreement between the two approaches is found. Comparing the FLHA with a full correlation treatment in the framework of quasidegenerate variational perturbation theory reveals that the leading contributions in the two approaches are identical. In the same way it could be shown that a much less demanding self-consistent field (SCF) calculation around a frozen local hole fully recovers, up to first order, all the leading single excitation contributions. Thus, both the FLHA and the above SCF approximation are well justified and provide a very promising and efficient alternative to fully correlated wave-function-based treatments of the valence and conduction bands in extended systems.     

X-ray spectra and core hole energy curves of some diatomic molecules

Soft X-ray emission spectra of the molecules CO, N₂, NO and O₂ are examined for the purpose of deriving information on their core hole energy curves. Molecular force constant and equilibrium bond lengths are determined for the core hole species C*O and N₂*, and a qualitative analysis is made for CO*, N*O, NO* and O₂*. The results show that differences of equilibrium geometries between the core hole states and the ground states are very well reproduced (better than 1 pm) by SCF calculations within the Hartree-Fock formalism. Inclusion of anharmonicity in the Franck-Condon analysis gives a small but significant effect on the best fitted value for the core hole state bond lengths (about 0.5 pm). Oxygen is binding energies determined from the X-ray spectra are shown to agree with ESCA data, in most cases within a few tenths of an eV. Calculated ΔSCF transition energies reproduce the experimental data within a few eV.     

Modification of the Roothaan equations to exclude BSSE from molecular interaction calculations

The Roothaan equations have been modified to compute molecular interactions between weakly bonded systems at the SCF level of theory without the basis set superposition error (BSSE). The increase in complication with respect to the usual SCF algorithm is negligible. Calculation of the SCF energy on large systems, such as nucleic acid pairs, does not pose any computational problem. At the same time, it is shown that a modest change in basis-set quality from 3-21G to 6-31G changes the binding energy by about 50% when computed according to standard SCF “supermolecule” techniques, while remaining practically constant when computed without introducing BSSE. Bader analysis shows that the amount of charge transferred between the interacting units is of the same order of magnitude when performed on standard SCF wave functions and those computed using the new method. The large difference between the corresponding computed energies is thus ascribed to the BSSE. © 1996 John Wiley & Sons, Inc.     

Core molecular orbital contribution to N2O isomerization as studied using theoretical electron momentum spectroscopy

Core molecular orbital contribution to the electronic structure of N2O isomers has been studied using quantum mechanical density functional theory combined with a plane wave impulse approximation method. Momentum distributions of wave functions for inner shell molecular orbitals of the linear NNO, cyclic and linear NON isomers of N2O are calculated through the (e, 2e) differential cross sections in momentum space. This is possible because this momentum distribution is directly proportional to the modulus squared of the momentum space wave function for the molecular orbital in question. While the momentum distributions of the NNO and cyclic N2O isomers demonstrate strong atomic orbital characteristics in their core space, the outer core molecular orbitals of the linear NON isomer exhibit configuration interactions between them and the valence molecular orbitals. It is suggested that the frozen core approximation breaks down in the prediction of the electronic structure of such an isomer. Core molecular orbital contributions to the electronic structure can alter the order of total energies of the isomers and lead to incorrect conclusions of the stability among the isomers. As a result, full electron calculations should be employed in the study of N2O isomerization.     

Methods of composite molecular wave functions. I. Variational principle and multiconfiguration SCF theory

The Silverstone–Stuebing variational principle for the discontinuous wave functions of one-electron systems is generalized for many-electron systems. The variational functional of energy takes real or complex value. The condition that it is real is given. Using the generalized variational principle, a multiconfiguration SCF theory for the composite molecular wave function is formulated. According to the theory, we may divide the whole space into space-filling cells, solve the SCF equations in each cell and build up the wave functions of the system by gathering the wave functions obtained in the cells. For use in the basis-set expansion method, the SCF equations are rewritten as matrix forms in which only one- and two-center integrals appear if an expansion center is located in each cell.     

Remarks on the use of the average occupation numbers in the SCF process for open shell systems

It is shown that the utilization of an average occupation number for open shell orbitals, having different occupation numbers in a degenerate wave function, tantamounts to disregard first order contributions in a CI -type expansion of the wave function. These contributions are taken into account in a SCF process that derives the differently occupied open shell orbitals as eigenfunctions of different Fock-type operators and accounts for the total symmetry.     

A theoretical investigation of electronic reorganizations accompanying core and valence ionization in simple hydrogen bonded systems

Non-empirical LCAO MO SCF calculations have been carried out on the ground state and core ionized states of some hydrogen bonded dimers, and in the particular case of H₂O the trimer has also been investigated. Comparison of absolute and relative binding energies and relaxation energies with respect to the corresponding monomers reveals that substantial changes occur in going to the associated species. The relaxation energies for a given core hole are shown to increase on going from monomer to dimer indicating that intermolecular contributions to relaxation energies are of the same sign irrespective of the sign for the shift in core binding energy. Creation of a core hole in the dimer species is shown to give rise to substantial changes in hydrogen bond energies compared with the neutral species. In the particular case of valence holes dominantly of 2s and 2p character it is shown that trends in shifts and relaxation energies parallel those for the core hole states.     

Ordinary field-theoretic methods for self-consistent wave functions which describe bond formation and dissociation. II. Commutative coupling case

Hartree–Bogoliubov–Valatin (HBV ) theory may be implemented with Lipkin Hamiltonians to obtain self-consistent BCS wave functions which describe bond formation and dissociation. These wave functions are in turn vacuua for Nambu's representation of Feynman–Dyson–Goldstone diagrammatic perturbation theory, and hence provide suitable references for the many-body treatment of correlation. Exact SCF solutions of the HBV equations are equivalent to special even-replacement MC –SCF solutions. The latter are similar to generalized valence bond theory, and require one Fock operator for each one-particle shell. The commutative coupling case of HBV theory is realized when the number-conserving renormalized one-body and number-nonconserving pairing operators commute. In this case, a set of orbital equations which involves a single Fock operator may be solved. Since this could prove to be a significant simplification for large systems, the commutative coupling and exact solutions are compared here for the fragmentation of H₂ and F₂. Results suggest that commutative coupling orbitals will be useful for the aforementioned many-body theory.     

Molecular states of atoms. II

Orbital energy parameters, previously obtained from atomic valence state energies, are used in calculating approximate wave functions for their orbitals. The radial factors of these wave functions are expressed as linear combinations of three Gaussian type orbitals with selected exponents, the coefficients being determined by normalisation and reproduction of the kinetic energy and interelectron repulsion parameters. Wave functions of universal form are obtained for the non-transition elements up to xenon. Each calculated s orbital wave function (except 1s) has a radial node, as is appropriate if there is a p orbital in the same shell with none.     

A non-emperical LCAO MO SCF investigation of the ground and core ionized states of tetrahedrane and cyclobutadiene

Non-empirical LCAO MO SCF calculations have been carried out on the ground and localised core hole states of tetrahedrane and cyclobutadiene to investigate the suitability of ESCA as a technique for distinguishing between structurally related, isomeric hydrocarbons. Differences in energy are found to be considerably magnified in the core hole state manifold and shifts in core binding energies and shake-up structures are shown to be distinctive.