A density functional study of Rh13

A density functional study was performed for the Rh13 cluster using the linear combination of Gaussian-type orbitals density functional theory (LCGTO-DFT) approach. The calculations employed both the local density approximation (LDA) as well as the generalized gradient approximation (GGA) in combination with a quasi-relativistic effective core potential (QECP). Initial structures for the geometry optimization were taken along Born–Oppenheimer molecular dynamics (BOMD) trajectories. The BOMD trajectories were performed at different temperatures and considered different potential energy surfaces (PES). As a result, several hundred isomers of the Rh13 cluster in different spin multiplicities were optimized with the aim to determine the lowest energy structures. All geometry optimizations were performed without any symmetry restriction. A vibrational analysis was performed to characterize these isomers. Structural parameters, relative stability energy, harmonic frequencies, binding energy, and most relevant Kohn–Sham (KS) molecular orbitals are reported. The obtained results are compared with available data from the literature. This study predicts a low symmetry biplanarlike structure as the ground-state structure of Rh13 with 11 unpaired electrons. This isomer was first noticed by inspection of first-principle Born–Oppenheimer molecular dynamics (BOMD) simulations between 300 and 600 K. This represents the most extensive theoretical study on the ground-state structure of the Rh13 cluster and underlines the importance of BOMD simulations to fully explore the PES landscapes of complicated systems.