M10 Angular momentum is conserved or not

Manual 10 for MD++
Angular momentum is conserved or not?

Keonwook Kang and Wei Cai

Jul 06 , 2009

---

In molecular dynamics simulations, positions and velocities of atoms are updated following Newton's equations of motion. Naturally, we would expect that linear and angular momenta are conserved during the time integration, or

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \begin{array}{rcl} \mathbf{L} & \equiv & \sum_i m_i\mathbf{v}_ i = \mathrm{Const.} \\ \mathbf{P} & \equiv & \sum_i \mathbf{r}_i \times m_i\mathbf{v}_i = \mathrm{Const.} \end{array}, }

where L is linear momentum and P is angular mumentum. Usually, we subtract ceter-of-mass velocity from the velocity of each atom so that the whole system can not drift and the linear momentume becomes zero.

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \begin{array}{rcl} \mathbf{v}_i & := & \mathbf{v}_ i - \mathbf{v}_{\mathrm{CM}} \\ \mathbf{L} & = & \sum_i m_i\mathbf{v}_i - \mathbf{v}_{\mathrm{CM}}\sum_i m_i= 0 \end{array}, }

where v_CM is defined as

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{v}_{\mathrm{CM}} \equiv \frac{\sum_i m_i \mathbf{v}_i}{\sum_i m_i} } .

Similarly, we can additionaly subtract velocity component contributing rotation so that the whole system can not rotate and the angular momentum becomes zero.

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{v}_i := \mathbf{v}_ i - \mathbf{\Omega}\times\mathbf{r}_i }