Researcher profile

Keola Wierschem

Keola Wierschem contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
15works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2016arXiv

Detection of symmetry-protected topological order in AKLT states by exact evaluation of the strange correlator

The strange correlator [Phys. Rev. Lett. 112, 247202 (2014)] has been proposed as a measure of symmetry protected topological order in one- and two-dimensional systems. It takes the form of a spin-spin correlation function, computed as a mixed overlap between the state of interest and a trivial local product state. We demonstrate that it can be computed exactly (asymptotically, in the Monte Carlo sense) for various Affleck-Kennedy-Lieb-Tasaki states by direct evaluation of the wave function within the valence bond loop gas framework. We present results for lattices with chain, square, honeycomb, cube, diamond, and hyperhoneycomb geometries. In each case, the spin quantum number S is varied such that 2S (the number of valence bonds emerging from each site) achieves various integer multiples of the lattice coordination number. We introduce the concept of strange correlator loop winding number and point to its utility in testing for the presence of symmetry protected topological order.

preprint2015arXiv

Characterizing the Haldane phase in quasi-one-dimensional spin-1 Heisenberg antiferromagnets

We review the basic properties of the Haldane phase in spin-1 Heisenberg antiferromagnetic chains, including its persistence in quasi-one-dimensional geometries. Using large-scale numerical simulations, we map out the phase diagram for a realistic model applicable to experimental Haldane compounds. We also investigate the effect of different chain coupling geometries and confirm a general mean field universality of the critical coupling times the coordination number of the lattice. Inspired by the recent development of characterization of symmetry protected topological states, of which the Haldane phase of spin-1 Heisenberg antiferromagnetic chain is a preeminent example, we provide direct evidence that the quasi-one-dimensional Haldane phase is indeed a non-trivial symmetry protected topological state.

preprint2015arXiv

Origin of modulated phases and magnetic hysteresis in TmB_4

We investigate the low temperature magnetic phases in TmB_4, a metallic quantum magnet on the geometrically frustrated Shastry-Sutherland lattice, using co-ordinated experimental and theoretical studies. Our results provide an explanation for the appearance of the intriguing fractional plateau in TmB_4 and accompanying magnetic hysteresis. Together with observation of the bump in the half-plateau, our results support the picture that the magnetization plateau structure in TmB_4 is strongly influenced by the zero-field modulated phases. We present a phenomenological model to explain the appearance of the modulated phases and a microscopic Hamiltonian that captures the complete magnetic behavior of TmB_4.

preprint2014arXiv

Magnetic phases in the S=1 Shastry-Sutherland model with uniaxial anisotropy

We explore the field induced magnetic phases of an $S=1$ $XXZ$ model with single-ion anisotropy and large Ising-like anisotropy on a Shastry Sutherland lattice over a wide range of Hamiltonian parameters and applied magnetic field. The multitude of ground state phases are characterized in detail in terms of their thermodynamic properties and the underlying classical (Ising limit) spin arrangements for the plateau phases are identified by calculating the static structure factors. The enlarged local Hilbert space of the $S=1$ spins results in several new ground state phases that are not realized for $S=1/2$ spins. These include the quantum paramagnetic state that is ubiquitous to $S=1$ spins with single ion anisotropy, two different spin supersolid phases (with distinct longitudinal ordering) and a magnetization plateau that arises as a direct descendant of the 1/3 plateau due to quantum fluctuations that are not possible for $S=1/2$ spins. We predict the same mechanism will lead to plateaus at smaller fractions of 1/3 for higher spins. The full momentum dependence of the longitudinal and transverse components of the static structure factor is calculated in the spin supersolid phase to demonstrate the simultaneous existence of diagonal and off-diagonal long-range order as well as the different longitudinal orderings.

preprint2014arXiv

Quenching the Haldane gap in spin-1 Heisenberg antiferromagnets

We consider a quasi-one-dimensional system of spin-1 Heisenberg antiferromagnetic chains in 2D and 3D hypercubic lattices with interchain coupling $J$ and uniaxial single-ion anisotropy $D$. Using large scale numerical simulations, we map out the $J-D$ phase diagram and investigate the low lying excitations of the Haldane phase in the $J\ll1$ limit. We also provide direct evidence that the Haldane phase remains a non-trivial symmetry protected topological state for small but finite $J$.

preprint2014arXiv

Strange correlations in spin-1 Heisenberg antiferromagnets

We study the behavior of the recently proposed "strange correlator" [Phys. Rev. Lett. {\bf 112}, 247202 (2014)] in spin-1 Heisenberg antiferromagnetic chains with uniaxial single-ion anisotropy. Using projective quantum Monte Carlo, we are able to directly access the strange correlator in a variety of phases, as well as to examine its critical behavior at the quantum phase transition between trivial and non-trivial symmetry protected topological phases. After finding the expected long-range behavior in these two symmetry conserving phases, we go on to verify the topological nature of two-leg and three-leg spin-1 Heisenberg antiferromagnetic ladders. This demonstrates the power of the strange correlator in distinguishing between trivial and non-trivial symmetry protected topological phases.

preprint2013arXiv

Columnar antiferromagnetic order and spin supersolid phase on the extended Shastry-Sutherland lattice

We use large scale quantum Monte Carlo simulations to study an extended version of the canonical Shastry-Sutherland model -- including additional interactions and exchange anisotropy -- over a wide range of interaction parameters and an applied magnetic field. The model is appropriate for describing the low energy properties of some members of the rare earth tetraborides. Working in the limit of large Ising-like exchange anisotropy, we demonstrate the stabilization of columnar antiferromagnetic order in the ground state at zero field and an extended magnetization plateau at 1/2 the saturation magnetization in the presence of an applied longitudinal magnetic field -- qualitatively similar to experimentally observed low-temperature phases in ErB$_4$. Our results show that for an optimal range of exchange parameters, a spin supersolid ground state is realized over a finite range of applied field between the columnar antiferromagnetic phase and the magnetization plateau. The full momentum dependence of the longitudinal and transverse components of the static structure factor is calculated in the spin supersolid phase to demonstrate the simultaneous existence of diagonal and off-diagonal long-range order. Our results will provide crucial guidance in designing further experiments to search for the interesting spin supersolid phase in ErB$_4$.

preprint2013arXiv

Dimensional Crossover in Quasi-one-dimensional Spin-1 Heisenberg Antiferromagnets

We study the quenching of the Haldane gap in quasi-one-dimensional systems of weakly coupled spin-1 antiferromagnetic Heisenberg chains. The critical interchain coupling Jc required to stabilize long range magnetic order can be accurately determined from large scale quantum Monte Carlo calculations. Several different geometries of coupled chains are studied, illustrating the dependence of Jc on the coordination of chains. For bipartite geometries, ferromagnetically coupled chains yield similar magnitudes for Jc.

preprint2013arXiv

Magnetization plateaus in generalized Shastry-Sutherland models

We study an anisotropic Heisenberg antiferromagnet with ferromagnetic transverse spin exchange using exact quantum Monte Carlo methods. Such a model is relevant to a class of rare earth tetraboride materials that display a range of magnetization plateaus under applied magnetic field. The layered arrangement of magnetic ions in these materials is topologically equivalent to the Shastry-Sutherland lattice. In this frustrated geometry, we study the interplay of next-nearest neighbor interactions in stabilizing a plateau at half the saturation magnetization (or 1/2 plateau). We also show hysteresis-like behavior at the onset of the 1/3 plateau.

preprint2013arXiv

Neel to spin-Peierls transition in a quasi-1D Heisenberg model coupled to bond phonons

The zero and finite temperature spin-Peierls transitions in a quasi-one-dimensional spin-1/2 Heisenberg model coupled to adiabatic bond phonons is investigated using the Stochastic Series Expansion (SSE) Quantum Monte Carlo (QMC) method. The quantum phase transition from a gapless Neel state to a spin-gapped Peierls state is studied in the parameter space spanned by spatial anisotropy, inter-chain coupling strength and spin-lattice coupling strength. It is found that for any finite inter-chain coupling, the transition to a dimerized Peierls ground state only occurs when the spin-lattice coupling exceeds a finite, non-zero critical value. This is in contrast to the pure 1D model (zero inter-chain coupling), where adiabatic/classical phonons lead to a dimerized ground state for any non-zero spin-phonon interaction. The phase diagram in the parameter space shows that for a strong inter-chain coupling, the relation between the inter-chain coupling and the critical value of the spin-phonon interaction is linear whereas for weak inter-chain coupling, this behavior is found to have a natural logarithm-like relation. No region was found to have a long range magnetic order and dimerization occurring simultaneously. Instead, the Neel state order vanishes simultaneously with the setting in of the spin-Peierls state. For the thermal phase transition, a continuous heat capacity with a peak at the critical temperature, $T_{c}$, shows a second order phase transition. The variation of the equilibrium bond length distortion, $δ_{eq}$, with temperature showed a power law relation which decayed to zero as the temperature was increased to $T_{c}$, indicating a continuous transition from the dimerized phase to a paramagnetic phase with uniform bond length and zero antiferromagnetic susceptibility.

preprint2013arXiv

Phase Diagram and Magnetic Excitations of Anisotropic Spin-One Magnets

We use a generalized spin wave approach and large scale quantum Monte Carlo (QMC) simulations to study the quantum phase diagram and quasiparticle excitations of the S=1 Heisenberg model with an easy-plane single-ion anisotropy in dimensions d=2 and 3. We consider two alternative approximations for describing the quantum paramagnetic state: the standard Holstein-Primakoff approximation and a modified treatment in which the local constraint (finite dimension of the local Hilbert space) is enforced by introducing a Lagrange multiplier. While both approximations produce qualitatively similar results, the latter approach is the only one that is in good quantitative agreement with the quantum phase diagram and the quasiparticle dispersions obtained with QMC. This result is very important for low-temperature studies of quantum paramagnets in magnetic fields because it shows that a simple modification of the standard analytical approach should produce much better quantitative agreement between theory and experiment.

preprint2012arXiv

Magnetic and nematic orderings in spin-1 antiferromagnets with single-ion anisotropy

We study a spin-one Heisenberg model with exchange interaction, J, uniaxial single-ion exchange anisotropy, D, and Zeeman coupling to a magnetic field, B, parallel to the symmetry axis. We com- pute the (D/J,B/J) quantum phase diagram for square and simple cubic lattices by combining an- alytical and Quantum Monte Carlo approaches, and find a transition between XY-antiferromagnetic and ferronematic phases that spontaneously break the U(1) symmetry of the model. In the language of bosonic gases, this is a transition between a Bose-Einstein condensate (BEC) of single bosons and a BEC of pairs. Our work opens up new avenues for measuring this transition in real magnets.

preprint2011arXiv

Dimensional crossover in spin-1 Heisenberg antiferromagnets: a quantum Monte Carlo study

We present results of large scale simulations of the spin-1 Heisenberg antiferromagnet on a tetragonal lattice. The stochastic series expansion quantum Monte Carlo method is used to calculate equilibrium thermodynamic variables in the presence of an external magnetic field. In particular, the low temperature magnetization curve is investigated in the quasi-one-dimensional (Q1D), quasi-two-dimensional (Q2D), and three-dimensional (3D) limits. Starting from the 3D limit, the Q1D (Q2D) limit is achieved by reducing the in-plane (out-of-plane) spin coupling strength towards zero. In the Q1D limit, a Haldane gap appears in the magnetization curve at low magnetic field. Additionally, near the saturation field the slope of the magnetization curve increases substantially, approaching the infinite-slope behavior of a one-dimensional spin-1 chain. A similar (though less pronounced) effect is seen in the Q2D limit. We also study the effect of uniaxial single-ion anisotropy on the magnetization curve for Q1D and Q2D systems. Our results will be important in understanding the field-induced behavior of a class of low-dimensional Ni-based quantum magnets.

preprint2011arXiv

Simulation of melting of two dimensional Lennard-Jones solids

We study the nature of melting of a two dimensional (2D) Lennard-Jones solid using large scale Monte Carlo simulation. We use systems of up to 102,400 particles to capture the decay of the correlation functions associated with translational order (TO) as well as the bond-orientational (BO) order. We study the role of dislocations and disclinations and their distribution functions. We computed the temperature dependence of the second moment of the two order parameters, i.e., of the TO order parameter as well as of the order parameter associated with BO order and using finite-size scaling we determined the corresponding two anomalous dimension critical exponents eta. We also computed the temperature dependent distribution of these two order parameters on the complex plane which support a two stage melting with a hexatic phase as an intermediate phase. The analysis of our results leads to a consistent picture strongly supporting a two stage melting scenario as predicted by the Kosterlitz, Thouless, Halperin, Nelson, and Young (KTHNY) theory. We find that the two correlation lengths associated with translational and bond-orientational order have a distinctly different temperature dependence each diverging at different temperature and that their finite size scaling properties are consistent with the KTHNY theory. We also used the temperature dependence of etas and their theoretical bounds to provide estimates for the two critical temperatures, which are found to be in agreement with those determined from the divergent correlation lengths and from the Binder ratio.

preprint2010arXiv

Quantum simulation of 3He impurities and of 4He interstitials in solid 4He

We have studied the role of an atomic 3He impurity and an interstitial 4He atom in two- and three-dimensional solid 4He using path integral Monte Carlo (PIMC) simulation. We find that when a substitutional 3He impurity is introduced, the impurity becomes localized and occupies an ideal lattice site. When an interstitial 3He impurity is introduced in the 4He solid, we find that the impurity becomes localized at a substitutional position and, thus, promotes the extra 4He atom to the interstitial space. As a consequence we find that the one-body density matrix (OBDM) and the superfluid fraction, for the case of a 4He solid with an interstitial impurity, are very similar to those calculated for a 4He solid with a 4He interstitial atom. Namely, while the off-diagonal OBDM approaches zero exponentially with increasing particle displacement for the "pure" solid, an interstitial 4He atom or a 3He impurity appear to enhance it at long distances. Finally, the effective mass of the 3He impurity quasiparticle in 2D and 3D crystalline 4He is estimated.