Questions tagged [molecular-dynamics]
Molecular dynamics simulations aim at solving the equations of motion of the atoms belonging to a molecular system using explicit time propagation and taking into account the effect of temperature on the atomic velocities (as opposite to geometry optimization). The interatomic forces can be described in classical (MM), quantum (QM) or mixed (QM/MM) frameworks.
101
questions
0
votes
0
answers
22
views
why does the average temperature vary a lot from my setting temperature in nve?
I set my t_start and t_stop as 0.52 (lj unit), and I ran the nvt first. after it reached equilibration I ran the nve but found the temperature of the last 1000steps went up to 0.56 while my coworker ...
- 1
2
votes
0
answers
55
views
Symplectic (or alike) integrator for system with Coulomb singularity and time-dependent potentials
I am trying to calculate classical trajectories for a single positive ion and a single electron inside an RF trap. Therefore, I am dealing with a two-body system that possesses:
Coulomb potential ...
- 21
0
votes
0
answers
21
views
Calculating radial distribution function in a cylindrical pore
I have tried calculating radial distribution functions in a cylindrical pore geometry. The pore has z-axis periodicity and it does not have x and y periodicity. But I got some troubles on that. The ...
0
votes
0
answers
25
views
Unit comparison between two molecular dynamics simulations with different time step
I need some clarification regarding the time step used in the molecular dynamics simulations. In the first simulation the time step was 1 nanosecond and in the second simulation the time step in 10 ns....
- 11
0
votes
1
answer
376
views
Who can explain the minimum image convention in molecular dynamic simulations?
How to choose the cutoff radius so that the atoms do not interact with its periodic image? Especially when simulating macromolecules (proteins).
0
votes
1
answer
109
views
How error scales with numerical precision in molecular dynamics?
In terms of time-step, numerical error in molecular dynamics scales with square, i.e. $error \approx dt^2$. But how it look for numerical precision ? E.g. how much bigger will be numerical error when ...
- 123
0
votes
0
answers
55
views
Comparison of computational complexities of MD versus MC simulations
In my humble understanding MD simulations of systems with short-range(like LJ interactions) and long-range(electrostatic) has a computational complexity $O(N . log(N))$. What will be the computational ...
- 81
2
votes
0
answers
31
views
Neighbor lists in overdamped Langevin MD simulations
The following question is about the implementation of molecular dynamics simulations with neighbor lists. In cases where the paths of individual particles are (i) stochastic and (ii) somewhat rough, ...
- 161
4
votes
2
answers
1k
views
Time Reversibility of Velocity Verlet Algorithm
I'm very new to computational Physics and am finding conflicting statements on whether the velocity Verlet algorithm, defined as:
$\begin{align}
x_{n+1} &= x_n + v_n \Delta t + \frac{1}{2} a_n \...
- 141
0
votes
0
answers
42
views
Reference request: Textbook similar in structure to "Computational Materials Science: An Introduction" by June Lee for Computational Biology
In the Computational Materials Science: An Introduction by June Lee, he discusses molecular dynamics and density functional theory with examples from LAMMPS and QuantumEspresso, and explains LAMMPS ...
- 1
1
vote
1
answer
70
views
Which are the right configurations in the Markov chain of a Hamiltonian Monte Carlo algorithm?
I have a question about the Markov Chain Hamiltonian Monte Carlo (MCHMC).
Hamiltonian Monte Carlo is known as Hybrid Monte Carlo too. I'll describe the steps of the algorithm.
We have at the ...
- 113
1
vote
0
answers
34
views
Tackling multiscale problem in numerical simulation
In a dusty plasma system there are more than one component with different masses, i.e, electrons, ions,neutrals and dust grains. Accordingly, there are more than one temporal and spatial scales ...
- 19
0
votes
0
answers
74
views
Why this LJ molecular dynamics result doesn't converge?
I am doing a molecular dynamics simulation of Leonard Jones 6-12 potential. But instead of converging to a particular value. It always stays between -5.58 to -5.62. The standard value is -5.517. The ...
- 27
4
votes
1
answer
1k
views
The velocity Verlet method and variable time steps
Does the velocity Verlet handle variable time steps? I found controversial statements about it.
In the paper Skeel, R. D., "Variable Step Size Destabilizes the Stömer/Leapfrog/Verlet Method", BIT ...
- 143
3
votes
1
answer
2k
views
Calculation of Mean Square Displacement for Brownian dynamics system with Lennard Jones interactions in python3
I have a problem getting a sensible result for the Mean Square Displacement (MSD) for a simulation of $N$ particles under Brownian dynamics with Lennard-Jones interaction between them with or without ...
0
votes
0
answers
92
views
Super Computer cluster service for students [duplicate]
I would like to ask whether there is a free of charge supercomputer for students, in my current institution they don't have a supercomputer, I am running MD simulations, I would be grateful for ...
- 109
1
vote
1
answer
156
views
Langevin Thermostat and overdamped Langevin Equations
I'm having some difficulties understanding the Langevin thermostat (MD).
In my notes, there is written that the Langevin equation is
$$ m\dot{v} = F - m\gamma v + f_R, \tag{1} \label{1}$$
where $...
- 394
0
votes
1
answer
49
views
Can LINCS algorithm be used for colliding molecules?
Supposing that one molecule is static and one is dynamic,
can the dynamic one be solved with LINCS for its shape (angle, bond length) constraints and also keep collisions with static molecule off, ...
1
vote
2
answers
221
views
Custom exponents for Lennard Jones in LAMMPS
I am trying to run an MD simulation using this generalized version of Lennard Jones.
$$
U(r) = \left(\frac{r_0}{r}\right)^A -\frac{A}{B}\left( \frac{r_0}{r} \right)^B
$$
However, I do not know how ...
- 255
0
votes
2
answers
102
views
How to calculate $g_r$ when using the cell method algorithm
I have a self-written MD code in C++. I am simulating atoms using the Langevin dynamics model. I use the linked-cell method to speed up the simulation; as a result, atomic interactions are only ...
- 187
0
votes
1
answer
128
views
Calculate forces on atoms from potential energy of system and position of atoms
Background
I am using a neural network to calculate the potential energy of atoms in a configuration and then adding energy of all atoms to compare it with the true energy of the configuration(label) ...
- 131
2
votes
2
answers
559
views
Effective way to build the neighbor's list in MD
I'm trying to implement the following form of the cell/neighbor list method in my MD code. I have divided my simulation box into a fixed number of cells, and according to its positions, I have ...
- 187
2
votes
1
answer
1k
views
Interatomic distance-periodic boundary conditions-non cubic unit cell
I am trying to find interatomic distance considering periodic boundary conditions for hexagon cubic cells (graphite). I tried to follow the answers to these two questions here but am unable to get the ...
- 21
7
votes
2
answers
123
views
What good are hard-sphere event-driven molecular dynamics simulations in the face of chaos?
Simple hard-sphere dynamical systems can exhibit chaotic dynamics. Due to finite-precision arithmetic when implemented on a computer, the presence of chaos implies that for a given set of initial data,...
- 171
3
votes
1
answer
1k
views
MD Simulation: Reference for the Neighbor's List Method
With a rather basic knowledge in C++, I have written my own MD simulation code. Currently, I calculate forces in the most naive way: I go through all the atoms and account for their interactions. This ...
- 187
1
vote
2
answers
213
views
Discussing the energy spectrum of Langevin Dynamics simulation of many atoms
UPDATED
I've coded a multiparticle MD simulation in 3D. It is based on Langevin Dynamics, with random impulse and dissipation. I think the program works correctly now? I have attached the plots of ...
- 187
2
votes
2
answers
553
views
How to compute forces in multi-particle MD
Suppose we have a system of $N$ particles that interact via the Lennard-Jones potential $$V(r)=V_0\left[\left(\frac{r_0}{r}\right)^{12}-2\ \left(\frac{r_0}{r}\right)^{6}\right].$$ No other forces ...
- 187
1
vote
0
answers
22
views
Rigid and flexible ZIF-8 framework [closed]
What are the differences in parameters in the force field for MD simulation study of RIGID and FLEXIBLE ZIF-8 framework?Is the electrostatic interaction term present in the rigid one?
1
vote
1
answer
213
views
Oil/Water interface simulation using GROMACS
I would like to perform a decane/water(tip4p) interface simulation using GROMACS to analyse the fluctuation of the interface over a long period of time (50ps - 1ns).
Though I have a rough idea on ...
0
votes
2
answers
275
views
Is there any scalable machine learning tool like molecular dynamics software NAMD?
I work in Molecular dynamics of biomolecules. I am familiar with the software NAMD. We use NAMD in our CPU cluster with the aid of mpirun.
Without GPUs can I implement large scale deep learning?
...
- 81
1
vote
0
answers
43
views
How to implement boundary conditions in arbitrary geometries, in MD\Lagrangian FD simulations
Let's say I have a random arbitrary geometry, like this one https://ibb.co/dG8fUa (something I just made up). I make one of these in 3D in SolidWorks. I have molecules/particles coming in through the ...
- 149
0
votes
1
answer
101
views
Generate a set of random x,y,z numbers, with a minimum difference between them, between defined limits
So I am trying a molecular dynamics simulation, and trying to populate a 3D rectangular volume with random particles which are uniformly distributed throughout the volume. Each molecule has a fixed ...
- 149
3
votes
1
answer
161
views
Why are Hamiltonian dynamics used in MCMC?
In Hamiltonian Monte Carlo, Hamiltonian dynamics are used to generate new proposals from the current state. I understand why these dynamics are used as opposed to random walk behavior to generate ...
- 31
15
votes
1
answer
2k
views
Why is leapfrog integration symplectic and RK4 not, if the latter is more accurate?
In a system where energy theoretically should be conserved, the most accurate simulation would conserve energy (as well as giving accurate positions, velocities and etc). RK4 is more accurate than ...
- 265
1
vote
3
answers
2k
views
How to integrate numerically Nosé Hoover equation?
In NVT molecular dynamics, Nosé Hoover thermostat is a method defining an extended system. I can understand perfectly how the Nosé Hoover differential equations are derived (Frenkel&Smit's book). ...
- 21
0
votes
1
answer
59
views
Forcefield for studying a two phase system, consisting of noble gases and gaseous organic components
I am preparing to do a MD-simulation of a system which consists of two phases. The gaseous phase consists of $H_{2}$, $Ar$, $He$, $N_{2}$, $CH_{4}$ and $CO_{2}$. In this system, $p=0.355 ~bar$ and $T\...
- 151
1
vote
2
answers
1k
views
How to control temperature in NVE molecular dynamics
I want to simulate a molecular system at a certain temperature. It would be best if it could be implemented by NVT molecular dynamics. However, if one were to use a NVE simulation (like velocity ...
- 21
0
votes
2
answers
278
views
Fast Python implementation of short-range interacting particles under Metroplis algorithm
Can anyone write a Python implementation of a set of particles interacting in 2D according to a short-range particle-particle force and evolving in time under a Metropolis algorithm, which randomly ...
user21228
0
votes
1
answer
130
views
List of Tersoff potentials?
I would like to know some GNU libraries or lists of downloadable Tersoff potentials. I couldn't find that many of them, and I also would like to know if there is a reason as to why this is the case.
- 111
2
votes
0
answers
230
views
temperature set point maintenance with Langevin thermostat + shake algorithm
During the past ~4.5 years of my PhD, I have written an MD code, based largely on LAMMPS and NAMD, however reparallelized to suit my computational scenario better. One important thing to mention ...
- 21
3
votes
0
answers
337
views
Error in MD simulation with Lennard-Jones potential in C, resulting in a lack of conservation of energy
I'm working on a MD simulation using the Lennard-Jones potential and the Velocity Verlet algorithm. The code compiles and runs, but the energy for each time step is not conserved; it decreases ...
- 151
3
votes
0
answers
174
views
Beowulf Cluster For MD simulations
What would you recommend as an all-in-one cluster management software for a Beowulf cluster? We have about 7(most of them are 8 core machines and 1 is a 48 core machine) very good compute nodes in our ...
- 365
6
votes
1
answer
472
views
Is resampling more accurate than block average for statistical analysis of data?
I'm working in laboratories where molecular dynamics data are almost always analysed usign block average as stated in the famous Allen and Tildesley book.
We divide the datas in blocks of size $M$ on ...
- 323
0
votes
1
answer
39
views
Modeling Coupling: vibrational,rotational, etc
I am trying to figure out how to deal coupling. I am trying to model a polymer chain. I have to deal with the angle between adjacent angles, dihedral angles, and vibration. I am able to compute all ...
- 111
1
vote
3
answers
973
views
Choosing the correct Ewald sum parameters
I'm currently working with a system where the long range Coulomb interactions are described with the Ewald summation method.
I am using DL_POLY Classic. The Ewald summation method used is described ...
- 151
9
votes
1
answer
659
views
Ways to start ab initio MD from classical MD
I am running molecular dynamics simulations of water for testing purposes. The box is quite small, if you ask a guy running classical MD, and relatively large, if you ask a DFT guy: I have 58 water ...
- 145
2
votes
0
answers
101
views
Why not solve molecular dynamics trajectories exactly?
The basis of Molecular Dynamics simulation is the integration of newton's equations of motion:
$$\frac{F}{m} = \frac{d^2r}{dt^2}$$
I understand that there are various methods based on Taylor's ...
1
vote
1
answer
101
views
Creating dense random configuration in for molecular dynamics
I am creating a random configuration of particles for a molecular dynamics simulation, where I would like to guarantee a certain density.
The strategy is as follows:
choose randomly the positions of ...
- 11
1
vote
0
answers
197
views
Maxwellian distribution of velocities with Shake algorithm present
I am writing a code to perform hybrid monte carlo molecular dynamics. To do this, I need to have a code to initialize the velocities of all particles according to a maxwell distribution. The code is ...
- 291
2
votes
0
answers
235
views
Effect of Periodic Boundary condition in Pair Distribution function
When we simulate a bulk molecular system using molecular dynamics in (say) NVE ensemble (microcanonical), we apply some periodic boundary conditions (PBC) to avoid surface effects. What effect does ...
- 201