Dear LAMMPS developers and users，I want to do a simulation similar to the active Brownian motion, which is achieved through overdamped Langevin. In the case of a small mass of 1, the frictional resistance is set very high to achieve this,a frictional force Ft = −ζtv and a torque Tr =−ζrω, ζt=6060, ζr=2020.
in LAMMPS
fix 1 mol_all langevin 1.0 1.0 0.000165 699483 (ζt=m/damp,damp=0.000165)
fix 2 mol_all viscous/sphere 2020
fix 3 core propel/self dipole 364
But when I run the program，ERROR on proc 0: Non-numeric atom coords - simulation unstable (src/OPENMP/domain_omp.cpp:58)
Can fix langevin operate under high frictional resistance？
I tried to reduce the frictional resistance
fix 1 mol_all langevin 1.0 1.0 0.0165 699483 (ζt=m/damp,damp=0.000165)
fix 2 mol_all viscous/sphere 20.2
It can run
I have already pointed out that debugging such cases is your job. Speculation won’t get you to a viable solution, you first need to understand why you are getting the error and identify its source at the core. Then you can start thinking and discussing about pathways to get more meaningful physics.
It is next to impossible to give advice from remote with limited context information and without being familiar with (let alone being an expert in) the kind of research you are doing.