Source author record

Pedro Bicudo

Pedro Bicudo appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

56works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

56 published item(s)

preprint2022arXiv

Bottomonium resonances from lattice QCD static-static-light-light potentials

We study $I=0$ quarkonium resonances decaying into pairs of heavy-light mesons using static-static-light-light potentials from lattice QCD. To this end, we solve a coupled channel Schrödinger equation with a confined quarkonium channel and channels with a heavy-light meson pair to compute phase shifts and $\mbox{T}$ matrix poles for the lightest decay channel. We discuss our results for $S$, $P$, $D$ and $F$ wave states in the context of corresponding experimental results, in particular for $Υ(10753)$ and $Υ(10860)$.

preprint2022arXiv

Computation of the quarkonium and meson-meson composition of the $Υ(nS)$ states and of the new $Υ(10753)$ Belle resonance from lattice QCD static potentials

We compute the composition of the bottomonium $Υ(nS)$ states (including $Υ(10860)$) and the new $Υ(10753)$ resonance reported by Belle in terms of quarkonium and meson-meson components. We use a recently developed novel approach utilizing lattice QCD string breaking potentials for the study of resonances. This approach is based on the Born Oppenheimer approximation and the unitary emergent wave method and allows to compute the poles of the $\mbox{S}$ matrix. We focus on $I=0$ bottomonium $S$ wave bound states and resonances, where the Schrödinger equation is a set of coupled differential equations. One of the channels corresponds to a confined heavy quark-antiquark pair $\bar b b$, the others to pairs of heavy-light mesons. In a previous study only one meson-meson channel $\bar{B}^{(\ast)} B^{(\ast)}$ was considered. Now we also include the closed strangeness channel $\bar{B}_s^{(\ast)} B_s^{(\ast)}$ extending our formalism significantly to have a more realistic description of bottomonium. We confirm the new Belle resonance $Υ(10753)$ as a dynamical meson-meson resonance with around $85 \%$ meson-meson content. Moreover, we identify $Υ(4S)$ and $Υ(10860)$ as states with both sizable quarkonium and meson-meson contents. With these results we contribute to the clarification of ongoing controversies in the vector bottomonium spectrum.

preprint2021arXiv

Highly excited pure gauge SU(3) flux tubes

Flux tube spectra are expected to have full towers of levels due to the quantization of the string vibrations. We study a spectrum of flux tubes with static quark and antiquark sources with pure gauge $SU(3)$ lattice QCD in 3+1 dimensions up to a significant number of excitations. To go high in the spectrum, we specialize in the most symmetric case $Σ_g^+$, use a large set of operators, solve the generalized eigenvalue and compare different lattice QCD gauge actions and anisotropies.

preprint2020arXiv

Bottomonium resonances with $I = 0$ from lattice QCD correlation functions with static and light quarks

We discuss, how to study $I = 0$ quarkonium resonances decaying into pairs of heavy-light mesons using static potentials from lattice QCD. These static potentials can be obtained from a set of correlation functions containing both static and light quarks. As a proof of concept we focus on bottomonium with relative orbital angular momentum $L = 0$ of the $\bar{b} b$ pair corresponding to $J^{P C} = 0^{- +}$ and $J^{P C} = 1^{- -}$. We use static potentials from an existing lattice QCD string breaking study and compute phase shifts and $\mbox{T}$ matrix poles for the lightest heavy-light meson-meson decay channel. We discuss our results in the context of corresponding experimental results, in particular for $Υ(10860)$ and $Υ(11020)$.

preprint2020arXiv

From string breaking to quarkonium spectrum

We present a preliminary computation of potentials between two static quarks in $n_f=2$ QCD with O(a) improved Wilson fermions. We explore different smearing choices (HYP, HYP2 and APE) and their effect on the signal to noise ratio in the computed static potentials. This is a part of a larger effort concerning, at first, a precise computation of the QCD string breaking parameters and their subsequent utilization for the recent approach based on Born-Oppenheimer approximation (Bicudo et al. 2020 \cite{Bicudo:2019ymo}) to study quarkonium resonances and bound states.

preprint2016arXiv

Flux tubes at Finite Temperature

We show the flux tubes produced by static quark-antiquark, quark-quark and quark-gluon charges at finite temperature. The sources are placed in the lattice with fundamental and adjoint Polyakov loops. We compute the square densities of the chromomagnetic and chromoelectric fields above and below the phase transition. Our results are gauge invariant and produced in pure gauge SU(3). The codes are written in CUDA and the computations are performed with GPUs.

preprint2016arXiv

Including heavy spin effects in a lattice QCD study of static-static-light-light tetraquarks

In previous works we predicted the existence of a $\bar b \bar b u d$ tetraquark with quantum numbers $I(J^P) = 0(1^+)$ using the static approximation for the $\bar b$ quarks and neglecting heavy spin effects. Since the binding energy is of the same order as expected for these heavy spin effects, it is essential to include them in the computation. Here we present a corresponding method and show evidence that binding is only slightly weakened and that the $\bar b \bar b u d$ tetraquark persists.

preprint2016arXiv

Investigation of $B\bar B$ four-quark systems using lattice QCD

We investigate $B \bar B$ systems by computing potentials of two static quarks in the presence of two quarks of finite mass using lattice QCD. By solving the Schrödinger equation we check whether these potentials are sufficiently attractive to host bound states. Particular focus is put on the experimentally most promising bottomonium-like tetraquark candidate $Z_b^\pm$ with quantum numbers $I(J^P)=1(1^+)$.

preprint2016arXiv

Lattice QCD study of heavy-heavy-light-light tetraquark candidates

We investigate heavy-light four-quark systems $ud\bar b \bar b$ with bottom quarks of finite mass which are treated in the framework of NRQCD. We focus on $I(J^P)=0(1^+)$, where we recently found evidence for the existence of a tetraquark state using static bottom quarks. Furthermore, we report on an investigation of the $u \bar d b \bar b$ four-quark system with quantum numbers $I(J^P)=1(1^+)$ again using static bottom quarks.

preprint2016arXiv

Mapping chiral symmetry breaking in the excited baryon spectrum

We study the conjectured "Insensitivity to Chiral Symmetry Breaking" in the highly excited light baryon spectrum. While the experimental spectrum is being measured at JLab and CBELSA/TAPS, this insensitivity remains to be computed theoretically in detail. As the only existing option to have both confinement, highly excited states and chiral symmetry, we adopt the truncated Coulomb gauge formulation of QCD, considering a linearly confining Coulomb term. Adopting a systematic and numerically intensive variational treatment up to 12 harmonic oscillator shells we are able to access several angular and radial excitations. We compute both the excited spectra of $I=1/2$ and $I=3/2$ baryons, up to large spin $J=13/2$, and study in detail the proposed chiral multiplets. While the static-light and light-light spectra clearly show chiral symmetry restoration high in the spectrum, the realization of chiral symmetry is more complicated in the baryon spectrum than earlier expected.

preprint2016arXiv

Tighter bounds on a hypothetical graviton screening mass from the gravitational wave observation GW150914 at LIGO

While quantum gravity is not solved yet, a screening mass for the graviton remains theoretically possible. If such a mass would screen gravity at distances of the order of the cluster galaxy radius, it could account for the universe expansion. The modified Newtonian dynamics model also could be related to a screening graviton mass at inter-galactic scales. Moreover, massive spin-2 theories constitute a very active theoretical topic. We briefly show how the very recent LIGO gravitational wave observation GW150914, emitted by a binary black hole merger distant $\sim 1.3 \times 10^9$ ly from the Earth, tightens the phenomenological bound on a massive graviton or on the screening of gravity.

preprint2015arXiv

BB interactions with static bottom quarks from Lattice QCD

The isospin, spin and parity dependent potential of a pair of $B$ mesons is computed using Wilson twisted mass lattice QCD with two flavours of degenerate dynamical quarks. The $B$ meson is addressed in the static-light approximation, i.e.\ the $b$ quarks are infinitely heavy. From the results of the $B\,B$ meson-meson potentials, a simple rule can be deduced stating which isospin, spin and parity combinations correspond to attractive and which to repulsive forces. We provide fits to the ground state potentials in the attractive channels and discuss the potentials in the repulsive and excited channels. The attractive channels are most important since they can possibly lead to a bound four-quark state, i.e.\ a $\bar{b}\bar{b}ud$ tetraquark. Using these attractive potentials in the Schrödinger equation, we find indication for such a tetraquark state of two static bottom antiquarks and two light $u/d$ quarks with mass extrapolated down to the physical value.

preprint2015arXiv

Evidence for the existence of $u d \bar{b} \bar{b}$ and the non-existence of $s s \bar{b} \bar{b}$ and $c c \bar{b} \bar{b}$ tetraquarks from lattice QCD

We combine lattice QCD results for the potential of two static antiquarks in the presence of two quarks $q q$ of finite mass and quark model techniques to study possibly existing $q q \bar{b} \bar{b}$ tetraquarks. While there is strong indication for a bound four-quark state for $q q = (ud-du) / \sqrt{2}$, i.e. isospin $I=0$, we find clear evidence against the existence of corresponding tetraquarks with $q q \in \{ uu , (ud+du) / \sqrt{2} , dd \}$, i.e. isospin $I=1$, $q q = s s$ and $q q = c c$.

preprint2015arXiv

Exploring possibly existing $q q \bar b \bar b$ tetraquark states with $q q = ud, ss, cc$

We compute potentials of two static antiquarks in the presence of two quarks $qq$ of finite mass using lattice QCD. In a second step we solve the Schrödinger equation, to determine, whether the resulting potentials are sufficiently attractive to host a bound state, which would indicate the existence of a stable $q q \bar b \bar b$ tetraquark. We find a bound state for $qq=(ud-du)/\sqrt{2}$ with corresponding quantum numbers $I(J^ P)=0(1^+)$ and evidence against the existence of bound states with isospin $I=1$ or $qq \in \{cc,ss \}$

preprint2015arXiv

Gauge fixing and the gluon propagator in renormalizable xi gauges

Covariant $R_ξ$ gauge fixing is notoriously difficult for large lattice volumes, large $ξ$ and small $N_c$. We thoroughly test different convergence techniques, which allows the gauge fixing of lattice configurations with a total volume of (3.25 fm)$^4$, up to $ξ=0.5$. We are able to study the gluon propagator in the infrared region and its dependence on the gauge fixing parameter $ξ$. As expected, the longitudinal gluon dressing functions stay constant at their tree-level value $ξ$. Similar to the Landau gauge, the transverse $R_ξ$ gauge gluon propagators saturate at a non-vanishing value in the deep infrared for all values of $ξ$ studied. We compare with very recent continuum studies and perform a simple analysis of the found saturation with a dynamically generated effective gluon mass.

preprint2015arXiv

Refined lattice/model investigation of $u d \bar{b} \bar{b}$ tetraquark candidates with heavy spin effects taken into account

We investigate four-quark systems consisting of two heavy anti-bottom quarks and two light up/down quarks. We propose to solve a coupled Schrödinger equation for the anti-bottom-anti-bottom separation using potentials computed via lattice QCD in the limit of static anti-bottom quarks. This coupled Schrödinger equation allows to incorporate effects due to the heavy anti-bottom spins. First exploratory numerical tests are discussed.

preprint2014arXiv

An unitarized model for tetraquarks with a color flip-flip potential

In this work, a color structure dependent flip-flop potential is developed for the two quarks and two antiquarks system. Then, this potential is applied to a microscopic quark model which, by integrating the internal degrees of freedom, is transformed into a model of mesons with non-local interactions. With this, the T matrix for the system is constructed and meson-meson scattering is studied. Tetraquarks states, interpreted as poles of of the T matrix, both bound states and resonances, are found. Special emphasis is given to the truly exotic $qq\bar{Q}\bar{Q}$ system, but some results for the crypto-exotic $qQ\bar{q}\bar{Q}$ are also presented.

preprint2014arXiv

Aspects of gluon propagation in Landau gauge: spectral densities, and mass scales at finite temperature

We discuss a method to extract the Källén-Lehmann spectral density of a particle (be it elementary or bound state) propagator and apply it to compute gluon spectral densities from lattice data. Furthermore, we also consider the interpretation of the Landau-gauge gluon propagator at finite temperature as a massive-type bosonic propagator.

preprint2014arXiv

Many faces of the Landau gauge gluon propagator at zero and finite temperature: positivity violation, spectral density and mass scales

We address several aspects of gluon propagation at zero and finite temperature. In particular, we study the violation of spectral positivity, we discuss a method to extract the Källén-Lehmann spectral density of a particle (be it elementary or bound state) propagator and apply it to compute gluon spectral densities from lattice data. Furthermore, we also consider the interpretation of the Landau gauge gluon propagator at finite temperature as a massive type bosonic propagator.

preprint2014arXiv

Matrix model for deconfinement in a SU(Nc) gauge theory in 2+1 dimensions

We use an effective matrix model to study deconfinement in a pure SU(Nc) gauge theory, without quarks, in d=2+1 dimensions. Expanding about a constant background A0 field we construct an effective potential for the eigenvalues of the thermal Wilson line for general Nc and in the large-Nc limit. The numerical results are presented using one, two, and four free parameters, which are determined by fitting directly to the lattice data for the pressure. The matrix model shows a good agreement with numerical lattice simulations for the pressure and the interaction measure, starting from the perturbative limit up to the critical temperature. For the pressure, the details of A0- dependent nonperturbative terms are relevant only in a narrow transition region below ~1.2 Td. This is also the range where the Polyakov loop deviates notably from one. In accordance with the lattice results we find that, up to a trivial factor Nc^2-1, there is only a mild dependence on the number of colors.

preprint2013arXiv

Lattice QCD signal for a bottom-bottom tetraquark

Utilizing lattice QCD results for the potential of two static antiquarks and two dynamical quarks as well as quark model techniques for the dynamics of two heavy antiquarks in a cloud of two light quarks, we are provided with an accurate framework for the study of possibly existing heavy-heavy-light-light tetraquarks. Among the possible quantum numbers of such a system, we find binding in only one channel, the scalar isosinglet. Solving the Schroedinger equation for the displacement of the heavy antiquarks and taking systematic errors into account, we find an antibottom-antibottom-light-light bound state with a confidence level of around 1.8 sigma ... 3.0 sigma and binding energy of approximately 30 MeV ... 57 MeV.

preprint2013arXiv

Matrix model for deconfinement in an SU(2) gauge theory in 2+1 dimensions

We use matrix models to characterize deconfinement at a nonzero temperature T for an SU(2) gauge theory in three spacetime dimensions. At one loop order, the potential for a constant vector potential A0 is ~T^3 times a trilogarithm function of A0/T. In addition, we add various nonperturbative terms to model deconfinement. The parameters of the model are adjusted by fitting the lattice results for the pressure. The nonperturbative terms are dominated by a constant term ~T^2Td, where Td is the temperature for deconfinement. Besides this constant, we add terms which are nontrivial functions of A0/T, both ~T^2Td and ~TTd^2. There is only a mild sensitivity to the details of these nonconstant terms. Overall we find a good agreement with the lattice results. For the pressure, the conformal anomaly, and the Polyakov loop the nonconstant terms are relevant only in a narrow region below ~1.2Td. We also compute the 't Hooft loop, and find that the details of the nonconstant terms enter in a much wider region, up to ~4Td.

preprint2012arXiv

Colour Fields of the Static Pentaquark System Computed in SU(3) Lattice QCD

We compute the colour fields of SU(3) lattice QCD created by static pentaquark systems, in a $24^3\times 48$ lattice at $β=6.2$ corresponding to a lattice spacing $a=0.07261(85)$ fm. We find that the pentaquark colour fields are well described by a multi-Y-type shaped flux tube. The flux tube junction points are compatible with Fermat-Steiner points minimizing the total flux tube length. We also compare the pentaquark flux tube profile with diquark-diantiquark central flux tube profile in the tetraquark and the quark-antiquark fundamental flux tube profile in the meson, and they match, thus showing that the pentaquark flux tubes are composed of fundamental flux tubes.

preprint2012arXiv

Colour flux-tubes in static Pentaquark and Tetraquark systems

The colour fields created by the static tetraquark and pentaquark systems are computed in quenched SU(3) lattice QCD, with gauge invariant lattice operators, in a 24^3 x 48 lattice at beta=6.2. We generate our quenched configurations with GPUs, and detail the respective benchmanrks in different SU(N) groups. While at smaller distances the coulomb potential is expected to dominate, at larger distances it is expected that fundamental flux tubes, similar to the flux-tube between a quark and an antiquark, emerge and confine the quarks. In order to minimize the potential the fundamental flux tubes should connect at 120o angles. We compute the square of the colour fields utilizing plaquettes, and locate the static sources with generalized Wilson loops and with APE smearing. The tetraquark system is well described by a double-Y-shaped flux-tube, with two Steiner points, but when quark-antiquark pairs are close enough the two junctions collapse and we have an X-shaped flux-tube, with one Steiner point. The pentaquark system is well described by a three-Y-shaped flux-tube where the three flux the junctions are Steiner points.

preprint2012arXiv

CuBA - a CUDA implementation of BAMPS

Using CUDA as programming language, we create a code named CuBA which is based on the CPU code "Boltzmann Approach for Many Parton Scattering (BAMPS)" developed in Frankfurt in order to study a system of many colliding particles resulting from heavy ion collisions. Furthermore, we benchmark our code with the Riemann Problem and compare the results with BAMPS. They demonstrate an improvement of the computational runtime, by one order of magitude.

preprint2012arXiv

Detailed study of quark-antiquark flux tubes and flux tube recombination

In this work we compute the color fields in the mediator plane between a static quark and a static antiquark using quenched lattice QCD. In special, we see the effect of the quark-antiquark distance on the flux tube. To obtain this results an improved multihit technique is developed and an extend smearing technique is used. Then, we also discuss the flux-tubes in a system composed of two quarks and two antiquarks. The ground and first excited states fields are studied for different dispositions of the system.

preprint2012arXiv

Generating SU(Nc) pure gauge lattice QCD configurations on GPUs with CUDA

The starting point of any lattice QCD computation is the generation of a Markov chain of gauge field configurations. Due to the large number of lattice links and due to the matrix multiplications, generating SU(Nc) lattice QCD configurations is a highly demanding computational task, requiring advanced computer parallel architectures such as clusters of several Central Processing Units (CPUs) or Graphics Processing Units (GPUs). In this paper we present and explore the performance of CUDA codes for NVIDIA GPUs to generate SU(Nc) lattice QCD pure gauge configurations. Our implementation in one GPU uses CUDA and in multiple GPUs uses OpenMP and CUDA. We present optimized CUDA codes SU(2), SU(3) and SU(4). We also show a generic SU(Nc) code for Nc$\,\geq 4$ and compare it with the optimized version of SU(4). Our codes are publicly available for free use by the lattice QCD community.

preprint2012arXiv

Landau Gauge Fixing on GPUs and String Tension

We explore the performance of CUDA in performing Landau gauge fixing in Lattice QCD, using the steepest descent method with Fourier acceleration. The code performance was tested in a Tesla C2070, Fermi architecture. We also present a study of the string tension at finite temperature in the confined phase. The string tension is extracted from the color averaged free energy and from the color singlet using Landau gauge fixing.

preprint2012arXiv

Study of compact U(1) flux tubes in 3+1 dimensions in lattice gauge theory using GPU's

We utilize Polyakov loop correlations to study (3+1)D compact U(1) flux tubes and the static electron-positron potential in lattice gauge theory. By using field operators it is possible in U(1) lattice gauge theory to probe directly the electric and magnetic fields. In order to improve the signal-to-noise ratio in the confinement phase, we apply the Lüscher-Weiss multilevel algorithm. Our code is written in CUDA, and we run it in NVIDIA FERMI generation GPU's, in order to achieve the necessary performance for our computations.

preprint2012arXiv

Tetraquark and the flux tube recombination

Here we study the static potential for the two quarks and two antiquarks system. First this is done using the tetraquark operator, which has been previously used to calculate the static potential. This is found to give good results in the region where the tetraquark is expected to be the ground state, however failing outside it. To repair this, we resort to a variational method. This let us study the first excited state besides the ground state of the system in two different particle dispositions, one where the quarks are on the same side of a rectangle and the other where they are at opposite sides. Results for the field components and for the lagrangian density are presented.

preprint2011arXiv

Colour flux-tubes in static Pentaquark systems

The colour fields created by the static tetraquark and pentaquark systems are computed in quenched SU(3) lattice QCD, with gauge invariant lattice operators, in a 24^3 x 48 lattice at beta=6.2 . We generate our quenched configurations with GPUs, and detail the respective benchmanrks in different SU(N) groups. While at smaller distances the coulomb potential is expected to dominate, at larger distances it is expected that fundamental flux tubes, similar to the flux-tube between a quark and an antiquark, emerge and confine the quarks. In order to minimize the potential the fundamental flux tubes should connect at 120o angles. We compute the square of the colour fields utilizing plaquettes, and locate the static sources with generalized Wilson loops and with APE smearing. The tetraquark system is well described by a double-Y-shaped flux-tube, with two Steiner points, but when quark-antiquark pairs are close enough the two junctions collapse and we have an X-shaped flux-tube, with one Steiner point. The pentaquark system is well described by a three-Y-shaped flux-tube where the three flux the junctions are Steiner points.

preprint2011arXiv

Excited charmonium states from Bethe-Salpeter equation

We solve the Bethe-Salpeter equation for a system of a heavy quark-antiquark pair interacting with a screened linear confining potential. First we show the spinless QFT model is inadequate and fail to describe even gross feature of the quarkonia spectrum. In order to get reliable description the spine degrees of freedom has to be considered. Within the approximation employed we reasonably reproduce known radial excitation of vector charmonium. The BSE favors relatively large string breaking scale $μ\simeq 350MeV$ . Using free charm quark propagators we observe that $J/Ψ$ is the only charmonium left bellow naive quark-antiquark threshold $2m_c$, while the all excited states are situated above this threshold. Within the numerical method we overcome obstacles related with threshold singularity and discuss the consequences of the use of free propagators for calculation of excited states above the threshold.

preprint2011arXiv

Finite temperature lattice QCD with GPUs

Graphics Processing Units (GPUs) are being used in many areas of physics, since the performance versus cost is very attractive. The GPUs can be addressed by CUDA which is a NVIDIA's parallel computing architecture. It enables dramatic increases in computing performance by harnessing the power of the GPU. We present a performance comparison between the GPU and CPU with single precision and double precision in generating lattice SU(2) configurations. Analyses with single and multiple GPUs, using CUDA and OPENMP, are also presented. We also present SU(2) results for the renormalized Polyakov loop, colour averaged free energy and the string tension as a function of the temperature.

preprint2011arXiv

String tension at finite temperature Lattice QCD

The critical curve of string tension as a function of the temperature is computed in SU(3) Lattice QCD. We present the results for the string tension utilizing a pair of Polyakov loop and antiloop, with two different techniques. We compare the colour averaged loop-antiloop which is gauge invariant but is only adequate to study the string tension, and the colour singlet loop-antiloop using the Landau gauge fixing which also enables to compute the coulomb part of the free energy.

preprint2011arXiv

SU(2) Lattice Gauge Theory Simulations on Fermi GPUs

In this work we explore the performance of CUDA in quenched lattice SU(2) simulations. CUDA, NVIDIA Compute Unified Device Architecture, is a hardware and software architecture developed by NVIDIA for computing on the GPU. We present an analysis and performance comparison between the GPU and CPU in single and double precision. Analyses with multiple GPUs and two different architectures (G200 and Fermi architectures) are also presented. In order to obtain a high performance, the code must be optimized for the GPU architecture, i.e., an implementation that exploits the memory hierarchy of the CUDA programming model. We produce codes for the Monte Carlo generation of SU(2) lattice gauge configurations, for the mean plaquette, for the Polyakov Loop at finite T and for the Wilson loop. We also present results for the potential using many configurations ($50\ 000$) without smearing and almost $2\ 000$ configurations with APE smearing. With two Fermi GPUs we have achieved an excellent performance of $200 \times$ the speed over one CPU, in single precision, around 110 Gflops/s. We also find that, using the Fermi architecture, double precision computations for the static quark-antiquark potential are not much slower (less than $2 \times$ slower) than single precision computations.

preprint2011arXiv

The chiral crossover, static-light and light-light meson spectra, and the deconfinement crossover

We study the chiral crossover, the spectra of light-light and of static-light mesons and the deconfinement crossover at finite temperature T. Our framework is the confining and chiral invariant quark model, related to truncated Coulomb gauge QCD. Since we are dealing with light quarks, where the linear potential dominates the quark condensate and the spectrum, we only specialize in the linear confining potential for the quark-antiquark interaction. We utilize T dependent string tensions previously fitted from lattice QCD data, and a fit of previously computed dynamically generated constituent quark masses. We scan the T effects on the constituent quark mass, on the meson spectra and on the polyakov loop.

preprint2010arXiv

Chromoelectric and chromomagnetic fields for the static gluon-quark-antiquark system

The chromoelectric and chromomagnetic fields, created by a static gluon-quark-antiquark system, are computed in the quenched approximation of lattice QCD, in a $24^3\times 48$ lattice at $β=6.2$. We study two geometries, one with a U shape and another with an L shape. The degenerate case of the two gluon glueball is also studied. This is relevant to understand the microscopic structure of hadrons, in particular of hybrids. This also contributes to understand confinement with flux tubes of the chromoelectric field, and to discriminate between the models of fundamental or adjoint tubes.

preprint2010arXiv

Colour fields in gauge invariant quenched SU(3) Lattice QCD

he colour fields, created by static sources belonging to different SU(3) representations, from the 3 to the 27, are computed in quenched SU(3) lattice QCD, in a 24^3 x 48 lattice at beta=6.2 and a=0.07261(85) fm. We utilize the technique of generalized Wilson Loops to localize the sources, correlated with plaquettes to measure the respective colour fields. We investigate the Casimir scaling of the fields, measured in the static potentials by Bali. We also study the coherence length, comparing with the dual Ginzburg-Landau approach. With the penetration and coherence lengths we determined the Ginzburg-Landau dimensionless parameter, this result is consistent with a type II superconductor picture, and with an effective dual gluon mass of 0.905 +- 0.163 GeV.

preprint2010arXiv

Colour fields of the static hybrid gluon-quark-antiquark system

The colour fields, created by a static gluon-quark-antiquark system, are computed in quenched SU(3) lattice QCD, in a $24^3\times 48$ lattice at $β=6.2$ and $a=0.07261(85)\,fm$. We study two geometries, one with a U shape and another with an L shape. The particular cases of the two gluon glueball and quark-antiquark are also studied, and the Casimir scaling is investigated in a microscopic perspective. This also contributes to understand confinement with flux tubes and to discriminate between the models of fundamental versus adjoint confining strings, analogous to type-II and type-I superconductivity.

preprint2010arXiv

Confinement and chiral symmetry breaking, one or two critical points?

We study the QCD phase diagram, in particular we study the critical points of the two main QCD phase transitions, confinement and chiral symmetry breaking. Confinement drives chiral symmetry breaking, and, due to the finite quark mass, at small density both transitions are a crossover, while they are a first or second order phase transition in large density. We study the QCD phase diagram with a quark potential model including both confinement and chiral symmetry. This formalism, in the Coulomb gauge hamiltonian formalism of QCD, is presently the only one able to microscopically include both a quark-antiquark confining potential and a vacuum condensate of quark-antiquark pairs. Our order parameters are the Polyakov loop and the quark mass gap. The confining potential is extracted from the Lattice QCD data of the Bielefeld group. We scan the QCD phase diagram for different quark masses, in order to address how the quark masses affect the critical point location in the phase diagram.

preprint2010arXiv

Confinement and chiral symmetry crossovers

We study the QCD phase diagram, in particular we study the critical points of the two main QCD phase transitions, confinement and chiral symmetry breaking. Confinement drives chiral symmetry breaking, and, due to the finite quark mass, at small density both transitions are a crossover, while they are a first or second order phase transition in large density. We study the QCD phase diagram with a quark potential model including both confinement and chiral symmetry. This formalism, in the Coulomb gauge hamiltonian formalism of QCD, is presently the only one able to microscopically include both a quark-antiquark confining potential and a vacuum condensate of quark-antiquark pairs. This model is able to address all the excited hadrons, and chiral symmetry breaking, at the same token. Our order parameters are the Polyakov loop and the quark mass gap. The confining potential is extracted from the Lattice QCD data of the Bielefeld group. We address how the quark masses affect the critical point location in the phase diagram.

preprint2010arXiv

Coulomb Gauge QCD and the Excited Hadron Spectrum

We discuss progress in understanding the light and heavy quark excited hadron spectrum from Coulomb gauge QCD. For light quark systems we highlight the insensitivity to spontaneous chiral symmetry breaking, which predicts Wigner parity-degeneracy in the highly excited hadron spectrum and allows the quark mass momentum dependence to be experimentally probed. For heavy quark meson decays we invoke the Franck-Condon principle, a consequence of small velocity changes for heavy quarks, to extract qualitative, but model independent, structure insight from the momentum distribution of the decay products.

preprint2010arXiv

Lattice SU(2) on GPU's

We discuss the CUDA approach to the simulation of pure gauge Lattice SU(2). CUDA is a hardware and software architecture developed by NVIDIA for computing on the GPU. We present an analysis and performance comparison between the GPU and CPU with single precision. Analysis with single and multiple GPU's, using CUDA and OPENMP, are also presented. In order to obtain a high performance, the code must be optimized for the GPU architecture, i.e., an implementation that exploits the memory hierarchy of the CUDA programming model. Using GPU texture memory and minimizing the data transfers between CPU and GPU, we achieve a speedup of $200\times$ using 2 NVIDIA 295 GTX cards relative to a serial CPU, which demonstrates that GPU's can serve as an efficient platform for scientific computing. With multi-GPU's we are able, in one day computation, to generate 1 000 000 gauge configurations in a $48^4$ lattice with $β=6.0$ and calculate the mean average plaquette. We present results for the mean average plaquette in several lattice sizes for different $β$. Finally we present results for the mean average Polyakov loop at finite temperature.

preprint2010arXiv

On confinement and chiral symmetry critical points

We study the QCD phase diagram, in particular we study the critical points of the two main QCD phase transitions, confinement and chiral symmetry breaking. Confinement drives chiral symmetry breaking, and, due to the finite quark mass, at small density both transitions are a crossover, while they are a first or second order phase transition in large density. We study the QCD phase diagram with a quark potential model including both confinement and chiral symmetry. This formalism, in the Coulomb gauge hamiltonian formalism of QCD, is presently the only one able to microscopically include both a quark-antiquark confining potential and a vacuum condensate of quark-antiquark pairs. This model is able to address all the excited hadrons, and chiral symmetry breaking, at the same token. Our order parameters are the Polyakov loop and the quark mass gap. The confining potential is extracted from the Lattice QCD data of the Bielefeld group. We address how the quark masses affect the critical point location in the phase diagram.

preprint2010arXiv

Tetraquark resonances with the triple flip-flop potential, decays in the cherry in a broken glass approximation

We develop a unitarized formalism to study tetraquarks using the triple flip-flop potential, which includes two meson-meson potentials and the tetraquark four-body potential. This can be related to the Jaffe-Wilczek and to the Karliner-Lipkin tetraquark models, where we also consider the possible open channels, since the four quarks and antiquarks may at any time escape to a pair of mesons. Here we study a simplified two-variable toy model and explore the analogy with a cherry in a glass, but a broken one where the cherry may escape from. It is quite interesting to have our system confined or compact in one variable and infinite in the other variable. In this framework we solve the two-variable Schrödinger equation in configuration space. With the finite difference method, we compute the spectrum, we search for localized states and we attempt to compute phase shifts. We then apply the outgoing spherical wave method to compute in detail the phase shifts and and to determine the decay widths. We explore the model in the equal mass case, and we find narrow resonances. In particular the existence of two commuting angular momenta is responsible for our small decay widths.

preprint2006arXiv

The BES f_0(1810): a new glueball candidate

We analyze the f_0(1810) state recently observed by the BES collaboration via radiative J/ψdecay to a resonant ϕωspectrum and confront it with DM2 data and glueball theory. The DM2 group only measured ωωdecays and reported a pseudoscalar but no scalar resonance in this mass region. A rescattering mechanism from the open flavored KKbar decay channel is considered to explain why the resonance is only seen in the flavor asymmetric ωϕbranch along with a discussion of positive C parity charmonia decays to strengthen the case for preferred open flavor glueball decays. We also calculate the total glueball decay width to be roughly 100 MeV, in agreement with the narrow, newly found f_0, and smaller than the expected estimate of 200-400 MeV. We conclude that this discovered scalar hadron is a solid glueball candidate and deserves further experimental investigation, especially in the K-Kbar channel. Finally we comment on other, but less likely, possible assignments for this state.

preprint2006arXiv

Vorticity and magnetic shielding in a type-II superconductor

We study in detail, solving the Bogoliubov-de Gennes equations, the magnetic field, supercurrent and order parameter profiles originated by a solenoid or magnetic whisker inserted in a type-II superconductor. We consider solutions of different vorticities, n, in the various cases. The results confirm the connection between the vorticity, the internal currents and the boundstates in a self-consistent way. The number of boundstates is given by the vorticity of the phase of the gap function as in the case with no external solenoid. In the limiting case of an infinitely thin solenoid, like a Dirac string, the solution is qualitatively different. The quasiparticle spectrum and wave functions are a function of n-n_ext, where n_ext is the vorticity of the solenoid. The flux is in all cases determined by the vorticity of the gap function.

preprint1999arXiv

A new Perspective on the Scalar meson Puzzle, from Spontaneous Chiral Symmetry Breaking Beyond BCS

We introduce coupled channels of Bethe-Salpeter mesons both in the boundstate equation for mesons and in the mass gap equation for chiral symmetry. Consistency is insured by the Ward Identities for axial currents, which preserve the Goldstone boson nature of the pion and prevents a systematic shift of the hadron spectrum. We study the decay of a scalar meson coupled to a pair of pseudoscalars. We also show that coupled channels reduce the breaking of chiral symmetry, with the same Feynman diagrams that appear in the coupling of a scalar meson to a pair of pseudoscalar mesons. Exact calculations are performed in a particular confining quark model, where we find that the groundstate $I=0, ^3P_0 q \bar q$ meson is the f_0(980) with a partial decay width of 40MeV. We also find a 30% reduction of the chiral condensate due to coupled channels.