This article is devoted to the study of a fully discrete A‐
ϕ
finite element method to solve nonlinear Maxwell's equations based on backward Euler discretization in time and nodal finite elements in space. The nonlinearity is owing to a field‐dependent conductivity with the power‐law form
|
E
|
α
−
1
,
0
<
α
<
1
. We design a nonlinear time‐discrete scheme for approximation in suitable function spaces. We show the well‐posedness of the problem, prove the convergence of the semidiscrete scheme based on the boundedness of the second derivative in the dual space and derive its error estimate. The Minty–Browder technique is introduced to obtain the convergence of the nonlinear term. Finally, we discuss the error estimate for the fully discretized problem and support the theoretical result by two numerical experiments. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 2083–2108, 2014