Some Comments on Multigrid Methods for Computing Propagators

I make three conceptual points regarding multigrid methods for computing propagators in lattice gauge theory: 1) The class of operators handled by the algorithm must be stable under coarsening. 2) Problems related by symmetry should have solution methods related by symmetry. 3) It is crucial to distinguish the vector space $V$ from its dual space $V^*$. All the existing algorithms violate one or more of these principles.Comment: 16 pages, LaTeX plus subeqnarray.sty (included at end), NYU-TH-93/07/0

Measuring the Decorrelation Times of Fourier Modes in Simulations

We describe a method to study the rate at which modes decorrelate in numerical simulations. We study the XY model updated with the Metropolis and Wolff dynamics respectively and compute the rate at which each eigenvector of the dynamics decorrelates. Our method allows us to identify the decorrelation time for each mode separately. We find that the autocorrelation function of the various modes is markedly different for the `local' Metropolis compared to the `non-local' Wolff dynamics. Equipped with this new insight, it may be possible to devise highly efficient algorithms.Comment: 16 pp (LaTeX), PUPT-1378 , IASSNS-HEP-93/

Online fabrication and characterization of capsule populations with a flow-focusing microfluidic system

This paper was presented at the 3rd Micro and Nano Flows Conference (MNF2011), which was held at the Makedonia Palace Hotel, Thessaloniki in Greece. The conference was organised by Brunel University and supported by the Italian Union of Thermofluiddynamics, Aristotle University of Thessaloniki, University of Thessaly, IPEM, the Process Intensification Network, the Institution of Mechanical Engineers, the Heat Transfer Society, HEXAG - the Heat Exchange Action Group, and the Energy Institute.We have designed a microfluidic system that combines a double flow-focusing setup for calibrated capsule fabrication with a microchannel for the characterization of their mechanical properties. The double flow-focusing system consists of a first Y junction to create the microdroplets and of a second Y junction to introduce the cross-linking agent allowing the membrane formation. The human serum albumin (HSA) aqueous solution for the dispersed solution, hydrophobic phase for the continuous solution and cross-linking agent solution are introduced by means of syringe pumps. A wavy channel after the second junction allows to control the reticulation time. A cylindrical microchannel then enables to deform and characterize the capsules formed. The mechanical properties of the capsule membrane are obtained by inverse analysis (Chu et al. 2011). The results show that the drop size increases with the flow rate ratio between the central and lateral channels and does not change much regardless of the flow rate of the reticulation phase. The mean shear modulus of the capsules fabricated after 23 s of reticulation is of the order of the surface tension of HSA solution with Dragoxat indicating that the reticulation time is too short to form an elastic membrane around the droplet. When the reticulation time is increased to 60 s, the membrane shear modulus is multiplied by a factor of 3 confirming that a solid membrane has formed around the drop

Anomalous diffusion in polymers: long-time behaviour

We study the Dirichlet boundary value problem for viscoelastic diffusion in polymers. We show that its weak solutions generate a dissipative semiflow. We construct the minimal trajectory attractor and the global attractor for this problem.Comment: 13 page

Dgsos on DTRS

We perform simulations of a discrete gaussian solid on solid (DGSOS) model on dynamical $\phi^3$ graphs, which is equivalent to coupling the model to 2d quantum gravity, using the cluster algorithms recently developed by Evertz et.al.for use on fixed lattices. We find evidence from the growth of the width-squared in the rough phase of KT-like behaviour, which is consistent with theoretical expectations. We also investigate the cluster statistics, dynamical critical exponent and lattice properties, and compare these with the dual XY model.Comment: 9 pages, COLO-HEP-32

Multigrid Monte Carlo Algorithms for SU(2) Lattice Gauge Theory: Two versus Four Dimensions

We study a multigrid method for nonabelian lattice gauge theory, the time slice blocking, in two and four dimensions. For SU(2) gauge fields in two dimensions, critical slowing down is almost completely eliminated by this method. This result is in accordance with theoretical arguments based on the analysis of the scale dependence of acceptance rates for nonlocal Metropolis updates. The generalization of the time slice blocking to SU(2) in four dimensions is investigated analytically and by numerical simulations. Compared to two dimensions, the local disorder in the four dimensional gauge field leads to kinematical problems.Comment: 24 pages, PostScript file (compressed and uuencoded), preprint MS-TPI-94-

Kinematics of Multigrid Monte Carlo

We study the kinematics of multigrid Monte Carlo algorithms by means of acceptance rates for nonlocal Metropolis update proposals. An approximation formula for acceptance rates is derived. We present a comparison of different coarse-to-fine interpolation schemes in free field theory, where the formula is exact. The predictions of the approximation formula for several interacting models are well confirmed by Monte Carlo simulations. The following rule is found: For a critical model with fundamental Hamiltonian H(phi), absence of critical slowing down can only be expected if the expansion of in terms of the shift psi contains no relevant (mass) term. We also introduce a multigrid update procedure for nonabelian lattice gauge theory and study the acceptance rates for gauge group SU(2) in four dimensions.Comment: 28 pages, 8 ps-figures, DESY 92-09

Multigrid Methods in Lattice Field Computations

The multigrid methodology is reviewed. By integrating numerical processes at all scales of a problem, it seeks to perform various computational tasks at a cost that rises as slowly as possible as a function of $n$, the number of degrees of freedom in the problem. Current and potential benefits for lattice field computations are outlined. They include: $O(n)$ solution of Dirac equations; just $O(1)$ operations in updating the solution (upon any local change of data, including the gauge field); similar efficiency in gauge fixing and updating; $O(1)$ operations in updating the inverse matrix and in calculating the change in the logarithm of its determinant; $O(n)$ operations per producing each independent configuration in statistical simulations (eliminating CSD), and, more important, effectively just $O(1)$ operations per each independent measurement (eliminating the volume factor as well). These potential capabilities have been demonstrated on simple model problems. Extensions to real life are explored.Comment: 4

Quenching 2D Quantum Gravity

We simulate the Ising model on a set of fixed random $\phi^3$ graphs, which corresponds to a {\it quenched} coupling to 2D gravity rather than the annealed coupling that is usually considered. We investigate the critical exponents in such a quenched ensemble and compare them with measurements on dynamical $\phi^3$ graphs, flat lattices and a single fixed $\phi^3$ graph.Comment: 8 page

On some discrete boundary value problems in canonical domains

We study some discrete boundary value problems for discrete elliptic pseudo-differential equations in a half-space. These statements are related with a special periodic factorization of an elliptic symbol and a number of boundary conditions depends on an index of periodic factorizatio