An implicit split-operator algorithm for the nonlinear time-dependent Schrödinger equation.

The explicit split-operator algorithm is often used for solving the linear and nonlinear time-dependent Schrödinger equations. However, when applied to certain nonlinear time-dependent Schrödinger equations, this algorithm loses time reversibility and second-order accuracy, which makes it very inefficient. Here, we propose to overcome the limitations of the explicit split-operator algorithm by abandoning its explicit nature. We describe a family of high-order implicit split-operator algorithms that are norm-conserving, time-reversible, and very efficient. The geometric properties of the integrators are proven analytically and demonstrated numerically on the local control of a two-dimensional model of retinal. Although they are only applicable to separable Hamiltonians, the implicit split-operator algorithms are, in this setting, more efficient than the recently proposed integrators based on the implicit midpoint method.

Time-reversible and norm-conserving high-order integrators for the nonlinear time-dependent Schrödinger equation: Application to local control theory.

The explicit split-operator algorithm has been extensively used for solving not only linear but also nonlinear time-dependent Schrödinger equations. When applied to the nonlinear Gross-Pitaevskii equation, the method remains time-reversible, norm-conserving, and retains its second-order accuracy in the time step. However, this algorithm is not suitable for all types of nonlinear Schrödinger equations. Indeed, we demonstrate that local control theory, a technique for the quantum control of a molecular state, translates into a nonlinear Schrödinger equation with a more general nonlinearity, for which the explicit split-operator algorithm loses time reversibility and efficiency (because it only has first-order accuracy). Similarly, the trapezoidal rule (the Crank-Nicolson method), while time-reversible, does not conserve the norm of the state propagated by a nonlinear Schrödinger equation. To overcome these issues, we present high-order geometric integrators suitable for general time-dependent nonlinear Schrödinger equations and also applicable to nonseparable Hamiltonians. These integrators, based on the symmetric compositions of the implicit midpoint method, are both norm-conserving and time-reversible. The geometric properties of the integrators are proven analytically and demonstrated numerically on the local control of a two-dimensional model of retinal. For highly accurate calculations, the higher-order integrators are more efficient. For example, for a wavefunction error of 10-9, using the eighth-order algorithm yields a 48-fold speedup over the second-order implicit midpoint method and trapezoidal rule, and a 400 000-fold speedup over the explicit split-operator algorithm.

Efficient geometric integrators for nonadiabatic quantum dynamics. II. The diabatic representation.

Exact nonadiabatic quantum evolution preserves many geometric properties of the molecular Hilbert space. In the first paper of this series ["Paper I," S. Choi and J. Vanícek, J. Chem. Phys. 150, 204112 (2019)], we presented numerical integrators of arbitrary-order of accuracy that preserve these geometric properties exactly even in the adiabatic representation, in which the molecular Hamiltonian is not separable into kinetic and potential terms. Here, we focus on the separable Hamiltonian in diabatic representation, where the split-operator algorithm provides a popular alternative because it is explicit and easy to implement, while preserving most geometric invariants. Whereas the standard version has only second-order accuracy, we implemented, in an automated fashion, its recursive symmetric compositions, using the same schemes as in Paper I, and obtained integrators of arbitrary even order that still preserve the geometric properties exactly. Because the automatically generated splitting coefficients are redundant, we reduce the computational cost by pruning these coefficients and lower memory requirements by identifying unique coefficients. The order of convergence and preservation of geometric properties are justified analytically and confirmed numerically on a one-dimensional two-surface model of NaI and a three-dimensional three-surface model of pyrazine. As for efficiency, we find that to reach a convergence error of 10-10, a 600-fold speedup in the case of NaI and a 900-fold speedup in the case of pyrazine are obtained with the higher-order compositions instead of the second-order split-operator algorithm. The pyrazine results suggest that the efficiency gain survives in higher dimensions.

On-the-fly ab initio semiclassical evaluation of time-resolved electronic spectra.

We present a methodology for computing vibrationally and time-resolved pump-probe spectra, which takes into account all vibrational degrees of freedom and is based on the combination of the thawed Gaussian approximation with on-the-fly ab initio evaluation of the electronic structure. The method is applied to the phenyl radical and compared with two more approximate approaches based on the global harmonic approximation-the global harmonic method expands both the ground- and excited-state potential energy surfaces to the second order about the corresponding minima, while the combined global harmonic/on-the-fly method retains the on-the-fly scheme for the excited-state wavepacket propagation. We also compare the spectra by considering their means and widths, and show analytically how these measures are related to the properties of the semiclassical wavepacket. We find that the combined approach is better than the global harmonic one in describing the vibrational structure, while the global harmonic approximation estimates better the overall means and widths of the spectra due to a partial cancellation of errors. Although the full-dimensional on-the-fly ab initio result seems to reflect the dynamics of only one mode, we show, by performing exact quantum calculations, that this simple structure cannot be recovered using a one-dimensional model. Yet, the agreement between the quantum and semiclassical spectra in this simple, but anharmonic model lends additional support for the full-dimensional ab initio thawed Gaussian calculation of the phenyl radical spectra. We conclude that the thawed Gaussian approximation provides a viable alternative to the expensive or unfeasible exact quantum calculations in cases, where low-dimensional models are not sufficiently accurate to represent the full system.

Semiclassical Approach to Photophysics Beyond Kasha's Rule and Vibronic Spectroscopy Beyond the Condon Approximation. The Case of Azulene.

Azulene is a prototypical molecule with an anomalous fluorescence from the second excited electronic state, thus violating Kasha's rule, and with an emission spectrum that cannot be understood within the Condon approximation. To better understand the photophysics and spectroscopy of azulene and other nonconventional molecules, we developed a systematic, general, and efficient computational approach combining the semiclassical dynamics of nuclei with ab initio electronic structure. First, to analyze the nonadiabatic effects, we complement the standard population dynamics by a rigorous measure of adiabaticity, estimated with the multiple-surface dephasing representation. Second, we propose a new semiclassical method for simulating non-Condon spectra, which combines the extended thawed Gaussian approximation with the efficient single-Hessian approach. S1 ← S0 and S2 ← S0 absorption and S2 → S0 emission spectra of azulene, recorded in a new set of experiments, agree very well with our calculations. We find that accuracy of the evaluated spectra requires the treatment of anharmonicity, Herzberg-Teller, and mode-mixing effects.