ezDyson

NO dimer


ezDyson is a C++ code that calculates absolute photodetachment/ photoionization cross sections, photoelectron angular distributions (PADs), and anisotropy parameters (beta) using Dyson orbitals computed by an ab initio program. The calculation is based on the following approximations: (1) weak field limit, i.e., the photodetachment or photoionization is treated in a perturbative regime; (2) dipole approximation, i.e., assuming that the wavelength of the radiation field is longer than the size of the molecule and that the process is one-photon (first order in perturbation theory); (3) sudden approximation, i.e., we ignore the interactions between the ejected photoelectron and the remaining core electrons; (4) strong orthogonality condition, i.e., the continuum orbital is orthogonal to all states of initial system. We use plane waves or Coulomb waves to describe the continuum state of the photoelectron.



Dyson orbitals are defined as the overlap between an N electron molecular wavefunction and the N-1 electron wave function of the ionized system:

eq1

They contain all the information about the ionized system required for calculating the differential cross-section:

eq2

where DkIF is the photoelectron dipole matrix element connecting the initial and the final wave functions:

eq3

The ejected photoelectron is described by a plane wave or Coulomb wave.

The code requires information from a quantum chemical calculation (geometry, atomic basis set), details of the experiment (laser polarization, ionization energy, molecular orientation, range of electron kinetic energies), the coefficients of the Dyson orbital in the AO basis, and the norm of the Dyson orbital. The ezDyson input is in xml format. It can be prepared manually, or a Python interface can be used to extract data from an electronic structure program. A Python interface is currently provided only for the Q-Chem outputs. Dyson orbital calculations are implemented within the equation-of-motion coupled-cluster (EOM-CC) suite of methods in Q-Chem. The details of the theoretical model employed as well as benchmark results can be found in the following publications:

[1] C.M. Oana and A.I. Krylov, J. Chem. Phys. 127, 234106 (2007)
[2] C.M. Oana and A.I. Krylov, J. Chem. Phys. 131, 124114 (2009)
[3] S. Gozem, A.O. Gunina, T. Ichino, D.L. Osborn, J.F. Stanton, and A.I. Krylov, J. Phys. Chem. Lett. 6, 4532-4540 (2015)

Please refer to the manual for detailed theoretical background, instructions on running ezDyson, an overview of recent updates, and a description of ezDyson input and output files.

Download ezDyson: ezDyson_v5.0.tar.gz

Download manual: ezDysonManual.pdf

Watch webinar describing ezSpectra suite by Professor Samer Gozem

To install, execute:
tar -xzf ezDyson_v5.0.tar.gz

Precompiled Linux and Mac OS executables are provided. However, the code can be easily recompiled for other platforms by using an appropriate makefile.

The original version of the code was developed by Dr. Melania Oana. Dr. Evgeny Epifanovsky and Dr. Ilya Kaliman have helped to parallelize the code. Dr. Kadir Diri and Dr. Anastasia Gunina contributed to developing the manual. Sahil Gulania wrote a section of how to plot Dyson orbitals. The ccman2 implementation of Dyson orbitals was developed by Anastasia Gunina and Tomek Kuš, with the help of Dr. Evgeny Epifanovsky. Dr. Takatoshi Ichino and Prof. John Stanton helped to validate the implementation. This work has been funded by the National Science Foundation and by the Department of Energy.

To acknowledge ezDyson use the following citation:
S. Gozem, A. I. Krylov, "The ezSpectra suite: An easy-to-use toolkit for spectroscopy modeling", WIRES CMS, e1546 (2021).

and

S. Gozem, A.O. Gunina, T. Ichino, D.L. Osborn, J.F. Stanton, and A.I. Krylov, J. Phys. Chem. Lett. 6, 4532-4540 (2015).