scispace - formally typeset
Search or ask a question
Journal ArticleDOI

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
Journal ArticleDOI
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

Journal ArticleDOI
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

Journal ArticleDOI
20 Jun 2016-ACS Nano
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

Journal ArticleDOI
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

Journal ArticleDOI
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
Journal ArticleDOI
TL;DR: These methods, capable of handling electron attachment to open-shell molecules, have been applied to the electron affinities of NH and C2, the excitation energies of CH, and the spectroscopic constants of all these molecules with the errors due to basis sets of finite sizes removed by extrapolation.
Abstract: High-order equation-of-motion coupled-cluster methods for electron attachment (EA-EOM-CC) have been implemented with the aid of the symbolic algebra program TCE into parallel computer programs Two types of size-extensive truncation have been applied to the electron-attachment and cluster excitation operators: (1) the electron-attachment operator truncated after the 2p-1h, 3p-2h, or 4p-3h level in combination with the cluster excitation operator after doubles, triples, or quadruples, respectively, defining EA-EOM-CCSD, EA-EOM-CCSDT, or EA-EOM-CCSDTQ; (2) the combination of up to the 3p-2h electron-attachment operator and up to the double cluster excitation operator [EA-EOM-CCSD(3p-2h)] or up to 4p-3h and triples [EA-EOM-CCSDT(4p-3h)] These methods, capable of handling electron attachment to open-shell molecules, have been applied to the electron affinities of NH and C2, the excitation energies of CH, and the spectroscopic constants of all these molecules with the errors due to basis sets of finite sizes removed by extrapolation The differences in the electron affinities or excitation energies between EA-EOM-CCSD and experiment are frequently in excess of 2 eV for these molecules, which have severe multideterminant wave functions Including higher-order operators, the EA-EOM-CC methods predict these quantities accurate to within 001 eV of experimental values In particular, the 3p-2h electron-attachment and triple cluster excitation operators are significant for achieving this accuracy

50 citations

Journal ArticleDOI
TL;DR: It is demonstrated that reliable excitation energies can be computed by EOM-DIP using a stabilized resonance wave function instead of the lowest energy solution corresponding to the neutral + free electron(s) state of the system, and benchmark calculations of the singlet-triplet energy gaps in several diradicals show a remarkable improvement of the E OM-Dip accuracy in problematic cases.
Abstract: The charge-stabilization method is applied to double ionization potential equation-of-motion (EOM-DIP) calculations to stabilize unstable dianion reference functions. The auto-ionizing character of the dianionic reference states spoils the numeric performance of EOM-DIP limiting applications of this method. We demonstrate that reliable excitation energies can be computed by EOM-DIP using a stabilized resonance wave function instead of the lowest energy solution corresponding to the neutral + free electron(s) state of the system. The details of charge-stabilization procedure are discussed and illustrated by examples. The choice of optimal stabilizing Coulomb potential, which is strong enough to stabilize the dianion reference, yet, minimally perturbs the target states of the neutral, is the crux of the approach. Two algorithms of choosing optimal parameters of the stabilization potential are presented. One is based on the orbital energies, and another – on the basis set dependence of the total Hartree-Fock energy of the reference. Our benchmark calculations of the singlet-triplet energy gaps in several diradicals show a remarkable improvement of the EOM-DIP accuracy in problematic cases. Overall, the excitation energies in diradicals computed using the stabilized EOM-DIP are within 0.2 eV from the reference EOM spin-flip values.

48 citations

Journal ArticleDOI
TL;DR: The current work generates the potential energy curves for the ground and multiple excited states of the Li2 and Na2 molecules and calculates spectroscopic parameters of the systems studied.
Abstract: The recently developed method [M. Musial, J. Chem. Phys. 136, 134111 (2012)] to study double electron attached states has been applied to the description of the ground and excited state potential energy curves of the alkali metal dimers. The method is based on the multireference coupled cluster scheme formulated within the Fock space formalism for the (2,0) sector. Due to the use of the efficient intermediate Hamiltonian formulation, the approach is free from the intruder states problem. The description of the neutral alkali metal dimers is accomplished via attaching two electrons to the corresponding doubly ionized system. This way is particularly advantageous when a closed shell molecule dissociates into open shell subunits while its doubly positive cation generates the closed shell fragments. In the current work, we generate the potential energy curves for the ground and multiple excited states of the Li2 and Na2 molecules. In all cases the potential energy curves are smooth for the entire range of interatomic distances (from the equilibrium point to the dissociation limit). Based on the calculated potential energy curves, we are able to compute spectroscopic parameters of the systems studied.

42 citations

Journal ArticleDOI
TL;DR: A novel parallel algorithm for noniterative multireference coupled cluster (MRCC) theories, which merges recently introduced reference-level parallelism (RLP) with the possibility of accelerating numerical calculations using graphics processing units (GPUs) is presented.
Abstract: A novel parallel algorithm for noniterative multireference coupled cluster (MRCC) theories, which merges recently introduced reference-level parallelism (RLP) [Bhaskaran-Nair, K.; Brabec, J.; Apra, E.; van Dam, H. J. J.; Pittner, J.; Kowalski, K. J. Chem. Phys.2012, 137, 094112] with the possibility of accelerating numerical calculations using graphics processing units (GPUs) is presented. We discuss the performance of this approach applied to the MRCCSD(T) method (iterative singles and doubles and perturbative triples), where the corrections due to triples are added to the diagonal elements of the MRCCSD effective Hamiltonian matrix. The performance of the combined RLP/GPU algorithm is illustrated on the example of the Brillouin-Wigner (BW) and Mukherjee (Mk) state-specific MRCCSD(T) formulations.

39 citations