Questions tagged [intermolecular-forces]

Refers to forces that mediate interaction between molecules, including forces of attraction or repulsion which act between molecules and other types of neighboring particles, e.g. atoms or ions.

Filter by
Sorted by
Tagged with
3 votes
1 answer
45 views

Energy balance in a molecular dynamics simulation with ghost atoms and periodic boundary conditions?

I am performing a simulation in LAMMPS with granules (DEM) and I had a question regarding the conservation of energy. If an interaction happens where one granule is in the simulation box and the other ...
6 votes
0 answers
41 views

Is it possible to evaluate rotational and vibrational partition functions for diffusion in the Eyring Theory of Absolute Rates?

Let us consider the usual expression for Diffusion processes and diffusion coefficients in the Eyring's theory: $$D=\lambda^2k = \lambda^2 \frac{k_BT}{h}\frac{\Omega_\ddagger}{\Omega}e^{-\frac{\...
2 votes
1 answer
46 views

Rotational invariance of anisotropic polarization models

After writing out this question, I actually figured out the answer, but I am surprised I haven't seen this detail mentioned in almost any paper on the topic, so I will answer my own question here in ...
  • 3,249
8 votes
1 answer
136 views

Developing a machine-learned interatomic potential for molecular dynamics

There are few questions in this SE about the advantages/disadvantages of machine-learned interatomic potentials in molecular dynamics, but there is not much discussion on how to make them. I thought ...
  • 271
7 votes
2 answers
144 views

Which equation makes it possible to calculate the forces of attraction between a surface and molecules?

Which equation makes it possible to calculate the forces of attraction between a surface and molecules? Is there any free software that implements such a calculation? [Clarification 1] Between a ...
  • 601
4 votes
1 answer
27 views

What's a quick way to get a molecular bond length using a Newtonian model?

Can anyone spare me a clue on this? I'm a self-taught student of Molecular Geometry trying to learn more about methods in modeling. After some extensive reading and study, it felt like it should be ...
  • 53
4 votes
0 answers
70 views

Periodic boundary conditions or spherical potential? [closed]

I was wondering whether you know if there is any sort of literature discussing whether periodic boundary conditions or spherical potentials could be better or worse for QM/MM or MM simulations, their ...
7 votes
0 answers
92 views

Sigma Profile in COSMO - RS model [closed]

As far as I understand, sigma profile in COSMO - RS is a probability distribution of screening charge density for certain molecular surface segment. In COSMO - RS, molecular surface area is segmented ...
4 votes
0 answers
43 views

Gromacs generate forces using .pdb topol.top and .mdp file [closed]

Given the start configuration (.pdb), topology file (topol.top) and run file(water.mdp). I want to generate the forces on atoms. I tried, using: ...
  • 1,083
7 votes
0 answers
154 views

What does the shape of the radial distribution function(RDF) indicate? [closed]

From my understanding, generally in the RDF plot for liquids, the value of $g(r)$ goes to a maximum near the first solvation shell and then goes to a minimum and then slowly tapers off to the value of ...
9 votes
1 answer
84 views

How might inclusion of a Grimme dispersion correction influence a simulation of a dense liquid that is not water?

I adapted an example first-principles water simulation to do simulations of liquid SiO2, but did not know what I was doing and did not remove a Grimme dispersion correction term. It would be great to ...
  • 709
14 votes
2 answers
541 views

How do you incorporate hydrogen bonding into molecular simulations?

How do computer simulations of liquid hydrogen-bonding systems incorporate hydrogen bonding? Since I think hydrogen bonding has mixed covalent and electrostatic character, wouldn't you need something ...
  • 1,453
6 votes
1 answer
717 views

What is the polarization catastrophe?

I heard about the polarization catastrophe in a talk about polarizable force fields. The speaker talked about how the dipole induction is damped at short range to avoid the “polarization catastrophe&...
  • 7,869
6 votes
1 answer
281 views

Should I use a different cut-off scheme for non-bonded interactions when simulating a system with MD in vacuum?

I am currently doing a simulation of the behaviour of various sizes of water droplets (~100-1000 water molecules) in vacuum. I am using NAMD with the spherical boundary conditions to keep the droplet ...
  • 6,216
7 votes
0 answers
74 views

Molecular dynamics to reproduce dispersion interactions [closed]

I'm trying to simulate the behavior of chromophore molecules incorporated into a polymer under an external electric field. I've tried to use classical molecular dynamics (GROMACS) but it's unable to ...
  • 2,333
6 votes
0 answers
58 views

References for Electrostatics Calculations Used in Molecular Dynamics

I'm looking for classic references that describe the usual electrostatic contributions which have to be calculated in molecular dynamics simulations. I find it very hard to apprehend what different ...
  • 3,249
10 votes
1 answer
218 views

Are there good tutorials that explain how to use experimental data to customise the force field in AMDock Vina?

I am not a chemist, but I know how to use AMDock to predict the binding affinity of a small molecule to a protein. There is a particular target for which AMDock Vina, with its default force fields, ...
8 votes
2 answers
125 views

Can the NCI (noncovalent interaction) index be generated starting with psi4?

Can NCI[1] index plots be generated starting from psi4? I'd be happy to just generate the NCI index values and mesh the isosurface myself. With a long-term goal to visualise various medicinal-...
  • 876
8 votes
1 answer
141 views

Should I add vdW interaction to estimate chemical potential of metallic elements?

I am considering layered oxides $\ce{LiCoO2}$. Previous papers suggested that vdW interaction should be added to get lattice parameters closer to experimental data. I would like to dope other elements ...
  • 1,935
17 votes
1 answer
268 views

Empirically determining thermostat damping factor

Is there a way to empirically determine an appropriate thermostat damping factor given a timestep size and a numerical integration method? For example, I would like one for the surface and the ...
  • 445
11 votes
1 answer
345 views

How does one compute the boiling point of a liquid made of a particular molecule?

This question is in relation to Anomalous boiling point of "iso-" alkanes on Chemistry SE. I want to compute the boiling point of the different isomers of an n-carbon alkane to check what ...
12 votes
1 answer
176 views

Rule-of-thumb for Morse potential cutoff in molecular dynamics?

Is there a rule-of-thumb for selecting a cutoff for the Morse potential in molecular dynamics simulation? For instance, is it typically some multiple of the equilibrium distance?
  • 985
23 votes
3 answers
834 views

Simulating breaking bonds in molecular dynamics

How does one introduce the possibility of breaking a bond in a molecule during MD simulation? I only found the cases when we just introduce the harmonic potential $U_{ij} = \frac{1}{2}k(r_i - r_j)^2$ ...
15 votes
1 answer
207 views

DFT modeling intermolecular interactions

Density Functional Theory based codes evaluate quite successfully intramolecular forces responsible of the atomic bonding. Within this technique one is able to obtain molecular relaxed structures that ...
  • 1,560
22 votes
1 answer
233 views

Benchmarking data for Ewald summation algorithms

Are there good benchmarks available for testing homemade Ewald summation code? Optimally it would be very nice if there was several types of systems. NIST has a very nice benchmark for energies of ...
  • 4,266
27 votes
3 answers
2k views

What is the appropriate way of determining a value for the Hubbard-like potential U for LDA+U / GGA+U calculations?

It is known that neither LDA or GGA correctly account for the strong on-site Coulomb interaction of localized electrons[1]. A cost-effective way of correcting this is by using a Hubbard-like term $U$. ...
24 votes
2 answers
342 views

What are the different ways of calculating dispersion constants?

There are many different dispersion corrections out there. The most famous is D3 [1] (and the new D4 [2]), but there's probably other approaches too. The dispersion energy can be written as: $$ E_{\...
36 votes
3 answers
899 views

In molecular mechanics, how are van der Waals forces modelled?

In terms of energy, how are van der Waals forces modelled (are there formulas/laws that govern these)?