Evaluation of the continuum limit of the axial anomaly and index is sketched for the staggered overlap Dirac operator. There are new complications compared to the usual overlap case due to the distribution of the spin and flavor components around lat
tice hypercubes in the staggered formalism. The index is found to correctly reproduce the continuum index, but for the axial anomaly this is only true after averaging over the sites of a lattice hypercube.
Results on the computational efficiency of 2-flavor staggered Wilson fermions compared to usual Wilson fermions in a quenched lattice QCD simulation on $16^3times32$ lattice at $beta=6$ are reported. We compare the cost of inverting the Dirac matrix
on a source by the conjugate gradient (CG) method for both of these fermion formulations, at the same pion masses, and without preconditioning. We find that the number of CG iterations required for convergence, averaged over the ensemble, is less by a factor of almost 2 for staggered Wilson fermions, with only a mild dependence on the pion mass. We also compute the condition number of the fermion matrix and find that it is less by a factor of 4 for staggered Wilson fermions. The cost per CG iteration, dominated by the cost of matrix-vector multiplication for the Dirac matrix, is known from previous work to be less by a factor 2-3 for staggered Wilson compared to usual Wilson fermions. Thus we conclude that staggered Wilson fermions are 4-6 times cheaper for inverting the Dirac matrix on a source in the quenched backgrounds of our study.
A new formulation of chiral fermions on the lattice is presented. It is a version of overlap fermions, but built from the computationally efficient staggered fermions rather than the previously used Wilson fermions. The construction reduces the four
quark flavors described by the staggered fermion to two quark flavors; this pair can be taken as the up and down quarks in Lattice QCD. The exact flavored chiral symmetry of the staggered fermion gets converted into an unflavored Ginsparg-Wilson chiral symmetry of the new overlap fermion, which also has pairs of exact chiral zero-modes satisfying the Index Theorem. Stability under radiative corrections is checked. A domain wall formulation giving a truncation of this overlap construction is also outlined.
We introduce a new Forward-Flux Sampling in Time (FFST) algorithm to efficiently measure transition times in rare-event processes in non-equilibrium systems, and apply it to study the first-order (discontinuous) kinetic transition in the Ziff-Gulari-
Barshad model of catalytic surface reaction. The average time for the transition to take place, as well as both the spinodal and transition points, are clearly found by this method.
In many dynamical systems there is a large separation of time scales between typical events and rare events which can be the cases of interest. Rare-event rates are quite difficult to compute numerically, but they are of considerable practical import
ance in many fields: for example transition times in chemical physics and extinction times in epidemiology can be very long, but are quite important. We present a very fast numerical technique that can be used to find long transition times (very small rates) in low-dimensional systems, even if they lack detailed balance. We illustrate the method for a bistable non-equilibrium system introduced by Maier and Stein and a two-dimensional (in parameter space) epidemiology model.
A way to identify the would-be zero-modes of staggered lattice fermions away from the continuum limit is presented. Our approach also identifies the chiralities of these modes, and their index is seen to be determined by gauge field topology in accor
dance with the Index Theorem. The key idea is to consider the spectral flow of a certain hermitian version of the staggered Dirac operator. The staggered fermion index thus obtained can be used as a new way to assign the topological charge of lattice gauge fields. In a numerical study in U(1) backgrounds in 2 dimensions it is found to perform as well as the Wilson index while being computationally more efficient. It can also be expressed as the index of an overlap Dirac operator with a new staggered fermion kernel.
We present a technique, which we call etching, which we use to study the harmonic measure of Fortuin-Kasteleyn clusters in the Q-state Potts model for Q=1-4. The harmonic measure is the probability distribution of random walkers diffusing onto the pe
rimeter of a cluster. We use etching to study regions of clusters which are extremely unlikely to be hit by random walkers, having hitting probabilities down to 10^(-4600). We find good agreement between the theoretical predictions of Duplantier and our numerical results for the generalized dimension D(q), including regions of small and negative q.
We obtain the harmonic measure of diffusion-limited aggregate (DLA) clusters using a biased random-walk sampling technique which allows us to measure probabilities of random walkers hitting sections of clusters with unprecedented accuracy; our result
s include probabilities as small as 10^(-80). We find the multifractal D(q) spectrum including regions of small and negative q. Our algorithm allows us to obtain the harmonic measure for clusters more than an order of magnitude larger than those achieved using the method of iterative conformal maps, which is the previous best method. We find a phase transition in the singularity spectrum f(alpha) at alpha approximately equal to 14 and also find a minimum q of D(q), q_{min} = 0.9 plus or minus 0.05.
The totally asymmetric simple exclusion process (TASEP) is a well studied example of far-from-equilibrium dynamics. Here, we consider a TASEP with open boundaries but impose a global constraint on the total number of particles. In other words, the bo
undary reservoirs and the system must share a finite supply of particles. Using simulations and analytic arguments, we obtain the average particle density and current of the system, as a function of the boundary rates and the total number of particles. Our findings are relevant to biological transport problems if the availability of molecular motors becomes a rate-limiting factor.
As a solvable and broadly applicable model system, the totally asymmetric exclusion process enjoys iconic status in the theory of non-equilibrium phase transitions. Here, we focus on the time dependence of the total number of particles on a 1-dimensi
onal open lattice, and its power spectrum. Using both Monte Carlo simulations and analytic methods, we explore its behavior in different characteristic regimes. In the maximal current phase and on the coexistence line (between high/low density phases), the power spectrum displays algebraic decay, with exponents -1.62 and -2.00, respectively. Deep within the high/low density phases, we find pronounced emph{oscillations}, which damp into power laws. This behavior can be understood in terms of driven biased diffusion with conserved noise in the bulk.
