Twisted valence quarks and hadron interactions on the lattice

We consider QCD with valence and sea quarks obeying different boundary conditions. We point out that the energy of low lying two hadron states do not depend on the boundary condition of the sea quarks (up to exponentially small corrections). Thus, the advantages in using twisted boundary conditions on the lattice QCD extraction of nucleon-nucleon phase shifts can be gained without the need of new gauge configurations, even in fully unquenched calculations.Comment: 11 pages, 5 figure

Geometrical and topological issues in octree based automatic meshing

Finite element meshes derived automatically from solid models through recursive spatial subdivision schemes (octrees) can be made to inherit the hierarchical structure and the spatial addressability intrinsic to the underlying grid. These two properties, together with the geometric regularity that can also be built into the mesh, make octree based meshes ideally suited for efficient analysis and self-adaptive remeshing and reanalysis. The element decomposition of the octal cells that intersect the boundary of the domain is discussed. The problem, central to octree based meshing, is solved by combining template mapping and element extraction into a procedure that utilizes both constructive solid geometry and boundary representation techniques. Boundary cells that are not intersected by the edge of the domain boundary are easily mapped to predefined element topology. Cells containing edges (and vertices) are first transformed into a planar polyhedron and then triangulated via element extractor. The modeling environments required for the derivation of planar polyhedra and for element extraction are analyzed

Correcting curvature-density effects in the Hamilton-Jacobi skeleton

The Hainilton-Jacobi approach has proven to be a powerful and elegant method for extracting the skeleton of two-dimensional (2-D) shapes. The approach is based on the observation that the normalized flux associated with the inward evolution of the object boundary at nonskeletal points tends to zero as the size of the integration area tends to zero, while the flux is negative at the locations of skeletal points. Nonetheless, the error in calculating the flux on the image lattice is both limited by the pixel resolution and also proportional to the curvature of the boundary evolution front and, hence, unbounded near endpoints. This makes the exact location of endpoints difficult and renders the performance of the skeleton extraction algorithm dependent on a threshold parameter. This problem can be overcome by using interpolation techniques to calculate the flux with subpixel precision. However, here, we develop a method for 2-D skeleton extraction that circumvents the problem by eliminating the curvature contribution to the error. This is done by taking into account variations of density due to boundary curvature. This yields a skeletonization algorithm that gives both better localization and less susceptibility to boundary noise and parameter choice than the Hamilton-Jacobi method

Characteristic Evolution and Matching

I review the development of numerical evolution codes for general relativity based upon the characteristic initial value problem. Progress in characteristic evolution is traced from the early stage of 1D feasibility studies to 2D axisymmetric codes that accurately simulate the oscillations and gravitational collapse of relativistic stars and to current 3D codes that provide pieces of a binary black hole spacetime. Cauchy codes have now been successful at simulating all aspects of the binary black hole problem inside an artificially constructed outer boundary. A prime application of characteristic evolution is to extend such simulations to null infinity where the waveform from the binary inspiral and merger can be unambiguously computed. This has now been accomplished by Cauchy-characteristic extraction, where data for the characteristic evolution is supplied by Cauchy data on an extraction worldtube inside the artificial outer boundary. The ultimate application of characteristic evolution is to eliminate the role of this outer boundary by constructing a global solution via Cauchy-characteristic matching. Progress in this direction is discussed.Comment: New version to appear in Living Reviews 2012. arXiv admin note: updated version of arXiv:gr-qc/050809

On the extraction of spectral quantities with open boundary conditions

We discuss methods to extract decay constants, meson masses and gluonic observables in the presence of open boundary conditions. The ensembles have been generated by the CLS effort and have 2+1 flavors of O(a)-improved Wilson fermions with a small twisted-mass term as proposed by L\"uscher and Palombi. We analyse the effect of the associated reweighting factors on the computation of different observables.Comment: 7 pages, talk presented at the 32nd International Symposium on Lattice Field Theory - Lattice 2014, Columbia University New Yor

Two body scattering length of Yukawa model on a lattice

The extraction of scattering parameters from Euclidean simulations of a Yukawa model in a finite volume with periodic boundary conditions is analyzed both in non relativistic quantum mechanics and in quantum field theory.Comment: 4 pages, talk at "18th International IUPAP conference on Few Body Problems in Physics" (Sao Paulo, August 2006