A theoretical-computational procedure, recently proposed for modelling Vibrational Energy Relaxation (VER) processes of a molecule (Quantum Center, QC) embedded in a complex atomic-molecular system, is extended and applied for analyzing in detail the features of the QC density matrix (DM) temporal evolution. The results, obtained using aqueous azide ion as a case study, show the total lack of coherence in the DM, when the system is prepared to be initially in a pure vibrational eigenstate. This finding is fully in line with the statistical interpretation of the process typically adopted also in the experimental studies where the relaxation processes are all described within the typical schemes of chemical kinetics. Consistently, when the initial vibrational state corresponds to an eigenstate mixture, although initially coherent, the DM relaxes to a fully incoherent condition with a mean lifetime related to the one of the diagonal elements relaxation. These specific DM features turn out to be essentially governed by the thermal equilibrium condition of the atomic-molecular classical coordinates which drive the ensemble of the quantum-trajectories toward the observed statistical regime. Finally, from the analysis of a single long timescale quantum vibrational trajectory it also clearly emerges its ergodic behaviour.
A theoretical-computational procedure, recently proposed for modelling Vibrational Energy Relaxation (VER) processes of a molecule (Quantum Center, QC) embedded in a complex atomic-molecular system, is extended and applied for analyzing in detail the features of the QC density matrix (DM) temporal evolution. The results, obtained using aqueous azide ion as a case study, show the total lack of coherence in the DM, when the system is prepared to be initially in a pure vibrational eigenstate. This finding is fully in line with the statistical interpretation of the process typically adopted also in the experimental studies where the relaxation processes are all described within the typical schemes of chemical kinetics. Consistently, when the initial vibrational state corresponds to an eigenstate mixture, although initially coherent, the DM relaxes to a fully incoherent condition with a mean lifetime related to the one of the diagonal elements relaxation. These specific DM features turn out to be essentially governed by the thermal equilibrium condition of the atomic-molecular classical coordinates which drive the ensemble of the quantum-trajectories toward the observed statistical regime. Finally, from the analysis of a single long timescale quantum vibrational trajectory it also clearly emerges its ergodic behaviour.
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.
The coherent dynamics of many body quantum system is nowadays an experimental reality: by means of the cold atoms in optical lattices, many Hamiltonians and time-dependent perturbations can be engineered. In this Thesis we discuss what happens in these systems when a periodic perturbation is applied. Thanks to Floquet theory, we can see that -- if the Floquet spectrum obeys certain continuity conditions possible in the thermodynamic limit-- dephasing among Floquet quasi-energies makes local observables relax to a periodic steady regime described by an effective density matrix: the Floquet diagonal ensemble (FDE). By means of numerical examples on the Quantum Ising Chain and the Lipkin model, we discuss the properties of the FDE focusing on the difference among ergodic and regular quantum dynamics and on how this reflects on the thermal properties (T=∞) of the asymptotic condition. We verify thermalization in the classically ergodic Lipkin model and we demonstrate that this effect is induced by the Floquet states being delocalized and obeying Eigenstate Thermalization Hypothesis.We discuss also, in the Ising chain case, the work probability distribution, whose asymptotic condition is not described by the form (Generalized Gibbs Ensemble) that FDE acquires for local obserbvables because of integrability. Dephasing makes some correlations invisible in the local observables, but they are still present in the system. We consider also the linear response limit: when the amplitude of the perturbation is vanishingly small, the Floquet diagonal ensemble is not sufficient to describe the asymptotic condition given by LRT. For every small but finite amplitude, there are quasi-degeneracies in the Floquet spectrum giving rise to pre-relaxation to the condition predicted by Linear Response; these phenomena are strictly related to energy absorption and boundedness of the spectrum.
This review emphasizes the physical principles underlying the dephasing of nonstationary vibrational states of molecules with multiple vibrational degrees of freedom. The motion of atoms within molecules can be described to a very good approximation by a many-level system of coupled anharmonic quantized oscillators. The separation of electronic and nuclear timescales, combined with weak symmetry breaking inherent in molecular structures, guarantees that the Hamiltonian can always be cast in a form that is local in the quantum state space of the molecular vibrations. The resulting nonexponential dephasing dynamics can be controlled with external fields to stabilize nonstationary quantum states, and both quantum–classical and fully quantum control formalisms are described. Coupled molecular vibrations interacting with external fields also offer prospects for quantum computing because vibrational level spacings can be made very large compared to thermal noise, and relaxation mechanisms such as infrared fluorescence are many orders of magnitude slower than the timescales required for coherence transfer. Finally, a toy model that provides insights into the interaction of vibrational degrees of freedom with 'solvent' modes is also discussed, and exhibits nonexponential dynamics in certain regimes.
Bohm theory is a formulation of Quantum Mechanics that characterises the state of a quantum system according to both the wave function, as in the conventional formulation, and the coordinates (positions) of all the particles that evolve in time drawing quantum continuous trajectories. Furthermore, a statistical ensemble of all the possible trajectories, raising from the impossibility to know the initial position of all the particles, establishes the exact correspondence with the traditional Quantum Mechanics. From a computational point of view, Bohm theory has found many applications in Chemical Physics especially to develop new methodologies for solving the Schrödinger equation and to address semi-classical approximations of Quantum Mechanics. From a theoretical point of view, the most appealing feature of Bohm theory is its capability to supply a conceptual map between the quantum formalism and our representation of what a chemical system is. Chemical systems are composed of molecules, but the same idea of molecule requires a specific arrangement in the space of particles, i.e., the nuclei of the atoms. The statistical description of conventional Quantum Mechanics on the basis of wave function alone is insufficient to establish a clear correspondence with such a picture of molecules. Indeed, chemists employ usually Classical Mechanics in order to overcome this drawback of the standard quantum theory. On the other hand, if the particles position is included in the quantum formalism, as Bohm theory does, the map can be defined in a self-consistent way. In other words, Bohm theory appears to be the suitable quantum framework to represent molecules and their motion. The chemical representation of molecular systems finds a natural correspondence with a single Bohm trajectory, since it is always implicitly assumed that molecular components have specific spatial position independently of our knowledge about it. Consequently, we develop a quantum method whose fundamental assumption is that a single Bohm trajectory, i.e., a quantum molecular trajectory, describes the molecular systems and the molecular motion correctly. First of all, we examine the correspondence between a single Bohm trajectory and the conventional Quantum Mechanics, without using the ensemble of trajectories. We verify that such a correspondence exists through numerical simulations and we prove formally that the statistical properties of a single Bohm trajectory explain the probabilistic description of Quantum Mechanics. Once the consistency of this original approach has been established, we investigate the predicted properties. For instance, we take into account the constants of motion (such as the energy) corresponding to the time evolution of the coordinates and the behaviour of simple chemical systems, e.g., the vibrational motion of single molecules interacting with a resonant field. In this way, unexpected features of the molecular motion are found. Secondly, we tackle the challenge of describing many components systems (like the chemical systems in ordinary conditions). As a matter of fact, the computation of the Bohm trajectory and of the wave function is extremely demanding. However, the statistical properties of the Bohm trajectory allow the derivation of stochastic theories for examining the dynamics of open quantum systems, i.e., few molecules (or few degrees of freedom) interacting with their environment (the other molecules). One of the developed stochastic methods correlates the dynamics of the reduced density matrix, for the degrees of freedom of interest, to the evolution of the corresponding Bohm coordinates. In other words, the Bohm equation, determining the set of all the particles velocities according to the full wave function, is replaced with a stochastic one that approximates the velocity of a subset of coordinates according to the reduced density matrix. In such a way, the quantum fluctuations induced by the environment are taken into account. The advantage of this method concerns its capability of describing quantum systems, including open quantum systems, in terms of a quantum trajectory. This could allow the understanding of the molecular motion during a spectroscopical experiment. The possibility of investigating reactive systems, such as conformational changes, is particularly interesting. As a matter of fact, chemical reactions can be completely characterised only through the particles motion and we define the suit- able quantum methodology providing a self-consistent description of the molecular motion.
The formalism of collective (Dicke-like) states and operators is used to describe the relaxation of the initially prepared state of a material in the presence of a laser field. The material states are shown to be expressible as minimum uncertainty wave packets formed from those special collective states that couple with the radiation field. It is shown that these coherent states of matter exhibit all the features that are usually attributed to ensemble averages of two-level molecules in a conventional statistical picture. However, there are new features in the quantum mechanical model when there are correlations among the molecules in the ensemble. The loss of minimum uncertainty characteristics of the wave packet correspnds to the transverse (T2) relaxation of the statistical model. Different types of T2 processes can be defined corresponding to the decay of different orders of coherence of the wave packet. Only for the case that the system decays in all orders of coherence do the quantum and statistical results coincide. There is expected to be a transverse relaxation that preserves one or more orders of coherence. Experiments are suggested that can test these ideas. A particular case of the formalism of collective states applies to molecular crystals and it provides a quantum mechanical picture of the dephasing of molecular excitons in terms of orders of coherence analogous to the usual manner of describing states of the radiation field.
We investigate the effect of inter-mode coupling on the vibrational relaxation dynamics of molecules in weak dissipative environments. The simulations are performed within the reduced density matrix formalism in the Markovian regime, assuming a Lindblad form for the system-bath interaction. The prototypical two-dimensional model system representing two CO molecules approaching a Cu(100) surface is adapted from an ab initio potential, while the diatom-diatom vibrational coupling strength is systematically varied. In the weak system-bath coupling limit and at low temperatures, only first order non-adiabatic uni-modal coupling terms contribute to surface-mediated vibrational relaxation. Since dissipative dynamics is non-unitary, the choice of representation will affect the evolution of the reduced density matrix. Two alternative representations for computing the relaxation rates and the associated operators are thus compared: the fully coupled spectral basis, and a factorizable ansatz. The former is well-established and serves as a benchmark for the solution of Liouville-von Neumann equation. In the latter, a contracted grid basis of potential-optimized discrete variable representation is tailored to incorporate most of the inter-mode coupling, while the Lindblad operators are represented as tensor products of one-dimensional operators, for consistency. This procedure results in a marked reduction of the grid size and in a much more advantageous scaling of the computational cost with respect to the increase of the dimensionality of the system. The factorizable method is found to provide an accurate description of the dissipative quantum dynamics of the model system, specifically of the time evolution of the state populations and of the probability density distribution of the molecular wave packet. The influence of intra-molecular vibrational energy redistribution appears to be properly taken into account by the new model on the whole range of coupling strengths. It demontrates that most of the mode mixing during relaxation is due to the potential part of the Hamiltonian and not to the coupling among relaxation operators.