Source author record

Yoshitaka Tanimura

Yoshitaka Tanimura appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

23works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

23 published item(s)

preprint2022arXiv

Numerically "exact" simulations of a quantum Carnot cycle: Analysis using thermodynamic work diagrams

We investigate the efficiency of a quantum Carnot engine based on open quantum dynamics theory. The model includes time-dependent external fields for the subsystems controlling the isothermal and isentropic processes and for the system--bath (SB) interactions controlling the transition between these processes. Numerical simulations are conducted in a nonperturbative and non-Markovian SB coupling regime using the hierarchical equations of motion under these fields at different cycle frequencies. The work applied to the total system and the heat exchanged with the baths are rigorously evaluated. In addition, by regarding quasi-static work as free energy, we compute the quantum thermodynamic variables and analyze the simulation results using thermodynamic work diagrams for the first time. Analysis of these diagrams indicates that, in the strong SB coupling region, the fields for the SB interactions are major sources of work, while in other regions, the field for the subsystem is a source of work. We find that the maximum efficiency is achieved in the quasi-static case and is determined solely by the bath temperatures, regardless of the SB coupling strength, which is a numerical manifestation of Carnot's theorem.

preprint2022arXiv

The laws of thermodynamics for quantum dissipative systems: A quasi-equilibrium Helmholtz energy approach

Using the quasi-equilibrium Helmholtz energy (qHE), defined as the thermodynamic work in a quasi-static process, we investigate the thermal properties of both an isothermal process and a transition process between the adiabatic and isothermal states (adiabatic transition). Here, the work is defined by the change in energy from a steady-state to another state under a time-dependent perturbation. In particular, the work for a quasi-static change is regarded as thermodynamic work. We employ a system--bath model that involves time-dependent perturbations in both the system and the system--bath interaction. We conduct numerical experiments for a three-stroke heat machine (a Kelvin-Planck cycle). For this purpose, we employ the hierarchical equations of motion (HEOM) approach. These experiments involve an adiabatic transition field that describes the operation of an adiabatic wall between the system and the bath. Thermodynamic--work diagrams for external fields and their conjugate variables, similar to the $P$--$V$ diagram, are introduced to analyze the work done for the system in the cycle. We find that the thermodynamic efficiency of this machine is zero because the field for the isothermal processes acts as a refrigerator, whereas that for the adiabatic wall acts as a heat engine. This is a numerical manifestation of the Kelvin--Planck statement, which states that it is impossible to derive mechanical effects from a single heat source. These HEOM simulations serve as a rigorous test of thermodynamic formulations because the second law of thermodynamics is only valid when the work involved in the operation of the adiabatic wall is treated accurately.

preprint2021arXiv

Exciton Transfer in Organic Photovoltaic Cells: A Role of Local and Nonlocal Electron-Phonon Interactions in a Donor Domain

We theoretically investigate an exciton transfer process in a donor domain of organic photovoltaic cells focusing on the roles of local and nonlocal electron-phonon interactions. Our model consists of a three-level system described by the Holstein-Peierls Hamiltonian coupled to multiple heat baths for local and nonlocal molecular modes characterized by Brownian spectral distribution functions. We chose tetracene as a reference donor molecule, where the spectral distribution functions of the local and nonlocal modes exist. We then employ the reduced hierarchy equations of motion (HEOM) approach to simulating the dynamics of the system under the influence of the environment as a function of the electron-phonon coupling strength and temperature. We rigorously calculate the reduced density matrix elements to explain the timescale of dynamics under the influence of the dissipative local and nonlocal modes. The results indicate that the strong nonlocal electron-phonon interaction under high temperature conditions favors the exciton transfer process and enhances the efficiency of organic photovoltaic materials, while the lifetime of the exciton becomes shorter due to a low frequency local mode.

preprint2021arXiv

Open quantum dynamics theory for a complex subenvironment system with a quantum thermostat: Application to a spin heat bath

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 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 of the anisotropy and the ferromagnetic or antiferromagnetic orders of the spins.

preprint2021arXiv

Open Quantum Dynamics Theory for Non-Equilibrium Work: Hierarchical Equations of Motion Approach

A system--bath (SB) model is considered to examine the Jarzynski equality in the fully quantum regime. In our previous paper [J. Chem. Phys. 153 (2020) 234107], we carried out "exact" numerical experiments using hierarchical equations of motion (HEOM) in which we demonstrated that the SB model describes behavior that is consistent with the first and second laws of thermodynamics and that the dynamics of the total system are time irreversible. The distinctive quantity in the Jarzynski equality is the "work characteristic function (WCF)", $\langle \exp[-βW] \rangle$, where $W$ is the work performed on the system and $β$ is the inverse temperature. In the present investigation, we consider the definitions based on the partition function (PF) and on the path, and numerically evaluate the WCF using the HEOM to determine a method for extending the Jarzynski equality to the fully quantum regime. We show that using the PF-based definition of the WCF, we obtain a result that is entirely inconsistent with the Jarzynski equality, while if we use the path-based definition, we obtain a result that approximates the Jarzynski equality, but may not be consistent with it.

preprint2020arXiv

Modeling intermolecular and intramolecular modes of liquid water using multiple heat baths: Machine learning approach

The vibrational motion of molecules in dissipative environments, such as solvation and protein molecules, is composed of contributions from both intermolecular and intramolecular modes. The existence of these collective modes introduces difficulty into quantum simulations of chemical and biological processes. In order to describe the complex molecular motion of the environment in a simple manner, we introduce a system-bath model in which the intramolecular modes with anharmonic mode-mode couplings are described by a system Hamiltonian, while the other degrees of freedom, arising from the environmental molecules, are described by heat bath. Employing a machine-learning based approach, we determine not only the system parameters of the intramolecular modes but also the spectral distribution of the system-bath coupling to describe the intermolecular modes, using the atomic trajectories obtained from molecular dynamics (MD) simulations. The capabilities of the present approach are demonstrated for liquid water using MD trajectories calculated from the SPC/E model and the polarizable water model for intramolecular and intermolecular vibrational spectroscopies (POLI2VS) by determining the system parameters describing the symmetric-stretch, asymmetric-stretch and bend modes with intramolecular interactions and the bath spectral distribution functions for each intramolecular mode representing the interaction with the intra-molecular modes. From these results, we were able to elucidate the energy relaxation pathway between the intramolecular modes and the intermolecular modes in a non-intuitive manner.

preprint2020arXiv

Open Quantum Dynamics Theory of Spin Relaxation: Application to $μ$SR and Low-Field NMR Spectroscopies

An open quantum system refers to a system, which is in turn 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 an 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, to help further the development of these experimental means.

preprint2020arXiv

Perspective: Numerically "exact" approach to open quantum dynamics: The hierarchical equations of motion (HEOM)

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 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.

preprint2020arXiv

Proton Tunneling in a Two-Dimensional Potential Energy Surface with a Non-linear System-Bath Interaction: Thermal Suppression of Reaction Rate

We consider a proton-transfer (PT) system described by a proton-transfer reaction (PTR) coordinate and a rate promoting vibrational (RPV) coordinate interacting with a non-Markovian heat-bath. While dynamics of PT processes has been widely discussed using two-dimensional (2D) potential energy surfaces (PES), the role of the heat-bath, in particular, for a realistic form of the system-bath interaction has not been well explored. Previous studies are largely based on one-dimensional model and linear-linear (LL) system-bath interaction. In the present study, we introduce an exponential-linear (EL) system-bath interaction, which is derived from the analysis of a PTR-PRV system in a realistic situation. This interaction mainly causes vibrational dephasing in the PTR mode and population relaxation in the RPV mode. Numerical simulations were carried out using hierarchy equations of motion approach. We analyze the role of the heat-bath interaction on the chemical reaction rate as a function of the system-bath coupling strength at different temperature and for different values of the bath correlation time. A prominent feature of the present result is that while the reaction rate predicted from classical and quantum Kramers theory increases as the temperature increases, the present EL interaction model exhibits opposite temperature dependence. Kramers turn-over profile of the reaction rate as a function of the system-bath coupling is also suppressed in the present EL model turning into a plateau-like curve for larger system-bath interaction strength. Such features arise from the interplay of the vibrational dephasing process in the PTR mode and the population relaxation process in the RPV mode.

preprint2016arXiv

Exploring complete positivity in hierarchy equations of motion

We derive a purely algebraic framework for the identification of hierarchy equations of motion that induce completely positive dynamics and demonstrate the applicability of our approach with several examples. We find bounds on the violation of complete positivity for microscopically derived hierarchy equations of motion and construct well-behaved phenomenological models with strongly non-Markovian revivals of quantum coherence.

preprint2016arXiv

Full molecular dynamics simulations of liquid water and carbon tetrachloride for two-dimensional Raman spectroscopy in the frequency domain

Frequency-domain two-dimensional Raman signals, which are equivalent to coherent two-dimensional Raman scattering (COTRAS) signals, for liquid water and carbon tetrachloride were calculated using an equilibrium-nonequilibrium hybrid MD simulation algorithm. We elucidate mechanisms governing the 2D signal profiles involving anharmonic mode-mode coupling and the nonlinearities of the polarizability for the intermolecular and intramolecular vibrational modes. The predicted signal profiles and intensities can be utilized to analyze recently developed single-beam 2D spectra, whose signals are generated from a coherently controlled pulse, allowing the single-beam measurement to be carried out more efficiently.

preprint2015arXiv

Analysis of 2D THz-Raman spectroscopy using a non-Markovian Brownian oscillator model with nonlinear system-bath interactions

We explore and describe the roles of inter-molecular vibrations employing a Brownian oscillator (BO) model with linear-linear (LL) and square-linear (SL) system-bath interactions, which we use to analyze two-dimensional (2D) THz-Raman spectra obtained by means of molecular dynamics (MD) simulations. In addition to linear absorption (1D IR), we calculated 2D Raman-THz-THz, THz-Raman-THz, and THz-THz-Raman signals for liquid formamide, water, and methanol using an equilibrium non-equilibrium hybrid MD simulation. The calculated 1D IR and 2D THz-Raman signals are compared with results obtained from the LL+SL BO model applied through use of hierarchal Fokker-Planck equations with non-perturbative and non-Markovian noise. We find that all of the qualitative features of the 2D profiles of the signals obtained from the MD simulations are reproduced with the LL+SL BO model, indicating that this model captures the essential features of the inter-molecular motion. We analyze the fitted 2D profiles in terms of anharmonicity, nonlinear polarizability, and dephasing time. The origins of the echo peaks of the librational motion and the elongated peaks parallel to the probe direction are elucidated using optical Liouville paths.

preprint2015arXiv

Linear, third- and fifth-order nonlinear spectroscopy of a charge transfer system coupled to an underdamped vibration

We study hole, electron and exciton transport in a charge transfer system in the presence of underdamped vibrational motion. We analyze the signature of these processes in the linear and third-, and fifth-order nonlinear electronic spectra. Calculations are performed with a numerically exact hierarchical equations of motion method for an underdamped Brownian oscillator spectral density. We find that combining electron, hole and exciton transfer can lead to non-trivial spectra with more structure than with excitonic coupling alone. Traces taken during the waiting time of a two-dimensional spectrum are dominated by vibrational motion and do not reflect the electron, hole, and exciton dynamics directly. We find that the fifth-order nonlinear response is particularly sensitive to the charge transfer process. While third-order 2D spectroscopy detects the correlation between two coherences, fifth-order 2D spectroscopy (2D population spectroscopy) is here designed to detect correlations between the excited states during two different time periods.

preprint2015arXiv

Real-Time and Imaginary-Time Quantum Hierarchal Fokker-Planck Equations

We consider a quantum mechanical system represented in phase space (referred to hereafter as "Wigner space"), coupled to a harmonic oscillator bath. We derive quantum hierarchal Fokker-Planck (QHFP) equations not only in real time, but also in imaginary time, which represents an inverse temperature. This is an extension of a previous work, in which we studied a spin-boson system, to a Brownian system. It is shown that the QHFP in real time obtained from a correlated thermal equilibrium state of the total system possess the same form as those obtained from a factorized initial state. A modified terminator for the hierarchal equations of motion is introduced to treat the non-Markovian case more efficiently. Using the imaginary-time QHFP, numerous thermodynamic quantities, including the free energy, entropy, internal energy, heat capacity, and susceptibility can be evaluated for any potential. These equations allow us to treat non-Markovian, non-perturbative system-bath interactions at finite temperature. Through numerical integration of the real-time QHFP for a harmonic system, we obtain the equilibrium distributions, the auto-correlation function, and the first- and second-order response functions. These results are compared with analytically exact results for the same quantities. This provides a critical test of the formalism for a non-factorized thermal state, and elucidates the roles of fluctuation, dissipation, non-Markovian effects, and system-bath coherence. Employing numerical solutions of the imaginary-time QHFP, we demonstrate the capability of this method to obtain thermodynamic quantities for any potential surface. It is shown that both types of QHFP equations can produce numerical results of any desired accuracy.

preprint2015arXiv

Reduced hierarchical equations of motion in real and imaginary time: Correlated initial states and thermodynamic quantities

For a system strongly coupled to a heat bath, the quantum coherence of the system and the heat bath plays an important role in the system dynamics. This is particularly true in the case of non-Markovian noise. We rigorously investigate the influence of system-bath coherence by deriving the reduced hierarchal equations of motion (HEOM), not only in real time, but also in imaginary time, which represents an inverse temperature. It is shown that the HEOM in real time obtained when we include the system-bath coherence of the initial thermal equilibrium state possess the same form as those obtained from a factorized initial state. We find that the difference in behavior of systems treated in these two manners results from the difference in initial conditions of the HEOM elements, which are defined in path integral form. We also derive HEOM along the imaginary time path to obtain the thermal equilibrium state of a system strongly coupled to a non-Markovian bath. Then, we show that the steady state hierarchy elements calculated from the real-time HEOM can be expressed in terms of the hierarchy elements calculated from the imaginary-time HEOM. Moreover, we find that the imaginary-time HEOM allow us to evaluate a number of thermodynamic variables, including the free energy, entropy, internal energy, heat capacity, and susceptibility. The expectation values of the system energy and system-bath interaction energy in the thermal equilibrium state are also evaluated.

preprint2013arXiv

An approach to quantum transport based on reduced hierarchy equations of motion: Application to a resonant tunneling diode

The quantum dissipative dynamics of a tunneling process through double barrier structures is investigated on the basis of a rigorous treatment for the first time. We employ a Caldeira-Leggett Hamiltonian with an effective potential calculated self-consistently, accounting for the electron distribution. With this Hamiltonian, we use the reduced hierarchy equations of motion in the Wigner space representation to study the effects of non-Markovian thermal fluctuations and dissipation at finite temperature in a rigorous manner. Hysteresis, double plateau-like behavior, and self-excited current oscillation are observed in a negative differential resistance (NDR) region of the current-voltage curve. We find that while most of the current oscillations decay in time in the NDR region, there is a steady oscillation characterized by a tornado-like rotation in the Wigner space in the upper plateau of the NDR region.

preprint2013arXiv

Self-excited current oscillations in a resonant tunneling diode described by a model based on the Caldeira-Leggett Hamiltonian

The quantum dissipative dynamics of a tunneling process through double barrier structures is investigated on the basis of non-perturbative and non-Markovian treatment. We employ a Caldeira-Leggett Hamiltonian with an effective potential calculated self-consistently, accounting for the electron distribution. With this Hamiltonian, we use the reduced hierarchy equations of motion in the Wigner space representation to study non-Markovian and non-perturbative thermal effects at finite temperature in a rigorous manner. We study current variation in time and the current-voltage (I-V) relation of the resonant tunneling diode for several widths of the contact region, which consists of doped GaAs. Hysteresis and both single and double plateau-like behavior are observed in the negative differential resistance (NDR) region. While all of the current oscillations decay in time in the NDR region in the case of a strong system-bath coupling, there exist self-excited high-frequency current oscillations in some parts of the plateau in the NDR region in the case of weak coupling. We find that the effective potential in the oscillating case possesses a basin-like form on the emitter side (emitter basin) and that the current oscillation results from tunneling between the emitter basin and the quantum well in the barriers. We find two distinct types of current oscillations, with large and small oscillation amplitudes, respectively. These two types of oscillation appear differently in the Wigner space, with one exhibiting tornado-like motion and the other exhibiting a two piston engine-like motion.

preprint2012arXiv

Dephasing by a Continuous-Time Random Walk Process

Stochastic treatments of magnetic resonance spectroscopy and optical spectroscopy require evaluations of functions like <exp(i int_0^t Q(s)ds)>, where t is time, Q(s) is the value of a stochastic process at time s, and the angular brackets denote ensemble averaging. This paper gives an exact evaluation of these functions for the case where Q is a continuous-time random walk process. The continuous time random walk describes an environment that undergoes slow, step-like changes in time. It also has a well-defined Gaussian limit, and so allows for non-Gaussian and Gaussian stochastic dynamics to be studied within a single framework. We apply the results to extract qubit-lattice interaction parameters from dephasing data of P-doped Si semiconductors (data collected elsewhere), and to calculate the two-dimensional spectrum of a three level harmonic oscillator undergoing random frequency modulations.

preprint2012arXiv

Reduced hierarchy equations of motion approach with Drude plus Brownian spectral distribution: Probing electron transfer processes by means of two- dimensionalcorrelation spectroscopy

We theoretically investigate an electron transfer (ET) process in a dissipative environment by means of two-dimensional (2D) correlation spectroscopy. We extend the reduced hierarchy equations of motion approach to include both overdamped Drude and underdamped Brownian modes. While the overdamped mode describes the inhomogeneity of a system in the slow modulation limit, the underdamped mode expresses the primary vibrational mode coupled with the electronic states. We outline a procedure for calculating 2D correlation spectrum that incorporates the ET processes. The present approach has the capability of dealing with system-bath coherence under an external perturbation, which is important to calculate nonlinear response functions for non-Markovian noise. The calculated 2D spectrum exhibits the effects of the ET processes through the presence of ET transition peaks along the $Ω_1$ axis, as well as the decay of echo signals.

preprint2011arXiv

Non-Markovianity: initial correlations and nonlinear optical measurements

By extending the response function approach developed in nonlinear optics, we analytically derive an expression for the non-Markovianity [Laine, et al, Phys. Rev. A 81, 062115 (2010)] in the time evolution of a system in contact with a quantum mechanical bath, and find a close connection with the directly observable nonlinear optical response. The result indicates that memory in the bath induced fluctuations rather than in the dissipation causes non-Markovianity. Initial correlations between states of the system and the bath are shown to be essential for a correct understanding of the non-Markovianity.

preprint2010arXiv

Correlated fluctuations in the exciton dynamics and spectroscopy of DNA

The absorption of ultraviolet light creates excitations in DNA, which subsequently start moving in the helix. Their fate is important for an understanding of photo damage, and is determined by the interplay of electronic couplings between bases and the structure of the DNA environment. We model the effect of dynamical fluctuations in the environment and study correlation, which is present when multiple base pairs interact with the same mode in the environment. We find that the correlations strongly affect the exciton dynamics, and show how they are observed in the decay of the anisotropy as a function of a coherence and a population time in a non-linear optical experiment.

preprint2010arXiv

Non-Markovian entanglement dynamics in the presence of system-bath coherence

A complete treatment of the entanglement of two-level systems, which evolves through the contact with a thermal bath, must include the fact that the system and the bath are not fully separable. Therefore, quantum coherent superpositions of system and bath states, which are almost never fully included in theoretical models, are invariably present when an entangled state is prepared experimentally. We show their importance for the time evolution of the entanglement of two qubits coupled to independent baths. In addition, our treatment is able to handle slow and low-temperature thermal baths.

preprint2002arXiv

Energy-level diagrams and their contribution to fifth-order Raman and second-order infrared responses: Distinction between relaxation mechanisms by two-dimensional spectroscopy

We develop a Feynman rule for energy-level diagrams emphasizing their connections to the double-sided Feynman diagrams and physical processes in the Liouville space. Thereby we completely identify such diagrams and processes contributing to the two-dimensional response function in the Brownian oscillator model. We classify such diagrams or processes in quartet and numerically present signals separately from each quartet of diagrams or Liouville-space processes. We find that the signal from each quartet is distinctly different from the others; we can identify each peaks in frequency domain with a certain quartet. This offers the basis for analyzing and assigning actual two-dimensional peaks and suggests the possibility of Liouville-space-path selective spectroscopy. As an application we demonstrate an example in which two familiar homogeneous mechanisms of relaxation are distinguished by existence or non-existence of certain peaks on the two-dimensional map; appearance or disappearance of certain peak is sensitive to the coupling mechanism. We also point out some confusion in the literature with regard to inclusion of relaxation effects.