Green Function Monte Carlo with Stochastic Reconfiguration

A new method for the stabilization of the sign problem in the Green Function Monte Carlo technique is proposed. The method is devised for real lattice Hamiltonians and is based on an iterative ''stochastic reconfiguration'' scheme which introduces some bias but allows a stable simulation with constant sign. The systematic reduction of this bias is in principle possible. The method is applied to the frustrated J1-J2 Heisenberg model, and tested against exact diagonalization data. Evidence of a finite spin gap for J2/J1 >~ 0.4 is found in the thermodynamic limit.Comment: 13 pages, RevTeX + 3 encapsulated postscript figure

Direct magneto-optical compression of an effusive atomic beam for high-resolution focused ion beam application

An atomic rubidium beam formed in a 70 mm long two-dimensional magneto-optical trap (2D MOT), directly loaded from a collimated Knudsen source, is analyzed using laser-induced fluorescence. The longitudinal velocity distribution, the transverse temperature and the flux of the atomic beam are reported. The equivalent transverse reduced brightness of an ion beam with similar properties as the atomic beam is calculated because the beam is developed to be photoionized and applied in a focused ion beam. In a single two-dimensional magneto-optical trapping step an equivalent transverse reduced brightness of $(1.0\substack{+0.8-0.4})$ $\times 10^6$ A/(m$^2$ sr eV) was achieved with a beam flux equivalent to $(0.6\substack{+0.3-0.2})$ nA. The temperature of the beam is further reduced with an optical molasses after the 2D MOT. This increased the equivalent brightness to $(6\substack{+5-2})$$\times 10^6$ A/(m$^2$ sr eV). For currents below 10 pA, for which disorder-induced heating can be suppressed, this number is also a good estimate of the ion beam brightness that can be expected. Such an ion beam brightness would be a six times improvement over the liquid metal ion source and could improve the resolution in focused ion beam nanofabrication.Comment: 10 pages, 8 figures, 1 tabl

Optimization of Gutzwiller Wavefunctions in Quantum Monte Carlo

Gutzwiller functions are popular variational wavefunctions for correlated electrons in Hubbard models. Following the variational principle, we are interested in the Gutzwiller parameters that minimize e.g. the expectation value of the energy. Rewriting the expectation value as a rational function in the Gutzwiller parameters, we find a very efficient way for performing that minimization. The method can be used to optimize general Gutzwiller-type wavefunctions both, in variational and in fixed-node diffusion Monte Carlo.Comment: 9 pages RevTeX with 10 eps figure

Sign Rules for Anisotropic Quantum Spin Systems

We present new and exact ``sign rules'' for various spin-s anisotropic spin-lattice models. It is shown that, after a simple transformation which utilizes these sign rules, the ground-state wave function of the transformed Hamiltonian is positive-definite. Using these results exact statements for various expectation values of off-diagonal operators are presented, and transitions in the behavior of these expectation values are observed at particular values of the anisotropy. Furthermore, the effects of sign rules in variational calculations and quantum Monte Carlo calculations are considered. They are illustrated by a simple variational treatment of a one-dimensional anisotropic spin model.Comment: 4 pages, 1 ps-figur

Charge fluctuations close to phase separation in the two dimensional t-J model

We have studied the t-J model using the Green Function Monte Carlo technique. We have obtained accurate energies well converged in the thermodynamic limit, by performing simulations up to 242 lattice sites. By studying the energy as a function of hole doping we conclude that there is no phase separation in the physical region, relevant for HTc superconductors. This finding is further supported by the hole-hole correlation function calculation. Remarkably, by approaching the phase separation instability, for $J_c/t\sim 0.5$,this function displays enhanced fluctuations at incommensurate wavevectors, scaling linearly with the doping, in agreement with experimental findings.Comment: To appear on Phys. Rev. Let

Green Function Monte Carlo with Stochastic Reconfiguration: an effective remedy for the sign problem disease

A recent technique, proposed to alleviate the ``sign problem disease'', is discussed in details. As well known the ground state of a given Hamiltonian $H$ can be obtained by applying the imaginary time propagator $e^{-H \tau}$ to a given trial state $\psi_T$ for large imaginary time $\tau$ and sampling statistically the propagated state $\psi_{\tau} = e^{-H \tau} \psi_T$. However the so called ``sign problem'' may appear in the simulation and such statistical propagation would be practically impossible without employing some approximation such as the well known ``fixed node'' approximation (FN). This method allows to improve the FN dynamic with a systematic correction scheme. This is possible by the simple requirement that, after a short imaginary time propagation via the FN dynamic, a number $p$ of correlation functions can be further constrained to be {\em exact} by small perturbation of the FN propagated state, which is free of the sign problem. By iterating this scheme the Monte Carlo average sign, which is almost zero when there is sign problem, remains stable and finite even for large $\tau$. The proposed algorithm is tested against the exact diagonalization results available on finite lattice. It is also shown in few test cases that the dependence of the results upon the few parameters entering the stochastic technique can be very easily controlled, unless for exceptional cases.Comment: 44 pages, RevTeX + 5 encaplulated postscript figure

From antiferromagnetism to d-wave superconductivity in the 2D t-J model

We have found that the two dimensional t-J model, for the physical parameter range J/t = 0.4 reproduces the main experimental qualitative features of High-Tc copper oxide superconductors: d-wave superconducting correlations are strongly enhanced upon small doping and clear evidence of off diagonal long range order is found at the optimal doping \delta ~ 0.15. On the other hand antiferromagnetic long range order, clearly present at zero hole doping, is suppressed at small hole density with clear absence of antiferromagnetism at \delta >~ 0.1.Comment: 4 pages, 5 figure

Screening, Coulomb pseudopotential, and superconductivity in alkali-doped Fullerenes

We study the static screening in a Hubbard-like model using quantum Monte Carlo. We find that the random phase approximation is surprisingly accurate almost up to the Mott transition. We argue that in alkali-doped Fullerenes the Coulomb pseudopotential $\mu^\ast$ is not very much reduced by retardation effects. Therefore efficient screening is important in reducing $\mu^{\ast}$ sufficiently to allow for an electron-phonon driven superconductivity. In this way the Fullerides differ from the conventional picture, where retardation effects play a major role in reducing the electron-electron repulsion.Comment: 4 pages RevTeX with 2 eps figures, additional material available at http://www.mpi-stuttgart.mpg.de/docs/ANDERSEN/fullerene

Random Exchange Quantum Heisenberg Chains

The one-dimensional quantum Heisenberg model with random $\pm J$ bonds is studied for $S=\frac{1}{2}$ and $S=1$. The specific heat and the zero-field susceptibility are calculated by using high-temperature series expansions and quantum transfer matrix method. The susceptibility shows a Curie-like temperature dependence at low temperatures as well as at high temperatures. The numerical results for the specific heat suggest that there are anomalously many low-lying excitations. The qualitative nature of these excitations is discussed based on the exact diagonalization of finite size systems.Comment: 13 pages, RevTex, 12 figures available on request ([email protected]

Spontaneous plaquette dimerization in the J1−J2J_1-J_2 Heisenberg model

We investigate the non magnetic phase of the spin-half frustrated Heisenberg antiferromagnet on the square lattice using exact diagonalization (up to 36 sites) and quantum Monte Carlo techniques (up to 144 sites). The spin gap and the susceptibilities for the most important crystal symmetry breaking operators are computed. A genuine and somehow unexpected `plaquette RVB', with spontaneously broken translation symmetry and no broken rotation symmetry, comes out from our numerical simulations as the most plausible ground state for $J_2/J_1 \simeq 0.5$.Comment: 4 pages, 5 postscript figure