• Chinese Physics B
  • Vol. 29, Issue 10, (2020)
Chuanbiao Zhang1 and Xin Zhou2,†
Author Affiliations
  • 1College of Physics and Electronic Engineering, Heze University, Heze 27405, China
  • 2School of Physical Sciences, University of Chinese Academy of Sciences, Beijing 100190, China
  • show less
    DOI: 10.1088/1674-1056/abad24 Cite this Article
    Chuanbiao Zhang, Xin Zhou. Find slow dynamic modes via analyzing molecular dynamics simulation trajectories[J]. Chinese Physics B, 2020, 29(10): Copy Citation Text show less

    Abstract

    It is a central issue to find the slow dynamic modes of biological macromolecules via analyzing the large-scale data of molecular dynamics simulation (MD). While the MD data are high-dimensional time-successive series involving all-atomic details and sub-picosecond time resolution, a few collective variables which characterizing the motions in longer than nanoseconds are needed to be chosen for an intuitive understanding of the dynamics of the system. The trajectory map (TM) was presented in our previous works to provide an efficient method to find the low-dimensional slow dynamic collective-motion modes from high-dimensional time series. In this paper, we present a more straight understanding about the principle of TM via the slow-mode linear space of the conformational probability distribution functions of MD trajectories and more clearly discuss the relation between the TM and the current other similar methods in finding slow modes.
    γd\boldsymbolq=U\boldsymbolqdt+2kBTγdtd\boldsymbolG,(1)

    View in Article

    tP(q,t)=\boldsymbolLFPP(q,t),(2)

    View in Article

    (1)

    View in Article

    tΨ(q,t)=\boldsymbolHΨ(q,t),(3)

    View in Article

    (1)

    View in Article

    P(q,t)=Peq(q)n=0an(0)eλntΦ^n(q),(4)

    View in Article

    P(q,t)Peq(q)1+n=1Nan(t)Φ^n(q),(5)

    View in Article

    Φ^n(q)=ibn,iPi(q)Peq(q),(6)

    View in Article

    Φ^n(q)=αcnαΘα(q).(7)

    View in Article

    Pi(q)Peq(q)=αciαΘα(q),(8)

    View in Article

    Pi(q)Peq(q)μpi,μAμ(q).(9)

    View in Article

    Aμ(q)Aν(q)eq=δμν,(10)

    View in Article

    \boldsymbolpi=pi,μA^μ(q),(11)

    View in Article

    B^α(q)=bαμA^μ(q),(12)

    View in Article

    B^α(t)=1Δttt+ΔtB^α(q(t))dt.(13)

    View in Article

    G(q,t|q0,0)=Peq(q)[1+n=1NeλntΦ^n(q)Φ^n(q0)],(14)

    View in Article

    Gνμ(t,0)=dqdqAμ(q)G(q,t|q,0)Aν(q)P(q,0)Aμ(q(t))Aν(q(0))0,(15)

    View in Article

    Gνμ(t)=1Z^μθμ(q(t))θν(q(0))eq=nμν(t,0)nν(0),(16)

    View in Article

    (1)

    View in Article

    Σμν=1τ0τdt(1tτ)[Cμν(t)+Cνμ(t)],(17)

    View in Article

    Cμν(t)=1τt0τtdt1Gμν(t+t1,t1).(18)

    View in Article

    Chuanbiao Zhang, Xin Zhou. Find slow dynamic modes via analyzing molecular dynamics simulation trajectories[J]. Chinese Physics B, 2020, 29(10):
    Download Citation