Phase stability in multicomponent solid solutions depends on configurational entropy beyond the ideal mixing limit, but capturing it together with vibrational entropy within the same atomistic framework remains challenging. Here, we extend non-equilibrium thermodynamic integration to composition-dependent transformations through an alchemical interpolation of the interactions, combined with Monte Carlo identity exchange moves and molecular dynamics that sample the vibrational and non-ideal configurational entropy along the integration path. We apply the framework to the Au-Cu binary alloy using Atomic Cluster Expansion potentials trained on density functional theory data using the LDA, PBE, and r2SCAN functionals, and construct composition-temperature phase diagrams directly from atomistic free energies. We find that explicit configurational sampling lowers the AuCu order-disorder transition temperature predicted by the ACE potential trained on LDA data from approximately 810 K to 710 K, closer to the experimental value of 683 K, and substantially widens the stability range of the solid solution. At the same time, the much larger sensitivity to the exchange-correlation functional shows that this level of agreement should not be interpreted as general predictive accuracy. Non-ideal configurational entropy must therefore be sampled explicitly, alongside a careful choice of functional, for a reliable atomistic description of binary phase diagrams.
Sarath Menon, M. Poul, T. Hickel et al.· 0 citations
Aromatic organic solutes in water exhibit a delicate balance between hydrophobic solvation and directional O-H$\cdots \pi$ hydrogen bonds, yet widely used force fields and state-of-the-art density functional approaches struggle to provide a consistent picture of these pivotal interactions. We introduce a data-efficient upfitting strategy to train a machine learning interatomic potential (MLIP) based on the graph atomic cluster expansion for aqueous aromatic molecules with CCSD(T) accuracy for condensed phase simulations, using only finite molecular clusters. We apply our method to aqueous toluene (C$_6$H$_5$CH$_3$). The resulting CCSD(T)-quality MLIP reproduces coupled cluster energies and forces in bulk and reveals that commonly employed methods do not capture the crucial balance between hydrophilic and hydrophobic solvation, distorting the interactions of aromatic molecules with their environment. Representative biomolecular force fields substantially understructure the hydrophobic solvation shell and misorient interfacial water, while overestimating $\pi$-contacts, yielding an inconsistent solvation balance. Even hybrid DFT and MP2 overestimate barriers to breaking of water-$\pi$ hydrogen bonds. Our workflow provides a practical, general route to CCSD(T)-quality condensed-phase simulations of aqueous solutions, and thus constructed interaction potentials now open the door to consistent, highly accurate benchmark studies of $\pi$-contacts and hydrophobic effects in biomolecular contexts such as solvation of proteins and DNA in aqueous environments.
N. Stolte, H. Forbert, Yu. Lysogorskiy et al.· 0 citations
