A methodology based on graph theory to optimize the operating temperatures of cryogenic computing systems is presented in this chapter. Each circuit unit is modeled as an edge in a graph with power and delay characteristics, enabling a search for the most energy-efficient temperature assignments under delay constraints. By including thermal interactions and interconnections, the approach identifies the optimal system configurations, as demonstrated with a quantum computing case study.