We determine the topological susceptibility of the gauge configurations generated by lattice simulations using 2 flavors of optimal domain-wall fermion on the 16 3 32 lattice with length 16 in the fifth dimension, at the inverse lattice spacing a -1 1.9 GeV. Using the adaptive thick-restart Lanczos algorithm, we project the low-lying eigenmodes of the overlap Dirac operator, and obtain the topological charge Q t of each gauge configuration in eight ensembles with pion masses in the range 230 -580 MeV. From our result of Q t , we compute the topological susceptibility and the second normalized cumulant. Our result of the topological susceptibility agrees with the sea-quark mass dependence predicted by the next-to-leading order chiral perturbation theory, and provides a determination of both the chiral condensate and the pion decay constant.