The emergence of experiments capable of probing quantum dynamics at the single-molecule level requires the development of new theoretical tools capable of simulating and analyzing these dynamics beyond an ensemble-averaged description. In this article, we present an efficient method for sampling and simulating the dynamics of the individual quantum systems that make up an ensemble and apply it to study the nonequilibrium dynamics of the ubiquitous spin-boson model. We generate an ensemble of single-system trajectories, and we analyze this trajectory ensemble using tools from classical statistical mechanics. Our results demonstrate that the dynamics of quantum coherence is highly heterogeneous at the single-system level due to variations in the initial bath configuration, which significantly affects the transient exchange of coherence between the system and its bath. We observe that single systems tend to retain coherence over time scales longer than that of the ensemble. We also compute a novel thermodynamic entanglement entropy that quantifies a thermodynamic driving force favoring system-bath entanglement.
The emergence of experiments capable of probing quantum dynamics at the single-molecule level requires the development of new theoretical tools capable of simulating and analyzing these dynamics beyond an ensemble-averaged description. In this article, we present an efficient method for sampling and simulating the dynamics of the individual quantum systems that make up an ensemble and apply it to study the nonequilibrium dynamics of the ubiquitous spin-boson model. We generate an ensemble of single-system trajectories, and we analyze this trajectory ensemble using tools from classical statistical mechanics. Our results demonstrate that the dynamics of quantum coherence is highly heterogeneous at the single-system level due to variations in the initial bath configuration, which significantly affects the transient exchange of coherence between the system and its bath. We observe that single systems tend to retain coherence over time scales longer than that of the ensemble. We also compute a novel thermodynamic entanglement entropy that quantifies a thermodynamic driving force favoring system-bath entanglement.
Recently, the possibility of investigating single molecule, or single spin observables, as well as the necessity of a better understanding of the mechanisms underlying quantum dynamics in order to obtain nanoscale devices and nanostructered materials suitable for quantum computing tasks, have revived the interest in foundational aspects of quantum statistical mechanics. This thesis aims to give a contribution to this field by re-considering the statistical characterization of a quantum system at the light of some paradigmatic changes in our understanding of quantum theory which have taken place in the last two decades. In particular the impressive development of quantum information theory has changed the perceptions of quantum entanglement: for a long time it has been considered a somewhat paradoxical property of the matter at the atomic scale, but now it is regarded as an essential and ubiquitous phenomenon whose consequences are affecting the very macroscopic world that we experience. Still the decoherence program has brought out the importance of considering a quantum system together with its environment in order to clarify some key aspects of quantum dynamics. Thus, we start from the idea that quantum correlations are ubiquitous and somewhat uncontrollable in systems with many degrees of freedom which are typically considered in statistical mechanics. As a consequence we assume the standpoint that quantum statistical mechanics has not to be based on the underlying idea of a collection of many, independent quantum systems but rather it has to emerge at the level of a global wavefunction (pure state) which describes the system as well as its environment as a whole. In order to investigate the consequences of these assumptions we study the equilibrium distribution of an isolated quantum system. This is defined, in analogy with the ergodic foundations of classical statistical mechanics, on the basis of the time evolution of the quantum state. Then, we study the emergence of thermodynamic properties in a quantum system by studying the probability distribution of some function of interests, as the entropy and the equilibrium state of a subsystem, on Ensembles of Pure States. Such a probability distribution is derived from the geometry of the Hilbert space, and the theoretical tools suitable for its characterization are developed. On the one hand we perform a numerical sampling of the ensemble distributions by employing Monte Carlo techniques, on the other hand simpler analytical approximation of the geometrical distributions are derived by means of a maximum entropy principle. Model systems composed of an ensemble of spins are chosen to illustrate the salient features which emerge from the developed theoretical framework: the main point is that the Ensemble Distributions of “thermodynamic observables” (entropy or equilibrium state of a subsystem) are sharply peaked around a typical value. From the analysis it emerges that each of the overwhelming majority of the wavefunctions which has appreciable weight in the considered ensemble, is characterized by the same value of the “macroscopic” functions. This is a striking evidence of the “typicality” of these properties. In the essence, our impossibility to know the state of the system in detail does not matter, just for the remarkable fact that almost all quantum states behave essentially in the same way. By virtue of this typicality the study of the behaviour of the typical values of the thermodynamic function become meaningful. Notably, under certain conditions, one recovers the results of standard statistical mechanics, that is, the equilibrium average of the state of a subsystem can be cast in the Boltzmann canonical form at the temperature given by the usual thermodynamical relation . In the second part of the thesis we consider the dynamical aspects of the equilibrium state of a subsystem interacting with its environment. The fluctuations around the equilibrium average critically depends on the entanglement between the system and the environment and on the form of the interaction Hamiltonian. The connection between the dynamics of the fluctuations of an observable at the equilibrium and the relaxation toward the equilibrium from a “non typical” initial value is also investigated with the aid of simple model systems. The study presented in this thesis was partly motivated by a critical analysis of the statistical methods available for the theoretical modelling of magnetic resonance experiments. One of these, the Stochastic Liouville Equation, has been employed in a work completed during the first year of my Ph.D. program in order to interpret some feature of a two dimensional electron spin resonance experiment, [Fresch B., Frezzato D., Moro G. J., Kothe G., Freed J. H.; J. Phys. Chem. B., 110, 24238, (2006)].
In this Letter, we present a new method for simulating quantum processes in the context of classical molecular dynamics simulations. The approach is based on solving numerically the quantum Liouville equation in the Wigner representation using ensembles of classical trajectories. Quantum effects arise in this formulation as a breakdown of the statistical independence of the ensemble. New interaction forces between ensemble members are derived, which require the trajectory ensemble representing the state to evolve as an entangled, unified whole. The method is applied to the simulation of quantum tunneling in a one-dimensional model system, yielding excellent agreement with exact quantum calculations.
A hybrid method is proposed to propagate system-bath quantum dynamics that use both basis functions and coupled quantum trajectories. In it, the bath is represented with an ensemble of Bohmian trajectories while the system degrees of freedom are accounted through reduced density matrices. By retaining the Hilbert space structure for the system, the method is able to capture interference processes that are challenging to describe in Bohmian dynamics due to singularities that these processes introduce in the quantum potential. By adopting quantum trajectories to represent the bath, the method beats the exponential scaling of the computational cost with the bath size. This combination makes the method suitable for large-scale ground and excited state fully quantum molecular dynamics simulations. Equations of motion for the quantum trajectories and reduced density matrices are derived from the Schrödinger equation and a computational algorithm to solve these equations is proposed. Through computations in two-dimensional model systems, the method is shown to offer an accurate description of subsystem observables and of quantum decoherence, which is difficult to obtain when the quantum nature of the bath is ignored. The scaling of the method is demonstrated using a model with 21 degrees of freedom. The limit of independent trajectories is recovered when the mass of bath degrees of freedom is much larger than the one of the system, in agreement with mixed quantum-classical descriptions.
In this paper, we describe a new method for simulating quantum processes using classical-like molecular dynamics. The approach is based on solving the quantum Liouville equation in the Wigner representation using ensembles of classical trajectories in phase space. The nonlocality of quantum mechanics is incorporated in the trajectory representation as nonclassical interactions between the members of the ensemble, leading to an entanglement of their evolution. The statistical independence of the individual trajectories making up an ensemble in the classical limit is lost when quantum effects are included, and the entire state of the system must be propagated as a unified whole. We develop the formalism and its numerical implementation, and illustrate its application on two model problems of quantum mechanical tunneling: escape from a metastable well and wave packet penetration of the Eckhart barrier.
A hybrid method is proposed to propagate system-bath quantum dynamics that use both basis functions and coupled quantum trajectories. In it, the bath is represented with an ensemble of Bohmian trajectories while the system degrees of freedom are accounted through reduced density matrices. By retaining the Hilbert space structure for the system, the method is able to capture interference processes that are challenging to describe in Bohmian dynamics due to singularities that these processes introduce in the quantum potential. By adopting quantum trajectories to represent the bath, the method beats the exponential scaling of the computational cost with the bath size. This combination makes the method suitable for large-scale ground and excited state fully quantum molecular dynamics simulations. Equations of motion for the quantum trajectories and reduced density matrices are derived from the Schrodinger equation and a computational algorithm to solve these equations is proposed. Through computations in two-dimensional model systems, the method is shown to offer an accurate description of subsystem observables and of quantum decoherence, which is difficult to obtain when the quantum nature of the bath is ignored. The scaling of the method is demonstrated using a model with 21 degrees of freedom. The limit of independent trajectories is recovered when the mass of bath degrees of freedom is much larger than the one of the system, in agreement with mixed quantum-classical descriptions.
Despite early developments on the foundations of quantum mechanics concern the wave function, quantum statistics has been developed with the density matrix formalism, leading to very important results in explaining molecular observations. Only recently, several authors argued a new interpretation by focusing on the wave function representing the quantum state of an isolated system, showing how a single wave function can exhibit statistical properties and generate the same results expected in the standard quantum statistical framework. Starting from these results, investigation on the foundations of quantum statistical mechanics has gained recently a renewed interest. As a matter of fact, the possibility of studying single molecule properties as well as the need of a better understanding of the effect of quantum dynamics, in order to develop new nanoscaled materials suitable to quantum computing tasks, have opened new intriguing questions leading to quantum statistical approaches far from being well understood and accepted. In this framework the behaviour of a single realization of quantum systems has gained a central role in the description of molecular systems. Furthermore, in recent years, an increasing number of studies has been presented on quantum dynamics through the numerical solution of the Schrödinger equation for systems of interacting components. These studies demonstrates that Quantum Dynamics Simulations could be a practicable route. In order to study phenomena such as dissipation, relaxation and thermalization, the focus has to be moved from isolated molecules to modular systems made of mutual interacting components, with model Hamiltonians possessing a sufficiently low dimensional representation. An important issue concerns the rules to be employed for the choice of the initial quantum state for the simulation of isolated systems. As long as one considers molecular degrees of freedom interacting with a (model) environment, there are no reasons to select a particular quantum state for the overall system and, therefore, a random choice has to be performed amongst a well defined statistical ensemble of pure states. Furthermore, one would like to operate a choice assuring that the simulation of the system is in a well defined thermal state with given temperature. This necessarily calls for a statistical description like for classical systems. A useful parametrization of the wave function will be presented in order to highlight the most appropriate variable for the statistical analysis. In particular some of them, called phases, retain all the dynamical information whereas the others, called populations, are the constants of motion. The latter, in particular are very important because they describe the equilibrium properties strictly related to the thermodynamic description. Since the dynamics of wave function does not supply any information about populations, the definition of a probability distribution on these variables is required. Different probability distributions on populations have been proposed, only on the bases of reasonable assumptions and they validation has been performed only with a posteriori considerations. In particular Fresch and Moro have demonstrated how the agreement with thermodynamics can be employed to discriminate different probability distributions on pure states. This had led the uniform ensembles to be, up to now, the most self-consistent models for quantum pure states. However, in this thesis I will highlight a drawback of the uniform distribution ensemble that can be described as follow: if we bring into contact two systems, even through a perturbative interaction, we are not able to describe the equilibrium properties after the interaction within the uniform distribution statistics, since the uniform character is lost. It represents a severe shortcoming of the statistical ensemble from a methodological point of view, since closed systems can be always considered as the result of interaction among previously isolated systems. On the other hand this drawback introduces a further requirement of a different nature that can be used for the definition of a new statistical ensemble. In this work I intend to find and characterize a statistical ensemble for populations that overcomes the drawbacks of the uniform distribution of pure states. The invariance of the thermal state in the coupling of identical systems will be used as a guideline in the definition af a new probability distribution on populations. Such an ensemble for pure states, called Thermalization Resilient Ensemble, provides a convenient framework for treating the interactions between quantum systems, as long as the structure of the statistical distribution is preserved and the identification of thermodynamic properties is assured. In perspective it should be the privileged statistical ensemble to implement Quantum Dynamics Simulations. Once the average properties of the Thermalization Resilient Ensemble have been introduced, I will obtain a probability distribution on pure states with the use of a geometrical analysis on the Hilbert space. The surface elements of an ellipsoidal manifold will be related to the probability density on populations. As a matter of fact the explicit form of the probability distribution is a prerequisite in order to perform Quantum Dynamical Simulations. However the results obtained through the geometrical analysis cannot be easily extended to systems with unbounded energy spectrum and an alternative strategy has been developed. A scaling algorithm on the basis of the uniform statistical ensemble will be described and this allows a well defined sampling of a probability distribution with desired averages. In this framework I demonstrate the emergence of thermodynamic behavior in the limit of macroscopic systems. In the last part of the thesis I consider the dynamical features of the thermalization experiment. Two identical systems, initially at different temperature, will be brought in interaction and the analysis of the final equilibrium state will be performed for two different generic forms of the interaction Hamiltonian, highlighting how the statistical approach can be very useful in the definition of the equilibrium in complex quantum systems.