Nonminimal Einstein–Maxwell–Vlasov-axion model