Computational materials science

In molecular dynamics, you want to integrate the equations of motion with a time step of 2 fs to calculate atomic motion over 1 ns. How many time steps are needed?

  1. A5×10²
  2. B5×10⁸
  3. C5×10⁵
Show answer and explanation

Answer: C. 5×10⁵

Since 1 ns = 10⁻⁹ s and 2 fs = 2×10⁻¹⁵ s, the number of steps is 10⁻⁹/(2×10⁻¹⁵) = 5×10⁵. The time step must be small enough to divide the period of atomic vibration (around 100 fs) into several tens of steps, and for systems containing light hydrogen it is set to 1 fs or less. If the step is too large, errors in numerical integration such as the Verlet method destroy energy conservation and the calculation breaks down. 5×10² takes fs as 10⁻¹² s (ps), and 5×10⁸ takes ns as 10⁻⁶ s (µs). Because forces on all atoms are calculated at every step, nanoseconds to microseconds is the rough limit of the time that molecular dynamics can follow.

Report an error

Practice this topic Back to topic

More Computational materials science questions

  1. Which statement about the interatomic potentials used in molecular dynamics is correct?
  2. Argon is described by the Lennard-Jones potential V(r) = 4ε[(σ/r)¹² − (σ/r)⁶] with σ = 0.34 nm. Which is closest to the interatomic distance r₀ at which the energy is a minimum?
  3. In a molecular dynamics calculation of a monatomic system, the average kinetic energy per atom was 0.0517 eV. Which is closest to the temperature of the system obtained from the equipartition theorem? Take the Boltzmann constant as k = 8.617×10⁻⁵ eV/K.
  4. In a Monte Carlo calculation by the Metropolis method, a trial move of an atom raises the energy of the system by 0.10 eV. At a temperature of 1000 K, which is closest to the probability of accepting this move? Take k = 8.617×10⁻⁵ eV/K.
  5. Which statement about the kinetic Monte Carlo (kMC) method is correct?