Researcher profile

Zhigang Shuai

Zhigang Shuai contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

6 published item(s)

preprint2022arXiv

Hybrid Quantum-Classical Boson Sampling Algorithm for Molecular Vibrationally Resolved Electronic Spectroscopy with Duschinsky Rotation and Anharmonicity

Using a photonic quantum computer for boson sampling has been demonstrated a tremendous advantage over classical supercomputers. It is highly desirable to develop boson sampling algorithms for realistic scientific problems. In this work, we propose a hybrid quantum-classical sampling (HQCS) algorithm to calculate the optical spectrum for complex molecules considering anharmonicity and Duschinsky rotation (DR) effects. The classical sum-over-state method for this problem has a computational complexity that exponentially increases with system size. In the HQCS algorithm, an intermediate harmonic potential energy surface (PES) is created, bridging the initial and final PESs. The magnitude and sign (-1 or +1) of the overlap between the initial state and the intermediate state are estimated by quantum boson sampling and by classical algorithms respectively, achieving an exponential speed-up. Additionally, the overlap between the intermediate state and the final state is efficiently evaluated by classical algorithms. The feasibility of HQCS is demonstrated in calculations of the emission spectrum of a Morse model as well as pyridine molecule by comparison with the nearly exact time-dependent density matrix renormalization group solutions. A near-term quantum advantage for realistic molecular spectroscopy simulation is proposed.

preprint2022arXiv

Unified definition of exciton coherence length for exciton-phonon coupled molecular aggregates

Exciton coherence length (ECL) is an essential concept to characterize the nature of exciton in molecular aggregates for photosynthesis, organic photovoltaics, and light-emitting diodes. ECL has been defined in a number of ways through the variance or purity of the electronic reduced density matrix. However, we find that these definitions fail to present a monotonic relationship with respect to the exciton radiative decay efficiency as it should be when exciton-phonon couplings are taken into accounts. We propose a unified definition of ECL by virtue of sum rule of oscillator strengths. Using the numerically accurate time-dependent matrix product states formalism applied to Frenkel-Holstein models for molecular aggregates in both one- and two-dimensional system, we find our ECL definition exhibits a monotonic relationship with respect to the radiative efficiency and can serve as an efficient and unified description for exciton coherence. We further predict that the two-dimensional aggregates can display maximum superradiance enhancement (SRE) at finite temperature, different from the previous knowledge of SRE-$1/T$ behavior.

preprint2021arXiv

Non-Markovian Stochastic Schrödinger Equation: Matrix Product State Approach to the Hierarchy of Pure States

We derive a stochastic hierarchy of matrix product states (HOMPS) for non-Markovian dynamics in open quantum system at finite temperature, which is numerically exact and efficient. HOMPS is obtained from the recently developed stochastic hierarchy of pure states (HOPS) by expressing HOPS in terms of formal creation and annihilation operators. The resulting stochastic first order differential equation is then formulated in terms of matrix product states and matrix product operators. In this way the exponential complexity of HOPS can be reduced to scale polynomial with the number of particles. The validity and efficiency of HOMPS is demonstrated for the spin-boson model and long chains where each site is coupled to a structured, strongly non-Markovian environment.

preprint2020arXiv

A General Automatic Method for Optimal Construction of Matrix Product Operators Using Bipartite Graph Theory

Constructing matrix product operators (MPO) is at the core of the modern density matrix renormalization group (DMRG) and its time dependent formulation. For DMRG to be conveniently used in different problems described by different Hamiltonians, in this work we propose a new generic algorithm to construct the MPO of an arbitrary operator with a sum-of-products form based on the bipartite graph theory. We show that the method has the following advantages: (i) It is automatic in that only the definition of the operator is required; (ii) It is symbolic thus free of any numerical error; (iii) The complementary operator technique can be fully employed so that the resulting MPO is globally optimal for any given order of degrees of freedom; (iv) The symmetry of the system could be fully employed to reduce the dimension of MPO. To demonstrate the effectiveness of the new algorithm, the MPOs of Hamiltonians ranging from the prototypical spin-boson model and Holstein model to the more complicated ab initio electronic Hamiltonian and the anharmonic vibrational Hamiltonian with sextic force field are constructed. It is found that for the former three cases, our automatic algorithm can reproduce exactly the same MPOs as the optimally hand-crafted ones already known in the literature.

preprint2020arXiv

Finite Temperature TD-DMRG for the Carrier Mobility of Organic Semiconductors

A large number of non-adiabatic dynamical studies have been applied to reveal the nature of carrier transport in organic semiconductors with different approximations. We present here a &#34;nearly exact&#34; graphical process unit (GPU) based finite temperature TD-DMRG method to evaluate the carrier mobility in organic semiconductors as described by electron-phonon model, in particular to rubrene crystal, one of the prototypical organic semiconductors, with parameters derived from first-principles. We find that (i) TD-DMRG is a general and robust method that can bridge the gap between hopping and band picture covering a wide range of electronic coupling strength; and (ii) with realistic parameters, TD-DMRG is able to account for the experimentally observed &#34;band-like&#34; transport behavior ($\partial μ/ \partial T < 0$) in rubrene. We further study the long-standing puzzle of isotope effect for charge transport and unambiguously demonstrate that the negative isotope effect ($\partial μ/ \partial m < 0$ where $m$ is the atomic mass) should be universal.

preprint2019arXiv

Numerical Assessment for Accuracy and GPU Acceleration of TD-DMRG Time Evolution Schemes

Time dependent density matrix renormalization group (TD-DMRG) has become one of the cutting edge methods of quantum dynamics for complex systems. In this paper, we comparatively study the accuracy of three time evolution schemes in TD-DMRG, the global propagation and compression method with Runge-Kutta algorithm (P&C-RK), the time dependent variational principle based methods with matrix unfolding algorithm (TDVPMU) and with projector-splitting algorithm (TDVP-PS), by performing benchmarks on the exciton dynamics of Fenna-Matthews-Olson (FMO) complex. We show that TDVP-MU and TDVP-PS yield the same result when the time step size is converged and they are more accurate than P&C-RK4, while TDVP-PS tolerates a larger time step size than TDVP-MU. We further adopt the graphical processing units (GPU) to accelerate the heavy tensor contractions in TD-DMRG and it is able to speed up the TDVP-MU and TDVP-PS schemes by up to 73 times.