The surface energy, particularly for nanoparticles, is one of the most important quantities in understanding the thermodynamics of particles. Therefore, it is astonishing that there is still great uncertainty about its value. The uncertainty increases if one questions its dependence on particle size. Different approaches, such as classical thermodynamics calculations, molecular dynamics simulations, and ab initio calculations, exist to predict this quantity. Generally, considerations based on classical thermodynamics lead to the prediction of decreasing values of the surface energy with decreasing particle size. This phenomenon is caused by the reduced number of next neighbors of surface atoms with decreasing particle size, a phenomenon that is partly compensated by the reduction of the binding energy between the atoms with decreasing particle size. Furthermore, this compensating effect may be expected by the formation of a disordered or quasi-liquid layer at the surface. The atomistic approach, based either on molecular dynamics simulations or ab initio calculations, generally leads to values with an opposite tendency. However, it is shown that this result is based on an insufficient definition of the particle size. A more realistic definition of the particle size is possible only by a detailed analysis of the electronic structure obtained from initio calculations. Except for minor variations caused by changes in the structure, only a minor dependence of the surface energy on the particle size is found. The main conclusion of this work is that surface energy values for the equivalent bulk materials should be used if detailed data for nanoparticles are not available.

Motivated by often contradictory literature reports on the dependence of the surface energy of gold nanoparticles on the variety of its size and shape, we performed an atomistic study combining molecular mechanics and ab initio calculations. We show that, in the case of Au nanocubes, their surface energy converges to the value for $\left(001\right)$ facets of bulk crystals. A fast convergence to a single valued surface energy is predicted also for nanospheres. However, the value of the surface energy is larger in this case than that of any low-index surface facet of bulk Au crystal. This fact can be explained by the complex structure of the surface with an extensive number of broken bonds due to edge and corner atoms. A similar trend was obtained also for the case of cuboctahedrons. Since the exact surface area of the nanoparticles is an ill-defined quantity, we have introduced the surface-induced excess energy and discuss this quantity as a function of (i) number of atoms forming the nano-object or (ii) characteristic size of the nano-object. In case (i), a universal power-law behaviour was obtained independent of the nanoparticle shape. Importantly, we show that the size-dependence of the surface energy is hugely reduced, if the surface area correction is considered due to its expansion by the electronic cloud, a phenomenon specifically important for small nanoparticles.