Universal Scaling Functions of the Gr{ü}neisen Ratio near Quantum Critical Points

The Grüneisen ratio, defined as Γg ≡ (1/T) (∂ T/∂ g)S, serves as a highly sensitive probe for detecting quantum critical points (QCPs) driven by an external feild g and for characterizing the magnetocaloric effect (MCE). Near a QCP, the Grüneisen ratio displays a universal divergence which is governed by a universality-class-dependent scaling function stemming from the scale invariance. In this work, we systematically investigate the universal scaling functions of Grüneisen ratio in both one-dimensional (1D) and two-dimensional (2D) quantum spin systems, including the transverse-field Ising model, the spin-1/2 Heisenberg model, the quantum q-state Potts model ($q=3,4$) and the J₁-J₂ columnar dimer model. Our approach employs the thermal tensor-network method for infinite-size 1D systems and the stochastic series expansion quantum Monte Carlo (SSE QMC) simulations for 2D systems, enabling precise calculations of the Grüneisen ratio near QCPs. Through data collapse analysis, we extract the corresponding scaling functions, which establish quantitative frameworks to interpret magnetocaloric experiments and guide the development of ultralow-temperature refrigeration.