Molecular dynamics nearly integrable system

The study of liquids near solid surfaces using microscopic (atomistic-based) descriptions of liquid molecules is relatively new. Given a potential energy function for the interaction between liquid molecules and between the liquid molecules and the solid surface, the integral equation for the liquid density profile and the liquid molecules orientation can be solved approximately, or the molecular dynamics method can be used to calculate these and many other structural and dynamic properties. In applying these methods to water near a metal surface, care must be taken to include additional features that are unique to this system (see later discussion). 

In a hybrid method, molecules are displaced in time according to conventional molecular dynamics (MD) algorithms, specifically, by integrating Newton s equations of motion for the system of interest. Once the initial coordinates and momenta of the particles are specified, motion is deterministic (i.e., one can determine with machine precision where the system will be in the near future). In the context of Eq. (2.1), the probability of proposing a transition from a state 0 to a state 1 is determined by the probability with which the initial velocities of the particles are assigned from that point on, motion is deterministic (it occurs with probability one). If initial velocities are sampled at random from a Maxwellian distribution at the temperature of interest, then the transition probability function required by Eq. 

---