The Theoretical Approaches to Complexity of Phase Transitions in Condensed Matter: From Van Der Waals Theory to Spectral Mean-Field Method

Abstract

The theoretical description of phase transitions in condensed matter represents one of the most profound and enduring challenges in modern physics, spanning more than a century of intellectual development. From the phenomenological models of van der Waals in 1873 and Landau in 1937, through the statistical-mechanical rigor of integral equation theories of Ornstein and Zernike in 1914, Percus and Yevick in 1958, to Wilson’s renormalization group in 1971 and contemporary spectral methods such as the Self-Consistent Spectral Potential Mean-Field (SC-SPMF) method of Villani, the theoretical landscape has evolved through successive refinements that have progressively deepened our understanding of critical phenomena, metastability, and universality. Despite these remarkable advances, the prediction of phase transition temperatures, the characterization of critical behavior, and the understanding of metastability and hysteresis remain active areas of research, particularly for first-order transitions such as melting and crystallization. This review develops a central thesis: the difficulty of describing phase transitions in condensed matter arises from the intrinsically irreducible complexity of N-body systems. The behavior of a system of N interacting particles cannot, in general, be reduced to the behavior of N independent two-body systems. This is the essence of the N-body problem: the interactions between particles create collective effects that cannot be understood by considering pairwise interactions in isolation. The irreducible nature of this problem manifests itself in the hierarchy of correlation functions, where the n-body distribution function depends on the (n + 1)-body distribution function, leading to an infinite hierarchy that cannot be closed without approximations. The SC-SPMF method resolves the complexity of the problem through a spectral reduction that transforms the irreducible N-body problem into a self-consistent two-body problem in a mean field. The spectral criterion for phase transitions introduced by the SC-SPMF method represents a fundamental departure from traditional approaches: the transition is marked by a change in the curvature of the first eigenvalue μ1 (T). This criterion does not rely on arbitrary thresholds but on a qualitative change in the spectral response of the system, reflecting a fundamental change in the nature of correlations. The review is organized into three main parts. The first establishes the historical and conceptual foundations, tracing the evolution of theoretical approaches from van der Waals to the renormalization group, examining the successes and limitations of each approach, and identifying the common thread of approximate closures and functionals that has motivated the development of alternative methods. The second addresses the fundamental problem of complexity in N-body systems, exploring the BBGKY hierarchy, the irreducibility of many-body correlations, and the failure of pairwise reduction; it then presents the spectral mean-field approach as a rigorous solution to this problem, demonstrating how the eigenvalue spectrum of the Fisher operator encodes the collective behavior of the system. The third provides a comprehensive account of the SC-SPMF method. The SC-SPMF method offers several distinct advantages over existing approaches: it requires no approximate closures (unlike integral equation theories), no excess free energy functionals (unlike density functional theory), and no extensive computational resources (unlike molecular dynamics or Monte Carlo simulations). The coordination shell structure emerges self-consistently from the spectral decomposition, and the method naturally captures hysteresis and metastability because the self-consistent cycle provides stable solutions as long as the Fisher equation yields negative eigenvalues. The predictive power of the method is demonstrated by its ability to reproduce experimental structural parameters and to reveal the path-dependent universality that has not been observed in previous theoretical studies.