Phase Frustration Induced Intrinsic Bose Glass in the Kitaev-Bose-Hubbard Model

We report an intrinsic "Bubble Phase" in the two-dimensional Kitaev-Bose-Hubbard model, driven purely by phase frustration between complex hopping and anisotropic pairing. By combining Inhomogeneous Gutzwiller Mean-Field Theory with a Bogoliubov-de Gennes stability analysis augmented by a novel Energy Penalty Method, we demonstrate that this phase spontaneously fragments into coherent islands, exhibiting the hallmark Bose glass signature of finite compressibility without global superfluidity. Notably, we propose a unified framework linking disorder-driven localization to deterministic phase frustration, identifying the Bubble Phase as a pristine, disorder-free archetype of the Bose glass. Our results provide a theoretical blueprint for realizing glassy dynamics in clean quantum simulators.