While this is not a direct answer, it might hopefully offer some inspiration from a neighbouring field.
It turns out that there is something called DNA melting. The DNA melting temperature is defined as the temperature at which half of the DNA strands are in the random coil or single-stranded state. While not at all a true melting temperature in the conventional sense, it bears a couple of key resemblances. First, like the melting point of solids (metals, in your case), this is a thermodynamical compromise between entropy and enthalpy. Second, like the melting point of solids, these terms can be estimated computationally by considering interatomic interactions, and this is non-trivial. So maybe there is something to be learned of how people from that field have been dealing with this problem that, for them, is of great practical importance.
Here is where (hopefully) the inspiration comes in: Leber et al, "A fractional programming approach to efficient DNA melting temperature calculation", Bioinformatics 2005, 21, 2375–2382.
Results: As the melting temperature can be expressed as a fraction in terms of enthalpy and entropy differences of the corresponding annealing reaction, we propose to use a fractional programming algorithm, the Dinkelbach algorithm, to solve the problem. To calculate the required differences of enthalpy and entropy, the Nearest Neighbor model is applied. Using this model, the substeps of the Dinkelbach algorithm in our problem setting turn out to be calculations of alignments which optimize an additive score function. Thus, the usual dynamic programming techniques can be applied. The result is an efficient algorithm to determine melting temperatures of two DNA strands, suitable for large-scale applications such as primer or probe design.
We can also see at least one instance of people employing LAMMPS to deal with DNA melting: Svaneborg, "LAMMPS framework for dynamic bonding and an application modeling DNA", Computer Physics Communications 2012, 183, 1793-1802 So maybe there is a path forward? If you have a way of estimating enthalpy and entropy for small chunks of your fcc metal as a function of pressure, perhaps this approach can lead you somewhere.
 For more on the Dinkelbach's algorithm: You et al, "Dinkelbach's algorithm as an efficient method to solve a class of MINLP models for large-scale cyclic scheduling problems" Computers & Chemical Engineering 2009, 33, 1879-1889.