ABSTRACT
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.
ABSTRACT
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.
ABSTRACT
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.
Subject(s)
Brain Ischemia , Ischemic Stroke , Mechanical Thrombolysis , Stroke , Humans , Stroke/diagnostic imaging , Stroke/therapy , Ischemic Stroke/diagnostic imaging , Ischemic Stroke/therapy , Thrombolytic Therapy/adverse effects , Brain Ischemia/diagnostic imaging , Brain Ischemia/therapy , Thrombectomy/adverse effects , Clopidogrel/adverse effects , Treatment Outcome , Intracranial Hemorrhages/chemically induced , Aspirin/adverse effects , Mechanical Thrombolysis/adverse effectsABSTRACT
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.
ABSTRACT
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.
Subject(s)
Algorithms , Quantum Theory , Physical Phenomena , Motion , Normal DistributionABSTRACT
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.
ABSTRACT
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.
ABSTRACT
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.
ABSTRACT
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.
ABSTRACT
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.
ABSTRACT
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.
ABSTRACT
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.
ABSTRACT
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.
ABSTRACT
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.
ABSTRACT
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.
ABSTRACT
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.
ABSTRACT
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.
ABSTRACT
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.
ABSTRACT
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.
ABSTRACT
Geometric integrators of the Schrödinger equation conserve exactly many invariants of the exact solution. Among these integrators, the split-operator algorithm is explicit and easy to implement but, unfortunately, is restricted to systems whose Hamiltonian is separable into kinetic and potential terms. Here, we describe several implicit geometric integrators applicable to both separable and nonseparable Hamiltonians and, in particular, to the nonadiabatic molecular Hamiltonian in the adiabatic representation. These integrators combine the dynamic Fourier method with the recursive symmetric composition of the trapezoidal rule or implicit midpoint method, which results in an arbitrary order of accuracy in the time step. Moreover, these integrators are exactly unitary, symplectic, symmetric, time-reversible, and stable and, in contrast to the split-operator algorithm, conserve energy exactly, regardless of the accuracy of the solution. The order of convergence and conservation of geometric properties are proven analytically and demonstrated numerically on a two-surface NaI model in the adiabatic representation. Although each step of the higher order integrators is more costly, these algorithms become the most efficient ones if higher accuracy is desired; a thousand-fold speedup compared to the second-order trapezoidal rule (the Crank-Nicolson method) was observed for a wavefunction convergence error of 10-10. In a companion paper [J. Roulet, S. Choi, and J. Vanícek, J. Chem. Phys. 150, 204113 (2019)], we discuss analogous, arbitrary-order compositions of the split-operator algorithm and apply both types of geometric integrators to a higher-dimensional system in the diabatic representation.