Equation of motion coupled cluster methods for electron attachment and ionization potential in polyacenes
TL;DR: In this paper, the authors employed the EA/IP equation of motion coupled cluster singles and doubles (EA/IP-EOMCCSD) method for linear polyacenes and their electron affinity (EA) and ionization potential (IP) properties.
About: This article is published in Chemical Physics Letters.The article was published on 2015-11-16 and is currently open access. It has received 12 citations till now. The article focuses on the topics: Hexacene & Coupled cluster.
Citations
More filters
TL;DR: The unusual intrinsic electronic structure of the [B12X12]2- MCAs provides the basis for a molecular level understanding of their observed unique physical and chemical properties and a new paradigm for understanding the properties of these MCAs with not well-separated charges that departs from the prevailing model used for spatially separated charges.
Abstract: The stability and electron loss process of numerous multiply charged anions (MCAs) have been traditionally explained in terms of the classical Coulomb interaction between spatially separated charged groups. An understanding of these processes in MCAs with not well-separated excess charges is still lacking. We report the surprising properties and physical behavior of [B12X12]2–, X = F, Cl, Br, I, At, which are MCAs with not well-separated excess charges and cannot be described by the prevailing classical picture. In this series of MCAs, comprising a “boron core” surrounded by a “halogen shell”, the sign of the total charge in these two regions changes along the halogen series from X = F–At. With the aid of experimental photoelectron spectroscopy and highly correlated ab initio electronic structure calculations, we demonstrate that the trend in the electronic stability of these MCAs is determined by the interplay between the Coulomb (de)stabilization originating from the “boron core” and “halogen shell” and...
51 citations
TL;DR: A key aspect of this work is the determination of the frequency dependent self-energy, Σ(ω), which reveals a diagonally dominate block structure where as the non-local correlation increases, the block size increases proportionally.
Abstract: In this paper, we report on the development of a parallel implementation of the coupled-cluster (CC) Green function formulation (GFCC) employing single and double excitations in the cluster operator (GFCCSD). A key aspect of this work is the determination of the frequency dependent self-energy, Σ(ω). The detailed description of the underlying algorithm is provided, including approximations used that preserve the pole structure of the full GFCCSD method, thereby reducing the computational costs while maintaining an accurate character of methodology. Furthermore, for systems with strong local correlation, our formulation reveals a diagonally dominate block structure where as the non-local correlation increases, the block size increases proportionally. To demonstrate the accuracy of our approach, several examples including calculations of ionization potentials for benchmark systems are presented and compared against experiment.
46 citations
TL;DR: The morphological evolution of poly(3-hexylthiophene) (P3HT) and phenyl-C61-butyric acid methyl ester blends that are commensurate with the size of a BHJ thin film is examined and nanoscale geometries of the P3HT oligomers, LUMO and HOMO energy levels calculated by density functional theory are found to be invariant across the donor/acceptor interface.
Abstract: Understanding how additives interact and segregate within bulk heterojunction (BHJ) thin films is critical for exercising control over structure at multiple length scales and delivering improvements in photovoltaic performance. The morphological evolution of poly(3-hexylthiophene) (P3HT) and phenyl-C61-butyric acid methyl ester (PCBM) blends that are commensurate with the size of a BHJ thin film is examined using petascale coarse-grained molecular dynamics simulations. Comparisons between two-component and three-component systems containing short P3HT chains as additives undergoing thermal annealing demonstrate that the short chains alter the morphology in apparently useful ways: they efficiently migrate to the P3HT/PCBM interface, increasing the P3HT domain size and interfacial area. Simulation results agree with depth profiles determined from neutron reflectometry measurements that reveal PCBM enrichment near substrate and air interfaces but a decrease in that PCBM enrichment when a small amount of shor...
28 citations
TL;DR: It is demonstrated that the electronegativity of C20nH20 systems remains, to a large extent, independent of the nanotube length.
Abstract: In this paper, we apply equation-of-motion coupled cluster (EOM-CC) methods in the studies of the vertical ionization potentials (IPs) and electron affinities (EAs) for a series of single-walled carbon nanotubes (SWCNT). The EOM-CC formulations for IPs and EAs employing excitation manifolds spanned by single and double excitations (IP/EA-EOM-CCSD) are used to study the IPs and EAs of the SWCNTs as a function of the nanotube length. Several armchair nanotubes corresponding to C20nH20 models with n = 2–6 have been used in benchmark calculations. In agreement with previous studies, we demonstrate that the electronegativity of C20nH20 systems remains, to a large extent, independent of the nanotube length. We also compare IP/EA-EOM-CCSD results with those obtained with coupled cluster models with single and double excitations corrected by perturbative triples, CCSD(T), and density functional theory (DFT) using global and range-separated hybrid exchange–correlation functionals.
9 citations
TL;DR: In this article , a subspace projection scheme was proposed to evaluate resonance positions and widths at the overall cost of a single electronic structure calculation for the complex absorbing potential equation-of-motion coupled-cluster (CAP-EOM-CC) method.
Abstract: The complex absorbing potential equation-of-motion coupled-cluster (CAP-EOM-CC) method is routinely used to investigate metastable electronic states in small molecules. However, the requirement of evaluating eigenvalue trajectories presents a barrier to larger simulations, as each point corresponding to a different value of the CAP strength parameter requires a unique eigenvalue calculation. Here, we present a new implementation of CAP-EOM-CCSD that uses a subspace projection scheme to evaluate resonance positions and widths at the overall cost of a single electronic structure calculation. We analyze the performance of the projected CAP-EOM-CC scheme against the conventional scheme, where the CAP is incorporated starting from the Hartree-Fock level, for various small and medium sized molecules, and investigate its sensitivity to various parameters. Finally, we report resonance parameters for a set of molecules commonly used for benchmarking CAP-based methods, and we report estimates of resonance energies and widths for 1- and 2-cyanonaphtalene, molecules that were recently detected in the interstellar medium.
7 citations
References
More filters
TL;DR: In this paper, a detailed study of correlation effects in the oxygen atom was conducted, and it was shown that primitive basis sets of primitive Gaussian functions effectively and efficiently describe correlation effects.
Abstract: In the past, basis sets for use in correlated molecular calculations have largely been taken from single configuration calculations. Recently, Almlof, Taylor, and co‐workers have found that basis sets of natural orbitals derived from correlated atomic calculations (ANOs) provide an excellent description of molecular correlation effects. We report here a careful study of correlation effects in the oxygen atom, establishing that compact sets of primitive Gaussian functions effectively and efficiently describe correlation effects i f the exponents of the functions are optimized in atomic correlated calculations, although the primitive (s p) functions for describing correlation effects can be taken from atomic Hartree–Fock calculations i f the appropriate primitive set is used. Test calculations on oxygen‐containing molecules indicate that these primitive basis sets describe molecular correlation effects as well as the ANO sets of Almlof and Taylor. Guided by the calculations on oxygen, basis sets for use in correlated atomic and molecular calculations were developed for all of the first row atoms from boron through neon and for hydrogen. As in the oxygen atom calculations, it was found that the incremental energy lowerings due to the addition of correlating functions fall into distinct groups. This leads to the concept of c o r r e l a t i o n c o n s i s t e n t b a s i s s e t s, i.e., sets which include all functions in a given group as well as all functions in any higher groups. Correlation consistent sets are given for all of the atoms considered. The most accurate sets determined in this way, [5s4p3d2f1g], consistently yield 99% of the correlation energy obtained with the corresponding ANO sets, even though the latter contains 50% more primitive functions and twice as many primitive polarization functions. It is estimated that this set yields 94%–97% of the total (HF+1+2) correlation energy for the atoms neon through boron.
26,705 citations
TL;DR: The coupled cluster singles and doubles model (CCSD) as discussed by the authors is derived algebraically, presenting the full set of equations for a general reference function explicitly in spin-orbital form, and the computational implementation of the CCSD model, which involves cubic and quartic terms, is discussed and results are compared with full CI calculations for H2O and BeH2.
Abstract: The coupled‐cluster singles and doubles model (CCSD) is derived algebraically, presenting the full set of equations for a general reference function explicitly in spin–orbital form. The computational implementation of the CCSD model, which involves cubic and quartic terms, is discussed and results are reported and compared with full CI calculations for H2O and BeH2. We demonstrate that the CCSD exponential ansatz sums higher‐order correlation effects efficiently even for BeH2, near its transition state geometry where quasidegeneracy efforts are quite large, recovering 98% of the full CI correlation energy. For H2O, CCSD plus the fourth‐order triple excitation correction agrees with the full CI energy to 0.5 kcal/mol. Comparisons with low‐order models provide estimates of the effect of the higher‐order terms T1T2, T21T2, T31, and T41 on the correlation energy.
5,603 citations
TL;DR: An overview of NWChem is provided focusing primarily on the core theoretical modules provided by the code and their parallel performance, as well as Scalable parallel implementations and modular software design enable efficient utilization of current computational architectures.
4,666 citations
TL;DR: In this paper, a comprehensive overview of the equation of motion coupled-cluster (EOM•CC) method and its application to molecular systems is presented by exploiting the biorthogonal nature of the theory, it is shown that excited state properties and transition strengths can be evaluated via a generalized expectation value approach that incorporates both the bra and ket state wave functions.
Abstract: A comprehensive overview of the equation of motion coupled‐cluster (EOM‐CC) method and its application to molecular systems is presented. By exploiting the biorthogonal nature of the theory, it is shown that excited state properties and transition strengths can be evaluated via a generalized expectation value approach that incorporates both the bra and ket state wave functions. Reduced density matrices defined by this procedure are given by closed form expressions. For the root of the EOM‐CC effective Hamiltonian that corresponds to the ground state, the resulting equations are equivalent to the usual expressions for normal single‐reference CC density matrices. Thus, the method described in this paper provides a universal definition of coupled‐cluster density matrices, providing a link between EOM‐CC and traditional ground state CC theory.Excitation energy,oscillator strength, and property calculations are illustrated by means of several numerical examples, including comparisons with full configuration interaction calculations and a detailed study of the ten lowest electronically excited states of the cyclic isomer of C4.
2,171 citations