Xi-Min Liang, Yuan-Zhu Wang, Tao Zhu, Wen Zhao, Xin Zhang

Gravitational wave (GW) observations provide an unprecedented laboratory for testing general relativity (GR) in the strong-field, highly dynamic, and relativistic regimes. Within the parameterized post-Newtonian (PN) formalisms, waveform generation tests have conventionally been limited to constraining inspiral coefficients up to the 3.5PN order. Leveraging the recent theoretical breakthrough that extended the analytical compact binary phasing to the 4.5PN order, we present the first observational constraints on these higher-order effects. Our analysis utilizes two exceptional events detected by the LIGO-Virgo-KAGRA (LVK) network: GW250114\_082203, which boasts the highest signal-to-noise ratio (SNR) recorded to date, and GW230627\_015337, which features a uniquely prolonged inspiral phase and the highest inspiral phase SNR to date. By performing Bayesian inference on the dimensionless deviation parameters ($δφ_i$) associated with the 4PN and 4.5PN coefficients, we find that our results are fully consistent with the predictions of GR. While the current 90\% credible intervals for the four deviation parameters are of order $\mathcal{O}(1) \text{-} \mathcal{O}(10)$, the general relativistic null values ($δ\hatφ_a= 0$) are entirely encapsulated within the bounds. This investigation establishes the first empirical baseline for 4PN and 4.5PN inspiral tests of GR, paving the way for high-precision null tests of GR with current and next-generation GW detectors.