Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userTime-dependent complete-active-space self-consistent-field method for atoms: Application to high-harmonic generation
82
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present the numerical implementation of the time-dependent complete-active-space self-consistent-field (TD-CASSCF) method [Phys. Rev. A, 88, 023402 (2013)] for atoms driven by a strong linearly polarized laser pulse. The present implementation treats the problem in its full dimensionality and introduces a gauge-invariant frozen-core approximation, an efficient evaluation of the Coulomb mean field scaling linearly with the number of basis functions, and a split-operator method specifically designed for stable propagation of stiff spatial derivative operators. We apply this method to high-harmonic generation in helium, beryllium, and neon and explore the role of electron correlations.
rate research
Read More
Fragmentation methods applied to multireference wave functions constitute a road towards the application of highly accurate ab initio wave function calculations to large molecules and solids. However, it is important for reproducibility and transferability that a fragmentation scheme be well-defined with minimal dependence on initial orbital guesses or user-designed ad hoc fragmentation schemes. One way to improve this sort of robustness is to ensure the energy obeys a variational principle; i.e., that the active orbitals and active space wave functions minimize the electronic energy in a certain ansatz for the molecular wave function. We extended the theory of the localized active space self-consistent field, LASSCF, method (JCTC 2019, 15, 972) to fully minimize the energy with respect to all orbital rotations, rendering it truly variational. The new method, called vLASSCF, substantially improves the robustness and reproducibility of the LAS wave function compared to LASSCF. We analyze the storage and operation cost scaling of vLASSCF compared to orbital optimization using a standard CASSCF approach and we show results of vLASSCF calculations on some simple test systems. We show that vLASSCF is energetically equivalent to CASSCF in the limit of one active subspace, and that vLASSCF significantly improves upon the reliability of LASSCF energy differences, allowing for more meaningful and subtle analysis of potential energy curves of dissociating molecules. We also show that all forms of LASSCF have a lower operation cost scaling than the orbital-optimization part of CASSCF.
We report an efficient algorithm using density fitting for the relativistic complete active space self-consistent field (CASSCF) method, which is significantly more stable than the algorithm previously reported by one of the authors [J. E. Bates and T. Shiozaki, J. Chem. Phys. 142, 044112 (2015)]. Our algorithm is based on the second-order orbital update scheme with an iterative augmented Hessian procedure, in which the density-fitted orbital Hessian is directly contracted to the trial vectors. Using this scheme, each microiteration is made less time consuming than one Dirac-Hartree-Fock iteration, and macroiterations converge quadratically. In addition, we show that the CASSCF calculations with the Gaunt and full Breit interactions can be efficiently performed by means of approximate orbital Hessians computed with the Dirac-Coulomb Hamiltonian. It is demonstrated that our algorithm can also be applied to systems under an external magnetic field, for which all of the molecular integrals are computed using gauge-including atomic orbitals.
The full-dimensional time-dependent Schrodinger equation for the electronic dynamics of single-electron systems in intense external fields is solved directly using a discrete method. Our approach combines the finite-difference and Lagrange mesh methods. The method is applied to calculate the quasienergies and ionization probabilities of atomic and molecular systems in intense static and dynamic electric fields. The gauge invariance and accuracy of the method is established. Applications to multiphoton ionization of positronium and hydrogen atoms and molecules are presented. At very high intensity above saturation threshold, we extend the method using a scaling technique to estimate the quasienergies of metastable states of the hydrogen molecular ion. The results are in good agreement with recent experiments.
We present an implementation of a time-dependent multiconfiguration self-consistent-field (TD-MCSCF) method [R. Anzaki et al., Phys. Chem. Chem. Phys. 19, 22008 (2017)] with the full configuration interaction expansion for coupled electron-nuclear dynamics in diatomic molecules subject to a strong laser field. In this method, the total wave function is expressed as a superposition of different configurations constructed from time-dependent electronic Slater determinants and time-dependent orthonormal nuclear basis functions. The primitive basis functions of nuclei and electrons are strictly independent of each other without invoking the Born-Oppenheimer approximation. Our implementation treats the electronic motion in its full dimensionality on curvilinear coordinates, while the nuclear wave function is propagated on a one-dimensional stretching coordinate with rotational nuclear motion neglected. We apply the present implementation to high-harmonic generation and dissociative ionization of a hydrogen molecule and discuss the role of electron-nuclear correlation.
Bichromatic extreme-ultraviolet pulses from a seeded free-electron laser enable us to measure photoelectron angular distribution (PAD) as a function of the relative phase between the different wavelength components. The time-dependent multiconfiguration self-consistent-field (TD-MCSCF) methods are powerful multielectron computation methods to accurately simulate such photoionization dynamics from the first principles. Here we propose a method to evaluate the amplitude and phase of each ionization path, which completely determines the photoionization processes, using TD-MCSCF simulation results. The idea is to exploit the capability of TD-MCSCF to calculate the partial wave amplitudes specified by the azimuthal and magnetic angular momenta (l,m) and the m-resolved PAD. The phases of the ionization paths as well as the amplitudes of the paths resulting in the same (l,m) are obtained through global fitting of the expression of the asymmetry parameters to the calculated m-resolved PAD, which depends on the relative phase of the bichromatic field. We apply the present method to ionization of Ne by combined fundamental and second-harmonic XUV pulses, demonstrating that the extracted amplitudes and phases excellently reproduce the asymmetry parameters.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
