-
Structure of Periodic Crystals and Quasicrystals in Ultrathin Film Ba-Ti-O
Authors:
Eric Cockayne,
Marek Mihalkovic,
Christopher L. Henley
Abstract:
We model the remarkable thin-film Ba-Ti-O structures formed by heat treatment of an initial perovskite BaTiO$_3$ thin film on a Pt(111) surface. All structures contain a rumpled Ti-O network with all Ti threefold coordinated with O, and with Ba occupying the larger. mainly Ti$_7$O$_7$, pores. The quasicrystal structure is a simple decoration of three types of tiles: square, triangle and 30 degree…
▽ More
We model the remarkable thin-film Ba-Ti-O structures formed by heat treatment of an initial perovskite BaTiO$_3$ thin film on a Pt(111) surface. All structures contain a rumpled Ti-O network with all Ti threefold coordinated with O, and with Ba occupying the larger. mainly Ti$_7$O$_7$, pores. The quasicrystal structure is a simple decoration of three types of tiles: square, triangle and 30 degree rhombus, with edge lengths 6.85 Angstroms, joined edge-to-edge in a quasicrystalline pattern; observed crystalline Ba-Ti-O structures are built from these and other tiles. Simulated STM images reproduce the patterns seen experimentally, and identify the bright protrusions as Ba atoms. The models are consistent with all experimental observations.
△ Less
Submitted 7 January, 2016; v1 submitted 21 October, 2015;
originally announced October 2015.
-
Schwinger boson mean field perspective on emergent spins in diluted Heisenberg antiferromagnets
Authors:
Shivam Ghosh,
Hitesh J. Changlani,
Christopher L. Henley
Abstract:
Using an adaptation of Schwinger boson mean field theory (SBMFT) for nonuniform systems, we study the nature of low-energy spin excitations on the square and Bethe lattice at their percolation threshold. The optimal SBMFT parameters are interpreted as on-site potentials and pairing amplitudes, which enables an explanation of why emergent local moments develop in this system on dilution [L. Wang an…
▽ More
Using an adaptation of Schwinger boson mean field theory (SBMFT) for nonuniform systems, we study the nature of low-energy spin excitations on the square and Bethe lattice at their percolation threshold. The optimal SBMFT parameters are interpreted as on-site potentials and pairing amplitudes, which enables an explanation of why emergent local moments develop in this system on dilution [L. Wang and A. W. Sandvik, Phys. Rev. Lett. 97, 117204 (2006); H.J. Changlani et al., ibid 111, 157201 (2013)] and why the corresponding single particle frequencies are driven to anomalously low values. We discuss how our mean field calculations suggest the strong link between the presence of sublattice imbalance and long range antiferromagnetic order, and why linear spin wave theory is inadequate for capturing this relation. Within the SBMFT framework, we also extract an energy scale for the interaction between emergent moments, which show qualitative agreement with many-body calculations.
△ Less
Submitted 3 August, 2015; v1 submitted 21 May, 2015;
originally announced May 2015.
-
Phase diagram of the Kondo lattice model on the Kagome lattice
Authors:
Shivam Ghosh,
Patrick O' Brien,
Christopher L. Henley,
Michael J. Lawler
Abstract:
We consider the potential for novel forms of magnetism arising from the subtle interplay between electrons and spins in the under-screened kagome Kondo lattice model. At weak coupling, we show that incommensurate non-coplanar multi-wave vector magnetic orders arise at nearly all fillings and that this results from Fermi surface effects that introduces competing interactions between the spins. At s…
▽ More
We consider the potential for novel forms of magnetism arising from the subtle interplay between electrons and spins in the under-screened kagome Kondo lattice model. At weak coupling, we show that incommensurate non-coplanar multi-wave vector magnetic orders arise at nearly all fillings and that this results from Fermi surface effects that introduces competing interactions between the spins. At strong coupling, we find that such complex order survives near half filling despite the presence of ferromagnetism at all other fillings. We show this arises due to state selection among a massive degeneracy of states at infinite coupling. Finally, we show that at intermediate filling, only commensurate orders seem to survive. But these orders still include non-coplanar magnetism. So, the mere presence of both local moments and itinerant electrons enables complex orders to form unlike any currently observed in kagome materials.
△ Less
Submitted 3 August, 2015; v1 submitted 20 July, 2014;
originally announced July 2014.
-
Density-matrix based numerical methods for discovering order and correlations in interacting systems
Authors:
Christopher L. Henley,
Hitesh J. Changlani
Abstract:
We review recently introduced numerical methods for the unbiased detection of the order parameter and/or dominant correlations, in many-body interacting systems, by using reduced density matrices. Most of the paper is devoted to the "quasi-degenerate density matrix" (QDDM) which is rooted in Anderson's observation that the degenerate symmetry-broken states valid in the thermodynamic limit, are man…
▽ More
We review recently introduced numerical methods for the unbiased detection of the order parameter and/or dominant correlations, in many-body interacting systems, by using reduced density matrices. Most of the paper is devoted to the "quasi-degenerate density matrix" (QDDM) which is rooted in Anderson's observation that the degenerate symmetry-broken states valid in the thermodynamic limit, are manifested in finite systems as a set of low-energy "quasi-degenerate" states (in addition to the ground state). This method, its original form due to Furukawa et al.[Phys. Rev. Lett. 96, 047211 (2006)], is given a number of improvements here, above all the extension from two-fold symmetry breaking to arbitrary cases. This is applied to two test cases (1) interacting spinless hardcore bosons on the triangular lattice and (2) a spin-1/2 antiferromagnetic system at the percolation threshold. In addition, we survey a different method called the "correlation density matrix", which detects (possibly long-range) correlations only from the ground state, but using the reduced density matrix from a cluster consisting of two spatially separated regions.
△ Less
Submitted 15 July, 2014;
originally announced July 2014.
-
Semi-classical spin dynamics of the antiferromagnetic Heisenberg model on the kagome lattice
Authors:
Mathieu Taillefumier,
Julien Robert,
Christopher L. Henley,
Roderich Moessner,
Benjamin Canals
Abstract:
We investigate the dynamical properties of the classical antiferromagnetic Heisenberg model on the kagome lattice using a combination of Monte Carlo and molecular dynamics simulations. We find that frustration induces a distribution of timescales in the cooperative paramagnetic regime (i.e. far above the onset of coplanarity), as recently reported experimentally in deuterium jarosite. At lower tem…
▽ More
We investigate the dynamical properties of the classical antiferromagnetic Heisenberg model on the kagome lattice using a combination of Monte Carlo and molecular dynamics simulations. We find that frustration induces a distribution of timescales in the cooperative paramagnetic regime (i.e. far above the onset of coplanarity), as recently reported experimentally in deuterium jarosite. At lower temperature, when the coplanar correlations are well established, we show that the weath- ervane loop fluctuations control the system relaxation : the time distribution observed at higher temperatures splits into two distinct timescales associated with fluctuations in the plane and out of the plane of coplanarity. The temperature and wave vector dependences of these two components are qualitatively consistent with loops diffusing in the entropically dominated free energy landscape. Numerical results are discussed and compared with the $O(N)$ model and recent experiments for both classical and quantum realizations of the kagome magnets.
△ Less
Submitted 31 March, 2014;
originally announced March 2014.
-
Orientational interaction and ordering of Cd_4 tetrahedra in a quasicrystal approximant
Authors:
Woosong Choi,
Christopher L. Henley,
Marek Mihalkovic
Abstract:
We model the quasicrystal-related structure CaCd$_6$, a bcc packing of icosahedral clusters containing tetrahedra which undergo orientational orderings at T<100 K. We use general schemes to evaluate an effective Hamltonian for inter-tetrahedron orientations, based on all-atom relaxations, either in terms of discrete cluster orientations, or of continuous rotation angles. The effective Hamiltonian…
▽ More
We model the quasicrystal-related structure CaCd$_6$, a bcc packing of icosahedral clusters containing tetrahedra which undergo orientational orderings at T<100 K. We use general schemes to evaluate an effective Hamltonian for inter-tetrahedron orientations, based on all-atom relaxations, either in terms of discrete cluster orientations, or of continuous rotation angles. The effective Hamiltonian is used in Monte Carlo simulations to find the (complex) ground state ordering pattern as a function of pressure. A preliminary investigation of thermal transitions found (in part of the pressure range) two different first-order transitions occurring below 100 K.
△ Less
Submitted 21 December, 2012;
originally announced December 2012.
-
NMR relaxation in spin ice due to diffusing emergent monopoles
Authors:
Christopher L. Henley
Abstract:
At low temperatures, spin dynamics in ideal spin ice is due mainly to dilute, thermally excited magnetic monopole excitations. I consider how these will affect the dynamics of a nuclear spin (the same theory applies to muon spin resonance if implanted muons do not diffuse). Up to the time scale for nearby monopoles to be rearranged, a stretched-exponential form of the relaxation functions is expec…
▽ More
At low temperatures, spin dynamics in ideal spin ice is due mainly to dilute, thermally excited magnetic monopole excitations. I consider how these will affect the dynamics of a nuclear spin (the same theory applies to muon spin resonance if implanted muons do not diffuse). Up to the time scale for nearby monopoles to be rearranged, a stretched-exponential form of the relaxation functions is expected. I work out the expected exponent in that exponential and the formulas for the $T_1$ (longitudinal) and $T_2$ (dephasing) relaxations, as a function of the monopole density. Experimental NMR is incompatible with the predictions and I suggest is due to magnetic impurities.
△ Less
Submitted 30 October, 2012;
originally announced October 2012.
-
Classical Topological Order in Abelian and Non-Abelian Generalized Height Models
Authors:
R. Zach Lamberty,
Stefanos Papanikolaou,
Christopher L. Henley
Abstract:
We present Monte Carlo simulations on a new class of lattice models in which the degrees of freedom are elements of an abelian or non-abelian finite symmetry group G, placed on directed edges of a two-dimensional lattice. The plaquette group product is constrained to be the group identity. In contrast to discrete gauge models (but similar to past work on height models) only elements of symmetry-re…
▽ More
We present Monte Carlo simulations on a new class of lattice models in which the degrees of freedom are elements of an abelian or non-abelian finite symmetry group G, placed on directed edges of a two-dimensional lattice. The plaquette group product is constrained to be the group identity. In contrast to discrete gauge models (but similar to past work on height models) only elements of symmetry-related subsets S of G are allowed on edges. These models have topological sectors labeled by group products along topologically non-trivial loops. Measurement of relative sector probabilities and the distribution of distance between defect pairs are done to characterize the types of order (topological or quasi-LRO) exhibited by these models. We present particular models in which fully local non-abelian constraints lead to global topological liquid properties.
△ Less
Submitted 30 May, 2013; v1 submitted 30 October, 2012;
originally announced October 2012.
-
Boson features in STM spectra of cuprate superconductors: Weak-coupling phenomenology
Authors:
Sumiran Pujari,
C. L. Henley
Abstract:
We derive the shape of the high-energy features due to a weakly coupled boson in cuprate superconductors, as seen experimentally in Bi_2 Sr_2 Ca_1 Cu_2 O_8+x (BSCCO) by Lee et al. [Nature (London) 442, 546 (2006)]. A simplified model is used of d-wave Bogoliubov quasiparticles coupled to Einstein oscillators with a momentum-independent electron-boson coupling and an analytic fitting form is derive…
▽ More
We derive the shape of the high-energy features due to a weakly coupled boson in cuprate superconductors, as seen experimentally in Bi_2 Sr_2 Ca_1 Cu_2 O_8+x (BSCCO) by Lee et al. [Nature (London) 442, 546 (2006)]. A simplified model is used of d-wave Bogoliubov quasiparticles coupled to Einstein oscillators with a momentum-independent electron-boson coupling and an analytic fitting form is derived, which allows us (a) to extract the boson mode's frequency and (b) to estimate the electron-boson coupling strength. We further calculate the maximum possible superconducting gap due to an Einstein oscillator with the extracted electron-boson coupling strength, which is found to be less than 0.2 times of the observed gap indicating at the observed boson's non-dominant role in the superconductivity's mechanism. The extracted momentum-independent electron-boson coupling parameter (that we show a posteriori to indeed be in the weak-coupling regime) is then to be interpreted as an (band-structure detail dependent weighted) average over the Brillouin zone of the actual momentum-dependent electron-boson coupling in BSCCO.
△ Less
Submitted 17 October, 2014; v1 submitted 28 October, 2012;
originally announced October 2012.
-
Ground States of the Classical Antiferromagnet on the Pyrochlore Lattice
Authors:
Matthew F. Lapa,
Christopher L. Henley
Abstract:
We study the classical ground states of the exchange-coupled Heisenberg antiferromagnet on the Pyrochlore lattice, a non-Bravais lattice made of corner-sharing tetrahedra. In particular, we map out the entire phase diagram for the case of first and second nearest neighbor interactions. In this phase diagram we find {\it four} complex non-coplanar ground states based on different ordering modes. Th…
▽ More
We study the classical ground states of the exchange-coupled Heisenberg antiferromagnet on the Pyrochlore lattice, a non-Bravais lattice made of corner-sharing tetrahedra. In particular, we map out the entire phase diagram for the case of first and second nearest neighbor interactions. In this phase diagram we find {\it four} complex non-coplanar ground states based on different ordering modes. These are the Cuboctahedral Stack state, a $<111>$ stacking of Cuboctahedral states, and three families of commensurate spiral: the Kawamura states, constructed from three different combinations of $\{\tfrac{3}{4}\tfrac{3}{4}0\}$ modes, the Double-Twist state, also constructed from $\{\tfrac{3}{4}\tfrac{3}{4}0\}$ modes, and the Multiply-Modulated Commensurate Spiral state, constructed from $\{\tfrac{3}{4}\tfrac{1}{4}\tfrac{1}{2}\}$ modes. We also briefly look at states involving the two kinds of third nearest neighbor interactions on the Pyrochlore lattice. In this region of parameter space we again find the Cuboctahedral Stack state, and we also find another non-coplanar state in the form of a new kind of Alternating Conic Spiral.
△ Less
Submitted 25 October, 2012;
originally announced October 2012.
-
Emergent Spin Excitations in a Bethe Lattice at Percolation
Authors:
Hitesh J. Changlani,
Shivam Ghosh,
Sumiran Pujari,
Christopher L. Henley
Abstract:
We study the spin 1/2 quantum Heisenberg antiferromagnet on a Bethe lattice diluted to the percolation threshold. Dilution creates areas of even/odd sublattice imbalance resulting in "dangling spins" (L. Wang and A. W. Sandvik, Phys. Rev. Lett. 97, 117204 (2006); Phys. Rev. B 81, 054417 (2010)). These collectively act as "emergent" spin 1/2 degrees of freedom and are responsible for the creation o…
▽ More
We study the spin 1/2 quantum Heisenberg antiferromagnet on a Bethe lattice diluted to the percolation threshold. Dilution creates areas of even/odd sublattice imbalance resulting in "dangling spins" (L. Wang and A. W. Sandvik, Phys. Rev. Lett. 97, 117204 (2006); Phys. Rev. B 81, 054417 (2010)). These collectively act as "emergent" spin 1/2 degrees of freedom and are responsible for the creation of a set of low lying "quasidegenerate states". Using Density Matrix Renormalization Group (DMRG) calculations, we detect the presence and location of these emergent spins. We find an effective Hamiltonian of these emergent spins, with Heisenberg interactions that decay exponentially with the distance between them.
△ Less
Submitted 8 October, 2013; v1 submitted 16 October, 2012;
originally announced October 2012.
-
Nonplanar ground states of frustrated antiferromagnets on an octahedral lattice
Authors:
Sophia R. Sklan,
Christopher L. Henley
Abstract:
We consider methods to identify the classical ground state for an exchange-coupled Heisenberg antiferromagnet on a non-Bravais lattice with interactions $J_{ij}$ to several neighbor distances. Here we apply this to the unusual "octahedral" lattice in which spins sit on the edge midpoints of a simple cubic lattice. Our approach is informed by the eigenvectors of $J_{ij}$ with largest eigenvalues. W…
▽ More
We consider methods to identify the classical ground state for an exchange-coupled Heisenberg antiferromagnet on a non-Bravais lattice with interactions $J_{ij}$ to several neighbor distances. Here we apply this to the unusual "octahedral" lattice in which spins sit on the edge midpoints of a simple cubic lattice. Our approach is informed by the eigenvectors of $J_{ij}$ with largest eigenvalues. We discovered two families of non-coplanar states: (i) two kinds of commensurate state with cubic symmetry, each having twelve sublattices with spins pointing in (1,1,0) directions in spin space (modulo a global rotation); (ii) varieties of incommensurate conic spiral. The latter family is addressed by projecting the three-dimensional lattice to a one-dimensional chain, with a basis of two (or more) sites per unit cell.
△ Less
Submitted 13 September, 2012; v1 submitted 6 September, 2012;
originally announced September 2012.
-
Heisenberg antiferromagnet on Cayley trees: low-energy spectrum and even/odd site imbalance
Authors:
Hitesh J. Changlani,
Shivam Ghosh,
Christopher L. Henley,
Andreas M. Läuchli
Abstract:
To understand the role of local sublattice imbalance in low-energy spectra of s=1/2 quantum antiferromagnets, we study the s=1/2 quantum nearest neighbor Heisenberg antiferromagnet on the coordination 3 Cayley tree. We perform many-body calculations using an implementation of the density matrix renormalization group (DMRG) technique for generic tree graphs. We discover that the bond-centered Cayle…
▽ More
To understand the role of local sublattice imbalance in low-energy spectra of s=1/2 quantum antiferromagnets, we study the s=1/2 quantum nearest neighbor Heisenberg antiferromagnet on the coordination 3 Cayley tree. We perform many-body calculations using an implementation of the density matrix renormalization group (DMRG) technique for generic tree graphs. We discover that the bond-centered Cayley tree has a quasidegenerate set of a low-lying tower of states and an "anomalous" singlet-triplet finite-size gap scaling. For understanding the construction of the first excited state from the many-body ground state, we consider a wave function ansatz given by the single-mode approximation, which yields a high overlap with the DMRG wave function. Observing the ground-state entanglement spectrum leads us to a picture of the low-energy degrees of freedom being "giant spins" arising out of sublattice imbalance, which helps us analytically understand the scaling of the finite-size spin gap. The Schwinger-boson mean-field theory has been generalized to nonuniform lattices, and ground states have been found which are spatially inhomogeneous in the mean-field parameters.
△ Less
Submitted 11 March, 2013; v1 submitted 8 August, 2012;
originally announced August 2012.
-
Order and supersymmetry at high filling zero energy states on the triangular lattice
Authors:
D. Galanakis,
C. L. Henley,
S. Papanikolaou
Abstract:
We perform exact diagonalization studies in d=2 dimensions for the Fendley and Schoutens model of hard-core and nearest-neighbor excluding fermions that displays an exact non-relativistic supersymmetry. Using clusters of all possible shapes up to 46 sites, we systematically study the behavior of the ground state phase diagram as a function of filling. We focus on the highly degenerate zero-energy…
▽ More
We perform exact diagonalization studies in d=2 dimensions for the Fendley and Schoutens model of hard-core and nearest-neighbor excluding fermions that displays an exact non-relativistic supersymmetry. Using clusters of all possible shapes up to 46 sites, we systematically study the behavior of the ground state phase diagram as a function of filling. We focus on the highly degenerate zero-energy states found at fillings between 1/7 and ~1/5. At the lower end of that interval, at filling 1/7, we explicitly show that the ground states are gapped crystals. Consistent with previous suggestions, we find that the extensive entropy of zero states peaks at a filling of ~0.178. At the higher end of the interval, we find zero energy ground states at fillings above 1/5, in contradiction to previous numerical studies and analytical suggestions; these display non-trivial amplitude degeneracies.
△ Less
Submitted 25 July, 2012; v1 submitted 10 June, 2012;
originally announced June 2012.
-
Ab-initio tiling and atomic structure for decagonal ZnMgDy quasicrystal
Authors:
M. Mihalkovic,
C. L. Henley,
J. Richmond-Decker,
M. Oxborrow
Abstract:
We discover the detailed atomic structure of $d$-MgZnY, a stable decagonal quasicrystal alloy of the layered Frank-Kasper type, and related phases, using the "tiling and decoration" approach. The atoms have invariable sites in the rectangle and triangle tiles of a 10-fold-symmetric planar tiling. To discover the lowest-energy structures, we combine the methods of density functional theory (DFT) to…
▽ More
We discover the detailed atomic structure of $d$-MgZnY, a stable decagonal quasicrystal alloy of the layered Frank-Kasper type, and related phases, using the "tiling and decoration" approach. The atoms have invariable sites in the rectangle and triangle tiles of a 10-fold-symmetric planar tiling. To discover the lowest-energy structures, we combine the methods of density functional theory (DFT) total energy calculations, empirical oscillating pair potentials (fitted to DFT), fitting effective Hamiltonians for tilings, and discovering optimum tiling structures using a nonlocal tile-reshuffling algorithm. We find a family of practically stable compounds with varying composition, including the decagonal quasicrystal and the known Mg$_4$Zn$_7$ phase these are more stable than competing icosahedral structures by a small margi
△ Less
Submitted 16 December, 2011;
originally announced December 2011.
-
Caged clusters in Al$_{11}$Ir$_4$: structural transition and insulating phase
Authors:
Marek Mihalkovič,
C. L. Henley
Abstract:
Using pair potentials fitted to an ab--initio database, combined with replica--exchange simulated annealing, we show the complex, quasicrystal-related Al$_{11}$Ir$_4$ compound contains a new version of the "pseudo-Mackay" icosahedral cluster, with non-icosahedral inner Al$_{10}$Ir and/or Al$_9$Ir clusters that exist in various orientations and account for partial occupancies in the reported struct…
▽ More
Using pair potentials fitted to an ab--initio database, combined with replica--exchange simulated annealing, we show the complex, quasicrystal-related Al$_{11}$Ir$_4$ compound contains a new version of the "pseudo-Mackay" icosahedral cluster, with non-icosahedral inner Al$_{10}$Ir and/or Al$_9$Ir clusters that exist in various orientations and account for partial occupancies in the reported structure. Two different compositions show first-order transitions to orientationally ordered phases doubling the (cubic) unit cell, which are respectively metallic and insulating.
△ Less
Submitted 23 June, 2012; v1 submitted 16 December, 2011;
originally announced December 2011.
-
Propagating left/right asymmetry in the zebrafish embryo: one-dimensional model
Authors:
Hanrong Chen,
C. L. Henley,
B. Xu
Abstract:
During embryonic development in vertebrates, left-right (L/R) asymmetry is reliably generated by a conserved mechanism: a L/R asymmetric signal is transmitted from the embryonic node to other parts of the embryo by the L/R asymmetric expression and diffusion of the TGF-$β$ related proteins Nodal and Lefty via propagating gene expression fronts in the lateral plate mesoderm (LPM) and midline. In ze…
▽ More
During embryonic development in vertebrates, left-right (L/R) asymmetry is reliably generated by a conserved mechanism: a L/R asymmetric signal is transmitted from the embryonic node to other parts of the embryo by the L/R asymmetric expression and diffusion of the TGF-$β$ related proteins Nodal and Lefty via propagating gene expression fronts in the lateral plate mesoderm (LPM) and midline. In zebrafish embryos, Nodal and Lefty expression can only occur along 3 narrow stripes that express the co-receptor \emph{one-eyed pinhead} (oep): Nodal along stripes in the left and right LPM, and Lefty along the midline. In wild-type embryos, Nodal is only expressed in the left LPM but not the right, because of inhibition by Lefty from the midline; however, bilateral Nodal expression occurs in loss-of-handedness mutants. A two-dimensional model of the zebrafish embryo predicts this loss of L/R asymmetry in oep mutants \cite{henley-xu-burdine}. In this paper, we simplify this two-dimensional picture to a one-dimensional model of Nodal and Lefty front propagation along the oep-expressing stripes. We represent Nodal and Lefty production by step functions that turn on when a linear function of Nodal and Lefty densities crosses a threshold. We do a parameter exploration of front propagation behavior, and find the existence of \emph{pinned} intervals, along which the linear function underlying production is pinned to the threshold. Finally, we find parameter regimes for which spatially uniform oscillating solutions are possible.
△ Less
Submitted 15 December, 2011;
originally announced December 2011.
-
Possible origins of macroscopic left-right asymmetry in organisms
Authors:
Christopher L. Henley
Abstract:
I consider the microscopic mechanisms by which a particular left-right (L/R) asymmetry is generated at the organism level from the microscopic handedness of cytoskeletal molecules. In light of a fundamental symmetry principle, the typical pattern-formation mechanisms of diffusion plus regulation cannot implement the "right-hand rule"; at the microscopic level, the cell's cytoskeleton of chiral fil…
▽ More
I consider the microscopic mechanisms by which a particular left-right (L/R) asymmetry is generated at the organism level from the microscopic handedness of cytoskeletal molecules. In light of a fundamental symmetry principle, the typical pattern-formation mechanisms of diffusion plus regulation cannot implement the "right-hand rule"; at the microscopic level, the cell's cytoskeleton of chiral filaments seems always to be involved, usually in collective states driven by polymerization forces or molecular motors. It seems particularly easy for handedness to emerge in a shear or rotation in the background of an effectively two-dimensional system, such as the cell membrane or a layer of cells, as this requires no pre-existing axis apart from the layer normal. I detail a scenario involving actin/myosin layers in snails and in C. elegans, and also one about the microtubule layer in plant cells. I also survey the other examples that I am aware of, such as the emergence of handedness such as the emergence of handedness in neurons, in eukaryote cell motility, and in non-flagellated bacteria.
△ Less
Submitted 2 May, 2012; v1 submitted 10 December, 2011;
originally announced December 2011.
-
Empirical oscillating potentials for alloys from ab-initio fits and the prediction of quasicrystal-related structures in the Al-Cu-Sc system
Authors:
M. Mihalkovic,
C. L. Henley
Abstract:
By fitting to a database of ab-initio forces and energies, we can extract pair potentials for alloys, with a simple six-parameter analytic form including Friedel oscillations, which give a remarkably faithful account of many complex intermetallic compounds. Furthermore, such potentials are combined with a method of discovering complex zero-temperature structures with hundreds of atoms per cell, gi…
▽ More
By fitting to a database of ab-initio forces and energies, we can extract pair potentials for alloys, with a simple six-parameter analytic form including Friedel oscillations, which give a remarkably faithful account of many complex intermetallic compounds. Furthermore, such potentials are combined with a method of discovering complex zero-temperature structures with hundreds of atoms per cell, given only the composition and the con- straint of known lattice parameters, using molecular-dynamics quenches. We apply this approach to structure prediction in the Al-Cu-Sc quasicrystal-related system.
△ Less
Submitted 30 September, 2011;
originally announced September 2011.
-
Properties of Resonating-Valence-Bond Spin Liquids and Critical Dimer Models
Authors:
Ying Tang,
Anders W. Sandvik,
Christopher L. Henley
Abstract:
We use Monte Carlo simulations to study properties of Anderson's resonating-valence-bond (RVB) spin-liquid state on the square lattice (i.e., the equal superposition of all pairing of spins into nearest-neighbor singlet pairs) and compare with the classical dimer model (CDM). The latter system also corresponds to the ground state of the Rokhsar-Kivelson quantum dimer model at its critical point. W…
▽ More
We use Monte Carlo simulations to study properties of Anderson's resonating-valence-bond (RVB) spin-liquid state on the square lattice (i.e., the equal superposition of all pairing of spins into nearest-neighbor singlet pairs) and compare with the classical dimer model (CDM). The latter system also corresponds to the ground state of the Rokhsar-Kivelson quantum dimer model at its critical point. We find that although spin-spin correlations decay exponentially in the RVB, four-spin valence-bond-solid (VBS) correlations are critical, qualitatively like the well-known dimer-dimer correlations of the CDM, but decaying more slowly (as $1/r^a$ with $a \approx 1.20$, compared with $a=2$ for the CDM). We also compute the distribution of monomer (defect) pair separations, which decay by a larger exponent in the RVB than in the CDM. We further study both models in their different winding number sectors and evaluate the relative weights of different sectors. Like the CDM, all the observed RVB behaviors can be understood in the framework of a mapping to a "height" model characterized by a gradient-squared stiffness constant $K$. Four independent measurements consistently show a value $K_{RVB} \approx 1.6 K_{CDM}$, with the same kinds of numerical evaluations of $K_{CDM}$ give results in agreement with the rigorously known value $K_{CDM}=π/16$. The background of a nonzero winding number gradient $W/L$ introduces spatial anisotropies and an increase in the effective K, both of which can be understood as a consequence of anharmonic terms in the height-model free energy, which are of relevance to the recently proposed scenario of "Cantor deconfinement" in extended quantum dimer models. We also study ensembles in which fourth-neighbor (bipartite) bonds are allowed, at a density controlled by a tunable fugacity, resulting (as expected) in a smooth reduction of K.
△ Less
Submitted 22 November, 2011; v1 submitted 29 October, 2010;
originally announced October 2010.
-
The "Coulomb phase" in frustrated systems
Authors:
Christopher L. Henley
Abstract:
The "Coulomb phase" is an emergent state for lattice models (particularly highly frustrated antiferromagnets) which have local constraints that can be mapped to a divergence-free "flux". The coarse-grained version of this flux or polarization behave analogously to electric or magnetic fields; in particular, defects at which the local constraint is violated behave as effective charges with Coulom…
▽ More
The "Coulomb phase" is an emergent state for lattice models (particularly highly frustrated antiferromagnets) which have local constraints that can be mapped to a divergence-free "flux". The coarse-grained version of this flux or polarization behave analogously to electric or magnetic fields; in particular, defects at which the local constraint is violated behave as effective charges with Coulomb interactions. I survey the derivation of the characteristic power-law correlation functions and the pinch-points in reciprocal space plots of diffuse scattering, as well as applications to magnetic relaxation, quantum-mechanical generalizations, phase transitions to long-range-ordered states, and the effects of disorder.
△ Less
Submitted 22 December, 2009;
originally announced December 2009.
-
Echolocation by Quasiparticles
Authors:
Sumiran Pujari,
C. L. Henley
Abstract:
It is shown that the local density of states (LDOS), measured in an Scanning Tunneling Microscopy (STM) experiment, at a single tip position contains oscillations as a function of Energy, due to quasiparticle interference, which is related to the positions of nearby scatterers. We propose a method of STM data analysis based on this idea, which can be used to locate the scatterers. In the case of…
▽ More
It is shown that the local density of states (LDOS), measured in an Scanning Tunneling Microscopy (STM) experiment, at a single tip position contains oscillations as a function of Energy, due to quasiparticle interference, which is related to the positions of nearby scatterers. We propose a method of STM data analysis based on this idea, which can be used to locate the scatterers. In the case of a superconductor, the method can potentially distinguish the nature of the scattering by a particular impurity.
△ Less
Submitted 12 November, 2009; v1 submitted 23 October, 2009;
originally announced October 2009.
-
Classical height models with topological order
Authors:
Christopher L. Henley
Abstract:
I discuss a family of statistical-mechanics models in which (some classes of) elements of a finite group $G$ occupy the (directed) edges of a lattice; the product around any plaquette is constrained to be the group identity $e$. Such a model may possess topological order, i.e. its equilibrium ensemble has distinct, symmetry-related thermodynamic components that cannot be distinguished by any local…
▽ More
I discuss a family of statistical-mechanics models in which (some classes of) elements of a finite group $G$ occupy the (directed) edges of a lattice; the product around any plaquette is constrained to be the group identity $e$. Such a model may possess topological order, i.e. its equilibrium ensemble has distinct, symmetry-related thermodynamic components that cannot be distinguished by any local order parameter. In particular, if $G$ is a non-abelian group, the topological order may be non-abelian. Criteria are given for the viability of particular models, in particular for Monte Carlo updates.
△ Less
Submitted 17 March, 2011; v1 submitted 23 October, 2009;
originally announced October 2009.
-
Correlation density matrices for 1- dimensional quantum chains based on the density matrix renormalization group
Authors:
W. Münder,
A. Weichselbaum,
A. Holzner,
J. von Delft,
C. L. Henley
Abstract:
A useful concept for finding numerically the dominant correlations of a given ground state in an interacting quantum lattice system in an unbiased way is the correlation density matrix. For two disjoint, separated clusters, it is defined to be the density matrix of their union minus the direct product of their individual density matrices and contains all correlations between the two clusters. We…
▽ More
A useful concept for finding numerically the dominant correlations of a given ground state in an interacting quantum lattice system in an unbiased way is the correlation density matrix. For two disjoint, separated clusters, it is defined to be the density matrix of their union minus the direct product of their individual density matrices and contains all correlations between the two clusters. We show how to extract from the correlation density matrix a general overview of the correlations as well as detailed information on the operators carrying long-range correlations and the spatial dependence of their correlation functions. To determine the correlation density matrix, we calculate the ground state for a class of spinless extended Hubbard models using the density matrix renormalization group. This numerical method is based on matrix product states for which the correlation density matrix can be obtained straightforwardly. In an appendix, we give a detailed tutorial introduction to our variational matrix product state approach for ground state calculations for 1- dimensional quantum chain models. We show in detail how matrix product states overcome the problem of large Hilbert space dimensions in these models and describe all techniques which are needed for handling them in practice.
△ Less
Submitted 5 October, 2009;
originally announced October 2009.
-
Exact ground states and correlation functions of chain and ladder models of interacting hardcore bosons or spinless fermions
Authors:
Siew-Ann Cheong,
Christopher L. Henley
Abstract:
By removing one empty site between two occupied sites, we map the ground states of chains of hardcore bosons and spinless fermions with infinite nearest-neighbor repulsion to ground states of chains of hardcore bosons and spinless fermions without nearest-neighbor repulsion respectively, and ultimately in terms of the one-dimensional Fermi sea. We then introduce the intervening-particle expansio…
▽ More
By removing one empty site between two occupied sites, we map the ground states of chains of hardcore bosons and spinless fermions with infinite nearest-neighbor repulsion to ground states of chains of hardcore bosons and spinless fermions without nearest-neighbor repulsion respectively, and ultimately in terms of the one-dimensional Fermi sea. We then introduce the intervening-particle expansion, where we write correlation functions in such ground states as a systematic sum over conditional expectations, each of which can be ultimately mapped to a one-dimensional Fermi-sea expectation. Various ground-state correlation functions are calculated for the bosonic and fermionic chains with infinite nearest-neighbor repulsion, as well as for a ladder model of spinless fermions with infinite nearest-neighbor repulsion and correlated hopping in three limiting cases. We find that the decay of these correlation functions are governed by surprising power-law exponents.
△ Less
Submitted 24 July, 2009;
originally announced July 2009.
-
Coarse-grained protein-protein stiffnesses and dynamics from all-atom simulations
Authors:
Stephen D. Hicks,
C. L. Henley
Abstract:
Large protein assemblies, such as virus capsids, may be coarse-grained as a set of rigid domains linked by generalized (rotational and stretching) harmonic springs. We present a method to obtain the elastic parameters and overdamped dynamics for these springs from all-atom molecular dynamics simulations of one pair of domains at a time. The computed relaxation times of this pair give a consisten…
▽ More
Large protein assemblies, such as virus capsids, may be coarse-grained as a set of rigid domains linked by generalized (rotational and stretching) harmonic springs. We present a method to obtain the elastic parameters and overdamped dynamics for these springs from all-atom molecular dynamics simulations of one pair of domains at a time. The computed relaxation times of this pair give a consistency check for the simulation, and (using a fluctuation-dissipation relationship) we find the corrective force needed to null systematic drifts. As a first application we predict the stiffness of an HIV capsid layer and the relaxation time for its breathing mode.
△ Less
Submitted 24 June, 2009; v1 submitted 11 March, 2009;
originally announced March 2009.
-
Effective Hamiltonians for state selection in Heisenberg antiferromagnets
Authors:
Brond E. Larson,
Christopher L. Henley
Abstract:
In frustrated antiferromagnets with isotropic exchange interactions, there is typically a manifold of degenerate classical ground states. This degeneracy is broken by the (free) energy of quantum or thermal fluctuations, or the uniform effects of bond disorder. We derive effective Hamiltonians to express each kind of selection effect, in both exact forms and convenient approximate forms. It is a…
▽ More
In frustrated antiferromagnets with isotropic exchange interactions, there is typically a manifold of degenerate classical ground states. This degeneracy is broken by the (free) energy of quantum or thermal fluctuations, or the uniform effects of bond disorder. We derive effective Hamiltonians to express each kind of selection effect, in both exact forms and convenient approximate forms. It is argued that biquadratic terms, representing the collinear-selecting effects of quantum fluctuations, should be included in classical simulations of large-$S$ frustrated magnets at low temperatures
△ Less
Submitted 28 September, 2009; v1 submitted 6 November, 2008;
originally announced November 2008.
-
Spontaneous Currents in Spinless Fermion Lattice Models at the Strong-Coupling Limit
Authors:
Sumiran Pujari,
C. L. Henley
Abstract:
What kind of lattice Hamiltonian manifestly has an ordered state with spontaneous orbital currents? We consider interacting spinless fermions on an array of square plaquettes, connected by weak hopping; the array geometry may be a 2 x 2L ladder, a 2 x 2 x 2L "tube", or a 2L x 2L square grid. At half filling, we derive an effective Hamiltonian in terms of pseudospins, of which one component repre…
▽ More
What kind of lattice Hamiltonian manifestly has an ordered state with spontaneous orbital currents? We consider interacting spinless fermions on an array of square plaquettes, connected by weak hopping; the array geometry may be a 2 x 2L ladder, a 2 x 2 x 2L "tube", or a 2L x 2L square grid. At half filling, we derive an effective Hamiltonian in terms of pseudospins, of which one component represents orbital currents, and find the conditions sufficient for orbital current long-range order. We consider spinfull variants of the aforesaid spinless models and make contact with other spinfull models in the literature purported to possess spontaneous currents.
△ Less
Submitted 24 November, 2008; v1 submitted 4 November, 2008;
originally announced November 2008.
-
Anharmonic Ground state selection in the pyrochlore antiferromagnet
Authors:
Uzi Hizi,
Christopher L. Henley
Abstract:
In the pyrochlore lattice Heisenberg antiferromagnet, for large spin length $S$, the massive classical ground state degeneracy is partly lifted by the zero-point energy of quantum fluctuations at harmonic order in spin-waves. However, there remains an infinite manifold of degenerate collinear ground states, related by a gaugelike symmetry. We have extended the spin-wave calculation to quartic or…
▽ More
In the pyrochlore lattice Heisenberg antiferromagnet, for large spin length $S$, the massive classical ground state degeneracy is partly lifted by the zero-point energy of quantum fluctuations at harmonic order in spin-waves. However, there remains an infinite manifold of degenerate collinear ground states, related by a gaugelike symmetry. We have extended the spin-wave calculation to quartic order, assuming a Gaussian variational wavefunction (equivalent to Hartree-Fock approximation). Quartic calculations \emph{do} break the harmonic-order degeneracy of periodic ground states. The form of the effective Hamiltonian describing this splitting, which depends on loops, was fitted numerically and also rationalized analytically. We find a family of states that are still almost degenerate, being split by the term from loops of length 26. We also calculated the anharmonic terms for the checkerboard lattice, and discuss why it (as well as the kagomé lattice) behave differently than the pyrochlore at anharmonic orders.
△ Less
Submitted 31 May, 2009; v1 submitted 3 November, 2008;
originally announced November 2008.
-
Possible mechanisms for initiating macroscopic left-right asymmetry in developing organisms
Authors:
Christopher L. Henley
Abstract:
How might systematic left-right (L/R) asymmetry of the body plan originate in multicellular animals (and plants)? Somehow, the microscopic handedness of biological molecules must be brought up to macroscopic scales. Basic symmetry principles suggest that the usual "biological" mechanisms -- diffusion and gene regulation -- are insufficient to implement the "right-hand rule" defining a third body…
▽ More
How might systematic left-right (L/R) asymmetry of the body plan originate in multicellular animals (and plants)? Somehow, the microscopic handedness of biological molecules must be brought up to macroscopic scales. Basic symmetry principles suggest that the usual "biological" mechanisms -- diffusion and gene regulation -- are insufficient to implement the "right-hand rule" defining a third body axis from the other two. Instead, on the cellular level, "physical" mechanisms (forces and collective dynamic states) are needed involving the long stiff fibers of the cytoskeleton. I discuss some possible scenarios; only in the case of vertebrate internal organs is the answer currently known (and even that is in dispute).
△ Less
Submitted 29 November, 2008; v1 submitted 31 October, 2008;
originally announced November 2008.
-
Long range order in the classical kagome antiferromagnet: effective Hamiltonian approach
Authors:
Christopher L. Henley
Abstract:
Following Huse and Rutenberg [Phys. Rev. B 45, 7536 (1992)], I argue the classical Heisenberg antiferromagnet on the kagomé lattice has long-range spin order of the $\sqrt{3}\times\sqrt{3}$ type (modulo gradual orientation fluctuations of the spins' plane). I start from the effective quartic Hamiltonian for the soft (out of plane) spin fluctuation modes, and treat as a perturbation those terms w…
▽ More
Following Huse and Rutenberg [Phys. Rev. B 45, 7536 (1992)], I argue the classical Heisenberg antiferromagnet on the kagomé lattice has long-range spin order of the $\sqrt{3}\times\sqrt{3}$ type (modulo gradual orientation fluctuations of the spins' plane). I start from the effective quartic Hamiltonian for the soft (out of plane) spin fluctuation modes, and treat as a perturbation those terms which depend on the discrete coplanar state. Soft mode correlations, which become the coefficients of a discrete effective Hamiltonian, are estimated analytically.
△ Less
Submitted 9 October, 2009; v1 submitted 31 October, 2008;
originally announced November 2008.
-
Matching rules from Al-Co potentials in an almost realistic model
Authors:
Sejoon Lim,
M. Mihalkovic,
C. L. Henley
Abstract:
We consider a model decagonal quasicrystal of composition Al$_{80.1}$Co$_{19.9}$ -- closely related to actual structures, and using realistic pair potentials -- on a quasilattice of candidate sites. Its ground state, according to simulations, is a Hexagon-Boat-Star tiling satisfying Penrose's matching rules. In this note, we rationalize these results in terms of the potentials; the Al-Co second-…
▽ More
We consider a model decagonal quasicrystal of composition Al$_{80.1}$Co$_{19.9}$ -- closely related to actual structures, and using realistic pair potentials -- on a quasilattice of candidate sites. Its ground state, according to simulations, is a Hexagon-Boat-Star tiling satisfying Penrose's matching rules. In this note, we rationalize these results in terms of the potentials; the Al-Co second-neighbor potential well is crucial.
△ Less
Submitted 30 August, 2008;
originally announced September 2008.
-
Penrose Matching Rules from Realistic Potentials in a Model System
Authors:
Sejoon Lim,
M. Mihalkovic,
C. L. Henley
Abstract:
We exhibit a toy model of a binary decagonal Al-Co quasicrystal -- closely related to actual structures -- in which realistic pair potentials yield a ground state which appears to perfectly implement Penrose's matching rules, for Hexagon-Boat-Star (HBS) tiles of edge 2.45 A. The second minimum of the potentials is crucial for this result.
We exhibit a toy model of a binary decagonal Al-Co quasicrystal -- closely related to actual structures -- in which realistic pair potentials yield a ground state which appears to perfectly implement Penrose's matching rules, for Hexagon-Boat-Star (HBS) tiles of edge 2.45 A. The second minimum of the potentials is crucial for this result.
△ Less
Submitted 30 August, 2008;
originally announced September 2008.
-
Depleted pyrochlore antiferromagnets
Authors:
Christopher L. Henley
Abstract:
I consider the class of "depleted pyrochlore" lattices of corner-sharing triangles, made by removing spins from a pyrochlore lattice such that every tetrahedron loses exactly one. Previously known examples are the "hyperkagome" and "kagome staircase". I give criteria in terms of loops for whether a given depleted lattice can order analogous to the kagome \sqrt{3} \times \sqrt{three} state, and a…
▽ More
I consider the class of "depleted pyrochlore" lattices of corner-sharing triangles, made by removing spins from a pyrochlore lattice such that every tetrahedron loses exactly one. Previously known examples are the "hyperkagome" and "kagome staircase". I give criteria in terms of loops for whether a given depleted lattice can order analogous to the kagome \sqrt{3} \times \sqrt{three} state, and also show how the pseudo-dipolar correlations (due to local constraints) generalize to even the random depleted case.
△ Less
Submitted 30 August, 2008;
originally announced September 2008.
-
Correlation density matrix: an unbiased analysis of exact diagonalizations
Authors:
Siew-Ann Cheong,
C. L. Henley
Abstract:
Given the ground state wavefunction for an interacting lattice model, we define a "correlation density matrix"(CDM) for two disjoint, separated clusters $A$ and $B$, to be the density matrix of their union, minus the direct product of their respective density matrices. The CDM can be decomposed systematically by a numerical singular value decomposition, to provide a systematic and unbiased way t…
▽ More
Given the ground state wavefunction for an interacting lattice model, we define a "correlation density matrix"(CDM) for two disjoint, separated clusters $A$ and $B$, to be the density matrix of their union, minus the direct product of their respective density matrices. The CDM can be decomposed systematically by a numerical singular value decomposition, to provide a systematic and unbiased way to identify the operator(s) dominating the correlations, even unexpected ones.
△ Less
Submitted 30 August, 2008;
originally announced September 2008.
-
Empirical oscillating potentials for alloys from ab-initio fits
Authors:
M. Mihalkovič,
C. L. Henley,
M. Widom,
P. Ganesh
Abstract:
By fitting to a database of ab-initio forces and energies, we can extract pair potentials for alloys, with a simple six-parameter analytic form including Friedel oscillations, which give a remarkably faithful account of many complex intermetallic compounds. As examples we show results for (crystal or quasicrystal) structure prediction and phonon spectrum for three systems: Fe--B, Al--Mg--Zn, and…
▽ More
By fitting to a database of ab-initio forces and energies, we can extract pair potentials for alloys, with a simple six-parameter analytic form including Friedel oscillations, which give a remarkably faithful account of many complex intermetallic compounds. As examples we show results for (crystal or quasicrystal) structure prediction and phonon spectrum for three systems: Fe--B, Al--Mg--Zn, and Al--Cu--Fe.
△ Less
Submitted 20 February, 2008; v1 submitted 20 February, 2008;
originally announced February 2008.
-
Frustrated minority spins in GeNi2O4
Authors:
M. Matsuda,
J. -H. Chung,
S. Park,
T. J. Sato,
K. Matsuno,
H. Aruga-Katori,
H. Takagi,
K. Kakurai,
K. Kamazawa,
Y. Tsunoda,
I. Kagomiya,
C. L. Henley,
S. -H. Lee
Abstract:
Recently, two consecutive phase transitions were observed, upon cooling, in an antiferromagnetic spinel GeNi$_2$O$_4$ at $T_{N1}=12.1$ K and $T_{N2}=11.4$ K, respectively \cite{matsuno, crawford}. Using unpolarized and polarized elastic neutron scattering we show that the two transitions are due to the existence of frustrated minority spins in this compound. Upon cooling, at $T_{N1}$ the spins o…
▽ More
Recently, two consecutive phase transitions were observed, upon cooling, in an antiferromagnetic spinel GeNi$_2$O$_4$ at $T_{N1}=12.1$ K and $T_{N2}=11.4$ K, respectively \cite{matsuno, crawford}. Using unpolarized and polarized elastic neutron scattering we show that the two transitions are due to the existence of frustrated minority spins in this compound. Upon cooling, at $T_{N1}$ the spins on the $<111>$ \kagome planes order ferromagnetically in the plane and antiferromagnetically between the planes (phase I), leaving the spins on the $<111>$ triangular planes that separate the \kagome planes frustrated and disordered. At the lower $T_{N2}$, the triangular spins also order in the $<111>$ plane (phase II). We also present a scenario involving exchange interactions that qualitatively explains the origin of the two purely magnetic phase transitions.
△ Less
Submitted 23 August, 2007;
originally announced August 2007.
-
Effective Hamiltonians for large-S pyrochlore antiferromagnet
Authors:
Uzi Hizi,
Christopher L. Henley
Abstract:
The pyrochlore lattice Heisenberg antiferromagnet has a massive classical ground state degeneracy. We summarize three approximation schemes, valid for large spin length $S$, to capture the (partial) lifting of this degeneracy when zero-point quantum fluctuations are taken into account; all three are related to analytic loop expansions. The first is harmonic order spin waves; at this order, there…
▽ More
The pyrochlore lattice Heisenberg antiferromagnet has a massive classical ground state degeneracy. We summarize three approximation schemes, valid for large spin length $S$, to capture the (partial) lifting of this degeneracy when zero-point quantum fluctuations are taken into account; all three are related to analytic loop expansions. The first is harmonic order spin waves; at this order, there remains an infinite manifold of degenerate collinear ground states, related by a gauge-like symmetry. The second is anharmonic (quartic order) spin waves, using a self-consistent approximation; the harmonic-order degeneracy is split, but (within numerical precision) some degeneracy may remain, with entropy still of order $L$ in a system of $L^3$ sites. The third is a large-N approximation, a standard and convenient approach for frustrated antiferromagnets; however, the large-N result contradicts the harmonic order at $O(S)$ hence must be wrong (for large $S$).
△ Less
Submitted 17 November, 2006;
originally announced November 2006.
-
Energy-based Structure Prediction for d(Al70Co20Ni10)
Authors:
Nan Gu,
M. Mihalkovic,
C. L. Henley
Abstract:
We use energy minimization principles to predict the structure of a decagonal quasicrystal - d(AlCoNi) - in the Cobalt-rich phase. Monte Carlo methods are then used to explore configurations while relaxation and molecular dynamics are used to obtain a more realistic structure once a low energy configuration has been found. We find five-fold symmetric decagons 12.8 A in diameter as the characteri…
▽ More
We use energy minimization principles to predict the structure of a decagonal quasicrystal - d(AlCoNi) - in the Cobalt-rich phase. Monte Carlo methods are then used to explore configurations while relaxation and molecular dynamics are used to obtain a more realistic structure once a low energy configuration has been found. We find five-fold symmetric decagons 12.8 A in diameter as the characteristic formation of this composition, along with smaller pseudo-five-fold symmetric clusters filling the spaces between the decagons. We use our method to make comparisons with a recent experimental approximant structure model from Sugiyama et al (2002).
△ Less
Submitted 7 March, 2006; v1 submitted 7 March, 2006;
originally announced March 2006.
-
Structure prediction for cobalt-rich decagonal AlCoNi from pair potentials
Authors:
Nan Gu,
M. Mihalkovic,
C. L. Henley
Abstract:
A systematic, decoration-based technique to discover the atomic structure of a decagonal quasicrystal, given pair potentials and experimentally measured lattice constants, is applied to the ``basic'' cobalt-rich decagonal Al-Co-Ni quasicrystal. First lattice-gas Monte Carlo simulations are performed, assuming the atomic sites are vertices of a rhombus tiling with edge 2.45 A. This phase is found…
▽ More
A systematic, decoration-based technique to discover the atomic structure of a decagonal quasicrystal, given pair potentials and experimentally measured lattice constants, is applied to the ``basic'' cobalt-rich decagonal Al-Co-Ni quasicrystal. First lattice-gas Monte Carlo simulations are performed, assuming the atomic sites are vertices of a rhombus tiling with edge 2.45 A. This phase is found to be dominated by 13 A diameter decagon-shaped clusters, each with a pentagon of Co atoms at the center. These, and another smaller cluster, decorate vertices of a ``binary tiling'' with rhombus edge 10.4 A. Further simulations with a restricted site list show that Al arrangements on the borders of the 13 A decagon cluster form Hexagon, Boat, and Star tiles with edge 2.45 A; they indicate specific sites for Co versus Ni atoms, and how the structure adapts to small composition changes. In the second half of the paper, relaxation (augmented by molecular dynamics annealing) is used to obtain realistic structures. The dominant new feature is a set of linear ``channels'' attractive to Al atoms and running transverse to the layers. Each is typically occupied by three atoms in four layers, implying puckering and a spontaneous period doubling to c $\approx$ 8 A. Puckering favors pentagonal long range order of the cluster orientations. Our simulation captures most features of the related $W$-AlCoNi crystal, except for its pentagonal bipyramid motif.
△ Less
Submitted 3 February, 2006;
originally announced February 2006.
-
An irreversible growth model for virus capsid assembly
Authors:
Stephen D. Hicks,
C. L. Henley
Abstract:
We model the spontaneous assembly of a capsid (a virus's closed outer shell) from many copies of identical units, using entirely irreversible steps and only information local to the growing edge. Our model is formulated in terms of (i) an elastic Hamiltonian with stretching and bending stiffness and a spontaneous curvature, and (ii) a set of rate constants for addition of new units or bonds. An…
▽ More
We model the spontaneous assembly of a capsid (a virus's closed outer shell) from many copies of identical units, using entirely irreversible steps and only information local to the growing edge. Our model is formulated in terms of (i) an elastic Hamiltonian with stretching and bending stiffness and a spontaneous curvature, and (ii) a set of rate constants for addition of new units or bonds. An ensemble of highly irregular capsids is generated, unlike the well-known icosahedrally symmetric viruses, but (we argue) plausible as a way to model the irregular capsids of retroviruses such as HIV. We found that (i) the probability of successful capsid completion decays exponentially with capsid size; (ii) capsid size depends strongly on spontaneous curvature and weakly on the ratio of the bending and stretching elastic stiffnesses of the shell; (iii) the degree of localization of Gaussian curvature (a measure of facetedness) depends heavily on the ratio of elastic stiffnesses.
△ Less
Submitted 8 August, 2006; v1 submitted 19 December, 2005;
originally announced December 2005.
-
Discussion on clusters, phasons, and quasicrystal stabilisation
Authors:
C. L. Henley,
M. de Boissieu,
W. Steurer
Abstract:
This paper summarises a two-hour discussion at the Ninth International Conference on Quasicrystals, including nearly 20 written comments sent afterwards, concerning (i) the meaning [if any] of clusters in quasicrystals; (ii) phason elasticity, and (iii) thermodynamic stabilisation of quasicrystals.
This paper summarises a two-hour discussion at the Ninth International Conference on Quasicrystals, including nearly 20 written comments sent afterwards, concerning (i) the meaning [if any] of clusters in quasicrystals; (ii) phason elasticity, and (iii) thermodynamic stabilisation of quasicrystals.
△ Less
Submitted 10 December, 2005;
originally announced December 2005.
-
Clusters, phason elasticity, and entropic stabilisation: a theory perspective
Authors:
C. L. Henley
Abstract:
Personal comments are made about the title subjects, including: the relation of Friedel oscillations to Hume-Rothery stabilisation; how calculations may resolve the random-tiling versus ideal pictures of quasicrystals; and the role of entropies apart from tile-configurational.
Personal comments are made about the title subjects, including: the relation of Friedel oscillations to Hume-Rothery stabilisation; how calculations may resolve the random-tiling versus ideal pictures of quasicrystals; and the role of entropies apart from tile-configurational.
△ Less
Submitted 10 December, 2005;
originally announced December 2005.
-
Co-rich decagonal Al-Co-Ni: predicting structure, orientational order, and puckering
Authors:
Nan Gu,
C. L. Henley,
M. Mihalkovic
Abstract:
We apply systematic methods previously used by Mihalkovic et al. to predict the structure of the `basic' Co-rich modification of the decagonal Al70 Co20 Ni10 layered quasicrystal, based on known lattice constants and previously calculated pair potentials. The modelling is based on Penrose tile decoration and uses Monte Carlo annealing to discover the dominant motifs, which are converted into rul…
▽ More
We apply systematic methods previously used by Mihalkovic et al. to predict the structure of the `basic' Co-rich modification of the decagonal Al70 Co20 Ni10 layered quasicrystal, based on known lattice constants and previously calculated pair potentials. The modelling is based on Penrose tile decoration and uses Monte Carlo annealing to discover the dominant motifs, which are converted into rules for another level of description. The result is a network of edge-sharing large decagons on a binary tiling of edge 10.5 A. A detailed analysis is given of the instability of a four-layer structure towards $c$-doubling and puckering of the atoms out of the layers, which is applied to explain the (pentagonal) orientational order.
△ Less
Submitted 10 December, 2005;
originally announced December 2005.
-
Effective Hamiltonian for the Pyrochlore antiferromagnet: semiclassical derivation and degeneracy
Authors:
U. Hizi,
C. L. Henley
Abstract:
In the classical pyrochlore lattice Heisenberg antiferromagnet, there is a macroscopic continuous ground state degeneracy. We study semiclassical limit of large spin length $S$, keeping only the lowest order (in 1/S) correction to the classical Hamiltonian. We perform a detailed analysis of the spin-wave modes, and using a real-space loop expansion, we produce an effective Hamiltonian, in which…
▽ More
In the classical pyrochlore lattice Heisenberg antiferromagnet, there is a macroscopic continuous ground state degeneracy. We study semiclassical limit of large spin length $S$, keeping only the lowest order (in 1/S) correction to the classical Hamiltonian. We perform a detailed analysis of the spin-wave modes, and using a real-space loop expansion, we produce an effective Hamiltonian, in which the degrees of freedom are Ising variables representing fluxes through loops in the lattice. We find a family of degenerate collinear ground states, related by gauge-like $Z_2$ transformations and provide bounds for the order of the degeneracy. We further show that the theory can readily be applied to determine the ground states of the Heisenberg Hamiltonian on related lattices, and to field-induced collinear magnetization plateau states.
△ Less
Submitted 16 December, 2005; v1 submitted 31 August, 2005;
originally announced September 2005.
-
Order by disorder and gauge-like degeneracy in quantum pyrochlore antiferromagnet
Authors:
Christopher L. Henley
Abstract:
The (three-dimensional) pyrochlore lattice antiferromagnet with Heisenberg spins of large spin length $S$ is a highly frustrated model with an macroscopic degeneracy of classical ground states. The zero-point energy of (harmonic order) spin wave fluctuations distinguishes a subset of these states. I derive an approximate but illuminating {\it effective Hamiltonian}, acting within the subspace of…
▽ More
The (three-dimensional) pyrochlore lattice antiferromagnet with Heisenberg spins of large spin length $S$ is a highly frustrated model with an macroscopic degeneracy of classical ground states. The zero-point energy of (harmonic order) spin wave fluctuations distinguishes a subset of these states. I derive an approximate but illuminating {\it effective Hamiltonian}, acting within the subspace of Ising spin configurations representing the {\it collinear} ground states. It consists of products of Ising spins around loops, i.e has the form of a $Z_2$ lattice gauge theory. The remaining ground state entropy is still infinite but not extensive, being $O(L)$ for system size $O(L^3)$. All these ground states have unit cells bigger than those considered previously.
△ Less
Submitted 31 August, 2005;
originally announced September 2005.
-
Many-Body Density Matrices On a Two-Dimensional Square Lattice: Noninteracting and Strongly Interacting Spinless Fermions
Authors:
Siew-Ann Cheong,
Christopher L. Henley
Abstract:
The reduced density matrix of an interacting system can be used as the basis for a truncation scheme, or in an unbiased method to discover the strongest kind of correlation in the ground state. In this paper, we investigate the structure of the many-body fermion density matrix of a small cluster in a square lattice. The cluster density matrix is evaluated numerically over a set of finite systems…
▽ More
The reduced density matrix of an interacting system can be used as the basis for a truncation scheme, or in an unbiased method to discover the strongest kind of correlation in the ground state. In this paper, we investigate the structure of the many-body fermion density matrix of a small cluster in a square lattice. The cluster density matrix is evaluated numerically over a set of finite systems, subject to non-square periodic boundary conditions given by the lattice vectors $\bR_1 \equiv (R_{1x}, R_{1y})$ and $\bR_2 \equiv (R_{2x}, R_{2y})$. We then approximate the infinite-system cluster density-matrix spectrum, by averaging the finite-system cluster density matrix (i) over degeneracies in the ground state, and orientations of the system relative to the cluster, to ensure it has the proper point-group symmetry; and (ii) over various twist boundary conditions to reduce finite size effects. We then compare the eigenvalue structure of the averaged cluster density matrix for noninteracting and strongly-interacting spinless fermions, as a function of the filling fraction $\nbar$, and discuss whether it can be approximated as being built up from a truncated set of single-particle operators.
△ Less
Submitted 22 July, 2006; v1 submitted 31 August, 2005;
originally announced August 2005.
-
Semiclassical ordering in the large-N pyrochlore antiferromagnet
Authors:
U. Hizi,
Prashant Sharma,
C. L. Henley
Abstract:
We study the semiclassical limit of the $Sp(N)$ generalization of the pyrochlore lattice Heisenberg antiferromagnet by expanding about the $N \to \infty$ saddlepoint in powers of a generalized inverse spin. To leading order, we write down an effective Hamiltonian as a series in loops on the lattice. Using this as a formula for calculating the energy of any classical ground state, we perform Mont…
▽ More
We study the semiclassical limit of the $Sp(N)$ generalization of the pyrochlore lattice Heisenberg antiferromagnet by expanding about the $N \to \infty$ saddlepoint in powers of a generalized inverse spin. To leading order, we write down an effective Hamiltonian as a series in loops on the lattice. Using this as a formula for calculating the energy of any classical ground state, we perform Monte-Carlo simulations and find a unique collinear ground state. This state is not a ground state of linear spin-wave theory, and can therefore not be a physical (N=1) semiclassical ground state.
△ Less
Submitted 12 October, 2005; v1 submitted 22 June, 2005;
originally announced June 2005.
-
Power-law spin correlations in pyrochlore antiferromagnets
Authors:
C. L. Henley
Abstract:
The ground state ensemble of the highly frustrated pyrochlore-lattice antiferromagnet can be mapped to a coarse-grained ``polarization'' field satisfying a zero-divergence condition From this it follows that the correlations of this field, as well as the actual spin correlations, decay with separation like a dipole-dipole interaction ($1/|R|^3$). Furthermore, a lattice version of the derivation…
▽ More
The ground state ensemble of the highly frustrated pyrochlore-lattice antiferromagnet can be mapped to a coarse-grained ``polarization'' field satisfying a zero-divergence condition From this it follows that the correlations of this field, as well as the actual spin correlations, decay with separation like a dipole-dipole interaction ($1/|R|^3$). Furthermore, a lattice version of the derivation gives an approximate formula for spin correlations, with several features that agree well with simulations and neutron-diffraction measurements of diffuse scattering, in particular the pinch-point (pseudo-dipolar) singularities at reciprocal lattice vectors. This system is compared to others in which constraints also imply diffraction singularities, and other possible applications of the coarse-grained polarization are discussed.
△ Less
Submitted 30 June, 2004;
originally announced July 2004.
-
Combined energy -- diffraction data refinement of decagonal AlNiCo
Authors:
M. Mihalkovič,
C. L. Henley,
M. Widom
Abstract:
We incorporate realistic pair potential energies directly into a non-linear least-square fit of diffraction data to quantitatively compare structure models with experiment for the Ni-rich $d$(AlNiCo) quasicrystal. The initial structure models are derived from a few {\it a priori} assumptions (gross features of the Patterson function) and the pair potentials. In place of the common hyperspace app…
▽ More
We incorporate realistic pair potential energies directly into a non-linear least-square fit of diffraction data to quantitatively compare structure models with experiment for the Ni-rich $d$(AlNiCo) quasicrystal. The initial structure models are derived from a few {\it a priori} assumptions (gross features of the Patterson function) and the pair potentials. In place of the common hyperspace approach to the structure refinement of quasicrystals, we use a real-space tile decoration scheme, which does not rely on strict quasiperiodicity, and makes it easy to enforce sensible local arrangements of the atoms. Inclusion of the energies provides information complementary to the diffraction data and protects the fit procedure from converging on spurious solutions. The method pinpoints sites which are likely to break the symmetry of their local environment.
△ Less
Submitted 26 November, 2003;
originally announced November 2003.