A cooling schedule for simulated annealing specifies successive effective temperatures and how many proposals to attempt at each. For example, with , using the geometric sequence, can define stages of equal chosen length. Changing changes the transition matrix: if the state probability distribution is , then , which generally differs from . Time spent at one temperature must be compared with its mixing time before equilibrium intuition is justified.
In the five-state barrier example for simulated annealing, fixed gives a finite irreducible Markov chain, so is eventually visited with probability one. That unbounded-time statement gives no fixed finite deadline. Over any finite number of trials starting at , the probability of staying at throughout is
Each factor is the rejection probability at . This is an explicit failure event even when every is positive: the best visited objective function value remains , above the global minimum value . More exploratory proposals can help cross a barrier; cooling before those moves occur can preserve a poor local minimum. Successful asymptotic schedules need additional analysis of the allowed moves and barriers.
Historically, Scott Kirkpatrick, Charles Daniel Gelatt Jr. and Mario P. Vecchi connected statistical mechanics with mathematical optimization in 1983. They adapted the Metropolis–Hastings algorithm to cooling searches with symmetric proposals and illustrated computer-design problems.[ref]
Mario P. Vecchi Created 2026-10-18 Updated 2026-11-05
Mario P. Vecchi is an electrical engineer whose work has included semiconductor research, broadband communications and computational optimization. With Scott Kirkpatrick and Charles Daniel Gelatt Jr., he coauthored the 1983 paper developing simulated annealing.[ref]
The paper connects statistical mechanics with component placement and wiring: the cost of a proposed configuration plays the role of energy, and controlled random changes explore alternatives. His later work also included optical switching and wavelength-division multiplexing in broadband networks, illustrating a career spanning physical devices and communication systems.[ref]
Scott Kirkpatrick Created 2026-10-18 Updated 2026-11-05
The method borrows a rule from the Metropolis algorithm: a proposed increase in cost can be accepted with probability , where is an adjustable control parameter in the same units as the cost. Allowing occasional uphill changes can escape a local minimum; reducing progressively makes them less likely. The original work applied this physical analogy to electronic component placement and wiring.[ref]
Figure 1.
Portrait of Scott Kirkpatrick
. Photograph from Hebrew University of Jerusalem's official researcher profile. Source.
The Sherrington-Kirkpatrick model couples every distinct pair of Ising spins using independent standard normal random variables for , with . Here and have units of energy. David Sherrington and Scott Kirkpatrick introduced this scaling in 1975; their original replica-symmetric calculation gave an unphysical negative low-temperature entropy. The model and that initial proposed solution must be distinguished.[ref]