We present a calculation of the change of free energy of a solid surface upon bending of the solid. It is based on extracting the surface stress through a molecular dynamics simulation of a bent slab by using a generalized stress theorem formula, and subsequent integration of the stress with respect to strain as a function of bending curvature. The method is exemplified by obtaining and comparing free energy changes with curvature of various reconstructed Au(001) surfaces.
The equations for calculating diffraction profiles for bent crystals are revisited for both meridional and sagittal bending. Two approximated methods for computing diffraction profiles are treated: multilamellar and Penning-Polder. A common treatment of crystal anisotropy is included in these models. The formulation presented is implemented into the XOP package, completing and updating the crystal module that simulates diffraction profiles for perfect, mosaic and now distorted crystals by elastic bending.
The total energy differences between various SiC polytypes (3C, 6H, 4H, 2H, 15R and 9R) were calculated using the full-potential linear muffin-tin orbital method using the Perdew-Wang-(91) generalized gradient approximation to the exchange-correlation functional in the density functional method. Numerical convergence versus k-point sampling and basis set completeness are demonstrated to be better than 1 meV/atom. The parameters of several generalized anisotropic next-nearest-neighbor Ising models are extracted and their significance and consequences for epitaxial growth are discussed.
We devise automated workflows for the calculation of Helmholtz and Gibbs free energies and their temperature and pressure dependence and provide the corresponding computational tools. We employ non-equilibrium thermodynamics for evaluating the free energy of solid and liquid phases at a given temperature and reversible scaling for computing free energies over a wide range of temperatures, including the direct integration of $P$-$T$ coexistence lines. By changing the chemistry and the interatomic potential, alchemical and upscaling free energy calculations are possible. Several examples illustrate the accuracy and efficiency of our implementation.
We study the excess free energy due to phase coexistence of fluids by Monte Carlo simulations using successive umbrella sampling in finite LxLxL boxes with periodic boundary conditions. Both the vapor-liquid phase coexistence of a simple Lennard-Jones fluid and the coexistence between A-rich and B-rich phases of a symmetric binary (AB) Lennard-Jones mixture are studied, varying the density rho in the simple fluid or the relative concentration x_A of A in the binary mixture, respectively. The character of phase coexistence changes from a spherical droplet (or bubble) of the minority phase (near the coexistence curve) to a cylindrical droplet (or bubble) and finally (in the center of the miscibility gap) to a slab-like configuration of two parallel flat interfaces. Extending the analysis of M. Schrader, P. Virnau, and K. Binder [Phys. Rev. E 79, 061104 (2009)], we extract the surface free energy gamma (R) of both spherical and cylindrical droplets and bubbles in the vapor-liquid case, and present evidence that for R -> Infinity the leading order (Tolman) correction for droplets has sign opposite to the case of bubbles, consistent with the Tolman length being independent on the sign of curvature. For the symmetric binary mixture the expected non-existence of the Tolman length is confirmed. In all cases {and for a range of radii} R relevant for nucleation theory, gamma(R) deviates strongly from gamma (Infinity) which can be accounted for by a term of order gamma(Infinity)/gamma(R)-1 ~ 1/R^2. Our results for the simple Lennard-Jones fluid are also compared to results from density functional theory and we find qualitative agreement in the behavior of gamma(R) as well as in the sign and magnitude of the Tolman length.
Top-down fabricated GaN nanowires, 250 nm in diameter and with various heights, have been used to experimentally determine the evolution of strain along the vertical direction of 1-dimensional objects. X-ray diffraction and photoluminescence techniques have been used to obtain the strain profile inside the nanowires from their base to their top facet for both initial compressive and tensile strains. The relaxation behaviors derived from optical and structural characterizations perfectly match the numerical results of calculations based on a continuous media approach. By monitoring the elastic relaxation enabled by the lateral free-surfaces, the height from which the nanowires can be considered strain-free has been estimated. Based on this result, NWs sufficiently high to be strain-free have been coalesced to form a continuous GaN layer. X-ray diffraction, photoluminescence and cathodoluminescence clearly show that despite the initial strain-free nanowires template, the final GaN layer is strained.