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

129
votes
3answers
16k views

How does a knife cut things at the atomic level?

As the title says. It is common sense that sharp things cut, but how do they work at the atomical level?
26
votes
6answers
3k views

How can fast moving particles gain energy from slow moving ones?

Imagine a large diameter piston filled with water connected to a small funnel. When you press on the piston slowly but with considerable force the water will move very quickly from the funnel in form ...
26
votes
6answers
5k views

When I stretch a rubber band, it breaks. When I hold the broken ends together, why doesn't it join again?

The question is simple. When we join the two broken surfaces, what is it that keeps the surfaces from connecting with each other, while earlier they were attached to each other? Also, would the two ...
14
votes
5answers
4k views

Does there exist a free good molecule / atom simulation software?

I'm looking for a software or software package (for example C/C++) that can simulate a lot (say thousands at least) of molecules in action (ie. in movement or attached to say static walls). I have ...
14
votes
1answer
577 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 ...
12
votes
1answer
252 views

Do Franck-Condon oscillations have natural lineshapes?

I recently found a paper (for the curious, this one) that talks about observing the motion of a nuclear wavepacket in H2O, as initiated by tunnel ionization. This wavepacket should be thought of as a ...
8
votes
4answers
5k 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., ...
7
votes
1answer
178 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 ...
7
votes
2answers
1k 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? ...
6
votes
1answer
475 views

Born-Oppenheimer separation in Dirac bra-ket notation

Most derivations I have seen of the Born-Oppenheimer approximation are made using wave-functions. To understand it better, I was trying to write a derivation using Dirac notation, but I am stuck. I am ...
6
votes
1answer
69 views

FCC-to-BCC phase transition in NaCl, Buckingham or Lennard-Jones potentials?

Background The transformation from B1 (face centered cubic (FCC) type) to B2 (body centered cubic (body centered cubic (BCC) type) structures is one of the best documented high pressure phase ...
6
votes
0answers
561 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
1k 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 ...
4
votes
1answer
491 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 ...
4
votes
1answer
780 views

Moment of inertia of rotating particles in center of mass frame?

I am trying to simulate a collision between two molecules. I know the energy for every position/orientation, from which I can calculate the forces. The treatment is classical and the molecules are ...
3
votes
2answers
106 views

Gases Expanding and Bonding in the Vacuum? Frozen Clouds

Why do gasses expand in space rather than attract forming a liquid? Do some gasses attract or bond when they get colder and what would be the characteristics of it in space? Would it be like powder ...
3
votes
1answer
158 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 ...
3
votes
2answers
107 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, ...
3
votes
1answer
46 views

Symmetry in program for ewald summation

The formula for Ewald summation as given in Allen and Tildesley - $$ U = U^{(r)} + U^{(k)} + U^{(bc)} + U^{self} $$ where the k-space contribution of potential is given by $$ U^{(k)} = \frac{1}{\pi ...
3
votes
1answer
163 views

Calculation of long-range forces in Molecular Dynamics - Ewald summation

I am trying to write a code to calculate the potential and forces, for the same using ewald summation.For this purpose, the formula for potential and force I have used is : $$ U = U^{(r)} + U^{(k)} ...
3
votes
1answer
3k views

Understanding the mean square displacement in molecular dynamics

In a Molecular Dynamics (MD) simulation, the mean square displacement $\text{MSD}$ is given by $$\text{MSD}(\delta t) = \left\langle\left|\vec{r}(\delta t)-\vec{r}(0)\right|^2\right\rangle,$$ where ...
3
votes
0answers
34 views

Nose-Hoover Barostat

Much can be found about the Nose-Hoover Thermostat. However I seem to be having difficulty finding out details about the Nose-Hoover Barostat, and how it is implemented. Would anyone be able to give ...
3
votes
1answer
237 views

Force field for calcium ions, potassium ions and water in molecular dynamics

I want to adapt my existing MD simulation so that it can handle these three species. For the interactions of the ions I thought that I would use the coulomb potential, but I don't know what sort of ...
2
votes
1answer
168 views

Time length molecular dynamics

What limits time length of simulation in molecular dynamics to femto and in coarse graining to nano or micro second?
2
votes
1answer
245 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 ...
2
votes
1answer
90 views

Why is modelling water so difficult?

There are several models of water that have come into existence so far. Nevertheless, there is no one particular model which effectively gives all the properties of water correctly. A model like TIP3P ...
2
votes
3answers
84 views

Do molecular bounded systems shield or reduce neutron cross-sections?

When talking about neutron cross-sections, literature is usually investigating isolated cases of Neutron + Atom. Here, the abundance of hydrogen is dominating neutron fluxes through material. I ...
2
votes
1answer
106 views

Estimate the difference between two sets of atoms

I've been working on amorphous structures derived from a crystalline one (using MD) containing $N$ atoms. I want to prove that these structures are different and to quantify their "differentness". One ...
2
votes
0answers
38 views

For molecular orbital, LCAO approximation is trusty?

I'm curious how the LCAO approximation is trusty. I'm majoring in physics and usually treat solid using tight binding method. Using tight binding method, We can think the periodic molecular orbital ...
2
votes
0answers
63 views

How can the fictitious mass in the Car-Parrinello method reproduce the “real” dynamics?

In the Car-Parrinello method, to solve simultaneously the classical equations of motion for the atoms and the Kohn-Sham equations for the electrons, the following effective Lagrangian is used: $$ ...
2
votes
1answer
51 views

Mass scaling in molecular dynamics

I've noticed some authors scale the mass of particles in molecular dynamics simulations while leaving the force field parameters the same in order to achieve materials of different densities. Does ...
2
votes
1answer
136 views

Thermionic emission, delayed emission and predissociation

In molecular photodissociation, the thermionic emission, delayed emission and predissociation are the same? Otherwise, what is the difference between them? My question is not about the solids, but I ...
2
votes
0answers
196 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 ...
2
votes
1answer
182 views

Reaction coordinate as a function of atomic positions

I'm going over some (molecular dynamics) related literature - specifically the derivation of the Weighted Histogram Analysis Method (WHAM). As a quick backdrop WHAM is a method for stitching ...
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 ...
1
vote
3answers
138 views

Langevin Equation - Stochastic Differential Equation. What are the subtleties?

I am trying to find out the motion of a particle in 3D governed by the Langevin equation, numerically. Anyway, the Langevin equation is given by $$m \ddot{x} = -(6\pi a\nu) \dot{x} + F_b $$ where ...
1
vote
1answer
32 views

Thermostat in Dissipative Particle Dynamics

I'm doing a computational experiment about dynamics of a solvated line polymer, with dissipative particle dynamics(DPD). However when I try to measure kinetic temperature of the whole system(monomers ...
1
vote
1answer
83 views

Is there any physical simulator precise to atomic levels and quantum effects? [duplicate]

There are tons and tons of physics simulators for classical/large scale physics. I'm interested in a simulator in which you feed an input as atoms and positions and it simulates the evolution of the ...
1
vote
1answer
160 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
3answers
2k views

Franck Condon Principle and Born Oppenheimer approximation

My question here is purely fundamental. I am confused with the concept in Franck Condon (FC) principle and Born Oppenheimer (BO) approximation. The FC principle is in accordance with the BO ...
1
vote
1answer
75 views

Good book for learning fluid statistics

I'm currently using Gray and Gubbins Theory of Molecular Fluids to learn about the statistical physics of fluids. It may be a fine reference text, but I'm not impressed with it as an introduction to ...
1
vote
3answers
61 views

Why is intermolecular potential energy given in 'per mole'?

I was just wondering that if intermolecular potential is the interaction potential energy of two atoms/molecules, then why do we give its value in J/mole? I don't understand why 'per mole' ...
1
vote
1answer
98 views

Molecular simulation applied to astrobiology

There is some research group which use molecular dynamics simulations to study the origin of life according to the hypothesis of abiogenesis? I'm just beginning to know this subject so I don't know ...
1
vote
1answer
425 views

Coherent intermediate scattering function from simulations

I want to calculate the coherent intermediate scattering function, $S_{\text{coh}} (\mathbf{Q},t)$, from a molecular dynamics trajectory, based on its definition. The definition of the function is: ...
1
vote
2answers
193 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 ...
1
vote
0answers
24 views

Atomic Stress (Virial Stress) and Harmonic Springs

The general context here is polymer physics. The stress tensor in a bead-spring simulation can be easily computed using the Virial Stress (which may be named the Kirkwood formula, or the Doi-Edwards ...
1
vote
0answers
19 views

Relation between solvent accessibility and brownian motion

Assume one has a molecule (made of nodes) inside a solvent. If one tries to model the average effect of the interaction between the molecule and the solvent, one has two effects: 1- Friction term on ...
1
vote
1answer
100 views

What actually makes boiling happen?

Boiling can be breezed over easily with a few rudimentary diagrams and a couple equations, but I seek a deeper explanation. The definition of boiling is that the vapor pressure in the liquid is ...
1
vote
1answer
42 views

What is thermophoresis?

I read wikipedia article and I saw a bad youtube presentation on thermophoresis, however I don't have a clear insight about the subject. I assume the forces are basically Brownian molecular forces ...
1
vote
0answers
40 views

Using reflective optics concepts to focus gas pressure in a vacuum

At very low vacuum pressures the bulk behaviour of molecules can be approximated by non interacting particles bouncing off of the chamber walls: For an monotonic molecule such as He with perfect ...