We have estimated the sensitivity to the anomalous couplings of the $\gamma\gamma\gamma Z$ vertex in the $\gamma\gamma\rightarrow\gamma Z$ scattering of the Compton backscattered photons at the CLIC. Both polarized and unpolarized collisions at the $e^+e^-$ energies 1500 GeV and 3000 GeV are addressed, and anomalous contributions to helicity amplitudes are derived. The differential and total cross sections are calculated. We have obtained 95\% C.L. exclusion limits on the anomalous quartic gauge couplings (QGCs). They are compared with corresponding bounds derived for the $\gamma\gamma\gamma Z$ couplings via $\gamma Z$ production at the LHC. The constraints on the anomalous QGCs are one to two orders of magnitude more stringent that at the HL-LHC. The partial-wave unitarity constraints on the anomalous couplings are examined. It is shown that the unitarity is not violated in the region of the anomalous QGCs studied in the paper.

Comment: 29 pages, 7 figures, 3 tables; minor corrections to the text, 4 references added