All simulations, including geometry optimization and vibration analysis under quasi-hydrostatic pressure, were conducted using the CP2K code. The PBE function and the Goldk-Teter-Hutter (GTH) pseudopotential with the double-ζ-with-polarization (DZVP) basis set were employed. Additionally, the semiempirical dispersion correction D3 (BJ) was used to account for the intermolecular dispersion force. The simulations were performed under periodic boundary conditions, with a finest grid level cutoff set to 400 Ry for the Γ point. The Broyden-Flecher-Goldfarb-Shanno (BFGS) optimization method was used, with tighter convergence criteria applied for the subsequent calculations of infrared spectra. Finally, the wave-function analysis, such as non-covalent interaction and LBO, was conducted using the Multiwfn and VMD packages.