This paper is concerned with the following Maxwell–Dirac system
class="formula" id="fd000005">
where
class="mathmlsrc">class="formulatext stixSupport mathImg" data-mathURL="/science?_ob=MathURL&_method=retrieve&_eid=1-s2.0-S0362546X15002412&_mathId=si2.gif&_user=111111111&_pii=S0362546X15002412&_rdoc=1&_issn=0362546X&md5=ed83e0463a374f314f0f12672a7be0f6" title="Click to view the MathML source">M(x)class="mathContainer hidden">class="mathCode"> is a external potential and
class="mathmlsrc">class="formulatext stixSupport mathImg" data-mathURL="/science?_ob=MathURL&_method=retrieve&_eid=1-s2.0-S0362546X15002412&_mathId=si3.gif&_user=111111111&_pii=S0362546X15002412&_rdoc=1&_issn=0362546X&md5=5bfc4241dbf658fd71c8454739370851" title="Click to view the MathML source">F(x,u)class="mathContainer hidden">class="mathCode"> is an asymptotically quadratic nonlinearity modeling various types of interactions. In view of the effects of the nonlocal term, we use some special techniques to deal with the nonlocal term. Moreover, existence and multiplicity of stationary solutions are obtained for system without any periodicity assumption via variational methods.