This code is designed to efficiently simulate the wormlike chain polymer model using various coarse-grainings where applicable.

For very stiff polymers, the usual wormlike chain is simulated.

For relatively more flexible polymers, the "stretchable, shearable" chain is used.

For VERY stretchable polymers, a purely Gaussian chain is used.

Simply typing make in the top level directory will build the simulator from source. The executable created (wlcsim.exe) will use the parameters in the file input/input and write its output to the data directory. Descriptions of the available parameters can be found at their definitions and where they are read in src/wlcsim/params.f03. Example input files are usually more useful, and can be found in the input directory. To force a rerun without having to manually delete all the old output files, you can also simply type make run at any time.

There are several ways to easily visualize simulation output. There are PyMol scripts in the vizualization directory, python -m wlcsim.plot_wlcsim from the repo's top level directory will launch a GUI designed to visualize BD simulations, and one can of course simply use the output in the data directory, which contains rank two arrays of shape num_beads*num_polymers-by-3, with one file per time point. By default, specifying multiple polymers just simulates them in parallel in the same reaction volume, no interactions are assumed.

To scan parameters, the Python script scan_wlcsim.py should be used. It takes care of saving the current git commit_hash, all inputs, etc. into a unique directory, and preventing race conditions even on shared filesystems, among other things.

This codebase is internal to the Spakowitz lab and is not guaranteed to be bug-free at any point. For battle-tested versions of our software, please see the links in the relevant papers.