[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[ESPResSo-users] Bond Breaking in Polymer Networks - United Atom Forcefi
[ESPResSo-users] Bond Breaking in Polymer Networks - United Atom Forcefield
Tue, 13 Nov 2012 12:02:13 +0200
I'm interested in performing MD simulation of polymer networks under tension,
where some chains could break as elongation proceeds.
I use a United Atom (UA) forcefield.
The bonded interactions of my model include:
(i) a bond stretching potential
(ii) an angle potential (bond bending)
(iii) a dihedral (torsional) potential.
The nonbonded, pairwise interactions are based on a 6-12 Lennard-Jones (LJ) potential,
which acts between all possible pairs (inter- and intra-chain),
except from intra-chain pairs that are separated by 1, 2, or 3 sequential bonds,
(the so called (1-5) intra-chain LJ potential).
I would like to ask you if ESPResSo is able to break bonds 'on the fly', as simulation proceeds,
according to a distance criterion for the bond length.
In terms of potentials, a bond scission event means that
(a) For the UAs making this bond, the bond potential will not be evaluated (or will be turned off) after breakage
(b) The angle and dihedral potentials related with this bond will cease to exist (or be evaluated)
(c) A LJ pairwise interaction between the UAs of the broken bond will start to apply.
Thank you very much for your help,
- [ESPResSo-users] Bond Breaking in Polymer Networks - United Atom Forcefield,
Stefanos Anogiannakis <=