In this work, we extend recent inverse statistical-mechanical methods developed for many-particle systems to the case of spin systems. For simplicity, we focus in this initial study on the two-state Ising model with radial spin-spin interactions of finite range (i.e., extending beyond nearest-neighbor sites) on the square lattice under periodic boundary conditions. Our interest herein is to find the optimal set of shortest-range pair interactions within this family of Hamiltonians, whose corresponding ground state is a targeted spin configuration such that the difference in energies between the energetically closest competitor and the target is maximized. For an exhaustive list of competitors, this optimization problem is solved exactly using linear programming. The possible outcomes for a given target configuration can be organized into the following three solution classes: unique (nondegenerate) ground state (class I), degenerate ground states (class II), and solutions not contained in the previous two classes (class III). We have chosen to study a general family of striped-phase spin configurations comprised of alternating parallel bands of up and down spins of varying thicknesses and a general family of rectangular block checkerboard spin configurations with variable block size, which is a generalization of the classic antiferromagnetic Ising model. Our findings demonstrate that the structurally anisotropic striped phases, in which the thicknesses of up- and down-spin bands are equal, are unique ground states for isotropic short-ranged interactions. By contrast, virtually all of the block checkerboard targets are either degenerate or fall within class III solutions. The degenerate class II spin configurations are identified up to a certain block size. We also consider other target spin configurations with different degrees of global symmetries and order. Our investigation reveals that the solution class to which a target belongs depends sensitively on the nature of the target radial spin-spin correlation function. In the future, it will be interesting to explore whether such inverse statistical-mechanical techniques could be employed to design materials with desired spin properties.
In this work, we extend recent inverse statistical-mechanical methods developed for many-particle systems to the case of spin systems. For simplicity, we focus in this initial study on the two-state Ising model with radial spin-spin interactions of finite range (i.e., extending beyond nearest-neighbor sites) on the square lattice under periodic boundary conditions. Our interest herein is to find the optimal set of shortest-range pair interactions within this family of Hamiltonians, whose corresponding ground state is a targeted spin configuration such that the difference in energies between the energetically closest competitor and the target is maximized. For an exhaustive list of competitors, this optimization problem is solved exactly using linear programming. The possible outcomes for a given target configuration can be organized into the following three solution classes: unique (nondegenerate) ground state (class I), degenerate ground states (class II), and solutions not contained in the previous two classes (class III). We have chosen to study a general family of striped-phase spin configurations comprised of alternating parallel bands of up and down spins of varying thicknesses and a general family of rectangular block checkerboard spin configurations with variable block size, which is a generalization of the classic antiferromagnetic Ising model. Our findings demonstrate that the structurally anisotropic striped phases, in which the thicknesses of up- and down-spin bands are equal, are unique ground states for isotropic short-ranged interactions. By contrast, virtually all of the block checkerboard targets are either degenerate or fall within class III solutions. The degenerate class II spin configurations are identified up to a certain block size. We also consider other target spin configurations with different degrees of global symmetries and order. Our investigation reveals that the solution class to which a target belongs depends sensitively on the nature of the target radial spin-spin correlation function. In the future, it will be interesting to explore whether such inverse statistical-mechanical techniques could be employed to design materials with desired spin properties.
In the first paper of this series [DiStasio, Jr., Marcotte, Car, Stillinger, and Torquato, Phys. Rev. B 88, 134104 (2013)], we applied inverse statistical-mechanical techniques to study the extent to which targeted spin configurations on the square lattice can be ground states of finite-ranged radial spin-spin interactions. In this sequel, we enumerate all of the spin configurations within a unit cell on the one-dimensional integer lattice and the two-dimensional square lattice up to some modest size under periodic boundary conditions. We then classify these spin configurations into those that can or cannot be unique classical ground states of the aforementioned radial pair spin interactions and found the relative occurrences of these ground-state solution classes for different system sizes. As a result, we also determined the minimal radial extent of the spin-spin interaction potentials required to stabilize those configurations whose ground states are either unique or degenerate (i.e., those sharing the same radial spin-spin correlation function). This enumeration study has established that unique ground states are not limited to simple target configurations. However, we also found that many simple target spin configurations cannot be unique ground states.
Chapter 1: The statistical mechanics of particles with shapes on a one-dimensional lattice is investigated in the context of the s = 1 Ising chain with uniform nearest-neighbor coupling, quadratic single-site potential, and magnetic field, which supports four distinct ground states: : ↑↓↑↓ ···⟩, : ∘ ∘ ···⟩, : ↑↑ ··· ⟩, : ↑ ∘ ↑ ∘ ··· ⟩. The complete spectrum is generated from each ground state by particles from a different set of six or seven species. Particles and elements of pseudo-vacuum are characterized by motifs (patterns of several consecutive site variables). Particles are floating objects that can be placed into open slots on the lattice. Open slots are recognized as permissible links between motifs. The energy of a particle varies between species but is independent of where it is placed. Placement of one particle changes the open-slot configuration for particles of all species. This statistical interaction is encoded in a generalized Pauli principle, from which the multiplicity of states for a given particle combination is determined and used for the exact statistical mechanical analysis. Particles from all species belong to one of four categories: compacts, hosts, tags, or hybrids. Compacts and hosts find open slots in segments of pseudo-vacuum. Tags find open slots inside hosts. Hybrids are tags with hosting capability. In the taxonomy of particles proposed here, 'species' is indicative of structure and 'category' indicative of function. The hosting function splits the Pauli principle into exclusion and accommodation parts. Near phase boundaries, the state of the Ising chain at low temperature is akin to that of miscible or immiscible liquids with particles from one species acting as surfactant molecules. Chapter 2: The s = 1/2 Ising chain with uniform nearest-neighbor and next-nearest-neighbor coupling is used to construct a system of floating particles characterized by motifs of up to six consecutive local spins. The spin couplings cause the assembly of particles which, in turn, remain free of interaction energies even at high density. All microstates are configurations of particles from one of three different sets, excited from pseudo-vacua associated with ground states of periodicities one, two, and four. The motifs of particles and elements of pseudo-vacuum interlink in two shared site variables. The statistical interaction between particles is encoded in a generalized Pauli principle, describing how the placement of one particle modifies the options for placing further particles. In the statistical mechanical analysis arbitrary energies can be assigned to all particle species. The entropy is a function of the particle populations. The statistical interaction specifications are transparently built into that expression. The energies and structures of the particles alone govern the ordering at low temperature. Under special circumstances the particles can be replaced by more fundamental particles with shorter motifs that interlink in only one shared site variable. Structures emerge from interactions on two levels: particles with shapes from coupled spins and long-range ordering tendencies from statistically interacting particles with shapes. Chapter 3: The s = [special characters omitted] Ising spin chain with uniform nearest-neighbor coupling, quadratic single-site potential, and magnetic field is shown to be equivalent to a system of 17 species of particles with internal structure. The same set of particles (with different energies) is shown to generate the spectrum of the s = ½ Ising chain with dimerized nearest-neighbor coupling. The particles are free of interaction energies even at high densities. The mutual exclusion statistics of particles from all species is determined by their internal structure and encoded in a generalized Pauli principle. The exact statistical mechanical analysis can be performed for thermodynamically open or closed systems and with arbitrary energies assigned to all particle species. Special circumstances make it possible to merge two or more species into a single species. All traits that distinguish the original species become ignorable. The particles from the merged species are effectively indistinguishable and obey modified exclusion statistics. Different mergers may yield the same end-product, implying that the inverse process (splitting any species into subspecies) is not unique. In a macroscopic system of two merged species at thermal equilibrium, the concentrations of the original species satisfy a functional relation governed by their mutual statistical interaction. That relation is derivable from an extremum principle. In the Ising context the system is open and the particle energies depend on the Hamiltonian parameters. Simple models of polymerization and solitonic paramagnetism each represent a closed system of two species that can transform into each other. Here they represent distinguishable traits with different energies of the same physical particle.
In spin systems, geometrical frustration describes the impossibility of minimizing simultaneously all the interactions in a Hamiltonian, often giving rise to macroscopic ground-state degeneracies and emergent low-temperature physics. In this thesis, combining tensor network (TN) methods to Monte Carlo (MC) methods and ground-state energy lower bound approaches, we study two-dimensional frustrated classical Ising models. In particular, we focus on the determination of the residual entropy in the presence of farther-neighbor interactions in kagome lattice Ising antiferromagnets (KIAFM). In general, using MC to determine the residual entropy is a significant challenge requiring ad-hoc updates, a precise evaluation of the energy at all temperatures to allow for thermodynamic integration, and a good control of the finite-size scaling behavior. As an alternative, we turn to TNs; however, we argue that, in the presence of frustration and macroscopic ground-state degeneracy, standard algorithms fail to converge at low temperatures on the usual TN formulation of partition functions. Inspired by methods for constructing ground-state energy lower bounds, we propose a systematic way to find the ground-state local rule using linear programming. Characterizing the rules as tiles that can be tessellated to form ground states of the model gives rise to a natural contractible TN formulation of the partition function. This method provides a direct access to the ground-state properties of frustrated models and, in particular, allows an extremely precise determination of their residual entropy. We then study two models inspired by artificial spin systems on the kagome lattice with out-of-plane (OOP) anisotropy. The first model is motivated by experiments on an array of chirally coupled nanomagnets. We argue that the farther-neighbor to nearest-neighbor couplings ratios in this system are much smaller than in the dipolar case, J2/J1 being of the order of 2%. A comparison of the experimental correlations with the results of extensive TN and MC simulations shows that (1) the experimental second- and third-neighbor correlations are inverted as compared to those of a pure nearest-neighbor model at equilibrium (even with a magnetic field), and (2) second-neighbor couplings as small as 1% of the nearest-neighbor couplings will affect the spin-spin correlations even at fairly high temperatures. Motivated by dipolar coupled artificial spin systems, we turn to the progressive lifting of the ground-state degeneracy of the KIAFM. We provide a detailed study of the ground-state phases of this model with up to third neighbor interactions, for arbitrary J2, J3 such that J1 >> J2, J3, obtaining exact results for the ground-state energies. When all couplings are antiferromagnetic, we exhibit three macroscopically degenerate ground-state phases and establish their residual entropy using our TN approach. Furthermore, in the phase corresponding to the dipolar KIAFM truncated to third neighbors, we use the ground-state tiles to establish the existence of a mapping to the ground-state manifold of the triangular Ising antiferromagnet.
This is a study of finite lattice methods applied to quantum spin models on one and two dimensional lattices.The motivation for the use of these techniques is provided by placing them in the context of the wider field of approximation methods in the statistical mechanics of spin models.After reviewing previous finite lattice calculations on quantum spin models, the theory needed to efficiently generate and analyse finite lattice data is given.Finite lattice calculations on three quantum spin models are then presented.The first calculation provides an estimate of the transverse susceptibility of the one dimensional antiferromagnetic XY model, by extrapolation from calculations on rings of up to 10 sites.Good convergence is found for KB T ~ 0•4 and these results are compared with recent c -axis susceptibility measurements on the compound Cs 2 CoC1 4 .Chapters 5 and 6 study the quantum hamiltonian analog of the two dimensional axial next nearest neighbour Ising model, A perturbation analysis using non degenerate Rayleigh-Schrodingerperturbation theory is followed by a mass gap scaling analysis.These analyses provide accurate phase boundaries for, and ellucidate several interesting features of, the phase diagram of this model.In particular, the model is shown to exhibit a phase whose modulation varies continuously with a change in the interaction parameters in the hamiltonian.This "modulated" phase region lies between a normal ordered phase and a paramagnetic phase.Evid~nce is given that the transition from the paramagnetic to modulated phase is exponential in nature.The third model studied is a spin 1 biquadratic Heisenberg model with hamiltonian,<ij> (Jl s .. s. +J2(~ . .~.)2). 1 J 1 J Accurate estimates of the ground state energy in both one and two dimensions are found by calculating this quantity on 4,6,8 and 10 site lattices.Comment is made about the ground state phase diagram, and about the existence of intermediate spin states and states with a complicated magnetic ordering.The thesis concludes by commenting on the performance of finite lattice methods in the three models studied and by outlining certain directions for future advance., ACKNOWLEVGEMENTS I would like to express
Classical and quantum spin systems are widely used models in both experimental and theoretical condensed matter physics. In many materials, the electronic interactions can be difficult to model exactly. However, in some insulators (Mott insulators), the magnetic (spin) interactions can be captured well with spin-only models. Several models are studied in this thesis. First, I report the solution of a long-standing issue in a classical frustrated spin model (i.e., where the quantum effects are neglected and the complexity is due to competing interactions): the nature of the thermal phase transition to a stripe state in the two dimensional (2D) J1-J2 model. Here J1 and J2 are nearest- and next-nearest-neighbor couplings. Monte Carlo simulations with single-spin updates are used for the calculations, and an extended-ensemble method, a generalization of the Wang-Landau algorithm, is also developed and tested. I focus on the study of "weak universality" behavior (continuously varying critical exponents, with one of the exponents staying fixed), for which I show a correspondence with a known class of conformal field theories with charge c = 1. Next, moving to quantum spins, to shed light on magnetic systems studied experimentally and to investigate new types of quantum states of interest in developing theories of quantum magnetism, I study the S = 1/2 Heisenberg model on the 2D square lattice with added six-spin interactions (the so-called J-Q3 model) as well as a set of 3D quantum antiferromagnets on dimerized lattices. Here I use the stochastic series expansion quantum Monte Carlo method. In the study of the J-Q3 model, I report on a similar weak-universality behavior as in the classical J1-J2 model, but with a mapping to a different known class of c = 1 conformal field theories. The critical behavior of the system again shows continuously changing exponents, in a way which corresponds to a gradual weakening of Z4 symmetry-breaking to an emergent U(1) symmetry. In the study of dimerized antiferromagnets, I report on a universal behavior of the Néel temperature TN, which can be related to ground state parameters independently of the microscopic interaction details in several different models.
For generalized 2D Ising model in an external magnetic field with the interaction of nearest neighbors, next nearest neighbors, all kinds of triple interactions and the quadruple interaction the formulas for finding free energy per lattice site in the thermodynamic limit were derived on a certain set of exact disordered solutions depending on seven parameters. Lattice models are considered with boundary conditions with a shift (similar to helical ones), and with cyclic closure of the set of all points in natural ordering. The elementary transfer matrix with nonnegative matrix elements are constructed. On the set of disorder solutions the largest eigenvalue of the transfer matrix is constant for every size of considering planar lattice, and, in particular, in the thermodynamic limit. Free energy per lattice site in the thermodynamic limit is expressed through the natural logarithm of the largest eigenvalue of transfer matrix. This largest eigenvalue can be found for a special form of eigenvector with positive components. The numerical example show the existence of nontrivial solutions of the resulting systems of equations. The system of equations and the value of free energy in the thermodynamic limit will remain the same for 2D generalized Ising models with Hamiltonians, in which the values of two (out of four) neighboring maximal spins in the natural ordering are replaced by the values of the spins at any other two lattice points adjacent in the natural ordering, this significantly expands the set of models having disordered exact solutions. The high degree of symmetry and inductive construction of the components of the eigenvectors, which disappear when going beyond the framework of the obtained set of exact solutions, is an occasion to search for phase transitions in the vicinity of this set of disordered solutions.