Constrained Orthogonal Polynomials

We define sets of orthogonal polynomials satisfying the additional constraint of a vanishing average. These are of interest, for example, for the study of the Hohenberg-Kohn functional for electronic or nucleonic densities and for the study of density fluctuations in centrifuges. We give explicit properties of such polynomial sets, generalizing Laguerre and Legendre polynomials. The nature of the dimension 1 subspace completing such sets is described. A numerical example illustrates the use of such polynomials.Comment: 11 pages, 10 figure

On the energy momentum dispersion in the lattice regularization

For a free scalar boson field and for U(1) gauge theory finite volume (infrared) and other corrections to the energy-momentum dispersion in the lattice regularization are investigated calculating energy eigenstates from the fall off behavior of two-point correlation functions. For small lattices the squared dispersion energy defined by $E_{\rm dis}^2=E_{\vec{k}}^2-E_0^2-4\sum_{i=1}^{d-1}\sin(k_i/2)^2$ is in both cases negative ($d$ is the Euclidean space-time dimension and $E_{\vec{k}}$ the energy of momentum $\vec{k}$ eigenstates). Observation of $E_{\rm dis}^2=0$ has been an accepted method to demonstrate the existence of a massless photon ($E_0=0$) in 4D lattice gauge theory, which we supplement here by a study of its finite size corrections. A surprise from the lattice regularization of the free field is that infrared corrections do {\it not} eliminate a difference between the groundstate energy $E_0$ and the mass parameter $M$ of the free scalar lattice action. Instead, the relation $E_0=\cosh^{-1} (1+M^2/2)$ is derived independently of the spatial lattice size.Comment: 9 pages, 2 figures. Parts of the paper have been rewritten and expanded to clarify the result

Monte Carlo simulation and global optimization without parameters

We propose a new ensemble for Monte Carlo simulations, in which each state is assigned a statistical weight $1/k$, where $k$ is the number of states with smaller or equal energy. This ensemble has robust ergodicity properties and gives significant weight to the ground state, making it effective for hard optimization problems. It can be used to find free energies at all temperatures and picks up aspects of critical behaviour (if present) without any parameter tuning. We test it on the travelling salesperson problem, the Edwards-Anderson spin glass and the triangular antiferromagnet.Comment: 10 pages with 3 Postscript figures, to appear in Phys. Rev. Lett

Particle Dispersion on Rapidly Folding Random Hetero-Polymers

We investigate the dynamics of a particle moving randomly along a disordered hetero-polymer subjected to rapid conformational changes which induce superdiffusive motion in chemical coordinates. We study the antagonistic interplay between the enhanced diffusion and the quenched disorder. The dispersion speed exhibits universal behavior independent of the folding statistics. On the other hand it is strongly affected by the structure of the disordered potential. The results may serve as a reference point for a number of translocation phenomena observed in biological cells, such as protein dynamics on DNA strands.Comment: 4 pages, 4 figure

Drops with non-circular footprints

In this paper we study the morphology of drops formed on partially wetting substrates, whose footprint is not circular. This type of drops is a consequence of the breakup processes occurring in thin films when anisotropic contact line motions take place. The anisotropy is basically due to hysteresis effects of the contact angle since some parts of the contact line are wetting, while others are dewetting. Here, we obtain a peculiar drop shape from the rupture of a long liquid filament sitting on a solid substrate, and analyze its shape and contact angles by means of goniometric and refractive techniques. We also find a non--trivial steady state solution for the drop shape within the long wave approximation (lubrication theory), and compare most of its features with experimental data. This solution is presented both in Cartesian and polar coordinates, whose constants must be determined by a certain group of measured parameters. Besides, we obtain the dynamics of the drop generation from numerical simulations of the full Navier--Stokes equation, where we emulate the hysteretic effects with an appropriate spatial distribution of the static contact angle over the substrate

A Decision Support Tool for Seed Mixture Calculations

Grassland species are normally seeded in mixtures rather than monocultures. In theory, seeding rates for mixtures are simply a sum of the amount of pure live seed (PLS) of each seed lot in the mix, an amount sufficient to ensure establishment and survival of each species. Mixtures can be complex because of the number of species used (especially in conservation and reclamation programs) and variations in seed purity and seed size. Soil limitations and seeding equipment settings need to be considered and in Canada, a metric conversion may be required. All these conditions make by-hand calculations of mixtures containing more than 3 species tedious and complicated. Thus, in practice, agronomists and growers use simple rules to set rates. The easiest rule is to estimate the mixture’s components as a percentage by weight of a standardized total weight of the seed required (e.g. 10% of 10 kg/ha). The resulting errors can be observed in the predominance of thin stands, the unexpected dominance of small seeded species and the added costs of interseeding to compete with weeds and fertilizer to increase yield. The objective of this project was to develop a decision support tool, a seed mixture calculator to simplify conversion and improve the estimates of seed required for individual seeding projects

Particles held by springs in a linear shear flow exhibit oscillatory motion

The dynamics of small spheres, which are held by linear springs in a low Reynolds number shear flow at neighboring locations is investigated. The flow elongates the beads and the interplay of the shear gradient with the nonlinear behavior of the hydrodynamic interaction among the spheres causes in a large range of parameters a bifurcation to a surprising oscillatory bead motion. The parameter ranges, wherein this bifurcation is either super- or subcritical, are determined.Comment: 4 pages, 5 figure

Multicanonical Recursions

The problem of calculating multicanonical parameters recursively is discussed. I describe in detail a computational implementation which has worked reasonably well in practice.Comment: 23 pages, latex, 4 postscript figures included (uuencoded Z-compressed .tar file created by uufiles), figure file corrected