Second-order shakeup terms in electron propagator calculations on F2 and H2O2

A new approximation of the electron propagator includes second-order terms in the superoperator Hamiltonian matrix that are not present in the nondiagonal, renormalized method (NR2) of J. Chem. Phys. 108 , 1008 (1998). These terms are not difficult to calculate and resemble expressions obtained with various propagator and coupled-cluster formalisms. Their inclusion produces larger, more accurate vertical ionization energies of F₂ and H₂O₂. Cancellations between the second-order terms considered here and the first-order terms present in the NR2 method reduce the importance of shakeup operators in describing the first few cationic states of these molecules. Dyson orbitals for ²A final states of H₂O₂ exhibit qualitatively important mixings between canonical orbitals. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 69: 175–182, 1998     

Investigation of ethynylpyridines using the electron propagator theory

The optimized geometry, dipole moment, and HOMO–LUMO gap for three monoethynylpyridines and six diethynylpyridines have been computed using DFT/B3LYP/6‐311++g(3df,3p) level of theory, and the first 11 vertical ionization energies and electron affinity of these compounds have been calculated using various electron propagator decouplings. The outer valence Green's function approximant of the electron propagator theory offers closest agreement with experimental photoelectron spectrum, and the results for structural suitability, dipole moment, HOMO–LUMO gap, ionization energies, and electron affinity indicate that 2‐ethynylpyridine among monoethynylpyridines and 2,6‐diethynylpyridine among diethynylpyridines may be useful precursors for the preparation of conducting polymers. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem 112:426–439, 2012     

Acidity predictions in an electron propagator approach

An alternative way to calculate vertical ionization potentials (VIP) and vertical electron affinity (VEA) is the application of Koopman's theorem, using the electron propagator theory. In the present work, the results of the application of this theorem using the electron propagator formalism have been compared with the experiment in order to validate different basis set. Using the basis set with the best performance, the acidity tendencies in some substituted acetic acid molecules have been analyzed by correlating the proton affinity (PA) with molecular electronegativity (χ) and hardness (η); these last indexes were obtained from the calculated VIP and VEA considering the finite difference approximation. The above correlations were compared with equivalent correlations using the energy of the frontier Hartree–Fock orbitals and the corresponding Kohn–Sham orbitals, which were calculated with the B3LYP‐DFT procedure. The results indicate that the electron propagator theory could be an interesting alternative to evaluate reactivity indexes, since this theory gives reliable values of VIP and VEA. It was also found that (i) the VIP values are very close to experiment, with only a 0.38% of error; (ii) acceptable results are inferred for VEA; (iii) a triple zeta quality function works quite well in these calculations, and particularly the 6‐311G(d,p) basis set is the best, as it had been reported; and (iv) using the depronation energy (DPE), good results were obtained in the correlations δDPE‐VEA and δDPE‐χ. The results tested that P3 approximation in the electron propagator approach can be a new and interesting alternative in predicting VIP, VEA, and some reactivity indexes, such as χ and η, at least for the compounds studied. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Possible DNA damage by oxidation products of guanine: A density functional and electron propagator theoretical study

The stability and reactivity of seven guanine oxidation products (GOP), which contain 8‐oxo‐7,8‐dihydroguanine (8‐oxoG) have been studied and compared with the four nucleobases, such as adenine, cytosine, guanine, and thymine. It has been possible with the use of density functional theory and electron propagator theory (EPT), to evaluate their relation with certain ionization induced process, which produce damage to DNA. Using the application of Koopmans's theorem, EPT provides alternatively a reliable way to calculate the vertical ionization potential (VIP) and vertical electron affinity (VEA). This process has been used to obtain other reactivity indexes, such as: electronegativity and hardness. In the nucleobases and the GOP studies, we observed the following: guanine and 8‐oxo‐7,8‐dihydroguanine were the lowest VIP, and 8‐oxoG was the lowest hardness. This let us confirm that these species are the most susceptible to change their electron densities and transform in other GOP. Particularly, the GOP, Sp(R), Sp(S), and Z were the highest VIP. It allows us to say that they are the most stable. Z and Iz were the highest VEA; this suggests that they have a big capacity to accept electrons and form anionic centers in DNA. The GOP, which was considered in this study, showed hardness values between 9.1 and 10.4 because of π‐conjugation; therefore, these GOP could be good candidates to participate in DNA transversions. © 2012 Wiley Periodicals, Inc.     

A generalized any-particle propagator theory: Calculations of nucleon's binding energies

Generalized one-particle propagator calculations were performed for fermions in atoms: neutrons, protons, and electrons. For this purpose, multicomponent Hartree-Fock equations were implemented using Gaussian basis sets where, for nucleons, we consider a non-Coulombic interaction, through a two-term Yukawa scalar potential and the interaction between electrons and the electrons with positive charge (protons) through a Coulombic potential. The strategy for evaluating the required interaction integrals follows Obara-Saika and Head-Gordon recurrence relations combined with the generalized Boys function suggested by Ten-no. Calculations on the isotopes ²H, ³H, ³He, ⁴He, ⁶Li, ⁶Be, ⁷Li, and ⁸Be were realized to test the accuracy of Koopmans' approximation and a second-order generalized one-particle propagator. Yukawa potentials were parametrized to reproduce nuclear properties as kinetic energies and radial distributions of density. These potentials produced the reference nuclear Hartree-Fock calculations on which fully ab initio propagator calculations were performed for these non-Coulombic potentials. This allowed us to explore the electronic structure of isotopes in an extended nucleus context.     

Electron propagator theory study of N-/O-methylglycine conformers

The low-lying conformers of N-/O-methylglycine are studied by ab initio calculations at the B3LYP, MP3, and MP4(SDQ) levels of theory with the aug-cc-pVDZ basis set. The conformers having the intramolecular hydrogen bonds N-H...O=C or O-H...N are more stable than the others. Vertical ionization energies for the valence molecular orbitals of each conformer predicted with the electron propagator theory in the partial third-order quasiparticle approximation are in good agreement with the experimental data available in the literatures. The relative energies of the conformers and comparison between the simulated and the experimental photoelectron spectra demonstrate that there are at least three and two conformers of N- and O-methylglycine, respectively, in the gas-phase experiments. The intramolecular hydrogen bonding O-H...N effects on the molecular electronic structures are discussed for the glycine methyl derivatives, on the basis of the ab initio electronic structure calculations, natural orbital bond, and atoms-in-molecules analyses. The intramolecular hydrogen bonding O-H...N interactions hardly affect the electronic structures of the O-NH2-CH2-C(=O)-O-CH3 and alpha-methylated NH2-CH2-C(CH3)OOH conformers, while the similar intramolecular interactions lead to the significantly lower-energy levels of the highest occupied molecular orbitals for the N-(CH3-NH-CH2-COOH) and beta-methylated (NH2-CH2-CH2-COOH) conformers.     

Recursive intermediate factorization and complete computational linearization of the coupled-cluster single,double, triple,and quadruple excitation equations

Summary The nonlinear CCSDTQ equations are written in a fully linearized form, via the introduction of computationally convenient intermediates. An efficient formulation of the coupled cluster method is proposed. Due to a recursive method for the calculation of intermediates, all computational steps involve the multiplication of an intermediate with aT vertex. This property makes it possible to express the CC equations exclusively in terms of matrix products which can be directly transformed into a highly vectorized program.This work has been supported by the U.S. Air Force Office of Scientific Research, Grant No. 90-0079     

Treatment of Be+ ( 1 s? 1)2S Auger resonance with different decouplings of the dilated electron propagator

The diagonal 2ph-TDA and quasiparticle decouplings of the dilated electron propagator (based on an underlying bi-variational SCF) are utilized to calculate energy and width of the Be⁺(1_s ⁻¹)²S Auger resonance for the first time. Comparison with experimental and other theoretical results reveals that the renormalized infinite order diagonal 2ph-TDA decoupling seems to offer a less balanced approach to the treatment of resonances than the second-order decoupling. The diagonal quasiparticle approximation to the self energy is seen to offer an effective and economic alternative to the non-diagonal propagator calculations.     

Improved description of the orbital relaxation effect by practical use of the self‐energy

The self‐energy shift in the orbital relaxation (OR) term of the polarization propagator complete through the second‐order is presented. In combination with the optimal damping parameter in the OR term, the modified propagator produces the excitation energy of the coupled‐cluster with singles and doubles (CCSD) accuracy. The self‐energy shift requires the floating‐point operation of , where N refers to the magnitude of the molecular size. Because the second‐order polarization propagator requires the floating‐point operation of , the additional computational effort to construct the self‐energy is negligibly small. Numerical results are shown for several molecules including glycine, 2,3,5,6‐tetrafluorobenzene, and naphthalene, and promising agreements with those of CCSD are confirmed within less than 0.2 eV. The basis set dependence is also tested for the water molecule using aug‐cc‐pV NZ (N = D–7), where this newly developed approach mimics the behavior of the CCSD values. The self‐energy shifting for the second‐order response matrix in combination with the use of a dumping parameter is efficiently implemented for calculations of medium‐sized molecular systems, including glycine and naphthalene. The developed approach provides CCSD‐like accuracy at a more affordable computational expense. © 2014 Wiley Periodicals, Inc.     

Modification of optimal complete active space choices for the multiconfigurational spin‐tensor electron propagator method for ionization potentials

The multiconfigurational spin tensor electron propagator (MCSTEP) method was developed as an implementation of electron propagator/single particle Green's function methods. MCSTEP was specifically designed for open‐shell and highly correlated (nondynamically correlated) initial states. Ionization or electron attachment is always from a state of pure spin symmetry to a state of pure spin symmetry even if the initial state is open shell. MCSTEP can be used as well for molecules with initial states that can be accurately described by a single determinant‐based theory. The initial state that is used in MCSTEP is typically a small complete active space (CAS) multiconfigurational self‐consistent field (MCSCF) state. We previously examined different small CAS choices for MCSTEP initial states and have developed a generally workable scheme. This article further examines some different ways to choose the CAS for MCSTEP. With several logical CAS choices, we have calculated the low‐lying vertical MCSTEP ionization potentials (IPs) of C₂, N₂, linear H₂O, O₂, CH₂, and NH₂, comparing them with large multireference configuration interaction (MRCI) calculations. We conclude that generally a small modification and extension of our previous schemes for choosing the MCSTEP CAS gives IPs that most effectively mimics the results of large scale MRCI IPs in general. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Testing theory beyond molecular structure: Electron density distributions of complex molecules

The electron density distribution (EDD) of a molecular system can be determined experimentally from elaborate X‐ray diffraction measurements or calculated with quantum mechanical methods: This provides a unique opportunity for mutual validation of the experimental and theoretical methods—a validation that goes far beyond comparison of molecular structures. Two examples of complex molecular systems of biologic relevance are presented. The first is the cocrystallized complex of betaine, imidazole, and picric acid, 1, which is a 75‐atom molecular complex serving as a model for the active site in the serine proteases class of enzymes, the so‐called catalytic triad. For 1 the experimental charge density was determined by combined modeling of single crystal synchrotron X‐ray and neutron diffraction data measured at 28(1) K, and it is compared with ab initio theoretical calculations at the B3LYP/6‐311G(d,p) level of theory. Overall, the agreement is good, but in one strong N? H? O hydrogen bond clear differences are observed. The second example concerns the EDD of the mixed valence trinuclear oxo‐centered iron carboxylate, [Fe₃O(OOCC(CH₃)₃)₆(NC₅H₅)₃], 2. This molecule contains 133 atoms (542 electrons) including three open‐shell iron atoms, and the experimental investigation is based on synchrotron X‐ray diffraction data. Calculations in the experimental geometry at the commonly used UB3LYP/LanL2DZ level of theory are not able to reproduce a number of experimentally observed electron density features. In particular, the sp³‐like distribution on the central oxygen atom and the electron deformations on the iron centers are at variance with experiment. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

Equivalent orbitals for multiconfigurational spin‐tensor electron propagator method (MCSTEP): The vertical ionization potentials of B,NO, CF,and OF

The multiconfigurational spin tensor electron propagator method (MCSTEP) was developed as an implementation of electron propagator/single particle Green's function methods for ionization potentials (IPs) and electron affinities (EAs). MCSTEP was specifically designed for open shell and highly correlated (nondynamically correlated) initial states. For computational efficiency the initial state used in MCSTEP is typically a small complete active space (CAS) multiconfigurational self‐consistent field (MCSCF) state. If in a molecule there are some degenerate orbitals which are not fully or half occupied, usual MCSCF calculations will make these orbitals inequivalent, i.e., the occupied ones will be different from the nonoccupied ones, so that the degeneracy is broken. In this article, we use a state averaged MCSCF method to get equivalent orbitals for the initial state and import the integrals into the subsequent MCSTEP calculations. This gives, in general, more reliable MCSTEP vertical IPs. © 2008 Wiley Periodicals, Inc., 2008     

A variational coupled cluster theory for closed shells using a propagator modification procedure

A general partial summation method for including arbitrary classes of diagrams to all orders in the coupled cluster based size consistent energy functional for closed shell states is developed. Since the various reduced density matrices which appear in the energy functional are essentially the time-independent analogues of the corresponding many body Green functions, it is possible to derive Dyson-like equations for these quantities. By expanding the associated proper self energy parts in terms of the T-amplitudes, one can carry out partial summations in the reduced density matrices and thus in energy. At a higher level, higher order terms in a proper self energy can also be generated by renormalizing the internal propagators in it, and considering only the irreducible self-energy terms.     

Electron propagator theory calculations of molecular photoionization cross sections: the first-row hydrides

Together with ionization potentials, cross sections provide valuable information for the interpretation of photoelectron spectra. We have developed a program to perform ab initio calculations of photoionization cross sections within the electric dipole approximation using electron propagator theory. Applications to the first-row hydrides CH(4), NH(3), H(2)O, and HF, using several approximations for the propagator self-energy and the plane-wave and orthogonalized-plane-wave approximations to represent the photoelectron, as well as comparison to experimental data, are presented. This program is implemented within the quantum chemistry package GAUSSIAN.     

电子断层三维重构技术在材料微观结构分析中的应用

电子断层三维重构技术是在透射电镜基础上发展起来的,用以解析材料三维结构的一种技术。本文以美国FEI公司Tecnai G~2 F20透射电镜三维重构系统——Xplore 3D系统为例,探讨了样品制备与取向的选择、样品漂移问题的成因与校正、降低缺失锲存在造成的模型失真及空间分辨率的提高等问题,从以上四个方面详细介绍电子断层三维重构技术的要点及在材料微观结构方面的应用经验。     

Electron and positron pair emission by low energy positron impact on surfaces

The emission of electron pairs from surfaces has the power to reveal details about the electron–electron interaction in condensed matter. This process, stimulated by a primary electron or photon beam, has been studied both in experiment and theory over the last two decades. An additional pathway, namely positron–electron pair emission, holds the promise to provide additional information. It is based on the notion that the Pauli exclusion principle does not need to be considered for this process.We have commissioned a laboratory based positron source and performed a systematic study on a variety of solid surfaces. In a symmetric emission geometry we can explore the fact that positron and electron are distinguishable particles. Following fundamental symmetry arguments we have to expect that the available energy is shared unequally among positron and electron. Experimentally we observe such a behavior for all materials studied. We find an universal feature for all materials in the sense that on average the positron carries a larger fraction of the available energy. This is qualitatively accounted for by a simplified scattering model. Numerical results, which we obtained by a microscopic theory of positron–electron emission from surfaces, reveal however that there are also cases in which the electron carries more energy. Whether the positron or the electron is more energetic depends on details of the bound electron state and of the emission geometry. The coincidence intensity is strongly material dependent and there exists an almost monotonic relation between the singles and coincidence intensity. These results resemble the findings obtained in electron and photon stimulated electron pair emission. An additional reaction channel is the emission of an electron pair upon positron impact. We will discuss the energy distributions and the material dependence of the coincidence signal which shows similar features as those for positron–electron pairs.     

透射电镜三维重构确定金纳米粒子的立体结构

贵金属纳米颗粒具有局域表面等离激元这一特性使其具有丰富的光学性质,而这一特性受制于纳米颗粒所形成的立体几何形状,而透射电镜和扫描电镜的二维图像不能真切地观测和确定纳米颗粒所形成的立体几何结构。透射电镜三维重构技术可作为一种确定纳米颗粒立体结构的直观有效的方法。本文利用透射电镜的三维重构技术,选择合适的参数进行二维图像的采集、图像匹配对中及重构、立体模型的构建,从而通过构建的模型对两种金纳米颗粒样品的不同几何形状所产生的边界形态进行了确认和分析。     

Monte Carlo Simulation of Electron Transport and X-Ray Generation. I. Electron Elastic and Inelastic Scattering

We present different theoretical approaches to determine differential cross sections for elastic and inelastic interactions of electrons. These cross sections are the basic ingredients for accurate Monte Carlo simulation of electron transport in matter. The considered models range from simple analytical approximations employed in early calculations to purely numerical differential cross sections described by large databases calculated with state-of-the-art theory.     

Photo-Induced Intermolecular Electron Transfer-Effect of Acceptor Molecular Structures

Photo-induced electron transfer versus molecular structure of acceptors is investigated using ultrafast time-resolved transient grating spectroscopy. Typical laser dyes Rhodamine 101 (Rh101) and Rhodamine 6G (Rh6G) in electron donor solvent-aniline are adopted as the objects. The forward electron transfer time constant from aniline to the excited singlet state of two Rhodamine dyes and subsequent back electron transfer from two dyes to aniline are measured. The experimental results denote that Rh6G presents faster electron transfer rates with aniline in both forward electron transfer and back electron transfer processes. With chemical calculation and qualitative analysis, it is found that the flexible molecular geometry of Rh6G leads to stronger electron coupling with donor solvent and further gives rise to larger electron transfer rates.     

Electronic spectra of finite polyenes and polyacetylene obtained by electron and polarization propagator calculations

Summary Polyacetylene and the short polyenes are used as common test cases for demonstrating the capability of electron and polarization propagator theory to simulate various types of molecular electronic spectra.Dedicated to Jan Linderberg on his 60th birthday