Constructing feed-forward artificial neural networks to fit potential energy surfaces for molecular simulation of high-temperature gas flows.

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.

Multiple translational temperature model and its shock structure solution.

Within a shock layer, the translational motion of the gas is more energetic in the direction perpendicular to the shock front than in the direction parallel to the shock. To predict this translational nonequilibrium, a generalized gas-kinetic BGK model is proposed. With the adaptation of a continuum gas distribution function truncated up to the Navier-Stokes order and a generalized particle collision time, this newly constructed kinetic model is used in the monatomic gas shock structure calculations, where two translational temperatures inside the shock layer are well captured.