We develop a method of variational optimization of the infinite projected
entangled pair states on the honeycomb lattice. The method is based on the
automatic differentiation of the honeycomb lattice corner transfer matrix
renormalization group. We apply the approach to the antiferromagnetic
Heisenberg spin-1/2 model on the honeycomb lattice. The developed formalism
gives quantitatively accurate results for the main physical observables and has
a necessary potential for further extensions