Efficient method to include nuclear quantum effects in the determination of phase boundaries

We developed a methodology to assess nuclear quantum effects in phase boundaries calculations that is based on the dynamical integration of Clausius-Clapeyron equation using path integral simulations. The technique employs non-equilibrium simulations that are very efficient. The approach was applied to the calculation of the melting line of Ne in an interval of pressures ranging from 1 to 3366 bar. Our results show a very good agreement with both experimental findings and results from previous calculations. The methodology can be applied to solid and liquid phases, without limitations regarding anharmonicities. The method allows the computation of coexistence lines for wide intervals of pressure and temperature using, in principle, a single simulation.