In the previous emails, I see somebody posted that the different energies existed between LAMPS and Material Studio for the polymeric materials with the same force field. You said that the difference was caused by the setting of special_bonds. Would you please tell me the default setting of special_bonds for PCFF in materials studio, using which we can get the same result from LAMMPS and material studio?
Would you please tell me the default setting of special_bonds for PCFF in materials studio, using which we can get the >same result from LAMMPS and material studio?
I have a question on the command "deform". Does this command only simply change the cell box length and let the atoms insides adjust themselves, or rescale the positions of all the atoms directly every step?
I got a trouble when using ave/spatial to calculate the density of the bulk materials (sample :PE). I build the model in Material Studio with a initial density 1 g/cm^3, and then transferred it into LAMMPS, and run a short NVT. I got the density file on X Y Z diractions, and I found the value is about 0.5-0.8 g/cm^3, less than the initial 1 g/cm^3. I send you the attached file, would you please check this for me?
I can't debug people's (complicated) input scripts. It's a simple
calculation. Take a snapshot. Find which atoms are in a layer.
Add up their mass. Divide by volume of the layer. You can do
it by hand for a small system. If you think LAMMPS isn't
giving the same answer as you compute, then show me why.