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: The diagonalization manifold in similarity transformed equation-of-motion coupled cluster (STEOM-CC) theory is extended to include doubly excited determinants in this article, and the resulting extended STEOM approach accurate results are obtained for doubly-excited states in small model systems for which full configuration interaction (CI) benchmark results are available (∼0.1 eV errors).
Abstract: The diagonalization manifold in similarity transformed equation-of-motion coupled cluster (STEOM-CC) theory is extended to include doubly excited determinants. In the resulting extended-STEOM approach accurate results are obtained for doubly excited states in small model systems for which full configuration interaction (CI) benchmark results are available (∼0.1 eV errors). On the other hand, extended-STEOM results are found to be virtually identical (<0.1 eV shifts) to the original STEOM results for states that are dominated by single excitations, at least in prototypical organic molecules. The extended-STEOM method is also applied to the transition metal complexes TiCl4, Ni(CO)4, and MnO4−, and yields improved results compared to STEOM and EOM-CCSD. For highly correlated systems, like the permangenate anion, results are not yet fully satisfactory however. In these cases the dominant source of error appears to be the description of ground, ionized, and attached states that underly the similarity transform...
108 citations
TL;DR: The application of these methods to Koopmans and satellite ionizations of N2 and CO has shown that IP-EOM-CCSD frequently accounts for orbital relaxation inadequately and displays errors exceeding a couple of eV, but these errors can be systematically reduced to tenths or even hundredths of an eV by IP-CCSDTT or CCSDTQ.
Abstract: Compact algebraic equations defining the equation-of-motion coupled-cluster (EOM-CC) methods for ionization potentials (IP-EOM-CC) have been derived and computer implemented by virtue of a symbolic algebra system largely automating these processes. Models with connected cluster excitation operators truncated after double, triple, or quadruple level and with linear ionization operators truncated after two-hole-one-particle (2h1p), three-hole-two-particle (3h2p), or four-hole-three-particle (4h3p) level (abbreviated as IP-EOM-CCSD, CCSDT, and CCSDTQ, respectively) have been realized into parallel algorithms taking advantage of spin, spatial, and permutation symmetries with optimal size dependence of the computational costs. They are based on spin-orbital formalisms and can describe both α and β ionizations from open-shell (doublet, triplet, etc.) reference states into ionized states with various spin magnetic quantum numbers. The application of these methods to Koopmans and satellite ionizations of N2 and C...
106 citations
TL;DR: In this paper, the gas-phase electron attachment free energies for polycyclic aromatic hydrocarbons with the electron-transfer equilibria technique were measured with electron transmission spectroscopy.
Abstract: The gas-phase electron attachment free energies ΔG a o , where-ΔGa a o ≃EA, were measured for several polycyclic aromatic hydrocarbons with the electron-transfer equilibria technique. The EA values in kcal/mol are as follow: anthracene, 13.2; benz[a]pyrene, 18.2; perylene, 22.3; tetracene, 23.9; and pentacene, 31.2. These results, when combined with the literature value for naphthalene, -4.4 kcal/mol, based on electron transmission spectroscopy, lead to an approximately linear correlation with the one-electron (polarographic) reduction potentials determined in solution. However, the slope of the line is not equal to 1 as previously assumed
105 citations
TL;DR: In this article, the generalized Hellmann-Feynman theorem was used to derive energy derivatives for nonvariational (e.g., coupled-cluster) methods invoking the generalized HELLMAN-FEYNman theorem, where one can construct a functional which, besides the usual wave-function parameters, contains new ones.
Abstract: It is possible to derive energy derivatives for nonvariational (e.g., coupled-cluster) methods invoking the generalized Hellmann–Feynman theorem. In such a procedure, one constructs a functional which, besides the usual wave-function parameters, contains new ones. One set of stationary conditions will reproduce exactly the original equations of the method, while the others will determine the value of the new parameters. We applied this straightforward procedure to derive analytic energy derivatives for several coupled-cluster (CC) methods applicable to excited states such as the Hilbert-space CC method, two-determinetal (TD) CC method, Fock-space CC method, and equation-of-motion–CC (EOM–CC) method. Finally, we compared the computational requirements for the different methods. © 1995 John Wiley & Sons, Inc.
101 citations
TL;DR: A benchmark theoretical determination of the electron affinities of benzene and linear oligoacenes ranging from naphthalene to hexacene is presented, using the principles of a focal point analysis.
Abstract: A benchmark theoretical determination of the electron affinities of benzene and linear oligoacenes ranging from naphthalene to hexacene is presented, using the principles of a focal point analysis. These energy differences have been obtained from a series of single-point calculations at the Hartree-Fock, second-, third-, and partial fourth-order Moller-Plesset (MP2, MP3, and MP4SDQ) levels and from coupled cluster calculations including single and double excitations (CCSD) as well as perturbative estimates of connected triple excitations [CCSD(T)], using basis sets of improving quality, containing up to 1386, 1350, 1824, 1992, 1630, and 1910 basis functions in the computations, respectively. Studies of the convergence properties of these energy differences as a function of the size of the basis set and order attained in electronic correlation enable a determination of the vertical electron affinities of the four larger terms of the oligoacene (C(2+4n)H(2+2n)) series within chemical accuracy (0.04 eV). According to our best estimates, these amount to +0.28, +0.82, +1.21, and +1.47 eV when n=3, 4, 5, and 6. Adiabatic electron affinities have been further calculated by incorporating corrections for zero-point vibrational energies and for geometrical relaxations. The same procedure was applied to determine the vertical electron affinities of benzene and naphthalene, which are found to be markedly negative ( approximately -1.53 and approximately -0.48 eV, respectively). Highly quantitative insights into experiments employing electron transmission spectroscopy on these compounds were also amenable from such an approach, provided diffuse atomic functions are deliberately removed from the basis set, in order to enforce confinement in the molecular region and enable a determination of pseudoadiabatic electron affinities (with respect to the timescale of nuclear motions). Comparison was made with calculations employing density functional theory and especially designed models that exploit the integer discontinuity in the potential or incorporate a potential wall in the unrestricted Kohn-Sham orbital equation for the anion.
97 citations