Free-energy functional for freezing transitions: hard-sphere systems freezing into crystalline and amorphous structures

Revised Enskog theory and molecular dynamics simulations of the viscosities and thermal conductivity of the hard-sphere fluid and crystal

Hard-sphere (HS) shear, longitudinal, cross, and bulk viscosities and the thermal conductivity are obtained by molecular dynamics (MD) simulations, covering the entire density range from the dilute fluid to the solid crystal near close-packing. The transport coefficient data for the HS crystal are largely new and display, unlike for the fluid, a surprisingly simple behavior in that they can be represented well by a simple function of the density compressibility factor. In contrast to the other four transport coefficients (which diverge), the bulk viscosity in the solid is quite small and decreases rapidly with increasing density, tending to zero in the close-packed limit. The so-called cross viscosity exhibits a different behavior to the other viscosities, in being negative over the entire solid range, and changes sign from negative to positive on increasing the density in the fluid phase. The extent to which the viscosity tensor and thermal conductivity of the HS crystal can be represented by revised Enskog theory (RET) is investigated. The RET expressions are sums of an instantaneous (I), a kinetic (K), and a so-called α part. The I part of the transport coefficients evaluated directly by MD are statistically indistinguishable from those of the corresponding kinetic theory (Enskog and RET) expressions. For the K part the integral over the spatial two-particle distribution function at contact was determined and the α part was estimated using the direct correlation function and density functional theory approximations. All three parts were determined in this work which allowed the accuracy of RET for solid systems to be assessed rigorously. It is found that in the case of the thermal conductivity the predictions of RET are in excellent agreement with the MD results. Also, for the shear viscosity the agreement over the entire solid phase is quite good but is considerably worse for the three remaining viscosities in the solid phase.

Super-Arrhenius behavior of molecular glass formers

A theory is developed to calculate values of the potential-energy barriers to structural relaxation in molecular glass formers from the data of static pair-correlation function. The barrier height is shown to increase due to an increase in the number of "stable bonds" a particle forms with its neighbors and the energy of each bond as liquids move deeper into the supercooled (supercompressed) region. We present results for a system of hard spheres and compare calculated values of the structural relaxation time with experimental and simulation results.

Dependence of the configurational entropy on amorphous structures of a hard-sphere fluid

The free energy of a hard-sphere fluid for which the average energy is trivial signifies how its entropy changes with packing. The packing η_{f} at which the free energy of the crystalline state becomes lower than that of the disordered fluid state marks the freezing point. For packing fractions η>η_{f} of the hard-sphere fluid, we use the modified weighted density functional approximation to identify metastable free energy minima intermediate between uniform fluid and crystalline states. The distribution of the sharply localized density profiles, i.e., the inhomogeneous density field ρ(x) characterizing the metastable state is primarily described by a pair function g_{s}(η/η_{0}). η_{0} is a structural parameter such that for η=η_{0} the pair function is identical to that for the Bernal random structure. The configurational entropy S_{c} of the metastable hard-sphere fluid is calculated by subtracting the corresponding vibrational entropy from the total entropy. The extrapolated S_{c} vanishes as η→η_{K} and η_{K} is in agreement with other works. The dependence of η_{K} on the structural parameter η_{0} is obtained.

Density-functional theory for fluid-solid and solid-solid phase transitions

We develop a theory to describe solid-solid phase transitions. The density functional formalism of classical statistical mechanics is used to find an exact expression for the difference in the grand thermodynamic potentials of the two coexisting phases. The expression involves both the symmetry conserving and the symmetry broken parts of the direct pair correlation function. The theory is used to calculate phase diagram of systems of soft spheres interacting via inverse power potentials u(r)=ε(σ/r)^{n}, where parameter n measures softness of the potential. We find that for 1/n<0.154 systems freeze into the face centered cubic (fcc) structure while for 1/n≥0.154 the body-centred-cubic (bcc) structure is preferred. The bcc structure transforms into the fcc structure upon increasing the density. The calculated phase diagram is in good agreement with the one found from molecular simulations.

Fluid-solid transition in simple systems using density functional theory

A free energy functional for a crystal which contains both the symmetry-conserved and symmetry-broken parts of the direct pair correlation function has been used to investigate the fluid-solid transition in systems interacting via purely repulsive Weeks-Chandler-Anderson Lennard-Jones potential and the full Lennard-Jones potential. The results found for freezing parameters for the fluid-face centred cubic crystal transition are in very good agreement with simulation results. It is shown that although the contribution made by the symmetry broken part to the grand thermodynamic potential at the freezing point is small compared to that of the symmetry conserving part, its role is crucial in stabilizing the crystalline structure and on values of the freezing parameters.

Correlation functions in liquids and crystals: free-energy functional and liquid-to-crystal transition

A free-energy functional for a crystal that contains both the symmetry-conserved and symmetry-broken parts of the direct pair-correlation function has been used to investigate the crystallization of fluids in three dimensions. The symmetry-broken part of the direct pair-correlation function has been calculated using a series in ascending powers of the order parameters and which contains three- and higher-body direct correlation functions of the isotropic phase. It is shown that a very accurate description of freezing transitions for a wide class of potentials is found by considering the first two terms of this series. The results found for freezing parameters including the structure of the frozen phase for fluids interacting via the inverse power potential u(r)=ε(σ/r)(n) for n ranging from 4 to ∞ are in very good agreement with simulation results. It is found that for n>6.5 the fluid freezes into a face-centered cubic (fcc) structure while for n≤6 the body-centered cubic (bcc) structure is preferred. The fluid-bcc-fcc triple point is found to be at 1/n=0.158, which is in good agreement with simulation result.

Freezing of a two-dimensional fluid into a crystalline phase: density functional approach

A free-energy functional for a crystal proposed by Singh and Singh [Europhys. Lett. 88, 16005 (2009)] which contains both the symmetry conserved and symmetry broken parts of the direct pair correlation function has been used to investigate the crystallization of a two-dimensional fluid. The results found for fluids interacting via the inverse power potential u(r)=ε(σ/r)(n) for n=3,6, and 12 are in good agreement with experimental and simulation results. The contribution made by the symmetry broken part to the grand thermodynamic potential at the freezing point is found to increase with the softness of the potential. Our results explain why the Ramakrishnan-Yussouff [Phys. Rev. B 19, 2775 (1979)] free-energy functional gave good account of freezing transitions of hard-core potentials but failed for potentials that have soft core and/or attractive tail.

Description of hard-sphere crystals and crystal-fluid interfaces: a comparison between density functional approaches and a phase-field crystal model

In materials science the phase-field crystal approach has become popular to model crystallization processes. Phase-field crystal models are in essence Landau-Ginzburg-type models, which should be derivable from the underlying microscopic description of the system in question. We present a study on classical density functional theory in three stages of approximation leading to a specific phase-field crystal model, and we discuss the limits of applicability of the models that result from these approximations. As a test system we have chosen the three-dimensional suspension of monodisperse hard spheres. The levels of density functional theory that we discuss are fundamental measure theory, a second-order Taylor expansion thereof, and a minimal phase-field crystal model. We have computed coexistence densities, vacancy concentrations in the crystalline phase, interfacial tensions, and interfacial order parameter profiles, and we compare these quantities to simulation results. We also suggest a procedure to fit the free parameters of the phase-field crystal model. Thereby it turns out that the order parameter of the phase-field crystal model is more consistent with a smeared density field (shifted and rescaled) than with the shifted and rescaled density itself. In brief, we conclude that fundamental measure theory is very accurate and can serve as a benchmark for the other theories. Taylor expansion strongly affects free energies, surface tensions, and vacancy concentrations. Furthermore it is phenomenologically misleading to interpret the phase-field crystal model as stemming directly from Taylor-expanded density functional theory.