Complex Absorbing Potentials with Voronoi Isosurfaces Wrapping Perfectly around Molecules
- 2 September 2015
- journal article
- research article
- Published by American Chemical Society (ACS) in Journal of Chemical Theory and Computation
- Vol. 11 (10) , 4627-4633
Abstract
Complex absorbing potentials (CAPs) are imaginary potentials that are added to a Hamiltonian to change the boundary conditions of the problem from scattering to square-integrable. In other words, with a CAP, standard bound-state methods can be used in problems involving unbound states such as identifying resonance states and predicting their energies and lifetimes. Although in wave packet dynamics, many CAP forms are used, in electronic structure theory, the so-called box-CAP is used almost exclusively, because of the ease of evaluating its integrals in a Gaussian basis set. However, the box-CAP does has certain disadvantages. First, it will, e.g., break the symmetry of Cnv molecules if n is odd and the main axis is placed along the z-axis by the “standard orientation” of the electronic structure code. Second, it provides a CAP starting at the smallest box around the entire molecular system. For larger molecules or clusters, which do not fill the space efficiently, that implies that much “dead space” within the molecule will be left, where there is neither a CAP nor a sufficient description with basis functions. Here, two new CAP forms are introduced and systematically explored: first, a Voronoi-CAP (that is, a CAP defined in each atom’s Voronoi cell), and second, a smooth Voronoi-CAP (which is similar to the Voronoi-CAP; however, the noncontinuously differentiable behavior at the surfaces between the Voronoi cells is smoothed out). Both have isosurfaces that are similar to the cavities used in solvation modeling. An obvious disadvantage of these two CAPs is that the integrals cannot be obtained analytically, but must be computed numerically. However, Voronoi-CAPs share the advantage of having the same symmetry as the molecular system, and, more importantly, considerably facilitate the treatment of larger molecules with asymmetric side chains and of molecular clusters.Keywords
Funding Information
- American Chemical Society Petroleum Research Fund
This publication has 28 references indexed in Scilit:
- Geometry-dependent lifetime of Interatomic coulombic decay using equation-of-motion coupled cluster methodThe Journal of Chemical Physics, 2014
- Complex absorbing potentials within EOM-CC family of methods: Theory, implementation, and benchmarksThe Journal of Chemical Physics, 2014
- Coupled cluster methods for autoionisation resonancesInternational Reviews in Physical Chemistry, 2014
- Strong field ionization rates simulated with time-dependent configuration interaction and an absorbing potentialThe Journal of Chemical Physics, 2014
- Equation-of-motion coupled-cluster method for the study of shape resonanceThe Journal of Chemical Physics, 2012
- CAP/SAC-CI method for calculating resonance states of metastable anionsChemical Physics Letters, 2012
- Complex absorbing potentials in the framework of electron propagator theory. II. Application to temporary anionsThe Journal of Chemical Physics, 2003
- Non-Hermitian electronic theory and applications to clustersPhysics Reports, 2002
- Ab initio calculation of energies and lifetimes of metastable dianions: The C22− resonanceThe Journal of Chemical Physics, 2000
- Temporary anions - calculation of energy and lifetime by absorbing potentials: the resonanceJournal of Physics B: Atomic, Molecular and Optical Physics, 1998