Rigid top propagator
Jump to navigation
Jump to search
For a rigid three dimensional asymmetric top the kernel is given by ([1][2]) (Eq. 15) ):
The contribution to the rotational energy of the interactions between beads Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle t}
and Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle t+1}
is given by (Eq. 16):
- Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle e_{rot}^{t,t+1}={\frac {1}{\rho _{\mathrm {rot} }^{t,t+1}}}\sum _{JM{\hat {K}}}\left({\frac {2J+1}{8\pi ^{2}}}\right)A_{{\hat {K}}M}^{(JM)}E_{\hat {K}}^{(JM)}\exp \left(-{\frac {\beta }{P}}E_{\hat {K}}^{(JM)}\right)\sum _{K}A_{{\hat {K}}K}^{(JM)}d_{MK}^{J}({\tilde {\theta }}^{t+1})\cos(M{\tilde {\phi }}^{t+1}+K{\tilde {\chi }}^{t+1})}
References[edit]
- ↑ M. H. Müser and B. J. Berne "Path-Integral Monte Carlo Scheme for Rigid Tops: Application to the Quantum Rotator Phase Transition in Solid Methane", Physical Review Letters 77 pp. 2638-2641 (1996)
- ↑ Eva G. Noya, Carlos Vega, and Carl McBride "A quantum propagator for path-integral simulations of rigid molecules", Journal of Chemical Physics 134 054117 (2011)
- Related reading