Lattice parameters for bulk rutile TiO2 are optimized through iterative minimization of lattice stress and atomic displacements until the external pressure is reduced to below 0.5 kbar. Electronic and ionic relaxations continue until residual changes in total energy fall below 10 -5 and 10 -4 eV, respectively. Under these conditions, the atomic forces of the converged bulk TiO2 structure are all below 0.01 eV/Å. We employed Crystal Orbital Hamilton Population (COHP) analysis using the LOBSTER code to understand the chemical bonding between Rh and the surrounding atoms in the support.