[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

From: |
Axel Kohlmeyer <akohlmey@...24...> |

Date: |
Wed, 20 Sep 2017 21:22:14 -0400 |

On Wed, Sep 20, 2017 at 4:31 PM, A. M.M <alaa.murar@...24...> wrote: > بسم الله الرحمن الرحيم > > Hello, > > I have some questions which are: > > 1. I am using fene/expand without LJ terms ( set epsillon and sigma equal to > zero ) as bond_style in my simulation, but there are errors which are too > long fene bond and bad fene bond. I know that fene potential has a > boundaries that the bonded atoms must not exceed them, but in my simulation > there are a two pair_style which are coul/debye and morse and these pairs > produce large repulsive ( same charge in coul ) forces when the distance > between the atoms becomes small and small and this leads to increase the > distance between atoms especially the bonded ones. therefore producing the > problem of the fene potential. > > How I can solve this problem ? this sounds like an issue with your force field or simulation parameters. either they are incorrect/unsuitable for your system. there are far too many details that matter, so it is next to impossible to help you from remote with this unless there is a big glaring error. your setup is complex and very specific, so you may need to approach this in stages and perhaps get some basic experience with a simpler system first. it is usually best to start with reproducing well documented examples from literature. talk to your adviser/supervisor. you need to consult with a local person with sufficient experience in MD. if that is not available, you need to join a group with such experience for a month or two. > 2. I am using fix langevin with fix nve to perform a Brownian dynamics and I > set Tstart equal to Tstop equal to 298.15 ( 25 C ) and the damp equal to > 39.4, but when the simulation run the temp. varies and I need it to be > fixed, so how can I fix it ? it looks like you need a better understanding statistical thermodynamics. on the microscopic scale, the _instantaneous_ temperature *always* fluctuates. that the most extreme example: a harmonic oscillator. at its turning point, all energy is potential energy, kinetic energy is 0. at equilibrium distance you'll have the maximum kinetic energy 1/2*m*v^2. on average, it will be half of that. in a typical simulation system at the microscopic level you have the equivalent of multiple of those oscillators and they are coupled, too, so they can transfer energy from one to the other. overall this leads to a situation similar to having multiple harmonic oscillators at different force constants, with varying kinetic energy and different phase. => there will be oscillations and those will be the larger, the smaller your system is, or get smaller, the larger your system gets. a "constant" temperature, you will only have in the limit of an infinite number of particles. however, for a system in equilibrium, you need not only average over particles, but you can also average over time. this should give you enough basic information to look up the details in a text book on statistical thermodynamics and learn why asking for an exactly constant temperature in an MD simulation makes no sense. > 3. for the damp parameter in fix langevin I do the follwing to include the > viscosity of the water as an implicit solvent; please note, that this mailing list is not a classroom and we are not your teaching assistants or your supervisors tasked with checking and grading your work. axel. > the coeff of the drag focre is m / damp and in BD the coeff is equal to > 3*pi*eta*diameter ( in other program i found 6 instead of 3 what the > difference ? ) > > I set m / damp = 3*pi*eta*d, therefore damp = m/3*pi*eta*d > > I put m ( mass of the beads in my simulation = 0.113 kg/mole > eta = 0.00089 kg/m.s for water > d ( for the bead ) = 5.68e-10 m > > the unit of the answer is s/mole > > I divide the answer by Avogadro's number and get 39.4 fs ( unit of time for > real ) > > is this correct for BD in water ? > > Thanks > > > ------------------------------------------------------------------------------ > Check out the vibrant tech community on one of the world's most > engaging tech sites, Slashdot.org! http://sdm.link/slashdot > _______________________________________________ > lammps-users mailing list > lammps-users@lists.sourceforge.net > https://lists.sourceforge.net/lists/listinfo/lammps-users > -- Dr. Axel Kohlmeyer akohlmey@...24... http://goo.gl/1wk0 College of Science & Technology, Temple University, Philadelphia PA, USA International Centre for Theoretical Physics, Trieste. Italy.

**References**:**[lammps-users] Some questions***From:*"A. M.M" <alaa.murar@...24...>

- Prev by Date:
**Re: [lammps-users] Using pair style hybrid with GPU** - Next by Date:
**Re: [lammps-users] ERROR: Invalid fix style** - Previous by thread:
**[lammps-users] Some questions** - Next by thread:
**[lammps-users] accepted force field deviation** - Index(es):