根据粒子运动规律，对系统每个粒子运动进行数值求解，得到每个粒子相空间的运动轨迹，再利用统计方法得到系统的静态和动态特征。

• 多粒子系统的时间演化动力学行为
• 相空间采样技术→ 观测量统计值

Program MD
Call init
t = 0
do while (t < T_max)
	call force(f, en)
	call integrate(f, en)
	t = t + dt
	call sample
enddo
end MD

• 1 设定物理模型
• 2 计算力和能量
• 3 积分运动方程
• 4 计算物理性质

边界条件

积分运动方程

力场

Ex 等温过程