Lecture 6: MD2 Flashcards
What happens at a boundary?
It can be a repulsive wall or a periodic boundary
What are periodic boundary conditions?
-they mimic bulk conditions
-as particle goes out the box, it’s periodic image comes in the other side
-system must be large enough that a particle doesn’t interact with its mirror image (as least 2x potential cutoff)
Why are periodic boundary conditions good?
Ensemble properties can be calculated with fewer particles
What ensemble do pure molecular dynamics occur in
The micro canonical (NVE) ensemble
What might be a disadvantage of using a microcanonical ensemble?
Not useful to compare with experiments which are usually at a constant temperature
What can be used to control the temp of the system
Thermostats
What is équipartition theory?
KE total=3N/2kB T
What is a drawback of the Velocity Rescaling Thermostat?
It does not generate the canonical ensemble
Consider a limiting case where we rescale the temp every step, which would lead to zero temp fluctuations
How does the Berendsen thermostat work?
Heats system with heat bath
But this suppresses the KE fluctuations not consistent with ensemble
So a heat decay of Tau is added
What is the Nosé-Hoover thermostat?
Introduces thermal reservoir and friction coefficient into equations of motion and so allows KE changes in the system
Allows canonical ensemble correctly
What is a common approach of thermostats in simulations
Use a weak coupling method to equilibrate the system, then use Nosé-Hoover thermostat for the final simulations
What is a barostat?
Couple the system to a pressure bath
What is the berendsen barostat?
Rescaled coordinates every step
Doesn’t generate proper thermodynamic ensemble
What is the Parrinello-Rahman barostat?
-similar to Nosé-Hoover thermostat where particle equations of motion are updated every step
-correct thermodynamic ensemble generated
How do you know when you’re ready to start sampling?
When the system has reached equilibrium