This work presents a dynamical version of the LCAO method which includ
es extension variables, such as temperature and pressure. As an illust
rative example, the method is applied to the temperature response of s
mall silicon clusters. A comparison is made with molecular dynamics ca
lculations using standard interatomic potentials.