Orbital-Free Molecular Dynamics Simulations of Melting in Na8 and Na20: Melting in Steps

Abstract

The melting-like transitions of Na8 and Na20 are investigated by ab initio constant energy molecular dynamics simulations, using a variant of the Car-Parrinello method which employs an explicit electronic kinetic energy functional of the density, thus avoiding the use of one-particle orbitals. Several melting indicators are evaluated in order to determine the nature of the various transitions, and compared with other simulations. Both Na8 and Na20 melt over a wide temperature range. For Na8, a transition is observed to begin at approx. 110 K, between a rigid phase and a phase involving isomerizations between the different permutational isomers of the ground state structure. The ``liquid'' phase is completely established at approx. 220 K. For Na20, two transitions are observed: the first, at approx. 110 K, is associated with isomerization transitions between those permutational isomers of the ground state structure which are obtained by interchanging the positions of the surface-like atoms; the second, at approx. 160 K, involves a structural transition from the ground state isomer to a new set of isomers with the surface molten. The cluster is completely ``liquid'' at approx. 220 K.

Revised version, accepted for publication in J. Chem. Phys. The changes include longer simulations for the Na20 microcluster, a more complete comparison to previous theoretical results, and the discussion of some technical details of the method applied

References in corpus (1)

• Ab initio molecular dynamics using density based energy functionals: application to ground state geometries of some small clusters

Cited by in corpus (14)

• Dynamical properties of liquid Al near melting. An orbital-free molecular dynamics study
• Negative heat capacity of sodium clusters
• Ionic and electronic structure of sodium clusters up to N=59
• Thermal expansion in small metal clusters and its impact on the electric polarizability
• First principles investigation of finite-temperature behavior in small sodium clusters
• Orbital-Free DFT Correctly Models Quantum Dots When Asymptotics, Nonlocality and Nonhomogeneity Are Accounted For
• On the premelting features in sodium clusters
• Impurity effects on the melting of Ni clusters
• Finite temperature behavior of impurity doped Lithium cluster {\em viz} Li${}_6$Sn
• Nonlocal Kinetic Energy Density Functionals for Isolated Systems via Local Density Approximation Kernels
• Thermodynamics of Na_8 and Na_{20} clusters studied with ab-initio electronic structure methods
• An orbital-free molecular dynamics study of melting in K_20, K_55, K_92, K_142, Rb_55 and Cs_55 clusters
• Exchange Monte Carlo for Molecular Simulations with Monoelectronic Hamiltonians
• The electronic structures, the equilibrium geometries and finite temperature properties of Na_n (n=39-55)