An open quantum system is a system coupled to an environment that can describe time-irreversible dynamics through which the system evolves toward the thermal equilibrium state.We present a quantum mechanically rigorous theory in order to help the analysis of spectra obtained from the advanced nuclear magnetic resonance (NMR) and muon spin rotation, relaxation or resonance (μSR) techniques.Our approach is based on the numerically "Exact" hierarchical equations of motion (HEOM) approach, which allows us to study the reduced system dynamics for non-perturbative and non-Markovian system-bath interactions at finite temperature even under strong time-dependent perturbations.We demonstrate the present theory to analyze μSR and low-field NMR spectra, as an extension of the Kubo-Toyabe theory focusing on the effects of temperature and anisotropy of a local magnetic field on spectra, to help further the development of these experimental means.
An open quantum system is a system coupled to an environment that can describe time-irreversible dynamics through which the system evolves toward the thermal equilibrium state.We present a quantum mechanically rigorous theory in order to help the analysis of spectra obtained from the advanced nuclear magnetic resonance (NMR) and muon spin rotation, relaxation or resonance (μSR) techniques.Our approach is based on the numerically "Exact" hierarchical equations of motion (HEOM) approach, which allows us to study the reduced system dynamics for non-perturbative and non-Markovian system-bath interactions at finite temperature even under strong time-dependent perturbations.We demonstrate the present theory to analyze μSR and low-field NMR spectra, as an extension of the Kubo-Toyabe theory focusing on the effects of temperature and anisotropy of a local magnetic field on spectra, to help further the development of these experimental means.
An open quantum system refers to a system that is further coupled to a bath system consisting of surrounding radiation fields, atoms, molecules, or proteins. The bath system is typically modeled by an infinite number of harmonic oscillators. This system-bath model can describe the time-irreversible dynamics through which the system evolves toward a thermal equilibrium state at finite temperature. In nuclear magnetic resonance and atomic spectroscopy, dynamics can be studied easily by using simple quantum master equations under the assumption that the system-bath interaction is weak (perturbative approximation) and the bath fluctuations are very fast (Markovian approximation). However, such approximations cannot be applied in chemical physics and biochemical physics problems, where environmental materials are complex and strongly coupled with environments. The hierarchical equations of motion (HEOM) can describe the numerically "exact" dynamics of a reduced system under nonperturbative and non-Markovian system-bath interactions, which has been verified on the basis of exact analytical solutions (non-Markovian tests) with any desired numerical accuracy. The HEOM theory has been used to treat systems of practical interest, in particular, to account for various linear and nonlinear spectra in molecular and solid state materials, to evaluate charge and exciton transfer rates in biological systems, to simulate resonant tunneling and quantum ratchet processes in nanodevices, and to explore quantum entanglement states in quantum information theories. This article presents an overview of the HEOM theory, focusing on its theoretical background and applications, to help further the development of the study of open quantum dynamics.
We investigate the quantum dynamics of Coulomb potential systems in thermal baths. We study these systems within the framework of open quantum dynamics theory, focusing on preserving the rotational symmetry of the entire system, including the baths. Thus, we employ a three-dimensional rotationally invariant system-bath model to derive numerically "exact" hierarchical equations of motion for atomic orbitals that enable a non-perturbative and non-Markovian treatment of system-bath interactions at finite temperatures. To assess the formalism, we calculated the linear absorption spectrum of an atomic system under an isotropic thermal environment, with systematic variation of system-bath coupling strength and temperature.
In this study, we provide a way to describe the dynamics of quantum tomography in an open system with a generalized master equation, considering a case where the relevant system under tomographic measurement is influenced by the environment. We apply this to spin tomography because such situations typically occur in μSR (muon spin rotation/relaxation/resonance) experiments where microscopic features of the material are investigated by injecting muons as probes. As a typical example to describe the interaction between muons and a sample material, we use a spin-boson model where the relevant spin interacts with a bosonic environment. We describe the dynamics of a spin tomogram using a time-convolutionless type of generalized master equation that enables us to describe short time scales and/or low-temperature regions. Through numerical evaluation for the case of Ohmic spectral density with an exponential cutoff, a clear interdependency is found between the time evolution of elements of the density operator and a spin tomogram. The formulation in this paper may provide important fundamental information for the analysis of results from, for example, μSR experiments on short time scales and/or in low-temperature regions using spin tomography.
Explanation of decoherence and quasi-equilibrium in systems with few degrees of freedom demands a deep theoretical analysis that considers the observed system as an open quantum system. In this work, we study the problem of decoherence of an observed system of quantum interacting particles, coupled to a quantum lattice. Our strategy is based on treating the environment and the system-environment Hamiltonians fully quantum mechanically, which yields a representation of the time evolution operator useful for disentangling the different time scales underlying in the observed system dynamics. To describe the possible different stages of the dynamics of the observed system, we introduce quantum mechanical definitions of essentially isolated, essentially adiabatic, and thermal-contact system-environment interactions. This general approach is then applied to the study of decoherence and quasi-equilibrium in proton nuclear magnetic resonance ((1)H NMR) of nematic liquid crystals. A summary of the original results of this work is as follows. We calculate the decoherence function and apply it to describe the evolution of a coherent spin state, induced by the coupling with the molecular environment, in absence of spin-lattice relaxation. By assuming quantum energy conserving or non-demolition interactions, we identify an intermediate time scale, between those controlled by self-interactions and thermalization, where coherence decays irreversibly. This treatment is also adequate for explaining the buildup of quasi-equilibrium of the proton spin system, via the process we called eigen-selectivity. By analyzing a hypothetical time reversal experiment, we identify two sources of coherence loss which are of a very different nature and give rise to distinct time scales of the spin dynamics: (a) reversible or adiabatic quantum decoherence and (b) irreversible or essentially adiabatic quantum decoherence. Local irreversibility arises as a consequence of the uncertainty introduced by the coupling with an infinite quantum environment. The reversible part can be represented by a semiclassical model, similar to standard line-shape adiabatic models. By exploiting the separation existing between the time scales of the spin coherences and the irreversible decoherence, we present a novel technique to obtain the orientational molecular distribution function for a nematic liquid crystal. The procedure is based on the comparison of the observed coherence time evolution and numerical calculation under the adiabatic quantum decoherence approach. As an example, it is used the experimental free induction decay from a nematic PAA(d6) sample to extract such an orientational distribution. This is the first theoretical description of the experimental liquid crystal NMR signal in the time domain. On the contrary, the irreversible decoherence is intrinsically full-quantum mechanical, as it is governed by the commutation properties of the environment and the spin-lattice Hamiltonians. Consistently, it depends on the molecular correlation in a decisive way, since it vanishes under a mean-field model for the molecular dynamics. The results of this work can contribute to the understanding of the open question of the applicability of the spin-temperature concept in spin systems with few degrees of freedom.
Complex environments, such as molecular matrices and biological material, play a fundamental role in many important dynamic processes in condensed phases. Because it is extremely difficult to conduct full quantum dynamics simulations on such environments due to their many degrees of freedom, here, we treat in detail the environment only around the main system of interest (the subenvironment), while the other degrees of freedom needed to maintain the equilibrium temperature are described by a simple harmonic bath, which we call a quantum thermostat. The noise generated by the subenvironment is spatially non-local and non-Gaussian and cannot be characterized by the fluctuation-dissipation theorem. We describe this model by simulating the dynamics of a two-level system (TLS) that interacts with a subenvironment consisting of a one-dimensional XXZ spin chain. The hierarchical Schrödinger equations of motion are employed to describe the quantum thermostat, allowing for time-irreversible simulations of the dynamics at arbitrary temperature. To see the effects of a quantum phase transition of the subenvironment, we investigate the decoherence and relaxation processes of the TLS at zero and finite temperatures for various values of the spin anisotropy. We observed the decoherence of the TLS at finite temperature even when the anisotropy of the XXZ model is enormous. We also found that the population-relaxation dynamics of the TLS changed in a complex manner with the change in the anisotropy and the ferromagnetic or antiferromagnetic orders of spins.
This thesis is concerned with the design of numerical methods for quantum simulation and the development of improved models for quantum relaxation. Analysis is presented for the treatment of quantum systems using the density matrix formalism. This approach has been developed from the early days of quantum mechanics as a tool to describe from a statistical point of view a large number of identical quantum ensembles. Traditional methods are well established and reliable, but they perform poorly for practical simulation as the system size is scaled up. Ad hoc schemes for nuclear spin dynamics appearing in the literature can be shown to fail in certain situations. The challenge is therefore to identify efficient reduction methods for the quantum system which are also based on a rigorous foundation. The method presented in the thesis, for the time–independent Hamiltonian case, combines a quantum density matrix formalism with a procedure based on Chebyshev polynomials; application of the method to Nuclear Magnetic Resonance (NMR) spectroscopy is considered, and it is shown that the new technique outperforms existing alternatives in term of computational costs. The case of a time–dependent Hamiltonian in NMR simulation is studied as well and some splitting methods are presented. To the author’s knowledge this is the first time such methods have been applied within the NMR framework, and the numerical results show a better error–to–cost rate than traditional methods. In a separate strand of research, formulations for open quantum systems are studied and new dynamical systems approaches are considered for this problem. Motivations This thesis work is mainly focused on nuclear spin dynamics. Nuclear spin dynamics constitutes the basis for NMR, which is a very powerful spectroscopy technique that exploits the interaction between nuclear spins and magnetic fields. The same technique is used to reveal the presence of hydrogen atoms in the blood for Magnetic Resonance Imaging (MRI). Within this framework the role of simulations is extremely important, as it provides a benchmark for studies of new materials, and the development of new magnetic fields. The main computational issue is that with current software for NMR simulation it is extremely expensive to deal with systems made of more than few (7–10) spins. There is therefore a strong need to develop new algorithms capable of simulating larger systems. In recent years NMR simulations have been found to be one of the most favorable candidates for quantum computing. There are two reasons for this: nuclear quantum states maintain extremely long coherences, and it is possible to attain a very strong control on the quantum state via the application of sequences of pulses. In order to develop a proper quantum computer it is fundamental to understand how the entangled states lose coherence and relax back to equilibrium by means of external interactions. This process is described as relaxation in an open quantum system. The theory for such systems has been available for 50 years but there are still substantial limitations in the two main approaches. There are also relatively few numerical approaches for the simulation of such systems, for this reason it is important to develop numerical alternatives for the description of open quantum systems. Thesis Outline The thesis is organized as follow: the first two chapters provide background material to familiarize the reader with fundamental concepts of both quantum mechanics and nuclear spin dynamics; in this part of the thesis no new results are presented. The first chapter introduces the concept of quantum systems and the mathematical environment with which we describe those systems. We also present the main equations we need to solve to determine the dynamics of a quantum system in a statistical framework. In the second chapter we introduce the nuclear spin system, that is the physical system that has been the main reference frame in this work, for both tests and practical applications of the new algorithms. We describe how nuclear spin systems are at the basis of very important applications like NMR spectroscopy and MRI. We present in some detail the physical features of the NMR technique and the equations we need to solve to describe the dynamics of a spin system; we also focus on the relevance of numerical simulations for these systems, and consequently which must be the interest in developing new algorithms, and the major obstacles which must be overcome. In the third chapter we investigate the numerical challenges that arise in simulation of quantum systems, we describe some of the methods that have been developed in the literature, focusing on the performances and the computational costs of them, setting the new developments of this thesis in the proper research frame. We discuss one of the major issues: the evaluation of the matrix exponential. We also present the analysis we have done of a recent method called Zero Track Elimination (ZTE) that has been developed specifically for NMR simulations. This analysis shows the limitations of this method but also gives a mathematical explanation of why–and in which cases–it works. In the fourth chapter we present the main result of the thesis, the development of a new method that directly evaluates the expectation values for a quantum simulation via a different application of the well known Chebyshev expansion. We have proved that this new method can provide an excellent boost in terms of performance, with computational costs that can be reduced by a factor ten in common cases. (The results of this chapter and the new method have been presented in international conferences and recently they have been submitted for publication). We also present some attempts we have made in the application of splitting methods for the evolution of the system in a time dependent environment. To our knowledge this is the first time splitting methods have been used for NMR simulations. The results of this approach are as follows: for a particular splitting technique combined with a Lanczos iteration method it is possible to speed up the calculation by a third if compared with a Lanczos type method whilst keeping the error below a critical threshold. This last approach is still a work in progress especially in terms of developing clever ways to split the Hamiltonian. The last chapter of this thesis deals with simulation of quantum systems interacting with an external environment. After presenting the main theoretical approaches for the description of such systems we then survey several the techniques that are currently used for the numerical implementation of such theories. As a work in progress we present a considerably different new approach we have been developing aiming to overcome some of the issues that arise when treating this kind of system within usual frameworks. This is somewhat speculative work that gives rise to some new directions in the development of a numerical description for open quantum systems. We also present some numerical results. (The main core of this chapter has been presented in international conferences).