Related Experiment Video
Updated: Jun 22, 2026
10:11
Temperature-Controlled Assembly and Characterization of a Droplet Interface Bilayer
Published on: April 19, 2021
Langevin thermostat for rigid body dynamics
Ruslan L Davidchack1, Richard Handel, M V Tretyakov
1Department of Mathematics, University of Leicester, Leicester LE1 7RH, United Kingdom. rld8@mcs.le.ac.uk
The Journal of Chemical Physics
|June 25, 2009
Summary
We developed a new simulation method for rigid molecules using quaternion and Langevin dynamics. This approach accurately samples molecular distributions for improved isothermal simulations.
Area of Science:
- Computational chemistry
- Molecular dynamics
- Statistical mechanics
Background:
- Isothermal simulations are crucial for understanding molecular behavior.
- Accurate sampling of canonical distribution is essential for rigid molecules.
- Existing methods may have limitations in efficiency or accuracy.
Purpose of the Study:
- To introduce a novel method for isothermal rigid body simulations.
- To enhance the sampling of canonical distribution for rigid molecules.
- To evaluate the performance of proposed numerical integrators.
Main Methods:
- Utilizing the quaternion representation for rotational degrees of freedom.
- Employing Langevin dynamics for rigid body simulations.
- Developing and testing quasisymplectic second-order numerical integrators.
- Applying the method to the TIP4P water model.
Main Results:
- The proposed method correctly samples the canonical distribution.
- Quasisymplectic integrators demonstrate good performance.
- Investigated optimal thermostat parameters for simulations.
Conclusions:
- The new quaternion-based Langevin dynamics method is effective for isothermal rigid body simulations.
- The developed numerical integrators offer a simple and accurate approach.
- The study provides insights into optimizing simulation parameters.
Related Concept Videos
01:21Rigid Body Equilibrium Problems - II
A rigid body is in static equilibrium when the net force and the net torque acting on the system are equal to zero.
Consider two children sitting on a seesaw, which has negligible mass. The first child has a mass (m1) of 26 kg and sits at point A, which is 1.6 meters (r1) from the pivot point B; the second child has a mass (m2) of 32 kg and sits at point C. How far from the pivot point B should the second child sit (r2) to balance the seesaw?
Consider two children sitting on a seesaw, which has negligible mass. The first child has a mass (m1) of 26 kg and sits at point A, which is 1.6 meters (r1) from the pivot point B; the second child has a mass (m2) of 32 kg and sits at point C. How far from the pivot point B should the second child sit (r2) to balance the seesaw?
00:49Rigid Body Equilibrium Problems - I
A rigid body is said to be in static equilibrium when the net force and the net torque acting on the system is equal to zero. To solve for rigid body equilibrium problems, do the following steps.
01:11Angular Momentum: Rigid Body
The total angular momentum of a rigid body can be calculated using the summation of the angular momentum of all the tiny particles rotating in the same plane. Considering all the tiny particles rotating in the x-y plane, the direction of angular momentum of all such particles and that of the rigid body would be perpendicular to the plane of the rotation along the z-axis.
This calculation can get complicated when tiny particles within the rigid body are not rotating in the same plane but have...
This calculation can get complicated when tiny particles within the rigid body are not rotating in the same plane but have...
01:12Equation of Motion for a Rigid Body
The movement of a rigid object can be understood through the equations that explain both translational and rotational motion about the center of mass of the object, point G. This center of mass is the point where the equation of motion for translational motion comes into play, as per Newton's Second Law.
The combined moments generated about the center of mass of the object are equal to the rate of change of the angular momentum of the body. An external force, when applied at a different point...
The combined moments generated about the center of mass of the object are equal to the rate of change of the angular momentum of the body. An external force, when applied at a different point...
01:13Kinetic Energy for a Rigid Body
Imagine a solid object involved in a general planar movement, with its center of mass pinpointed at a spot labeled G. The object's kinetic energy relative to an arbitrary point A can be quantified for each of its particles - the ith particle in this case. This measurement is achieved through the employment of the relative velocity definition. The position vector, known as rA, extends from point A to the mass element i.
01:06Virtual Work for a System of Connected Rigid Bodies
Virtual work is a powerful method used to solve problems involving several connected rigid bodies. When the system is in equilibrium, virtual work is zero. This allows the calculation of the resulting forces when a system undergoes a virtual displacement. When attempting to analyze such a system, first, use a free-body diagram, where an independent coordinate represents the configuration of the links, and mark its deflected position resulting from the positive virtual displacement.
Next,...
Next,...