The following figures and data are our benchmarking results (for detailed discussion please see ref1).
Figure 1. Results of the 24-mode LVCM for pyrazine. Panels (a1)-(a4): Population of the state 2. Panels (b1)-(b4): The average dimensionless coordinate of the nuclear normal mode v6a. Panels (c1)-(c4): The average dimensionless momentum of the nuclear normal mode . In panels (a1), (b1) and (c1), the green, blue and cyan solid lines represent the results of NaF-cc, NaF-cx and NaF-TW, respectively. In panels (a2), (b2) and (c2), the cyan, pink, orange and magenta solid lines represent the results of NaF-TW, NaF-TW2, NaF-HW and NaF-HW2, respectively. In panels (a3), (b3) and (c3), the orange dashed line, purple dashed line, brown dashed line and cyan solid line denote the results of SH-1, SH-2, SH-3 and NaF-TW, respectively. In panels (a4), (b4) and (c4), the blue dashed line, green dashed line and cyan solid line denote the results of CMM ( = 0.366), CMM ( = 0.5) and NaF-TW, respectively. The numerically exact results produced by MCTDH2 are demonstrated by black solid lines with black points in each panel. data
Figure 2. Results of the 3-mode LVCM for pyrazine. Panels (a1)-(a4): Population of the state 2. Panels (b1)-(b4): The average dimensionless coordinate of the nuclear normal mode v6a. Panels (c1)-(c4): The average dimensionless momentum of the nuclear normal mode . In panels (a1), (b1) and (c1), the green, blue and cyan solid lines represent the results of NaF-cc, NaF-cx and NaF-TW, respectively. In panels (a2), (b2) and (c2), the cyan, pink, orange and magenta solid lines represent the results of NaF-TW, NaF-TW2, NaF-HW and NaF-HW2, respectively. In panels (a3), (b3) and (c3), the orange dashed line, purple dashed line, brown dashed line and cyan solid line denote the results of SH-1, SH-2, SH-3 and NaF-TW, respectively. In panels (a4), (b4) and (c4), the blue dashed line, green dashed line and cyan solid line denote the results of CMM ( = 0.366), CMM ( = 0.5) and NaF-TW, respectively. The numerically exact results produced by MCTDH2 are demonstrated by black solid lines with black points in each panel. data
Figure 3. Results of the 2-mode LVCM for the molecule. Panel (a1)-(a10): Population dynamics, where the black, red and blue solid lines represent the population of states 1−3, respectively, and the numerically exact results (taken from ref 3) are plotted by dashed lines with corresponding colors. Panel (a1): NaF-cc. Panel (a2): NaF-cx. Panel (a3): NaF-TW. Panel (a4): NaF-TW2. Panel (a5): NaF-HW. Panel (a6): NaF-HW2. Panel (a7): SH-1. Panel (a8): SH-2. Panel (a9): CMM ( = 0.333). Panel (a10): CMM ( = 0.5). Panels (b1)-(b4): The average dimensionless coordinate of the second nuclear normal mode. In panel (b1), the green, blue and cyan solid lines represent the results of NaF-cc, NaF-cx and NaF-TW, respectively. In panel (b2), the cyan, pink, orange and magenta solid lines represent the results of NaF-TW, NaFTW2, NaF-HW and NaF-HW2, respectively. In panel (b3), the orange dashed line, purple dashed line and cyan solid line denote the results of SH-1, SH-2 and NaF-TW, respectively. In panel (b4), the blue dashed line, green dashed line and cyan solid line denote the results of CMM ( = 0.333), CMM ( = 0.5) and NaF-TW, respectively. Note that SH-3 is not applicable for this 3-state model. The numerically exact results produced by MCTDH (taken from ref 3) are demonstrated by black solid lines with black points in panels (b1)-(b4). data
Figure 4. Results of the 39-mode LVCM for Thymine parametrized from the CAM-B3LYP calculations. In each panel, population dynamics of 3 states are demonstrated, where the blue, green and orange solid lines represent the population of the state, the state and the state, respectively. Panel (a): NaF-cc. Panel (b): NaF-cx. Panel (c): NaF-TW. Panel (d): NaF-TW2. Panel (e): NaF-HW. Panel (f): NaF-HW2. Panel (g): SH-1. Panel (h): SH-2. Panel (i): CMM ( = 0.261). Panel (j): CMM ( = 0.5). Note that SH-3 is not applicable for this 7-state model. The numerically exact results produced by ML-MCTDH (taken from ref 4) are demonstrated by dashed lines with corresponding colors in each panel. data
1B. Wu, B. Li, X. He, X. Cheng, J. Ren and J. Liu, "Nonadiabatic field: A conceptually novel approach for nonadiabatic quantum molecular dynamics ", J. Chem. Theory Comput. 21, 3775-3813 (2025).↩
2G. A. Worth, M. H. Beck, A. Jackle and H.-D. Meyer, The MCTDH Package, Version 8.2, (2000). H.-D. Meyer Version 8.3 (2002), Version 8.4 (2007). O. Vendrell and H.-D. Meyer Version 8.5 (2013). Version 8.5 contains the ML-MCTDH algorithm. See http://mctdh.uni-hd.de. (accessed on November 1st, 2023) Used version: 8.5.14.↩↩
3G. A. Worth, G. Welch, and M. J. Paterson, “Wavepacket dynamics study of Cr(CO)5 after formation by photodissociation: Relaxation through an (E ⊕ A) ⊗ e Jahn–Teller conical intersection”, Mol. Phys. 104, 1095–1105 (2006).↩↩
4J. A. Green, M. Y. Jouybari, D. Aranda, R. Improta, and F. Santoro, “Nonadiabatic absorption spectra and ultrafast dynamics of DNA and RNA photoexcited nucleobases”, Molecules 26, 1743 (2021).↩