Finite-temperature vibronic spectra from the split-operator coherence thermofield dynamics.

We present a numerically exact approach for evaluating vibrationally resolved electronic spectra at finite temperatures using the coherence thermofield dynamics. In this method, which avoids implementing an algorithm for solving the von Neumann equation for coherence, the thermal vibrational ensemble is first mapped to a pure-state wavepacket in an augmented space, and this wavepacket is then propagated by solving the standard, zero-temperature Schrödinger equation with the split-operator Fourier method. We show that the finite-temperature spectra obtained with the coherence thermofield dynamics in a Morse potential agree exactly with those computed by Boltzmann-averaging the spectra of individual vibrational levels. Because the split-operator thermofield dynamics on a full tensor-product grid is restricted to low-dimensional systems, we briefly discuss how the accessible dimensionality can be increased by various techniques developed for the zero-temperature split-operator Fourier method.

High-order geometric integrators for the local cubic variational Gaussian wavepacket dynamics.

Gaussian wavepacket dynamics has proven to be a useful semiclassical approximation for quantum simulations of high-dimensional systems with low anharmonicity. Compared to Heller's original local harmonic method, the variational Gaussian wavepacket dynamics is more accurate, but much more difficult to apply in practice because it requires evaluating the expectation values of the potential energy, gradient, and Hessian. If the variational approach is applied to the local cubic approximation of the potential, these expectation values can be evaluated analytically, but they still require the costly third derivative of the potential. To reduce the cost of the resulting local cubic variational Gaussian wavepacket dynamics, we describe efficient high-order geometric integrators, which are symplectic, time-reversible, and norm-conserving. For small time steps, they also conserve the effective energy. We demonstrate the efficiency and geometric properties of these integrators numerically on a multidimensional, nonseparable coupled Morse potential.

Safety and Efficacy of Baseline Antiplatelet Treatment in Patients Undergoing Mechanical Thrombectomy for Ischemic Stroke: Antiplatelets Before Mechanical Thrombectomy.

PURPOSE: To investigate the safety and efficacy of baseline antiplatelet treatment in patients with acute ischemic stroke (AIS) undergoing mechanical thrombectomy (MT). MATERIALS AND METHODS: Baseline use of antiplatelet medication before MT for (AIS) may provide benefit on reperfusion and clinical outcome but could also carry an increased risk of intracranial hemorrhage (ICH). All consecutive patients with AIS and treated with MT with and without intravenous thrombolysis (IVT) between January 2012 and December 2019 in all centers performing MT nationwide were reviewed. Data were prospectively collected in national registries (eg, SITS-TBY and RES-Q). Primary outcome was functional independence (modified Rankin Scale 0-2) at 3 months; secondary outcome was ICH. RESULTS: Of the 4,351 patients who underwent MT, 1,750 (40%) and 666 (15%) were excluded owing to missing data from the functional independence and ICH outcome cohorts, respectively. In the functional independence cohort (n = 2,601), 771 (30%) patients received antiplatelets before MT. Favorable outcome did not differ in any antiplatelet, aspirin, and clopidogrel groups when compared with that in the no-antiplatelet group: odds ratio (OR), 1.00 (95% CI, 0.84-1.20); OR, 1.05 (95% CI, 0.86-1.27); and OR, 0.88 (95% CI, 0.55-1.41), respectively. In the ICH cohort (n = 3,685), 1095 (30%) patients received antiplatelets before MT. The rates of ICH did not increase in any treatment options (any antiplatelet, aspirin, clopidogrel, and dual antiplatelet groups) when compared with those in the no-antiplatelet group: OR, 1.03 (95% CI, 0.87-1.21); OR, 0.99 (95% CI, 0.83-1.18); OR, 1.10 (95% CI, 0.82-1.47); and OR, 1.43 (95% CI, 0.87-2.33), respectively. CONCLUSIONS: Antiplatelet monotherapy before MT did not improve functional independence or increase the risk of ICH.

Isotope Effects on the Electronic Spectra of Ammonia from Ab Initio Semiclassical Dynamics.

Despite its simplicity, the single-trajectory thawed Gaussian approximation has proven useful for calculating the vibrationally resolved electronic spectra of molecules with weakly anharmonic potential energy surfaces. Here, we show that the thawed Gaussian approximation can capture surprisingly well even more subtle observables, such as the isotope effects in the absorption spectra, and we demonstrate it on the four isotopologues of ammonia (NH3, NDH2, ND2H, and ND3). The differences in their computed spectra are due to the differences in the semiclassical trajectories followed by the four isotopologues, and the isotope effects─narrowing of the transition band and reduction of the peak spacing─are accurately described by this semiclassical method. In contrast, the adiabatic harmonic model shows a double progression instead of the single progression seen in the experimental spectra. The vertical harmonic model correctly shows only a single progression but fails to describe the anharmonic peak spacing. Analysis of the normal-mode activation upon excitation provides insight into the elusiveness of the symmetric stretching progression in the spectra.

Family of Gaussian wavepacket dynamics methods from the perspective of a nonlinear Schrödinger equation.

Many approximate solutions of the time-dependent Schrödinger equation can be formulated as exact solutions of a nonlinear Schrödinger equation with an effective Hamiltonian operator depending on the state of the system. We show that Heller's thawed Gaussian approximation, Coalson and Karplus's variational Gaussian approximation, and other Gaussian wavepacket dynamics methods fit into this framework if the effective potential is a quadratic polynomial with state-dependent coefficients. We study such a nonlinear Schrödinger equation in full generality: we derive general equations of motion for the Gaussian's parameters, demonstrate time reversibility and norm conservation, and analyze conservation of energy, effective energy, and symplectic structure. We also describe efficient, high-order geometric integrators for the numerical solution of this nonlinear Schrödinger equation. The general theory is illustrated by examples of this family of Gaussian wavepacket dynamics, including the variational and nonvariational thawed and frozen Gaussian approximations and their special limits based on the global harmonic, local harmonic, single-Hessian, local cubic, and local quartic approximations for the potential energy. We also propose a new method by augmenting the local cubic approximation with a single fourth derivative. Without substantially increasing the cost, the proposed "single-quartic" variational Gaussian approximation improves the accuracy over the local cubic approximation and, at the same time, conserves both the effective energy and symplectic structure, unlike the much more expensive local quartic approximation. Most results are presented in both Heller's and Hagedorn's parametrizations of the Gaussian wavepacket.

High-order geometric integrators for the variational Gaussian approximation.

Among the single-trajectory Gaussian-based methods for solving the time-dependent Schrödinger equation, the variational Gaussian approximation is the most accurate one. In contrast to Heller's original thawed Gaussian approximation, it is symplectic, conserves energy exactly, and may partially account for tunneling. However, the variational method is also much more expensive. To improve its efficiency, we symmetrically compose the second-order symplectic integrator of Faou and Lubich and obtain geometric integrators that can achieve an arbitrary even order of convergence in the time step. We demonstrate that the high-order integrators can speed up convergence drastically compared to the second-order algorithm and, in contrast to the popular fourth-order Runge-Kutta method, are time-reversible and conserve the norm and the symplectic structure exactly, regardless of the time step. To show that the method is not restricted to low-dimensional systems, we perform most of the analysis on a non-separable twenty-dimensional model of coupled Morse oscillators. We also show that the variational method may capture tunneling and, in general, improves accuracy over the non-variational thawed Gaussian approximation.

How to Find Molecules with Long-lasting Charge Migration?

Under certain conditions, the ionization of a molecule may create a superposition of electronic states, leading to ultrafast electron dynamics. If controlled, this motion could be used in attochemistry applications, but it has been shown that the decoherence induced by the nuclear motion typically happens in just a few femtoseconds. We recently developed an efficient algorithm for finding molecules exhibiting long-lasting electronic coherence and charge migration across the molecular structure after valence ionization. Here, we first explain why the but-3-ynal molecule is a promising candidate to study this type of ultrafast electron dynamics. Then, we use the 3-oxopropanenitrile molecule, which does not induce long-lasting charge migration in any of three different ionization scenarios, as an example demonstrating that several different properties must be fulfilled simultaneously to make the attochemistry applications possible.

Search for long-lasting electronic coherence using on-the-fly ab initio semiclassical dynamics.

Using a combination of high-level ab initio electronic structure methods with efficient on-the-fly semiclassical evaluation of nuclear dynamics, we performed a massive scan of small polyatomic molecules searching for a long-lasting oscillatory dynamics of the electron density triggered by the outer-valence ionization. We observed that in most of the studied molecules, either the sudden removal of an electron from the system does not lead to the appearance of the electronic coherence or the created coherences become damped by the nuclear rearrangement on a time scale of a few femtoseconds. However, we report several so far unexplored molecules with the electronic coherences lasting up to 10 fs, which can be good candidates for experimental studies. In addition, we present the full-dimensional simulations of the electronic coherences coupled to nuclear motion in several molecules which were studied previously only in the fixed nuclei approximation.

Efficient Semiclassical Evaluation of Electronic Coherences in Polyatomic Molecules.

Exposing a molecule to intense light pulses may bring this molecule to a nonstationary quantum state, thus launching correlated dynamics of electronic and nuclear subsystems. Although much had been achieved in the understanding of fundamental physics behind the electron-nuclear interactions and dynamics, accurate numerical simulations of light-induced processes taking place in polyatomic molecules remain a formidable challenge. Here, we review a recently developed theoretical approach for evaluating electronic coherences in molecules, in which the ultrafast electronic dynamics is coupled to nuclear motion. The presented technique, which combines accurate ab initio on-the-fly simulations of electronic structure with efficient semiclassical procedure to compute the dynamics of nuclear wave packets, is not only computationally efficient, but also can help shed light on the underlying physical mechanisms of decoherence and revival of the electronic coherences driven by nuclear rearrangement.

Core-Valence Attosecond Transient Absorption Spectroscopy of Polyatomic Molecules.

Tracing ultrafast processes induced by interaction of light with matter is often very challenging. In molecular systems, the initially created electronic coherence becomes damped by the slow nuclear rearrangement on a femtosecond timescale which makes real-time observations of electron dynamics in molecules particularly difficult. In this work, we report an extension of the theory underlying the attosecond transient absorption spectroscopy (ATAS) for the case of molecules, including a full account for the coupled electron-nuclear dynamics in the initially created wave packet, and apply it to probe the oscillations of the positive charge created after outer-valence ionization of the propiolic acid molecule. By taking advantage of element-specific core-to-valence transitions induced by x-ray radiation, we show that the resolution of ATAS makes it possible to trace the dynamics of electron density with atomic spatial resolution.

How important are the residual nonadiabatic couplings for an accurate simulation of nonadiabatic quantum dynamics in a quasidiabatic representation?

Diabatization of the molecular Hamiltonian is a standard approach to remove the singularities of nonadiabatic couplings at conical intersections of adiabatic potential energy surfaces. In general, it is impossible to eliminate the nonadiabatic couplings entirely-the resulting "quasidiabatic" states are still coupled by smaller but nonvanishing residual nonadiabatic couplings, which are typically neglected. Here, we propose a general method for assessing the validity of this potentially drastic approximation by comparing quantum dynamics simulated either with or without the residual couplings. To make the numerical errors negligible to the errors due to neglecting the residual couplings, we use the highly accurate and general eighth-order composition of the implicit midpoint method. The usefulness of the proposed method is demonstrated on nonadiabatic simulations in the cubic Jahn-Teller model of nitrogen trioxide and in the induced Renner-Teller model of hydrogen cyanide. We find that, depending on the system, initial state, and employed quasidiabatization scheme, neglecting the residual couplings can result in wrong dynamics. In contrast, simulations with the exact quasidiabatic Hamiltonian, which contains the residual couplings, always yield accurate results.

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.

High-order geometric integrators for representation-free Ehrenfest dynamics.

Ehrenfest dynamics is a useful approximation for ab initio mixed quantum-classical molecular dynamics that can treat electronically nonadiabatic effects. Although a severe approximation to the exact solution of the molecular time-dependent Schrödinger equation, Ehrenfest dynamics is symplectic, is time-reversible, and conserves exactly the total molecular energy as well as the norm of the electronic wavefunction. Here, we surpass apparent complications due to the coupling of classical nuclear and quantum electronic motions and present efficient geometric integrators for "representation-free" Ehrenfest dynamics, which do not rely on a diabatic or adiabatic representation of electronic states and are of arbitrary even orders of accuracy in the time step. These numerical integrators, obtained by symmetrically composing the second-order splitting method and exactly solving the kinetic and potential propagation steps, are norm-conserving, symplectic, and time-reversible regardless of the time step used. Using a nonadiabatic simulation in the region of a conical intersection as an example, we demonstrate that these integrators preserve the geometric properties exactly and, if highly accurate solutions are desired, can be even more efficient than the most popular non-geometric integrators.

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 Semiclassical Dynamics for Vibronic Spectroscopy beyond Harmonic, Condon, and Zero-Temperature Approximations.

Understanding light-induced processes in biological and human-made molecular systems is one of the main goals of physical chemistry. It has been known for years that the photoinduced dynamics of atomic nuclei can be studied by looking at the vibrational substructure of electronic absorption and emission spectra. However, theoretical simulation is needed to understand how dynamics translates into the spectral features. Here, we review several recent developments in the computation of vibrationally resolved electronic spectra (sometimes simply called 'vibronic' spectra). We present a theoretical approach for computing such spectra beyond the commonly used zero-temperature, Condon, and harmonic approximations. More specifically, we show how the on-the-fly ab initio thawed Gaussian approximation, which partially includes anharmonicity effects, can be combined with the thermo-field dynamics to treat non-zero temperature and with the Herzberg-Teller correction to include non-Condon effects. The combined method, which can treat all three effects, is applied to compute the S1 ← S0 and S2 ← S0 absorption spectra of azulene.

On-the-Fly Ab Initio Semiclassical Evaluation of Electronic Coherences in Polyatomic Molecules Reveals a Simple Mechanism of Decoherence.

Irradiation of a molecular system by an intense laser field can trigger dynamics of both electronic and nuclear subsystems. The lighter electrons usually move on much faster, attosecond timescale but the slow nuclear rearrangement damps ultrafast electronic oscillations, leading to the decoherence of the electronic dynamics within a few femtoseconds. We show that a simple, single-trajectory semiclassical scheme can evaluate the electronic coherence time in polyatomic molecules accurately by demonstrating an excellent agreement with full-dimensional quantum calculations. In contrast to numerical quantum methods, the semiclassical one reveals the physical mechanism of decoherence beyond the general blame on nuclear motion. In the propiolic acid, the rate of decoherence and the large deviation from the static frequency of electronic oscillations are quantitatively described with just two semiclassical parameters-the phase space distance and signed area between the trajectories moving on two electronic surfaces. Because it evaluates the electronic structure on the fly, the semiclassical technique avoids the "curse of dimensionality" and should be useful for preselecting molecules for experimental studies.

On-the-fly ab initio semiclassical evaluation of vibronic spectra at finite temperature.

To compute and analyze vibrationally resolved electronic spectra at zero temperature, we have recently implemented the on-the-fly ab initio extended thawed Gaussian approximation [A. Patoz et al., J. Phys. Chem. Lett. 9, 2367 (2018)], which accounts for anharmonicity, mode-mode coupling, and Herzberg-Teller effects. Here, we generalize this method in order to evaluate spectra at non-zero temperature. In line with thermo-field dynamics, we transform the von Neumann evolution of the coherence component of the density matrix to the Schrödinger evolution of a wavefunction in an augmented space with twice as many degrees of freedom. Due to the efficiency of the extended thawed Gaussian approximation, this increase in the number of coordinates results in nearly no additional computational cost. More specifically, compared to the original, zero-temperature approach, the finite-temperature method requires no additional ab initio electronic structure calculations. At the same time, the new approach allows for a clear distinction among finite-temperature, anharmonicity, and Herzberg-Teller effects on spectra. We show, on a model Morse system, the advantages of the finite-temperature thawed Gaussian approximation over the commonly used global harmonic methods and apply it to evaluate the symmetry-forbidden absorption spectrum of benzene, where all of the aforementioned effects contribute.

Which form of the molecular Hamiltonian is the most suitable for simulating the nonadiabatic quantum dynamics at a conical intersection?

Choosing an appropriate representation of the molecular Hamiltonian is one of the challenges faced by simulations of the nonadiabatic quantum dynamics around a conical intersection. The adiabatic, exact quasidiabatic, and strictly diabatic representations are exact and unitary transforms of each other, whereas the approximate quasidiabatic Hamiltonian ignores the residual nonadiabatic couplings in the exact quasidiabatic Hamiltonian. A rigorous numerical comparison of the four different representations is difficult because of the exceptional nature of systems where the four representations can be defined exactly and the necessity of an exceedingly accurate numerical algorithm that avoids mixing numerical errors with errors due to the different forms of the Hamiltonian. Using the quadratic Jahn-Teller model and high-order geometric integrators, we are able to perform this comparison and find that only the rarely employed exact quasidiabatic Hamiltonian yields nearly identical results to the benchmark results of the strictly diabatic Hamiltonian, which is not available in general. In this Jahn-Teller model and with the same Fourier grid, the commonly employed approximate quasidiabatic Hamiltonian led to inaccurate wavepacket dynamics, while the Hamiltonian in the adiabatic basis was the least accurate, due to the singular nonadiabatic couplings at the conical intersection.

On-the-fly ab initio semiclassical evaluation of third-order response functions for two-dimensional electronic spectroscopy.

Ab initio computation of two-dimensional electronic spectra is an expanding field, whose goal is improving upon simple, few-dimensional models often employed to explain experiments. Here, we propose an accurate and computationally affordable approach, based on the single-trajectory semiclassical thawed Gaussian approximation, to evaluate two-dimensional electronic spectra. Importantly, the method is exact for arbitrary harmonic potentials with mode displacement, changes in the mode frequencies, and inter-mode coupling (Duschinsky effect), but can also account partially for the anharmonicity of the involved potential energy surfaces. We test its accuracy on a set of model Morse potentials and use it to study anharmonicity and Duschinsky effects on the linear and two-dimensional electronic spectra of phenol. We find that in this molecule, the anharmonicity effects are weak, whereas the Duschinsky rotation and the changes in the mode frequencies must be included in accurate simulations. In contrast, the widely used displaced harmonic oscillator model captures only the basic physics of the problem but fails to reproduce the correct vibronic lineshape.

A time-reversible integrator for the time-dependent Schrödinger equation on an adaptive grid.

One of the most accurate methods for solving the time-dependent Schrödinger equation uses a combination of the dynamic Fourier method with the split-operator algorithm on a tensor-product grid. To reduce the number of required grid points, we let the grid move together with the wavepacket but find that the naïve algorithm based on an alternate evolution of the wavefunction and grid destroys the time reversibility of the exact evolution. Yet, we show that the time reversibility is recovered if the wavefunction and grid are evolved simultaneously during each kinetic or potential step; this is achieved by using the Ehrenfest theorem together with the splitting method. The proposed algorithm is conditionally stable, symmetric, and time-reversible and conserves the norm of the wavefunction. The preservation of these geometric properties is shown analytically and demonstrated numerically on a three-dimensional harmonic model and collinear model of He-H2 scattering. We also show that the proposed algorithm can be symmetrically composed to obtain time-reversible integrators of an arbitrary even order. We observed 10 000-fold speedup by using the tenth-order instead of the second-order method to obtain a solution with a time discretization error below 10-9. Moreover, using the adaptive grid instead of the fixed grid resulted in a 64-fold reduction in the required number of grid points in the harmonic system and made it possible to simulate the He-H2 scattering for six times longer while maintaining reasonable accuracy. Applicability of the algorithm to high-dimensional quantum dynamics is demonstrated using the strongly anharmonic eight-dimensional Hénon-Heiles model.