Kinetic rates for thermochemical nonequilibrium models are generally computed from quasiclassical trajectory (QCT) calculations on accurate ab initio potential energy surfaces (PES). In this article, we use a feed-forward artificial neural network (ANN) to fit existing single-point energies for N_{2}+N_{2} interactions [Bender et al., J. Chem. Phys. 143, 054304 (2015)JCPSA60021-960610.1063/1.4927571] to construct a PES suitable for molecular simulation of high-temperature gas flows. We then perform detailed comparisons with a widely used N_{4} PES that was built using the permutation invariant polynomials (PIP) method. Specific physical considerations in the construction of the ANN for this application are detailed. Translation, rotation, and permutation invariance are precisely satisfied by mapping the interatomic distances onto a set of permutation invariant inputs, known as fundamental invariants (FI) that generate the permutation invariant polynomial ring. The diatomic energy is imposed by decomposing the total potential energy into a sum of a two-body and a many-body energy contribution. To obtain the correct dynamical behavior with the most basic, yet computationally efficient ANN, spurious long-distance interactions must be removed to avoid incorrect physical behavior at the dissociation threshold. We use a simple apodization function to smoothly taper off to zero any residual many-body interaction at large separations. Both accuracy and performance of the FI-ANN PES are assessed. QCT calculations are used to compute dissociation probabilities and vibrational energy distributions at various equilibrium temperatures. Excellent agreement with the results obtained from the PIP PES is found. For our test case, the ANN PES is also significantly more computationally efficient than the PIP PES at comparable root-mean-square error levels.