Ab initio Hartree–Fock self‐consistent‐field molecular dynamics study of structure and dynamics of Li8

FIRST PRINCIPLES MOLECULAR DYNAMICS INVOLVING EXCITED STATES AND NONADIABATIC TRANSITIONS

Extensions of traditional molecular dynamics to excited electronic states and non-Born–Oppenheimer dynamics are reviewed focusing on applicability to chemical reactions of large molecules, possibly in condensed phases. The latter imposes restrictions on both the level of accuracy of the underlying electronic structure theory and the treatment of nonadiabaticity. This review, therefore, exclusively deals with ab initio "on the fly" molecular dynamics methods. For the same reason, mainly mixed quantum-classical approaches to nonadiabatic dynamics are considered.

Quantum Effects on the Caloric Curves of Lithium Clusters

A tight-binding quantum Hamiltonian and an empirical embedded-atom model (EAM) potential are used to get insight into the finite-temperature behavior of small Lin clusters, n = 8, 20, and 55. Exchange Monte Carlo simulations provide an extensive sampling of configuration space, including the putative global minimum and many relevant isomers. The heat capacities obtained from the classical simulations are corrected for low-temperature quantum delocalization using the Pitzer-Gwinn approximation. Alternatively, the caloric curves are estimated from the database of local minima using the quantum harmonic superposition approximation. While the two atomistic models predict qualitatively similar features, including some premelting effects in Li20 but none in Li55, strong variations are observed in the melting temperatures, the EAM potential giving unexpectedly low values.

Specific heat and Lindemann-like parameter of metallic clusters: Mono- and polyvalent metals

The Brownian-type molecular dynamics simulation is revisited and applied to study the thermal and geometric properties of four mono- and two polyvalent metallic clusters. For the thermal property, we report the specific heat at constant volume CV and study the solid-liquid-like transition by scrutinizing its characteristic. For the geometric property, we calculate the root mean square relative bond-length fluctuation delta as a function of increasing temperature. The thermal change in delta reflects the movement of atoms and hence is a relevant parameter in understanding the phase transition in clusters. The simulated results for the CV of alkali and aluminum clusters whose ground state structures exhibit icosahedral symmetry generally show one phase transition. In contrast, the tetravalent lead is quite often seen to exhibit two phase transitions, a premelting process followed by a progressive melting. In connection with the premelting scenario, it is found here that those (magic number) clusters identified to be of lesser stability (among other stable ones) according to the second energy difference are clusters showing a greater possibility of undergoing premelting process. This energy criterion applies to aluminum clusters nAl=28 and 38. To delve further into the thermal behavior of clusters, we have analyzed also the thermal variation of deltaT and attempted to correlate it with CV(T). It turns out that the premelting (if exist) and melting temperatures of the smaller size clusters (n less, similar 50) extracted from CV do not always agree quantitatively with that deduced from delta.