Instabilities and Non-Reversibility of Molecular Dynamics Trajectories

The theoretical justification of the Hybrid Monte Carlo algorithm depends upon the molecular dynamics trajectories within it being exactly reversible. If computations were carried out with exact arithmetic then it would be easy to ensure such reversibility, but the use of approximate floating point arithmetic inevitably introduces violations of reversibility. In the absence of evidence to the contrary, we are usually prepared to accept that such rounding errors can be made small enough to be innocuous, but in certain circumstances they are exponentially amplified and lead to blatantly erroneous results. We show that there are two types of instability of the molecular dynamics trajectories which lead to this behavior, instabilities due to insufficiently accurate numerical integration of Hamilton's equations, and intrinsic chaos in the underlying continuous fictitious time equations of motion themselves. We analyze the former for free field theory, and show that it is essentially a finite volume effect. For the latter we propose a hypothesis as to how the Liapunov exponent describing the chaotic behavior of the fictitious time equations of motion for an asymptotically free quantum field theory behaves as the system is taken to its continuum limit, and explain why this means that instabilities in molecular dynamics trajectories are not a significant problem for Hybrid Monte Carlo computations. We present data for pure $SU(3)$ gauge theory and for QCD with dynamical fermions on small lattices to illustrate and confirm some of our results.Comment: 28 pages latex with 19 color postscript figures included by eps

Chiral properties of domain-wall fermions from eigenvalues of 4 dimensional Wilson-Dirac operator

We investigate chiral properties of the domain-wall fermion (DWF) system by using the four-dimensional hermitian Wilson-Dirac operator. We first derive a formula which connects a chiral symmetry breaking term in the five dimensional DWF Ward-Takahashi identity with the four dimensional Wilson-Dirac operator, and simplify the formula in terms of only the eigenvalues of the operator, using an ansatz for the form of the eigenvectors. For a given distribution of the eigenvalues, we then discuss the behavior of the chiral symmetry breaking term as a function of the fifth dimensional length. We finally argue the chiral property of the DWF formulation in the limit of the infinite fifth dimensional length, in connection with spectra of the hermitian Wilson-Dirac operator in the infinite volume limit as well as in the finite volume.Comment: Added a reference and modified the acknowledgmen

Electronic transport through ballistic chaotic cavities: reflection symmetry, direct processes, and symmetry breaking

We extend previous studies on transport through ballistic chaotic cavities with spatial left-right (LR) reflection symmetry to include the presence of direct processes. We first analyze fully LR-symmetric systems in the presence of direct processes and compare the distribution w(T) of the transmission coefficient T with that for an asymmetric cavity with the same "optical" S matrix. We then study the problem of "external mixing" of the symmetry caused by an asymmetric coupling of the cavity to the outside. We first consider the case where symmetry breaking arises because two symmetrically positioned waveguides are coupled to the cavity by means of asymmetric tunnel barriers. Although this system is asymmetric with respect to the LR operation, it has a striking memory of the symmetry of the cavity it was constructed from. Secondly, we break LR symmetry in the absence of direct proceses by asymmetrically positioning the two waveguides and compare the results with those for the completely asymmetric case.Comment: 15 pages, 8 Postscript figures, submitted to Phys. Rev.

Moments of Structure Functions in Full QCD

Moments of the quark density distribution, moments of the quark helicity distribution, and the tensor charge are calculated in full QCD. Calculations of matrix elements of operators from the operator product expansion have been performed on $16^3 \times 32$ lattices for Wilson fermions at $\beta = 5.6$ using configurations from the SESAM collaboration and at $\beta = 5.5$ using configurations from SCRI. One-loop perturbative renormalization corrections are included. Selected results are compared with corresponding quenched calculations and with calculations using cooled configurations.Comment: Lattice 2000 (Hadronic Matrix Elements), 4 pages, 5 figure

Novel gonadal characteristics in an aged bovine freemartin

The gonads from a five year old freemartin Holstein animal were subjected to morphological analysis and to immunohistochemistry using antibodies against developmental and functional markers. We demonstrate, for the first time, the retention of anti-mullerian hormone (AMH) producing intratubular cells (Sertoli cells) in the context of abundant steroidogenic interstitial cells, and structures consistent with clusters of luteal cells. This novel report describes the clinical, gross and histological findings accompanying this newly described gonadal immunophenotype, and its implication in the understanding of freemartin development

On topological charge carried by nexuses and center vortices

In this paper we further explore the question of topological charge in the center vortex-nexus picture of gauge theories. Generally, this charge is locally fractionalized in units of 1/N for gauge group SU(N), but globally quantized in integral units. We show explicitly that in d=4 global topological charge is a linkage number of the closed two-surface of a center vortex with a nexus world line, and relate this linkage to the Hopf fibration, with homotopy $\Pi_3(S^3)\simeq Z$; this homotopy insures integrality of the global topological charge. We show that a standard nexus form used earlier, when linked to a center vortex, gives rise naturally to a homotopy $\Pi_2(S^2)\simeq Z$, a homotopy usually associated with 't Hooft-Polyakov monopoles and similar objects which exist by virtue of the presence of an adjoint scalar field which gives rise to spontaneous symmetry breaking. We show that certain integrals related to monopole or topological charge in gauge theories with adjoint scalars also appear in the center vortex-nexus picture, but with a different physical interpretation. We find a new type of nexus which can carry topological charge by linking to vortices or carry d=3 Chern-Simons number without center vortices present; the Chern-Simons number is connected with twisting and writhing of field lines, as the author had suggested earlier. In general, no topological charge in d=4 arises from these specific static configurations, since the charge is the difference of two (equal) Chern-Simons number, but it can arise through dynamic reconnection processes. We complete earlier vortex-nexus work to show explicitly how to express globally-integral topological charge as composed of essentially independent units of charge 1/N.Comment: Revtex4; 3 .eps figures; 18 page

Proposal for the numerical solution of planar QCD

Using quenched reduction, we propose a method for the numerical calculation of meson correlation functions in the planar limit of QCD. General features of the approach are outlined, and an example is given in the context of two-dimensional QCD.Comment: 31 pages, 10 figures, uses axodraw.sty, To appear in Physical Review

Monopole characteristics in various Abelian gauges

Renormalization group (RG) smoothing is employed on the lattice to investigate and to compare the monopole structure of the SU(2) vacuum as seen in different gauges (maximally Abelian (MAG), Polyakov loop (PG) and Laplacian gauge (LG)). Physically relevant types of monopoles (LG and MAG) are distinguished by their behavior near the deconfining phase transition. For the LG, Abelian projection reproduces well the gauge independent monopole structure encoded in an auxiliary Higgs field. Density and localization properties of monopoles, their non-Abelian action and topological charge are studied. Results are presented confirming the Abelian dominance with respect to the non-perturbative static potential for all gauges considered.Comment: 36 pages, 12 figure

Confinement and Chiral Symmetry Breaking via Domain-Like Structures in the QCD Vacuum

A qualitative mechanism for the emergence of domain structured background gluon fields due to singularities in gauge field configurations is considered, and a model displaying a type of mean field approximation to the QCD partition function based on this mechanism is formulated. Estimation of the vacuum parameters (gluon condensate, topological susceptibility, string constant and quark condensate) indicates that domain-like structures lead to an area law for the Wilson loop, nonzero topological susceptibility and spontaneous breakdown of chiral symmetry. Gluon and ghost propagators in the presence of domains are calculated explicitly and their analytical properties are discussed. The Fourier transforms of the propagators are entire functions and thus describe confined dynamical fields.Comment: RevTeX, 48 pages (32 pages + Appendices A-E), new references added [1,2,4,5] and minor formulae corrected for typographical error