On the determination of partial differential cross sections for photodetachment and photoionization processes producing polyatomic molecules with electronic states coupled by conical intersections

A formalism is derived for the computation of partial differential cross sections for electron photodetachment and photoionization processes that leave the residual or target molecule in electronic states that are strongly coupled by conical intersections. Because the electronic states of the target are nonadiabatically coupled, the standard adiabatic states approach of solving the electronic Schro?dinger equation for the detached electron at fixed nuclear geometries and then vibrationally averaging must be fundamentally modified. We use a Lippmann-Schwinger equation based approach, which leads naturally to a partitioning of the transition amplitude into a Dyson orbital like part plus a scattering correction. The requisite Green's function is that developed in our previous paper for the direct determination of total integral cross sections. The method takes proper account of electron exchange, possible nonorthogonality of the orbital describing the detached electron, and nonadiabatic effects in the product molecule. The Green's function is constructed in an L(2) basis using complex scaling techniques. The accurate treatment of nonadiabatic effects in the residual molecule is accomplished using the multimode vibronic coupling model. For photodetachment, an approximate approach, which is less computationally demanding, is suggested.     

Virtual state scattering with cold electrons: para-xylene and para-difluorobenzene

The scattering of electrons with kinetic energies down to a few meV by para-xylene and para-difluorobenzene has been observed experimentally with an electron beam energy resolution of 0.95 to 1.5 meV (full width half maximum). At low electron energies the collisions can be considered as cold scattering events because the de Broglie wavelength of the electron is considerably larger than the target dimensions. The scattering cross sections measured rise rapidly at low energy due to virtual state scattering. The nature of this scattering process is discussed using s- and p-wave phase shifts derived from the experimental data. Scattering lengths are derived of, respectively, -9.5+/-0.5 and -8.0+/-0.5 a.u. for para-xylene and para-difluorobenzene. The virtual state effect is interpreted in terms of nuclear diabatic and partially adiabatic models, involving the electronic and vibronic symmetries of the unoccupied orbitals in the target species. The concept of direct and indirect virtual state scattering is introduced, through which the present species, in common with carbon dioxide and benzene, scatter through an indirect virtual state process, whereas other species, such as perfluorobenzene, scatter through a direct process.     

Isomorphic electron orbitals for vibronic flexibility in a cyclopropenyl radical molecular device

Summary The vibronic character of this molecular device has been studied using isomorphic electron orbitals. The leading role of the softest vibrational mode for the electron transport process is stressed by the quantum mechanical treatment of the rearrangement operator. The theory was used to investigate the possible function of the soliton valve, which has been suggested as a switching tip. The electronic flexibility of the cyclopropenyl radical with respect to molecular vibrations, which is important for the function of the molecular device, is well characterized by the hardness and softness of the electron structure in terms of the orbital energy-occupation number correlation diagram.     

An ab initio study of the vibronic, spin-orbit, and magnetic hyperfine structure in the X2Pi electronic state of NCO

In the present study we give the results of the ab initio calculations on the vibronic, spin-orbit, and magnetic hyperfine structure in the X (2)Pi electronic state of the NCO radical. The calculations of the potential surfaces and the electronic mean values of the hyperfine coupling constants are carried out by means of the density functional theory approach (B3LYP functional combined with an atomic orbital basis set suitable for calculations of the hyperfine structure). The vibronic levels, spin-orbit splitting, and the vibronic mean values of the components of the hyperfine tensor in the vibronic species are calculated using a variational method. The results of the calculations are in good agreement with the available experimental data.     

Møller-Plesset (MP2) perturbation theory for large molecules

Summary A novel formulation of MP2 theory is presented which starts from the Laplace transform MP2 ansatz, and subsequently moves from a molecular orbital (MO) representation to an atomic orbital (AO) representation. Consequently, the new formulation is denoted AO-MP2. As in traditional MP2 approaches electron repulsion integrals still need to be transformed. Strict bounds on the individual MP2 energy contribution of each intermediate four-index quantity allow to screen off numerically insignificant integrals with a single threshold parameter. Implicit in our formulation is a bound to two-particle density matrix elements. For small molecules the computational cost for AO-MP2 calculations is about a factor of 100 higher than for traditional MO-based approaches, but due to screening the computational effort in larger systems will only grow with the fourth power of the size of the system (or less) as is demonstrated both in theory and in application. MP2 calculations on (non-metallic) crystalline systems seem to be a feasible extension of the Laplace transform approach. In large molecules the AO-MP2 ansatz allows massively parallel MP2 calculations without input/output of four-index quantities provided that each processor has in-core memory for a limited number of two-index quantities. Energy gradient formulas for the AO-MP2 approach are derived.Dedicated to Prof. W. Kutzelnigg whose books on theoretical chemistry aroused my interest in this field     

Empirical correlation methods for temporary anions

A temporary anion is a short-lived radical anion that decays through electron autodetachment into a neutral molecule and a free electron. The energies of these metastable species are often predicted using empirical correlation methods because ab initio predictions are computationally very expensive. Empirical correlation methods can be justified in the framework of Weisskopf-Fano-Feshbach theory but tend to work well only within closely related families of molecules or within a restricted energy range. The reason for this behavior can be understood using an alternative theoretical justification in the framework of the Hazi-Taylor stabilization method, which suggests that the empirical parameters do not so much correct for the coupling of the computed state to the continuum but for electron correlation effects and that therefore empirical correlation methods can be improved by using more accurate electronic structure methods to compute the energy of the confined electron. This idea is tested by choosing a heterogeneous reference set of temporary states and comparing empirical correlation schemes based on Hartree-Fock orbital energies, Kohn-Sham orbital energies, and attachment energies computed with the equation-of-motion coupled-cluster method. The results show that using more reliable energies for the confined electron indeed enhances the predictive power of empirical correlation schemes and that useful correlations can be established beyond closely related families of molecules. Certain types of σ* states are still problematic, and the reasons for this behavior are analyzed. On the other hand, preliminary results suggest that the new scheme can even be useful for predicting energies of bound anions at a fraction of the computational cost of reliable ab initio calculations. It is then used to make predictions for bound and temporary states of the furantrione and croconic acid radical anions.     

Second order explicitly correlated R12 theory revisited: a second quantization framework for treatment of the operators' partitionings

Second order R12 theory is presented and derived alternatively using the second quantized hole-particle formalism. We have shown that in order to ensure the strong orthogonality between the R12 and the conventional part of the wave function, the explicit use of projection operators can be easily avoided by an appropriate partitioning of the involved operators to parts which are fully describable within the computational orbital basis and complementary parts that involve imaginary orbitals from the complete orbital basis. Various Hamiltonian splittings are discussed and computationally investigated for a set of nine molecules and their atomization energies. If no generalized Brillouin condition is assumed, with all relevant partitionings the one-particle contribution arising in the explicitly correlated part of the first order wave function has to be considered and has a significant role when smaller atomic orbital basis sets are used. The most appropriate Hamiltonian splitting results if one follows the conventional perturbation theory for a general non-Hartree-Fock reference. Then, no couplings between the R12 part and the conventional part arise within the first order wave function. The computationally most favorable splitting when the whole complementary part of the Hamiltonian is treated as a perturbation fails badly. These conclusions also apply to MP2-F12 approaches with different correlation factors.     

Vibronic coupling simulations for linear and nonlinear optical processes: theory

A comprehensive vibronic coupling model based on the time-dependent wavepacket approach is derived to simulate linear optical processes, such as one-photon absorbance and resonance Raman scattering, and nonlinear optical processes, such as two-photon absorbance and resonance hyper-Raman scattering. This approach is particularly well suited for combination with first-principles calculations. Expressions for the Franck-Condon terms, and non-Condon effects via the Herzberg-Teller coupling approach in the independent-mode displaced harmonic oscillator model are presented. The significance of each contribution to the different spectral types is discussed briefly.     

Nonadiabatic scattering of NO off Au3 clusters: A simple and robust diabatic state manifold generation method for multiconfigurational wavefunctions

The presence of nonadiabatic effects during the interaction of small molecules with metals has been observed experimentally for the last decades. Specially remarkable are the effects found for NO/Au, where experiments have suggested the presence of very strong vibronic coupling during the molecular scattering. However, the accurate inclusion of the nonadiabatic effects in periodic boundary conditions (PBC) theoretical methods remain an unapproachable challenge. Here, aiming to give some theoretical insight to the strong vibronic coupling, we have adopted a pragmatic point of view, taking use of an auxiliary simplified system, NO/Au₃. We show the importance of nonadiabatic coupling, during the scattering of NO from a Au₃ cluster, using a diabatic representation of 12 electronic states of the system, including a few charge-transfer states. Our diabatic representation is obtained by rotating the orbital and configuration interaction (CI) vectors of a restricted active space (RAS) wavefunction. We present a strategy for extracting the best effective manifold of states relevant to the system, below some prescribed energy, directly from the RAS CI vectors. This scheme is able to disentangle a large dense manifold of adiabatic states with strong coupling and crossings. This approach is also shown to work for multireference configuration interaction (MRCI). By performing quantum propagations, we observed an increase in vibrational redistribution with increasing initial vibrational or translational energies. We suggest that these nonadiabatic effects should also be present at smaller energies in larger clusters. © 2018 Wiley Periodicals, Inc.     

Independent particle theory with electron correlation

We formulate an effective independent particle model where the effective Hamiltonian is composed of the Fock operator and a correlation potential. Within the model the kinetic energy and the exchange energy can be expressed exactly leaving the correlation energy functional as the remaining unknown. Our efforts concentrate on finding a correlation potential such that exact ionization potentials and electron affinities can be reproduced as orbital energies. The equation-of-motion coupled-cluster approach enables us to define an effective Hamiltonian from which a correlation potential can be extracted. We also make the connection to electron propagator theory. The disadvantage of the latter is the inherit energy dependence of the potential resulting in a different Hamiltonian for each orbital. Alternatively, the Fock space coupled-cluster approach employs an effective Hamiltonian which is energy independent and universal for all orbitals. A correlation potential is extracted which yields the exact ionization potentials and electron affinities and a set of associated molecular orbitals. We also describe the close relationship to Brueckner theory.     

Negative electron binding energies observed in a triply charged anion: photoelectron spectroscopy of 1-hydroxy-3,6,8-pyrene-trisulfonate

We report the observation of negative electron binding energies (BEs) in a triply charged anion, 1-hydroxy-3,6,8-pyrene-trisulfonate (HPTS(3-)). Low-temperature photoelectron spectra were obtained for HPTS(3-) at several photon energies, revealing three detachment features below 0 electron BE. The HPTS(3-) trianion was measured to possess a negative BE of -0.66 eV. Despite the relatively high excess energy stored in HPTS(3-), it was observed to be a long-lived anion due to its high repulsive Coulomb barrier (RCB) ( approximately 3.3 eV), which prevents spontaneous electron emission. Theoretical calculations were carried out, which confirmed the negative electron BEs observed. The calculations further showed that the highest occupied molecular orbital in HPTS(3-) is an antibonding pi orbital on the pyrene rings, followed by lone pair electrons in the peripheral -SO(3) (-) groups. Negative electron BE is a unique feature of multiply charged anions due to the presence of the RCB. Such metastable species may be good models to study electron-electron and vibronic interactions in complex molecules.     

Exponential perturbation theories for inelastic scattering

An Exponential Perturbation Theory (EPT) is derived whereby one calculates a phase-shift matrix by an nth order perturbation theory and then exponentiates it to obtain the scattering matrix. The theory has been developed to include high-order terms, closed channels and resonances. The radial wavefunctions used are WKB solutions which are generalized to cases where there are multiple turning points. The orbital angular momentum may be treated exactly or in the classical or sudden limits. Calculations are done for the rotationally inelastic scattering in He + H₂, Ar + N₂ and Ar + HCl. The first two systems give fair to good agreement with accurate calculations; the last case gives poor agreement. The first-order EPT is very much better than the first-order distorted-wave approximation.     

Determining partial differential cross sections for low-energy electron photodetachment involving conical intersections using the solution of a Lippmann-Schwinger equation constructed with standard electronic structure techniques

A method for obtaining partial differential cross sections for low energy electron photodetachment in which the electronic states of the residual molecule are strongly coupled by conical intersections is reported. The method is based on the iterative solution to a Lippmann-Schwinger equation, using a zeroth order Hamiltonian consisting of the bound nonadiabatically coupled residual molecule and a free electron. The solution to the Lippmann-Schwinger equation involves only standard electronic structure techniques and a standard three-dimensional free particle Green's function quadrature for which fast techniques exist. The transition dipole moment for electron photodetachment, is a sum of matrix elements each involving one nonorthogonal orbital obtained from the solution to the Lippmann-Schwinger equation. An expression for the electron photodetachment transition dipole matrix element in terms of Dyson orbitals, which does not make the usual orthogonality assumptions, is derived.     

The estimation of electron affinities fromab initio 1s orbital energies

Summary It has been found that the electron affinities of alkoxy-radicals can be estimated using a correlation with the 1s orbital energy of the oxygen on the associated alkoxy-anion, EA=–0.64503 * (1s orbital energy) –351.58. The method assumes that the species of interest accepts the electron into an orbital which is localized on the oxygen.     

Reaction mechanism in the mechanochemical synthesis of dibenzophenazine: application of vibronic coupling density analysis

The reaction mechanism for mechanochemical synthesis of dibenzophenazine was theoretically investigated in terms of the vibronic coupling density, which describes the interactions between electrons and nuclear motions. The concept theoretically indicates experimentally observed reactive sites that cannot be explained by the conventional frontier orbital theory. The results of vibronic coupling density analysis suggested the difference between reaction mechanisms under thermal and mechanochemical conditions.     

Local approach to coupled cluster evaluation of polarizabilities for long conjugated molecules

Electron correlation effect in conjugated polymer is a long-standing problem, especially for the nonlinear optical properties. We have implemented a spin-adapted Coupled Cluster with singles and doubles excitation with local molecular orbital approach. As a first application, we evaluate the static polarizability of conjugated polyene chains with finite field approach. It is found that the local molecular orbital approach can tremendously reduce the computational costs at sufficiently high accuracy. It is also found that the electron correlation can largely reduce the molecular polarizability with respect to the Hartree-Fock mean field results.     

Vibronic coupling in benzene cation and anion: vibronic coupling and frontier electron density in Jahn-Teller molecules

Vibronic coupling constants of Jahn-Teller molecules, benzene radical cation and anion, are computed as matrix elements of the electronic part of the vibronic coupling operator using the electronic wave functions calculated by generalized restricted Hartree-Fock and state-averaged complete active space self-consistent-field methods. The calculated vibronic coupling constants for benzene cation agree well with the experimental and theoretical values. Vibronic coupling density analysis, which illustrates the local properties of the coupling, is performed in order to explain the order of magnitude of the coupling constant from view of the electronic and vibrational structures. This analysis reveals that the couplings of the e2g2 and e2g3 modes in which the large displacements locate on C-C bonds are strong in the cation. On the other hand, they are greatly weakened in the anion because of the decrease of electron density in the region of the C-C bonds, which originates from the antibonding nature of the singly occupied molecular orbital of the anion. However, the difference of the electronic structure has a little influence on the vibronic coupling of the e2g4 mode. These results indicate that the vibronic coupling depends not only on the direction of the nuclear displacement but also on the frontier electron density.     

Analysis of the Electronic Structure of Aqueous Urea and Its Derivatives: A Systematic Soft X‐Ray–TD‐DFT Approach

Soft X‐ray emission (XE), absorption (XA), and resonant inelastic scattering (RIXS) experiments have been conducted at the nitrogen K‐edge of urea and its derivatives in aqueous solution and were compared with density functional theory and time‐dependent density functional theory calculations. This comprehensive study provides detailed information on the occupied and unoccupied molecular orbitals of urea, thiourea, acetamide, dimethylurea, and biuret at valence levels. By identifying the electronic transitions that contribute to the experimental spectral features, the energy gap between the highest occupied and the lowest unoccupied molecular orbital of each molecule is determined. Moreover, a theoretical approach is introduced to simulate resonant inelastic X‐ray scattering spectra by adding an extra electron to the lowest unoccupied molecular orbital, thereby mimicking the real initial state of the core‐electron absorption before the subsequent relaxation process.