Gluon flux-tube distribution and linear confinement in baryons

We have observed the formation of gluon flux-tubes within baryons using lattice QCD techniques. A high-statistics approach, based on translational and rotational symmetries of the four-dimensional lattice, enables us to observe correlations between vacuum action density and quark positions in a completely gauge independent manner. This contrasts with earlier studies which used gauge-dependent smoothing techniques. We used 200 O(a^2) improved quenched QCD gauge-field configurations on a 16^3x32 lattice with a lattice spacing of 0.123 fm. In the presence of static quarks flux tubes representing the suppression of gluon-field fluctuations are observed. We have analyzed 11 L-shapes and 8 T and Y shapes of varying sizes in order to explore a variety of flux-tube topologies, including the ground state. At large separations, Y-shape flux-tube formation is observed. T-shaped paths are observed to relax towards a Y-shaped topology, whereas L-shaped paths give rise to a large potential energy. We do not find any evidence for the formation of a Delta-shaped flux-tube (empty triangle) distribution. However, at small quark separations, we observe an expulsion of gluon-field fluctuations in the shape of a filled triangle with maximal expulsion at the centre of the triangle. Having identified the precise geometry of the flux distribution, we are able to perform quantitative comparison between the length of the flux-tube and the associated static quark potential. For every source configuration considered we find a universal string tension, and conclude that, for large quark separations, the ground state potential is that which minimizes the length of the flux-tube. The flux tube radius of the baryonic ground state potential is found to be 0.38 \pm 0.03 fm, with vacuum fluctuations suppressed by 7.2 \pm 0.6 %.Comment: 16 pages, final version as accepted for publication in Physical review D1. Abstract, text, references and some figures have been revise

Finding the Minimum-Weight k-Path

Given a weighted $n$-vertex graph $G$ with integer edge-weights taken from a range $[-M,M]$, we show that the minimum-weight simple path visiting $k$ vertices can be found in time \tilde{O}(2^k \poly(k) M n^\omega) = O^*(2^k M). If the weights are reals in $[1,M]$, we provide a $(1+\varepsilon)$-approximation which has a running time of \tilde{O}(2^k \poly(k) n^\omega(\log\log M + 1/\varepsilon)). For the more general problem of $k$-tree, in which we wish to find a minimum-weight copy of a $k$-node tree $T$ in a given weighted graph $G$, under the same restrictions on edge weights respectively, we give an exact solution of running time \tilde{O}(2^k \poly(k) M n^3) and a $(1+\varepsilon)$-approximate solution of running time \tilde{O}(2^k \poly(k) n^3(\log\log M + 1/\varepsilon)). All of the above algorithms are randomized with a polynomially-small error probability.Comment: To appear at WADS 201

Nonlinear cellular instabilities of planar premixed flames: numerical simulations of the Reactive Navier-Stokes equations

Two-dimensional compressible Reactive Navier-Stokes numerical simulations of intrinsic planar, premixed flame instabilities are performed. The initial growth of a sinusoidally perturbed planar flame is first compared with the predictions of a recent exact linear stability analysis, and it is shown the analysis provides a necessary but not sufficient test problem for validating numerical schemes intended for flame simulations. The long-time nonlinear evolution up to the final nonlinear stationary cellular flame is then examined for numerical domains of increasing width. It is shown that for routinely computationally affordable domain widths, the evolution and final state is, in general, entirely dependent on the width of the domain and choice of numerical boundary conditions. It is also shown that the linear analysis has no relevance to the final nonlinear cell size. When both hydrodynamic and thermal-diffusive effects are important, the evolution consists of a number of symmetry breaking cell splitting and re-merging processes which results in a stationary state of a single very asymmetric cell in the domain, a flame shape which is not predicted by weakly nonlinear evolution equations. Resolution studies are performed and it is found that lower numerical resolutions, typical of those used in previous works, do not give even the qualitatively correct solution in wide domains. We also show that the long-time evolution, including whether or not a stationary state is ever achieved, depends on the choice of the numerical boundary conditions at the inflow and outflow boundaries, and on the numerical domain length and flame Mach number for the types of boundary conditions used in some previous works

Evaporation and growth of crystals - propagation of step density compression waves at vicinal surfaces

We studied the step dynamics during crystal sublimation and growth in the limit of fast surface diffusion and slow kinetics of atom attachment-detachment at the steps. For this limit we formulate a model free of the quasi-static approximation in the calculation of the adatom concentration on the terraces at the crystal surface. Such a model provides a relatively simple way to study the linear stability of a step train in a presence of step-step repulsion and an absence of destabilizing factors (as Schwoebel effect, surface electromigration etc.). The central result is that a critical velocity of the steps in the train exists which separates the stability and instability regimes. When the step velocity exceeds its critical value the plot of these trajectories manifests clear space and time periodicity (step density compression waves propagate on the vicinal surface). This ordered motion of the steps is preceded by a relatively short transition period of disordered step dynamics.Comment: 18 pages, 6 figure

Glueball enhancements in p(gamma,VV)p through vector meson dominance

Double vector meson photoproduction, p(gamma, G -> VV)p, mediated by a scalar glueball G is investigated. Using vector meson dominance (VMD) and Regge/pomeron phenomenology, a measureable glueball enhancement is predicted in the invariant VV = rho rho and omega omega mass spectra. The scalar glueball is assumed to be the lightest physical state on the daughter pomeron trajectory governing diffractive vector meson photoproduction. In addition to cross sections, calculations for hadronic and electromagnetic glueball decays, G -> V V' (V,V'= rho, omega, phi, gamma), and gamma_v V -> G transition form factors are presented based upon flavor universality, VMD and phenomenological couplings from phi photoproduction analyses. The predicted glueball decay widths are similar to an independent theoretical study. A novel signature for glueball detection is also discussed

Carrier drift velocity and edge magnetoplasmons in graphene

We investigate electron dynamics at the graphene edge by studying the propagation of collective edge magnetoplasmon (EMP) excitations. By timing the travel of narrow wave-packets on picosecond time scales around exfoliated samples, we find chiral propagation with low attenuation at a velocity which is quantized on Hall plateaus. We extract the carrier drift contribution from the EMP propagation and find it to be slightly less than the Fermi velocity, as expected for an abrupt edge. We also extract the characteristic length for Coulomb interaction at the edge and find it to be smaller than for soft, depletion edge systems.Comment: 5 pages, 3 figures of main text and 6 pages, 6 figures of supplemental materia

Potential and current distribution in strongly anisotropic Bi(2)Sr(2) CaCu(2)O(8) single crystals at current breakdown

Experiments on potential differences in the low-temperature vortex solid phase of monocrystalline platelets of superconducting Bi(2)Sr(2)CaCu(2)O(8) (BSCCO) subjected to currents driven either through an "ab" surface or from one such surface to another show evidence of a resistive/nonresistive front moving progressively out from the current contacts as the current increases. The depth of the resistive region has been measured by a novel in-depth voltage probe contact. The position of the front associated with an injection point appears to depend only on the current magnitude and not on its withdrawal point. It is argued that enhanced nonresistive superconducting anisotropy limits current penetration to less than the London length and results in a flat rectangular resistive region with simultaneous "ab" and "c" current breakdown which moves progressively out from the injection point with increasing current. Measurements in "ab" or "c" configurations are seen to give the same information, involving both ab-plane and c-axis conduction properties.Comment: 9 pages, 13 figures, typo error corrected, last section was refine