[lammps-users] fix thermal/conductivity

Dear all,
I’m calculating thermal conductivity of molecular crystals using fix thermal/conductivity. For flexible molecules, the atom exchange method works, but for molecules with constraints, the swap of velocities has to be done between center-of-mass of molecules. I wonder if anyone has amended this molecular swap method in LAMMPS. If not, I would probably give it a try, but I’m not an expert of C++, nor do I know LAMMPS code, any input from you is welcome.

Dong Wang

what kind of constraints are you topic about ?

like bond length constraints using SHAKE.

This would be non-trivial for large molecules, which span processors.
Why not do the calculation with flexible molecules, which are generally
better models anyway?