Prometheus: Unsupervised Discovery of Phase Transitions and Order Parameters in the Two-Dimensional Ising Model Using Variational Autoencoders

We present an unsupervised machine learning framework for discovering phase transitions and order parameters in the two-dimensional Ising model without prior knowledge of the underlying physics. Using convolutional beta-variational autoencoders (β-VAE), we learn compressed representations that naturally separate ordered and disordered phases. Our physics-informed training procedures incorporate data augmentation respecting Ising symmetries, cosine annealing optimization, and progressive temperature curriculum learning. Experimental validation on Monte Carlo simulated configurations demonstrates automatic order parameter discovery with correlation ρ = 0.85 ± 0.04 to theoretical magnetization (95% CI [0.82, 0.88]) and critical temperature detection within 0.27% of the Onsager solution (p < 0.001). The method achieves 89% improvement over principal component analysis and 124% improvement over t-distributed stochastic neighbor embedding. Latent representations exhibit clear phase separation with physics consistency score S = 0.88 ± 0.03. Ablation studies reveal optimal performance at β = 1.0, balancing reconstruction fidelity with disentanglement. This framework enables automated identification of emergent phenomena in complex many-body systems, advancing machine learning applications in condensed matter physics.