8.1 Quantum Mechanics and Molecular Dynamics

Quantum mechanics forms the foundation of molecular modeling, describing matter's behavior at atomic scales. It introduces wave-particle duality and the uncertainty principle, with the Schrödinger equation as its cornerstone. These concepts are crucial for understanding electronic structures and potential energy surfaces in molecules.

The Born-Oppenheimer approximation simplifies molecular calculations by separating electronic and nuclear motions. This approach allows for the creation of potential energy surfaces, which guide our understanding of molecular geometry, vibrations, and chemical reactions. These principles are essential for accurate molecular simulations.

Quantum Mechanics in Molecular Modeling

Fundamentals of quantum mechanics

Top images from around the web for Fundamentals of quantum mechanics

• 

  Schrödinger Equation [The Physics Travel Guide] View original

  Is this image relevant?

• 

  quantum mechanics - Schrodinger equation to find general wave function - Mathematics Stack Exchange View original

  Is this image relevant?

• 

  quantum mechanics - How to "read" Schrodinger's equation? - Physics Stack Exchange View original

  Is this image relevant?

• 

  Schrödinger Equation [The Physics Travel Guide] View original

  Is this image relevant?

• 

  quantum mechanics - Schrodinger equation to find general wave function - Mathematics Stack Exchange View original

  Is this image relevant?

1 of 3

Top images from around the web for Fundamentals of quantum mechanics

• 

  Schrödinger Equation [The Physics Travel Guide] View original

  Is this image relevant?

• 

  quantum mechanics - Schrodinger equation to find general wave function - Mathematics Stack Exchange View original

  Is this image relevant?

• 

  quantum mechanics - How to "read" Schrodinger's equation? - Physics Stack Exchange View original

  Is this image relevant?

• 

  Schrödinger Equation [The Physics Travel Guide] View original

  Is this image relevant?

• 

  quantum mechanics - Schrodinger equation to find general wave function - Mathematics Stack Exchange View original

  Is this image relevant?

1 of 3

• Quantum mechanics describes the behavior of matter at the atomic and subatomic scales
  • Wave-particle duality: Particles exhibit both wave-like (interference, diffraction) and particle-like (discrete energy levels) properties
  • Uncertainty principle: The position and momentum of a particle cannot be simultaneously determined with arbitrary precision (Heisenberg's uncertainty principle)
• Schrödinger equation: Fundamental equation in quantum mechanics that describes the time evolution of a quantum system
  • Time-dependent Schrödinger equation: $i\hbar\frac{\partial}{\partial t}\Psi(\mathbf{r},t) = \hat{H}\Psi(\mathbf{r},t)$
  • Time-independent Schrödinger equation: $\hat{H}\Psi(\mathbf{r}) = E\Psi(\mathbf{r})$
  • $\hat{H}$: Hamiltonian operator, representing the total energy of the system (kinetic and potential energy)
  • $\Psi(\mathbf{r},t)$: Wave function, containing all information about the quantum system (probability amplitude)
  • $E$: Energy eigenvalue, the allowed energy levels of the system
• Applying quantum mechanics to molecular modeling
  • Electronic structure calculations: Solving the Schrödinger equation for electrons in a molecule to determine molecular properties (energy, geometry, dipole moment)
  • Potential energy surfaces: Mapping the energy of a molecular system as a function of its geometry (bond lengths, angles, dihedral angles)
  • Quantum chemistry methods: Ab initio (Hartree-Fock, coupled cluster), density functional theory (DFT), and semi-empirical methods (AM1, PM3) used to solve the electronic Schrödinger equation

Born-Oppenheimer approximation in molecules

• Born-Oppenheimer approximation: Decouples electronic and nuclear motions in molecules based on the large difference in their masses
  • Assumes that electrons adjust instantaneously to changes in nuclear positions due to their much smaller mass (1800 times lighter than protons)
• Electronic Schrödinger equation: $\hat{H}_\text{elec}\Psi_\text{elec}(\mathbf{r};\mathbf{R}) = E_\text{elec}(\mathbf{R})\Psi_\text{elec}(\mathbf{r};\mathbf{R})$
  • $\hat{H}_\text{elec}$: Electronic Hamiltonian, including kinetic energy of electrons and potential energy from electron-electron and electron-nucleus interactions
  • $\Psi_\text{elec}(\mathbf{r};\mathbf{R})$: Electronic wave function, dependent on electronic coordinates $\mathbf{r}$ and parametrically dependent on nuclear coordinates $\mathbf{R}$
  • $E_\text{elec}(\mathbf{R})$: Electronic energy, a function of nuclear coordinates, represents the potential energy surface for nuclear motion
• Potential energy surface: $E_\text{tot}(\mathbf{R}) = E_\text{elec}(\mathbf{R}) + V_\text{nn}(\mathbf{R})$
  • $V_\text{nn}(\mathbf{R})$: Nuclear-nuclear repulsion energy, a function of nuclear coordinates
  • $E_\text{tot}(\mathbf{R})$: Total potential energy surface, governing the motion of nuclei (vibrations, rotations, conformational changes)
• Adiabatic approximation: Assumes that the system remains in the same electronic state during nuclear motion, valid for most ground-state processes (exceptions: photochemistry, electron transfer)

Molecular Dynamics Simulations

Implementation of molecular dynamics simulations

• Molecular dynamics (MD) simulations: Computational method to study the time-dependent behavior of molecular systems by numerically solving Newton's equations of motion for a system of interacting particles
  • Classical mechanics: Particles obey Newton's laws of motion ($\mathbf{F} = m\mathbf{a}$)
  • Quantum effects: Usually neglected, but can be included through ab initio MD or quantum-classical hybrid methods (QM/MM)
• Force field: Mathematical description of the potential energy of a system as a function of particle positions
  • Bonded interactions: Bond stretching (harmonic potential), angle bending (harmonic potential), and torsional terms (cosine series)
  • Non-bonded interactions: Van der Waals (Lennard-Jones potential) and electrostatic (Coulomb potential) terms
  • Parametrization: Force field parameters (force constants, equilibrium values) obtained from experimental data or quantum chemical calculations
• Numerical integration: Algorithms for propagating particle positions and velocities over time
  • Verlet algorithm: $\mathbf{r}(t+\Delta t) = 2\mathbf{r}(t) - \mathbf{r}(t-\Delta t) + \frac{\mathbf{F}(t)}{m}\Delta t^2$
    1. Calculate new positions using current positions, previous positions, and forces
    2. Update forces using new positions
    3. Repeat steps 1-2 for the desired number of time steps
  • Velocity Verlet algorithm: $\mathbf{r}(t+\Delta t) = \mathbf{r}(t) + \mathbf{v}(t)\Delta t + \frac{1}{2}\frac{\mathbf{F}(t)}{m}\Delta t^2$ and $\mathbf{v}(t+\Delta t) = \mathbf{v}(t) + \frac{1}{2}\left(\frac{\mathbf{F}(t)}{m} + \frac{\mathbf{F}(t+\Delta t)}{m}\right)\Delta t$
    1. Calculate new positions using current positions, velocities, and forces
    2. Update forces using new positions
    3. Calculate new velocities using current velocities and the average of old and new forces
    4. Repeat steps 1-3 for the desired number of time steps
• Periodic boundary conditions: Simulate bulk properties by replicating the simulation box in all directions, eliminating surface effects
• Thermostats and barostats: Control temperature and pressure in MD simulations to sample different statistical ensembles
  • Nosé-Hoover thermostat: Introduces a fictitious heat bath with a coupling parameter to maintain constant temperature
  • Parrinello-Rahman barostat: Allows the simulation box to change shape and size to maintain constant pressure, coupled with a thermostat

Analysis of molecular dynamics results

• Ensemble averages: Calculate macroscopic properties as averages over the microscopic states sampled during the simulation
  • Microcanonical ensemble (NVE): Constant number of particles, volume, and energy, corresponding to an isolated system
  • Canonical ensemble (NVT): Constant number of particles, volume, and temperature, corresponding to a system in contact with a heat bath
  • Isothermal-isobaric ensemble (NPT): Constant number of particles, pressure, and temperature, corresponding to a system in contact with a heat bath and a pressure bath
• Radial distribution function (RDF): Probability of finding a particle at a distance $r$ from another particle, normalized by the average density
  • Provides insight into the local structure and ordering of the system (coordination numbers, solvation shells)
  • Can be compared with experimental data from X-ray or neutron scattering
• Mean square displacement (MSD): Average squared distance traveled by particles over time, $\text{MSD}(t) = \langle |\mathbf{r}(t) - \mathbf{r}(0)|^2 \rangle$
  • Relates to diffusion coefficient $D$ through the Einstein relation: $\text{MSD}(t) = 6Dt$ for 3D systems
  • Can be used to study transport properties (ionic conductivity, viscosity) and phase transitions (solid-liquid, glass transition)
• Velocity autocorrelation function (VACF): Correlation between particle velocities at different times, $\text{VACF}(t) = \langle \mathbf{v}(0) \cdot \mathbf{v}(t) \rangle$
  • Provides information about the dynamics and vibrational spectrum of the system (power spectrum, density of states)
  • Related to the diffusion coefficient through the Green-Kubo relation: $D = \frac{1}{3} \int_0^\infty \text{VACF}(t) dt$
• Free energy calculations: Techniques to compute free energy differences and barriers between different states or along a reaction coordinate
  • Thermodynamic integration: $\Delta A = \int_0^1 \left\langle \frac{\partial U(\lambda)}{\partial \lambda} \right\rangle_\lambda d\lambda$, where $\lambda$ is a coupling parameter that connects the initial and final states
  • Umbrella sampling: Applies a biasing potential to sample high-energy regions of the configuration space, with subsequent unbiasing to obtain the true free energy profile