Light hadrons with improved staggered quarks: approaching the continuum limit

We have extended our program of QCD simulations with an improved Kogut-Susskind quark action to a smaller lattice spacing, approximately 0.09 fm. Also, the simulations with a approximately 0.12 fm have been extended to smaller quark masses. In this paper we describe the new simulations and computations of the static quark potential and light hadron spectrum. These results give information about the remaining dependences on the lattice spacing. We examine the dependence of computed quantities on the spatial size of the lattice, on the numerical precision in the computations, and on the step size used in the numerical integrations. We examine the effects of autocorrelations in "simulation time" on the potential and spectrum. We see effects of decays, or coupling to two-meson states, in the 0++, 1+, and 0- meson propagators, and we make a preliminary mass computation for a radially excited 0- meson.Comment: 43 pages, 16 figure

Modeling two-language competition dynamics

During the last decade, much attention has been paid to language competition in the complex systems community, that is, how the fractions of speakers of several competing languages evolve in time. In this paper we review recent advances in this direction and focus on three aspects. First we consider the shift from two-state models to three state models that include the possibility of bilingual individuals. The understanding of the role played by bilingualism is essential in sociolinguistics. In particular, the question addressed is whether bilingualism facilitates the coexistence of languages. Second, we will analyze the effect of social interaction networks and physical barriers. Finally, we will show how to analyze the issue of bilingualism from a game theoretical perspective.Comment: 15 pages, 5 figures; published in the Special Issue of Advances in Complex Systems "Language Dynamics

A length scale for the superconducting Nernst signal above Tc_{c} in Nb0.15_{0.15}Si0.85_{0.85}

We present a study of the Nernst effect in amorphous superconducting thin films of Nb$_{0.15}$Si$_{0.85}$. The field dependence of the Nernst coefficient above T$_{c}$ displays two distinct regimes separated by a field scale set by the Ginzburg-Landau correlation length. A single function $F(\xi)$, with the correlation length as its unique argument set either by the zero-field correlation length (in the low magnetic field limit) or by the magnetic length (in the opposite limit), describes the Nernst coefficient. We conclude that the Nernst signal observed on a wide temperature ($30 \times T_c$) and field ($4 \times B_{c2}$) range is exclusively generated by short-lived Cooper pairs.Comment: 4 pages, 4 figure

The scaling dimension of low lying Dirac eigenmodes and of the topological charge density

As a quantitative measure of localization, the inverse participation ratio of low lying Dirac eigenmodes and topological charge density is calculated on quenched lattices over a wide range of lattice spacings and volumes. Since different topological objects (instantons, vortices, monopoles, and artifacts) have different co-dimension, scaling analysis provides information on the amount of each present and their correlation with the localization of low lying eigenmodes.Comment: Lattice2004(topology), Fermilab, June 21 - 26, 2004; 3 pages, 3 figure

High Energy Physics from High Performance Computing

We discuss Quantum Chromodynamics calculations using the lattice regulator. The theory of the strong force is a cornerstone of the Standard Model of particle physics. We present USQCD collaboration results obtained on Argonne National Lab's Intrepid supercomputer that deepen our understanding of these fundamental theories of Nature and provide critical support to frontier particle physics experiments and phenomenology.Comment: Proceedings of invited plenary talk given at SciDAC 2009, San Diego, June 14-18, 2009, on behalf of the USQCD collaboratio

Inner approximated reachability analysis

International audienceComputing a tight inner approximation of the range of a function over some set is notoriously di cult, way beyond obtaining outer approximations. We propose here a new method to compute a tight inner approximation of the set of reachable states of non-linear dynamical systems on a bounded time interval. This approach involves a ne forms and Kaucher arithmetic, plus a number of extra ingredients from set-based methods. An implementation of the method is discussed, and illustrated on representative numerical schemes, discrete-time and continuous-time dynamical systems

Sharp constants in weighted trace inequalities on Riemannian manifolds

We establish some sharp weighted trace inequalities W^{1,2}(\rho^{1-2\sigma}, M)\hookrightarrow L^{\frac{2n}{n-2\sigma}}(\pa M) on $n+1$ dimensional compact smooth manifolds with smooth boundaries, where $\rho$ is a defining function of $M$ and $\sigma\in (0,1)$. This is stimulated by some recent work on fractional (conformal) Laplacians and related problems in conformal geometry, and also motivated by a conjecture of Aubin.Comment: 34 page

Light pseudoscalar decay constants, quark masses, and low energy constants from three-flavor lattice QCD

As part of our program of lattice simulations of three flavor QCD with improved staggered quarks, we have calculated pseudoscalar meson masses and decay constants for a range of valence quark masses and sea quark masses on lattices with lattice spacings of about 0.125 fm and 0.09 fm. We fit the lattice data to forms computed with staggered chiral perturbation theory. Our results provide a sensitive test of the lattice simulations, and especially of the chiral behavior, including the effects of chiral logarithms. We find: f_\pi=129.5(0.9)(3.5)MeV, f_K=156.6(1.0)(3.6)MeV, and f_K/f_\pi=1.210(4)(13), where the errors are statistical and systematic. Following a recent paper by Marciano, our value of f_K/f_\pi implies |V_{us}|=0.2219(26). Further, we obtain m_u/m_d= 0.43(0)(1)(8), where the errors are from statistics, simulation systematics, and electromagnetic effects, respectively. The data can also be used to determine several of the constants of the low energy effective Lagrangian: in particular we find 2L_8-L_5=-0.2(1)(2) 10^{-3} at chiral scale m_\eta. This provides an alternative (though not independent) way of estimating m_u; 2L_8-L_5 is far outside the range that would allow m_u=0. Results for m_s^\msbar, \hat m^\msbar, and m_s/\hat m can be obtained from the same lattice data and chiral fits, and have been presented previously in joint work with the HPQCD and UKQCD collaborations. Using the perturbative mass renormalization reported in that work, we obtain m_u^\msbar=1.7(0)(1)(2)(2)MeV and m_d^\msbar=3.9(0)(1)(4)(2)MeV at scale 2 GeV, with errors from statistics, simulation, perturbation theory, and electromagnetic effects, respectively.Comment: 86 pages, 22 figures. v3: Remarks about m_u=0 and the strong CP problem modified; reference added. Figs 5--8 modified for clarity. Version to be published in Phys. Rev. D. v2: Expanded discussion of finite volume effects, normalization in Table I fixed, typos and minor errors correcte