E-CAM velocities_resolve_EVC module¶
velocities_resolve_EVC is a module that resolves the excluded volume constraint with a velocity formulation (no potential
applied between two overlapped bonds).
velocities_resolve_EVC uses the module E-CAM minDist2segments_KKT module to find the minimal distance between two bonds.
To study the long term memory of the initial conformation of a highly entangled polymer we need to preserve the topology. It means that two bonds cannot cross. It is of great importance for the study of post-mitotic chromosome unfolding. To resolve the excluded volume constraints you could use a soft or hard potential between the two points associated to the minimal distance. Here, we propose to change the relative velocity between overlapped bonds to resolve the excluded volume constraint in one time-step of molecular dynamics.
- Polymer simulation.
- To resolve the excluded volume constraints.
- It is used in a scientific collaboration.
- Publications: not currently available.
We would use the present module to avoid topology violation in an entangled polymer system. The present module uses the E-CAM module E-CAM minDist2segments_KKT module.
This module is a part of a pilot project (E-CAM post-doc). We would use it to avoid topology violation in an entangled polymer system.
You can find a PDF file with a detailed derivation of the velocity-based method we use to resolve the excluded volume constraint in one time-step of molecular dynamics on the velocities_resolve_EVC GitLab repository.
I provide a simple Makefile you can find at the same location as the source code.
You need C++11 in order to use pseudo-random number generator.
Before the compilation you can clean the previous build with “make mrproper” command.
The purpose of the module is to resolve excluded volume constraints.
Therefore, we provide a simple example of a system of
N bonds with volume interactions.
We test every
n iterations the average overlap pairwise.