Make life simple: unleash the full power of the parallel tempering algorithm

We introduce a new update scheme to systematically improve the efficiency of parallel tempering simulations. We show that by adapting the number of sweeps between replica exchanges to the canonical autocorrelation time, the average round-trip time of a replica in temperature space can be significantly decreased. The temperatures are not dynamically adjusted as in previous attempts but chosen to yield a 50% exchange rate of adjacent replicas. We illustrate the new algorithm with results for the Ising model in two and the Edwards-Anderson Ising spin glass in three dimensionsComment: 4 pages, 5 figure

Monte Carlo Study of Cluster-Diameter Distribution: A New Observable to Estimate Correlation Lengths

We report numerical simulations of two-dimensional $q$-state Potts models with emphasis on a new quantity for the computation of spatial correlation lengths. This quantity is the cluster-diameter distribution function $G_{diam}(x)$, which measures the distribution of the diameter of stochastically defined cluster. Theoretically it is predicted to fall off exponentially for large diameter $x$, $G_{diam} \propto \exp(-x/\xi)$, where $\xi$ is the correlation length as usually defined through the large-distance behavior of two-point correlation functions. The results of our extensive Monte Carlo study in the disordered phase of the models with $q=10$, 15, and $20$ on large square lattices of size $300 \times 300$, $120 \times 120$, and $80 \times 80$, respectively, clearly confirm the theoretically predicted behavior. Moreover, using this observable we are able to verify an exact formula for the correlation length $\xi_d(\beta_t)$ in the disordered phase at the first-order transition point $\beta_t$ with an accuracy of about $1$ for all considered values of $q$. This is a considerable improvement over estimates derived from the large-distance behavior of standard (projected) two-point correlation functions, which are also discussed for comparison.Comment: 20 pages, LaTeX + 13 postscript figures. See also http://www.cond-mat.physik.uni-mainz.de/~janke/doc/home_janke.htm

Fluctuation Pressure of a Stack of Membranes

We calculate the universal pressure constants of a stack of N membranes between walls by strong-coupling theory. The results are in very good agreement with values from Monte-Carlo simulations.Comment: Author Information under http://www.physik.fu-berlin.de/~kleinert/institution.html Latest update of paper also at http://www.physik.fu-berlin.de/~kleinert/31

Lattice Models of Quantum Gravity

Standard Regge Calculus provides an interesting method to explore quantum gravity in a non-perturbative fashion but turns out to be a CPU-time demanding enterprise. One therefore seeks for suitable approximations which retain most of its universal features. The $Z_2$-Regge model could be such a desired simplification. Here the quadratic edge lengths $q$ of the simplicial complexes are restricted to only two possible values $q=1+\epsilon\sigma$, with $\sigma=\pm 1$, in close analogy to the ancestor of all lattice theories, the Ising model. To test whether this simpler model still contains the essential qualities of the standard Regge Calculus, we study both models in two dimensions and determine several observables on the same lattice size. In order to compare expectation values, e.g. of the average curvature or the Liouville field susceptibility, we employ in both models the same functional integration measure. The phase structure is under current investigation using mean field theory and numerical simulation.Comment: 4 pages, 1 figure

Simplified Transfer Matrix Approach in the Two-Dimensional Ising Model with Various Boundary Conditions

A recent simplified transfer matrix solution of the two-dimensional Ising model on a square lattice with periodic boundary conditions is generalized to periodic-antiperiodic, antiperiodic-periodic and antiperiodic-antiperiodic boundary conditions. It is suggested to employ linear combinations of the resulting partition functions to investigate finite-size scaling. An exact relation of such a combination to the partition function corresponding to Brascamp-Kunz boundary conditions is found.Comment: Phys.Rev.E, to be publishe

Application of Multicanonical Multigrid Monte Carlo Method to the Two-Dimensional ϕ4\phi^4-Model: Autocorrelations and Interface Tension

We discuss the recently proposed multicanonical multigrid Monte Carlo method and apply it to the scalar $\phi^4$-model on a square lattice. To investigate the performance of the new algorithm at the field-driven first-order phase transitions between the two ordered phases we carefully analyze the autocorrelations of the Monte Carlo process. Compared with standard multicanonical simulations a real-time improvement of about one order of magnitude is established. The interface tension between the two ordered phases is extracted from high-statistics histograms of the magnetization applying histogram reweighting techniques.Comment: 49 pp. Latex incl. 14 figures (Fig.7 not included, sorry) as uuencoded compressed tar fil

THE PROFITABILITY OF SUSTAINABLE AGRICULTURE ON A REPRESENTATIVE GRAIN FARM IN THE MID-ATLANTIC REGION, 1981-89

A long-term whole-farm analysis compared conventional and low-input farming systems. Data from a nine-year agronomic study at the Rodale Research Farm, Kutztown, Pennsylvania, were used to analyze profitability, liquidity, solvency, and risk on a representative commercial grain farm. Conventional and low-input farms participating in government programs are the most profitable scenarios, followed by conventional and low-input farms not participating in government programs. All farms increased their net worth. The low-input approach is advantageous for risk-averse farmers using a safety-first criterion.Crop Production/Industries, Environmental Economics and Policy,

Critical exponents of the classical Heisenberg ferromagnet

In a recent letter, R.G. Brown and M. Ciftan (Phys. Rev. Lett. 76, 1352, 1996) reported high precision Monte Carlo (MC) estimates of the static critical exponents of the classical 3D Heisenberg model, which stand in sharp contrast to values obtained by four independent approaches, namely by other recent high statistics MC simulations, high-temperature series analyses, field theoretical methods, and experimental studies. In reply to the above cited work we submitted this paper as a comment to Phys. Rev. Lett