Construction of complex STO-NG basis sets by the method of least squares and their applications

Time-dependentab initioapproaches for high-harmonic generation spectroscopy

High-harmonic generation (HHG) is a nonlinear physical process used for the production of ultrashort pulses in XUV region, which are then used for investigating ultrafast phenomena in time-resolved spectroscopies. Moreover, HHG signal itself encodes information on electronic structure and dynamics of the target, possibly coupled to the nuclear degrees of freedom. Investigating HHG signal leads to HHG spectroscopy, which is applied to atoms, molecules, solids and recently also to liquids. Analysing the number of generated harmonics, their intensity and shape gives a detailed insight of, e.g., ionisation and recombination channels occurring in the strong-field dynamics. A number of valuable theoretical models has been developed over the years to explain and interpret HHG features, with the three-step model being the most known one. Originally, these models neglect the complexity of the propagating electronic wavefunction, by only using an approximated formulation of ground and continuum states. Many effects unravelled by HHG spectroscopy are instead due to electron correlation effects, quantum interference, and Rydberg-state contributions, which are all properly captured by anab initioelectronic-structure approach. In this review we have collected recent advances in modelling HHG by means ofab initiotime-dependent approaches relying on the propagation of the time-dependent Schrödinger equation (or derived equations) in presence of a very intense electromagnetic field. We limit ourselves to gas-phase atomic and molecular targets, and to solids. We focus on the various levels of theory employed for describing the electronic structure of the target, coupled with strong-field dynamics and ionisation approaches, and on the basis used to represent electronic states. Selected applications and perspectives for future developments are also given.

Fitting continuum wavefunctions with complex Gaussians: Computation of ionization cross sections

We implement a full nonlinear optimization method to fit continuum states with complex Gaussians. The application to a set of regular scattering Coulomb functions allows us to validate the numerical feasibility, to explore the range of convergence of the approach, and to demonstrate the relative superiority of complex over real Gaussian expansions. We then consider the photoionization of atomic hydrogen, and ionization by electron impact in the first Born approximation, for which the closed form cross sections serve as a solid benchmark. Using the proposed complex Gaussian representation of the continuum combined with a real Gaussian expansion for the initial bound state, all necessary matrix elements within a partial wave approach become analytical. The successful numerical comparison illustrates that the proposed all-Gaussian approach works efficiently for ionization processes of one-center targets.

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.

Optimization of complex slater-type functions with analytic derivative methods for describing photoionization differential cross sections

The complex basis function (CBF) method applied to various atomic and molecular photoionization problems can be interpreted as an L2 method to solve the driven-type (inhomogeneous) Schrödinger equation, whose driven term being dipole operator times the initial state wave function. However, efficient basis functions for representing the solution have not fully been studied. Moreover, the relation between their solution and that of the ordinary Schrödinger equation has been unclear. For these reasons, most previous applications have been limited to total cross sections. To examine the applicability of the CBF method to differential cross sections and asymmetry parameters, we show that the complex valued solution to the driven-type Schrödinger equation can be variationally obtained by optimizing the complex trial functions for the frequency dependent polarizability. In the test calculations made for the hydrogen photoionization problem with five or six complex Slater-type orbitals (cSTOs), their complex valued expansion coefficients and the orbital exponents have been optimized with the analytic derivative method. Both the real and imaginary parts of the solution have been obtained accurately in a wide region covering typical molecular regions. Their phase shifts and asymmetry parameters are successfully obtained by extrapolating the CBF solution from the inner matching region to the asymptotic region using WKB method. The distribution of the optimized orbital exponents in the complex plane is explained based on the close connection between the CBF method and the driven-type equation method. The obtained information is essential to constructing the appropriate basis sets in future molecular applications. © 2017 Wiley Periodicals, Inc.