Search

The purpose of this document is to describe the solution and implementation of the time-independent and time-dependent Schr\"odinger using pseudospectral methods. Currently, the description is for single particle systems interacting with a classical electromagnetic field in spherical coordinates.

We present a numerical approach to magnetic optical rotation based on real-time time-dependent electronic-structure theory. Not relying on perturbation expansions in the magnetic-field strength, the formulation allows us to test the range of validity of the linear relation between the rotation angle per unit path length and the magnetic-field strength that was established empirically by Verdet 160 years ago. Results obtained from time-dependent coupled-cluster and time-dependent current density-functional theory are presented for the closed-shell molecules H2, HF, and CO in magnetic fields up to 55 kT at standard temperature and pressure conditions. We find that Verdet's linearity remains valid up to roughly 10-20 kT, above which significant deviations from linearity are observed. Among the three current density-functional approximations tested in this work, the current-dependent Tao-Perdew-Staroverov-Scuseria hybrid functional performs the best in comparison with time-dependent coupled-cluster singles and doubles results for the magnetic optical rotation.

Based on the work done by an electromagnetic field on an atomic or molecular electronic system, a general gauge invariant formulation of transient absorption spectroscopy is presented within the semi-classical approximation. Avoiding multipole expansions, a computationally viable expression for the spectral response function is derived from the minimal-coupling Hamiltonian of an electronic system interacting with one or more laser pulses described by a source-free, enveloped electromagnetic vector potential. With a fixed-basis expansion of the electronic wave function, the computational cost of simulations of laser-driven electron dynamics beyond the dipole approximation is the same as simulations adopting the dipole approximation. We illustrate the theory by time-dependent configuration interaction and coupled-cluster simulations of core-level absorption and circular dichroism spectra.

Real-time simulations of laser-driven electron dynamics contain information about molecular optical properties through all orders in response theory. These properties can be extracted by assuming convergence of the power series expansion of induced electric and magnetic multipole moments. However, the accuracy relative to analytical results from response theory quickly deteriorates for higher-order responses due to the presence of high-frequency oscillations in the induced multipole moment in the time domain. This problem has been ascribed to missing higher-order corrections. We here demonstrate that the deviations are caused by nonadiabatic effects arising from the finite-time ramping from zero to full strength of the external laser field. Three different approaches, two using a ramped wave and one using a pulsed wave, for extracting electrical properties from real-time time-dependent electronic-structure simulations are investigated. The standard linear ramp is compared to a quadratic ramp, which is found to yield highly accurate results for polarizabilities, and first and second hyperpolarizabilities, at roughly half the computational cost. Results for the third hyperpolarizability are presented along with a simple, computable measure of reliability.

In this work we present a derivation of the real-time time-dependent orbital-optimized M{\o}ller-Plesser TDOMP2 and its biorthogonal companion, time-dependent non-orthogonal OMP2 (TDNOMP2), theory starting from the time-dependent bivariational principle and a parametrization based on the exponential orbital-rotation operator formulation commonly used in time-independent molecular electronic-structure theory. We apply the TDOMP2 method to extract absorption spectra and frequency-dependent polarizabilities and first hyperpolarizabilities from real-time simulations, comparing the results with those obtained from conventional time-dependent coupled-cluster singles and doubles (TDCCSD) simulations and from its second-order approximation TDCC2. We also compare with results from CCSD and CC2 linear and quadratic response theory. Our results indicate that while TDOMP2 absorption spectra are of the same quality as TDCC2 spectra, frequency-dependent polarizabilities and hyperpolarizabilties from TDOMP2 simulations are significantly closer to TDCCSD results than those from TDCC2 simulations.