Source author record

Shan-Ho Tsai

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

Catalog footprint

Research overview

11works
4topics
4close collaborators

Published work

11 publications

preprint2015arXiv

Exploring Replica-Exchange Wang-Landau sampling in higher-dimensional parameter space

We considered a higher-dimensional extension for the replica-exchange Wang-Landau algorithm to perform a random walk in the energy and magnetization space of the two-dimensional Ising model. This hybrid scheme combines the advantages of Wang-Landau and Replica-Exchange algorithms, and the one-dimensional version of this approach has been shown to be very efficient and to scale well, up to several thousands of computing cores. This approach allows us to split the parameter space of the system to be simulated into several pieces and still perform a random walk over the entire parameter range, ensuring the ergodicity of the simulation. Previous work, in which a similar scheme of parallel simulation was implemented without using replica exchange and with a different way to combine the result from the pieces, led to discontinuities in the final density of states over the entire range of parameters. From our simulations, it appears that the replica-exchange Wang-Landau algorithm is able to overcome this difficulty, allowing exploration of higher parameter phase space by keeping track of the joint density of states.

preprint2004arXiv

Phase transitions in ferro-antiferromagnetic bilayers with a stepped interface

We have studied magnetic ordering in ferro/antiferromagnetic (F/AF) bilayers using Monte Carlo simulations of classical Heisenberg spins. For both flat and stepped interfaces we observed order in the AF above the Neel temperature, with the AF spins aligning collinearly with the F moments. In the case of the stepped interface there is a transition from collinear to perpendicular alignment of the F and AF spins at a lower temperature.

preprint2003arXiv

Dynamic critical behavior of the classical anisotropic BCC Heisenberg antiferromagnet

Using a recently implemented integration method [Krech et. al.] based on an iterative second-order Suzuki-Trotter decomposition scheme, we have performed spin dynamics simulations to study the critical dynamics of the BCC Heisenberg antiferromagnet with uniaxial anisotropy. This technique allowed us to probe the narrow asymptotic critical region of the model and estimate the dynamic critical exponent z=2.25±0.08z=2.25 \pm 0.08. Comparisons with competing theories and experimental results are presented.

preprint1999arXiv

Improved Spin Dynamics Simulations of Magnetic Excitations

Using Suzuki-Trotter decompositions of exponential operators we describe new algorithms for the numerical integration of the equations of motion for classical spin systems. These techniques conserve spin length exactly and, in special cases, also conserve the energy and maintain time reversibility. We investigate integration schemes of up to eighth order and show that these new algorithms can be used with much larger time steps than a well established predictor-corrector method. These methods may lead to a substantial speedup of spin dynamics simulations, however, the choice of which order method to use is not always straightforward.

preprint1999arXiv

Spin dynamics simulations of excitations and critical dynamics in RbMnF_3

Spin-dynamics simulations have been used to investigate the dynamic behavior of RbMnF_3, treating it as a classical Heisenberg antiferromagnet on a simple cubic lattice. Time-evolutions of spin configurations were determined numerically from coupled equations of motion for individual spins using a new algorithm which is based on Suzuki-Trotter decompositions of exponential operators. The dynamic structure factor was calculated from the space- and time-displaced spin-spin correlation function. The crossover from hydrodynamic to critical behavior of the dispersion curve and spin-wave half-width was studied as the temperature was increased towards the critical value. The dynamic critical exponent was estimated to be z=(1.43\pm 0.03), which is slightly lower than the dynamic scaling prediction, but in good agreement with a recent experimental value. Comparisons are made of both the dispersion curve and the lineshapes obtained from our simulations with very recent experimental results for RbMnF_3 are presented.

preprint1998arXiv

Chromatic Polynomials for J(∏H)IJ(\prod H)I Strip Graphs and their Asymptotic Limits

We calculate the chromatic polynomials PP for nn-vertex strip graphs of the form J(∏ℓ=1mH)IJ(\prod_{\ell=1}^m H)I, where JJ and II are various subgraphs on the left and right ends of the strip, whose bulk is comprised of mm-fold repetitions of a subgraph HH. The strips have free boundary conditions in the longitudinal direction and free or periodic boundary conditions in the transverse direction. This extends our earlier calculations for strip graphs of the form (∏ℓ=1mH)I(\prod_{\ell=1}^m H)I. We use a generating function method. From these results we compute the asymptotic limiting function W=lim⁡n→∞P1/nW=\lim_{n \to \infty}P^{1/n}; for q∈Z+q \in {\mathbb Z}_+ this has physical significance as the ground state degeneracy per site (exponent of the ground state entropy) of the qq-state Potts antiferromagnet on the given strip. In the complex qq plane, WW is an analytic function except on a certain continuous locus B{\cal B}. In contrast to the (∏ℓ=1mH)I(\prod_{\ell=1}^m H)I strip graphs, where B{\cal B} (i) is independent of II, and (ii) consists of arcs and possible line segments that do not enclose any regions in the qq plane, we find that for some J(∏ℓ=1mH)IJ(\prod_{\ell=1}^m H)I strip graphs, B{\cal B} (i) does depend on II and JJ, and (ii) can enclose regions in the qq plane. Our study elucidates the effects of different end subgraphs II and JJ and of boundary conditions on the infinite-length limit of the strip graphs.

preprint1998arXiv

End Graph Effects on Chromatic Polynomials for Strip Graphs of Lattices and their Asymptotic Limits

We report exact calculations of the ground state degeneracy per site (exponent of the ground state entropy) W({G},q)W(\{G\},q) of the qq-state Potts antiferromagnet on infinitely long strips with specified end graphs for free boundary conditions in the longitudinal direction and free and periodic boundary conditions in the transverse direction. This is equivalent to calculating the chromatic polynomials and their asymptotic limits for these graphs. Making the generalization from q∈Z+q \in {\mathbb Z}_+ to q∈Cq \in {\mathbb C}, we determine the full locus B{\cal B} on which WW is nonanalytic in the complex qq plane. We report the first example for this class of strip graphs in which B{\cal B} encloses regions even for planar end graphs. The bulk of the specific strip graph that exhibits this property is a part of the (33⋅42)(3^3 \cdot 4^2) Archimedean lattice.

preprint1998arXiv

Fourth-Order Cumulants to Characterize the Phase Transitions of a Spin-1 Ising Model

Fourth-order cumulants of physical quantities have been used to characterize the nature of a phase transition. In this paper we report some Monte Carlo simulations to illustrate the behavior of fourth-order cumulants of magnetization and energy across second and first-order transitions in the phase diagram of a well known spin-1 Ising model. Simple ideas from the theory of thermodynamic fluctuations are used to account for the behavior of these cumulants.

preprint1998arXiv

Ground State Entropy of Potts Antiferromagnets on Homeomorphic Families of Strip Graphs

We present exact calculations of the zero-temperature partition function, and ground-state degeneracy (per site), WW, for the qq-state Potts antiferromagnet on a variety of homeomorphic families of planar strip graphs G=(Ch)k1,k2,Σ,k,mG = (Ch)_{k_1,k_2,Σ,k,m}, where k1k_1, k2k_2, ΣΣ, and kk describe the homeomorphic structure, and mm denotes the length of the strip. Several different ways of taking the total number of vertices to infinity, by sending (i) m→∞m \to \infty with k1k_1, k2k_2, and kk fixed; (ii) k1k_1 and/or k2→∞k_2 \to \infty with mm, and kk fixed; and (iii) k→∞k \to \infty with mm and p=k1+k2p=k_1+k_2 fixed are studied and the respective loci of points B{\cal B} where WW is nonanalytic in the complex qq plane are determined. The B{\cal B}'s for limit (i) are comprised of arcs which do not enclose regions in the qq plane and, for many values of pp and kk, include support for Re(q)<0Re(q) < 0. The B{\cal B} for limits (ii) and (iii) is the unit circle ∣q−1∣=1|q-1|=1.

preprint1998arXiv

Ground State Entropy of Potts Antiferromagnets: Homeomorphic Classes with Noncompact W Boundaries

We present exact calculations of the zero-temperature partition function Z(G,q,T=0)Z(G,q,T=0) and ground-state degeneracy W({G},q)W(\{G\},q) for the qq-state Potts antiferromagnet on a number of families of graphs GG for which (generalizing qq from Z+{\mathbb Z}_+ to C{\mathbb C}) the boundary B{\cal B} of regions of analyticity of WW in the complex qq plane is noncompact, passing through z=1/q=0z=1/q=0. For these types of graphs, since the reduced function Wred.=q−1WW_{red.}=q^{-1}W is nonanalytic at z=0z=0, there is no large--qq Taylor series expansion of Wred.W_{red.}. The study of these graphs thus gives insight into the conditions for the validity of the large--qq expansions. It is shown how such (families of) graphs can be generated from known families by homeomorphic expansion.

preprint1997arXiv

Chromatic Polynomials for Families of Strip Graphs and their Asymptotic Limits

We calculate the chromatic polynomials P((Gs)m,q)P((G_s)_m,q) and, from these, the asymptotic limiting functions W({Gs},q)=lim⁡n→∞P(Gs,q)1/nW(\{G_s\},q)=\lim_{n \to \infty}P(G_s,q)^{1/n} for families of nn-vertex graphs (Gs)m(G_s)_m comprised of mm repeated subgraphs HH adjoined to an initial graph II. These calculations of W({Gs},q)W(\{G_s\},q) for infinitely long strips of varying widths yield important insights into properties of W(Λ,q)W(Λ,q) for two-dimensional lattices ΛΛ. In turn, these results connect with statistical mechanics, since W(Λ,q)W(Λ,q) is the ground state degeneracy of the qq-state Potts model on the lattice ΛΛ. For our calculations, we develop and use a generating function method, which enables us to determine both the chromatic polynomials of finite strip graphs and the resultant W({Gs},q)W(\{G_s\},q) function in the limit n→∞n \to \infty. From this, we obtain the exact continuous locus of points B{\cal B} where W({Gs},q)W(\{G_s\},q) is nonanalytic in the complex qq plane. This locus is shown to consist of arcs which do not separate the qq plane into disconnected regions. Zeros of chromatic polynomials are computed for finite strips and compared with the exact locus of singularities B{\cal B}. We find that as the width of the infinitely long strips is increased, the arcs comprising B{\cal B} elongate and move toward each other, which enables one to understand the origin of closed regions that result for the (infinite) 2D lattice.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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