|
From: | address@hidden |
Subject: | Re: [ESPResSo-users] FENE+LJ bond |
Date: | Tue, 21 May 2019 21:11:19 +0800 |
Actually, they did not specify in the paper for which pairs of particles they used the usual LJ potential (Eq.1). They only say that the fene + LJ (WCA) with cut=2^1/6 sigma) is used for polymer-polymer interactions. The text suggests that they used Eq.1 for all particle pairs including the polymer but that would be a very wierd setting.If you need to model exactly the same system, then I think that you should ask the authors for clarification. Otherwise, I would use WCA potential for polymer-polymer and LJ with cut=2.5sigma fro all other pairs.Greetings,peterOn Tue, May 21, 2019 at 2:54 PM Dusan Racko <address@hidden> wrote:Hi, isn't that a cutoff for the attractive potential to model magnetic particles only? Sorry, I don't know the papers, just suggestingDear Peter,Actually, in the paper, the cutoff distance of LJ potential between nonbonded beads is 2.5[sigma] while the cutoff distance of LJ potential between bonded beads is 2^(1/6)[sigma]. This can't be done by your way.Thanks anyway!Regards!Ricky
From: Peter KošovanDate: 2019-05-21 20:44To: address@hiddenCC: address@hiddenSubject: Re: [ESPResSo-users] FENE+LJ bondHi Ricky,Just set up a FENE potential between the given pair of particles, and an LJ potential between the given particle types (probably the same type) and you get exactly what you want. This is the standard setup which is just described in a weird way in the paper that you cited.The fene potential acts just between the two particles with a given id. The LJ potential acts between any pair of particles. If your pair of particles interacts both with fene and LJ, then you get the fene+LJ interaction.See also the script minimal-polymer.py in the samples directory.peterDear Sebastian,Thanks for your reply at first!I want to model the hard-soft block copolymers with magnetic particles and the copolymer model I use is from this paperIn this paper, the bond potential is FENE+LJ. Beyond that, I can also find paper using FENE+LJ potentials in https://aip.scitation.org/doi/10.1063/1.458541. What's more, in LAMMPS, it also uses the FENE+LJ type fene bond.I will use the FENE+LJ bond with bonded beads and the LJ potential with nonbonded beads.Regards!Ricky
--www.natur.cuni.cz/chemistry/fyzchem/Přírodovědecká fakulta Univerzity Karlovy v PrazeKatedra fyzikální a makromolekulární chemieFaculty of Science, Charles University in Prague, Czech RepublicDr. Peter Košovan
Department of Physical and Macromolecular Chemistry
Tel. +420221951029Fax +420224919752--_____________________Polymer Institute of the Slovak Academy of SciencesDubravska cesta 3845 41 Bratislava, Slovak Republic [‡]tel: +421 2 3229 4321
--www.natur.cuni.cz/chemistry/fyzchem/Přírodovědecká fakulta Univerzity Karlovy v PrazeKatedra fyzikální a makromolekulární chemieFaculty of Science, Charles University in Prague, Czech RepublicDr. Peter Košovan
Department of Physical and Macromolecular Chemistry
Tel. +420221951029Fax +420224919752
[Prev in Thread] | Current Thread | [Next in Thread] |