N-Mode Quantized Anharmonic Vibronic Hamiltonians for Matrix Product State Dynamics
Pith reviewed 2026-05-18 01:44 UTC · model grok-4.3
The pith
N-mode quantization of anharmonic vibronic Hamiltonians supports accurate matrix product state dynamics for photochemical processes.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The n-mode quantization of all vibronic Hamiltonian terms comprised of general high-dimensional model representations results in a second-quantized framework for accurate vibronic calculations employing the density matrix renormalization group algorithm, as demonstrated by calculating the excited state quantum dynamics of maleimide with analysis of convergence and parameters.
What carries the argument
The n-mode quantization procedure applied to anharmonic vibronic Hamiltonian terms and off-diagonal nonadiabatic coupling terms, which converts high-dimensional model representations into a second-quantized form compatible with matrix product state representations for dynamics.
If this is right
- Excited state quantum dynamics calculations become feasible for molecules exhibiting complex anharmonic features using tensor network methods such as the density matrix renormalization group algorithm.
- The time-dependent density matrix renormalization group algorithm can be applied directly to these n-mode quantized Hamiltonians, with convergence properties that can be systematically analyzed.
- Photochemical processes can be modeled while retaining both anharmonic effects in the potential energy surfaces and the full nonadiabatic coupling terms.
- Parameter choices for the underlying time-dependent density matrix renormalization group algorithm can be optimized to achieve reliable results for such quantized vibronic systems.
Where Pith is reading between the lines
- This quantization approach might extend naturally to other photochemical systems where anharmonicity influences reaction pathways or spectra.
- The framework could interface with additional tensor-network or machine-learning representations of potential energy surfaces to handle even larger molecules.
- Testing the method on benchmark systems with known exact solutions would clarify the trade-off between n-mode truncation order and long-time dynamical accuracy.
Load-bearing premise
The n-mode quantization of general high-dimensional model representations is assumed to preserve sufficient accuracy for anharmonic effects and nonadiabatic couplings without introducing uncontrolled approximations that would invalidate the dynamics results.
What would settle it
A direct comparison of the computed excited-state dynamics and spectra for maleimide against experimental time-resolved spectroscopy or higher-accuracy reference calculations on the same system would test whether the quantization preserves the necessary accuracy.
Figures
read the original abstract
Theoretical predictions of photochemical processes are essential for interpreting and understanding spectral features. Reliable quantum dynamics calculations of vibronic systems require precise modeling of anharmonic effects in the potential energy surfaces and off-diagonal nonadiabatic coupling terms. In this work, we present the n-mode quantization of all vibronic Hamiltonian terms comprised of general high-dimensional model representations. This results in a second-quantized framework for accurate vibronic calculations employing the density matrix renormalization group algorithm. We demonstrate the accuracy and reliability of this approach by calculating the excited state quantum dynamics of maleimide. We analyze convergence and the choice of parameters of the underlying time-dependent density matrix renormalization group algorithm for the n-mode vibronic Hamiltonian, demonstrating that it enables accurate calculations of complex photochemical dynamics.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript introduces an n-mode quantization procedure to convert general high-dimensional anharmonic vibronic Hamiltonians into a second-quantized form amenable to matrix product state (MPS) representations. This enables the use of the density matrix renormalization group (DMRG) algorithm for time-dependent quantum dynamics simulations. The approach is demonstrated through calculations of the excited-state dynamics of maleimide, with accompanying analysis of convergence with respect to DMRG parameters such as bond dimension and time step.
Significance. If the n-mode quantization accurately preserves the anharmonic potential features and nonadiabatic couplings, this method could significantly advance the simulation of photochemical processes in complex molecular systems by providing a scalable, accurate framework for vibronic dynamics using tensor network techniques. The convergence analysis in the maleimide example is a positive step toward establishing reliability.
major comments (2)
- [§III] §III (N-Mode Quantization): the procedure for quantizing high-dimensional anharmonic terms and off-diagonal couplings lacks any explicit error metric or direct numerical comparison of the resulting second-quantized Hamiltonian against the original model representation for the maleimide surfaces; this is load-bearing for the central accuracy claim in the dynamics.
- [§V] §V (Results for Maleimide): convergence is reported only versus DMRG bond dimension and time step, with no propagation analysis of quantization truncation errors into the time-dependent observables or benchmark against dynamics on the unquantized Hamiltonian.
minor comments (2)
- [§II] Notation for the original high-dimensional representation versus the n-mode quantized form is occasionally ambiguous in the Hamiltonian definitions; a consistent symbol or subscript would improve clarity.
- [§V] Figure captions for the convergence plots could explicitly state the observable being plotted (e.g., population or coherence) to aid immediate interpretation.
Simulated Author's Rebuttal
We thank the referee for their careful reading of the manuscript and for the constructive comments, which have helped us strengthen the validation of the n-mode quantization procedure. We address each major comment below and describe the revisions planned for the next version of the manuscript.
read point-by-point responses
-
Referee: [§III] §III (N-Mode Quantization): the procedure for quantizing high-dimensional anharmonic terms and off-diagonal couplings lacks any explicit error metric or direct numerical comparison of the resulting second-quantized Hamiltonian against the original model representation for the maleimide surfaces; this is load-bearing for the central accuracy claim in the dynamics.
Authors: We agree that an explicit error metric and direct numerical comparison would provide stronger support for the accuracy of the quantization. The n-mode procedure is constructed to represent the anharmonic potentials and couplings exactly within the chosen mode expansion, but we acknowledge that a quantitative assessment against the original surfaces for maleimide is valuable. In the revised manuscript we will add a dedicated error analysis subsection in §III that reports root-mean-square deviations of the quantized potential and coupling terms relative to the original model at sampled geometries along the maleimide surfaces. This comparison will be used to confirm that quantization errors remain below the threshold relevant for the reported dynamics. revision: yes
-
Referee: [§V] §V (Results for Maleimide): convergence is reported only versus DMRG bond dimension and time step, with no propagation analysis of quantization truncation errors into the time-dependent observables or benchmark against dynamics on the unquantized Hamiltonian.
Authors: We concur that examining the propagation of quantization truncation errors into observables would improve the manuscript. A direct benchmark against the unquantized Hamiltonian is not computationally tractable for the full-dimensional maleimide system, which is the motivation for introducing the second-quantized representation that enables efficient MPS propagation. To address the concern, the revised §V will include additional calculations that vary the n-mode truncation level and demonstrate convergence of the time-dependent excited-state populations and coherences with respect to this parameter. These results will show that the key dynamical features stabilize once a sufficient number of modes is retained, thereby indicating limited propagation of quantization errors into the reported observables. revision: partial
Circularity Check
No circularity: derivation builds on general high-dimensional representations without self-referential reduction
full rationale
The paper introduces an n-mode quantization procedure applied to general high-dimensional vibronic model representations to obtain a second-quantized Hamiltonian suitable for DMRG/MPS dynamics. This construction is presented as a direct mapping from the input model surfaces (including anharmonic and coupling terms) rather than a fit or redefinition that presupposes the output. The maleimide demonstration relies on separate convergence checks with respect to DMRG bond dimension and time-step parameters, which are independent of the quantization step itself. No equations or claims reduce a prediction to a quantity defined by the authors' own prior equations or self-citations in a load-bearing way. The framework remains self-contained against external benchmarks for the quantization procedure.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption General high-dimensional model representations can be quantized in n-mode form while retaining anharmonic and coupling information.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
n-mode expansion of an arbitrary function F(Q) ... F(Q) = Σ F[i]1(Qi) + Σ F[ij]2(Qi,Qj) + ... ; second-quantized form Hnmode = Σ H[i]ki,hi b†ki bhi + ...
-
IndisputableMonolith/Foundation/AbsoluteFloorClosure.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
TD-DMRG propagation with bond-dimension truncation and local Hilbert-space size Nmax
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
Blanchet, V., Zgierski, M. Z., Seideman, T. & Stolow, A. Discerning vibronic molecular dynamics using time-resolved photoelectron spectroscopy.Nature401, 52–54 (1999)
work page 1999
-
[2]
Doppagne, B.et al.Vibronic spectroscopy with submolecular resolution from stm-induced electroluminescence.Phys. Rev. Lett.118, 127401 (2017)
work page 2017
-
[3]
Matselyukh, D. T., Despr´ e, V., Golubev, N. V., Kuleff, A. I. & W¨ orner, H. J. Decoherence and revival in attosecond charge migration driven by non-adiabatic dynamics.Nat. Phys.18, 1206–1213 (2022)
work page 2022
-
[4]
Bloino, J., Biczysko, M., Crescenzi, O. & Barone, V. Integrated computational approach to vibrationally resolved electronic spectra: Anisole as a test case.J. Chem. Phys.128(2008)
work page 2008
-
[5]
Baiardi, A., Bloino, J. & Barone, V. General time dependent approach to vi- bronic spectroscopy including franck–condon, herzberg–teller, and duschinsky effects.J. Chem. Theory Comput.9, 4097–4115 (2013)
work page 2013
-
[6]
Eng, J., Gourlaouen, C., Gindensperger, E. & Daniel, C. Spin-vibronic quantum dynamics for ultrafast excited-state processes.Acc. Chem. Res.48, 809–817 (2015)
work page 2015
-
[7]
J., Chen, T., Candian, A., Lee, T
Mackie, C. J., Chen, T., Candian, A., Lee, T. J. & Tielens, A. G. Fully anhar- monic infrared cascade spectra of polycyclic aromatic hydrocarbons.J. Chem. Phys.149(2018)
work page 2018
-
[8]
McCoy, A. B. & Duncan, M. A. Evidence of anharmonicity in the vibrational spectrum of protonated ethylene.J. Mol. Spectrosc.389, 111704 (2022)
work page 2022
-
[9]
& Bonaˇ ci´ c-Kouteck` y, V.Electronic Aspects of Organic Photochemistry (Wiley & Sons, 1990)
Michl, J. & Bonaˇ ci´ c-Kouteck` y, V.Electronic Aspects of Organic Photochemistry (Wiley & Sons, 1990)
work page 1990
-
[10]
Schoenlein, R. W., Peteanu, L. A., Mathies, R. A. & Shank, C. V. The first step in vision: femtosecond isomerization of rhodopsin.Science254, 412–415 (1991)
work page 1991
-
[11]
A second quantization formulation of multimode dynamics.J
Christiansen, O. A second quantization formulation of multimode dynamics.J. Chem. Phys.120, 2140–2148 (2004). 15
work page 2004
-
[12]
Verstraete, F.et al.Density matrix renormalization group, 30 years on.Nat. Rev. Phys.5, 273–276 (2023)
work page 2023
-
[13]
White, S. R. Density matrix formulation for quantum renormalization groups. Phys. Rev. Lett.69, 2863 (1992)
work page 1992
-
[14]
The density-matrix renormalization group in the age of matrix product states.Ann
Schollw¨ ock, U. The density-matrix renormalization group in the age of matrix product states.Ann. Phys.326, 96–192 (2011)
work page 2011
-
[15]
Baiardi, A. & Reiher, M. The density matrix renormalization group in chemistry and molecular physics: Recent developments and new challenges.J. Chem. Phys. 152(2020)
work page 2020
-
[16]
Ma, H., Schollw¨ ock, U. & Shuai, Z.Density Matrix Renormalization Group (DMRG)-Based Approaches in Computational Chemistry(Elsevier, 2022)
work page 2022
-
[17]
J., Kollath, C., Schollw¨ ock, U
Daley, A. J., Kollath, C., Schollw¨ ock, U. & Vidal, G. Time-dependent density- matrix renormalization-group using adaptive effective hilbertspaces.J. Stat. Mech: Theory Exp.2004, P04005 (2004)
work page 2004
- [18]
-
[19]
Wolf, F. A., Go, A., McCulloch, I. P., Millis, A. J. & Schollw¨ ock, U. Imaginary- time matrix product state impurity solver for dynamical mean-field theory.Phys. Rev. X5, 041032 (2015)
work page 2015
-
[20]
Ren, J., Shuai, Z. & Kin-Lic Chan, G. Time-dependent density matrix renor- malization group algorithms for nearly exact absorption and fluorescence spectra of molecular aggregates at both zero and finite temperature.J. Chem. Theory Comput.14, 5027–5039 (2018)
work page 2018
-
[21]
Baiardi, A. & Reiher, M. Large-scale quantum dynamics with matrix product states.J. Chem. Theory Comput.15, 3481–3498 (2019)
work page 2019
-
[22]
Electron dynamics with the time-dependent density matrix renor- malization group.J
Baiardi, A. Electron dynamics with the time-dependent density matrix renor- malization group.J. Chem. Theory Comput.17, 3320–3334 (2021)
work page 2021
-
[23]
Ren, J., Li, W., Jiang, T., Wang, Y. & Shuai, Z. Time-dependent density matrix renormalization group method for quantum dynamics in complex systems.Wiley Interdiscip. Rev.: Comput. Mol. Sci.12, e1614 (2022)
work page 2022
-
[24]
Efficient simulation of one-dimensional quantum many-body systems
Vidal, G. Efficient simulation of one-dimensional quantum many-body systems. Phys. Rev. Lett.93, 040502 (2004)
work page 2004
-
[25]
White, S. R. & Feiguin, A. E. Real-time evolution using the density matrix renormalization group.Phys. Rev. Lett.93, 076401 (2004). 16
work page 2004
-
[26]
Feiguin, A. E. & White, S. R. Time-step targeting methods for real-time dynam- ics using the density matrix renormalization group.Phys. Rev. B72, 020404 (2005)
work page 2005
-
[27]
Al-Hassanieh, K. A., Feiguin, A. E., Riera, J. A., B¨ usser, C. & Dagotto, E. Adap- tive time-dependent density-matrix renormalization-group technique for calcu- lating the conductance of strongly correlated nanostructures.Phys. Rev. B73, 195304 (2006)
work page 2006
-
[28]
Ronca, E., Li, Z., Jimenez-Hoyos, C. A. & Chan, G. K.-L. Time-step targeting time-dependent and dynamical density matrix renormalization group algorithms with ab initio hamiltonians.J. Chem. Theory Comput.13, 5560–5571 (2017)
work page 2017
-
[29]
Greene, S. M. & Batista, V. S. Tensor-train split-operator fourier transform (tt-soft) method: Multidimensional nonadiabatic quantum dynamics.J. Chem. Theory Comput.13, 4034–4042 (2017)
work page 2017
-
[30]
Time-dependent variational principle.Int
Moccia, R. Time-dependent variational principle.Int. J. Quantum Chem.7, 779–783 (1973)
work page 1973
-
[31]
Broeckhove, J., Lathouwers, L., Kesteloot, E. & Van Leuven, P. On the equiva- lence of time-dependent variational principles.Chem. Phys. Lett.149, 547–550 (1988)
work page 1988
-
[32]
Holtz, S., Rohwedder, T. & Schneider, R. The alternating linear scheme for tensor optimization in the tensor train format.SIAM J. Sci. Comput.34, A683– A713 (2012)
work page 2012
-
[33]
Lubich, C., Oseledets, I. V. & Vandereycken, B. Time integration of tensor trains.SIAM J. Numer. Anal.53, 917–941 (2015)
work page 2015
-
[34]
Haegeman, J., Lubich, C., Oseledets, I., Vandereycken, B. & Verstraete, F. Uni- fying time evolution and optimization with matrix product states.Phys. Rev. B 94, 165116 (2016)
work page 2016
-
[35]
De Chiara, G., Montangero, S., Calabrese, P. & Fazio, R. Entanglement entropy dynamics of heisenberg chains.J. Stat. Mech: Theory Exp.2006, P03001 (2006)
work page 2006
-
[36]
Bardarson, J. H., Pollmann, F. & Moore, J. E. Unbounded growth of entangle- ment in models of many-body localization.Phys. Rev. Lett.109, 017202 (2012)
work page 2012
-
[37]
Ren, J., Wang, Y., Li, W., Jiang, T. & Shuai, Z. Time-dependent density matrix renormalization group coupled with n-mode representation potentials for the excited state radiationless decay rate: Formalism and application to azulene. Chin. J. Chem. Phys.34, 565–582 (2021)
work page 2021
-
[38]
Wang, Y., Ren, J. & Shuai, Z. Evaluating the anharmonicity contributions to the molecular excited state internal conversion rates with finite temperature TD-DMRG.J. Chem. Phys.154(2021). 17
work page 2021
-
[39]
Seliskar, C. J. & McGlynn, S. P. Electronic spectroscopy of maleimide and its isoelectronic molecules. i. maleimide and n-alkylmaleimides.J. Chem. Phys.55, 4337–4342 (1971)
work page 1971
-
[40]
Lehr, A., G´ omez, S., Parkes, M. A. & Worth, G. A. The role of vibronic cou- pling in the electronic spectroscopy of maleimide: a multi-mode and multi-state quantum dynamics study.Phys. Chem. Chem. Phys.22, 25272–25283 (2020)
work page 2020
-
[41]
Meyer, H.-D., Manthe, U. & Cederbaum, L. S. The multi-configurational time- dependent hartree approach.Chem. Phys. Lett.165, 73–78 (1990)
work page 1990
-
[42]
Manthe, U., Meyer, H.-D. & Cederbaum, L. S. Wave-packet dynamics within the multiconfiguration hartree framework: General aspects and application to nocl.J. Chem. Phys.97, 3199–3213 (1992)
work page 1992
-
[43]
Meyer, H.-D. & Worth, G. A. Quantum molecular dynamics: propagating wavepackets and density operators using the multiconfiguration time-dependent hartree method.Theor. Chem. Acc.109, 251–267 (2003)
work page 2003
-
[44]
Wang, H. & Thoss, M. Multilayer formulation of the multiconfiguration time- dependent hartree theory.J. Chem. Phys.119, 1289–1299 (2003)
work page 2003
-
[45]
Manthe, U. A multilayer multiconfigurational time-dependent hartree approach for quantum dynamics on general potential energy surfaces.J. Chem. Phys.128 (2008)
work page 2008
-
[46]
Tully, J. C. Molecular dynamics with electronic transitions.J. Chem. Phys.93, 1061–1071 (1990)
work page 1990
-
[47]
Basics of surface hopping in mixed quantum/classical simulations
Drukker, K. Basics of surface hopping in mixed quantum/classical simulations. J. Comput. Phys.153, 225–272 (1999)
work page 1999
-
[48]
Tully, J. C. Perspective: Nonadiabatic dynamics theory.J. Chem. Phys.137 (2012)
work page 2012
-
[49]
Boguslawski, K., Tecmer, P., Legeza, O. & Reiher, M. Entanglement measures for single-and multireference correlation effects.J. Phys. Chem. Lett.3, 3129– 3135 (2012)
work page 2012
-
[50]
Boguslawski, K., Tecmer, P., Barcza, G., Legeza, O. & Reiher, M. Orbital entanglement in bond-formation processes.J. Chem. Theory Comput.9, 2959– 2973 (2013)
work page 2013
-
[51]
Stein, C. J. & Reiher, M. Automated selection of active orbital spaces.J. Chem. Theory Comput.12, 1760–1771 (2016)
work page 2016
-
[52]
Glaser, N., Baiardi, A., Lieberherr, A. Z. & Reiher, M. Vibrational entanglement through the lens of quantum information measures.J. Phys. Chem. Lett.15, 6958–6965 (2024). 18
work page 2024
-
[53]
Verstraete, F., Garcia-Ripoll, J. J. & Cirac, J. I. Matrix product density opera- tors: Simulation of finite-temperature and dissipative systems.Phys. Rev. Lett. 93, 207204 (2004)
work page 2004
-
[54]
Feiguin, A. E. & White, S. R. Finite-temperature density matrix renormalization using an enlarged hilbert space.Phys. Rev. B72, 220401 (2005)
work page 2005
-
[55]
Glaser, N., Baiardi, A. & Reiher, M. Flexible DMRG-based framework for an- harmonic vibrational calculations.J. Chem. Theory Comput.19, 9329–9343 (2023)
work page 2023
-
[56]
Rabitz, H. & Ali¸ s, ¨O. F. General foundations of high-dimensional model repre- sentations.J. Math. Chem.25, 197–233 (1999)
work page 1999
-
[57]
Alı¸ s,¨O. F. & Rabitz, H. Efficient implementation of high dimensional model representations.J. Math. Chem.29, 127–142 (2001)
work page 2001
-
[58]
Manzhos, S. & Carrington, T. A random-sampling high dimensional model representation neural network for building potential energy surfaces.J. Chem. Phys.125(2006)
work page 2006
- [59]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.