The QDPT CI approach for excitation and ionization spectra: A test on the CO molecule

A detailed analysis of the performance of the QDPT CI method, which is based on the application of quasi-degenerate perturbation theory (QDPT ) within the configuration interaction (CI ) approach, is presented by comparing the QDPT CI results and those obtained from the diagonalization of the perturbative space. Calculations of core and valence excitation and ionization spectra of the CO molecule, which is a strongly correlated system, are carried out using a DZ basis set and the accuracy of the eigenvalues and eigenfunctions is checked through the comparison of excitation and ionization energies and intensities. Excellent agreement is obtained for all values, indicating the computational advantages of QDPT CI which is able to explore wide energy ranges without requiring the diagonalization of large spaces. © 1995 John Wiley & Sons, Inc.     

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.     

Electron momentum spectroscopy of trisubstituted amines: The valence shell orbitals of triethylamine

The ionization potentials of the valence shell orbitals (up to 40 eV) of triethylamine have been measured by means of the binary (e,2e) technique. Satellite structure, due to transitions to ionic excited states, has been observed in the outer valence shell for binding energies larger than 15 eV. The electron momentum distributions of the valence orbitals have been measured on ionization peaks corresponding to main and satellite transitions. Results are compared with SCF calculations. The electron momentum distribution of the most external orbital, formed mostly by the N 2p lone pair, is discussed in detail.     

Investigation into the valence electronic structure of norbornene using electron momentum spectroscopy, Green's function, and density functional theories

Results of a study of the valence electronic structure of norbornene (C(7)H(10)), up to binding energies of 30 eV, are reported. Experimental electron momentum spectroscopy (EMS) and theoretical Green's function and density functional theory approaches were utilized in this investigation. A stringent comparison between the electron momentum spectroscopy and theoretical orbital momentum distributions found that, among the tested models, the combination of the Becke-Perdew functional and a polarized valence basis set of triple-zeta quality provides the best representation of the electron momentum distributions for all 19 valence orbitals of norbornene. This experimentally validated model was then used to extract other molecular properties of norbornene (geometry, infrared spectrum). When these calculated properties are compared to corresponding results from independent measurements, reasonable agreement is typically found. Due to the improved energy resolution, EMS is now at a stage to very finely image the effective topology of molecular orbitals at varying distances from the molecular center, and the way the individual atomic components interact with each other, often in excellent agreement with theory. This will be demonstrated here. Green's Function calculations employing the third-order algebraic diagrammatic construction scheme indicate that the orbital picture of ionization breaks down at binding energies larger than about 22 eV. Despite this complication, they enable insights within 0.2 eV accuracy into the available ultraviolet emission and newly presented (e,2e) ionization spectra. Finally, limitations inherent to calculations of momentum distributions based on Kohn-Sham orbitals and employing the vertical depiction of ionization processes are emphasized, in a formal discussion of EMS cross sections employing Dyson orbitals.     

A combined freeze-and-cut strategy for the description of large molecular systems using a localized orbitals approach

A technique to reduce the computational effort in calculating ab initio energies using a localized orbitals approach is presented. By exploiting freeze strategy at the self-consistent field (SCF) level and a cut of the unneeded atomic orbitals, it is possible to perform a localized complete active space (CAS-SCF) calculation on a reduced system. This will open the possibility to perform ab initio treatments on very large molecular systems, provided that the chemically important phenomena happen in a localized zone of the molecule. Two test cases are discussed, to illustrate the performance of the method: the cis-trans interconversion curves for the (7Z)-13 ammoniotridec-7-enoate, which demonstrates the ability of the method to reproduce the interactions between charged groups; and the cisoid-transoid energy barrier for the aldehydic group in the C13 polyenal molecule.     

Norbornane: an investigation into its valence electronic structure using electron momentum spectroscopy, and density functional and Green's function theories

We report on the results of an exhaustive study of the valence electronic structure of norbornane (C(7)H(12)), up to binding energies of 29 eV. Experimental electron momentum spectroscopy and theoretical Green's function and density functional theory approaches were all utilized in this investigation. A stringent comparison between the electron momentum spectroscopy and theoretical orbital momentum distributions found that, among all the tested models, the combination of the Becke-Perdew functional and a polarized valence basis set of triple-zeta quality provides the best representation of the electron momentum distributions for all of the 20 valence orbitals of norbornane. This experimentally validated quantum chemistry model was then used to extract some chemically important properties of norbornane. When these calculated properties are compared to corresponding results from other independent measurements, generally good agreement is found. Green's function calculations with the aid of the third-order algebraic diagrammatic construction scheme indicate that the orbital picture of ionization breaks down at binding energies larger than 22.5 eV. Despite this complication, they enable insights within 0.2 eV accuracy into the available ultraviolet photoemission and newly presented (e,2e) ionization spectra, except for the band associated with the 1a(2) (-1) one-hole state, which is probably subject to rather significant vibronic coupling effects, and a band at approximately 25 eV characterized by a momentum distribution of "s-type" symmetry, which Green's function calculations fail to reproduce. We note the vicinity of the vertical double ionization threshold at approximately 26 eV.     

Molecular MC SCF calculations by direct minimization

Calculations on the H₂CO ground state and excited states at the equilibrium planar geometry, using the single excitation MC SCF method, are described for two basis sets. The results obtained by a full calculation including all the electrons are compared with those obtained by a pseudopotential version of the method including only the valence electrons. The results agree quite well both with each other and with the observed values.     

A SCF Xα MO study of the optical and Ti2p core photoemission spectra of TiCl4

SCF Xα MO calculations on the ground state and optical excitation transition states of TiCl₄ accurately predict the energies of its UV absorption peaks. Calculations on the Ti2p core ion state and associated transition states indicate that the recently observed low energy (4.0 eV) Ti2p satellite arises from ligand to metal charge transfer excitations while the satellite at high energy (9.4 eV), similar to those previously observed in Ti(IV) compounds, can be attributed to transitions from the highest filled orbitals to empty orbitals with Cl3pTi4s. 4p antibonding character.     

An ab initio study of the ground and excited states of p-benzoquinone

The results of anab initio SCF calculation for the ground state and CI calculations for the excited states of p-benzoquinone are presented and discussed. A minimum basis set of Slater type orbitals was employed and the CI calculations were performed by considering single excitations from valence to virtual SCF molecular orbitals. The convergence of the calculated excitation energies is studied as a function of the number of orbitals used in the CI calculations. These calculations explain quite well the experimental results.     

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.     

Coupled-cluster studies of the electronic excitation spectra of silanes

The electronic excitation spectra of unsubstituted linear silanes (n-Si(m)H(2m+2), m = 1-6), cyclopentasilane (c-Si5H10), and neopentasilane (neo-Si5H12) have been studied at the coupled-cluster approximate singles and doubles (CC2) level using Dunning's quadruple-zeta basis sets augmented with diffuse functions (aug-cc-pVQZ). Comparisons with measured ultraviolet spectra for Si2H6 and n-Si3H8 show that CC2 calculations using these basis sets yield excitation energies in good agreement with experiment. The calculated excitation thresholds for Si2H6 and n-Si3H8 of 7.61 and 6.68 eV are only 0.05 eV larger than the gas-phase values of 7.56 and 6.63 eV, respectively. For n-Si4H10, n-Si5H12, and neo-Si5H12, the calculated excitation thresholds of 6.51, 6.14, and 6.87 eV for the lowest dipole-allowed transitions are about 0.4 eV larger than the corresponding liquid-phase data of 6.05, 5.77, and 6.53 eV; the discrepancy can mainly be attributed to solvent effects. The obtained excitation thresholds for n-Si6H14 is 5.85 eV, whereas no experimental data are available for its optical gap. Calculations using the Karlsruhe triple-zeta valence basis sets augmented with single and double sets of polarization functions show that very large basis sets augmented with diffuse functions are needed for obtaining accurate excitation energies. The optical gaps for silanes obtained using the triple-zeta polarization basis sets were found to be 0.4 and 0.2 eV larger than those obtained using Dunning's quadruple-zeta basis sets. Excitation thresholds calculated at density functional theory levels using generalized gradient approximation are 0.7-1.0 eV smaller than the experimental values and by employing hybrid functionals they are 0.3-0.4 eV below the experimental thresholds. By adding the present basis-set correction and environmental effects to the previously calculated CC2 value for the excitation threshold of the Si29H36 silicon nanocluster, the extrapolated absorption threshold is 4.0 eV as compared to the recently reported experimental value of 3.7 eV.     

Study of the photoelectron and electron momentum spectra of cyclopentene using benchmark Dyson orbital theories

A complete study of the valence electronic structure and related electronic excitation properties of cyclopentene in its C(s) ground state geometry is presented. Ionization spectra obtained from this compound by means of photoelectron spectroscopy (He I and He II) and electron momentum spectroscopy have been analyzed in details up to electron binding energies of 30 eV using one-particle Green's function (1p-GF) theory along with the outer-valence (OVGF) and the third-order algebraic diagrammatic construction [ADC(3)] schemes. The employed geometries derive from DFT/B3LYP calculations in conjunction with the aug-cc-pVTZ basis set, and closely approach the structures inferred from experiments employing microwave spectroscopy or electron diffraction in the gas phase. The 1p-GF/ADC(3) calculations indicate that the orbital picture of ionization breaks down at electron binding energies larger than approximately 17 eV in the inner-valence region, and that the outer-valence 7a' orbital is also subject to a significant dispersion of the ionization intensity over shake-up states. This study confirms further the rule that OVGF pole strengths smaller than 0.85 foretell a breakdown of the orbital picture of ionization at the ADC(3) level. Spherically averaged (e, 2e) electron momentum distributions at an electron impact energy of 1200 eV that were experimentally inferred from an angular analysis of EMS intensities have been interpreted by comparison with accurate simulations employing ADC(3) Dyson orbitals. Very significant discrepancies were observed with momentum distributions obtained from several outer-valence ionization bands using standard Kohn-Sham orbitals.     

Valence type vacant orbitals for configuration interaction calculations

A method is proposed to determine the valence type vacant orbitals, which are suitable for CI calculations and for the initial guess orbitals in MC SCF calculations. The method was applied to calculate the ionization energies of series of molecules and to draw the potential energy curves of various states of N₂ and N⁺₂.     

Calculation of excitation and photoionization spectra by quasi-degenerate perturbation theory

The ability of the QDPT CI method, in which the quasi-degenerate perturbation theory is applied within the configuration interaction (CI ) approach, in dealing with the calculation of excitation and photoionization spectra is shown through an overview of comparisons between the QDPT CI values and the full CI ones in the same basis set. A direct comparison with the experimental data is given for the core photoelectron spectrum of Ne and the valence photoelectron spectra of Ar and HCl. The quality of the results obtained is very satisfactory, both as regards the energies and the wave functions, indicating the validity and the flexibility of the method that can confidently be applied also in cases where strong correlation effects are present. © 1994 John Wiley & Sons, Inc.     

Threshold-photoelectron spectroscopic study of methyl-substituted hydrazine compounds

The valence shell electronic structures of methylhydrazine (CH(3)NHNH(2)), 1,1-dimethylhydrazine ((CH(3))(2)NNH(2)) and tetramethylhydrazine ((CH(3))(4)N(2)) have been studied by recording threshold and conventional (kinetic energy resolved) photoelectron spectra. Ab initio calculations have been performed on ammonia and the three methyl substituted hydrazines, with the structures being optimized at the B3-LYP/6-31+G(d) level of theory. The ionization energies of the valence molecular orbitals were calculated using the Green's function method, allowing the photoelectron bands to be assigned to specific molecular orbitals. The ground-state adiabatic and vertical ionization energies, as determined from the threshold photoelectron spectra, were IE(a) = 8.02 +/- 0.16 eV and IE(v) = 9.36 +/- 0.02 eV for methylhydrazine, IE(a) = 7.78 +/- 0.16 eV and IE(v) = 8.86 +/- 0.01 eV for 1,1-dimethylhydrazine and IE(a) = 7.26 +/- 0.16 eV and IE(v) = 8.38 +/- 0.01 eV for tetramethylhydrazine. Due to the large geometry change that occurs upon ionization, these IE(a) values are all higher than the true thresholds. New features have been observed in the inner valence region and these have been compared with similar structure in the spectrum of hydrazine. The effect of resonant autoionization on the threshold photoelectron yield is discussed. New heats of formation (Delta(f)H) are proposed for the three hydrazines on the basis of G3 calculations: 107, 94, and 95 kJ/mol for methylhydrazine, 1,1-dimethyhydrazine and tetramethylhydrazine, respectively. The previously reported Delta(f)H for tetramethylhydrazine is shown to be erroneous.     

Application of the MC SCF method to the π → π* excitation energies of ethylene

The MC SCF method is employed to calculate the N → T and N → V π → π^* vertical excitation energies of ethylene. To obtain accurate excitation energies it is found to be necessary to utilize an expanded valence space containing two π and two π^* orbitals. Relatively small MC SCF calculations, allowing at most one-electron excitations from the sigma space, are found to yield excitation energies and spatial extents of the excited states in excellent agreement with the predictions of large multi-reference or iterative-natural-orbital CI calculations. These results show that within an MC SCF framework σ-σ correlation is unimportant for describing the π → π^* processes. We also conclude that the neglect of the effects of unlinked cluster terms in some of the CI calculations may have introduced small, but important, errors in the excitation energies and predictions of the spatial extent of the V state.     

Valence electronic structure of CH3F and CH3Cl: electron momentum distributions and separation energies

Fluoromethane (CH₃F) has been studied by binary (e,2e) coincidence spectroscopy at 1200 eV using non-coplanar symmetric kinematics. Separation energy spectra have been determined in the energy range up to 50 eV at azimuthal angles of 0° and 9°. The separation energy spectra and electron momentum distributions measured for the valence orbitals of CH₃F and CH₃Cl are compared with the results of calculations employing SCF wavefunctions and outer valence as well as extended 2ph—TDA Green function methods. Electron density and momentum density maps have been generated for all valence orbitals of both molecules using the SCF wavefunctions and are used to explain differences in the bonding properties of the halomethanes investigated here.     

Low-lying quartet electronic states of nitrogen dioxide

The environmentally active molecule nitrogen dioxide (NO2) has been systematically studied using high level theoretical methods. The electronic ground state and the low-lying quartet states of NO2 have been investigated. Single reference restricted open-shell self-consistent field (SCF), complete active space SCF (CASSCF), spin-restricted (R) and spin-unrestricted (U) configuration interaction with single and double excitations (CISD), coupled cluster with single and double excitations (CCSD), CCSD with perturbative triple excitations [CCSD(T)], and internally contracted multireference configuration interaction (ICMRCI) methods along with Dunning's correlation consistent polarized valence cc-pVXZ and augmented cc-pVXZ (where X=T,Q,5) basis sets were used in this research. At the aug-cc-pV5Z/UCCSD(T) level the classical adiabatic excitation energies (Te values) of the three lowest-lying quartet excited states were predicted to be 83.3 kcalmol (3.61 eV, 29 200 cm(-1)) for the ? 4A2 state, 93.3 kcalmol (4.05 eV, 32 600 cm(-1)) for the b 4B2 state, and 100.8 kcalmol (4.37 eV, 35 300 cm(-1)) for the c 4A1 state. The quantum mechanical excitation energies (T 0 values) were determined to be 81.6 kcalmol (3.54 eV, 28 500 cm(-1)) for the a 4A2 state and 90.7 kcalmol (3.93 eV, 31 700 cm(-1)) for the b 4B2 state. The lowest quartet linear Renner-Teller 4Pi state gives rise to the a 4A2 state with 112.8 degrees and the b 4B2 state with 124.4 degrees <(ONO) bond angles upon bending. The b state shows some peculiar behavior. Although CASSCF, RCISD, UCISD, RCCSD, UCCSD, and RCCSD(T) methods predicted the presence of a Cs equilibrium geometry (a double minimum 4A' state), SCF, UCCSD(T), and ICMRCI wave functions predicted the C2v structure for the b 4B2 state. The importance of both dynamical and nondynamical correlation treatments for the energy difference between C2v and Cs structures of b state is highlighted in this context. The c 4A1 state is predicted to have a very small bond angle of 85.8 degrees . Potential energy diagrams with respect to the bond angles of the ground state and four quartet states are presented.