Molecular dynamics (MD) is an approach used for simulating dynamic processes through successive integration of Newton's equations of motion.

learn more… | top users | synonyms

6
votes
4answers
4k views

Why is the canonical ($NVT$) ensemble often used for (classical) molecular dynamics (MD) simulations?

Molecular dynamics (MD) simulation is a common approach to the (classical) many-body problem. It relies on integration of Newton's equations of motion to simulate the trajectories of many (e.g., ...
0
votes
2answers
183 views

Molecular mean free path probability

Let's let $p(\xi)$ be the probability a molecule travels at least $\xi$ between collisions, lets say $\xi=0.01$. When I think of this statement, I think it it is the probability the molecule is able ...
1
vote
0answers
123 views

Correct way to include constant external force in virial and pressure calculation

Halo, given a simulation cell with N particles where particles interact only with bond and pair potentials and periodic boundary conditions (minimum image convention) are used. On a subgroup of ...
0
votes
2answers
2k views

Does the Lennard-Jones force equation give its answer in Newtons?

I'm trying to do the dimensional analysis of the Lennard-Jones force to work out what units are being used in my MD simulation. The lennard Jones force is given as the negative derivative of the ...
7
votes
2answers
926 views

What is the “direction” of the transition dipole moment? (Understanding Eq. 9.29, Charge and Energy Transfer 3rd Ed, May & Kuhn)

For a real vector $\mathbf{r}$, the direction is given by: $\hat{\mathbf{n}}=\mathbf{r}/\left|\mathbf{r}\right|$. The transition dipole moment is a complex vector. How do you define its direction? ...
4
votes
1answer
338 views

Molecular dynamics: Thermostat via damping modulation

Molecular dynamics (MD) often use a thermostat so that they can control the temperature. Temperature tends to drift due to endo/exothermic reactions (that would normally absorb/emit heat to a large ...
1
vote
0answers
82 views

Any quadrupole approximation? Any example?

In atomic and molecular physics we quite often encounter with electric dipole approximation. The dipole approximation we do when the wave-length of the type of electromagnetic radiation which induces, ...
1
vote
2answers
161 views

Molecular Hamiltonian

I was reading some material on the Molecular Hamiltonian on Wiki. It said that, Almost all calculations of molecular wavefunctions are based on the separation of the Coulomb Hamiltonian first ...
0
votes
0answers
67 views

Can peridynamics give accurate quantitative results?

I'm applying for a position that is interested in studying crack growth using peridynamics. I have experience with MD on polymeric systems, but I had never knew of this method before. It seems very ...
6
votes
0answers
427 views

Surface tension between water (liquid) methane (gas)

In my project I should calculate surface tension between water and methane with lammps. I write the program.I use SPC model for water and $l_j$ for methane and water - methane. But now I don't know ...
5
votes
2answers
806 views

Deuterium density in seawater

I heard today that deuterium for nuclear fusion could be extracted from the Oceans, in which it is present in "heavy water." I heard this claim: "Deuterium is distributed uniformly with Ocean ...
0
votes
0answers
77 views

What is the physical property of metal nanoparticles?

I am a Math student but now I have to deal with gold nanoparticles in aqueous solution. Now I was wondering whether the physical properties of gold nanoparticles are the same as the properties of gold ...
3
votes
1answer
144 views

When is classical mechanics valid for describing motion of atoms?

In Molecular Dynamics simulations, the Newton's equation of motion is used to calculate the time evolution of system. Once, I read in an introductory text that when the thermal de Broglie wavelength ...
0
votes
0answers
99 views

Molecular dynamics simulation: fluctuating dipole model implementation

I'm conducting a molecular dynamics simulation for silica. Some time ago I turned to the fluctuating dipole model, and after much effort I'm still having problems implementing it. In short, all ...
2
votes
1answer
188 views

Is there a typo in this modified Lennard-Jones potential?

The standard 12-6 Lennard Jones potential is given by $$U(r_ij) = 4\epsilon\left[ \left(\frac{\sigma_{ij}}{r_{ij}}\right)^{12} - \left(\frac{\sigma_{ij}}{r_{ij}}\right)^{6} \right]$$ where ...
1
vote
1answer
136 views

Deriving the change in the Helmholtz free energy in the context of the free energy perturbation method

I am reading Free Energy Calculations: Theory and Applications in Chemistry and Biology by Chipot and Pohorille. At the beginning of the text (page 19, for example), the authors define the Helmholtz ...
1
vote
2answers
2k views

Molecular Dynamics (MD) Simulation: energy fluctuations in NVE ensemble

I'm writing my first MD simulation (ever) for liquid Argon. The code is up and running. I am supposed to do the calculations in the NVE ensemble. Having implemented a 4th order symplectic integrator ...
2
votes
0answers
132 views

Phase transition through Molecular Dynamics Studies

I am doing non-equilibrium molecular dynamics to see the phase transition of single crystal metals in extreme condition (shock compression). I am using NVE to conserve total energy, every thing going ...
0
votes
1answer
253 views

Thermionic emission and delayed emission

I want to understand the concepts behind the thermionic emission. In thermionic emission, the energy randomization occurs and the energy may be split to electronic or roto-vibrational states. If this ...
7
votes
1answer
162 views

How does one geometrically quantize the Bloch equations?

I've just now rated David Bar Moshe's post (below) as an "answer", for which appreciation and thanks are given. Nonetheless there's more to be said, and in hopes of stimulating further posts, I've ...
14
votes
1answer
490 views

Why isn't the Gear predictor-corrector algorithm for integration of the equations of motion symplectic?

Okumura et al., J. Chem. Phys. 2007 states that the Gear predictor-corrector integration scheme, used in particular in some molecular dynamics packages for the dynamics of rigid bodies using ...