General Algorithm For Improved Lattice Actions on Parallel Computing Architectures

Quantum field theories underlie all of our understanding of the fundamental forces of nature. The are relatively few first principles approaches to the study of quantum field theories [such as quantum chromodynamics (QCD) relevant to the strong interaction] away from the perturbative (i.e., weak-coupling) regime. Currently the most common method is the use of Monte Carlo methods on a hypercubic space-time lattice. These methods consume enormous computing power for large lattices and it is essential that increasingly efficient algorithms be developed to perform standard tasks in these lattice calculations. Here we present a general algorithm for QCD that allows one to put any planar improved gluonic lattice action onto a parallel computing architecture. High performance masks for specific actions (including non-planar actions) are also presented. These algorithms have been successfully employed by us in a variety of lattice QCD calculations using improved lattice actions on a 128 node Thinking Machines CM-5. {\underline{Keywords}}: quantum field theory; quantum chromodynamics; improved actions; parallel computing algorithms

The Potts-q random matrix model : loop equations, critical exponents, and rational case

In this article, we study the q-state Potts random matrix models extended to branched polymers, by the equations of motion method. We obtain a set of loop equations valid for any arbitrary value of q. We show that, for q=2-2 \cos {l \over r} \pi (l, r mutually prime integers with l < r), the resolvent satisfies an algebraic equation of degree 2 r -1 if l+r is odd and r-1 if l+r is even. This generalizes the presently-known cases of q=1, 2, 3. We then derive for any 0 \leq q \leq 4 the Potts-q critical exponents and string susceptibility.Comment: 7 pages, submitted to Phys. Letters

Comparison of fragment partitions production in peripheral and central collisions

Ensembles of single-source events, produced in peripheral and central collisions and correponding respectively to quasi-projectile and quasi-fusion sources, are analyzed. After selections on fragment kinematic properties, excitation energies of the sources are derived using the calorimetric method and the mean behaviour of fragments of the two ensembles are compared. Differences observed in their partitions, especially the charge asymmetry, can be related to collective energy deposited in the systems during the collisions.Comment: 7 pages, 2 figures, presented at the International Workshop on Multifragmentation and Related Topics, Caen France, 4-7th november 2007 (IWM2007

Sounding of finite solid bodies by way of topological derivative

This paper is concerned with an application of the concept of topological derivative to elastic-wave imaging of finite solid bodies containing cavities. Building on the approach originally proposed in the (elastostatic) theory of shape optimization, the topological derivative, which quantifies the sensitivity of a featured cost functional due to the creation of an infinitesimal hole in the cavity-free (reference) body, is used as a void indicator through an assembly of sampling points where it attains negative values. The computation of topological derivative is shown to involve an elastodynamic solution to a set of supplementary boundary-value problems for the reference body, which are here formulated as boundary integral equations. For a comprehensive treatment of the subject, formulas for topological sensitivity are obtained using three alternative methodologies, namely (i) direct differentiation approach, (ii) adjoint field method, and (iii) limiting form of the shape sensitivity analysis. The competing techniques are further shown to lead to distinct computational procedures. Methodologies (i) and (ii) are implemented within a BEM-based platform and validated against an analytical solution. A set of numerical results is included to illustrate the utility of topological derivative for 3D elastic-wave sounding of solid bodies; an approach that may perform best when used as a pre-conditioning tool for more accurate, gradient-based imaging algorithms. Despite the fact that the formulation and results presented in this investigation are established on the basis of a boundary integral solution, the proposed methodology is readily applicable to other computational platforms such as the finite element and finite difference techniques

Topological derivative for the inverse scattering of elastic waves

To establish an alternative analytical framework for the elastic-wave imaging of underground cavities, the focus of this study is an extension of the concept of topological derivative, rooted in elastostatics and shape optimization, to three-dimensional elastodynamics involving semi-infinite and infinite solids. The main result of the proposed boundary integral approach is a formula for topological derivative, explicit in terms of the elastodynamic fundamental solution, obtained by an asymptotic expansion of the misfit-type cost functional with respect to the creation of an infinitesimal hole in an otherwise intact (semi-infinite or infinite) elastic medium. Valid for an arbitrary shape of the infinitesimal cavity, the formula involves the solution of six canonical exterior elastostatic problems, and becomes fully explicit when the vanishing cavity is spherical. A set of numerical results is included to illustrate the potential of topological derivative as a computationally efficient tool for exposing an approximate cavity topology, location, and shape via a grid-type exploration of the host solid. For a comprehensive solution to three-dimensional inverse scattering problems involving elastic waves, the proposed approach can be used most effectively as a pre-conditioning tool for more refined, albeit computationally intensive minimization-based imaging algorithms. To the authors' knowledge, an application of topological derivative to inverse scattering problems has not been attempted before; the methodology proposed in this paper could also be extended to acoustic problems

Elastic-wave identification of penetrable obstacles using shape-material sensitivity framework

This study deals with elastic-wave identification of discrete heterogeneities (inclusions) in an otherwise homogeneous ``reference'' solid from limited-aperture waveform measurements taken on its surface. On adopting the boundary integral equation (BIE) framework for elastodynamic scattering, the inverse query is cast as a minimization problem involving experimental observations and their simulations for a trial inclusion that is defined through its boundary, elastic moduli, and mass density. For an optimal performance of the gradient-based search methods suited to solve the problem, explicit expressions for the shape (i.e. boundary) and material sensitivities of the misfit functional are obtained via the adjoint field approach and direct differentiation of the governing BIEs. Making use of the message-passing interface, the proposed sensitivity formulas are implemented in a data-parallel code and integrated into a nonlinear optimization framework based on the direct BIE method and an augmented Lagrangian whose inequality constraints are employed to avoid solving forward scattering problems for physically inadmissible (or overly distorted) trial inclusion configurations. Numerical results for the reconstruction of an ellipsoidal defect in a semi-infinite solid show the effectiveness of the proposed shape-material sensitivity formulation, which constitutes an essential computational component of the defect identification algorithm

Small-inclusion asymptotic of misfit functionals for inverse problems in acoustics

The aim of this study is an extension and employment of the concept of topological derivative as it pertains to the nucleation of infinitesimal inclusions in a reference (i.e. background) acoustic medium. The developments are motivated by the need to develop a preliminary indicator functional that would aid the solution of inverse scattering problems in terms of a rational initial 'guess' about the geometry and material characteristics of a hidden (finite) obstacle; an information that is often required by iterative minimization algorithms. To this end the customary definition of topological derivative, which quantifies the sensitivity of a given cost functional with respect to the creation of an infinitesimal hole, is adapted to permit the nucleation of a dissimilar acoustic medium. On employing the Green's function for the background domain, computation of topological sensitivity for the three-dimensional Helmholtz equation is reduced to the solution of a reference, Laplace transmission problem. Explicit formulae are given for the nucleating inclusions of spherical and ellipsoidal shapes. For generality the developments are also presented in an alternative, adjoint-field setting that permits nucleation of inclusions in an infinite, semi-infinite or finite background medium. Through numerical examples, it is shown that the featured topological sensitivity could be used, in the context of inverse scattering, as an effective obstacle indicator through an assembly of sampling points where it attains pronounced negative values. On varying a material characteristic (density) of the nucleating obstacle, it is also shown that the proposed methodology can be used as a preparatory tool for both geometric and material identification

Scaling Behavior of the Landau Gauge Overlap Quark Propagator

The properties of the momentum space quark propagator in Landau gauge are examined for the overlap quark action in quenched lattice QCD. Numerical calculations are done on three lattices with different lattice spacings and similar physical volumes to explore the approach of the quark propagator towards the continuum limit. We have calculated the nonperturbative momentum-dependent wavefunction renormalization function $Z(p^2)$ and the nonperturbative mass function $M(p^2)$ for a variety of bare quark masses and extrapolate to the chiral limit. We find the behavior of $Z(p^2)$ and $M(p^2)$ are in good agreement for the two finer lattices in the chiral limit. The quark condensate is also calculated.Comment: 3 pages, Lattice2003(Chiral fermions