The replica-permutation method (RPM) and the replica sub-permutation method (RSPM) have been proposed as improved alternatives to the replica-exchange method (REM). Here, we implemented the RPM and RSPM in the NVT and NPT ensembles into an open-source classical molecular dynamics package, LAMMPS. We applied the RPM and RSPM to a polyethylene chain in vacuum and an alanine dipeptide in explicit water to test the implemented codes. We demonstrated that the RPM and RSPM by our codes achieved higher transition ratios of temperatures and faster convergence of physical quantities than the REM. We also validated that the RPM and RSPM generate statistical ensembles correctly. Program summaryProgram Title: temper/permutation, temper/subpermutation, temper/permutation/npt, and temper/subpermutation/nptCPC Library link to program files:https://doi.org/10.17632/t8pbjy8zb7.1Developer's repository link:https://github.com/MasatakaYm/lammps-rpm-rspmLicensing provisions: GNU General Public License, version 3Programming language: C++/MPINature of problem: The replica-permutation and replica sub-permutation methods [1–3] require implementations to construct permutation or sub-permutation function, to divide replicas into subsets when many replicas are used, and to attempt permutation by the Suwa–Todo algorithm [4].Solution method: We have implemented the replica-permutation and replica sub-permutation methods as an extension of existing temper command of LAMMPS.