Hi there,
I am trying to generate a polymer by using espresso polymer module(version 4.0.2).
Here is my script:
import sys
import numpy as np
import espressomd
from espressomd.interactions import FeneBond
from espressomd import polymer
num_poly = 1
num_mono = 5
fene = FeneBond(k=30, d_r_max=2)
polymer.create_polymer(start_pos=[1, 1, 1], N_P=num_poly,bond_length=0.9, bond=fene,MPC=num_mono)
But it always raises a excetion that the bonded interaction passed as 'bond' keyword argument has not yet been added to the list of active bonds in Espresso.
Does anyone know how to fix this? Thank you very much.
Best,
Shawn