Gravitational form factors of the pion in the self-consistent light-front quark model

We present a self-consistent light-front quark model (LFQM) analysis of the pion's gravitational form factors (GFFs), incorporating the Bakamjian-Thomas (BT) construction consistently throughout the framework. By uniformly applying the BT formalism to both hadronic matrix elements and their associated Lorentz structures, we achieve a current-component-independent extraction of the pion GFFs $A_π(t)$ and $D_π(t)$, thereby eliminating the light-front zero-mode ambiguities that typically hinder conventional LFQM approaches. By tuning the model parameters, we identify an optimal set that successfully reproduces the decay constant and electromagnetic form factor of the pion, while yielding a $D$-term value $D_π(0) \approx -1$, consistent with predictions from chiral perturbation theory. The $D$-term emerges as a sensitive probe of the pion's internal dynamics, governing its mechanical radius -- the largest among the charge, mass, and mechanical radii. We further examine the pion's spatial structure via its associated two-dimensional light-front densities, including the momentum density, transverse pressure, and shear stress, all of which satisfy the required normalization and von Laue stability conditions. Our results reveal a detailed mechanical landscape: a centrally peaked momentum density that decreases monotonically; a repulsive pressure near the center (up to $x_\perp = 0.33$~fm) that transitions to attraction in the outer region; and a shear stress profile peaking at an intermediate distance ($x_\perp \approx 0.2$~fm).