In classical mechanics the state of a dynamical system at any particular time can be described by the values of a set of coordinates and their conjugate momenta, thus, if the system has n degrees of freedom, by 2n numbers. In quantum mechanics, on the other hand, we have to describe a state of the system by a wave function involving a set of coordinates, thus by a function of n variables. The quantum description is, therefore, much more complicated than the classical one. Let us consider, however, an ensemble of systems in Gibbs' sense, i.e. not a large number of actual systems which could, perhaps, interact with one another, but a large number of hypothetical systems which are introduced to describe one actual system of which our knowledge is only of a statistical nature. The basis of the quantum treatment of such an ensemble has been given by Neumann. The description obtained by Neumann of an ensemble on the quantum theory is no more complicated than the corresponding classical description. Thus the quantum theory, which appears to such a disadvantage on the score of complication when applied to individual systems, recovers its own when applied to an ensemble. It is the object of the present note to examine this question more closely and to show how complete the analogy is between the quantum and classical treatments of an ensemble.
In classical mechanics the state of a dynamical system at any particular time can be described by the values of a set of coordinates and their conjugate momenta, thus, if the system has n degrees of freedom, by 2n numbers. In quantum mechanics, on the other hand, we have to describe a state of the system by a wave function involving a set of coordinates, thus by a function of n variables. The quantum description is, therefore, much more complicated than the classical one. Let us consider, however, an ensemble of systems in Gibbs' sense, i.e. not a large number of actual systems which could, perhaps, interact with one another, but a large number of hypothetical systems which are introduced to describe one actual system of which our knowledge is only of a statistical nature. The basis of the quantum treatment of such an ensemble has been given by Neumann. The description obtained by Neumann of an ensemble on the quantum theory is no more complicated than the corresponding classical description. Thus the quantum theory, which appears to such a disadvantage on the score of complication when applied to individual systems, recovers its own when applied to an ensemble. It is the object of the present note to examine this question more closely and to show how complete the analogy is between the quantum and classical treatments of an ensemble.
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)].
For systems which obey quantum mechanics, the formulation of the problem of treating large numbers of particles is, of course, somewhat different than it is for classical systems. The microscopic description of the system is provided by a wave function which (in the absence of spin) is a function of the classical coordinates { q i }. The mathematical model is provided by a Hamiltonian operator H which is often obtained from the corresponding classical Hamiltonian by the replacement p i → ( ħ /i)(∂/∂ q i ). In other cases the form of the Hamiltonian operator is simply postulated. The microscopic dynamics are provided by the Schrödinger equation i ħ (∂Ψ/∂ t ) = H Ψ which requires as initial condition the knowledge of the wave function Ψ({ q i }, t ) at some initial time t 0 . (Boundary conditions on Ψ({ q i }, t ) must be specified as part of the description of the model as well.) The results of experiments in quantum mechanics are characterized by operators, usually obtained, like the Hamiltonian, from their classical forms and termed observables. Operators associated with observables must be Hermitian. In general, the various operators corresponding to observables do not commute with one another. It is possible to find sets of commuting operators whose mutual eigenstates span the Hilbert space in which the wave function is confined by the Schrödinger equation and the boundary conditions. A set of such (time independent) eigenstates, termed ψ ν ( q ), is a basis for the Hilbert space.
In this chapter the basic concepts of quantum theory are formulated via a self-consistent set of six postulates. These six postulates deal with the following: operator images of physical observables; properties of state vectors and wave functions; the connection with experimental measurements; coordinate-space forms of those operators having classical analogs; the primary dynamical equation of quantum mechanics; and the completeness of the eigenvectors of observables. The postulates are the framework on which the edifice of quantum mechanics is constructed. A separate section is devoted to the statement of each postulate and to a few examples chosen to illustrate or apply it. A comprehensive set of applications is presented in Chapters 6 and 7, where various 1-D systems are used as the vehicles for illustrating the postulates. The basic theory of this chapter is extended and developed in Chapters 8 and 9. Since probability is a key interpretational concept, a brief introduction to this topic is given in an appendix. We emphasize that our formulation of quantum theory (in this chapter and throughout this book) is restricted to the case of pure states. This is not a major restriction. It means that, when we consider ensembles of individual quantum systems (e.g., single particles, atoms, nuclei, etc.), all of the systems in the ensemble are assumed to be in the same quantum state, not necessarily an eigenstate of a quantal operator. Ensembles for which this assumption does not hold are described by a density operator or density matrix rather than by a pure state; treatments of such cases can be found in references cited later.
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.
It is shown that application of the time-dependent variational principle to trial wave functions of a particular form results in a description of the quantum system as a classical system consisting of two mutually interacting parts: the ordinary classical analog of the quantum system and a (generally finite) number of additional degrees of freedom. This holds out the prospect of simulating quantum dynamics by coupling the corresponding classical system to additional degrees of freedom, and it is pointed out that such calculations have, in effect, already been performed. Since the trial wave functions that qualify for this treatment may be arbitrarily accurate, there is no limitation to the accuracy of the resulting classical description. In a particular case, the canonical coordinates and momenta for all degrees of freedom of the overall classical system are identified and investigated. In general, certain parameters appearing in the trial wave functions treated in this work may become redundant at certain times, leading to singularities in the Hamiltonian equations of motion. It is shown, however, that in certain such cases, a classical formalism for the quantum dynamics remains valid if the Hamiltonian function is suitably modified. Calculations applying the present Hamiltonian formalism are carried out that demonstrate that it has better convergence properties than an alternative, nonvariational, Hamiltonian technique suggested by previous work. These numerical studies also suggest that the observed spreading of wave packets and development of complicated structure in wave functions may be interpreted classically in terms of energy transfer among the various degrees of freedom in the overall classical system.
Introduction In our treatment of quantum mechanics we discussed the notion of the state of a given system. The state was specified by giving the eigenvalues of a maximal set of mutually commuting observables. This represented the most information we could have about a given system in principle. In practice, of course, we do not really have all of the necessary information about the eigenvalues available. No experimentalist sets out with the primary intention of obtaining results which correspond to a maximal set of mutually compatible observables for any macroscopic system, such as a gas. Even a very small volume of gas contains an enormous number of particles. The quantum mechanical state would require us to know a great deal of information about the internal states of the electrons and nuclei of the gas molecules. The experimentalist may only be concerned with gross properties such as the pressure, internal energy, or specific heat for example. Most of the detailed information about the detailed state of the system is neither observed, nor of direct interest. Instead we must content ourselves with a description of the average, or most likely, behaviour of the system in the probabilistic sense. The usual technique for doing this consists of the introduction of an ensemble, meaning a very large collection of systems whose measured macroscopic properties are the same as those for the system we are interested in, but whose unobserved properties can vary.