Abstract
We analyse numerically thermal fluctuations of the static overlap between equilibrium configurations in a glass-forming liquid approaching the glass transition. We find that the emergence of slow dynamics near the onset temperature correlates with the development of non-Gaussian probability distributions of overlap fluctuations, measured using both annealed and quenched definitions. Below a critical temperature, a thermodynamic field conjugate to the overlap induces a first-order phase transition, whose existence we numerically demonstrate in the annealed case. These results establish that the approach to the glass transition is accompanied by profound changes in the nature of thermodynamic fluctuations, deconstructing the view that glassy dynamics occurs with little structural evolution.