Simulated annealing with adaptive cooling rates
Mariia Karabin1, Steven J Stuart1
1Department of Chemistry, Clemson University, Clemson, South Carolina 29634, USA.
Abstract:
As one of the most robust global optimization methods, simulated annealing has received considerable attention with many variations that attempt to improve the cooling schedule. This paper introduces a variant of molecular dynamics-based simulated annealing that is useful for optimizing atomistic structures, and makes use of the heat capacity of the system, determined on the fly during optimization, to adaptively control the cooling rate. This adaptive cooling approach is demonstrated to be more computationally efficient than classical simulated annealing when applied to Lennard-Jones clusters. The increase in efficiency is approximately a factor of two for clusters with 25-40 atoms, and improves as the size of the system increases.
Related Concept Videos
02:55Temperature Dependence on Reaction Rate
Atoms, molecules, or ions must collide before they can react with each other. Atoms must be close together to form chemical bonds. This premise is the basis for a theory that explains many observations regarding chemical kinetics, including factors affecting reaction rates.
The collision theory is based on the postulates that (i) the reaction rate is proportional to the rate of reactant collisions, (ii) the reacting species collide in an orientation allowing contact between...
01:27Thermal expansion and Thermal stress: Problem Solving
To solve the problem, first, identify the known and unknown quantities. The initial length (L) of the bridge is 1275 m, the coefficient of linear expansion (α) for steel is 12 x 10-6/°C, and the change in temperature (ΔT) is 55...
02:44Heating and Cooling Curves
For instance, the addition of heat raises the temperature of a solid; the amount of heat absorbed depends on the heat capacity of the solid (q = mcsolidΔT). According to thermochemistry, the relation between the amount of heat absorbed or released by a substance, q, and its...
01:20Maxwell-Boltzmann Distribution: Problem Solving
This distribution function f(v) is defined by saying that the expected number N (v1,v2) of particles with speeds between v1 and v2 is given by
01:18Adiabatic Processes for an Ideal Gas
01:29Mechanistic Models: Compartment Models in Algorithms for Numerical Problem Solving
In individual population analyses, different algorithms are employed, such as Cauchy's method, which uses a...