We compute the electric dipole moment d_n of the neutron from a fully
dynamical simulation of lattice QCD with 2+1 flavors of clover fermions and
nonvanishing theta term. The latter is rotated into the pseudoscalar density in
the fermionic action using the axial anomaly. To make the action real, the
vacuum angle theta is taken to be purely imaginary. The physical value of d_n
is obtained by analytic continuation. We find d_n = -3.8(2)(9) x 10^{-16}
[theta e cm], which, when combined with the experimental limit on d_n, leads to
the upper bound theta < 7.6 x 10^{-11}.Comment: 12 pages, 8 figures, matches PRL published versio