Source author record

Frank X. Lee

Frank X. Lee 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

23works
6topics
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

23 published item(s)

preprint2022arXiv

Pole position of the $a_1(1260)$ resonance in a three-body unitary framework

Masses, widths, and branching ratios of hadronic resonances are quantified by their pole positions and residues with respect to transition amplitudes on the Riemann sheets of the complex energy-plane. In this study we discuss the analytic structure in the physical energy region of three-body scattering amplitudes on such manifolds. As an application, we determine the pole position of the $a_1(1260)$ meson from the ALEPH experiment by allowing for $πρ$ coupled channels in S- and D-wave. We find it to be $\sqrt{s_0}=(1232^{+15+9}_{-0-11}-i266^{+0+15}_{-22-27})~\text{MeV}$.

preprint2021arXiv

Towards charged hadron polarizabilities from four-point functions in lattice QCD

We show how to compute electromagnetic polarizabilities of charged hadrons using four-point functions in lattice QCD. The low-energy behavior of Compton scattering amplitude is matched to matrix elements of current-current correlation functions on the lattice. Working in momentum space, formulas for electric polarizability ($α_E$) and magnetic polarizability ($β_M$) are derived for both charged pion and proton. Lattice four-point correlation functions are constructed from quark and gluon fields to be used in Monte-Carlo simulations. The content of the functions is assessed in detail and specific prescriptions are given to isolate the polarizabilities. The connected quark-line diagrams can be done today as a small lattice project. The disconnected diagrams are more challenging but are within reach of dedicated resources for medium to large lattice projects. We also draw attention to the potential of four-point functions as a multi-purpose tool for hadron structure.

preprint2020arXiv

Finite-volume energy spectrum of the $K^-K^-K^-$ system

The dynamics of multi-kaon systems are of relevance for several areas of nuclear physics. However, even the simplest systems, two and three kaons, are hard to prepare and study experimentally. Here we show how to extract this information using first-principle lattice QCD results. We (1) extend the relativistic three-body quantization condition to the strangeness sector, predicting for the first time the excited level finite-volume spectrum of three kaon systems at maximal isospin, and (2) present a first lattice QCD calculation of the excited levels of this system in a finite box. We compare our predictions with the lattice results reported here and with previous ground state calculations and find very good agreement.

preprint2020arXiv

Roper State from Overlap Fermions

The Roper state is extracted with valence overlap fermions on a $2+1$-flavor domain-wall fermion lattice (spacing $a = 0.114$ fm and $m_π = 330$ MeV) using both the Sequential Empirical Bayes (SEB) method and the variational method. The results are consistent, provided that a large smearing-size interpolation operator is included in the variational calculation to have better overlap with the lowest radial excitation. Similar calculations carried out for an anisotropic clover lattice with similar parameters find the Roper $\approx 280$ MeV higher than that of the overlap fermion. The fact that the prediction of the Roper state by overlap fermions is consistently lower than those of clover fermions, chirally improved fermions, and twisted-mass fermions over a wide range of pion masses has been dubbed a "Roper puzzle." To understand the origin of this difference, we study the hairpin $Z$-diagram in the isovector scalar meson ($a_0$) correlator in the quenched approximation. Comparing the $a_0$ correlators for clover and overlap fermions, at a pion mass of 290 MeV, we find that the spectral weight of the ghost state with clover fermions is smaller than that of the overlap at $a = 0.12$ fm and $0.09$ fm, whereas the whole $a_0$ correlators of clover and overlap at $a = 0.06$ fm coincide within errors. This suggests that chiral symmetry is restored for clover at $a \le 0.06$ fm and that the Roper should come down at and below this $a$. We conclude that this work supports a resolution of the "Roper puzzle" due to $Z$-graph type chiral dynamics. This entails coupling to higher components in the Fock space (e.g. $Nπ$, $Nππ$ states) to induce the effective flavor-spin interaction between quarks as prescribed in the chiral quark model, resulting in the parity-reversal pattern as observed in the experimental excited states of $N, Δ$ and $Λ$.

preprint2019arXiv

A cross-channel study of pion scattering from lattice QCD

We use a chiral model for pion interactions, in the inverse amplitude formalism, to perform a simultaneous analysis of lattice QCD results for pion-pion scattering in all three isospin channels. The input is the finite-volume two-pion spectrum computed using lattice QCD from six ensembles on lattices elongated in one of the spatial dimensions. A two-flavor dynamical lattice QCD action is used with two quark masses corresponding to a pion mass of 315 MeV and 224 MeV. The spectrum in the elastic region is subjected to a global fit which takes into account full correlations across isospin, pion mass and decay constant. The parameters from the fit are used to perform a chiral extrapolation to the physical point. The cross-channel fit results in a more precise determination of the parameters of the model when compared with single channel fits. We obtain $m_πa_0^{I=0}=0.2132(9)$, and $m_πa_0^{I=2}=0.0433(2)$ as well as $m_σ=443(3)-i221(6)$ MeV and $m_ρ=724(4)-i67(1)$ MeV. Several aspects of scale setting and consistency with previous analyses of lattice QCD results are discussed as well.

preprint2015arXiv

Charm and strange quark masses and $f_{D_s}$ from overlap fermions

We use overlap fermions as valence quarks to calculate meson masses in a wide quark mass range on the $2+1$-flavor domain-wall fermion gauge configurations generated by the RBC and UKQCD Collaborations. The well-defined quark masses in the overlap fermion formalism and the clear valence quark mass dependence of meson masses observed from the calculation facilitate a direct derivation of physical current quark masses through a global fit to the lattice data, which incorporates $O(a^2)$ and $O(m_c^4a^4)$ corrections, chiral extrapolation, and quark mass interpolation. Using the physical masses of $D_s$, $D_s^*$ and $J/ψ$ as inputs, Sommer's scale parameter $r_0$ and the masses of charm quark and strange quark in the $\overline{\rm MS}$ scheme are determined to be $r_0=0.465(4)(9)$ fm, $m_c^{\overline{\rm MS}}(2\,{\rm GeV})=1.118(6)(24)$ GeV (or $m_c^{\overline{\rm MS}}(m_c)=1.304(5)(20)$ GeV), and $m_s^{\overline{\rm MS}}(2\,{\rm GeV})=0.101(3)(6)\,{\rm GeV}$, respectively. Furthermore, we observe that the mass difference of the vector meson and the pseudoscalar meson with the same valence quark content is proportional to the reciprocal of the square root of the valence quark masses. The hyperfine splitting of charmonium, $M_{J/ψ}-M_{η_c}$, is determined to be 119(2)(7) MeV, which is in good agreement with the experimental value. We also predict the decay constant of $D_s$ to be $f_{D_s}=254(2)(4)$ MeV. The masses of charmonium $P$-wave states $χ_{c0}, χ_{c1}$ and $h_c$ are also in good agreement with experiments.

preprint2014arXiv

Sea contributions to the electric polarizability of the hadrons

We present a lattice QCD calculation of the polarizability of the neutron and other neutral hadrons that includes the effects of the background field on the sea quarks. This is done by perturbatively reweighting the charges of the sea quarks to couple them to the background field. The main challenge in such a calculation is stochastic estimation of the weight factors, and we discuss the difficulties in this estimation. Here we use an extremely aggressive dilution scheme to reduce the stochastic noise to a manageable level. The pion mass in our calculation is 306 MeV and the lattice size is 3 fm. For neutron, we find that $α_E = 2.70(55) * 10^{-4}fm^3$, which is the most precise lattice QCD determination of the polarizability to date that includes sea effects.

preprint2014arXiv

Update on the sea contributions to hadron polarizabilities via reweighting

We have made significant progress on extending lattice QCD calculation of the polarizability of the neutron and other hadrons to include the effects of charged dynamical quarks. This is done by perturbatively reweighting the charges of the sea quarks to couple them to the background field. The dominant challenge in such a calculation is stochastic estimation of the weight factors, and we discuss the difficulties in this estimation. Here we use an extremely aggressive dilution scheme with N = 124,416 sources per configuration to reduce the stochastic noise to a manageable level. We find that α_E = 2.70(55) * 10^-4 fm^3 for the neutron on one ensemble. We show that low-mode substitution can be used in tandem with dilution to construct an even better estimator, and introduce the offdiagonal matrix element mapping technique for predicting estimator quality.

preprint2013arXiv

Update on the Sea Contributions to Hadron Electric Polarizabilities through Reweighting

We present the results of a reweighting calculation to compute the contribution of the charged quark sea to the neutron electric polarizability. The chief difficulty is the stochastic estimation of weight factors, and we present a hopping parameter expansion-based technique for reducing the stochastic noise, along with a discussion of why this particular reweighting is so difficult. We used this technique to estimate weight factors for 300 configurations of nHYP-clover fermions and compute the neutron polarizability, but the reweighting greatly inflates the overall statistical error, driven by the stochastic noise in the weight factors.

preprint2012arXiv

A Problem in Particle Physics and Its Bayesian Analysis

There is a class of statistical problems that arises in several contexts, the Lattice QCD problem of particle physics being one that has attracted the most attention. In essence, the problem boils down to the estimation of an infinite number of parameters from a finite number of equations, each equation being an infinite sum of exponential functions. By introducing a latent parameter into the QCD system, we are able to identify a pattern which tantamounts to reducing the system to a telescopic series. A statistical model is then endowed on the series, and inference about the unknown parameters done via a Bayesian approach. A computationally intensive Markov Chain Monte Carlo (MCMC) algorithm is invoked to implement the approach. The algorithm shares some parallels with that used in the particle Kalman filter. The approach is validated against simulated as well as data generated by a physics code pertaining to the quark masses of protons. The value of our approach is that we are now able to answer questions that could not be readily answered using some standard approaches in particle physics. The structure of the Lattice QCD equations is not unique to physics. Such architectures also appear in mathematical biology, nuclear magnetic imaging, network analysis, ultracentrifuge, and a host of other relaxation and time decay phenomena. Thus, the methodology of this paper should have an appeal that transcends the Lattice QCD scenario which motivated us.

preprint2011arXiv

$ρ$ meson decay on asymmetrical lattices

We present a lattice QCD calculation of the characteristics of the $ρ$ meson decay. The study is carried out on spatially asymmetric boxes using nHYP-smeared clover fermions in the quenched approximation. The resonance mass and coupling constant are calculate using the P-wave scattering phaseshifts, of the isospin I=1 two-pion system. We use pion masses m_π= 418 MeV and m_π=312 MeV. In both cases, the $ρ$ decay is kinematically feasible. We work on lattice sizes N_z X 24^2 X 48 with lattice spacing a=0.1 fm and N_z=24,30,34,48.

preprint2011arXiv

Efficient implementation of the overlap operator on multi-GPUs

Lattice QCD calculations were one of the first applications to show the potential of GPUs in the area of high performance computing. Our interest is to find ways to effectively use GPUs for lattice calculations using the overlap operator. The large memory footprint of these codes requires the use of multiple GPUs in parallel. In this paper we show the methods we used to implement this operator efficiently. We run our codes both on a GPU cluster and a CPU cluster with similar interconnects. We find that to match performance the CPU cluster requires 20-30 times more CPU cores than GPUs.

preprint2011arXiv

Spin Polarizabilities on the Lattice

Spin polarizabilities provide information on the internal structure of hadrons in the presence of weak external electromagnetic fields, and are actively studied by Compton scattering experiments. They provide finer detail than the regular polarizabilities since they require space and time-varying fields. Using an effective action in the weak field limit, we have identified methods to isolate each of the physical quantities ($μ, α, β, γ_{E1}, γ_{M1}, γ_{E2}, γ_{M2}$) for spin-1/2 hadrons, both neutral and charged. We also perform a lattice QCD simulation to investigate the feasibility of the effective action approach.

preprint2011arXiv

Transition magnetic moment of $Λ$ $\to$ $Σ^0$ in QCD sum rules

The $Λ$ $\to$ $Σ^0$ transition magnetic moment is computed in the QCD sum rules approach. Three independent tensor structures are derived in the external field method using generalized interpolating fields. They are analyzed together with the $Λ$ and $Σ^0$ mass sum rules using a Monte-Carlo-based analysis, with attention to OPE convergence, ground-state dominance, and the role of the transitions in the intermediate states. Relations between sum rules for magnetic moments of $Λ$ and $Σ^0$ and sum rules for transition magnetic moment of $Λ$ $\to$ $Σ^0$ are also examined. Our best prediction for the transition magnetic moment is $μ_{Σ^0Λ}= 1.60\pm 0.07\; μ_N$. A comparison is made with other calculations in the literature.

preprint2010arXiv

Eta-nucleon coupling constant in QCD with SU(3) symmetry breaking

We study the $η$NN coupling constant using the method of QCD sum rules starting from the vacuum-to-eta correlation function of the interpolating fields of two nucleons. The matrix element of this correlation has been taken with respect to nucleon spinors to avoid unwanted pole contribution. The SU(3)-flavor symmetry breaking effects have been accounted for via the $η$-mass, s-quark mass and eta decay constant to leading order. Out of the four sum rules obtained by taking the ratios of the two sum rules in conjunction with the two sum rules in nucleon mass, three are found to give mutually consistent results. We find the SU(3) breaking effects significant, as large as 50% of the SU(3) symmetric part.

preprint2010arXiv

Magnetic Moments of Negative-Parity Baryons from Lattice QCD

We report preliminary results for the magnetic moments of negative-parity baryons extracted from mass shifts in the presence of static external magnetic fields. The calculations are done on $24^3\times 48$ quenched lattices using standard Wilson actions, with $β$=6.0 and pion mass down to about 520 MeV, and 1000 configurations. Reasonable signals for the negative-parity states are observed and the sign of their magnetic moments is established. The results are compared to model calculations.

preprint2009arXiv

Neutron electric polarizability

We use the background field method to extract the "connected" piece of the neutron electric polarizability. We present results for quenched simulations using both clover and Wilson fermions and discuss our experience in extracting the mass shifts and the challenges we encountered when we lowered the quark mass. For the neutron we find that as the pion mass is lowered below $500\MeV$, the polarizability starts rising in agreement with predictions from chiral perturbation theory. For our lowest pion mass, $m_π=320\MeV$, we find that $α_n = 3.8(1.3)\times 10^{-4}\fm^3$, which is still only one third of the experimental value. We also present results for the neutral pion; we find that its polarizability turns negative for pion masses smaller than $500\MeV$ which is puzzling.

preprint2008arXiv

The background field method on the lattice

The background field method has been used successfully to determine hadron electromagnetic polarizabilities. Recently questions have been raised regarding the proper way to deal with the electric field on the lattice. In this paper, we show how the presence of a background electric field affects the quark hopping matrix. We use this formulation to carry out simulations on quenched configurations, and we present the results for neutron electric polarizability for pion masses as low as $500\MeV$. We find that the polarizability is roughly constant for the quark masses considered, $α\sim 1.5 \times 10^{-4} \fm^3$. While the polarizability is positive, it is significantly smaller than the experimental value due to the fact that the quark masses used are too large.

preprint2004arXiv

Topological Charge Evolution in the Markov-Chain of QCD

The topological charge is studied on lattices of large physical volume and fine lattice spacing. We illustrate how a parity transformation on the SU(3) link-variables of lattice gauge configurations reverses the sign of the topological charge and leaves the action invariant. Random applications of the parity transformation are proposed to traverse from one topological charge sign to the other. The transformation provides an improved unbiased estimator of the ensemble average and is essential in improving the ergodicity of the Markov chain process.

preprint2002arXiv

Predictive Ability of QCD Sum Rules for Excited Baryons

The masses of octet baryons are calculated by the method of QCD sum rules. Using generalized interpolating fields, three independent sets of QCD sum rules are derived which allow the extraction of low-lying N* states with spin-parity 1/2+, 1/2- and 3/2- in both the non-strange and strange channels. The predictive ability of the sum rules is examined by a Monte-Carlo based analysis procedure in which the three phenomenological parameters (mass, coupling, threshold) are treated as free parameters simultaneously. Realistic uncertainties in these parameters are obtained by simultaneously exploring all uncertainties in the QCD input parameters. Those sum rules with good predictive power are identified and their predictions are compared with experiment where available.

preprint2000arXiv

N* Mass Spectrum from an Anisotropic Action

Results are reported for N* masses in the 1/2+ and 1/2- sectors on an anisotropic lattice. The gauge action is the usual plaquette plus rectangle type, the quark action is of the D234 type, both having tadpole-improved tree-level coefficients. Clear splittings from the nucleon ground state are observed with smeared operators and 500 configurations. The first excited states in each channel are further isolated, using a correlation matrix and the variational method. The basic pattern of these splittings is consistent with experimental observations.

preprint1997arXiv

Light Hadron Spectroscopy on Coarse Lattices with O(a^2) Mean-Field Improved Actions

The masses and dispersions of light hadrons are calculated in lattice QCD using an O(a^2) tadpole-improved gluon action and an O(a^2) tadpole-improved next-nearest-neighbor fermion action originally proposed by Hamber and Wu. Two lattices of constant volume with lattice spacings of approximately 0.40 fm and 0.24 fm are considered. The results reveal some scaling violations at the coarser lattice spacing on the order of 5%. At the finer lattice spacing, the calculated mass ratios reproduce state-of-the-art results using unimproved actions. Good dispersion and rotational invariance up to momenta of pa ~ 1 are also found. The relative merit of alternative choices for improvement operators is assessed through close comparisons with other plaquette-based tadpole-improved actions.