  This study reports the lowest energy structure of bare Cu$_{13}$ nanoclusters
as a pair of enantiomers for temperatures ranging from 20 to 1200 K. Moreover,
we compute the enantiomerization energy for the interconversion from
$\mathcal{M}$ to $\mathcal{P}$ structures in the chiral putative global minimum
for temperatures ranging from 20 to 1300 K. Additionally, employing statistical
thermodynamics and nanothermodynamics, we compute the Boltzmann Probability for
each particular isomer as a function of temperature. To achieve that, we
explore the free energy surface of the Cu$_{13}$ cluster, employing a genetic
algorithm coupled with density functional theory and statistical
thermodynamics. Moreover, we discuss the energetic ordering of structures
computed at the DFT level and compared to high level DLPNO-CCSD(T1) reference
energies, and we present chemical bonding analysis using the AdNDP method in
the chiral putative global minimum. Based on the computed relative abundances,
our results show that the chiral putative global minimum strongly dominates for
temperatures ranging from 20 to 1100 K.
