We employed radial (G 2 ) and angular (G 5 ) symmetry functions, totaling 70 symmetry functions (16 radial and 54 angular), as parameterized in previous work. For these symmetry functions, we chose a cutoff radius of 6 Å based on testing different values conducted in our work. Secondly, the appropriate level of network complexity to accurately represent the underlying physics without overfitting the data was determined by testing various number of hidden layers and node counts in each of them. It was found that a network configuration with two hidden layers, each containing 17 nodes, reduced the root mean square error (RMSE) to below 22 meV/atom. Finally, the output layer, comprising a single node, yields the energy of atom in the system. For our ANN with two hidden layers and 17 nodes in each hidden layer and 70 symmetry function in the input layer, the total energy of the atom is obtained as