Q1.

In molecular dynamics simulation of a liquid, explain how the Taylor-series based algorithm can
be used to compute the next step of positions and velocities. Highlight the advantages and
disadvantages of this algorithm.

In molecular dynamics (MD) simulations, the goal is to compute the time evolution of a system of particles based on
Newton's equations of motion. The Taylor-series based algorithm is one such method for propagating the positions and
velocities of particles forward in time. This approach expands particle positions and velocities as Taylor series in time.

Let the position r(t) and velocity v(t) of a particle be known at time t. We can expand them as Taylor series around time t:

In practical implementations, only the first few terms are used due to computational cost and error propagation.
Application in MD Algorithms: This expansion forms the basis of time integration methods like:

Velocity Verlet Algorithm (2nd-order Taylor):

This is derived by truncating the Taylor series at appropriate orders and using symmetries.
•Higher-order algorithms can retain more Taylor terms, giving better accuracy but at greater computational cost.
Advantages of Taylor-Series Based Algorithms

Predictability and Accuracy: When expanded to a suitable order, Taylor series methods can be very accurate for smooth,
conservative force fields.
Systematic Improvement: Accuracy can be improved systematically by including more terms.
Time Reversibility: Symmetric forms (like Velocity Verlet) preserve time-reversal symmetry, important for conservation
laws.

Disadvantages of Taylor-Series Based Algorithms

❖ Computational Cost: Higher-order derivatives (e.g., jerk, snap) require computing derivatives of forces, which is often
non-trivial and expensive.
❖ Stability Issues: For stiff systems or large time steps, these methods may become unstable unless sufficiently high-order
terms are included.
❖ Limited Utility in Complex Systems: For systems with discontinuities (e.g., collisions, reactive events), Taylor expansions
lose validity quickly.
❖ Not Adaptive: Taylor methods assume smooth, deterministic evolution and do not adapt to fast-changing dynamics
without adjustment.

Summary
Taylor-series based algorithms provide a fundamental basis for time integration in molecular dynamics. In practice, the
Velocity Verlet algorithm (a truncated Taylor-series approach) is commonly used due to its balance of efficiency and
accuracy. Full Taylor expansions are more academic than practical for most large-scale or long-time simulations.
Q2. In molecular dynamics simulation of a liquid write the steps involved in general
molecular dynamic simulations.

In a molecular dynamics (MD) simulation of a liquid, the goal is to model the motion of atoms or molecules over time based
on Newtonian mechanics. Here's a general outline of the steps involved in a typical MD simulation:

1. Define the System

•Choose the molecules: Select the type(s) of atoms/molecules (e.g., water, argon).
•Set the number of particles (N).
•Determine the simulation box size and shape, typically with periodic boundary conditions (PBC) to mimic a bulk liquid.

2. Select an Interaction Potential

•Choose a force field or potential function (e.g., Lennard-Jones, TIP3P for water, CHARMM, AMBER) to model interatomic
forces.
•This defines how particles interact (bonded and non-bonded interactions).

3. Initialize Positions and Velocities

•Positions: Usually assigned by:
• Placing atoms on a lattice or
• Randomly placing them (avoiding overlaps).
•Velocities: Initialized from a Maxwell-Boltzmann distribution at the desired temperature.
 4. Apply Boundary Conditions
•Periodic Boundary Conditions (PBC) are applied to simulate an infinite system and reduce surface effects.

5. Compute Initial Forces

•Based on initial positions, compute forces using the chosen potential.

6. Integrate Equations of Motion

•Use an integration algorithm (e.g., Velocity Verlet, Leapfrog, or Runge-Kutta) to update positions and velocities:
• r(t+Δt), v(t+Δt)

• Repeat for each time step.

7. Temperature and Pressure Control (if needed)
•Use thermostats (e.g., Berendsen, Nose-Hoover) to control temperature.
•Use barostats (e.g., Parrinello-Rahman) to control pressure for NPT ensembles.

8. Equilibration Phase
•Allow the system to reach equilibrium (steady temperature, pressure, energy).
•Data from this phase is usually discarded.
 9. Production Run
•Once equilibrated, run the simulation for the required time (e.g., nanoseconds to microseconds).
•Collect data: trajectories, energies, temperature, pressure, radial distribution functions, etc.

10. Analyze Results

•Post-process data to compute physical and structural properties:
• Radial distribution function (RDF)
• Diffusion coefficients
• Energy fluctuations
• Structural snapshots
• Thermodynamic averages

Summary of MD Simulation Workflow:

▪ Define system and model
▪ Initialize positions and velocities
▪ Apply boundary conditions
▪ Compute forces
▪ Integrate equations of motion
▪ Equilibrate the system
▪ Run production simulation
▪ Analyze the data
Q3: In molecular dynamics simulation of a liquid, 2048 molecules are considered in a
simulation box. While doing the simulation, discuss the periodic boundary condition and
minimum image convention.
In molecular dynamics (MD) simulations of liquids, Periodic Boundary Conditions (PBC) and the Minimum Image
Convention (MIC) are essential tools for simulating bulk-like behavior with a finite number of particles. Here's how they
apply when simulating a system of 2048 molecules in a simulation box:

1. Periodic Boundary Conditions (PBC)

Purpose:
•To eliminate edge effects by mimicking an infinite system using a finite simulation box.
•Ensures that molecules leaving one side of the box re-enter from the opposite side, maintaining constant density.

How It Works:
•The simulation box is surrounded by 26 identical copies (in 3D).
•If a molecule moves outside the box through one face, it appears on the opposite face.
•This allows molecules to interact as if they are in a bulk liquid without surface effects.

Example (2048 Molecules):

•Suppose your box is cubic with volume VVV, filled uniformly with 2048 molecules.
•Each time a molecule crosses the boundary at, say, x=Lx = Lx=L, it reappears at x=0x = 0x=0, and vice versa.
•Interactions consider not just the primary box but also the periodic images, which allows simulating continuous, infinite
matter.
 2. Minimum Image Convention (MIC)
Purpose:
•To ensure that each particle interacts only with the closest periodic image of every other particle.
•Reduces computation and avoids redundant force calculations.

How It Works:
•When computing the distance between two particles iii and jjj, the MIC ensures you measure the distance to the nearest
image of jjj relative to iii.
•For a cubic box of side LLL, if the vector between particles is:

•This formula brings the vector into the range [−L/2,L/2] for each coordinate, ensuring it's the shortest possible distance.

Condition:
•For MIC to work correctly, the cutoff radius rcr_crc for interactions must be less than or equal to L/2.
•This ensures no ambiguity in choosing the closest image.
 Combined Effect in MD Simulations:
•Together, PBC and MIC ensure that:
• The system behaves like an infinite, continuous liquid.
• Each particle "sees" its surroundings correctly and efficiently.
• Computational cost is minimized by avoiding unnecessary interaction calculations.

Key Advantages:
•Eliminates boundary artifacts.
•Enables modeling of realistic liquid behavior with a limited number of molecules (like 2048).
•Ensures computational efficiency and physical accuracy.

Consider box of size L, repeat box infinitely many times in all directions

• Each particle interacts (in principle) with all particles in all

boxes → problems for long-range interactions (infinite
re-summation necessary)

• short-range interactions: minimum image convention: at

most the closest of all images of a particle j can interact with
a given particle i → great simplification: pick the closest
image and use this to calculate V(rij)