Modus Operandi of a Pedalo-Type Molecular Switch: Insight from Dynamics and Theoretical Spectroscopy.

Molecular switches which can be triggered by light to interconvert between two or more well-defined conformation differing in their chemical or physical properties are fundamental for the development of materials with on-demand functionalities. Recently, a novel molecular switch based on a the azodicarboxamide core has been reported. It exhibits a volume-conserving conformational change upon excitation, making it a promising candidate for embedding in confined environments. In order to rationally implement and efficiently utilize the azodicarboxamide molecular switch, detailed insight into the coordinates governing the excited-state dynamics is needed. Here, we report a detailed comparative picture of the molecular motion at the atomic level in the presence and absence of explicit solvent. Our hybrid quantum mechanics/molecular mechanics (QM/MM) excited state simulations reveal that, although the energy landscape is slightly modulated by the solvation, the light-induced motion is dominated by a bending-assisted pedalo-type motion independent of the solvation. To support the predicted mechanism, we simulate time-resolved IR spectroscopy from first principles, thereby resolving fingerprints of the light-induced switching process. Our calculated time-resolved data are in good agreement with previously reported measured spectra.

Quantum-Classical Hybrid Systems and Ehrenfest's Theorem.

The conceptual analysis of quantum mechanics brings to light that a theory inherently consistent with observations should be able to describe both quantum and classical systems, i.e., quantum-classical hybrids. For example, the orthodox interpretation of measurements requires the transient creation of quantum-classical hybrids. Despite its limitations in defining the classical limit, Ehrenfest's theorem makes the simplest contact between quantum and classical mechanics. Here, we generalized the Ehrenfest theorem to bipartite quantum systems. To study quantum-classical hybrids, we employed a formalism based on operator-valued Wigner functions and quantum-classical brackets. We used this approach to derive the form of the Ehrenfest theorem for quantum-classical hybrids. We found that the time variation of the average energy of each component of the bipartite system is equal to the average of the symmetrized quantum dissipated power in both the quantum and the quantum-classical case. We expect that these theoretical results will be useful both to analyze quantum-classical hybrids and to develop self-consistent numerical algorithms for Ehrenfest-type simulations.

A Quantum-Classical Model of Brain Dynamics.

The study of the human psyche has elucidated a bipartite structure of logic reflecting the quantum-classical nature of the world. Accordingly, we posited an approach toward studying the brain by means of the quantum-classical dynamics of a mixed Weyl symbol. The mixed Weyl symbol can be used to describe brain processes at the microscopic level and, when averaged over an appropriate ensemble, can provide a link to the results of measurements made at the meso and macro scale. Within this approach, quantum variables (such as, for example, nuclear and electron spins, dipole momenta of particles or molecules, tunneling degrees of freedom, and so on) can be represented by spinors, whereas the electromagnetic fields and phonon modes can be treated either classically or semi-classically in phase space by also considering quantum zero-point fluctuations. Quantum zero-point effects can be incorporated into numerical simulations by controlling the temperature of each field mode via coupling to a dedicated Nosé-Hoover chain thermostat. The temperature of each thermostat was chosen in order to reproduce quantum statistics in the canonical ensemble. In this first paper, we introduce a general quantum-classical Hamiltonian model that can be tailored to study physical processes at the interface between the quantum and the classical world in the brain. While the approach is discussed in detail, numerical calculations are not reported in the present paper, but they are planned for future work. Our theory of brain dynamics subsumes some compatible aspects of three well-known quantum approaches to brain dynamics, namely the electromagnetic field theory approach, the orchestrated objective reduction theory, and the dissipative quantum model of the brain. All three models are reviewed.

PyUNIxMD: A Python-based excited state molecular dynamics package.

Theoretical/computational description of excited state molecular dynamics is nowadays a crucial tool for understanding light-matter interactions in many materials. Here we present an open-source Python-based nonadiabatic molecular dynamics program package, namely PyUNIxMD, to deal with mixed quantum-classical dynamics for correlated electron-nuclear propagation. The PyUNIxMD provides many interfaces for quantum chemical calculation methods with commercial and noncommercial ab initio and semiempirical quantum chemistry programs. In addition, the PyUNIxMD offers many nonadiabatic molecular dynamics algorithms such as fewest-switch surface hopping and its derivatives as well as decoherence-induced surface hopping based on the exact factorization (DISH-XF) and coupled-trajectory mixed quantum-classical dynamics (CTMQC) for general purposes. Detailed structures and flows of PyUNIxMD are explained for the further implementations by developers. We perform a nonadiabatic molecular dynamics simulation for a molecular motor system as a simple demonstration.

WavePacket: A Matlab package for numerical quantum dynamics. III. Quantum-classical simulations and surface hopping trajectories.

WavePacket is an open-source program package for numerical simulations in quantum dynamics. Building on the previous Part I (Schmidt and Lorenz, Comput. Phys. Commun. 2017, 213, 223] and Part II (Schmidt and Hartmann, Comput. Phys. Commun. 2018, 228, 229] which dealt with quantum dynamics of closed and open systems, respectively, the present Part III adds fully classical and mixed quantum-classical propagation techniques to WavePacket. There classical phase-space densities are sampled by trajectories which follow (diabatic or adiabatic) potential energy surfaces. In the vicinity of (genuine or avoided) intersections of those surfaces, trajectories may switch between them. To model these transitions, two classes of stochastic algorithms have been implemented: (1) Tully's fewest switches surface hopping and (2) Landau-Zener-based single switch surface hopping. The latter one offers the advantage of being based on adiabatic energy gaps only, thus not requiring nonadiabatic coupling information any more. The present work describes the MATLAB version of WavePacket 6.1.0, which is essentially an object-oriented rewrite of previous versions, allowing to perform fully classical, quantum-classical and quantum-mechanical simulations on an equal footing, that is, for the same physical system described by the same WavePacket input. The software package is hosted and further developed at the Sourceforge platform, where also extensive Wiki-documentation as well as numerous worked-out demonstration examples with animated graphics are available. © 2019 Wiley Periodicals, Inc.

Complex Fluid Models of Mixed Quantum-Classical Dynamics.

Several methods in nonadiabatic molecular dynamics are based on Madelung's hydrodynamic description of nuclear motion, while the electronic component is treated as a finite-dimensional quantum system. In this context, the quantum potential leads to severe computational challenges and one often seeks to neglect its contribution, thereby approximating nuclear motion as classical. The resulting model couples classical hydrodynamics for the nuclei to the quantum motion of the electronic component, leading to the structure of a complex fluid system. This type of mixed quantum-classical fluid models has also appeared in solvation dynamics to describe the coupling between liquid solvents and the quantum solute molecule. While these approaches represent a promising direction, their mathematical structure requires a certain care. In some cases, challenging higher-order gradients make these equations hardly tractable. In other cases, these models are based on phase-space formulations that suffer from well-known consistency issues. Here, we present a new complex fluid system that resolves these difficulties. Unlike common approaches, the current system is obtained by applying the fluid closure at the level of the action principle of the original phase-space model. As a result, the system inherits a Hamiltonian structure and retains energy/momentum balance. After discussing some of its structural properties and dynamical invariants, we illustrate the model in the case of pure-dephasing dynamics. We conclude by presenting some invariant planar models.

A Simple Solution to the Trivial Crossing Problem in Surface Hopping.

Surface hopping studies on supramolecular and nanoscale systems suffer severely from the trivial crossing problem, arising due to high density of adiabatic potential energy surfaces. We present a straightforward solution to the problem by introducing a self-consistency test to the well-known fewest switches surface hopping (FSSH) procedure. If the test is failed, the hopping probabilities are corrected with a simple procedure. The novel self-consistent fewest switches surface hopping (SC-FSSH) approach is applied to the Holstein Hamiltonian to study the time-dependence of the electron population. Already in the five-state system, SC-FSSH allows us to reduce the simulation time 10(4)-fold to achieve the FSSH accuracy. The reliable performance and simple formulation of SC-FSSH greatly expands the applicability range of the surface hopping method.

Accurate Calculations of Rotationally Inelastic Scattering Cross Sections Using Mixed Quantum/Classical Theory.

For computational treatment of rotationally inelastic scattering of molecules, we propose to use the mixed quantum/classical theory, MQCT. The old idea of treating translational motion classically, while quantum mechanics is used for rotational degrees of freedom, is developed to the new level and is applied to Na + N2 collisions in a broad range of energies. Comparison with full-quantum calculations shows that MQCT accurately reproduces all, even minor, features of energy dependence of cross sections, except scattering resonances at very low energies. The remarkable success of MQCT opens up wide opportunities for computational predictions of inelastic scattering cross sections at higher temperatures and/or for polyatomic molecules and heavier quenchers, which is computationally close to impossible within the full-quantum framework.