We present a lattice calculation of the neutron and proton electric dipole moments (EDM's) with N f = 2 + 1 flavors of domain-wall fermions. The neutron and proton EDM form factors are extracted from three-point functions at the next-to-leading order in the θ vacuum of QCD. In this computation, we use pion masses of 0.33 and 0.42 GeV and 2.7 fm 3 lattices with Iwasaki gauge action and a 0.17 GeV pion and 4.6 fm 3 lattice with I-DSDR gauge action, all generated by the RBC and UKQCD collaborations. The all-mode-averaging technique enables an efficient and high statistics calculation. Chiral behavior of lattice EDM's is discussed in the context of baryon chiral perturbation theory. In addition, we also show numerical evidence on the relationship of threeand two-point correlation functions with the local topological charge distribution.