Researcher profile

Jiajun Ren

Jiajun Ren contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
8works
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

8 published item(s)

preprint2026arXiv

Autonomous Quantum Simulation through Large Language Model Agents

We demonstrate that large language model (LLM) agents can autonomously perform tensor network simulations of quantum many-body systems, achieving approximately 90% success rate across representative benchmark tasks. Tensor network methods are powerful tools for quantum simulation, but their effective use requires expertise typically acquired through years of graduate training. By combining in-context learning with curated documentation and multi-agent decomposition, we create autonomous AI agents that can be trained in specialized computational domains within minutes. We benchmark three configurations (baseline, single-agent with in-context learning, and multi-agent with in-context learning) on problems spanning quantum phase transitions, open quantum system dynamics, and photochemical reactions. Systematic evaluation using DeepSeek-V3.2, Gemini 2.5 Pro, and Claude Opus 4.5 demonstrates that both in-context learning and multi-agent architecture are essential. Analysis of failure modes reveals characteristic patterns across models, with the multi-agent configuration substantially reducing implementation errors and hallucinations compared to simpler architectures.

preprint2025arXiv

Efficient Quantum Simulation of Non-Adiabatic Molecular Dynamics with Precise Electronic Structure

In the study of non-adiabatic chemical processes such as photocatalysis and photosynthesis, non-adiabatic molecular dynamics (NAMD) is an indispensable theoretical tool, which requires precise potential energy surfaces (PESs) of ground and excited states. Quantum computing offers promising potential for calculating PESs that are intractable for classical computers. However, its realistic application poses significant challenges to the development of quantum algorithms that are sufficiently general to enable efficient and precise PES calculations across chemical systems with diverse properties, as well as to seamlessly adapt existing NAMD theories to quantum computing. In this work, we introduce a quantum-adapted extension to the Landau-Zener-Surface-Hopping (LZSH) NAMD. This extension incorporates curvature-driven hopping corrections that protect the population evolution while maintaining the efficiency gained from avoiding the computation of non-adiabatic couplings (NACs), as well as preserving the trajectory independence that enables parallelization. Furthermore, to ensure the high-precision PESs required for surface hopping dynamics, we develop a sub-microhartree-accurate PES calculation protocol. This protocol supports active space selection, enables parallel acceleration either on quantum or classical clusters, and demonstrates adaptability to diverse chemical systems - including the charged H3+ ion and the C2H4 molecule, a prototypical multi-reference benchmark. This work paves the way for practical application of quantum computing in NAMD, showcasing the potential of parallel simulation on quantum-classical heterogeneous clusters for ab-initio computational chemistry.

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.