Convex Hull of Planar H-Polyhedra

Suppose are planar (convex) H-polyhedra, that is, $A_i \in \mathbb{R}^{n_i \times 2}$ and $\vec{c}_i \in \mathbb{R}^{n_i}$. Let $P_i = \{\vec{x} \in \mathbb{R}^2 \mid A_i\vec{x} \leq \vec{c}_i \}$ and $n = n_1 + n_2$. We present an $O(n \log n)$ algorithm for calculating an H-polyhedron with the smallest $P = \{\vec{x} \in \mathbb{R}^2 \mid A\vec{x} \leq \vec{c} \}$ such that $P_1 \cup P_2 \subseteq P$

A Multi-variate Discrimination Technique Based on Range-Searching

We present a fast and transparent multi-variate event classification technique, called PDE-RS, which is based on sampling the signal and background densities in a multi-dimensional phase space using range-searching. The employed algorithm is presented in detail and its behaviour is studied with simple toy examples representing basic patterns of problems often encountered in High Energy Physics data analyses. In addition an example relevant for the search for instanton-induced processes in deep-inelastic scattering at HERA is discussed. For all studied examples, the new presented method performs as good as artificial Neural Networks and has furthermore the advantage to need less computation time. This allows to carefully select the best combination of observables which optimally separate the signal and background and for which the simulations describe the data best. Moreover, the systematic and statistical uncertainties can be easily evaluated. The method is therefore a powerful tool to find a small number of signal events in the large data samples expected at future particle colliders.Comment: Submitted to NIM, 18 pages, 8 figure

Heat Conduction and Entropy Production in a One-Dimensional Hard-Particle Gas

We present large scale simulations for a one-dimensional chain of hard-point particles with alternating masses. We correct several claims in the recent literature based on much smaller simulations. Both for boundary conditions with two heat baths at different temperatures at both ends and from heat current autocorrelations in equilibrium we find heat conductivities kappa to diverge with the number N of particles. These depended very strongly on the mass ratios, and extrapolation to N -> infty resp. t -> infty is difficult due to very large finite-size and finite-time corrections. Nevertheless, our data seem compatible with a universal power law kappa ~ N^alpha with alpha approx 0.33. This suggests a relation to the Kardar-Parisi-Zhang model. We finally show that the hard-point gas with periodic boundary conditions is not chaotic in the usual sense and discuss why the system, when kept out of equilibrium, leads nevertheless to energy dissipation and entropy production.Comment: 4 pages (incl. 5 figures), RevTe

A Search for Instantons at HERA

A search for QCD instanton (I) induced events in deep-inelastic scattering (DIS) at HERA is presented in the kinematic range of low x and low Q^2. After cutting into three characteristic variables for I-induced events yielding a maximum suppression of standard DIS background to the 0.1% level while still preserving 10% of the I-induced events, 549 data events are found while 363^{+22}_{-26} (CDM) and 435^{+36}_{-22} (MEPS) standard DIS events are expected. More events than expected by the standard DIS Monte Carlo models are found in the data. However, the systematic uncertainty between the two different models is of the order of the expected signal, so that a discovery of instantons can not be claimed. An outlook is given on the prospect to search for QCD instanton events using a discriminant based on range searching in the kinematical region Q^2\gtrsim100\GeV^2 where the I-theory makes safer predictions and the QCD Monte Carlos are expected to better describe the inclusive data.Comment: Invited talk given at the Ringberg Workshop on HERA Physics on June 19th, 2001 on behalf of the H1 collaboratio

Staircase polygons: moments of diagonal lengths and column heights

We consider staircase polygons, counted by perimeter and sums of k-th powers of their diagonal lengths, k being a positive integer. We derive limit distributions for these parameters in the limit of large perimeter and compare the results to Monte-Carlo simulations of self-avoiding polygons. We also analyse staircase polygons, counted by width and sums of powers of their column heights, and we apply our methods to related models of directed walks.Comment: 24 pages, 7 figures; to appear in proceedings of Counting Complexity: An International Workshop On Statistical Mechanics And Combinatorics, 10-15 July 2005, Queensland, Australi

Phase Transition in the Aldous-Shields Model of Growing Trees

We study analytically the late time statistics of the number of particles in a growing tree model introduced by Aldous and Shields. In this model, a cluster grows in continuous time on a binary Cayley tree, starting from the root, by absorbing new particles at the empty perimeter sites at a rate proportional to c^{-l} where c is a positive parameter and l is the distance of the perimeter site from the root. For c=1, this model corresponds to random binary search trees and for c=2 it corresponds to digital search trees in computer science. By introducing a backward Fokker-Planck approach, we calculate the mean and the variance of the number of particles at large times and show that the variance undergoes a `phase transition' at a critical value c=sqrt{2}. While for c>sqrt{2} the variance is proportional to the mean and the distribution is normal, for c<sqrt{2} the variance is anomalously large and the distribution is non-Gaussian due to the appearance of extreme fluctuations. The model is generalized to one where growth occurs on a tree with $m$ branches and, in this more general case, we show that the critical point occurs at c=sqrt{m}.Comment: Latex 17 pages, 6 figure

Efficient Monte Carlo algorithm and high-precision results for percolation

We present a new Monte Carlo algorithm for studying site or bond percolation on any lattice. The algorithm allows us to calculate quantities such as the cluster size distribution or spanning probability over the entire range of site or bond occupation probabilities from zero to one in a single run which takes an amount of time scaling linearly with the number of sites on the lattice. We use our algorithm to determine that the percolation transition occurs at occupation probability 0.59274621(13) for site percolation on the square lattice and to provide clear numerical confirmation of the conjectured 4/3-power stretched-exponential tails in the spanning probability functions.Comment: 8 pages, including 3 postscript figures, minor corrections in this version, plus updated figures for the position of the percolation transitio

Two-Dimensional Quantum XY Model with Ring Exchange and External Field

We present the zero-temperature phase diagram of a square lattice quantum spin 1/2 XY model with four-site ring exchange in a uniform external magnetic field. Using quantum Monte Carlo techniques, we identify various quantum phase transitions between the XY-order, striped or valence bond solid, staggered Neel antiferromagnet and fully polarized ground states of the model. We find no evidence for a quantum spin liquid phase.Comment: 4 pages, 4 figure

Microscopic models for fractionalized phases in strongly correlated systems

We construct explicit examples of microscopic models that stabilize a variety of fractionalized phases of strongly correlated systems in spatial dimension bigger than one, and in zero external magnetic field. These include models of charge fractionalization in boson-only systems, and various kinds of spin-charge separation in electronic systems. We determine the excitation spectrum and show the consistency with that expected from field theoretic descriptions of fractionalization. Our results are further substantiated by direct numerical calculation of the phase diagram of one of the models.Comment: 10 pages, 4 figure

Topological Entanglement Entropy of a Bose-Hubbard Spin Liquid

The Landau paradigm of classifying phases by broken symmetries was demonstrated to be incomplete when it was realized that different quantum Hall states could only be distinguished by more subtle, topological properties. Today, the role of topology as an underlying description of order has branched out to include topological band insulators, and certain featureless gapped Mott insulators with a topological degeneracy in the groundstate wavefunction. Despite intense focus, very few candidates for these topologically ordered "spin liquids" exist. The main difficulty in finding systems that harbour spin liquid states is the very fact that they violate the Landau paradigm, making conventional order parameters non-existent. Here, we uncover a spin liquid phase in a Bose-Hubbard model on the kagome lattice, and measure its topological order directly via the topological entanglement entropy. This is the first smoking-gun demonstration of a non-trivial spin liquid, identified through its entanglement entropy as a gapped groundstate with emergent Z2 gauge symmetry.Comment: 4+ pages, 3 figure