PubMed Health⌕ Search

Biomedical subjects

Anders W Sandvik

Publications and source records attributed to Anders W Sandvik.

10 recordsLinked to original sources

Low-energy dynamics of the two-dimensional S=1/2 Heisenberg antiferromagnet on percolating clusters.

We investigate the quantum dynamics of site diluted S=1/2 Heisenberg antiferromagnetic clusters at the 2D percolation threshold. We use Lanczos diagonalization to calculate the lowest excitation gap Delta and, to reach larger sizes, use quantum Monte Carlo simulations to study an upper bound for Delta obtained from sum rules involving the staggered structure factor and susceptibility. Scaling the gap distribution with the cluster length L, Delta approximately L(-), we obtain a dynamic exponent z approximately 2D(f), where D(f)=91/48 is the fractal dimensionality of the percolating cluster. This is in contrast with previous expectations of z=D(f). We argue that the low-energy excitations are due to weakly coupled effective moments formed due to local imbalance in sublattice occupation.

Journal Article↗

Quantum criticality and percolation in dimer-diluted two-dimensional antiferromagnets.

The S = 1/2 Heisenberg model is considered on bilayer and single-layer square lattices with couplings J1, J2, with each spin belonging to one J2-coupled dimer. A transition from a Néel to disordered ground state occurs at a critical value of g = J2/J1. The systems are here studied at their dimer-dilution percolation points p*. The multicritical point (g*,p*) previously found for the bilayer is not reproduced for the single layer. Instead, there is a line of critical points (g < g*, p*) with continuously varying exponents. The uniform magnetic susceptibility diverges as T(-alpha) with alpha element of [1/2,1]. This unusual behavior is attributed to an effective free-moment density approximately T(1-alpha). The susceptibility of the bilayer is not divergent but exhibits remarkably robust quantum-critical scaling.

Journal Article↗

Ground state projection of quantum spin systems in the valence-bond basis.

A Monte Carlo method for quantum spin systems is formulated in the valence-bond basis. The nonorthogonality allows for an efficient importance-sampled projection of the ground state out of an arbitrary state. The method provides access to resonating valence-bond physics, enables a direct estimator for the singlet-triplet gap, and extends the class of models that can be studied without negative-sign problems. As a demonstration, the valence-bond distribution in the ground state of the 2D Heisenberg antiferromagnet is calculated. Generalizations of the method to fermion systems are also discussed.

Journal Article↗

Stochastic series expansion algorithm for the S = 1/2 XY model with four-site ring exchange.

We describe a stochastic series expansion quantum Monte Carlo method for a two-dimensional S = 1/2 XY model (or, equivalently, hard-core bosons at half filling) which in addition to the standard pair interaction J includes a four-site term K that flips spins on a square plaquette. The model has three ordered ground state phases; for K/J approximately < or = 8 it has long-range xy spin order (superfluid bosons), for K/J approximately > or = 15 it has staggered spin order in the z direction (charge-density wave), and between these phases it is in a state with columnar order in the bond and plaquette energy densities. We discuss an implementation of directed-loop updates for the SSE simulations of this model and also introduce a "multibranch" cluster update which significantly reduces the autocorrelation times for large K/J. In addition to the pure J-K model, which in the z basis has only off-diagonal terms, we also discuss modifications of the algorithm needed when various diagonal interactions are included.

Journal Article↗

Ground state phases of the half-filled one-dimensional extended hubbard model.

Using quantum Monte Carlo simulations, results of a strong-coupling expansion, and Luttinger liquid theory, we determine quantitatively the ground state phase diagram of the one-dimensional extended Hubbard model with on-site and nearest-neighbor repulsions U and V. We show that spin frustration stabilizes a bond-ordered (dimerized) state for U approximately V/2 up to U/t approximately 9, where t is the nearest-neighbor hopping. The transition from the dimerized state to the staggered charge-density-wave state for large V/U is continuous for U < or approximately 5.5 and first order for higher U.

Journal Article↗

Stochastic series expansion method for quantum Ising models with arbitrary interactions.

A quantum Monte Carlo algorithm for the transverse Ising model with arbitrary short- or long-range interactions is presented. The algorithm is based on sampling the diagonal matrix elements of the power-series expansion of the density matrix (stochastic series expansion), and avoids the interaction summations necessary in conventional methods. In the case of long-range interactions, the scaling of the computation time with the system size N is therefore reduced from N2 to N ln(N). The method is tested on a one-dimensional ferromagnet in a transverse field, with interactions decaying as 1/r(2).

Journal Article↗

Susceptibility of the 2D spin-1 / 2 Heisenberg antiferromagnet with an impurity.

We use a quantum Monte Carlo method (stochastic series expansion) to study the effects of a magnetic or nonmagnetic impurity on the magnetic susceptibility of the two-dimensional Heisenberg antiferromagnet. At low temperatures, we find a log-divergent contribution to the transverse susceptibility. We also introduce an effective few-spin model that can quantitatively capture the differences between magnetic and nonmagnetic impurities at high and intermediate temperatures.

Journal Article↗

Multicritical point in a diluted bilayer Heisenberg quantum antiferromagnet.

The S=1/2 Heisenberg bilayer antiferromagnet with randomly removed interlayer dimers is studied using quantum Monte Carlo simulations. A zero-temperature multicritical point (p(*),g(*)) at the classical percolation density p=p(*) and interlayer coupling g(*) approximately equal 0.16 is demonstrated. The quantum critical exponents of the percolating cluster are determined using finite-size scaling. It is argued that the associated finite-temperature quantum critical regime extends to zero interlayer coupling and could be relevant for antiferromagnetic cuprates doped with nonmagnetic impurities.

Journal Article↗

Quantum Monte Carlo with directed loops.

We introduce the concept of directed loops in stochastic series expansion and path-integral quantum Monte Carlo methods. Using the detailed balance rules for directed loops, we show that it is possible to smoothly connect generally applicable simulation schemes (in which it is necessary to include backtracking processes in the loop construction) to more restricted loop algorithms that can be constructed only for a limited range of Hamiltonians (where backtracking can be avoided). The "algorithmic discontinuities" between general and special points (or regions) in parameter space can hence be eliminated. As a specific example, we consider the anisotropic S=1/2 Heisenberg antiferromagnet in an external magnetic field. We show that directed-loop simulations are very efficient for the full range of magnetic fields (zero to the saturation point) and anisotropies. In particular, for weak fields and anisotropies, the autocorrelations are significantly reduced relative to those of previous approaches. The back-tracking probability vanishes continuously as the isotropic Heisenberg point is approached. For the XY model, we show that back tracking can be avoided for all fields extending up to the saturation field. The method is hence particularly efficient in this case. We use directed-loop simulations to study the magnetization process in the two-dimensional Heisenberg model at very low temperatures. For LxL lattices with L up to 64, we utilize the step structure in the magnetization curve to extract gaps between different spin sectors. Finite-size scaling of the gaps gives an accurate estimate of the transverse susceptibility in the thermodynamic limit: chi( perpendicular )=0.0659+/-0.0002.

Journal Article↗