Variational approach to continuum orbitals in a spline basis: An application to H2+ photoionization

Continuum Electronic States: The Tiresia Code

A multicenter (LCAO) B-spline basis is described in detail, and its capabilities concerning affording convergent solutions for electronic continuum states and wavepacket propagation are presented. It forms the core of the Tiresia code, which implements static-DFT and TDDFT hamiltonians, as well as single channel Dyson-DFT and Dyson-TDDFT descriptions to include correlation in the bound states. Together they afford accurate and computationally efficient descriptions of photoionization properties of complex systems, both in the single photon and strong field environments. A number of examples are provided.

Trends in angle-resolved molecular photoelectron spectroscopy

In this perspective article, main trends of angle-resolved molecular photoelectron spectroscopy in the laboratory up to the molecular frame, in different regimes of light-matter interactions, are highlighted with emphasis on foundations and most recent applications.

Calculation of photoionization differential cross sections using complex Gauss-type orbitals

Accurate theoretical calculation of photoelectron angular distributions for general molecules is becoming an important tool to image various chemical reactions in real time. We show in this article that not only photoionization total cross sections but also photoelectron angular distributions can be accurately calculated using complex Gauss-type orbital (cGTO) basis functions. Our method can be easily combined with existing quantum chemistry techniques including electron correlation effects, and applied to various molecules. The so-called two-potential formula is applied to represent the transition dipole moment from an initial bound state to a final continuum state in the molecular coordinate frame. The two required continuum functions, the zeroth-order final continuum state and the first-order wave function induced by the photon field, have been variationally obtained using the complex basis function method with a mixture of appropriate cGTOs and conventional real Gauss-type orbitals (GTOs) to represent the continuum orbitals as well as the remaining bound orbitals. The complex orbital exponents of the cGTOs are optimized by fitting to the outgoing Coulomb functions. The efficiency of the current method is demonstrated through the calculations of the asymmetry parameters and molecular-frame photoelectron angular distributions of H2+ and H2 . In the calculations of H2 , the static exchange and random phase approximations are employed, and the dependence of the results on the basis functions is discussed. © 2017 Wiley Periodicals, Inc.

Attosecond electronic and nuclear quantum photodynamics of ozone monitored with time and angle resolved photoelectron spectra

Recently we reported a series of numerical simulations proving that it is possible in principle to create an electronic wave packet and subsequent electronic motion in a neutral molecule photoexcited by a UV pump pulse within a few femtoseconds. We considered the ozone molecule: for this system the electronic wave packet leads to a dissociation process. In the present work, we investigate more specifically the time-resolved photoelectron angular distribution of the ozone molecule that provides a much more detailed description of the evolution of the electronic wave packet. We thus show that this experimental technique should be able to give access to observing in real time the creation of an electronic wave packet in a neutral molecule and its impact on a chemical process.

Time-dependent quantum simulation of coronene photoemission spectra

Photoelectron spectroscopy is usually described by a simple equation that relates the binding energy of the photoemitted electron, Ebinding, its kinetic energy, Ekinetic, the energy of the ionizing photon, Ephoton, and the work function of the spectrometer, ϕ, Ebinding = Ephoton - Ekinetic - ϕ. Behind this equation there is an extremely rich physics, which we describe here using as an example a relatively simple conjugated molecule, namely coronene. The theoretical analysis of valence band and C1s core level photoemission spectra showed that multiple excitations play an important role in determining the intensities of the final spectrum. An explicit, time-evolving model is applied, which is able to count all possible photo-excitations occurring during the photoemission process, showing that they evolve on a short time-scale, of about 10 fs. The method reveals itself to be a valid approach to reproduce photoemission spectra of polycyclic aromatic hydrocarbons (PAHs).

Time-dependent Dyson orbital theory

Probabilities of single ionization channels or double and multiple simultaneous ionization can be determined from exact master equations for time-dependent Dyson orbitals.

Conformational Effects on Circular Dichroism in the Photoelectron Angular Distribution

The B-spline density-functional method has been applied to the conformers of the (1R, 2R)-1,2-dibromo-1,2-dichloro-1,2-difluoroethane molecule. The cross section, asymmetry, and dichroic parameters relative to core and valence orbitals, which do not change their nature along the conformational curve, have been systematically studied. While the cross section and the asymmetry parameter are weakly affected, the dichroic parameter appears to be rather sensitive to the particular conformer of the molecule, suggesting that this dynamical property could be a useful tool for conformational analysis. The computational method has also been applied to methyl rotation in methyloxirane. Unexpected and dramatic sensitivity of the dichroic-parameter profile to the methyl rotation, both in the core and valence states, has been found. Boltzmann averaging over the conformers reproduces quite closely the profiles previously obtained for the minimum-energy conformation, which is in good agreement with the experimental results.

Simple method for obtaining electron scattering phase shifts from energies of an atom in a cavity

We present a simple method for obtaining elastic scattering phase shifts and cross sections from precise ab initio many-body perturbation theory energies of atoms in variable cavities. This method does not require calculations of wave functions of continuum states, can be generalized to many atoms and ions, and is extremely convenient because existing codes developed for energy calculations can be used without modification. The high precision of the method and close agreement with experiment are illustrated on examples of e-Ar and e-Kr scattering. Correlations as well as relativistic corrections are systematically considered.

Theoretical study on the circular dichroism in core and valence photoelectron angular distributions of camphor enantiomers

In the present work the photoelectron circular dichroism of camphor has been theoretically studied using B-spline and continuum multiple scattering-Xalpha methods, and comparisons are made with available experimental data. In general, rather large dichroism effects have been found for both valence and core (O 1s, C 1s) photoionizations. The agreement between the two calculations reported here and previous experimental measurements for core C 1s data is essentially quantitative. For valence ionization satisfactory agreement between theory and experiment has been obtained and the discrepancies have been attributed to both exchange-correlation potential limitations and the absence of response effects in the adopted formalism. The calculations predict, moreover, important features in the cross-section profiles, which have been discussed in terms of dipole-prepared continuum orbitals.

Branching ratio deviations from statistical behavior in core photoionization

Accurate calculations of carbon 1s photoionization cross sections have been performed at the density functional level with the B-spline linear combination of atomic orbitals approach. The molecules considered are FC[triple bond]CH, FC[triple bond]CCH3, FC[triple bond]CCN, F2C=CH2, CF3COOCH2CH3, and C3H5O. The variation of the branching ratios relative to inequivalent C 1s ionizations have been evaluated from threshold to about 100 eV photoelectron kinetic energy. Large deviations from the statistical ratios are observed at low energies, which remain often significant several tens of eV above threshold. The importance of taking into account core branching ratios for peak deconvolution and quantitative analysis, as well as an additional tool for structural information, is pointed out. Strong shape resonant effects are found to largely cancel in branching ratios. Their nature and variation along the series is analyzed in the framework of excitations into sigma* valence orbitals.

Time-dependent density-functional theory for molecular photoionization with noniterative algorithm and multicenter B-spline basis set: CS2 and C6H6 case studies

In this work a new direct (noniterative) algorithm to solve the time-dependent density-functional theory equations for molecular photoionization has been proposed and implemented, using a multicentric basis set expansion of B-spline functions and complete exploiting of the molecular point-group symmetry. The method has been applied to study the photoionization dynamics of CS2 and C6H6: the results confirmed the expectation of large screening effects in CS2. For C6H6 the screening effects have been found to play a minor role than in CS2, however, also in this case the quality of the final results is definitely improved. The method has proven suitable to study with confidence molecules of medium size, and there is still room for further improvement working on more elaborate treatment of the exchange-correlation functional.

Density functional study on the circular dichroism of photoelectron angular distribution from chiral derivatives of oxirane

The linear combination of atomic orbitals B-spline density functional method has been successfully applied to a series of four chiral derivatives of oxirane, to calculate the photoionization dynamical parameters, the circular dichroism in the angular distribution effect, and to identify trends along the series. The computational algorithm has proven numerically stable and computationally competitive. The photoionization cross section, asymmetry, and dichroic parameter profiles relative to valence orbitals have been systematically studied for the states which retain their nature along the series: the identified trends have been ascribed to the different electronic properties of the substituents. A rather unexpected sensitivity of the dichroic parameter to changes in the electronic structure has been found in many instances, making this dynamical property suitable to investigate the electronic structure of chiral compounds. The magnitude of the circular dichroism in the angular distribution effect does not seem to be associated with the initial state chirality, but rather to be governed by the ability of the delocalized photoelectron wave function to probe the asymmetry of the molecular effective potential.