Measure-preserving integrators for molecular dynamics in the isothermal-isobaric ensemble derived from the Liouville operator

Density distribution of n-decane at 333 K and 1 atm using NPT(iso).