Phase coexistence of cluster crystals: beyond the Gibbs phase rule

We report a study of the phase behavior of multiple-occupancy crystals through simulation. We argue that in order to reproduce the equilibrium behavior of such crystals it is essential to treat the number of lattice sites as a constraining thermodynamic variable. The resulting free-energy calculations thus differ considerably from schemes used for single-occupancy lattices. Using our approach, we obtain the phase diagram and the bulk modulus for a generalized exponential model that forms cluster crystals at high densities. We compare the simulation results with existing theoretical predictions. We also identify two types of density fluctuations that can lead to two sound modes and evaluate the corresponding elastic constants.

Finite size effects in the presence of a chemical potential: A study in the classical nonlinear O(2) sigma model

In the presence of a chemical potential, the physics of level crossings leads to singularities at zero temperature, even when the spatial volume is finite. These singularities are smoothed out at a finite temperature but leave behind nontrivial finite size effects which must be understood in order to extract thermodynamic quantities using Monte Carlo methods, particularly close to critical points. We illustrate some of these issues using the classical nonlinear O(2) sigma model with a coupling β and chemical potential μ on a 2+1-dimensional Euclidean lattice. In the conventional formulation this model suffers from a sign problem at nonzero chemical potential and hence cannot be studied with the Wolff cluster algorithm. However, when formulated in terms of the worldline of particles, the sign problem is absent, and the model can be studied efficiently with the "worm algorithm." Using this method we study the finite size effects that arise due to the chemical potential and develop an effective quantum mechanical approach to capture the effects. As a side result we obtain energy levels of up to four particles as a function of the box size and uncover a part of the phase diagram in the (β,μ) plane. © 2010 The American Physical Society.

Fermion Mass Generation without Spontaneous Symmetry Breaking

The conventional mechanism of fermion mass generation in the Standard Model involves Spontaneous Symmetry Breaking (SSB). In this thesis, we study an alternate mechanism for the generation of fermion masses that does not require SSB, in the context of lattice field theories. Being inherently strongly coupled, this mechanism requires a non-perturbative approach like the lattice approach.

In order to explore this mechanism, we study a simple lattice model with a four-fermion interaction that has massless fermions at weak couplings and massive fermions at strong couplings, but without any spontaneous symmetry breaking. Prior work on this type of mass generation mechanism in 4D, was done long ago using either mean-field theory or Monte-Carlo calculations on small lattices. In this thesis, we have developed a new computational approach that enables us to perform large scale quantum Monte-Carlo calculations to study the phase structure of this theory. In 4D, our results confirm prior results, but differ in some quantitative details of the phase diagram. In contrast, in 3D, we discover a new second order critical point using calculations on lattices up to size $ 60^3$. Such large scale calculations are unprecedented. The presence of the critical point implies the existence of an alternate mechanism of fermion mass generation without any SSB, that could be of interest in continuum quantum field theory.

Predictions of the Pt(8)Ti phase in unexpected systems.

The binary A(8)B phase (prototype Pt(8)Ti) has been experimentally observed in 11 systems. A high-throughput search over all the binary transition intermetallics, however, reveals 59 occurrences of the A(8)B phase: Au(8)Zn(dagger), Cd(8)Sc(dagger), Cu(8)Ni(dagger), Cu(8)Zn(dagger), Hg(8)La, Ir(8)Os(dagger), Ir(8)Re, Ir(8)Ru(dagger), Ir(8)Tc, Ir(8)W(dagger), Nb(8)Os(dagger), Nb(8)Rh(dagger), Nb(8)Ru(dagger), Nb(8)Ta(dagger), Ni(8)Fe, Ni(8)Mo(dagger)*, Ni(8)Nb(dagger)*, Ni(8)Ta*, Ni(8)V*, Ni(8)W, Pd(8)Al(dagger), Pd(8)Fe, Pd(8)Hf, Pd(8)Mn, Pd(8)Mo*, Pd(8)Nb, Pd(8)Sc, Pd(8)Ta, Pd(8)Ti, Pd(8)V*, Pd(8)W*, Pd(8)Zn, Pd(8)Zr, Pt(8)Al(dagger), Pt(8)Cr*, Pt(8)Hf, Pt(8)Mn, Pt(8)Mo, Pt(8)Nb, Pt(8)Rh(dagger), Pt(8)Sc, Pt(8)Ta, Pt(8)Ti*, Pt(8)V*, Pt(8)W, Pt(8)Zr*, Rh(8)Mo, Rh(8)W, Ta(8)Pd, Ta(8)Pt, Ta(8)Rh, V(8)Cr(dagger), V(8)Fe(dagger), V(8)Ir(dagger), V(8)Ni(dagger), V(8)Pd, V(8)Pt, V(8)Rh, and V(8)Ru(dagger) ((dagger) = metastable, * = experimentally observed). This is surprising for the wealth of new occurrences that are predicted, especially in well-characterized systems (e.g., Cu-Zn). By verifying all experimental results while offering additional predictions, our study serves as a striking demonstration of the power of the high-throughput approach. The practicality of the method is demonstrated in the Rh-W system. A cluster-expansion-based Monte Carlo model reveals a relatively high order-disorder transition temperature.