In this paper, we are concerned with the existence of multiple solutions to the critical magnetic Schrödinger equation (0.1) \((-{\rm i}\nabla-A(x))^2u+\lambda V(x)u=\mu \vert u\vert ^{p-2}u+\Big(\int_{\mathbb{R^N}}\frac{\vert u(y)\vert ^{2^*_\alpha}}{\vert x-y\vert ^\alpha}{\rm d}y\Big)\vert u\vert ^{2^*_\alpha-2}u\quad {\rm in}\ \mathbb{R}^N,\)
where \(N\geq4, 2\leq p<2^*, 2^*_\alpha=\frac{2N-\alpha}{N-2}\) with 0 < α < 4, λ > 0, μ ∈ ℝ, A(x) = (A1 (x), A2 (x), ⋯, AN (x)) is a real local Hölder continuous vector function, i is the imaginary unit, and V(x) is a real valued potential function on ℝN. Supposing that Ω = int V−1 (0) ⊂ ℝN is bounded, we show that problem (0.1) possesses at least catΩ(Ω) nontrivial solutions if λ is large.