Abstract: A series of detailed numerical simulations are used to investigate the properties ofmultipactor breakdown in circularwaveguidespropagating the TM01 mode. AMonte Carlo model is constructed to track the motion of the electrons, study the multipactor scenarios, and predict the multipactor thresholds. The theoretical and numerical analyses indicate that the product of the frequency and the gap (f . D) affects both the intensity of the ponderomotive force and its spatial distribution, which results from the nonuniformity of the radio frequency (RF) field and significantly influences the electrons' trajectoriesandmultipactor trends. The decrease in f . D results in a remarkable enhancement in the magnitude of the ponderomotive force, while the maximal intensity gradually moves toward the half radius R/2 area. Low values of f . D correspond to high ponderomotive potential, which sustains the short-range electrons and triggers the single-sidedmultipactor. In contrast, high values of f . D correspond to low ponderomotive potential, contributing to long-range electrons and exciting the double-sided multipactor. Fitting to the susceptibility diagram produces the border line and a modified f . D threshold of (f . D) th approximate to 338.4 GHz mm, which separates the susceptibility diagram into single-sided, double-sided, andmixed-sided zones. The initial electron energy influences their trajectories at high f . D and low RF power. This effect tends to dominate the multipactor behavior in the mixed-sided region.