Spectral profiles of two selected transitions for each O-bearing COM. The black curve shows the best-fit model with multiple kinematic components, and the colored dashed lines represent the contribution from each component. The blue, red, and gray bars at the bottom indicate the systemic velocities of VLA4A, VLA4B, and the circumbinary disk, respectively. The insets at the bottom of each panel show the residual from the best-fit model.