The descriptors of each atom flowed into a DNN of 5 hidden layers with decreasing numbers of neurons, (512, 256, 128, 32, 8), which mapped the descriptors into the output layer, which produces the atomic energy E i . The non-linear activation function was taken to be the hyperbolic tangent. The forces on each atom were computed as negative derivatives with respect to position. The loss function was taken as