Source author record

Jianshu Cao

Jianshu Cao 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

42works
11topics
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

42 published item(s)

preprint2023arXiv

Coherent spatial control of wave packet dynamics on quantum lattices

Quantum lattices are pivotal in the burgeoning fields of quantum materials and information science. Rapid developments in microscopy and quantum engineering allow for preparing and monitoring wave-packet dynamics on quantum lattices with increasing spatial and temporal resolution. Motivated by these emerging research interests, we present an analytical study of wave packet diffusivity and diffusion length on tight-binding quantum lattices subject to stochastic noise. Our analysis points to the crucial role of spatial coherence and predicts a set of novel phenomena: noise can enhance the transient diffusivity and diffusion length of sufficiently extended initial states; A smooth Gaussian initial state spreads slower than a localized initial state; A standing or traveling initial state with large momentum spreads faster than a localized initial state and exhibits a noise-induced peak in the transient diffusivity; The change in the time-dependent diffusivity and diffusion length relative to a localized initial state follows a universal dependence on the Gaussian width. These theoretical predictions and the underlying mechanism of spatial coherence suggest the possibility of controlling the wave packet dynamics on quantum lattices by spatial manipulations, which will have implications for materials science and quantum technologies.

preprint2022arXiv

Dynamical symmetries of periodically-driven quantum systems and their spectroscopic signatures

Spatial symmetries of quantum systems leads to important effects in spectroscopy, such as selection rules and dark states. Motivated by the increasing strength of light-matter interaction achieved in recent experiments, we investigate a set of dynamically-generalized symmetries for quantum systems, which are subject to a strong periodic driving. Based on Floquet response theory, we study rotational, particle-hole, chiral and time-reversal symmetries and their signatures in spectroscopy, including symmetry-protected dark states (spDS), a Floquet band selection rule (FBSR), and symmetry-induced transparency (siT). Specifically, a dynamical rotational symmetry establishes dark state conditions, as well as selection rules for inelastic light scattering processes; a particle-hole symmetry introduces dark states for symmetry related Floquet states and also a transparency effect at quasienergy crossings; chiral symmetry and time-reversal symmetry alone do not imply dark state conditions, but can be combined to the particle-hole symmetry. Our predictions reveal new physical phenomena when a quantum system reaches the strong light-matter coupling regime, important for superconducting qubits, atoms and molecules in optical or plasmonic field cavities, and optomechanical systems.

preprint2022arXiv

Generalized resonance energy transfer theory: Applications to vibrational energy flow in optical cavities

A general rate theory for resonance energy transfer is formulated to incorporate any degrees of freedom (e.g., rotation, vibration, exciton, and polariton) as well as coherently-coupled composite states. The compact rate expression allows us to establish useful relationships: (i) detailed balance condition when the donor and acceptor are at the same temperature; (ii) proportionality to the overlap between donor's emission and acceptor's absorption spectra; (iii) scaling with the effective coherent size, i.e., the number of coherently coupled molecules; (iv) spatial and orientational dependences as derived from the interaction potential. When applied to cavity-assisted vibrational energy transfer, the rate formalism provides an intuitive and quantitative explanation of intriguing phenomena such as cooperativity, resonance, and nonlinearity in the collective vibrational strong coupling regime, as demonstrated in recent simulations.

preprint2022arXiv

Higher-Order Photon Statistics as a New Tool to Reveal Hidden Excited States in a Plasmonic Cavity

Among the best known quantities obtainable from photon correlation measurements are the $g^{(m)}$~correlation functions. Here, we introduce a new procedure to evaluate these correlation functions based on higher-order factorial cumulants $C_{\text{F},m}$ which integrate over the time dependence of the correlation functions, i.e., summarize the available information at different time spans. In a systematic manner, the information content of higher-order correlation functions as well as the distribution of photon waiting times is taken into account. Our procedure greatly enhances the sensitivity for probing correlations and, moreover, is robust against a limited counting efficiency and time resolution in experiment. It can be applied even in case $g^{(m)}$ is not accessible at short time spans. We use the new evaluation scheme to analyze the photon emission of a plasmonic cavity coupled to a single quantum dot. We derive criteria which must hold if the system can be described by a generic Jaynes-Cummings model. A violation of the criteria can be explained by the presence of an additional excited quantum dot state.

preprint2022arXiv

Long-Range Non-Equilibrium Coherent Tunneling Induced by Fractional Vibronic Resonances

We study the influence of a linear energy bias on a non-equilibrium excitation on a chain of molecules coupled to local phonons (a tilted Holstein model) using both a random-walk rate kernel theory and a nonperturbative, massively parallelized adaptive-basis algorithm. We uncover structured and discrete vibronic resonance behavior fundamentally different from both linear response theory and homogeneous polaron dynamics. Remarkably, resonance between the phonon energy $\hbarω$ and the bias $δ_ε$ occurs not only at integer but also fractional ratios $δ_ε/(\hbarω) = \frac{m}{n}$, which effect long-range $n$-bond $m$-phonon tunneling. These observations are also reproduced in a model calculation of a recently demonstrated Cy3 system. Potential applications range from molecular electronics to optical lattices and artificial light harvesting via vibronic engineering of coherent quantum transport.

preprint2022arXiv

Unusual Dynamical Properties of Disordered Polaritons in Micocavities

The strong light-matter interaction in microcavities gives rise to intriguing phenomena, such as cavity-mediated transport that can potentially overcome the Anderson localization. Yet, an accurate theoretical treatment is challenging as the matter (e.g.,molecules) are subject to large energetic disorder. In this article, we develop the Green's function solution to the Fano-Anderson model and use the exact analytical solution to quantify the effects of energetic disorder on the spectral and transport properties in microcavities. Starting from microscopic equation of motions, we derive an effective non-Hermitian Hamiltonian and predict a set of scaling laws: (i) The complex eigen-energies of the effective Hamiltonian exhibit an exceptional point, which leads to underdamped coherent dynamics in the weak disorder regime, where the decay rate increases with disorder, and overdamped incoherent dynamics in the strong disorder regime, where the slow decay rate decreases with disorder. (ii) The total density of states of disordered ensembles can be exactly partitioned into the cavity, bright-state and dark-state local density of states, which are determined by the complex eigen solutions and can be measured via spectroscopy. (iii) The cavity-mediated relaxation and transport dynamics are intimately related such that the energy-resolved relaxation and transport rates are proportional to the cavity local density of states. The ratio of the disorder averaged relaxation and transport rates equals the molecule number, which can be interpreted as a result of a quantum random walk. (iv) A turnover in the rates as a function of disorder or molecule density can be explained in terms of the overlap of the disorder distribution function and the cavity local density of states. These findings reveal the significant impact of the dark states on the transport properties of disordered ensembles in cavities.

preprint2020arXiv

Magnetic field induced symmetry breaking in nonequilibrium quantum networks

We study the effect of an applied magnetic field on the nonequilibrium transport properties of a general cubic quantum network described by a tight-binding Hamiltonian with specially designed couplings to the leads that preserve open-system symmetries. We demonstrate that the symmetry of open systems can be manipulated by the direction of the magnetic field. Starting with all the symmetries preserved in absence of a field, the anisotropic and isotropic fields systematically break the symmetries, influencing all nonequilibrium properties. For simple cubic systems, we are able to identify the steady states that comprise of pure states, bath-dependent states (nonequilibrium steady states), and also nonphysical states. As an application, we show numerically for large cubic networks that the symmetry breaking can control nonequilibrium currents and that different environmental interactions can lead to novel features which can be engineered in artificial super-lattices and cold atoms.

preprint2020arXiv

Steady-State Analysis of Light-harvesting Energy Transfer Driven by Incoherent Light: From Dimers to Networks

The question of how quantum coherence facilitates energy transfer has been intensively debated in the scientific community. Since natural and artificial light-harvesting units operate under the stationary condition, we address this question via a non-equilibrium steady-state analysis of a molecular dimer irradiated by incoherent sunlight and then generalize the key predictions to arbitrarily-complex exciton networks. The central result of the steady-state analysis is the coherence-flux-efficiency relation:$η=c\sum_{i\neq j}F_{ij}κ_j=2c\sum_{i\neq j}J_{ij}{\rm Im}[ρ_{ij}]κ_j$ with $c$ the normalization constant. In this relation, the first equality indicates that energy transfer efficiency $η$ is uniquely determined by the trapping flux, which is the product of flux $F$ and branching ratio $κ$ for trapping at the reaction centers, and the second equality indicates that the energy transfer flux $F$ is equivalent to quantum coherence measured by the imaginary part of the off-diagonal density matrix, i.e., $F_{ij}=2J_{ij}{\rm Im}[ρ_{ij}]$. Consequently, maximal steady-state coherence gives rise to optimal efficiency. The coherence-flux-efficiency relation holds rigorously and generally for any exciton networks of arbitrary connectivity under the stationary condition and is not limited to incoherent radiation or incoherent pumping. For light-harvesting systems under incoherent light, non-equilibrium energy transfer flux (i.e. steady-state coherence) is driven by the breakdown of detailed balance and by the quantum interference of light-excitations and leads to the optimization of energy transfer efficiency. It should be noted that the steady-state coherence or, equivalently, efficiency is the combined result of light-induced transient coherence, inhomogeneous depletion, and system-bath correlation, and is thus not necessarily correlated with quantum beatings.

preprint2019arXiv

Generalized Kasha's Scheme for Classifying Two-Dimensional Excitonic Molecular Aggregates: Temperature Dependent Absorption Peak Frequency Shift

We propose a generalized theoretical framework for classifying two-dimensional (2D) excitonic molecular aggregates based on an analysis of temperature dependent spectra. In addition to the monomer-aggregate absorption peak shift, which defines the conventional J- and H-aggregates, we incorporate the peak shift associated with increasing temperature as a measure to characterize the exciton band structure. First we show that there is a one-to-one correspondence between the monomer-aggregate and the T-dependent peak shifts for Kasha's well-established model of 1D aggregates, where J-aggregates exhibit further redshift upon increasing temperature and H-aggregates exhibit further blueshift. On the contrary, 2D aggregate structures are capable of supporting the two other combinations: blueshifting J-aggregates and redshifting H-aggregates, owing to their more complex exciton band structures. Secondly, using spectral lineshape theory, the T-dependent shift is associated with the relative abundance of states on each side of the bright state. We further establish that the density of states can be connected to the microscopic packing condition leading to these four classes of aggregates by separately considering the short and long-range contribution to the excitonic couplings. In particular the T-dependent shift is shown to be an unambiguous signature for the sign of net short-range couplings: Aggregates with net negative (positive) short-range couplings redshift (blueshift) with increasing temperature. Lastly, comparison with experiments shows that our theory can be utilized to quantitatively account for the observed but previously unexplained T-dependent absorption lineshapes. Thus, our work provides a firm ground for elucidating the structure-function relationships for molecular aggregates and is fully compatible with existing experimental and theoretical structure characterization tools.

preprint2017arXiv

Light Adaptation in Phycobilisome antennas: Influence on the Rod Length and Structural Arrangement

Phycobilisomes, the light-harvesting antennas of cyanobacteria, can adapt to a wide range of environments thanks to a composition and function response to stress conditions. We study how structural changes influence excitation transfer in these super-complexes. Specifically, we show the influence of the rod length on the photon absorption and subsequent excitation transport to the core. Despite the fact that the efficiency of individual disks on the rod decreases with increasing rod length, we find an optimal length for which the average rod efficiency is maximal. Combining this study with experimental structural measurements, we propose models for the arrangement of the phycobiliproteins inside the thylakoid membranes, evaluate the importance of rod length, and predict the corresponding transport properties for different cyanobacterial species. This analysis, which links the functional and structural properties of full phycobilisome complexes, thus provides further rationals to help resolving their exact structure.

preprint2016arXiv

Dynamical signatures of molecular symmetries in nonequilibrium quantum transport

Symmetries play a crucial role in ubiquitous systems found in Nature. In this work, we propose an elegant approach to detect symmetries by measuring quantum currents. Our detection scheme relies on initiating the system in an anti-symmetric initial condition, with respect to the symmetric sites, and using a probe that acts like a local noise. Depending on the position of the probe the currents exhibit unique signatures such as a quasi-stationary plateau indicating the presence of meta-stability and multi-exponential decays in case of multiple symmetries. The signatures are sensitive to the probe and vanish completely when the timescale of the coherent system dynamics is much longer than the timescale of the probe. These results are demonstrated using a $4$-site model and an archetypal example of the para-benzene ring and are shown to be robust under a weak disorder.

preprint2016arXiv

Efficient simulation of non-Markovian system-environment interaction

In this work, we combine an established method for open quantum systems -- the time evolving density matrix using orthogonal polynomials algorithm (TEDOPA) -- with the transfer tensors formalism (TTM), a new tool for the analysis, compression and propagation of non-Markovian processes. This enables the investigation of previously inaccessible long-time dynamics, such as those ensuing from low temperature regimes with arbitrary, possibly highly structured, spectral densities. We briefly introduce both methods, followed by a benchmark to prove viability and combination synergies. Subsequently we illustrate the capabilities of this approach at the hand of specific examples and conclude our analysis by highlighting possible further applications of our method.

preprint2016arXiv

How two-dimensional brick layer J-aggregates differ from linear ones: excitonic properties and line broadening mechanisms

We study the excitonic coupling and homogeneous spectral line width of brick layer J-aggregate films. We begin by analysing the structural information revealed by the two-exciton states probed in two-dimensional spectra. Our first main result is that the relation between the excitonic couplings and the spectral shift in a two-dimensional structure is different (larger shift for the same nearest neighbour coupling) from that in a one-dimensional structure, which leads to an estimation of dipolar coupling in two-dimensional lattices. We next investigate the mechanisms of homogeneous broadening - population relaxation and pure dephasing - and evaluate their relative importance in linear and two-dimensional aggregates. Our second main result is that pure dephasing dominates the line width in two-dimensional systems up to a crossover temperature, which explains the linear temperature dependence of the homogeneous line width. This is directly related to the decreased density of states at the band edge when compared with linear aggregates, thus reducing the contribution of population relaxation to dephasing. Pump-probe experiments are suggested to directly measure the lifetime of the bright state and can therefore support the proposed model.

preprint2016arXiv

Quantum Diffusion on Molecular Tubes: Universal Scaling of the 1D to 2D Transition

The transport properties of disordered systems are known to depend critically on dimensionality. We study the diffusion coefficient of a quantum particle confined to a lattice on the surface of a tube, where it scales between the 1D and 2D limits. It is found that the scaling relation is universal and independent of the disorder and noise parameters, and the essential order parameter is the ratio between the localization length in 2D and the circumference of the tube. Phenomenological and quantitative expressions for transport properties as functions of disorder and noise are obtained and applied to real systems: In the natural chlorosomes found in light-harvesting bacteria the exciton transfer dynamics is predicted to be in the 2D limit, whereas a family of synthetic molecular aggregates is found to be in the homogeneous limit and is independent of dimensionality.

preprint2016arXiv

Quantum transport in d-dimensional lattices

We prove analytically that both fermionic and bosonic uniform d-dimensional lattices can be reduced to a set of independent one-dimensional modes. This reduction leads to the conclusion that the dynamics in uniform fermionic and bosonic lattices is always ballistic. By the use of the Jordan-Wigner transformation we extend our analysis to spin lattices, proving the existence of both ballistic and non-ballistic subspaces in any dimension and for any system size. We then relate the nature of transport with the number of excitations in the spin lattice, indicating that a single excitation propagates always ballistically and that the non-ballistic behavior of uniform spin lattices is a consequence of the interaction between different excitations.

preprint2015arXiv

A Continued Fraction Resummation Form of Bath Relaxation Effect in the Spin-Boson Model

In the spin-boson model, a continued fraction form is proposed to systematically resum high-order quantum kinetic expansion (QKE) rate kernels, accounting for the bath relaxation effect beyond the second-order perturbation. In particular, the analytical expression of the sixth-order QKE rate kernel is derived for resummation. With higher-order correction terms systematically extracted from higher-order rate kernels, the resummed quantum kinetic expansion (RQKE) approach in the continued fraction form extends the Pade approximation and can fully recover the exact quantum dynamics as the expansion order increases.

preprint2015arXiv

Effects of system-bath entanglement on the performance of light-harvesting systems: A quantum heat engine perspective

We explore energy transfer in a generic three-level system, which is coupled to three non-equilibrium baths. Built on the concept of quantum heat engine, our three-level model describes non-equilibrium quantum processes including light-harvesting energy transfer, nano-scale heat transfer, photo-induced isomerization, and photovoltaics in double quantum-dots. In the context of light-harvesting, the excitation energy is first pumped up by sunlight, then is transferred via two excited states which are coupled to a phonon bath, and finally decays to the ground state. The efficiency of this process is evaluated by steady state analysis via a polaron-transformed master equation; thus a wide range of the system-phonon coupling strength can be covered. We show that the coupling with the phonon bath not only modifies the steady state, resulting in population inversion, but also introduces a finite steady state coherence which optimizes the energy transfer flux and efficiency. In the strong coupling limit, the steady state coherence disappears and the efficiency approaches the heat engine limit given by Scovil and Schultz-Dubois in Phys. Rew. Lett. 2, 262 (1959).

preprint2015arXiv

Forster resonance energy transfer, absorption and emission spectra in multichromophoric systems: III. Exact stochastic path integral evaluation

A numerically exact path integral treatment of the absorption and emission spectra of open quantum systems is presented that requires only the straightforward solution of a stochastic differential equation. The approach converges rapidly enabling the calculation of spectra of large excitonic systems across the complete range of system parameters and for arbitrary bath spectral densities. With the numerically exact absorption and emission operators one can also immediately compute energy transfer rates using the multi-chromophoric Forster resonant energy transfer formalism. Benchmark calculations on the emission spectra of two level systems are presented demonstrating the efficacy of the stochastic approach. This is followed by calculations of the energy transfer rates between two weakly coupled dimer systems as a function of temperature and system-bath coupling strength. It is shown that the recently developed hybrid cumulant expansion is the only perturbative method capable of generating uniformly reliable energy transfer rates and spectra across a broad range of system parameters.

preprint2015arXiv

Non-Markovian Dynamical Maps: Numerical Processing of Open Quantum Trajectories

The initial stages of the evolution of an open quantum system encode the key information of its underlying dynamical correlations, which in turn can predict the trajectory at later stages. We propose a general approach based on non-Markovian dynamical maps to extract this information from the initial trajectories and compress it into non-Markovian transfer tensors. Assuming time-translational invariance, the tensors can be used to accurately and efficiently propagate the state of the system to arbitrarily long timescales. The non-Markovian Transfer Tensor Method (TTM) demonstrates the coherent-to-incoherent transition as a function of the strength of quantum dissipation and predicts the non-canonical equilibrium distribution due to the system-bath entanglement. TTM is equivalent to solving the Nakajima-Zwanzig equation, and therefore can be used to reconstruct the dynamical operators (the system Hamiltonian and memory kernel) from quantum trajectories obtained in simulations or experiments. The concept underlying the approach can be generalized to physical observables with the goal of learning and manipulating the trajectories of an open quantum system.

preprint2015arXiv

Stiffening of Red Blood Cells Induced by Disordered Cytoskeleton Structures: A Joint Theory-experiment Study

The functions and elasticities of the cell are largely related to the structures of the cytoskeletons underlying the lipid bi-layer. Among various cell types, the Red Blood Cell (RBC) possesses a relatively simple cytoskeletal structure. Underneath the membrane, the RBC cytoskeleton takes the form of a two dimensional triangular network, consisting of nodes of actins (and other proteins) and edges of spectrins. Recent experiments focusing on the malaria infected RBCs (iRBCs) showed that there is a correlation between the elongation of spectrins in the cytoskeletal network and the stiffening of the iRBCs. Here we rationalize the correlation between these two observations by combining the worm-like chain (WLC) model for single spectrins and the Effective Medium Theory (EMT) for the network elasticity. We specifically focus on how the disorders in the cytoskeletal network affect its macroscopic elasticity. Analytical and numerical solutions from our model reveal that the stiffness of the membrane increases with increasing end-to-end distances of spectrins, but has a non-monotonic dependence on the variance of the end-to-end distance distributions. These predictions are verified quantitively by our AFM and micropipette aspiration measurements of iRBCs. The model may, from a molecular level, provide guidelines for future identification of new treatment methods for RBC related diseases, such as malaria infection.

preprint2014arXiv

Can Natural Sunlight Induce Coherent Exciton Dynamics?

Excitation of a model photosynthetic molecular aggregate by incoherent sunlight is systematically examined. For a closed system, the excited state coherence induced by the sunlight oscillates with an average amplitude that is inversely proportional to the excitonic gap, and reaches a stationary amplitude that depends on the temperature and coherence time of the radiation field. For an open system, the light-induced dynamical coherence relaxes to a static coherence determined by the non-canonical thermal distribution resulting from the entanglement with the phonon bath. The decay of the excited state population to the common ground state establishes a non-equilibrium steady-state flux driven by the sunlight, and it defines a time window to observe the transition from dynamical to static coherence. For the parameters relevant to photosynthetic systems, the exciton dynamics initiated by the sunlight exhibits a non-negligible amount of dynamical coherence (quantum beats) on the sub-picosecond timescale; however, this sub-picosecond time-scale is long enough for light-harvesting systems to establish static coherence, which plays a crucial role in efficient energy transfer. Further, a relationship is established between the non-equilibrium steady-state induced by the sunlight and the coherent dynamics initiated from the ground state by a laser $δ$-pulse, thereby making a direct connection between incoherent sunlight excitation and ultrafast spectroscopy.

preprint2014arXiv

Coherent exciton dynamics in the presence of underdamped vibrations

Recent ultrafast optical experiments show that excitons in large biological light-harvesting complexes are coupled to molecular vibration modes. These high-frequency vibrations will not only affect the optical response, but also drive the exciton transport. Here, using a model dimer system, the frequency of the underdamped vibration is shown to have a strong effect on the exciton dynamics such that quantum coherent oscillations in the system can be present even in the case of strong noise. Two mechanisms are identified to be responsible for the enhanced transport efficiency: critical damping due to the tunable effective strength of the coupling to the bath, and resonance coupling where the vibrational frequency coincides with the energy gap in the system. The interplay of these two mechanisms determines parameters responsible for the most efficient transport, and these optimal control parameters are comparable to those in realistic light-harvesting complexes. Interestingly, oscillations in the excitonic coherence at resonance are suppressed in comparison to the case of an off-resonant vibration.

preprint2014arXiv

Coherent quantum transport in disordered systems: A unified polaron treatment of hopping and band-like transport

Quantum transport in disordered systems is studied using a polaron-based master equation. The polaron approach is capable of bridging the results from the coherent band-like transport regime governed by the Redfield equation to incoherent hopping transport in the classical regime. A non-monotonic dependence of the diffusion coefficient is observed both as a function of temperature and system-phonon coupling strength. In the band-like transport regime, the diffusion coefficient is shown to be linearly proportional to the system-phonon coupling strength, and vanishes at zero coupling due to Anderson localization. In the opposite classical hopping regime, we correctly recover that the dynamics are described by the Fermi's Golden Rule (FGR) and establish that the scaling of the diffusion coefficient depends on the phonon bath relaxation time. In both the hopping and band-like transport regimes, it is demonstrated that at low temperature the zero-point fluctuations of the bath lead to non-zero transport rates, and hence a finite diffusion constant. Application to rubrene and other organic semiconductor materials shows a good agreement with experimental mobility data.

preprint2014arXiv

Förster resonance energy transfer, absorption and emission spectra in multichromophoric systems: I. Cumulant expansions

We study the Förster resonant energy transfer (FRET) rate in multichromophoric systems. The multichromophoric FRET rate is determined by the overlap integral of the donor's emission and acceptor's absorption spectra, which are obtained via 2nd-order cumulant expansion techniques developed in this work. We calculate the spectra and multichromophoric FRET rate for both localized and delocalized systems. (i) The role of the initial entanglement between the donor and its bath is found to be crucial in both the emission spectrum and the multichromophoric FRET rate. (ii) The absorption spectra obtained by the cumulant expansion method are quite close to the exact one for both localized and delocalized systems, even when the system-bath coupling is far from the perturbative regime. (iii) For the emission spectra, the cumulant expansion can give very good results for the localized system, but fail to obtain reliable spectra of the high excitations of a delocalized system, when the system-bath coupling is large and the thermal energy is small. (iv) Even though, the multichromophoric FRET rate is good enough since it is determined by the overlap integral of the spectra.

preprint2014arXiv

Nonequilibrium Energy Transfer at Nanoscale: A Unified Theory from Weak to Strong Coupling

We investigate the microscopic mechanism of quantum energy transfer in the nonequilibrium spin-boson model. By developing a nonequilibrium polaron-transformed Redfield equation based on fluctuation decoupling, we dissect the energy transfer into multi-boson associated processes with even or odd parity. Based on this, we analytically evaluate the energy flux, which smoothly bridges the transfer dynamics from the weak spin-boson coupling regime to the strong-coupling one. Our analysis explains previous limiting predictions and provides a unified interpretation of several observations, including coherence-enhanced heat flux and absence of negative differential thermal conductance in the nonequilibrium spin-boson model. The results may find wide applications for the energy and information control in nanodevices.

preprint2014arXiv

Spectrins in Axonal Cytoskeletons: Dynamics Revealed by Extensions and Fluctuations

The macroscopic properties, the properties of individual components and how those components interact with each other are three important aspects of a composited structure. An understanding of the interplay between them is essential in the study of complex systems. Using axonal cytoskeleton as an example system, here we perform a theoretical study of slender structures that can be coarse-grained as a simple smooth 3-dimensional curve. We first present a generic model for such systems based on the fundamental theorem of curves. We use this generic model to demonstrate the applicability of the well-known worm-like chain (WLC) model to the network level and investigate the situation when the system is stretched by strong forces (weakly bending limit). We specifically studied recent experimental observations that revealed the hitherto unknown periodic cytoskeleton structure of axons and measured the longitudinal fluctuations. Instead of focusing on single molecules, we apply analytical results from the WLC model to both single molecule and network levels and focus on the relations between extensions and fluctuations. We show how this approach introduces constraints to possible local dynamics of the spectrin tetramers in the axonal cytoskeleton and finally suggests simple but self-consistent dynamics of spectrins in which the spectrins in one spatial period of axons fluctuate in-sync.

preprint2013arXiv

A hybrid stochastic hierarchy equations of motion approach to treat the low temperature dynamics of non-Markovian open quantum systems

The hierarchical equations of motion technique has found widespread success as a tool to generate the numerically exact dynamics of non-Markovian open quantum systems. However, its application to low temperature environments remains a serious challenge due to the need for a deep hierarchy that arises from the Matsubara expansion of the bath correlation function. Here we present a hybrid stochastic hierarchical equation of motion (sHEOM) approach that alleviates this bottleneck and leads to a numerical cost that is nearly independent of temperature. Additionally, the sHEOM method generally converges with fewer hierarchy tiers allowing for the treatment of larger systems. Benchmark calculations are presented on the dynamics of two level systems at both high and low temperatures to demonstrate the efficacy of the approach. Then the hybrid method is used to generate the exact dynamics of systems that are nearly impossible to treat by the standard hierarchy. First, exact energy transfer rates are calculated across a broad range of temperatures revealing the deviations from the Forster rates. This is followed by computations of the entanglement dynamics in a system of two qubits at low temperature spanning the weak to strong system-bath coupling regimes.

preprint2013arXiv

A Novel Construction of Complex-valued Gaussian Processes with Arbitrary Spectral Densities and its Application to Excitation Energy Transfer

The recent experimental discoveries about excitation energy transfer (EET) in light harvesting antenna (LHA) attract a lot of interest. As an open non-equilibrium quantum system, the EET demands more rigorous theoretical framework to understand the interaction between system and environment and therein the evolution of reduced density matrix. A phonon is often used to model the fluctuating environment and convolutes the reduced quantum system temporarily. In this paper, we propose a novel way to construct complex-valued Gaussian processes to describe thermal quantum phonon bath exactly by converting the convolution of influence functional into the time correlation of complex Gaussian random field. Based on the construction, we propose a rigorous and efficient computational method, the covariance decomposition (CD) and conditional propagation scheme, to simulate the temporarily entangled reduced system. The new method allows us to study the non-Markovian effect without perturbation under the influence of different spectral densities of the linear system-phonon coupling coefficients. Its application in the study of EET in the Fenna-Matthews-Olson (FMO) model Hamiltonian under four different spectral densities is discussed. Since the scaling of our algorithm is linear due to its Monte Carlo nature, the future application of the method for large LHA systems is attractive. In addition, this method can be used to study the effect of correlated initial condition on the reduced dynamics in the future.

preprint2013arXiv

Coherent quantum transport in disordered systems I: The influence of dephasing on the transport properties and absorption spectra on one-dimensional systems

Excitonic transport in static disordered one dimensional systems is studied in the presence of thermal fluctuations that are described by the Haken-Strobl-Reineker model. For short times, non-diffusive behavior is observed that can be characterized as the free-particle dynamics in the Anderson localized system. Over longer time scales, the environment-induced dephasing is sufficient to overcome the Anderson localization caused by the disorder and allow for transport to occur which is always seen to be diffusive. In the limiting regimes of weak and strong dephasing quantum master equations are developed, and their respective scaling relations imply the existence of a maximum in the diffusion constant as a function of the dephasing rate that is confirmed numerically. In the weak dephasing regime, it is demonstrated that the diffusion constant is proportional to the square of the localization length which leads to a significant enhancement of the transport rate over the classical prediction. Finally, the influence of noise and disorder on the absorption spectrum is presented and its relationship to the transport properties is discussed.

preprint2013arXiv

Correlated local bending of DNA double helix and its effect on the cyclization of short DNA fragments

We report a theoretical study of DNA flexibility and quantitatively predict the ring closure probability as a function of DNA contour length. Recent experimental studies show that the flexibility of short DNA fragments (as compared to the persistence length of DNA l_P~150 base pairs) cannot be described by the traditional worm-like chain (WLC) model, e.g., the observed ring closure probability is much higher than predicted. To explain these observations, DNA flexibility is investigated with explicit considerations of a new length scale l_D~10 base pairs, over which DNA local bend angles are correlated. In this correlated worm-like chain (C-WLC) model, a finite length correction term is analytically derived and the persistence length is found to be contour length dependent. While our model reduces to the traditional worm-like chain model when treating long DNA at length scales much larger than l_P, it predicts that DNA becomes much more flexible at shorter sizes, which helps explain recent cyclization measurements of short DNA fragments around 100 base pairs.

preprint2013arXiv

Higher-Order Kinetic Expansion of Quantum Dissipative Dynamics: Mapping Quantum Networks to Kinetic Networks

We apply a new formalism to derive the higher-order quantum kinetic expansion (QKE) for studying dissipative dynamics in a general quantum network coupled with an arbitrary thermal bath. The dynamics of system population is described by a time-convoluted kinetic equation, where the time-nonlocal rate kernel is systematically expanded on the order of off-diagonal elements of the system Hamiltonian. In the second order, the rate kernel recovers the expression of the noninteracting-blip approximation (NIBA) method. The higher-order corrections in the rate kernel account for the effects of the multi-site quantum coherence and the bath relaxation. In a quantum harmonic bath, the rate kernels of different orders are analytically derived. As demonstrated by four examples, the higher-order QKE can reliably predict quantum dissipative dynamics, comparing well with the hierarchic equation approach. More importantly, the higher-order rate kernels can distinguish and quantify distinct nontrivial quantum coherent effects, such as long-range energy transfer from quantum tunneling and quantum interference arising from the phase accumulation of interactions.

preprint2013arXiv

Optimal tunneling enhances quantum photovoltaic effect in double quantum dots

We investigate quantum photovoltaic effect in double quantum dots by applying nonequilibrium quantum master equation. The drastic suppression of the photovoltaic current is observed near the open circuit voltage, which leads to the large filling factor. We find that there always exists an optimal inter-dot tunneling that significantly enhances the photovoltaic current. The maximal output power will also be obtained around the optimal inter-dot tunneling. Moreover, the open circuit voltage approximately behaves as the product of the eigen-level gap and the Carnot efficiency. These results suggest a great potential for double quantum dots as efficient photovoltaic devices.

preprint2013arXiv

Shape transition of unstrained flattest single-walled carbon nanotubes under pressure

Single walled carbon nanotube's cross section can be flattened under hydrostatic pressure. One example is the cross section of a single walled carbon nanotube successively deforms from the original round shape to oval shape, then to peanut-like shape. At the transition point of reversible deformation between convex shape and concave shape, the side wall of nanotube is flattest. This flattest tube has many attractive properties. In the present work, an approximate approach is developed to determine the equilibrium shape of this unstrained flattest tube and the curvature distribution of this tube. Our results are in good agreement with recent numerical results, and can be applied to the study of pressure controlled electric properties of single walled carbon nanotubes. The present method can also be used to study other deformed inorganic and organic tube-like structures.

preprint2012arXiv

Accuracy of Second Order Perturbation theory in the Polaron and Variational Polaron Frames

In the study of open quantum systems, the polaron transformation has recently attracted a renewed interest as it offers the possibility to explore the strong system-bath coupling regime. Despite this interest, a clear and unambiguous analysis of the regimes of validity of the polaron transformation is still lacking. Here we provide such a benchmark, comparing second order perturbation theory results in the original untransformed frame, the polaron frame and the variational extension with numerically exact path integral calculations of the equilibrium reduced density matrix. Equilibrium quantities allow a direct comparison of the three methods without invoking any further approximations as is usually required in deriving master equations. It is found that the second order results in the original frame are accurate for weak system-bath coupling, the full polaron results are accurate in the opposite regime of strong coupling, and the variational method is capable of interpolating between these two extremes. As the bath becomes more non-Markovian (slow bath), all three approaches become less accurate.

preprint2012arXiv

An exact equilibrium reduced density matrix formulation I: The influence of noise, disorder, and temperature on localization in excitonic systems

An exact method to compute the entire equilibrium reduced density matrix for systems characterized by a system-bath Hamiltonian is presented. The approach is based upon a stochastic unraveling of the influence functional that appears in the imaginary time path integral formalism of quantum statistical mechanics. This method is then applied to study the effects of thermal noise, static disorder, and temperature on the coherence length in excitonic systems. As representative examples of biased and unbiased systems, attention is focused on the well-characterized light harvesting complexes of FMO and LH2, respectively. Due to the bias, FMO is completely localized in the site basis at low temperatures, whereas LH2 is completely delocalized. In the latter, the presence of static disorder leads to a plateau in the coherence length at low temperature that becomes increasingly pronounced with increasing strength of the disorder. The introduction of noise, however, precludes this effect. In biased systems, it is shown that the environment may increase the coherence length, but only decrease that of unbiased systems. Finally it is emphasized that for typical values of the environmental parameters in light harvesting systems, the system and bath are entangled at equilibrium in the single excitation manifold. That is, the density matrix cannot be described as a product state as is often assumed, even at room temperature. The reduced density matrix of LH2 is shown to be in precise agreement with the steady state limit of previous exact quantum dynamics calculations.

preprint2012arXiv

Efficient Energy Transfer in Light-Harvesting Systems, II: Quantum-Classical Comparison, Flux Network, and Robustness Analysis

Following the calculation of optimal energy transfer in thermal environment in our first paper (Wu et al., New J. Phys., 2010, 12, 105012), full quantum dynamics and leading-order `classical' hopping kinetics are compared in the seven-site Fenna-Matthews-Olson (FMO) protein complex. The difference between these two dynamic descriptions is due to higher-order quantum corrections. Two thermal bath models, classical white noise (the Haken-Strobl-Reineker model) and quantum Debye model, are considered. In the seven-site FMO model, we observe that higher-order corrections lead to negligible changes in the trapping time or in energy transfer efficiency around the optimal and physiological conditions (2% in the HSR model and 0.1% in the quantum Debye model for the initial site at BChl 1). However, using the concept of integrated flux, we can identify significant differences in branching probabilities of the energy transfer network between hopping kinetics and quantum dynamics (26% in the HSR model and 32% in the quantum Debye model for the initial site at BChl 1). This observation indicates that the quantum coherence can significantly change the distribution of energy transfer pathways in the flux network with the efficiency nearly the same. The quantum-classical comparison of the average trapping time with the removal of the bottleneck site, BChl 4, demonstrates the robustness of the efficient energy transfer by the mechanism of multi-site quantum coherence. To reconcile with the latest eight-site FMO model, the quantum-classical comparison with the flux network analysis is summarized in the appendix. The eight-site FMO model yields similar trapping time and network structure as the seven-site FMO model but leads to a more disperse distribution of energy transfer pathways.

preprint2012arXiv

Excitonic energy transfer in light-harvesting complexes in purple bacteria

Two distinct approaches, the Frenkel-Dirac time-dependent variation and the Haken-Strobl model, are adopted to study energy transfer dynamics in single-ring and double-ring light-harvesting systems in purple bacteria. It is found that inclusion of long-range dipolar interactions in the two methods results in significant increases in intra- or inter-ring exciton transfer efficiency. The dependence of exciton transfer efficiency on trapping positions on single rings of LH2 (B850) and LH1 is similar to that in toy models with nearest-neighbor coupling only. However, owing to the symmetry breaking caused by the dimerization of BChls and dipolar couplings, such dependence has been largely suppressed. In the studies of coupled-ring systems, both methods reveal interesting role of dipolar interaction in increasing energy transfer efficiency by introducing multiple intra/inter-ring transfer paths. Importantly, the time scale (~4ps) of inter-ring exciton transfer obtained from polaron dynamics is in good agreement with previous studies. In a double-ring LH2 system, dipole-induced symmetry breaking leads to global minima and local minima of the average trapping time when there is a finite value of non-zero dephasing rate, suggesting that environment plays a role in preserving quantum coherent energy transfer. In contrast, dephasing comes into play only when the perfect cylindrical symmetry in the hypothetic system is broken. This study has revealed that dipolar interaction between chromophores may play an important part in the high energy transfer efficiency in the LH2 system and many other natural photosynthetic systems.

preprint2012arXiv

Generic Mechanism of Optimal Energy Transfer Efficiency: A Scaling Theory of the Mean First Passage Time in Exciton Systems

An asymptotic scaling theory is presented using the conceptual basis of trapping-free subspace (i.e., orthogonal subspace) to establish the generic mechanism of optimal efficiency of excitation energy transfer (EET) in light-harvesting systems. Analogous to Kramers' turnover in classical rate theory, the enhanced efficiency in the weak damping limit and the suppressed efficiency in the strong damping limit define two asymptotic scaling regimes, which are interpolated to predict the functional form of optimal efficiency of the trapping-free subspace. In the presence of static disorder, the scaling law of transfer time with respect to dephasing rate changes from linear to square root, suggesting a weaker dependence on the environment. Though formulated in the context of EET, the analysis and conclusions apply in general to open quantum processes, including electron transfer, fluorescence emission, and heat conduction.

preprint2012arXiv

Non-Canonical Statistics of a Spin-Boson Model: Theory and Exact Monte-Carlo Simulations

Equilibrium canonical distribution in statistical mechanics assumes weak system-bath coupling (SBC). In real physical situations this assumption can be invalid and equilibrium quantum statistics of the system may be non-canonical. By exploiting both polaron transformation and perturbation theory in a spin-boson model, an analytical treatment is advocated to study non-canonical statistics of a two-level system at arbitrary temperature and for arbitrary SBC strength, yielding theoretical results in agreement with exact Monte-Carlo simulations. In particular, the eigen-representation of system's reduced density matrix is used to quantify non-canonical statistics as well as the quantumness of the open system. For example, it is found that irrespective of SBC strength, non-canonical statistics enhances as temperature decreases but vanishes at high temperature.

preprint2011arXiv

Efficient energy transfer in light-harvesting systems, III: The influence of the eighth bacteriochlorophyll on the dynamics and efficiency in FMO

The most recent crystal structure of the Fenna-Matthews-Olson (FMO) protein complex indicates that each subunit contains an additional eighth chromophore. It has been proposed that this extra site functions as a link between the chlorosome antenna complex and the remaining seven chromophores in FMO [Schmidt am Busch et al, J. Phys. Chem. Lett., {\bf 2}, 93 (2011)]. Here, we investigate the implications of this scenario through numerical calculations with the generalized Bloch-Redfield (GBR) equation and the non-interacting blip approximation (NIBA). Three key insights into the population dynamics and energy transfer efficiency in FMO are provided. First, it is shown that the oscillations that are often observed in the population relaxation of the dimer composed of sites one and two may be completely suppressed in the eight site model. The presence of the coherent oscillations is shown to depend upon the particular initial preparation of the dimer state. Secondly it is demonstrated that while the presence of the eighth chromophore does not cause a dramatic change in the energy transfer efficiency, it does however lead to a dominant energy transfer pathway which can be characterized by an effective three site system arranged in an equally spaced downhill configuration. Such a configuration leads to an optimal value of the site energy of the eighth chromophore which is shown to be near to its suggested value. Finally we confirm that the energy transfer process in the eight site FMO complex remains efficient and robust. The optimal values of the bath parameters are computed and shown to be closer to the experimentally fitted values than those calculated previously for the seven site system.

preprint2010arXiv

Efficient energy transfer in light-harvesting systems, I: optimal temperature, reorganization energy, and spatial-temporal correlations

Understanding the mechanisms of efficient and robust energy transfer in light-harvesting systems provides new insights for the optimal design of artificial systems. In this paper, we use the Fenna-Matthews-Olson (FMO) protein complex and phycocyanin 645 (PC 645) to explore the general dependence on physical parameters that help maximize the efficiency and maintain its stability. With the Haken-Strobl model, the maximal energy transfer efficiency (ETE) is achieved under an intermediate optimal value of dephasing rate. To avoid the infinite temperature assumption in the Haken-Strobl model and the failure of the Redfield equation in predicting the Forster rate behavior, we use the generalized Bloch-Redfield (GBR) equation approach to correctly describe dissipative exciton dynamics and find that maximal ETE can be achieved under various physical conditions, including temperature, reorganization energy, and spatial-temporal correlations in noise. We also identify regimes of reorganization energy where the ETE changes monotonically with temperature or spatial correlation and therefore cannot be optimized with respect to these two variables.

preprint2010arXiv

Optimal Efficiency of Self-Assembling Light-Harvesting Arrays

Using a classical master equation that describes energy transfer over a given lattice, we explore how energy transfer efficiency along with the photon capturing ability depends on network connectivity, on transfer rates, and on volume fractions - the numbers and relative ratio of fluorescence chromophore components, e.g., donor (D), acceptor (A), and bridge (B) chromophores. For a one-dimensional AD array, the exact analytical expression for efficiency shows a steep increase with a D-to-A transfer rate when a spontaneous decay is sufficiently slow. This result implies that the introduction of B chromophores can be a useful method for improving efficiency for a two-component AD system with inefficient D-to-A transfer and slow spontaneous decay. Analysis of this one-dimensional system can be extended to higher-dimensional systems with chromophores arranged in structures such as a helical or stacked-disk rod, which models the self-assembling monomers of the tobacco mosaic virus coat protein. For the stacked-disk rod, we observe the following: (1) With spacings between sites fixed, a staggered conformation is more efficient than an eclipsed conformation. (2) For a given ratio of A and D chromophores, the uniform distribution of acceptors that minimizes the mean first passage time to acceptors is a key point to designing the optimal network for a donor-acceptor system with a relatively small D-to-A transfer rate. (3) For a three-component ABD system with a large B-to-A transfer rate, a key design strategy is to increase the number of the pathways in accordance with the directional energy flow from D to B to A chromophores.