These atomic energies are summed to produce N total energies. QRNN is an example of an HDNNP and utilizes two separate NN for each atom type that is supported. The first such NN takes modified Behler-Parinello type symmetry functions as input and predicts atomic electronegativities. These then enter into a simplified charge equilibration equation that predicts atomic charges. These charges are determined recursively by using them as features into the same NN. We only perform two recursive steps and then use the final charges as features, in addition to geometric symmetry functions and charge-weighted radial symmetry functions, as input to a NN which predicts atomic energies. Finally we sum the atomic energies, together with a damped coulomb interaction based on the predicted atomic charges and Grimme's D2 dispersion correction to predict a total energy.