"Exact" Algorithm for Random-Bond Ising Models in 2D

We present an efficient algorithm for calculating the properties of Ising models in two dimensions, directly in the spin basis, without the need for mapping to fermion or dimer models. The algorithm gives numerically exact results for the partition function and correlation functions at a single temperature on any planar network of N Ising spins in O(N^{3/2}) time or less. The method can handle continuous or discrete bond disorder and is especially efficient in the case of bond or site dilution, where it executes in O(L^2 ln L) time near the percolation threshold. We demonstrate its feasibility on the ferromagnetic Ising model and the +/- J random-bond Ising model (RBIM) and discuss the regime of applicability in cases of full frustration such as the Ising antiferromagnet on a triangular lattice.Comment: 4.2 pages, 5 figures, accepted for publication in Phys. Rev. Let

Noise Predictions for STM in Systems with Local Electron Nematic Order

We propose that thermal noise in local stripe orientation should be readily detectable via STM on systems in which local stripe orientations are strongly affected by quenched disorder. Stripes, a unidirectional, nanoscale modulation of electronic charge, are strongly affected by quenched disorder in two-dimensional and quasi-two-dimensional systems. While stripe orientations tend to lock to major lattice directions, dopant disorder locally breaks rotational symmetry. In a host crystal with otherwise $C_4$ rotational symmetry, stripe orientations in the presence of quenched disorder map to the random field Ising model. While the low temperature state of such a system is generally a stripe glass in two dimensional or strongly layered systems, as the temperature is raised, stripe orientational fluctuations become more prevalent. We propose that these thermally excited fluctuations should be readily detectable in scanning tunneling spectroscopy as {\em telegraph noise} in the high voltage part of the local $I(V)$ curves. We predict the spatial, temporal, and thermal evolution of such noise, including the circumstances under which such noise is most likely to be observed. In addition, we propose an in-situ test, amenable to any local scanning probe, for assessing whether such noise is due to correlated fluctuations rather than independent switchers.Comment: 8 pages, 8 figure

Vlasov Description Of Dense Quark Matter

We discuss properties of quark matter at finite baryon densities and zero temperature in a Vlasov approach. We use a screened interquark Richardson's potential consistent with the indications of Lattice QCD calculations. We analyze the choices of the quark masses and the parameters entering the potential which reproduce the binding energy (B.E.) of infinite nuclear matter. There is a transition from nuclear to quark matter at densities 5 times above normal nuclear matter density. The transition could be revealed from the determination of the position of the shifted meson masses in dense baryonic matter. A scaling form of the meson masses in dense matter is given.Comment: 15 pages 4 figure

Dimers on the Triangular Kagome Lattice

We derive exact results for close-packed dimers on the triangular kagome lattice (TKL), formed by inserting triangles into the triangles of the kagome lattice. Because the TKL is a non-bipartite lattice, dimer-dimer correlations are short-ranged, so that the ground state at the Rokhsar-Kivelson (RK) point of the corresponding quantum dimer model on the same lattice is a short-ranged spin liquid. Using the Pfaffian method, we derive an exact form for the free energy, and we find that the entropy is 1/3 ln2 per site, regardless of the weights of the bonds. The occupation probability of every bond is 1/4 in the case of equal weights on every bond. Similar to the case of lattices formed by corner-sharing triangles (such as the kagome and squagome lattices), we find that the dimer-dimer correlation function is identically zero beyond a certain (short) distance. We find in addition that monomers are deconfined on the TKL, indicating that there is a short-ranged spin liquid phase at the RK point. We also find exact results for the ground state energy of the classical Heisenberg model. The ground state can be ferromagnetic, ferrimagnetic, locally coplanar, or locally canted, depending on the couplings. From the dimer model and the classical spin model, we derive upper bounds on the ground state energy of the quantum Heisenberg model on the TKL.Comment: 9 pages, 7 figures, http://www.physics.purdue.edu/~dyao

Practical solution to the Monte Carlo sign problem: Realistic calculations of 54Fe

We present a practical solution to the "sign problem" in the auxiliary field Monte Carlo approach to the nuclear shell model. The method is based on extrapolation from a continuous family of problem-free Hamiltonians. To demonstrate the resultant ability to treat large shell-model problems, we present results for 54Fe in the full fp-shell basis using the Brown-Richter interaction. We find the Gamow-Teller beta^+ strength to be quenched by 58% relative to the single-particle estimate, in better agreement with experiment than previous estimates based on truncated bases.Comment: 11 pages + 2 figures (not included

Bond-Propagation Algorithm for Thermodynamic Functions in General 2D Ising Models

Recently, we developed and implemented the bond propagation algorithm for calculating the partition function and correlation functions of random bond Ising models in two dimensions. The algorithm is the fastest available for calculating these quantities near the percolation threshold. In this paper, we show how to extend the bond propagation algorithm to directly calculate thermodynamic functions by applying the algorithm to derivatives of the partition function, and we derive explicit expressions for this transformation. We also discuss variations of the original bond propagation procedure within the larger context of Y-Delta-Y-reducibility and discuss the relation of this class of algorithm to other algorithms developed for Ising systems. We conclude with a discussion on the outlook for applying similar algorithms to other models.Comment: 12 pages, 10 figures; submitte

The structure of EAS at E 0.1 EeV

The ratio of extensive air showers (EAS) total shower energy in the electromagnetic channel (E em) to the size of the shower at maximum development (N max) from a direct measurement of shower longitudinal development using the air fluorescence technique was calculated. The values are not inconsistent with values based upon track length integrals of the Gaisser-Hillas formula for shower development or the known relation between shower energy and size at maximum for pure electromagnetic cascades. Using Linsley's estimates for undetected shower energy based on an analysis of a wide variety of cosmic ray data, the following relation for total shower energy E vs N max is obtained. The Gaisser Hillas implied undetected shower energy fractions

Algorithm for Linear Response Functions at Finite Temperatures: Application to ESR spectrum of s=1/2 Antiferromagnet Cu benzoate

We introduce an efficient and numerically stable method for calculating linear response functions $\chi(\vec{q},\omega)$ of quantum systems at finite temperatures. The method is a combination of numerical solution of the time-dependent Schroedinger equation, random vector representation of trace, and Chebyshev polynomial expansion of Boltzmann operator. This method should be very useful for a wide range of strongly correlated quantum systems at finite temperatures. We present an application to the ESR spectrum of s=1/2 antiferromagnet Cu benzoate.Comment: 4 pages, 4 figure

Photoemission spectra of many-polaron systems

The cross over from low to high carrier densities in a many-polaron system is studied in the framework of the one-dimensional spinless Holstein model, using unbiased numerical methods. Combining a novel quantum Monte Carlo approach and exact diagonalization, accurate results for the single-particle spectrum and the electronic kinetic energy on fairly large systems are obtained. A detailed investigation of the quality of the Monte Carlo data is presented. In the physically most important adiabatic intermediate electron-phonon coupling regime, for which no analytical results are available, we observe a dissociation of polarons with increasing band filling, leading to normal metallic behavior, while for parameters favoring small polarons, no such density-driven changes occur. The present work points towards the inadequacy of single-polaron theories for a number of polaronic materials such as the manganites.Comment: 15 pages, 13 figures; final version, accepted for publication in Phys. Rev.