Molecular dynamics sampling scheme realizing multiple distributions

We present a molecular-dynamics sampling scheme in which any summation of multiple arbitrarily given distributions can be realized deterministically by a certain theoretical guide. Our scheme thus provides a static view of the composition of the distributions, rather than a dynamic view in which some parameter is switched at a certain time in the simulation process. The proposed method induces the orbit to jump automatically among different phase space regions, without the use of any artificial timing parameter. In addition, the proposed method does not require an assumption that the density of states is smooth. We applied it to multiple Tsallis distributions and established a suitable series of parameter values for which the sum of the distributions allows broad sampling. Numerical simulations applied to fundamental models with multi-Tsallis distributions showed efficient sampling, characterized by an energy trajectory that was totally different from that associated with each single distribution.