|
From: | Peter Košovan |
Subject: | Re: Reactive polymer cross-linking |
Date: | Sun, 5 Jan 2020 23:54:52 +0100 |
Hi all,
I’m new to using espresso-md, so I’ve been pouring over the documentation and examples to better familiarize myself.
I’m starting with a somewhat simple example of poly-phenylene sulfide curing and cross-linking. So I have a few polymer chains with harmonic bonds, DPD style forces, and free oxygen ‘beads’.
I setup a reaction ensemble with the reaction to take a monomer location on two different molecules along with oxygen, and then bond them. I’m currently doing this now by introducing a new monomer type ‘linked’, which I then add a new bond type after seeing its been replaced in the reaction.
Long story short: the reaction ensemble is reacting over the entire simulation space, and picking any two monomer elements, reacting them and then forcing strange long distance bonds. Can I control the reaction volume with espressomd? Ideally I’d like a small radius around the reactants.
I also looked into the cool dynamic bonds thing, but it kept creating new bonds from ‘collisions’ of my monomer elements along a polymer chain.
Thanks!
Chris Buurma, Ph. D.
[Prev in Thread] | Current Thread | [Next in Thread] |