RESUMO
We present a Monte Carlo method for the direct evaluation of the difference between the free energies of two crystal structures. The method is built on a lattice-switch transformation that maps a configuration of one structure onto a candidate configuration of the other by "switching" one set of lattice vectors for the other, while keeping the displacements with respect to the lattice sites constant. The sampling of the displacement configurations is biased, multicanonically, to favor paths leading to gateway arrangements for which the Monte Carlo switch to the candidate configuration will be accepted. The configurations of both structures can then be efficiently sampled in a single process, and the difference between their free energies evaluated from their measured probabilities. We explore and exploit the method in the context of extensive studies of systems of hard spheres. We show that the efficiency of the method is controlled by the extent to which the switch conserves correlated microstructure. We also show how, microscopically, the procedure works: the system finds gateway arrangements which fulfill the sampling bias intelligently. We establish, with high precision, the differences between the free energies of the two close packed structures (fcc and hcp) in both the constant density and the constant pressure ensembles.
RESUMO
We show that the continuous phase space of a hard particle system can be mapped onto a discrete but infinite phase space. For three pointlike particles confined to a ring, the evolution of the system maps onto a chaotic walk on a hexagonal lattice. This facilitates direct measurement of the departure of the system from its original configuration. In special cases of mass ratios the phase space becomes closed and finite (nonergodic). There are qualitative differences between this chaotic walk and a random walk, in particular a more rapid sampling of phase space.
RESUMO
The structure of polydisperse hard sphere fluids, in the presence of a wall, is studied by the Rosenfeld density functional theory. Within this approach, the local excess free energy depends on only four combinations of the full set of density fields. The case of continuous polydispersity thereby becomes tractable. We predict, generically, an oscillatory size segregation close to the wall, and connect this, by a perturbation theory for narrow distributions, with the reversible work for changing the size of one particle in a monodisperse reference fluid.
RESUMO
Recent high-pressure work has suggested that elemental barium forms a high-pressure self-hosting structure (Ba IV) involving two "types" of barium atom. Uniquely among reported elemental structures it cannot be described by a single crystalline lattice, instead involving two interpenetrating incommensurate lattices. In this Letter we report pseudopotential calculations demonstrating the stability and the potentially disordered nature of the "guest" structure. Using band structures and nearly free electron theory we relate the appearance of Ba IV to an instability in the close-packed structure, demonstrate that it has a zero energy vibrational mode, and speculate about the structure's stability in other divalent elements.