Hi, all. Just a minor point. In bond_fene.cpp, the “equilibrium_distance” function returns .97*sigma. Actually, the minimum of the FENE potential for the standard parameters (k = 30, R_0 = 1.5, epsilon = 1.0, sigma = 1.0) is at .960897. So, while .97 may be the actual equilib. distance for the standard *temperature* T = 1.0 epsilon/k_B, this result is not general, and besides depends on the parameters k and R_0.

Thanks,

Rob