LVCM (Linear vibronic coupling models)

back to the main tutorial

The following figures and data are our benchmarking results (for detailed discussion please see ref1).

 

Figure1

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 R¯(t) of the nuclear normal mode v6a. Panels (c1)-(c4): The average dimensionless momentum P¯(t) of the nuclear normal mode ν6a. 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

 

Figure2

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 R¯(t) of the nuclear normal mode v6a. Panels (c1)-(c4): The average dimensionless momentum P¯(t) of the nuclear normal mode ν6a. 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

 

Figure3

Figure 3. Results of the 2-mode LVCM for the Cr(CO)5 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 R¯(t) 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

 

Figure4

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 nOπ1 state, the ππ1 state and the ππ2 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

 

 

 


1 B. 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).
2 G. 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.
3 G. 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).
4 J. 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).