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

An efficient scheme for numerical simulations of the spin-bath decoherence

We demonstrate that the Chebyshev expansion method is a very efficient numerical tool for studying spin-bath decoherence of quantum systems. We consider two typical problems arising in studying decoherence of quantum systems consisting of few coupled spins: (i) determining the pointer states of the system, and (ii) determining the temporal decay of quantum oscillations. As our results demonstrate, for determining the pointer states, the Chebyshev-based scheme is at least a factor of 8 faster than existing algorithms based on the Suzuki-Trotter decomposition. For the problems of second type, the Chebyshev-based approach has been 3--4 times faster than the Suzuki-Trotter-based schemes. This conclusion holds qualitatively for a wide spectrum of systems, with different spin baths and different Hamiltonians.Comment: 8 pages (RevTeX), 3 EPS figure

A Constrained Path Quantum Monte Carlo Method for Fermion Ground States

We propose a new quantum Monte Carlo algorithm to compute fermion ground-state properties. The ground state is projected from an initial wavefunction by a branching random walk in an over-complete basis space of Slater determinants. By constraining the determinants according to a trial wavefunction $|\Psi_T \rangle$, we remove the exponential decay of signal-to-noise ratio characteristic of the sign problem. The method is variational and is exact if $|\Psi_T\rangle$ is exact. We report results on the two-dimensional Hubbard model up to size $16\times 16$, for various electron fillings and interaction strengths.Comment: uuencoded compressed postscript file. 5 pages with 1 figure. accepted by PRL

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

Thermal transport in a granular metal array

We obtain the Kubo formula for the electronic thermal conductivity kappa(T) of a granular metal array at low temperatures for the Ambegaokar-Eckern-Schoen (AES) model and study the kinetic and potential contributions in the diamagnetic (local) and paramagnetic (current-current) terms. For small values of dimensionless intergrain tunneling conductance, g << 1, we show that inelastic cotunneling processes contribute to thermal conductivity due to non-cancellation of the diamagnetic and paramagnetic terms, unlike electrical conductivity. We find that the electrical conductivity obeys the Arrhenius law, sigma(T) ~ ge^{-E_c/T}, however kappa(T) decreases only algebraically, kappa(T) \~ g^2 T^3/E_c^2. At large values of intergrain coupling, g >> 1, we find it plausible that the Wiedemann-Franz law weakly deviates from the free-electron theory due to Coulomb effects.Comment: 5 pages RevTeX, to appear in Physical Review Letter