Full QCD calculation of neutron electric dipole moment with the external electric field method
Shintani, E.; Aoki, S.; Kuramashi, Y.
2008-07-01
We have calculated the neutron electric dipole moment (EDM) in the presence of the CP violating {theta} term in lattice QCD with two-flavor dynamical clover quarks, using the external electric field method. Accumulating a large number of statistics by the averages over 16 different source points and over forward and backward nucleon propagators, we have obtained nonzero signals of neutron and proton EDM beyond 1 standard deviation at each quark mass in full QCD. We have investigated the quark mass dependence of nucleon EDM in full QCD, and have found that nucleon EDM in full QCD does not decrease toward the chiral limit, as opposed to the theoretical expectation. We briefly discuss possible reasons for this behavior.
Nonperturbative QCD Calculations
NASA Astrophysics Data System (ADS)
Dellby, Niklas
1995-01-01
The research described in this thesis is an exact transformation of the Yang-Mills quantum chromodynamics (QCD) Lagrangrian into a form that is suitable for nonperturbative calculations. The conventional Yang-Mills Lagrangian has proven to be an excellent basis for perturbative calculations, but in nonperturbative calculations it is difficult to separate gauge problems from physical properties. To mitigate this problem, I develop a new equivalent Lagrangian that is not only expressed completely in terms of the field strengths ofthe gauge field but is also manifestly Lorentz and gauge invariant. The new Lagrangian is quadratic in derivatives, with non-linear local couplings, thus it is ideally suited for a numerical calculation. The field-strength Lagrangian is of such a form that it is possible to do a straightforward numerical stationary path expansion and find the fundamental QCD properties. This thesis examines several approximations analytically, investigating different ways to utilize the new Lagrangian. (Copies available exclusively from MIT Libraries, Rm. 14-0551, Cambridge, MA 02139-4307. Ph. 617-253-5668; Fax 617-253 -1690.).
Nucleon Axial Charge in Full Lattice QCD
Edwards, R.G.; Richards, D.G.; Fleming, G.T.; Haegler, Ph.; Negele, J.W.; Pochinsky, A.V.; Orginos, K.; Renner, D.B.; Schroers, W.
2006-02-10
The nucleon axial charge is calculated as a function of the pion mass in full QCD. Using domain wall valence quarks and improved staggered sea quarks, we present the first calculation with pion masses as light as 354 MeV and volumes as large as (3.5 fm){sup 3}. We show that finite volume effects are small for our volumes and that a constrained fit based on finite volume chiral perturbation theory agrees with experiment within 7% statistical errors.
Strange Baryon Physics in Full Lattice QCD
Huey-Wen Lin
2007-11-01
Strange baryon spectra and form factors are key probes to study excited nuclear matter. The use of lattice QCD allows us to test the strength of the Standard Model by calculating strange baryon quantities from first principles.
QCD Factorization and PDFs from Lattice QCD Calculation
NASA Astrophysics Data System (ADS)
Ma, Yan-Qing; Qiu, Jian-Wei
2015-02-01
In this talk, we review a QCD factorization based approach to extract parton distribution and correlation functions from lattice QCD calculation of single hadron matrix elements of quark-gluon operators. We argue that although the lattice QCD calculations are done in the Euclidean space, the nonperturbative collinear behavior of the matrix elements are the same as that in the Minkowski space, and could be systematically factorized into parton distribution functions with infrared safe matching coefficients. The matching coefficients can be calculated perturbatively by applying the factorization formalism on to asymptotic partonic states.
The {Lambda}(1405) in Full QCD
Menadue, Benjamin J.; Kamleh, Waseem; Leinweber, Derek B.; Mahbub, M. Selim
2011-12-14
At 1405.1 MeV, the lowest-lying negative-parity state of the {Lambda} baryon lies surprising low. Indeed, this is lower than the lowest negative-parity state of the nucleon, even though the {Lambda}(1405) possesses a valence strange quark. However, previous Lattice QCD studies have been unable to identify such a low-lying state. Using the PACS-CS (2+1)-flavour full-QCD ensembles, available through the ILDG, we utilise a variational analysis with source and sink smearing to isolate this elusive state. We find three low-lying odd-parity states, and for the first time reproduce the correct level ordering with respect to the nearby scattering thresholds.
The nucleon axial charge in full lattice QCD
Robert Edwards; George Fleming; P Hagler; John Negele; Konstantinos Orginos; Andrew Pochinsky; Dru Renner; David Richards; Wolfram Schroers
2005-10-13
The nucleon axial charge is calculated as a function of the pion mass in full QCD. Using domain wall valence quarks and improved staggered sea quarks, we present the first calculation with pion masses as light as 354 MeV and volumes as large as (3.5 fm)3. We show that finite volume effects are small for our volumes and that a constrained fit based on finite volume chiral perturbation theory agrees with experiment within 5% statistical errors.
Full QCD in external chromomagnetic field
NASA Astrophysics Data System (ADS)
Cosmai, Leonardo
2006-12-01
We investigate the deconfining phase transition in full QCD with two flavors of staggered fermions in presence of a constant abelian chromomagnetic field. We find that the deconfinement tempera- ture decreases and eventually goes to zero by increasing the strength of the chromomagnetic field. Moreover our results suggest that the chiral transition coincides with the deconfinement transition and therefore even the chiral critical temperature depends on the applied chromomagnetic field. We also find that the chiral condensate increases with the strength of the chromomagnetic field.
Calculation of the nucleon axial charge in lattice QCD
D. B. Renner; R. G. Edwards; G. Fleming; Ph. Hagler; J. W. Negele; K. Orginos; A. V. Pochinsky; D. G. Richards; W. Schroers
2006-09-01
Protons and neutrons have a rich structure in terms of their constituents, the quarks and gluons. Understanding this structure requires solving Quantum Chromodynamics (QCD). However QCD is extremely complicated, so we must numerically solve the equations of QCD using a method known as lattice QCD. Here we describe a typical lattice QCD calculation by examining our recent computation of the nucleon axial charge.
Consistent Perturbative Fixed Point Calculations in QCD and Supersymmetric QCD
NASA Astrophysics Data System (ADS)
Ryttov, Thomas A.
2016-08-01
We suggest how to consistently calculate the anomalous dimension γ* of the ψ ¯ ψ operator in finite order perturbation theory at an infrared fixed point for asymptotically free theories. If the n +1 loop beta function and n loop anomalous dimension are known, then γ* can be calculated exactly and fully scheme independently in a Banks-Zaks expansion through O (Δfn) , where Δf=N¯ f-Nf , Nf is the number of flavors, and N¯f is the number of flavors above which asymptotic freedom is lost. For a supersymmetric theory, the calculation preserves supersymmetry order by order in Δf. We then compute γ* through O (Δf2) for supersymmetric QCD in the dimensional reduction scheme and find that it matches the exact known result. We find that γ* is astonishingly well described in perturbation theory already at the few loops level throughout the entire conformal window. We finally compute γ* through O (Δf3) for QCD and a variety of other nonsupersymmetric fermionic gauge theories. Small values of γ* are observed for a large range of flavors.
Consistent Perturbative Fixed Point Calculations in QCD and Supersymmetric QCD.
Ryttov, Thomas A
2016-08-12
We suggest how to consistently calculate the anomalous dimension γ_{*} of the ψ[over ¯]ψ operator in finite order perturbation theory at an infrared fixed point for asymptotically free theories. If the n+1 loop beta function and n loop anomalous dimension are known, then γ_{*} can be calculated exactly and fully scheme independently in a Banks-Zaks expansion through O(Δ_{f}^{n}), where Δ_{f}=N[over ¯]_{f}-N_{f}, N_{f} is the number of flavors, and N[over ¯]_{f} is the number of flavors above which asymptotic freedom is lost. For a supersymmetric theory, the calculation preserves supersymmetry order by order in Δ_{f}. We then compute γ_{*} through O(Δ_{f}^{2}) for supersymmetric QCD in the dimensional reduction scheme and find that it matches the exact known result. We find that γ_{*} is astonishingly well described in perturbation theory already at the few loops level throughout the entire conformal window. We finally compute γ_{*} through O(Δ_{f}^{3}) for QCD and a variety of other nonsupersymmetric fermionic gauge theories. Small values of γ_{*} are observed for a large range of flavors. PMID:27563948
Full CKM matrix with lattice QCD
Okamoto, Masataka; /Fermilab
2004-12-01
The authors show that it is now possible to fully determine the CKM matrix, for the first time, using lattice QCD. |V{sub cd}|, |V{sub cs}|, |V{sub ub}|, |V{sub cb}| and |V{sub us}| are, respectively, directly determined with the lattice results for form factors of semileptonic D {yields} {pi}lv, D {yields} Klv, B {yields} {pi}lv, B {yields} Dlv and K {yields} {pi}lv decays. The error from the quenched approximation is removed by using the MILC unquenced lattice gauge configurations, where the effect of u, d and s quarks is included. The error from the ''chiral'' extrapolation (m{sub l} {yields} m{sub ud}) is greatly reduced by using improved staggered quarks. The accuracy is comparable to that of the Particle Data Group averages. In addition, |V{sub ud}|, |V{sub ts}|, |V{sub ts}| and |V{sub td}| are determined by using unitarity of the CKM matrix and the experimental result for sin (2{beta}). In this way, they obtain all 9 CKM matrix elements, where the only theoretical input is lattice QCD. They also obtain all the Wolfenstein parameters, for the first time, using lattice QCD.
Unlimited full configuration interaction calculations
NASA Astrophysics Data System (ADS)
Knowles, Peter J.; Handy, Nicholas C.
1989-08-01
In very large full configuration interaction (full CI), nearly all of the CI coefficients are very small. Calculations, using a newly developed algorithm which exploits this fact, on NH3 with a DZP basis are reported, involving 2×108 Slater determinants. Such calculations are impossible with other existing full CI codes. The new algorithm opens up the opportunity of full CI calculations which are unlimited in size.
A Framework for Lattice QCD Calculations on GPUs
Winter, Frank; Clark, M A; Edwards, Robert G; Joo, Balint
2014-08-01
Computing platforms equipped with accelerators like GPUs have proven to provide great computational power. However, exploiting such platforms for existing scientific applications is not a trivial task. Current GPU programming frameworks such as CUDA C/C++ require low-level programming from the developer in order to achieve high performance code. As a result porting of applications to GPUs is typically limited to time-dominant algorithms and routines, leaving the remainder not accelerated which can open a serious Amdahl's law issue. The lattice QCD application Chroma allows to explore a different porting strategy. The layered structure of the software architecture logically separates the data-parallel from the application layer. The QCD Data-Parallel software layer provides data types and expressions with stencil-like operations suitable for lattice field theory and Chroma implements algorithms in terms of this high-level interface. Thus by porting the low-level layer one can effectively move the whole application in one swing to a different platform. The QDP-JIT/PTX library, the reimplementation of the low-level layer, provides a framework for lattice QCD calculations for the CUDA architecture. The complete software interface is supported and thus applications can be run unaltered on GPU-based parallel computers. This reimplementation was possible due to the availability of a JIT compiler (part of the NVIDIA Linux kernel driver) which translates an assembly-like language (PTX) to GPU code. The expression template technique is used to build PTX code generators and a software cache manages the GPU memory. This reimplementation allows us to deploy an efficient implementation of the full gauge-generation program with dynamical fermions on large-scale GPU-based machines such as Titan and Blue Waters which accelerates the algorithm by more than an order of magnitude.
ColorFull: a C++ library for calculations in SU(Nc) color space
NASA Astrophysics Data System (ADS)
Sjödahl, Malin
2015-05-01
ColorFull, a C++ package for treating QCD color structure, is presented. ColorFull, which utilizes the trace basis approach, is intended for interfacing with event generators, but can also be used as a stand-alone package for squaring QCD amplitudes, calculating interferences, and describing the effect of gluon emission and gluon exchange.
Determination of the Lambda parameter from full lattice QCD
Goeckeler, M.; Horsley, R.; Irving, A.C.; Rakow, P.E.L.; Pleiter, D.; Schierholz, G.; Stueben, H.
2006-01-01
We present a determination of the QCD parameter {lambda} in the quenched approximation (n{sub f}=0) and for two flavors (n{sub f}=2) of light dynamical quarks. The calculations are performed on the lattice using O(a) improved Wilson fermions and include taking the continuum limit. We find {lambda}{sub n{sub f}=0}{sup MS}=259(1)(19) MeV and {lambda}{sub n{sub f}=2}{sup MS}=261(17)(26) MeV, using r{sub 0}=0.467 fm to set the scale. Extrapolating our results to five flavors, we obtain for the running coupling constant at the mass of the Z boson {alpha}{sub s}{sup MS}(m{sub Z})=0.112(1)(2)
Calculation of hadronic matrix elements using lattice QCD
Gupta, R.
1993-08-01
The author gives a brief introduction to the scope of lattice QCD calculations in his effort to extract the fundamental parameters of the standard model. This goal is illustrated by two examples. First the author discusses the extraction of CKM matrix elements from measurements of form factors for semileptonic decays of heavy-light pseudoscalar mesons such as D {yields} Ke{nu}. Second, he presents the status of results for the kaon B parameter relevant to CP violation. He concludes the talk with a short outline of his experiences with optimizing QCD codes on the CM5.
Uncertainty quantification in lattice QCD calculations for nuclear physics
Beane, Silas R.; Detmold, William; Orginos, Kostas; Savage, Martin J.
2015-02-05
The numerical technique of Lattice QCD holds the promise of connecting the nuclear forces, nuclei, the spectrum and structure of hadrons, and the properties of matter under extreme conditions with the underlying theory of the strong interactions, quantum chromodynamics. A distinguishing, and thus far unique, feature of this formulation is that all of the associated uncertainties, both statistical and systematic can, in principle, be systematically reduced to any desired precision with sufficient computational and human resources. As a result, we review the sources of uncertainty inherent in Lattice QCD calculations for nuclear physics, and discuss how each is quantified in current efforts.
Choice of nucleon current in QCD sum rule calculations
Indumathi, D.
1995-06-01
The author studies the correlator occurring in the calculation of nucleon structure functions in the QCD sum rule technique. The author shows+ that the property of gauge invariance of this four-point function in the spin dependent sector uniquely fixes the choice of nucleon current in such computations. No such constraint is obtainable from similar considerations in the spin-independent sector. 14 refs., 5 figs.
First-principles Calculation of Excited State Spectra in QCD
Jozef Dudek,Robert Edwards,Michael Peardon,David Richards,Christopher Thomas
2011-05-01
Recent progress at understanding the excited state spectra of mesons and baryons is described. I begin by outlining the application of the variational method to compute the spectrum of QCD, and then present results for the excited meson spectrum, with continuum quantum numbers of the states clearly delineated. I emphasise the need to extend the calculation to encompass multi-hadron contributions, and describe a recent calculation of the I=2 pion-pion energy-dependent phase shifts as a precursor to the study of channels with resonant behavior. I conclude with recent results for the low lying baryon spectrum, and the prospects for future calculations.
First-principles Calculation of Excited State Spectra in QCD
Dudek, Jozef J.; Edwards, Robert G.; Richards, David G.; Thomas, Christopher E.; Peardon, Michael J.
2011-05-24
Recent progress at understanding the excited state spectra of mesons and baryons is described. I begin by outlining the application of the variational method to compute the spectrum of QCD, and then present results for the excited meson spectrum, with continuum quantum numbers of the states clearly delineated. I emphasise the need to extend the calculation to encompass multi-hadron contributions, and describe a recent calculation of the I = 2{pi}{pi} energy-dependent phase shifts as a precursor to the study of channels with resonant behavior. I conclude with recent results for the low lying baryon spectrum, and the prospects for future calculations.
Nucleon Generalized Parton Distributions from Full Lattice QCD
Robert Edwards; Philipp Haegler; David Richards; John Negele; Konstantinos Orginos; Wolfram Schroers; Jonathan Bratt; Andrew Pochinsky; Michael Engelhardt; George Fleming; Bernhard Musch; Dru Renner
2007-07-03
We present a comprehensive study of the lowest moments of nucleon generalized parton distributions in N_f=2+1 lattice QCD using domain wall valence quarks and improved staggered sea quarks. Our investigation includes helicity dependent and independent generalized parton distributions for pion masses as low as 350 MeV and volumes as large as (3.5 fm)^3.
Nucleon generalized parton distributions from full lattice QCD
Haegler, Ph.; Musch, B.; Schroers, W.; Edwards, R. G.; Richards, D. G.; Engelhardt, M.; Fleming, G. T.; Orginos, K.; Renner, D. B.
2008-05-01
We present a comprehensive study of the lowest moments of nucleon generalized parton distributions in N{sub f}=2+1 lattice QCD using domain-wall valence quarks and improved staggered sea quarks. Our investigation includes helicity dependent and independent generalized parton distributions for pion masses as low as 350 MeV and volumes as large as (3.5 fm){sup 3}, for a lattice spacing of 0.124 fm. We use perturbative renormalization at one-loop level with an improvement based on the nonperturbative renormalization factor for the axial vector current, and only connected diagrams are included in the isosinglet channel.
Advancing QCD-based calculations of energy loss
NASA Astrophysics Data System (ADS)
Tywoniuk, Konrad
2013-08-01
We give a brief overview of the basics and current developments of QCD-based calculations of radiative processes in medium. We put an emphasis on the underlying physics concepts and discuss the theoretical uncertainties inherently associated with the fundamental parameters to be extracted from data. An important area of development is the study of the single-gluon emission in medium. Moreover, establishing the correct physical picture of multi-gluon emissions is imperative for comparison with data. We will report on progress made in both directions and discuss perspectives for the future.
Very large full configuration interaction calculations
NASA Astrophysics Data System (ADS)
Knowles, Peter J.
1989-03-01
The extreme sparsity of the solution of the full configuration interaction (full CI) secular equations is exploited in a new algorithm. For very large problems, the high speed memory, disk storage, and CPU requirements are reduced considerably, compared to previous techniques. This allows the possibility of full CI calculations with more than 10 8 Slater determinants. The power of the method is demonstrated in preliminary full CI calculations for the NH molecule, including up to 27901690 determinants.
Full CI benchmark calculations on CH3
NASA Technical Reports Server (NTRS)
Bauschlicher, Charles W., Jr.; Taylor, Peter R.
1987-01-01
Full CI calculations have been performed on the CH3 radical. The full CI results are compared to those obtained using CASSCF/multireference CI and coupled-pair functional methods, both at the equilibrium CH distance and at geometries with the three CH bonds extended. In general, the performance of the approximate methods is similar to that observed in calculations on other molecules in which one or two bonds were stretched.
Improved semileptonic form factor calculations in lattice QCD
Evans, Richard; Bali, Gunnar; Collins, Sara
2010-11-01
We investigate the computational efficiency of two stochastic based alternatives to the sequential propagator method used in lattice QCD calculations of heavy-light semileptonic form factors. In the first method, we replace the sequential propagator, which couples the calculation of two of the three propagators required for the calculation, with a stochastic propagator so that the calculations of all three propagators are independent. This method is more flexible than the sequential propagator method but introduces stochastic noise. We study the noise to determine when this method becomes competitive with the sequential propagator method, and find that for any practical calculation it is competitive with or superior to the sequential propagator method. We also examine a second stochastic method, the so-called 'one-end trick', concluding it is relatively inefficient in this context. The investigation is carried out on two gauge field ensembles, using the nonperturbatively improved Wilson-Sheikholeslami-Wohlert action with N{sub f}=2 mass-degenerate sea quarks. The two ensembles have similar lattice spacings but different sea-quark masses. We use the first stochastic method to extract O(a)-improved, matched lattice results for the semileptonic form factors on the ensemble with lighter sea quarks, extracting f{sub +}(0).
First Calculation of Hyperon Axial Couplings from Lattice QCD
Huey-Wen Lin; Konstantinos Orginos
2007-12-06
In this work, we report the first lattice calculation of the hyperon axial couplings, using the 2+1-flavor MILC configurations and domain-wall fermion valence quarks. Both the $\\Sigma$ and $\\Xi$ axial couplings are for the first time done in lattice QCD, and we find the numbers with greater precision than previous chiral perturbation theory and large-$N_c$ theory estimate: $g_{\\Sigma\\Sigma} = 0.450(21)_{\\rm stat}(22)_{\\rm syst}$ and $g_{\\Xi\\Xi} = -0.277(15)_{\\rm stat}(16)_{\\rm syst}$. As a side product, we also determine the low-energy chiral parameters $D$ and $F$ extracted from these coupling constants: $D=0.715(6)_{\\rm stat}(6)_{\\rm syst}$ and $F=0.453(5)_{\\rm stat}(5)_{\\rm syst}$.
Calculation of the equation of state of QCD at finite chemical and zero temperature
Zong Hongshi; Sun Weimin
2008-09-01
In this paper, we give a direct method for calculating the partition function, and hence the equation of state (EOS) of quantum chromodynamics (QCD) at finite chemical potential and zero temperature. In the EOS derived in this paper the pressure density is the sum of two terms: the first term P({mu})|{sub {mu}}{sub =0} (the pressure density at {mu}=0) is a {mu}-independent constant; the second term, which is totally determined by G{sub R}[{mu}](p) (the renormalized dressed quark propagator at finite {mu}), contains all the nontrivial {mu}-dependence. By applying a general result in the rainbow-ladder approximation of the Dyson-Schwinger approach obtained in our previous study [Phys. Rev. C 71, 015205 (2005)], G{sub R}[{mu}](p) is calculated from the meromorphic quark propagator proposed in [Phys. Rev. D 70, 014014 (2004)]. From this the full analytic expression of the EOS of QCD at finite {mu} and zero T is obtained (apart from the constant term P({mu})|{sub {mu}}{sub =0} which can in principle be calculated from the Cornwall-Jackiw-Tomboulis effective action). A comparison between our EOS and the cold, perturbative EOS of QCD of Fraga, Pisarski, and Schaffner-Bielich is made. It is expected that our EOS can provide a possible new approach for the study of neutron stars.
ϒ and ϒ' leptonic widths, aμb , and mb from full lattice QCD
NASA Astrophysics Data System (ADS)
Colquhoun, B.; Dowdall, R. J.; Davies, C. T. H.; Hornbostel, K.; Lepage, G. P.; Hpqcd Collaboration
2015-04-01
We determine the decay rate to leptons of the ground-state ϒ meson and its first radial excitation in lattice QCD for the first time. We use radiatively improved nonrelativistic QCD for the b quarks and include u , d , s and c quarks in the sea with u /d masses down to their physical values. We find Γ (ϒ →e+e-)=1.19 (11 ) keV and Γ (ϒ'→e+e-)=0.69 (9 ) keV , both in good agreement with experimental results. The decay constants we obtain are included in a summary plot of meson decay constants from lattice QCD given in the Conclusions. We also test time moments of the vector current-current correlator against values determined from the b -quark contribution to σ (e+e-→hadrons) and calculate the b -quark piece of the hadronic vacuum polarization contribution to the anomalous magnetic moment of the muon, aμb=0.271 (37 )×1 0-10 . Finally we determine the b -quark mass, obtaining in the MS ¯ scheme, m¯ b(m¯ b,nf=5 )=4.196 (23 ) GeV , the most accurate result from lattice QCD to date.
Bornyakov, V.G.
2005-06-01
Possibilities that are provided by a lattice regularization of QCD for studying nonperturbative properties of QCD are discussed. A review of some recent results obtained from computer calculations in lattice QCD is given. In particular, the results for the QCD vacuum structure, the hadron mass spectrum, and the strong coupling constant are considered.
NASA Astrophysics Data System (ADS)
Fleming, Sean
In this talk I review recent experimental and theoretical results in QCD. Since the topic is too vast to cover within given time constraints I choose to highlight some of the subjects that I find particularly exciting. On the experimental side I focus on measurements made at the Tevatron. Specifically jet production rates, and the cross section for B meson production. In addition I discuss an interesting measurement made by the Belle collaboration of double exclusive charmonium production. On the theory side I quickly review recent advances in computing hadronic cross sections at subleading order in perturbation theory. I then move on to soft-collinear effective theory. After a lightning review of the formalism I discuss recently published results on color-suppressed B → D decays.
The Excited-state Spectrum of QCD through Lattice Gauge Theory Calculations
David Richards
2012-12-01
I describe recent progress at understanding the excited state spectrum of QCD through lattice gauge calculations. I begin by outlining the evolution of the lattice effort at JLab. I detail the impact of recent lattice calculations on the present and upcoming experimental programs, and in particular that of the 12 GeV upgrade of Jefferson Laboratory. I conclude with the prospect for future calculations.
Calculation of equation of state of QCD at zero temperature and finite chemical potential
NASA Astrophysics Data System (ADS)
Jiang, Yu; Li, Ning; Sun, Wei-Min; Zong, Hong-Shi
2010-09-01
In this paper we calculate the equation of state (EOS) of QCD at zero temperature and finite chemical potential by using several models of quark propagators including the Dyson-Schwinger equations (DSEs) model, the hard-dense-loop (HDL) approximation and the quasi-particle model. The results are analyzed and compared with the known results in the literature.
Reduction of One-Loop Amplitudes at the Integrand Level -- NLO QCD Calculations
NASA Astrophysics Data System (ADS)
Ossola, G.; Papadopoulos, C. G.; Pittau, R.
2008-07-01
The recently proposed method (OPP) to extract the coefficients of the scalar one-loop integrals to any multi-particle (sub)-amplitude is described. Within this method no analytical information on the structure of the amplitude is needed, allowing for a purely numerical, but still algebraic, implementation of the algorithm. The algorithm can be used to automatically perform one-loop calculation both in QCD and in the EW Theory. As an application, we give QCD one-loop results for the process p p to ZZZ at the LHC.
A Review of Nucleon Spin Calculations in Lattice QCD
Huey-Wen Lin
2009-08-01
We review recent progress on lattice calculations of nucleon spin structure, including the parton distribution functions, form factors, generalization parton distributions, and recent developments in lattice techniques.
Calculation of neutral weak nucleon form factors with the AdS/QCD correspondence
NASA Astrophysics Data System (ADS)
Lohmann, Mark
The AdS/QCD (Anti-de Sitter/Quantum Chromodynamics) is a mathematical formalism applied to a theory based on the original AdS/CFT (Anti-de Sitter/ Conformal Field Theory) correspondence. The aim is to describe properties of the strong force in an essentially non-perturbative way. AdS/QCD theories break the conformal symmetry of the AdS metric (a sacrifice) to arrive at a boundary theory which is QCD-like (a payoff). This correspondence has been used to calculate well-known quantities in nucleon spectra and structure like Regge trajectories, form factors, and many others within an error of less than 20% from experiment. This is impressive considering that ordinary perturbation theory in QCD applied to the strongly interacting domain usually obtains an error of about 30%. In this thesis, the AdS/QCD correspondence method of light-front holography established by Brodsky and de Teramond is used in an attempt to calculate the Dirac and Pauli neutral weak form factors, FZ1 (Q2) and FZ2 (Q 2) respectively, for both the proton and the neutron. With this approach, we were able to determine the neutral weak Dirac form factor for both nucleons and the Pauli form factor for the proton, while the method did not succeed at determining the neutral weak Pauli form factor for the neutron. With these we were also able to extract the proton's strange electric and magnetic form factor, which addresses important questions in nucleon sub-structure that are currently being investigated through experiments at the Thomas Jefferson National Accelerator Facility.
Uncertainty quantification in lattice QCD calculations for nuclear physics
NASA Astrophysics Data System (ADS)
Beane, Silas R.; Detmold, William; Orginos, Kostas; Savage, Martin J.
2015-03-01
The numerical technique of lattice quantum chromodynamics (LQCD) holds the promise of connecting the nuclear forces, nuclei, the spectrum and structure of hadrons, and the properties of matter under extreme conditions with the underlying theory of the strong interactions, quantum chromodynamics. A distinguishing, and thus far unique, feature of this formulation is that all of the associated uncertainties, both statistical and systematic can, in principle, be systematically reduced to any desired precision with sufficient computational and human resources. We review the sources of uncertainty inherent in LQCD calculations for nuclear physics, and discuss how each is quantified in current efforts.
Variational perturbation theory and nonperturbative calculations in QCD
NASA Astrophysics Data System (ADS)
Solovtsova, O. P.
2013-10-01
A nonperturbative approach based on the variational perturbation theory in quantum chromodynamics is developed. The variational series is different from the conventional perturbative expansion and can be used to go beyond the weak-coupling regime. The approach suggested takes into account the summation of threshold singularities and the involvement of nonperturbative light quark masses. Phenomenological applications of this approach to describe physical quantities connected with the hadronic τ-decay data: the R τ ratio, the light-quark Adler function, and the smeared R Δ function are presented. The description of examined quantities includes an infrared region and, therefore, they cannot be directly calculated within the standard perturbation theory. It is shown that in spite of this fact the approach suggested gives a rather good result for these quantities down to the lowest energy scale.
Lattice QCD calculations of transverse momentum-dependent parton distributions (TMDs)
NASA Astrophysics Data System (ADS)
Engelhardt, M.; Musch, B.; Bhattacharya, T.; Green, J. R.; Gupta, R.; Hägler, P.; Krieg, S.; Negele, J.; Pochinsky, A.; Schäfer, A.; Syritsyn, S.; Yoon, B.
2016-03-01
An ongoing program of evaluating TMD observables within Lattice QCD is reviewed, summarizing recent progress with respect to several challenges faced by such calculations. These lattice calculations are based on a definition of TMDs through hadronic matrix elements of quark bilocal operators containing staple-shaped gauge connections. A parametrization of the matrix elements in terms of invariant amplitudes serves to cast them in the Lorentz frame preferred for a lattice calculation. Data on the naively T-odd Sivers and Boer-Mulders effects as well as the transversity TMD are presented.
Lattice QCD calculation of form factors describing the rare decays B→K*ℓ+ℓ- and Bs→ϕℓ+ℓ-
NASA Astrophysics Data System (ADS)
Horgan, Ronald R.; Liu, Zhaofeng; Meinel, Stefan; Wingate, Matthew
2014-05-01
The rare decays B0→K*0μ+μ- and Bs→ϕμ+μ- are now being observed with enough precision to test Standard Model predictions. A full understanding of these decays requires accurate determinations of the corresponding hadronic form factors. Here we present results of lattice QCD calculations of the B→K* and Bs→ϕ form factors. We also determine the form factors relevant for the decays Bs→K*ℓν and Bs→K¯*0ℓ+ℓ-. We use full-QCD configurations including 2+1 flavors of sea quarks using an improved staggered action, and we employ lattice nonrelativistic QCD to describe the bottom quark.
Lattice QCD for parallel computers
NASA Astrophysics Data System (ADS)
Quadling, Henley Sean
Lattice QCD is an important tool in the investigation of Quantum Chromodynamics (QCD). This is particularly true at lower energies where traditional perturbative techniques fail, and where other non-perturbative theoretical efforts are not entirely satisfactory. Important features of QCD such as confinement and the masses of the low lying hadronic states have been demonstrated and calculated in lattice QCD simulations. In calculations such as these, non-lattice techniques in QCD have failed. However, despite the incredible advances in computer technology, a full solution of lattice QCD may still be in the too-distant future. Much effort is being expended in the search for ways to reduce the computational burden so that an adequate solution of lattice QCD is possible in the near future. There has been considerable progress in recent years, especially in the research of improved lattice actions. In this thesis, a new approach to lattice QCD algorithms is introduced, which results in very significant efficiency improvements. The new approach is explained in detail, evaluated and verified by comparing physics results with current lattice QCD simulations. The new sub-lattice layout methodology has been specifically designed for current and future hardware. Together with concurrent research into improved lattice actions and more efficient numerical algorithms, the very significant efficiency improvements demonstrated in this thesis can play an important role in allowing lattice QCD researchers access to much more realistic simulations. The techniques presented in this thesis also allow ambitious QCD simulations to be performed on cheap clusters of commodity computers.
Thomas-Fermi Quark Model and Techniques to Improve Lattice QCD Calculation
NASA Astrophysics Data System (ADS)
Liu, Quan
Two topics are discussed separately in this thesis. In the first part a semiclassical quark model, called the Thomas-Fermi quark model, is reviewed. After a modified approach to spin in the model is introduced, I present the calculation of the spectra of octet and decuplet baryons. The six-quark doubly strange H-dibaryon state is also investigated. In the second part, two numerical techniques which improve latice QCD calculations are covered. The first one, which we call Polynomial-Preconditioned GMRES-DR(PP-GMRESDR), is used to speed up the calculation of large systems of linear equations in LQCD. The second one, called the Polynomial-Subtraction method, is used to help reduce the noise variance of the calculations for disconnected loops in LQCD.
Matching excluded-volume hadron-resonance gas models and perturbative QCD to lattice calculations
NASA Astrophysics Data System (ADS)
Albright, M.; Kapusta, J.; Young, C.
2014-08-01
We match three hadronic equations of state at low energy densities to a perturbatively computed equation of state of quarks and gluons at high energy densities. One of them includes all known hadrons treated as point particles, which approximates attractive interactions among hadrons. The other two include, in addition, repulsive interactions in the form of excluded volumes occupied by the hadrons. A switching function is employed to make the crossover transition from one phase to another without introducing a thermodynamic phase transition. A χ2 fit to accurate lattice calculations with temperature 100
Full orbit calculation for lost alpha particle measurement on ITER
Funaki, D.; Isobe, M.; Nishiura, M.; Sato, Y.; Okamoto, A.; Kobuchi, T.; Kitajima, S.; Sasao, M.
2008-10-15
An orbit following calculation code with full gyromotion under the ITER magnetic field configuration has been developed to investigate escaping alpha particle orbits in ITER and to determine the geometrical arrangement for alpha particle detection. The code contained the full geometrical information of the first wall panels. It was carefully investigated whether an alpha particle escaping from the plasma through the last closed flux surface does not touch or intersect the first wall boundary before reaching the detection point. Candidates of blanket module modification have been studied to achieve effective measurement geometry for escaping alpha particle detection. The calculations showed that direct orbit loss and banana diffusion can be detected with a probe head recessed from the first wall surface.
Full $\mathrm{QED}\mathbf{+}\mathrm{QCD}$ Low-Energy Constants through Reweighting
Ishikawa, Tomomi; Blum, Thomas; Hayakawa, Masashi; Izubuchi, Taku; Jung, Chulwoo; Zhou, Ran
2012-08-01
The effect of sea quark electromagnetic charge on meson masses is investigated, and first results for full QED+QCD low-energy constants are presented. The electromagnetic charge for sea quarks is incorporated in quenched QED+full QCD lattice simulations by a reweighting method. The reweighting factor, which connects quenched and unquenched QED, is estimated using a stochastic method on 2+1 flavor dynamical domain-wall quark ensembles.
MathQCDSR: a Mathematica Package for QCD Sum Rules Calculations
NASA Astrophysics Data System (ADS)
Wang, Lai; Lee, Frank X.
2010-02-01
We present a software package written in Mathematica for standard QCD sum rules calculations. Two examples are given to demonstrate how to use the package. One is for the mass spectrum of octet baryons from two-point correlation functions; the other for the magnetic moments of octet baryons in the external-field method. The software package FeynCalc is used to handle the gamma-matrix computations. The derived sum rules can be saved in various forms (native Mathematica, JAVA, Fortran, or Latex). Several notebooks are provided to carry out a Monte-Carlo-based numerical analysis of the sum rules, complete with in-line graphical display of sum rule matching, error distributions, scatter plots for correlations, etc. )
Quenched domain wall QCD with DBW2 gauge action toward nucleon decay matrix element calculation
NASA Astrophysics Data System (ADS)
Aoki, Yasumichi
2001-10-01
The domain wall fermion action is a promising way to control chiral symmetry in lattice gauge theory. By the good chiral symmetry of this approach even at finite lattice spacing, one is able to extract hadronic matrix elements, like kaon weak matrix elements, for which the symmetry is extremely important. Ordinary fermions with poor chiral symmetry make calculation difficult because of the large mixing of operators with different chiral structure. Even though the domain wall fermion action with the simple Wilson gauge action has a good chiral symmetry, one can further improve the symmetry by using a different gauge action. We take a non-perturbatively improved action, the DBW2 action of the QCD Taro group. Hadron masses are systematically examined for a range of parameters. Application to nucleon decay matrix element is also discussed.
FAPT: A Mathematica package for calculations in QCD Fractional Analytic Perturbation Theory
NASA Astrophysics Data System (ADS)
Bakulev, Alexander P.; Khandramai, Vyacheslav L.
2013-01-01
We provide here all the procedures in Mathematica which are needed for the computation of the analytic images of the strong coupling constant powers in Minkowski (A(s;nf) and Aνglob(s)) and Euclidean (A(Q2;nf) and Aνglob(Q2)) domains at arbitrary energy scales (s and Q2, correspondingly) for both schemes — with fixed number of active flavours nf=3,4,5,6 and the global one with taking into account all heavy-quark thresholds. These singularity-free couplings are inevitable elements of Analytic Perturbation Theory (APT) in QCD, proposed in [10,69,70], and its generalization — Fractional APT, suggested in [42,46,43], needed to apply the APT imperative for renormalization-group improved hadronic observables. Program summaryProgram title: FAPT Catalogue identifier: AENJ_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AENJ_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 1985 No. of bytes in distributed program, including test data, etc.: 1895776 Distribution format: tar.gz Programming language: Mathematica. Computer: Any work-station or PC where Mathematica is running. Operating system: Windows XP, Mathematica (versions 5 and 7). Classification: 11.5. Nature of problem: The values of analytic images A(Q2) and A(s) of the QCD running coupling powers αsν(Q2) in Euclidean and Minkowski regions, correspondingly, are determined through the spectral representation in the QCD Analytic Perturbation Theory (APT). In the program FAPT we collect all relevant formulas and various procedures which allow for a convenient evaluation of A(Q2) and A(s) using numerical integrations of the relevant spectral densities. Solution method: FAPT uses Mathematica functions to calculate different spectral densities and then performs numerical integration of these spectral integrals
A Molecular Full-Potential LMTO Calculation for Copper Clusters
NASA Astrophysics Data System (ADS)
Datta, Radhika Prosad; Banerjea, Amitava; Mookerjee, Abhijit; Bhattacharyya, A. K.
We study the electronic properties of small (10-20 atoms) copper clusters using the newly-developed molecular full-potential linearized muffin-tin orbital two-centre-fit (TCF) method of Methfessel and van Schilfgaarde. The geometric structures of the clusters had earlier been determined by us through simulated annealing using the Equivalent Crystal Theory to compute total energies. We report the variation of the binding energy, as obtained from the TCF calculations, with cluster size and compare these to the binding energies determined, for the same structures, from the ECT. We also show the variation of the HOMO-LUMO gap with cluster size, and the pseudo-density of states for select cluster sizes.
Full configuration interaction benchmark calculations for TiH
Bauschlicher, C.W. Jr.
1988-06-02
Full configuration interaction (FCI) calculations have been performed for the /sup 3/F and /sup 5/F states of Ti atom and the /sup 4/Phi and /sup 2/..delta.. states of TiH. The FCI calculations are compared to approximate treatments of the correlation problem; for the /sup 2/..delta.. state, the CASSCF/MRCI treatment agrees with the FCI results for r/sub e/, omega/sub e/, the dipole moment, and the dipole derivative. For the /sup 4/Phi state, the CASSCF/MRCI approach agrees well with the FCI for r/sub e/, omega/sub e/, and D/sub e/. However, the agreement between CASSCF/MRCI and FCI treatments for the dipole moment is not as good, even when the CASSCF and MRCI reference spaces contain up to 800 CSFs. Natural orbital iterations improve the dipole moment but have a smaller effect on the other spectroscopic parameters. The CPF and MCPF methods agree well with the FCI for the /sup 4/Phi state, which is reasonably well described by the SCF; this is true even for the dipole moment, where a natural orbital iteration must be performed for the CASSCF/MRCI treatment. The CPF and MCPF treatments do not agree as well for the /sup 2/..delta.. state, which is not as well described by the SCF.
Full waveform modelling and misfit calculation using the VERCE platform
NASA Astrophysics Data System (ADS)
Garth, Thomas; Spinuso, Alessandro; Casarotti, Emanuele; Magnoni, Federica; Krischner, Lion; Igel, Heiner; Schwichtenberg, Horst; Frank, Anton; Vilotte, Jean-Pierre; Rietbrock, Andreas
2016-04-01
simulated and recorded waveforms, enabling seismologists to specify and steer their misfit analyses using existing python tools and libraries such as Pyflex and the dispel4py data-intensive processing library. All these processes, including simulation, data access, pre-processing and misfit calculation, are presented to the users of the gateway as dedicated and interactive workspaces. The VERCE platform can also be used to produce animations of seismic wave propagation through the velocity model, and synthetic shake maps. We demonstrate the functionality of the VERCE platform with two case studies, using the pre-loaded velocity model and mesh for Chile and Northern Italy. It is envisioned that this tool will allow a much greater range of seismologists to access these full waveform inversion tools, and aid full waveform tomographic and source inversion, synthetic shake map production and other full waveform applications, in a wide range of tectonic settings.
Charmed spectroscopy from a nonperturbatively determined relativistic heavy quark action in full QCD
Huey-Wen Lin
2006-07-28
We present a preliminary calculation of the charmed meson spectrum using the 2+1 flavor domain wall fermion lattice configurations currently being generated by the RBC and UKQCD collaborations. The calculation is performed using the 3-parameter, relativistic heavy quark action with nonperturbatively determined coefficients. We will also demonstrate a step-scaling procedure for determining these coefficients nonperturbatively using a series of quenched, gauge field ensembles generated for three different lattice spacings.
Calculation of HELAS amplitudes for QCD processes using graphics processing unit (GPU)
NASA Astrophysics Data System (ADS)
Hagiwara, K.; Kanzaki, J.; Okamura, N.; Rainwater, D.; Stelzer, T.
2010-11-01
We use a graphics processing unit (GPU) for fast calculations of helicity amplitudes of quark and gluon scattering processes in massless QCD. New HEGET ( HELAS Evaluation with GPU Enhanced Technology) codes for gluon self-interactions are introduced, and a C++ program to convert the MadGraph generated FORTRAN codes into HEGET codes in CUDA (a C-platform for general purpose computing on GPU) is created. Because of the proliferation of the number of Feynman diagrams and the number of independent color amplitudes, the maximum number of final state jets we can evaluate on a GPU is limited to 4 for pure gluon processes ( gg→4 g), or 5 for processes with one or more quark lines such as qoverline{q}→ 5g and qq→ qq+3 g. Compared with the usual CPU-based programs, we obtain 60-100 times better performance on the GPU, except for 5-jet production processes and the gg→4 g processes for which the GPU gain over the CPU is about 20.
Calculation of the heavy-hadron axial couplings g1, g2, and g3 using lattice QCD
Will Detmold, David Lin, Stefan Meinel
2012-06-01
In a recent paper [arXiv:1109.2480] we have reported on a lattice QCD calculation of the heavy-hadron axial couplings g{sub 1}, g{sub 2}, and g{sub 3}. These quantities are low-energy constants of heavy-hadron chiral perturbation theory (HH{chi}PT) and are related to the B*B{pi}, {Sigma}{sub b}*{Sigma}{sub b}{pi}, and {Sigma}{sub b}{sup (*)}{Lambda}{sub b}{pi} couplings. In the following, we discuss important details of the calculation and give further results. To determine the axial couplings, we explicitly match the matrix elements of the axial current in QCD with the corresponding matrix elements in HH{chi}PT. We construct the ratios of correlation functions used to calculate the matrix elements in lattice QCD, and study the contributions from excited states. We present the complete numerical results and discuss the data analysis in depth. In particular, we demonstrate the convergence of SU(4|2) HH{chi}PT for the axial current matrix elements at pion masses up to about 400 MeV and show the impact of the nonanalytic loop contributions. Finally, we present additional predictions for strong and radiative decay widths of charm and bottom baryons.
fK /f{pi} in Full QCD with Domain Wall Valence Quarks
Silas Beane; Paulo Bedaque; Konstantinos Orginos; Martin Savage
2007-05-01
We compute the ratio of pseudoscalar decay constants f{sub K}/f{sub {pi}} using domain-wall valence quarks and rooted improved Kogut-Susskind sea quarks. By employing continuum chiral perturbation theory, we extract the Gasser-Leutwyler low-energy constant L{sub 5}, and extrapolate f{sub K}/f{sub {pi}} to the physical point. We find: f{sub K}/f{sub {pi}} = 1.218 {+-} 0.002{sub -0.024}{sup +0.011} where the first error is statistical and the second error is an estimate of the systematic due to chiral extrapolation and fitting procedures. This value agrees within the uncertainties with the determination by the MILC collaboration, calculated using Kogut-Susskind valence quarks, indicating that systematic errors arising from the choice of lattice valence quark are small.
A full variational calculation based on a tensor product decomposition
NASA Astrophysics Data System (ADS)
Senese, Frederick A.; Beattie, Christopher A.; Schug, John C.; Viers, Jimmy W.; Watson, Layne T.
1989-08-01
A new direct full variational approach exploits a tensor (Kronecker) product decomposition of the Hamiltonian. Explicit assembly and storage of the Hamiltonian matrix is avoided by using the Kronecker product structure to form matrix-vector products directly from the molecular integrals. Computation-intensive integral transformations and formula tapes are unnecessary. The wavefunction is expanded in terms of spin-free primitive kets rather than Slater determinants or configuration state functions, and the expansion is equivalent to a full configuration interaction expansion. The approach suggests compact storage schemes and algorithms which are naturally suited to parallel and pipelined machines.
A QCD sum rule calculation of the X± (5568) → Bs0 π± decay width
NASA Astrophysics Data System (ADS)
Dias, J. M.; Khemchandani, K. P.; Martínez Torres, A.; Nielsen, M.; Zanetti, C. M.
2016-07-01
To understand the nature of the X (5568), recently observed in the mass spectrum of the Bs0 π± system by the D0 Collaboration, we have investigated, in a previous work, a scalar tetraquark (diquak-antidiquark) structure for it, within the two-point QCD sum rules method. We found that it is possible to obtain a stable value of the mass compatible with the D0 result, although a rigorous QCD sum rule constrained analysis led to a higher value of mass. As a continuation of our investigation, we calculate the width of the tetraquark state with same quark content as X (5568), to the channel Bs0 π±, using the three-point QCD sum rule. We obtain a value of (20.4 ± 8.7) MeV for the mass ∼ 5568 MeV, which is compatible with the experimental value of 21.9 ± 6.4 (sta)-2.5+5.0 (syst) MeV /c2. We find that the decay width to Bs0 π± does not alter much for a higher mass state.
Full configuration interaction benchmark calculations for TiH
NASA Technical Reports Server (NTRS)
Bauschlicher, Charles W., Jr.
1988-01-01
Full configuration interaction (FCI) results for the 3F and 5F states of Ti and the 4Phi and 2Delta states of TiH are presented. While the coupled pair functional (CPF) or modified CPF approaches are found to work well for the 4Phi state of TiH, they do not perform as well for the 2Delta state. Although for mu, the CASSCF/MRCI methods do well for the 2Delta state, when the active space is expanded outside the normal valence definition, the dipole moment is only brought into agreement with the FCI results by natural orbital iterations.
Accelerating Full Configuration Interaction Calculations for Nuclear Structure
Yang, Chao; Sternberg, Philip; Maris, Pieter; Ng, Esmond; Sosonkina, Masha; Le, Hung Viet; Vary, James; Yang, Chao
2008-04-14
One of the emerging computational approaches in nuclear physics is the full configuration interaction (FCI) method for solving the many-body nuclear Hamiltonian in a sufficiently large single-particle basis space to obtain exact answers - either directly or by extrapolation. The lowest eigenvalues and correspondingeigenvectors for very large, sparse and unstructured nuclear Hamiltonian matrices are obtained and used to evaluate additional experimental quantities. These matrices pose a significant challenge to the design and implementation of efficient and scalable algorithms for obtaining solutions on massively parallel computer systems. In this paper, we describe the computational strategies employed in a state-of-the-art FCI code MFDn (Many Fermion Dynamics - nuclear) as well as techniques we recently developed to enhance the computational efficiency of MFDn. We will demonstrate the current capability of MFDn and report the latest performance improvement we have achieved. We will also outline our future research directions.
NASA Astrophysics Data System (ADS)
Vitev, I.
We review recent developments in the QCD description of jet production and modification in reactions with heavy nuclei at relativistic energies. Our goal is to formulate a perturbative expansion in the presence of nuclear matter that allows to systematically improve the accuracy of the theoretical predictions. As an example, we present calculations of inclusive jet cross sections at RHIC, Z(0/gamma^*) -tagged jet cross sections at the LHC, and jet shapes that include both next-to-leading order perturbative effects and the effects of the nuclear medium.
QCD sum rule calculation of quark-gluon three-body components in the B-meson wave function
Nishikawa, Tetsuo; Tanaka, Kazuhiro
2011-10-21
We discuss the QCD sum rule calculation of the heavy-quark effective theory parameters {lambda}{sub E} and {lambda}{sub H}, which represent quark-gluon three-body components in the B-meson wave function. We update the sum rules for {lambda}{sub E,H} calculating the new higher-order contributions to the operator product expansion for the corresponding correlator, i.e., the order {alpha}{sub s} radiative corrections to the Wilson coefficients associated with the dimension-5 quark-gluon mixed condensate, and the power corrections due to the dimension-6 vacuum condensates. We find that the new radiative corrections significantly improve stability of the corresponding Borel sum rules, modifying the values of {lambda}{sub E,H}.
Lattice QCD in rotating frames.
Yamamoto, Arata; Hirono, Yuji
2013-08-23
We formulate lattice QCD in rotating frames to study the physics of QCD matter under rotation. We construct the lattice QCD action with the rotational metric and apply it to the Monte Carlo simulation. As the first application, we calculate the angular momenta of gluons and quarks in the rotating QCD vacuum. This new framework is useful to analyze various rotation-related phenomena in QCD. PMID:24010426
A full CI treatment of Ne atom - A benchmark calculation performed on the NAS CRAY 2
NASA Technical Reports Server (NTRS)
Bauschlicher, C. W., Jr.; Langhoff, S. R.; Partridge, H.; Taylor, P. R.
1986-01-01
Full CI calculations are performed for Ne atom using Gaussian basis sets of up to triple-zeta plus double polarization quality. The total valence correlation energy through double, triple, quadruple and octuple excitations is compared for eight different basis sets. These results are expected to be an important benchmark for calibrating methods for estimating the importance of higher excitations.
Harrison, R.J.; Stahlberg, E.A.
1994-10-01
We describe an implementation of the benchmark ab initio electronic structure full configuration interaction model on the Intel Touchstone Delta. Its performance is demonstrated with several calculations, the largest of which (95 million configurations, 418 million determinants) is the largest full-CI calculation yet completed. The feasibility of calculations with over one billion configurations is discussed. A sustained computation rate in excess of 4 GFLOP/s on 512 processors is achieved, with an average aggregate communication rate of 155 Mbytes/s. Data-compression techniques and a modified diagonalization method were required to minimize I/O. The object-oriented design has increased portability and provides the distinction between local and non-local data essential for use of a distributed-data model.
Hadronic Resonances from Lattice QCD
Lichtl, Adam C.; Bulava, John; Morningstar, Colin; Edwards, Robert; Mathur, Nilmani; Richards, David; Fleming, George; Juge, K. Jimmy; Wallace, Stephen J.
2007-10-26
The determination of the pattern of hadronic resonances as predicted by Quantum Chromodynamics requires the use of non-perturbative techniques. Lattice QCD has emerged as the dominant tool for such calculations, and has produced many QCD predictions which can be directly compared to experiment. The concepts underlying lattice QCD are outlined, methods for calculating excited states are discussed, and results from an exploratory Nucleon and Delta baryon spectrum study are presented.
Hadronic Resonances from Lattice QCD
John Bulava; Robert Edwards; George Fleming; K. Jimmy Juge; Adam C. Lichtl; Nilmani Mathur; Colin Morningstar; David Richards; Stephen J. Wallace
2007-06-16
The determination of the pattern of hadronic resonances as predicted by Quantum Chromodynamics requires the use of non-perturbative techniques. Lattice QCD has emerged as the dominant tool for such calculations, and has produced many QCD predictions which can be directly compared to experiment. The concepts underlying lattice QCD are outlined, methods for calculating excited states are discussed, and results from an exploratory Nucleon and Delta baryon spectrum study are presented.
NASA Astrophysics Data System (ADS)
Taguchi, Masayoshi
1982-05-01
The neoclassical transport coefficients from the plateau regime to the Pfirsch-Schlüter regime are calculated by using the linearized full Fokker-Planck collision operator. The inverse aspect ratio is assumed to be small and the mass ratio of the electron to the ion is much less than unity. The lowest-order terms of these ratios are retained in this calculation. The results obtained give correction to those of Rawls et al. who employed a model collision operator. In the Pfirsch-Schlüter regime, the results of this paper agree fairly well with those given by Hazeltine and Hinton.
3D Neutron Transport PWR Full-core Calculation with RMC code
NASA Astrophysics Data System (ADS)
Qiu, Yishu; She, Ding; Fan, Xiao; Wang, Kan; Li, Zeguang; Liang, Jingang; Leroyer, Hadrien
2014-06-01
Nowadays, there are more and more interests in the use of Monte Carlo codes to calculate the detailed power density distributions in full-core reactors. With the Inspur TS1000 HPC Server of Tsinghua University, several calculations have been done based on the EDF 3D Neutron Transport PWR Full-core benchmark through large-scale parallelism. To investigate and compare the results of the deterministic method and Monte Carlo method, EDF R&D and Department of Engineering Physics of Tsinghua University are having a collaboration to make code to code verification. So in this paper, two codes are used. One is the code COCAGNE developed by the EDF R&D, a deterministic core code, and the other is the Monte Carlo code RMC developed by Department of Engineering Physics in Tsinghua University. First, the full-core model is described and a 26-group calculation was performed by these two codes using the same 26-group cross-section library provided by EDF R&D. Then the parallel and tally performance of RMC is discussed. RMC employs a novel algorithm which can cut down most of the communications. It can be seen clearly that the speedup ratio almost linearly increases with the nodes. Furthermore the cell-mapping method applied by RMC consumes little time to tally even millions of cells. The results of the codes COCAGNE and RMC are compared in three ways. The results of these two codes agree well with each other. It can be concluded that both COCAGNE and RMC are able to provide 3D-transport solutions associated with detailed power density distributions calculation in PWR full-core reactors. Finally, to investigate how many histories are needed to obtain a given standard deviation for a full 3D solution, the non-symmetrized condensed 2-group fluxes of RMC are discussed.
NASA Technical Reports Server (NTRS)
Shertzer, Janine; Temkin, Aaron
2004-01-01
The development of a practical method of accurately calculating the full scattering amplitude, without making a partial wave decomposition is continued. The method is developed in the context of electron-hydrogen scattering, and here exchange is dealt with by considering e-H scattering in the static exchange approximation. The Schroedinger equation in this approximation can be simplified to a set of coupled integro-differential equations. The equations are solved numerically for the full scattering wave function. The scattering amplitude can most accurately be calculated from an integral expression for the amplitude; that integral can be formally simplified, and then evaluated using the numerically determined wave function. The results are essentially identical to converged partial wave results.
Skands, Peter Z.; /Fermilab
2005-07-01
Recent developments in QCD phenomenology have spurred on several improved approaches to Monte Carlo event generation, relative to the post-LEP state of the art. In this brief review, the emphasis is placed on approaches for (1) consistently merging fixed-order matrix element calculations with parton shower descriptions of QCD radiation, (2) improving the parton shower algorithms themselves, and (3) improving the description of the underlying event in hadron collisions.
Status of lattice field theory calculations
Sharpe, S.R.
1990-01-01
This report briefly discusses the following topics: overview of all present calculation; reliability criteria for quenched calculation; quenched versus full QCD, and difficulties facing full QCD; results for the quenched pion wavefunction''; results for the quenched hadron spectrum; results for quenched B{sub K}; A new method for calculating the surface tension; the non-pertubative upper bound on the Higgs mass; and toward the TERAFLOP machine.
Full circuit calculation for electromagnetic pulse transmission in a high current facility
NASA Astrophysics Data System (ADS)
Zou, Wenkang; Guo, Fan; Chen, Lin; Song, Shengyi; Wang, Meng; Xie, Weiping; Deng, Jianjun
2014-11-01
We describe herein for the first time a full circuit model for electromagnetic pulse transmission in the Primary Test Stand (PTS)—the first TW class pulsed power driver in China. The PTS is designed to generate 8-10 MA current into a z -pinch load in nearly 90 ns rise time for inertial confinement fusion and other high energy density physics research. The PTS facility has four conical magnetic insulation transmission lines, in which electron current loss exists during the establishment of magnetic insulation. At the same time, equivalent resistance of switches and equivalent inductance of pinch changes with time. However, none of these models are included in a commercially developed circuit code so far. Therefore, in order to characterize the electromagnetic transmission process in the PTS, a full circuit model, in which switch resistance, magnetic insulation transmission line current loss and a time-dependent load can be taken into account, was developed. Circuit topology and an equivalent circuit model of the facility were introduced. Pulse transmission calculation of shot 0057 was demonstrated with the corresponding code FAST (full-circuit analysis and simulation tool) by setting controllable parameters the same as in the experiment. Preliminary full circuit simulation results for electromagnetic pulse transmission to the load are presented. Although divergences exist between calculated and experimentally obtained waveforms before the vacuum section, consistency with load current is satisfactory, especially at the rising edge.
Equation of state in two-, three-, and four-color QCD at nonzero temperature and density
NASA Astrophysics Data System (ADS)
Gorda, Tyler; Romatschke, Paul
2015-07-01
We calculate the equation of state at nonzero temperature and density from first principles in two-, three-, and four-color QCD with two fermion flavors in the fundamental and two-index, antisymmetric representation. By matching low-energy results (from a "hadron resonance gas") to high-energy results from (resummed) perturbative QCD, we obtain results for the pressure and trace anomaly that are in quantitative agreement with full lattice-QCD studies for three colors at zero chemical potential. Our results for nonzero chemical potential at zero temperature constitute predictions for the equation of state in QCD-like theories that can be tested by traditional lattice studies for two-color QCD with two fundamental fermions and four-color QCD with two two-index, antisymmetric fermions. We find that the speed of sound squared at zero temperature can exceed 1 /3 , which may be relevant for the phenomenology of high-mass neutron stars.
Constantia Alexandrou; Bojan Bistrovic; Robert Edwards; P de Forcrand; George Fleming; Philipp Haegler; John Negele; Konstantinos Orginos; Andrew Pochinsky; Dru Renner; David Richards; Wolfram Schroers; Antonios Tsapalis
2005-10-01
Lattice QCD is an essential complement to the current and anticipated DOE-supported experimental program in hadronic physics. In this poster we address several key questions central to our understanding of the building blocks of nuclear matter, nucleons and pions. Firstly, we describe progress at computing the electromagnetic form factors of the nucleon, describing the distribution of charge and current, before considering the role played by the strange quarks. We then describe the study of transition form factors to the Delta resonance. Finally, we present recent work to determine the pion form factor, complementary to the current JLab experimental determination and providing insight into the approach to asymptotic freedom.
Nakatsuji, Hiroshi
2011-12-15
The simplest iterative complement (SIC) calculations starting from Hartree-Fock and giving full configuration interaction (CI) at convergence were performed using regular and inverse Hamiltonians. Each iteration step is variational and involves only one variable. The convergence was slow when we used the regular Hamiltonian, but became very fast when we used the inverse Hamiltonian. This difference is due to the Coulomb singularity problem inherent in the regular Hamiltonian; the inverse Hamiltonian does not have such a problem. For this reason, the merit of the inverse Hamiltonian over the regular one becomes even more dramatic when we use a better-quality basis set. This was seen by comparing the calculations due to the minimal and double-{zeta} basis sets. Similar problematic situations exist in the Krylov sequence and in the Lanczos and Arnoldi methods.
Statistical error propagation in ab initio no-core full configuration calculations of light nuclei
NASA Astrophysics Data System (ADS)
Navarro Pérez, R.; Amaro, J. E.; Ruiz Arriola, E.; Maris, P.; Vary, J. P.
2015-12-01
We propagate the statistical uncertainty of experimental N N scattering data into the binding energy of 3H and 4He. We also study the sensitivity of the magnetic moment and proton radius of the 3H to changes in the N N interaction. The calculations are made with the no-core full configuration method in a sufficiently large harmonic oscillator basis. For those light nuclei we obtain Δ Estat(3H) =0.015 MeV and Δ Estat(4He) =0.055 MeV .
Full CI benchmark calculations for several states of the same symmetry
NASA Technical Reports Server (NTRS)
Bauschlicher, Charles W., Jr.; Taylor, Peter R.
1987-01-01
Full CI (FCI) wave functions are used to compute energies for several electronic states of the same symmetry for SiH2, CH2, and CH2(+). It is found that CASSCF/multireference CI wave functions yield results very similar to FCI, irrespective of whether the CASSCF MOs are optimized independently for each state or using an average of the CASSCF energies for all desired states. The ionization potentials and excitation energies obtained from the FCI calculations should help calibrate methods (such as Green's function approaches, equations of motion and propagator methods, and cluster expansions) in which energy differences are computed directly.
On the importance of full-dimensionality in low-energy molecular scattering calculations
Faure, Alexandre; Jankowski, Piotr; Stoecklin, Thierry; Szalewicz, Krzysztof
2016-01-01
Scattering of H2 on CO is of great importance in astrophysics and also is a benchmark system for comparing theory to experiment. We present here a new 6-dimensional potential energy surface for the ground electronic state of H2-CO with an estimated uncertainty of about 0.6 cm−1 in the global minimum region, several times smaller than achieved earlier. This potential has been used in nearly exact 6-dimensional quantum scattering calculations to compute state-to-state cross-sections measured in low-energy crossed-beam experiments. Excellent agreement between theory and experiment has been achieved in all cases. We also show that the fully 6-dimensional approach is not needed with the current accuracy of experimental data since an equally good agreement with experiment was obtained using only a 4-dimensional treatment, which validates the rigid-rotor approach widely used in scattering calculations. This finding, which disagrees with some literature statements, is important since for larger systems full-dimensional scattering calculations are currently not possible. PMID:27333870
On the importance of full-dimensionality in low-energy molecular scattering calculations.
Faure, Alexandre; Jankowski, Piotr; Stoecklin, Thierry; Szalewicz, Krzysztof
2016-01-01
Scattering of H2 on CO is of great importance in astrophysics and also is a benchmark system for comparing theory to experiment. We present here a new 6-dimensional potential energy surface for the ground electronic state of H2-CO with an estimated uncertainty of about 0.6 cm(-1) in the global minimum region, several times smaller than achieved earlier. This potential has been used in nearly exact 6-dimensional quantum scattering calculations to compute state-to-state cross-sections measured in low-energy crossed-beam experiments. Excellent agreement between theory and experiment has been achieved in all cases. We also show that the fully 6-dimensional approach is not needed with the current accuracy of experimental data since an equally good agreement with experiment was obtained using only a 4-dimensional treatment, which validates the rigid-rotor approach widely used in scattering calculations. This finding, which disagrees with some literature statements, is important since for larger systems full-dimensional scattering calculations are currently not possible. PMID:27333870
NASA Astrophysics Data System (ADS)
Ravindran, P.; Kjekshus, A.; Fjellvåg, H.; James, P.; Nordström, L.; Johansson, B.; Eriksson, O.
2001-04-01
The computational framework of this study is based on the local-spin-density approximation with first-principles full-potential linear muffin-tin orbital calculations including orbital polarization (OP) correction. We have studied the magnetic anisotropy for a series of bilayer CuAu(I)-type materials such as FeX, MnX (X=Ni,Pd,Pt), CoPt, NiPt, MnHg, and MnRh in a ferromagnetic state using experimental structural parameters to understand the microscopic origin of magnetic-anisotropy energy (MAE) in magnetic multilayers. Except for MnRh and MnHg, all these phases show perpendicular magnetization. We have analyzed our results in terms of angular momentum-, spin- and site-projected density of states, magnetic-angular-momentum-projected density of states, orbital-moment density of states, and total density of states. The orbital-moment number of states and the orbital-moment anisotropy for FeX (X=Ni,Pd,Pt) are calculated as a function of band filling to study its effect on MAE. The total and site-projected spin and orbital moments for all these systems are calculated with and without OP when the magnetization is along or perpendicular to the plane. The results are compared with available experimental as well as theoretical results. Our calculations show that OP always enhances the orbital moment in these phases and brings them closer to experimental values. The changes in MAE are analyzed in terms of exchange splitting, spin-orbit splitting, and tetragonal distortion/crystal-field splitting. The calculated MAE is found to be in good agreement with experimental values when the OP correction is included. Some of the materials considered here show large magnetic anisotropy of the order of meV. In particular we found that MnPt will have a very large MAE if it could be stabilized in a ferromagnetic configuration. Our analysis indicates that apart from large spin-orbit interaction and exchange interaction from at least one of the constituents, a large crystal-field splitting
QCD: Questions, challenges, and dilemmas
Bjorken, J.
1996-11-01
An introduction to some outstanding issues in QCD is presented, with emphasis on work by Diakonov and co-workers on the influence of the instanton vacuum on low-energy QCD observables. This includes the calculation of input valence-parton distributions for deep-inelastic scattering. 35 refs., 3 figs.
Full band calculations of the intrinsic lower limit of contact resistivity
NASA Astrophysics Data System (ADS)
Maassen, J.; Jeong, C.; Baraskar, A.; Rodwell, M.; Lundstrom, M.
2013-03-01
The intrinsic lower limit of contact resistivity (ρcLL) for InAs, In0.53Ga0.47As, GaSb, and Si is calculated using a full band ballistic quantum transport approach. Surprisingly, our results show that ρcLL is almost independent of the semiconductor. An analytical model, derived for 1D, 2D, and 3D, correctly reproduces the numerical results and explains why ρcLL is very similar in all cases. Our analysis sets a minimal carrier density required to meet the International Technology Roadmap for Semiconductors call for ρc=10-9 Ω-cm2 by 2023. Comparison with experiments shows there is room for improvement, which will come from optimizing interfacial properties.
Calculation of the Full Scattering Amplitude without Partial Wave Decomposition II
NASA Technical Reports Server (NTRS)
Shertzer, J.; Temkin, A.
2003-01-01
As is well known, the full scattering amplitude can be expressed as an integral involving the complete scattering wave function. We have shown that the integral can be simplified and used in a practical way. Initial application to electron-hydrogen scattering without exchange was highly successful. The Schrodinger equation (SE) can be reduced to a 2d partial differential equation (pde), and was solved using the finite element method. We have now included exchange by solving the resultant SE, in the static exchange approximation. The resultant equation can be reduced to a pair of coupled pde's, to which the finite element method can still be applied. The resultant scattering amplitudes, both singlet and triplet, as a function of angle can be calculated for various energies. The results are in excellent agreement with converged partial wave results.
NASA Technical Reports Server (NTRS)
Shertzer, Janine; Temkin, A.
2003-01-01
As is well known, the full scattering amplitude can be expressed as an integral involving the complete scattering wave function. We have shown that the integral can be simplified and used in a practical way. Initial application to electron-hydrogen scattering without exchange was highly successful. The Schrodinger equation (SE), which can be reduced to a 2d partial differential equation (pde), was solved using the finite element method. We have now included exchange by solving the resultant SE, in the static exchange approximation, which is reducible to a pair of coupled pde's. The resultant scattering amplitudes, both singlet and triplet, calculated as a function of energy are in excellent agreement with converged partial wave results.
Godunov-type schemes with an inertia term for unsteady full Mach number range flow calculations
NASA Astrophysics Data System (ADS)
Moguen, Yann; Delmas, Simon; Perrier, Vincent; Bruel, Pascal; Dick, Erik
2015-01-01
An inertia term is introduced in the AUSM+-up scheme. The resulting scheme, called AUSM-IT (IT for Inertia Term), is designed as an extension of the AUSM+-up scheme allowing for full Mach number range calculations of unsteady flows including acoustic features. In line with the continuous asymptotic analysis, the AUSM-IT scheme satisfies the conservation of the discrete linear acoustic energy at first order in the low Mach number limit. Its capability to properly handle low Mach number unsteady flows, that may include acoustic waves or discontinuities, is numerically illustrated. The approach for building the AUSM-IT scheme from the AUSM+-up scheme is applicable to any other Godunov-type scheme.
Transonic flow analysis for rotors. Part 2: Three-dimensional, unsteady, full-potential calculation
NASA Technical Reports Server (NTRS)
Chang, I. C.
1985-01-01
A numerical method is presented for calculating the three-dimensional unsteady, transonic flow past a helicopter rotor blade of arbitrary geometry. The method solves the full-potential equations in a blade-fixed frame of reference by a time-marching implicit scheme. At the far-field, a set of first-order radiation conditions is imposed, thus minimizing the reflection of outgoing wavelets from computational boundaries. Computed results are presented to highlight radial flow effects in three dimensions, to compare surface pressure distributions to quasi-steady predictions, and to predict the flow field on a swept-tip blade. The results agree well with experimental data for both straight- and swept-tip blade geometries.
Hybrid parallel code acceleration methods in full-core reactor physics calculations
Courau, T.; Plagne, L.; Ponicot, A.; Sjoden, G.
2012-07-01
When dealing with nuclear reactor calculation schemes, the need for three dimensional (3D) transport-based reference solutions is essential for both validation and optimization purposes. Considering a benchmark problem, this work investigates the potential of discrete ordinates (Sn) transport methods applied to 3D pressurized water reactor (PWR) full-core calculations. First, the benchmark problem is described. It involves a pin-by-pin description of a 3D PWR first core, and uses a 8-group cross-section library prepared with the DRAGON cell code. Then, a convergence analysis is performed using the PENTRAN parallel Sn Cartesian code. It discusses the spatial refinement and the associated angular quadrature required to properly describe the problem physics. It also shows that initializing the Sn solution with the EDF SPN solver COCAGNE reduces the number of iterations required to converge by nearly a factor of 6. Using a best estimate model, PENTRAN results are then compared to multigroup Monte Carlo results obtained with the MCNP5 code. Good consistency is observed between the two methods (Sn and Monte Carlo), with discrepancies that are less than 25 pcm for the k{sub eff}, and less than 2.1% and 1.6% for the flux at the pin-cell level and for the pin-power distribution, respectively. (authors)
J.W. Negele; R.C. Brower; P. Dreher; R. Edwards; G. Fleming; Ph. Hagler; U.M. Heller; Th. Lippert; A.V.Pochinsky; D.B. Renner; D. Richards; K. Schilling; W. Schroers
2004-04-01
This talk presents recent calculations in full QCD of the lowest three moments of generalized parton distributions and the insight they provide into the behavior of nucleon electromagnetic form factors, the origin of the nucleon spin, and the transverse structure of the nucleon. In addition, new exploratory calculations in the chiral regime of full QCD are discussed.
Torok, Aaron
2011-10-24
The {pi}{sup +}{Sigma}{sup +} and {pi}{sup +}{Xi}{sup 0} scattering lengths were calculated in mixed-action Lattice QCD with domain-wall valence quarks on the asqtad-improved coarse MILC configurations at four light-quark masses, and at two light-quark masses on the fine MILC configurations. Heavy Baryon Chiral Perturbation Theory with two and three flavors of light quarks was used to perform the chiral extrapolations. To NNLO in the three-flavor chiral expansion, the kaon-baryon processes that were investigated show no signs of convergence. Using the two-flavor chiral expansion for extrapolation, the pion-hyperon scattering lengths are found to be a{sub {pi}}{sup +}{sub {Sigma}}{sup +} = -0.197{+-}0.017 fm, and a{sub {pi}}{sup +}{sub {Xi}}{sup 0} = -0.098{+-}0.017 fm, where the comprehensive error includes statistical and systematic uncertainties.
Full potential calculation of electronics and thermoelectric properties of doped Mg{sub 2}Si
Poopanya, P.; Yangthaisong, A.
2013-12-04
We present the calculations of the electronic structure and transport properties on the anti-fluorite Mg{sub 2}Si using the full potential linearized augmented plane-wave (FP-LAPW) method and the semi-classical Boltzmann theory. The modified Becke-Johnson (mBJ) exchange potentials are used to derive energy gaps and correct band gaps according to experimental values. It is found that Mg{sub 2}Si is an indirect band gap (Γ→X) material with the gap of 0.56 eV which is in good agreement with the experimental observation. Note that the band structure of Mg{sub 2}Si is directly used in combination with the semi-classical Boltzmann theory to obtain the transport coefficients. It is found that the material is the n-type semiconductor with the lowest electron concentration of 3.03×10{sup 14} cm{sup −3} at 300 K. We have also calculated the thermoelectric properties of Mg{sub 2}Si based on the rigid band approximation by varying the p-type and n-type doping levels. At room temperature, the highest power factor for p-type and n-type dopants are obtained at the hole and electron concentration of 1.63×10{sup 20} cm{sup −3} and 1.15×1021 cm{sup −3}, respectively. From the electronic states, we also found that the n-type doping region is dominated by the Mg−2p{sup 6} 3s{sup 2} and Si−3p{sup 2} states, while the Mg−2p{sup 6} and Si−3p{sup 2} states are important in the p-type doped Mg{sub 2}Si.
Full potential calculations on the electron bandstructures of Sphalerite, Pyrite and Chalcopyrite
NASA Astrophysics Data System (ADS)
Edelbro, R.; Sandström, Å.; Paul, J.
2003-02-01
The bulk electronic structures of Sphalerite, Pyrite and Chalcopyrite have been calculated within an ab initio, full potential, density functional approach. The exchange term was approximated with the Dirac exchange functional, the Vosko-Wilk-Nusair parameterization of the Cepler-Alder free electron gas was used for correlation and linear combinations of Gaussian type orbitals were used as basis functions. The Sphalerite (zinc blende) band gap was calculated to be direct with a width of 2.23 eV. The Sphalerite valence band was 5.2 eV wide and composed of a mixture of sulfur and zinc orbitals. The band below the valence band located around -6.2 eV was mainly composed of Zn 3d orbitals. The S 3s orbitals gave rise to a band located around -12.3 eV. Pyrite was calculated to be a semiconductor with an indirect band gap of 0.51 eV, and a direct gap of 0.55 eV. The valence band was 1.25 eV wide and mainly composed of non-bonding Fe 3d orbitals. The band below the valence band was 4.9 eV wide and composed of a mixture of sulfur and iron orbitals. Due to the short inter-atomic distance between the sulfur dumbbells, the S 3s orbitals in Pyrite were split into a bonding and an anti-bonding range. Chalcopyrite was predicted to be a conductor, with no band-crossings at the Fermi level. The bands at -13.2 eV originate from the sulfur 3s orbitals and were quite similar to the sulfur 3s bands in Sphalerite, though somewhat shifted to lower energy. The top of the valence band consisted of a mixture of orbitals from all the atoms. The lower part of the same band showed metal character. Computational modeling as a tool for illuminating the flotation and leaching processes of Pyrite and Chalcopyrite, in connection with surface science experiments, is discussed.
NASA Astrophysics Data System (ADS)
Sboev, A. G.; Ilyashenko, A. S.; Vetrova, O. A.
1997-02-01
The method of bucking evaluation, realized in the MOnte Carlo code MCS, is described. This method was applied for calculational analysis of well known light water experiments TRX-1 and TRX-2. The analysis of this comparison shows, that there is no coincidence between Monte Carlo calculations, obtained by different ways: the MCS calculations with given experimental bucklings; the MCS calculations with given bucklings evaluated on base of full core MCS direct simulations; the full core MCNP and MCS direct simulations; the MCNP and MCS calculations, where the results of cell calculations are corrected by the coefficients taking into the account the leakage from the core. Also the buckling values evaluated by full core MCS calculations have differed from experimental ones, especially in the case of TRX-1, when this difference has corresponded to 0.5 percent increase of Keff value.
Plunkett, R.; The CDF Collaboration
1991-10-01
Results are presented for hadronic jet and direct photon production at {radical}{bar s} = 1800 GeV. The data are compared with next-to-leading QCD calculations. A new limit on the scale of possible composite structure of the quarks is also reported. 12 refs., 4 figs.
Andreas S. Kronfeld
2002-09-30
After reviewing some of the mathematical foundations and numerical difficulties facing lattice QCD, I review the status of several calculations relevant to experimental high-energy physics. The topics considered are moments of structure functions, which may prove relevant to search for new phenomena at the LHC, and several aspects of flavor physics, which are relevant to understanding CP and flavor violation.
NASA Astrophysics Data System (ADS)
Majumder, Moumita; Dawes, Richard; Wang, Xiao-Gang; Carrington, Tucker; Li, Jun; Guo, Hua; Manzhos, Sergei
2014-06-01
New potential energy surfaces for methane were constructed, represented as analytic fits to about 100,000 individual high-level ab initio data. Explicitly-correlated multireference data (MRCI-F12(AE)/CVQZ-F12) were computed using Molpro [1] and fit using multiple strategies. Fits with small to negligible errors were obtained using adaptations of the permutation-invariant-polynomials (PIP) approach [2,3] based on neural-networks (PIP-NN) [4,5] and the interpolative moving least squares (IMLS) fitting method [6] (PIP-IMLS). The PESs were used in full-dimensional vibrational calculations with an exact kinetic energy operator by representing the Hamiltonian in a basis of products of contracted bend and stretch functions and using a symmetry adapted Lanczos method to obtain eigenvalues and eigenvectors. Very close agreement with experiment was produced from the purely ab initio PESs. References 1- H.-J. Werner, P. J. Knowles, G. Knizia, 2012.1 ed. 2012, MOLPRO, a package of ab initio programs. see http://www.molpro.net. 2- Z. Xie and J. M. Bowman, J. Chem. Theory Comput 6, 26, 2010. 3- B. J. Braams and J. M. Bowman, Int. Rev. Phys. Chem. 28, 577, 2009. 4- J. Li, B. Jiang and Hua Guo, J. Chem. Phys. 139, 204103 (2013). 5- S Manzhos, X Wang, R Dawes and T Carrington, JPC A 110, 5295 (2006). 6- R. Dawes, X-G Wang, A.W. Jasper and T. Carrington Jr., J. Chem. Phys. 133, 134304 (2010).
Nucleon Structure from Dynamical Lattice QCD
Huey-Wen Lin
2007-06-01
We present lattice QCD numerical calculations of hadronic structure functions and form factors from full-QCD lattices, with a chirally symmetric fermion action, domain-wall fermions, for the sea and valence quarks. The lattice spacing is about 0.12 fm with physical volume approximately (2 fm)3 for RBC 2-flavor ensembles and (3 fm)3 for RBC/UKQCD 2+1-flavor dynamical ones. The lightest sea quark mass is about 1/2 the strange quark mass for the former ensembles and 1/4 for the latter ones. Our calculations include: isovector vector- and axial-charge form factors and the first few moments of the polarized and unpolarized structure functions of the nucleon. Nonperturbative renormalization in RI/MOM scheme is applied.
Nucleon Structure from Dynamical Lattice QCD
Lin, H.-W.
2007-06-13
We present lattice QCD numerical calculations of hadronic structure functions and form factors from full-QCD lattices, with a chirally symmetric fermion action, domain-wall fermions, for the sea and valence quarks. The lattice spacing is about 0.12 fm with physical volume approximately (2 fm)3 for RBC 2-flavor ensembles and (3 fm)3 for RBC/UKQCD 2+1-flavor dynamical ones. The lightest sea quark mass is about 1/2 the strange quark mass for the former ensembles and 1/4 for the latter ones. Our calculations include: isovector vector- and axial-charge form factors and the first few moments of the polarized and unpolarized structure functions of the nucleon. Nonperturbative renormalization in RI/MOM scheme is applied.
Kronfeld, A.S.; Allison, I.F.; Aubin, C.; Bernard, C.; Davies, C.T.H.; DeTar, C.; Di Pierro, M.; Freeland, E.D.; Gottlieb, Steven; Gray, A.; Gregor, E.; Heller, U.M.; Hetrick, J.E.; El-Khadra, Aida X.; Levkova, L.; Mackenzie, P.B.; Maresca, F.; Menscher, D.; Nobes, M.; Okamoto, M.; Oktay, M.B.; /Fermilab /Glasgow U. /Columbia U. /Washington U., St. Louis /Utah U. /DePaul U. /Art Inst. of Chicago /Indiana U. /Ohio State U. /Arizona U. /APS, New York /U. Pacific, Stockton /Illinois U., Urbana /Cornell U., LEPP /Simon Fraser U. /UC, Santa Barbara
2005-09-01
In the past year, we calculated with lattice QCD three quantities that were unknown or poorly known. They are the q{sup 2} dependence of the form factor in semileptonic D {yields} K/{nu} decay, the decay constant of the D meson, and the mass of the B{sub c} meson. In this talk, we summarize these calculations, with emphasis on their (subsequent) confirmation by experiments.
Clinical implementation of full Monte Carlo dose calculation in proton beam therapy.
Paganetti, Harald; Jiang, Hongyu; Parodi, Katia; Slopsema, Roelf; Engelsman, Martijn
2008-09-01
The goal of this work was to facilitate the clinical use of Monte Carlo proton dose calculation to support routine treatment planning and delivery. The Monte Carlo code Geant4 was used to simulate the treatment head setup, including a time-dependent simulation of modulator wheels (for broad beam modulation) and magnetic field settings (for beam scanning). Any patient-field-specific setup can be modeled according to the treatment control system of the facility. The code was benchmarked against phantom measurements. Using a simulation of the ionization chamber reading in the treatment head allows the Monte Carlo dose to be specified in absolute units (Gy per ionization chamber reading). Next, the capability of reading CT data information was implemented into the Monte Carlo code to model patient anatomy. To allow time-efficient dose calculation, the standard Geant4 tracking algorithm was modified. Finally, a software link of the Monte Carlo dose engine to the patient database and the commercial planning system was established to allow data exchange, thus completing the implementation of the proton Monte Carlo dose calculation engine ('DoC++'). Monte Carlo re-calculated plans are a valuable tool to revisit decisions in the planning process. Identification of clinically significant differences between Monte Carlo and pencil-beam-based dose calculations may also drive improvements of current pencil-beam methods. As an example, four patients (29 fields in total) with tumors in the head and neck regions were analyzed. Differences between the pencil-beam algorithm and Monte Carlo were identified in particular near the end of range, both due to dose degradation and overall differences in range prediction due to bony anatomy in the beam path. Further, the Monte Carlo reports dose-to-tissue as compared to dose-to-water by the planning system. Our implementation is tailored to a specific Monte Carlo code and the treatment planning system XiO (Computerized Medical Systems Inc
Lattice QCD and Nuclear Physics
Konstantinos Orginos
2007-03-01
A steady stream of developments in Lattice QCD have made it possible today to begin to address the question of how nuclear physics emerges from the underlying theory of strong interactions. Central role in this understanding play both the effective field theory description of nuclear forces and the ability to perform accurate non-perturbative calculations in lo w energy QCD. Here I present some recent results that attempt to extract important low energy constants of the effective field theory of nuclear forces from lattice QCD.
Bazavov, A; Bernard, C; Bouchard, C M; Detar, C; Du, Daping; El-Khadra, A X; Foley, J; Freeland, E D; Gámiz, E; Gottlieb, Steven; Heller, U M; Kim, Jongjeong; Kronfeld, A S; Laiho, J; Levkova, L; Mackenzie, P B; Neil, E T; Oktay, M B; Qiu, Si-Wei; Simone, J N; Sugar, R; Toussaint, D; Van de Water, R S; Zhou, Ran
2014-03-21
We calculate the kaon semileptonic form factor f+(0) from lattice QCD, working, for the first time, at the physical light-quark masses. We use gauge configurations generated by the MILC Collaboration with Nf = 2 + 1 + 1 flavors of sea quarks, which incorporate the effects of dynamical charm quarks as well as those of up, down, and strange. We employ data at three lattice spacings to extrapolate to the continuum limit. Our result, f+(0) = 0.9704(32), where the error is the total statistical plus systematic uncertainty added in quadrature, is the most precise determination to date. Combining our result with the latest experimental measurements of K semileptonic decays, one obtains the Cabibbo-Kobayashi-Maskawa matrix element |V(us)| = 0.22290(74)(52), where the first error is from f+(0) and the second one is from experiment. In the first-row test of Cabibbo-Kobayashi-Maskawa unitarity, the error stemming from |V(us)| is now comparable to that from |V(ud)|. PMID:24702353
Deur, Alexandre; Brodsky, Stanley J.; de Téramond, Guy F.
2016-05-09
Here, we review present knowledge onmore » $$\\alpha_{s}$$, the Quantum Chromodynamics (QCD) running coupling. The dependence of $$\\alpha_s(Q^2)$$ on momentum transfer $Q$ encodes the underlying dynamics of hadron physics --from color confinement in the infrared domain to asymptotic freedom at short distances. We will survey our present theoretical and empirical knowledge of $$\\alpha_s(Q^2)$$, including constraints at high $Q^2$ predicted by perturbative QCD, and constraints at small $Q^2$ based on models of nonperturbative dynamics. In the first, introductory, part of this review, we explain the phenomenological meaning of the coupling, the reason for its running, and the challenges facing a complete understanding of its analytic behavior in the infrared domain. In the second, more technical, part of the review, we discuss $$\\alpha_s(Q^2)$$ in the high momentum transfer domain of QCD. We review how $$\\alpha_s$$ is defined, including its renormalization scheme dependence, the definition of its renormalization scale, the utility of effective charges, as well as `` Commensurate Scale Relations" which connect the various definitions of the QCD coupling without renormalization scale ambiguity. We also report recent important experimental measurements and advanced theoretical analyses which have led to precise QCD predictions at high energy. As an example of an important optimization procedure, we discuss the ``Principle of Maximum Conformality" which enhances QCD's predictive power by removing the dependence of the predictions for physical observables on the choice of the gauge and renormalization scheme. In last part of the review, we discuss $$\\alpha_s(Q^2)$$ in the low momentum transfer domain, where there has been no consensus on how to define $$\\alpha_s(Q^2)$$ or its analytic behavior. We will discuss the various approaches used for low energy calculations. Among them, we will discuss the light-front holographic approach to QCD in the strongly coupled
Calculation of K →π π decay amplitudes with improved Wilson fermion action in lattice QCD
NASA Astrophysics Data System (ADS)
Ishizuka, N.; Ishikawa, K.-I.; Ukawa, A.; Yoshié, T.
2015-10-01
We present our result for the K →π π decay amplitudes for both the Δ I =1 /2 and 3 /2 processes with the improved Wilson fermion action. Expanding on the earlier works by Bernard et al. and by Donini et al., we show that mixings with four-fermion operators with wrong chirality are absent even for the Wilson fermion action for the parity odd process in both channels due to CPS symmetry. Therefore, after subtraction of an effect from the lower dimensional operator, a calculation of the decay amplitudes is possible without complications from operators with wrong chirality, as for the case with chirally symmetric lattice actions. As a first step to verify the possibility of calculations with the Wilson fermion action, we consider the decay amplitudes at an unphysical quark mass mK˜2 mπ . Our calculations are carried out with Nf=2 +1 gauge configurations generated with the Iwasaki gauge action and nonperturbatively O (a )-improved Wilson fermion action at a =0.091 fm , mπ=280 MeV , and mK=580 MeV on a 323×64 (L a =2.9 fm ) lattice. For the quark loops in the penguin and disconnected contributions in the I =0 channel, the combined hopping parameter expansion and truncated solver method work very well for variance reduction. We obtain, for the first time with a Wilson-type fermion action, that Re A0=60 (36 )×1 0-8 GeV and Im A0=-67 (56 )×1 0-12 GeV for a matching scale q*=1 /a . The dependence on the matching scale q* for these values is weak.
Cool QCD: Hadronic Physics and QCD in Nuclei
NASA Astrophysics Data System (ADS)
Cates, Gordon
2015-10-01
QCD is the only strongly-coupled theory given to us by Nature, and it gives rise to a host of striking phenomena. Two examples in hadronic physics include the dynamic generation of mass and the confinement of quarks. Indeed, the vast majority of the mass of visible matter is due to the kinetic and potential energy of the massless gluons and the essentially massless quarks. QCD also gives rise to the force that binds protons and neutrons into nuclei, including subtle effects that have historically been difficult to understand. Describing these phenomena in terms of QCD has represented a daunting task, but remarkable progress has been achieved in both theory and experiment. Both CEBAF at Jefferson Lab and RHIC at Brookhaven National Lab have provided unprecedented experimental tools for investigating QCD, and upgrades at both facilities promise even greater opportunities in the future. Also important are programs at FermiLab as well as the LHC at CERN. Looking further ahead, an electron ion collider (EIC) has the potential to answer whole new sets of questions regarding the role of gluons in nuclear matter, an issue that lies at the heart of the generation of mass. On the theoretical side, rapid progress in supercomputers is enabling stunning progress in Lattice QCD calculations, and approximate forms of QCD are also providing deep new physical insight. In this talk I will describe both recent advances in Cool QCD as well as the exciting scientific opportunities that exist for the future.
A numerical model for calculating vibration from a railway tunnel embedded in a full-space
NASA Astrophysics Data System (ADS)
Hussein, M. F. M.; Hunt, H. E. M.
2007-08-01
Vibration generated by underground railways transmits to nearby buildings causing annoyance to inhabitants and malfunctioning to sensitive equipment. Vibration can be isolated through countermeasures by reducing the stiffness of railpads, using floating-slab tracks and/or supporting buildings on springs. Modelling of vibration from underground railways has recently gained more importance on account of the need to evaluate accurately the performance of vibration countermeasures before these are implemented. This paper develops an existing model, reported by Forrest and Hunt, for calculating vibration from underground railways. The model, known as the Pipe-in-Pipe model, has been developed in this paper to account for anti-symmetrical inputs and therefore to model tangential forces at the tunnel wall. Moreover, three different arrangements of supports are considered for floating-slab tracks, one which can be used to model directly-fixed slabs. The paper also investigates the wave-guided solution of the track, the tunnel, the surrounding soil and the coupled system. It is shown that the dynamics of the track have significant effect on the results calculated in the wavenumber-frequency domain and therefore an important role on controlling vibration from underground railways.
Reconciling the light-cone and nonrelativistic QCD approaches to calculating e+e-→J/ψ+ηc
NASA Astrophysics Data System (ADS)
Bodwin, Geoffrey T.; Kang, Daekyoung; Lee, Jungil
2006-12-01
It has been suggested in Ref. [A. E. Bondar and V. L. Chernyak, Phys. Lett. B 612, 215 (2005)PYLBAJ0370-269310.1016/j.physletb.2005.03.021] that the disagreement between theoretical calculations and experimental observations for the rate for the process e+e-→J/ψ+ηc at the B factories might be resolved by using the light-cone method to take into account the relative momentum of the heavy-quark and antiquark in the quarkonia. The light-cone result for the production cross section in Ref. [A. E. Bondar and V. L. Chernyak, Phys. Lett. B 612, 215 (2005)PYLBAJ0370-269310.1016/j.physletb.2005.03.021] is almost an order-of-magnitude larger than existing NRQCD factorization results. We investigate this apparent theoretical discrepancy. We compute light-cone distribution functions by making use of quarkonium wave functions from the Cornell potential model. Our light-cone distribution functions are similar in shape to those of Ref. [A. E. Bondar and V. L. Chernyak, Phys. Lett. B 612, 215 (2005)PYLBAJ0370-269310.1016/j.physletb.2005.03.021] and yield a similar cross section. However, when we subtract parts of the light-cone distribution functions that correspond to corrections of relative-order αs in the NRQCD approach, we find that the cross section decreases by about a factor of three. When we set certain renormalization factors Zi in the light-cone calculation equal to unity, we find a further reduction in the cross section of about a factor of two. The resulting light-cone cross section is similar in magnitude to the NRQCD factorization cross sections and shows only a modest enhancement over the light-cone cross section in which the relative momentum of the heavy-quark and antiquark is neglected.
Hess, Peter O.
2006-09-25
A review is presented on the contributions of Mexican Scientists to QCD phenomenology. These contributions range from Constituent Quark model's (CQM) with a fixed number of quarks (antiquarks) to those where the number of quarks is not conserved. Also glueball spectra were treated with phenomenological models. Several other approaches are mentioned.
Neoclassical orbit calculations with a full-f code for tokamak edge plasmas
NASA Astrophysics Data System (ADS)
Rognlien, T. D.; Cohen, R. H.; Dorr, M.; Hittinger, J.; Xu, X. Q.; Collela, P.; Martin, D.
2008-11-01
Ion distribution function modifications are considered for the case of neoclassical orbit widths comparable to plasma radial-gradient scale-lengths. Implementation of proper boundary conditions at divertor plates in the continuum TEMPEST code, including the effect of drifts in determining the direction of total flow, enables such calculations in single-null divertor geometry, with and without an electrostatic potential. The resultant poloidal asymmetries in densities, temperatures, and flows are discussed. For long-time simulations, a slow numerical instability develops, even in simplified (circular) geometry with no endloss, which aids identification of the mixed treatment of parallel and radial convection terms as the cause. The new Edge Simulation Laboratory code, expected to be operational, has algorithmic refinements that should address the instability. We will present any available results from the new code on this problem as well as geodesic acoustic mode tests.
Full dimension Rb2He ground triplet potential energy surface and quantum scattering calculations
NASA Astrophysics Data System (ADS)
Guillon, Grégoire; Viel, Alexandra; Launay, Jean-Michel
2012-05-01
We have developed a three-dimensional potential energy surface for the lowest triplet state of the Rb2He complex. A global analytic fit is provided as in the supplementary material [see supplementary material at http://dx.doi.org/10.1063/1.4709433E-JCPSA6-136-034218 for the corresponding Fortran code]. This surface is used to perform quantum scattering calculations of 4He and 3He colliding with 87Rb2 in the partial wave J = 0 at low and ultralow energies. For the heavier helium isotope, the computed vibrational relaxation probabilities show a broad and strong shape resonance for a collisional energy of 0.15 K and a narrow Feshbach resonance at about 17 K for all initial Rb2 vibrational states studied. The broad resonance corresponds to an efficient relaxation mechanism that does not occur when 3He is the colliding partner. The Feshbach resonance observed at higher collisional energy is robust with respect to the isotopic substitution. However, its effect on the vibrational relaxation mechanism is faint for both isotopes.
Full dimension Rb2He ground triplet potential energy surface and quantum scattering calculations.
Guillon, Grégoire; Viel, Alexandra; Launay, Jean-Michel
2012-05-01
We have developed a three-dimensional potential energy surface for the lowest triplet state of the Rb(2)He complex. A global analytic fit is provided as in the supplementary material [see supplementary material at http://dx.doi.org/10.1063/1.4709433 for the corresponding Fortran code]. This surface is used to perform quantum scattering calculations of (4)He and (3)He colliding with (87)Rb(2) in the partial wave J = 0 at low and ultralow energies. For the heavier helium isotope, the computed vibrational relaxation probabilities show a broad and strong shape resonance for a collisional energy of 0.15 K and a narrow Feshbach resonance at about 17 K for all initial Rb(2) vibrational states studied. The broad resonance corresponds to an efficient relaxation mechanism that does not occur when (3)He is the colliding partner. The Feshbach resonance observed at higher collisional energy is robust with respect to the isotopic substitution. However, its effect on the vibrational relaxation mechanism is faint for both isotopes. PMID:22583230
Bret, A.
2014-02-15
The filamentation (Weibel) instability plays a key role in the formation of collisionless shocks which are thought to produce Gamma-Ray-Bursts and High-Energy-Cosmic-Rays in astrophysical environments. While it has been known for long that a flow-aligned magnetic field can completely quench the instability, it was recently proved in 2D that in the cold regime, such cancelation is possible if and only if the field is perfectly aligned. Here, this result is finally extended to a 3D geometry. Calculations are conducted for symmetric and asymmetric counter-streaming relativistic plasma shells. 2D results are retrieved in 3D: the instability can never be completely canceled for an oblique magnetic field. In addition, the maximum growth-rate is always larger for wave vectors lying in the plan defined by the flow and the oblique field. On the one hand, this bears consequences on the orientation of the generated filaments. On the other hand, it certifies 2D simulations of the problem can be performed without missing the most unstable filamentation modes.
Giannetti, P. )
1991-05-01
Recent analysis of jet data taken at the Fermilab Tevatron Collider at {radical}S = 1.8 Tev are presented. Inclusive jet, dijet, trijet and direct photon measurements are compared to QCD parton level calculations, at orders {alpha}{sub s}{sup 3} or {alpha}{sub s}{sup 2}. The large total transverse energy events are well described by the Herwig shower Montecarlo. 19 refs., 20 figs., 1 tab.
Roberts, C.D.
1994-09-01
The Dyson-Schwinger equations (DSEs) are a tower of coupled integral equations that relate the Green functions of QCD to one another. Solving these equations provides the solution of QCD. This tower of equations includes the equation for the quark self-energy, which is the analogue of the gap equation in superconductivity, and the Bethe-Salpeter equation, the solution of which is the quark-antiquark bound state amplitude in QCD. The application of this approach to solving Abelian and non-Abelian gauge theories is reviewed. The nonperturbative DSE approach is being developed as both: (1) a computationally less intensive alternative and; (2) a complement to numerical simulations of the lattice action of QCD. In recent years, significant progress has been made with the DSE approach so that it is now possible to make sensible and direct comparisons between quantities calculated using this approach and the results of numerical simulations of Abelian gauge theories. Herein the application of the DSE approach to the calculation of pion observables is described: the {pi}-{pi} scattering lengths (a{sub 0}{sup 0}, a{sub 0}{sup 2}, A{sub 1}{sup 1}, a{sub 2}{sup 2}) and associated partial wave amplitudes; the {pi}{sup 0} {yields} {gamma}{gamma} decay width; and the charged pion form factor, F{sub {pi}}(q{sup 2}). Since this approach provides a straightforward, microscopic description of dynamical chiral symmetry breaking (D{sub X}SB) and confinement, the calculation of pion observables is a simple and elegant illustrative example of its power and efficacy. The relevant DSEs are discussed in the calculation of pion observables and concluding remarks are presented.
Mixed augmented variational formulation (MAVF) for lower hybrid full-wave calculations
NASA Astrophysics Data System (ADS)
Peysson, Y.; Roche, J. R.; Bertrand, P.; Chatenet, J.-H.; Kirsch, C.; Mokrani, A.; Labrunie, S.
2009-11-01
In the continuation of the works led in cylindrical geometry [2], a full toroidal description for an arbitrary poloidal cross-section of the plasma has been developed. For simulation purpose a mixed augmented variational formulation (MAVF), which is particularly well suited for solving Maxwell equations, is considered [4]. The discretization of the MAVF is carried out using Taylor-Hood P2-iso-P1 finite elements. This formulation provide a natural implementation for parallel processing, a particularly important aspect when simulations for plasmas of large size must be considered. Details on the specific application of the MAVF to the LH problem are presented, as well as the structure of the corresponding matrices. A first application to a realistic small tokamak configuration is considered.
Numerical calculation of steady inviscid full potential compressible flow about wind turbine blades
NASA Technical Reports Server (NTRS)
Dulikravich, D. S.
1980-01-01
The air flow through a propeller-type wind turbine rotor is characterized by three-dimensional rotating cascade effects about the inner portions of the rotor blades and compressibility effects about the tip regions of the blades. In the case of large rotor diameter and/or increased rotor angular speed, the existence of small supersonic zones terminated by weak shocks is possible. An exact nonlinear mathematical model (called a steady Full Potential Equation - FPE) that accounts for the above phenomena has been rederived. An artificially time dependent version of FPE was iteratively solved by a finite volume technique involving an artificial viscosity and a three-level consecutive mesh refinement. The exact boundary conditions were applied by generating a boundary conforming periodic computation mesh.
Verley, R.; Ness, O.B.
1995-12-31
This paper presents the results of full scale bending tests on 16 in. and 20 in. diameter, concrete coated pipes with polyethene and asphalt corrosion coatings. Constant moment, four-point bending was applied to a pipe string consisting of one pipe joint welded between two half-length joints. The strain concentration factor (SCF) at the field joints (FJ), expressing the ratio between the strain in the FJ and the average strain for the pipe joint, was investigated and compared to predictions using an analytical model presented in an accompanying paper (Ness and Verley, 1995). Material tests on the pipe steel, the corrosion coating and the concrete were conducted. The analytical model is found to give a good prediction of the strain distribution along the pipe joint, for both the steel and the concrete, and therefore also of the SCF. The sliding of the concrete over the steel is also predicted reasonably well.
Nawa, Kanabu; Suganuma, Hideo; Kojo, Toru
2007-04-15
We study baryons in holographic QCD with D4/D8/D8 multi-D-brane system. In holographic QCD, the baryon appears as a topologically nontrivial chiral soliton in a four-dimensional effective theory of mesons. We call this topological soliton brane-induced Skyrmion. Some review of D4/D8/D8 holographic QCD is presented from the viewpoints of recent hadron physics and QCD phenomenologies. A four-dimensional effective theory with pions and {rho} mesons is uniquely derived from the non-Abelian Dirac-Born-Infeld (DBI) action of D8 brane with D4 supergravity background at the leading order of large N{sub c}, without small amplitude expansion of meson fields to discuss chiral solitons. For the hedgehog configuration of pion and {rho}-meson fields, we derive the energy functional and the Euler-Lagrange equation of brane-induced Skyrmion from the meson effective action induced by holographic QCD. Performing the numerical calculation, we obtain the soliton solution and figure out the pion profile F(r) and the {rho}-meson profile G-tilde(r) of the brane-induced Skyrmion with its total energy, energy density distribution, and root-mean-square radius. These results are compared with the experimental quantities of baryons and also with the profiles of standard Skyrmion without {rho} mesons. We analyze interaction terms of pions and {rho} mesons in brane-induced Skyrmion, and find a significant {rho}-meson component appearing in the core region of a baryon.
FOREWORD: Extreme QCD 2012 (xQCD)
NASA Astrophysics Data System (ADS)
Alexandru, Andrei; Bazavov, Alexei; Liu, Keh-Fei
2013-04-01
The Extreme QCD 2012 conference, held at the George Washington University in August 2012, celebrated the 10th event in the series. It has been held annually since 2003 at different locations: San Carlos (2011), Bad Honnef (2010), Seoul (2009), Raleigh (2008), Rome (2007), Brookhaven (2006), Swansea (2005), Argonne (2004), and Nara (2003). As usual, it was a very productive and inspiring meeting that brought together experts in the field of finite-temperature QCD, both theoretical and experimental. On the experimental side, we heard about recent results from major experiments, such as PHENIX and STAR at Brookhaven National Laboratory, ALICE and CMS at CERN, and also about the constraints on the QCD phase diagram coming from astronomical observations of one of the largest laboratories one can imagine, neutron stars. The theoretical contributions covered a wide range of topics, including QCD thermodynamics at zero and finite chemical potential, new ideas to overcome the sign problem in the latter case, fluctuations of conserved charges and how they allow one to connect calculations in lattice QCD with experimentally measured quantities, finite-temperature behavior of theories with many flavors of fermions, properties and the fate of heavy quarkonium states in the quark-gluon plasma, and many others. The participants took the time to write up and revise their contributions and submit them for publication in these proceedings. Thanks to their efforts, we have now a good record of the ideas presented and discussed during the workshop. We hope that this will serve both as a reminder and as a reference for the participants and for other researchers interested in the physics of nuclear matter at high temperatures and density. To preserve the atmosphere of the event the contributions are ordered in the same way as the talks at the conference. We are honored to have helped organize the 10th meeting in this series, a milestone that reflects the lasting interest in this
Brodsky, Stanley J.; Deshpande, Abhay L.; Gao, Haiyan; McKeown, Robert D.; Meyer, Curtis A.; Meziani, Zein-Eddine; Milner, Richard G.; Qiu, Jianwei; Richards, David G.; Roberts, Craig D.
2015-02-26
This White Paper presents the recommendations and scientific conclusions from the Town Meeting on QCD and Hadronic Physics that took place in the period 13-15 September 2014 at Temple University as part of the NSAC 2014 Long Range Planning process. The meeting was held in coordination with the Town Meeting on Phases of QCD and included a full day of joint plenary sessions of the two meetings. The goals of the meeting were to report and highlight progress in hadron physics in the seven years since the 2007 Long Range Plan (LRP07), and present a vision for the future by identifying the key questions and plausible paths to solutions which should define the next decade. The introductory summary details the recommendations and their supporting rationales, as determined at the Town Meeting on QCD and Hadron Physics, and the endorsements that were voted upon. The larger document is organized as follows. Section 2 highlights major progress since the 2007 LRP. It is followed, in Section 3, by a brief overview of the physics program planned for the immediate future. Finally, Section 4 provides an overview of the physics motivations and goals associated with the next QCD frontier: the Electron-Ion-Collider.
None
2011-10-06
Modern QCD - Lecture 4 We will consider some processes of interest at the LHC and will discuss the main elements of their cross-section calculations. We will also summarize the current status of higher order calculations.
Full-dimensional vibrational calculations for H5O2+ using an ab initio potential energy surface
NASA Astrophysics Data System (ADS)
McCoy, Anne B.; Huang, Xinchuan; Carter, Stuart; Landeweer, Marc Y.; Bowman, Joel M.
2005-02-01
We report quantum diffusion Monte Carlo (DMC) and variational calculations in full dimensionality for selected vibrational states of H5O2+ using a new ab initio potential energy surface [X. Huang, B. Braams, and J. M. Bowman, J. Chem. Phys. 122, 044308 (2005)]. The energy and properties of the zero-point state are focused on in the rigorous DMC calculations. OH-stretch fundamentals are also calculated using "fixed-node" DMC calculations and variationally using two versions of the code MULTIMODE. These results are compared with infrared multiphoton dissociation measurements of Yeh et al. [L. I. Yeh, M. Okumura, J. D. Myers, J. M. Price, and Y. T. Lee, J. Chem. Phys. 91, 7319 (1989)]. Some preliminary results for the energies of several modes of the shared hydrogen are also reported.
Full-CI calculations of alternant cyclic polyenes (CH) N, N=2, 4, 6, ƒ 18, in the PPP approximation
NASA Astrophysics Data System (ADS)
Bendazzoli, G. L.; Evangelisti, S.
1991-10-01
Full-CI calculations are performed on the PPP model of alternant cyclic polyenes using a recently developed direct full-CI algorithm based on an expansion in Slater determinants. By exploiting the cyclic symmetry and a partial spin adaptation, we have been able to compute a full-CI wavefunction and energy up to N = 18 on the CRAY YMP 8/432 of the Italian University Computing Centre CINECA. The wavefunctions and energies for the series N = 4ν and N = 4ν + 2 are discussed and their behaviour as functions of N is analyzed.
NASA Astrophysics Data System (ADS)
Dudek, Jozef J.
2016-03-01
I describe how hadron-hadron scattering amplitudes are related to the eigenstates of QCD in a finite cubic volume. The discrete spectrum of such eigenstates can be determined from correlation functions computed using lattice QCD, and the corresponding scattering amplitudes extracted. I review results from the Hadron Spectrum Collaboration who have used these finite volume methods to study ππ elastic scattering, including the ρ resonance, as well as coupled-channel πK, ηK scattering. The very recent extension to the case where an external current acts is also presented, considering the reaction πγ* → ππ, from which the unstable ρ → πγ transition form factor is extracted. Ongoing calculations are advertised and the outlook for finite volume approaches is presented.
Gupta, R.
1998-12-31
The goal of the lectures on lattice QCD (LQCD) is to provide an overview of both the technical issues and the progress made so far in obtaining phenomenologically useful numbers. The lectures consist of three parts. The author`s charter is to provide an introduction to LQCD and outline the scope of LQCD calculations. In the second set of lectures, Guido Martinelli will discuss the progress they have made so far in obtaining results, and their impact on Standard Model phenomenology. Finally, Martin Luescher will discuss the topical subjects of chiral symmetry, improved formulation of lattice QCD, and the impact these improvements will have on the quality of results expected from the next generation of simulations.
QCD measurements at the Tevatron
Bandurin, Dmitry; /Florida State U.
2011-12-01
Selected quantum chromodynamics (QCD) measurements performed at the Fermilab Run II Tevatron p{bar p} collider running at {radical}s = 1.96 TeV by CDF and D0 Collaborations are presented. The inclusive jet, dijet production and three-jet cross section measurements are used to test perturbative QCD calculations, constrain parton distribution function (PDF) determinations, and extract a precise value of the strong coupling constant, {alpha}{sub s}(m{sub Z}) = 0.1161{sub -0.0048}{sup +0.0041}. Inclusive photon production cross-section measurements reveal an inability of next-to-leading-order (NLO) perturbative QCD (pQCD) calculations to describe low-energy photons arising directly in the hard scatter. The diphoton production cross-sections check the validity of the NLO pQCD predictions, soft-gluon resummation methods implemented in theoretical calculations, and contributions from the parton-to-photon fragmentation diagrams. Events with W/Z+jets productions are used to measure many kinematic distributions allowing extensive tests and tunes of predictions from pQCD NLO and Monte-Carlo (MC) event generators. The charged-particle transverse momenta (p{sub T}) and multiplicity distributions in the inclusive minimum bias events are used to tune non-perturbative QCD models, including those describing the multiple parton interactions (MPI). Events with inclusive production of {gamma} and 2 or 3 jets are used to study increasingly important MPI phenomenon at high p{sub T}, measure an effective interaction cross section, {sigma}{sub eff} = 16.4 {+-} 2.3 mb, and limit existing MPI models.
NASA Technical Reports Server (NTRS)
Meitner, P. L.
1978-01-01
A computer program that calculates the coolant flow and the metal temperatures of a full-coverage-film-cooled vane or blade was developed. The analysis was based on compressible, one-dimensional fluid flow and on one-dimensional heat transfer and treats the vane or blade shell as a porous wall. The calculated temperatures are average values for the shell outer-surface area associated with each film-cooling hole row. A thermal-barrier coating may be specified on the shell outer surface, and centrifugal effects can be included for blade calculations. The program is written in FORTRAN 4 and is operational on a UNIVAC 1100/42 computer. The method of analysis, the program input, the program output, and two sample problems are provided.
QCD at nonzero chemical potential: Recent progress on the lattice
NASA Astrophysics Data System (ADS)
Aarts, Gert; Attanasio, Felipe; Jäger, Benjamin; Seiler, Erhard; Sexty, Dénes; Stamatescu, Ion-Olimpiu
2016-01-01
We summarise recent progress in simulating QCD at nonzero baryon density using complex Langevin dynamics. After a brief outline of the main idea, we discuss gauge cooling as a means to control the evolution. Subsequently we present a status report for heavy dense QCD and its phase structure, full QCD with staggered quarks, and full QCD with Wilson quarks, both directly and using the hopping parameter expansion to all orders.
Nikitin, A. V.; Rey, M.; Tyuterev, Vl. G.
2015-03-07
A simultaneous use of the full molecular symmetry and of an exact kinetic energy operator (KEO) is of key importance for accurate predictions of vibrational levels at a high energy range from a potential energy surface (PES). An efficient method that permits a fast convergence of variational calculations would allow iterative optimization of the PES parameters using experimental data. In this work, we propose such a method applied to tetrahedral AB{sub 4} molecules for which a use of high symmetry is crucial for vibrational calculations. A symmetry-adapted contracted angular basis set for six redundant angles is introduced. Simple formulas using this basis set for explicit calculation of the angular matrix elements of KEO and PES are reported. The symmetric form (six redundant angles) of vibrational KEO without the sin(q){sup −2} type singularity is derived. The efficient recursive algorithm based on the tensorial formalism is used for the calculation of vibrational matrix elements. A good basis set convergence for the calculations of vibrational levels of the CH{sub 4} molecule is demonstrated.
NASA Astrophysics Data System (ADS)
Nikitin, A. V.; Rey, M.; Tyuterev, Vl. G.
2015-03-01
A simultaneous use of the full molecular symmetry and of an exact kinetic energy operator (KEO) is of key importance for accurate predictions of vibrational levels at a high energy range from a potential energy surface (PES). An efficient method that permits a fast convergence of variational calculations would allow iterative optimization of the PES parameters using experimental data. In this work, we propose such a method applied to tetrahedral AB4 molecules for which a use of high symmetry is crucial for vibrational calculations. A symmetry-adapted contracted angular basis set for six redundant angles is introduced. Simple formulas using this basis set for explicit calculation of the angular matrix elements of KEO and PES are reported. The symmetric form (six redundant angles) of vibrational KEO without the sin(q)-2 type singularity is derived. The efficient recursive algorithm based on the tensorial formalism is used for the calculation of vibrational matrix elements. A good basis set convergence for the calculations of vibrational levels of the CH4 molecule is demonstrated.
Lattice and Phase Diagram in QCD
Lombardo, Maria Paola
2008-10-13
Model calculations have produced a number of very interesting expectations for the QCD Phase Diagram, and the task of a lattice calculations is to put these studies on a quantitative grounds. I will give an overview of the current status of the lattice analysis of the QCD phase diagram, from the quantitative results of mature calculations at zero and small baryochemical potential, to the exploratory studies of the colder, denser phase.
Yamamoto, Arata
2016-07-29
We propose the lattice QCD calculation of the Berry phase, which is defined by the ground state of a single fermion. We perform the ground-state projection of a single-fermion propagator, construct the Berry link variable on a momentum-space lattice, and calculate the Berry phase. As the first application, the first Chern number of the (2+1)-dimensional Wilson fermion is calculated by the Monte Carlo simulation. PMID:27517766
NASA Astrophysics Data System (ADS)
Diaz, Susana B.; Frederick, John E.; Lucas, Timothy; Booth, C. Rocky; Smolskaia, Irina
A spectroradiometer located at Ushuaia, Tierra del Fuego, Argentina obtained a high quality data set on the solar ultraviolet (UV) spectral irradiance for a full annual cycle during 1992. The unique aspect of the ground-based irradiances is their ability to characterize the effects of cloudiness. Measured irradiances at 340 nm can be used to characterize the attenuation provided by cloudy skies. When irradiances at shorter wavelengths, 302.5-320.0 nm, are adjusted for this attenuation, they show good agreement with radiative transfer calculations which assume clear skies and utilize data from the Total Ozone Mapping Spectrometer (TOMS) as inputs. The differences which exist are qualitatively consistent with the neglect of spherical geometry in the calculated irradiances.
NASA Astrophysics Data System (ADS)
Kokubun, Yasuo; Watanabe, Tatsuhiko; Kawata, Ryo; Morita, Kohei
2016-08-01
We propose and demonstrate a novel mode analysis method that can provide a full set of amplitudes, phases, and polarization states of guided modes including degenerate modes in few-mode fibers. The method is based on the calculation of amplitude components and phase differences from the intensity profiles that passed through a polarizer with angles of 0, 45, and 90°. The accuracies of calculation formulas are shown by simulation to be less than 10‑10 for amplitudes and 10‑3 rad for phase differences. The method was applied to a graded-index three-mode fiber, and the off-axis mode excitation ratio characteristics were compared with theoretical ones.
NASA Astrophysics Data System (ADS)
Okaya, D. A.; Van Avendonk, H. J.
2013-12-01
Recent anisotropy studies at scales ranging from crust to full mantle have recognized the importance of 3D anisotropy geometry and heterogeneity as well as variability in anisotropic symmetry and orientation (tilt) of the Earth. The strong relationship between seismic anisotropy and geodynamic processes highlights the need to construct realistic Earth models that can explain observations of anisotropy in modern seismic data sets. For example, ray paths through a mantle slab window or a mountain belt may show that the crust or mantle exhibits low-order anisotropy due to a history of deformation and the development of tectonic fabrics. Observed traveltimes might not be fit with simple Transverse Isotropy (TI), so realistic calculations require an Earth model that accurately describes the wave speeds of compressional and shear waves. We have developed an anisotropic traveltime solver that allows for full 3D heterogeneity of anisotropy tensors, degrees of symmetry, and arbitrary orientation. This traveltime solver is based on the robust shortest path method (SPM) and a ray-bending algorithm that were previously applied to isotropic media (e.g., Van Avendonk et al., 2001). Instead of using an isotropic description of the seismic wave velocity, we define the full elastic tensor at each location in the model. The directional seismic velocity can subsequently be extracted using solutions of the Christoffel equations. For computational efficiency, we calculate all directional seismic velocities at each model node before the start of ray tracing. As we calculate a new ray segment, this information is quickly retrieved. We use these directional velocity maps to separately describe the propagation of compressional (P) and shear (S) body waves in anisotropic media and to subsequently calculate their traveltimes. Patterns within the velocity maps represent tensor symmetries and tilts, allowing for the construction of discretized large-scale 3D LPO flow fields or fabric
HFOLD - A program package for calculating two-body MSSM Higgs decays at full one-loop level.
Frisch, W; Eberl, H; Hluchá, H
2011-10-01
HFOLD (Higgs Full One Loop Decays) is a Fortran program package for calculating all MSSM Higgs two-body decay widths and the corresponding branching ratios at full one-loop level. The package is done in the SUSY Parameter Analysis convention and supports the SUSY Les Houches Accord input and output format. PROGRAM SUMMARY: Program title: HFOLD Catalogue identifier: AEJG_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEJG_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 340 621 No. of bytes in distributed program, including test data, etc.: 1 760 051 Distribution format: tar.gz Programming language: Fortran 77 Computer: Workstation, PC Operating system: Linux RAM: 524 288 000 Bytes Classification: 11.1 External routines: LoopTools 2.2 (http://www.feynarts.de/looptools/), SLHALib 2.2 (http://www.feynarts.de/slha/). The LoopTools code is included in the distribution package. Nature of problem: A future high-energy e+e- linear collider will be the best environment for the precise measurements of masses, cross sections, branching ratios, etc. Experimental accuracies are expected at the per-cent down to the per-mile level. These must be matched from the theoretical side. Therefore higher order calculations are mandatory. Solution method: This program package calculates all MSSM Higgs two-body decay widths and the corresponding branching ratios at full one-loop level. The renormalization is done in the DR scheme following the SUSY Parameter Analysis convention. The program supports the SUSY Les Houches Accord input and output format. Running time: The example provided takes only a few seconds to run. PMID:21969735
Lattice QCD in Background Fields
William Detmold, Brian Tiburzi, Andre Walker-Loud
2009-06-01
Electromagnetic properties of hadrons can be computed by lattice simulations of QCD in background fields. We demonstrate new techniques for the investigation of charged hadron properties in electric fields. Our current calculations employ large electric fields, motivating us to analyze chiral dynamics in strong QED backgrounds, and subsequently uncover surprising non-perturbative effects present at finite volume.
Parton distributions from lattice QCD: an update
Detmold, W; Melnitchouk, W; Thomas, A W
2004-04-01
We review the extraction of parton distributions from their moments calculated in lattice QCD, focusing in particular on their extrapolation to the physical region. As examples, we consider both the unpolarized and polarized isovector parton distributions of the nucleon.
Dru Renner
2012-04-01
Precision computation of hadronic physics with lattice QCD is becoming feasible. The last decade has seen precent-level calculations of many simple properties of mesons, and the last few years have seen calculations of baryon masses, including the nucleon mass, accurate to a few percent. As computational power increases and algorithms advance, the precise calculation of a variety of more demanding hadronic properties will become realistic. With this in mind, I discuss the current lattice QCD calculations of generalized parton distributions with an emphasis on the prospects for well-controlled calculations for these observables as well. I will do this by way of several examples: the pion and nucleon form factors and moments of the nucleon parton and generalized-parton distributions.
Glueball decay in holographic QCD
Hashimoto, Koji; Tan, C.-I; Terashima, Seiji
2008-04-15
Using holographic QCD based on D4-branes and D8-anti-D8-branes, we have computed couplings of glueballs to light mesons. We describe glueball decay by explicitly calculating its decay widths and branching ratios. Interestingly, while glueballs remain less well understood both theoretically and experimentally, our results are found to be consistent with the experimental data for the scalar glueball candidate f{sub 0}(1500). More generally, holographic QCD predicts that decay of any glueball to 4{pi}{sup 0} is suppressed, and that mixing of the lightest glueball with qq mesons is small.
Neutron star structure from QCD
NASA Astrophysics Data System (ADS)
Fraga, Eduardo S.; Kurkela, Aleksi; Vuorinen, Aleksi
2016-03-01
In this review article, we argue that our current understanding of the thermodynamic properties of cold QCD matter, originating from first principles calculations at high and low densities, can be used to efficiently constrain the macroscopic properties of neutron stars. In particular, we demonstrate that combining state-of-the-art results from Chiral Effective Theory and perturbative QCD with the current bounds on neutron star masses, the Equation of State of neutron star matter can be obtained to an accuracy better than 30% at all densities.
QCD inequalities for hadron interactions.
Detmold, William
2015-06-01
We derive generalizations of the Weingarten-Witten QCD mass inequalities for particular multihadron systems. For systems of any number of identical pseudoscalar mesons of maximal isospin, these inequalities prove that near threshold interactions between the constituent mesons must be repulsive and that no bound states can form in these channels. Similar constraints in less symmetric systems are also extracted. These results are compatible with experimental results (where known) and recent lattice QCD calculations, and also lead to a more stringent bound on the nucleon mass than previously derived, m_{N}≥3/2m_{π}. PMID:26196617
QCD corrections to triboson production
NASA Astrophysics Data System (ADS)
Lazopoulos, Achilleas; Melnikov, Kirill; Petriello, Frank
2007-07-01
We present a computation of the next-to-leading order QCD corrections to the production of three Z bosons at the Large Hadron Collider. We calculate these corrections using a completely numerical method that combines sector decomposition to extract infrared singularities with contour deformation of the Feynman parameter integrals to avoid internal loop thresholds. The NLO QCD corrections to pp→ZZZ are approximately 50% and are badly underestimated by the leading order scale dependence. However, the kinematic dependence of the corrections is minimal in phase space regions accessible at leading order.
Precision lattice QCD: challenges and prospects
NASA Astrophysics Data System (ADS)
Hashimoto, Shoji
2013-04-01
With Peta-flops scale computational resources, lattice QCD simulation has recently reached one of its primary goals, i.e. reproducing the low-lying hadron spectrum starting from the QCD Lagrangian. Applications to various other phenomenological quantities, for which no other way of precise theoretical calculation is available, would become the next milestone. In this talk I will provide a brief overview of the field and summarize the remaining problems to be solved before achieving the precision calculations.
Lattice QCD and the Jefferson Laboratory Program
Jozef Dudek, Robert Edwards, David Richards, Konstantinos Orginos
2011-06-01
Lattice gauge theory provides our only means of performing \\textit{ab initio} calculations in the non-perturbative regime. It has thus become an increasing important component of the Jefferson Laboratory physics program. In this paper, we describe the contributions of lattice QCD to our understanding of hadronic and nuclear physics, focusing on the structure of hadrons, the calculation of the spectrum and properties of resonances, and finally on deriving an understanding of the QCD origin of nuclear forces.
Kenneth Wilson and Lattice QCD
NASA Astrophysics Data System (ADS)
Ukawa, Akira
2015-09-01
We discuss the physics and computation of lattice QCD, a space-time lattice formulation of quantum chromodynamics, and Kenneth Wilson's seminal role in its development. We start with the fundamental issue of confinement of quarks in the theory of the strong interactions, and discuss how lattice QCD provides a framework for understanding this phenomenon. A conceptual issue with lattice QCD is a conflict of space-time lattice with chiral symmetry of quarks. We discuss how this problem is resolved. Since lattice QCD is a non-linear quantum dynamical system with infinite degrees of freedom, quantities which are analytically calculable are limited. On the other hand, it provides an ideal case of massively parallel numerical computations. We review the long and distinguished history of parallel-architecture supercomputers designed and built for lattice QCD. We discuss algorithmic developments, in particular the difficulties posed by the fermionic nature of quarks, and their resolution. The triad of efforts toward better understanding of physics, better algorithms, and more powerful supercomputers have produced major breakthroughs in our understanding of the strong interactions. We review the salient results of this effort in understanding the hadron spectrum, the Cabibbo-Kobayashi-Maskawa matrix elements and CP violation, and quark-gluon plasma at high temperatures. We conclude with a brief summary and a future perspective.
Nuclear reactions from lattice QCD
Briceño, Raúl A.; Davoudi, Zohreh; Luu, Thomas C.
2015-01-13
In this study, one of the overarching goals of nuclear physics is to rigorously compute properties of hadronic systems directly from the fundamental theory of strong interactions, Quantum Chromodynamics (QCD). In particular, the hope is to perform reliable calculations of nuclear reactions which will impact our understanding of environments that occur during big bang nucleosynthesis, the evolution of stars and supernovae, and within nuclear reactors and high energy/density facilities. Such calculations, being truly ab initio, would include all two-nucleon and three- nucleon (and higher) interactions in a consistent manner. Currently, lattice QCD provides the only reliable option for performing calculations of some of the low-energy hadronic observables. With the aim of bridging the gap between lattice QCD and nuclear many-body physics, the Institute for Nuclear Theory held a workshop on Nuclear Reactions from Lattice QCD on March 2013. In this review article, we report on the topics discussed in this workshop and the path planned to move forward in the upcoming years.
Axion cosmology, lattice QCD and the dilute instanton gas
NASA Astrophysics Data System (ADS)
Borsanyi, Sz.; Dierigl, M.; Fodor, Z.; Katz, S. D.; Mages, S. W.; Nogradi, D.; Redondo, J.; Ringwald, A.; Szabo, K. K.
2016-01-01
Axions are one of the most attractive dark matter candidates. The evolution of their number density in the early universe can be determined by calculating the topological susceptibility χ (T) of QCD as a function of the temperature. Lattice QCD provides an ab initio technique to carry out such a calculation. A full result needs two ingredients: physical quark masses and a controlled continuum extrapolation from non-vanishing to zero lattice spacings. We determine χ (T) in the quenched framework (infinitely large quark masses) and extrapolate its values to the continuum limit. The results are compared with the prediction of the dilute instanton gas approximation (DIGA). A nice agreement is found for the temperature dependence, whereas the overall normalization of the DIGA result still differs from the non-perturbative continuum extrapolated lattice results by a factor of order ten. We discuss the consequences of our findings for the prediction of the amount of axion dark matter.
Innovations in Lattice QCD Algorithms
Konstantinos Orginos
2006-06-25
Lattice QCD calculations demand a substantial amount of computing power in order to achieve the high precision results needed to better understand the nature of strong interactions, assist experiment to discover new physics, and predict the behavior of a diverse set of physical systems ranging from the proton itself to astrophysical objects such as neutron stars. However, computer power alone is clearly not enough to tackle the calculations we need to be doing today. A steady stream of recent algorithmic developments has made an important impact on the kinds of calculations we can currently perform. In this talk I am reviewing these algorithms and their impact on the nature of lattice QCD calculations performed today.
NASA Astrophysics Data System (ADS)
Orlando, Roberto; De La Pierre, Marco; Zicovich-Wilson, Claudio M.; Erba, Alessandro; Dovesi, Roberto
2014-09-01
Use of symmetry can dramatically reduce the computational cost (running time and memory allocation) of self-consistent-field ab initio calculations for molecular and crystalline systems. Crucial for running time is symmetry exploitation in the evaluation of one- and two-electron integrals, diagonalization of the Fock matrix at selected points in reciprocal space, reconstruction of the density matrix. As regards memory allocation, full square matrices (overlap, Fock, and density) in the Atomic Orbital (AO) basis are avoided and a direct transformation from the packed AO to the symmetry adapted crystalline orbital basis is performed, so that the largest matrix to be handled has the size of the largest sub-block in the latter basis. Quantitative examples, referring to the implementation in the CRYSTAL code, are given for high symmetry families of compounds such as carbon fullerenes and nanotubes.
Orlando, Roberto Erba, Alessandro; Dovesi, Roberto; De La Pierre, Marco; Zicovich-Wilson, Claudio M.
2014-09-14
Use of symmetry can dramatically reduce the computational cost (running time and memory allocation) of self-consistent-field ab initio calculations for molecular and crystalline systems. Crucial for running time is symmetry exploitation in the evaluation of one- and two-electron integrals, diagonalization of the Fock matrix at selected points in reciprocal space, reconstruction of the density matrix. As regards memory allocation, full square matrices (overlap, Fock, and density) in the Atomic Orbital (AO) basis are avoided and a direct transformation from the packed AO to the symmetry adapted crystalline orbital basis is performed, so that the largest matrix to be handled has the size of the largest sub-block in the latter basis. Quantitative examples, referring to the implementation in the CRYSTAL code, are given for high symmetry families of compounds such as carbon fullerenes and nanotubes.
NASA Technical Reports Server (NTRS)
Bauschlicher, Charles W., Jr.; Langhoff, Stephen R.
1987-01-01
Full configuration interaction (CI) calculations on the ground states of N2, NO, and O2 using a DZP Gaussian basis are compared with single-reference SDCI and coupled pair approaches (CPF), as well as with CASSCF multireference CI approaches. The CASSCF/MRCI technique is found to describe multiple bonds as well as single bonds. Although the coupled pair functional approach gave chemical accuracy (1 kcal/mol) for bonds involving hydrogen, larger errors occur in the CPF approach for the multiple bonded systems considered here. CI studies on the 1Sigma(g +) state of N2, including all single, double, triple, and quadruple excitations show that triple excitations are very important for the multiple bond case, and accounts for most of the deficiency in the coupled pair functional methods.
NASA Technical Reports Server (NTRS)
Farrell, C. A.
1982-01-01
A fast, reliable computer code is described for calculating the flow field about a cascade of arbitrary two dimensional airfoils. The method approximates the three dimensional flow in a turbomachinery blade row by correcting for stream tube convergence and radius change in the throughflow direction. A fully conservative solution of the full potential equation is combined with the finite volume technique on a body-fitted periodic mesh, with an artificial density imposed in the transonic region to insure stability and the capture of shock waves. The instructions required to set up and use the code are included. The name of the code is QSONIC. A numerical example is also given to illustrate the output of the program.
NASA Technical Reports Server (NTRS)
Jutte, Christine V.; Ko, William L.; Stephens, Craig A.; Bakalyar, John A.; Richards, W. Lance
2011-01-01
A ground loads test of a full-scale wing (175-ft span) was conducted using a fiber optic strain-sensing system to obtain distributed surface strain data. These data were input into previously developed deformed shape equations to calculate the wing s bending and twist deformation. A photogrammetry system measured actual shape deformation. The wing deflections reached 100 percent of the positive design limit load (equivalent to 3 g) and 97 percent of the negative design limit load (equivalent to -1 g). The calculated wing bending results were in excellent agreement with the actual bending; tip deflections were within +/- 2.7 in. (out of 155-in. max deflection) for 91 percent of the load steps. Experimental testing revealed valuable opportunities for improving the deformed shape equations robustness to real world (not perfect) strain data, which previous analytical testing did not detect. These improvements, which include filtering methods developed in this work, minimize errors due to numerical anomalies discovered in the remaining 9 percent of the load steps. As a result, all load steps attained +/- 2.7 in. accuracy. Wing twist results were very sensitive to errors in bending and require further development. A sensitivity analysis and recommendations for fiber implementation practices, along with, effective filtering methods are included
Bai, D.; Levine, S.L. ); Luoma, J.; Mahgerefteh, M. )
1992-01-01
The Three Mile Island unit 1 core reloads have been designed using fast but accurate scoping codes, PSUI-LEOPARD and ADMARC. PSUI-LEOPARD has been normalized to EPRI-CPM2 results and used to calculate the two-group constants, whereas ADMARC is a modern two-dimensional, two-group diffusion theory nodal code. Problems in accuracy were encountered for cycles 8 and higher as the core lifetime was increased beyond 500 effective full-power days. This is because the heavier loaded cores in both {sup 235}U and {sup 10}B have harder neutron spectra, which produces a change in the transport effect in the baffle reflector region, and the burnable poison (BP) simulations were not accurate enough for the cores containing the increased amount of {sup 10}B required in the BP rods. In the authors study, a technique has been developed to take into account the change in the transport effect in the baffle region by modifying the fast neutron diffusion coefficient as a function of cycle length and core exposure or burnup. A more accurate BP simulation method is also developed, using integral transport theory and CPM2 data, to calculate the BP contribution to the equivalent fuel assembly (supercell) two-group constants. The net result is that the accuracy of the scoping codes is as good as that produced by CASMO/SIMULATE or CPM2/SIMULATE when comparing with measured data.
LATTICE QCD AT FINITE TEMPERATURE AND DENSITY.
BLUM,T.; CREUTZ,M.; PETRECZKY,P.
2004-02-24
With the operation of the RHIC heavy ion program, the theoretical understanding of QCD at finite temperature and density has become increasingly important. Though QCD at finite temperature has been extensively studied using lattice Monte-Carlo simulations over the past twenty years, most physical questions relevant for RHIC (and future) heavy ion experiments remain open. In lattice QCD at finite temperature and density there have been at least two major advances in recent years. First, for the first time calculations of real time quantities, like meson spectral functions have become available. Second, the lattice study of the QCD phase diagram and equation of state have been extended to finite baryon density by several groups. Both issues were extensively discussed in the course of the workshop. A real highlight was the study of the QCD phase diagram in (T, {mu})-plane by Z. Fodor and S. Katz and the determination of the critical end-point for the physical value of the pion mass. This was the first time such lattice calculations at, the physical pion mass have been performed. Results by Z Fodor and S. Katz were obtained using a multi-parameter re-weighting method. Other determinations of the critical end point were also presented, in particular using a Taylor expansion around {mu} = 0 (Bielefeld group, Ejiri et al.) and using analytic continuation from imaginary chemical potential (Ph. de Forcrand and O. Philipsen). The result based on Taylor expansion agrees within errors with the new prediction of Z. Fodor and S. Katz, while methods based on analytic continuation still predict a higher value for the critical baryon density. Most of the thermodynamics studies in full QCD (including those presented at this workshop) have been performed using quite coarse lattices, a = 0.2-0.3 fm. Therefore one may worry about cutoff effects in different thermodynamic quantities, like the transition temperature T{sub tr}. At the workshop U. Heller presented a study of the transition
Li, Bin; Bian, Wensheng
2008-07-14
Full-dimensional quantum calculations of vibrational states of C(2)H(2) and C(2)D(2) are performed in the high-energy region (above 20,400 cm(-1) relative to the acetylene minimum). The theoretical scheme is a combination of several methods. To exploit the full parity and permutation symmetry, the CC-HH diatom-diatom Jacobi coordinates are chosen; phase space optimization in combination with physical considerations is used to obtain an efficient radial discrete variable representation, whereas a basis contraction scheme is applied for angular coordinates. The preconditioned inexact spectral transform method combined with an efficient preconditioner is employed to compute eigenstates within a desired spectral window. The computation is efficient. More definite assignments on vinylidene states than previous studies are acquired using the normal mode projection; in particular, a consistent analysis of the nu(1) (symmetric CH stretch) state is provided. The computed vinylidene vibrational energy levels are in general good agreement with experiment, and several vinylidene states are reported for the first time. PMID:18624520
Vasilkov, Alexander P; Herman, Jay R; Ahmad, Ziauddin; Kahru, Mati; Mitchell, B Greg
2005-05-10
Quantitative assessment of the UV effects on aquatic ecosystems requires an estimate of the in-water radiation field. Actual ocean UV reflectances are needed for improving the total ozone retrievals from the total ozone mapping spectrometer (TOMS) and the ozone monitoring instrument (OMI) flown on NASA's Aura satellite. The estimate of underwater UV radiation can be done on the basis of measurements from the TOMS/OMI and full models of radiative transfer (RT) in the atmosphere-ocean system. The Hydrolight code, modified for extension to the UV, is used for the generation of look-up tables for in-water irradiances. A look-up table for surface radiances generated with a full RT code is input for the Hydrolight simulations. A model of seawater inherent optical properties (IOPs) is an extension of the Case 1 water model to the UV. A new element of the IOP model is parameterization of particulate matter absorption based on recent in situ data. A chlorophyll product from ocean color sensors is input for the IOP model. Verification of the in-water computational scheme shows that the calculated diffuse attenuation coefficient Kd is in good agreement with the measured Kd. PMID:15943340
Full-wave calculation of fast-wave current drive in tokamaks including k sub parallel variations
Jaeger, E.F; Batchelor, D.B.
1991-01-01
When fast waves propagate inward from the edge of a tokamak toward the plasma center, the k{perpendicular} spectrum produced by the antenna is not maintained but is shifted and deformed due to the presence of the finite poloidal magnetic field. This k{perpendicular} shift causes a variation in the parallel phase speed of the wave and can therefore have a strong effect on electron damping and current drive efficiency. In this paper, we include this effect in a new full-wave calculation (PICES) which represents the wave fields as a superposition of poloidal modes, thereby reducing k{perpendicular} to an algebraic operator. The wave equation is solved in general flux coordinates, including a full (non-perturbative) solution for E{perpendicular} and a reduced-order dielectric formulation to eliminate short-wavelength ion Bernstein modes. A simplified current drive model which includes particle trapping is used to estimate the effect of the k{perpendicular} shift on current drive efficiency in ITER and D3-D. Results suggest that when single-pass absorption is weak, reflected power may drive current nearly as efficiently as that absorbed on the first pass. 15 refs., 5 figs.
Experimenting with Langevin lattice QCD
Gavai, R.V.; Potvin, J.; Sanielevici, S.
1987-05-01
We report on the status of our investigations of the effects of systematic errors upon the practical merits of Langevin updating in full lattice QCD. We formulate some rules for the safe use of this updating procedure and some observations on problems which may be common to all approximate fermion algorithms.
Nuclear reactions from lattice QCD
Briceño, Raúl A.; Davoudi, Zohreh; Luu, Thomas C.
2015-01-13
In this study, one of the overarching goals of nuclear physics is to rigorously compute properties of hadronic systems directly from the fundamental theory of strong interactions, Quantum Chromodynamics (QCD). In particular, the hope is to perform reliable calculations of nuclear reactions which will impact our understanding of environments that occur during big bang nucleosynthesis, the evolution of stars and supernovae, and within nuclear reactors and high energy/density facilities. Such calculations, being truly ab initio, would include all two-nucleon and three- nucleon (and higher) interactions in a consistent manner. Currently, lattice QCD provides the only reliable option for performing calculationsmore » of some of the low-energy hadronic observables. With the aim of bridging the gap between lattice QCD and nuclear many-body physics, the Institute for Nuclear Theory held a workshop on Nuclear Reactions from Lattice QCD on March 2013. In this review article, we report on the topics discussed in this workshop and the path planned to move forward in the upcoming years.« less
Foundations of Perturbative QCD
NASA Astrophysics Data System (ADS)
Collins, John
2011-04-01
1. Introduction; 2. Why QCD?; 3. Basics of QCD; 4. Infra-red safety and non-safety; 5. Libby-Sterman analysis and power counting; 6. Parton model to parton theory I; 7. Parton model to parton theory II; 8. Factorization; 9. Corrections to the parton model in QCD; 10. Factorization and subtractions; 11. DIS in QCD; 12. Fragmentation; 13. TMD factorization; 14. Hadron-hadron collisions; 15. More advanced topics; Appendices; References; Index.
Foundations of Perturbative QCD
NASA Astrophysics Data System (ADS)
Collins, John
2013-11-01
1. Introduction; 2. Why QCD?; 3. Basics of QCD; 4. Infra-red safety and non-safety; 5. Libby-Sterman analysis and power counting; 6. Parton model to parton theory I; 7. Parton model to parton theory II; 8. Factorization; 9. Corrections to the parton model in QCD; 10. Factorization and subtractions; 11. DIS in QCD; 12. Fragmentation; 13. TMD factorization; 14. Hadron-hadron collisions; 15. More advanced topics; Appendices; References; Index.
Sudakov safety in perturbative QCD
NASA Astrophysics Data System (ADS)
Larkoski, Andrew J.; Marzani, Simone; Thaler, Jesse
2015-06-01
Traditional calculations in perturbative quantum chromodynamics (pQCD) are based on an order-by-order expansion in the strong coupling αs. Observables that are calculable in this way are known as "safe." Recently, a class of unsafe observables was discovered that do not have a valid αs expansion but are nevertheless calculable in pQCD using all-orders resummation. These observables are called "Sudakov safe" since singularities at each αs order are regulated by an all-orders Sudakov form factor. In this paper, we give a concrete definition of Sudakov safety based on conditional probability distributions, and we study a one-parameter family of momentum sharing observables that interpolate between the safe and unsafe regimes. The boundary between these regimes is particularly interesting, as the resulting distribution can be understood as the ultraviolet fixed point of a generalized fragmentation function, yielding a leading behavior that is independent of αs.
NASA Astrophysics Data System (ADS)
Page, Alexander; Uher, Ctirad; Poudeu, Pierre Ferdinand; Van der Ven, Anton
2015-11-01
Previous studies have indicated that the figure of merit (ZT ) of half-Heusler (HH) alloys with composition M NiSn (M =Ti , Zr, or Hf) is greatly enhanced when the alloys contain a nano-scale full-Heusler (FH) MN i2Sn second phase. However, the formation mechanism of the FHnanostructures in the HH matrix and their vibrational properties are still not well understood. We report on first-principles studies of thermodynamic phase equilibria in the MNiSn-MN i2Sn pseudobinary system as well as HH and FH vibrational properties. Thermodynamic phase diagrams as functions of temperature and Ni concentration were developed using density functional theory (DFT) combined with a cluster expansion and Monte Carlo simulations. The phase diagrams show very low excess Ni solubility in HH alloys even at high temperatures, which indicates that any Ni excess will decompose into a two-phase mixture of HH and FH compounds. Vibrational properties of HH and FH alloys are compared. Imaginary vibrational modes in the calculated phonon dispersion diagram of TiN i2Sn indicate a dynamical instability with respect to cubic [001] transverse acoustic modulations. Displacing atoms along unstable vibrational modes in cubic TiN i2Sn reveals lower-energy structures with monoclinic symmetry. The energy of the monoclinic structures is found to depend strongly on the lattice parameter. The origin of the instability in cubic TiN i2Sn and its absence in cubic ZrN i2Sn and HfN i2Sn is attributed to the small size of the Ti 3 d shells compared to those of Zr and Hf atoms. Lattice constants and heat capacities calculated by DFT agree well with experiment.
Excited light isoscalar mesons from lattice QCD
Christopher Thomas
2011-07-01
I report a recent lattice QCD calculation of an excited spectrum of light isoscalar mesons, something that has up to now proved challenging for lattice QCD. With novel techniques we extract an extensive spectrum with high statistical precision, including spin-four states and, for the first time, light isoscalars with exotic quantum numbers. In addition, the hidden flavour content of these mesons is determined, providing a window on annihilation dynamics in QCD. I comment on future prospects including applications to the study of resonances.
QCD thermodynamics and missing hadron states
NASA Astrophysics Data System (ADS)
Petreczky, Peter
2016-03-01
Equation of State and fluctuations of conserved charges in hot strongly interacting matter are being calculated with increasing accuracy in lattice QCD, and continuum results at physical quark masses become available. At sufficiently low temperature the thermodynamic quantities can be understood in terms of hadron resonance gas model that includes known hadrons and hadronic resonances from Particle Data Book. However, for some quantities it is necessary to include undiscovered hadronic resonances (missing states) that are, however, predicted by quark model and lattice QCD study of hadron spectrum. Thus, QCD thermodynamics can provide indications for the existence of yet undiscovered hadron states.
Death to perturbative QCD in exclusive processes?
Eckardt, R.; Hansper, J.; Gari, M.F.
1994-04-01
The authors discuss the question of whether perturbative QCD is applicable in calculations of exclusive processes at available momentum transfers. They show that the currently used method of determining hadronic quark distribution amplitudes from QCD sum rules yields wave functions which are completely undetermined because the polynomial expansion diverges. Because of the indeterminacy of the wave functions no statement can be made at present as to whether perturbative QCD is valid. The authors emphasize the necessity of a rigorous discussion of the subject and the importance of experimental data in the range of interest.
The QCD vacuum, hadrons and superdense matter
Shuryak, E.
1986-01-01
This is probably the only textbook available that gathers QCD, many-body theory and phase transitions in one volume. The presentation is pedagogical and readable. Contents: The QCD Vacuum: Introduction; QCD on the Lattice Topological Effects in Gauges Theories. Correlation Functions and Microscopic Excitations: Introduction; Operator Product Expansion; The Sum Rules beyond OPE; Nonpower Contributions to Correlators and Instantons; Hadronic Spectroscopy on the Lattice. Dense Matter: Hadronic Matter; Asymptotically Dense Quark-Gluon Plasma; Instantons in Matter; Lattice Calculations at Finite Temperature; Phase Transitions; Macroscopic Excitations and Experiments: General Properties of High Energy Collisions; ''Barometers'', ''Thermometers'', Interferometric ''Microscope''; Experimental Perspectives.
Excited light meson spectroscopy from lattice QCD
Christopher Thomas, Hadron Spectrum Collaboration
2012-04-01
I report on recent progress in calculating excited meson spectra using lattice QCD, emphasizing results and phenomenology. With novel techniques we can now extract extensive spectra of excited mesons with high statistical precision, including spin-four states and those with exotic quantum numbers. As well as isovector meson spectra, I will present new calculations of the spectrum of excited light isoscalar mesons, something that has up to now been a challenge for lattice QCD. I show determinations of the flavor content of these mesons, including the eta-eta' mixing angle, providing a window on annihilation dynamics in QCD. I will also discuss recent work on using lattice QCD to map out the energy-dependent phase shift in pi-pi scattering and future applications of the methodology to the study of resonances and decays.
Lattice QCD input for axion cosmology
NASA Astrophysics Data System (ADS)
Berkowitz, Evan; Buchoff, Michael I.; Rinaldi, Enrico
2015-08-01
One intriguing beyond-the-Standard-Model particle is the QCD axion, which could simultaneously provide a solution to the Strong C P Problem and account for some, if not all, of the dark matter density in the Universe. This particle is a pseudo-Nambu-Goldstone boson of the conjectured Peccei-Quinn symmetry of the Standard Model. Its mass and interactions are suppressed by a heavy symmetry-breaking scale, fa, the value of which is roughly greater than 109 GeV (or, conversely, the axion mass, ma, is roughly less than 104 μ eV ). The density of axions in the Universe, which cannot exceed the relic dark matter density and is a quantity of great interest in axion experiments like ADMX, is a result of the early Universe interplay between cosmological evolution and the axion mass as a function of temperature. The latter quantity is proportional to the second derivative of the temperature-dependent QCD free energy with respect to the C P -violating phase, θ . However, this quantity is generically nonperturbative, and previous calculations have only employed instanton models at the high temperatures of interest (roughly 1 GeV). In this and future works, we aim to calculate the temperature-dependent axion mass at small θ from first-principle lattice calculations, with controlled statistical and systematic errors. Once calculated, this temperature-dependent axion mass is input for the classical evolution equations of the axion density of the Universe, which is required to be less than or equal to the dark matter density. Due to a variety of lattice systematic effects at the very high temperatures required, we perform a calculation of the leading small-θ cumulant of the theta vacua on large volume lattices for SU(3) Yang-Mills with high statistics as a first proof of concept, before attempting a full QCD calculation in the future. From these pure glue results, the misalignment mechanism yields the axion mass bound ma≥(14.6 ±0.1 ) μ eV when Peccei-Quinn breaking occurs
Thermodynamics of Strong-Interaction Matter from Lattice QCD
NASA Astrophysics Data System (ADS)
Ding, Heng-Tong; Karsch, Frithjof; Mukherjee, Swagato
We review results from lattice QCD calculations on the thermodynamics of strong-interaction matter with emphasis on input these calculations can provide to the exploration of the phase diagram and properties of hot and dense matter created in heavy ion experiments. This review is organized in sections as follows: (1) Introduction, (2) QCD thermodynamics on the lattice, (3) QCD phase diagram at high temperature, (4) Bulk thermodynamics, (5) Fluctuations of conserved charges, (6) Transport properties, (7) Open heavy flavors and heavy quarkonia, (8) QCD in external magnetic fields, (9) Summary.
Vranas, P
2007-06-18
Quantum Chromodynamics is the theory of nuclear and sub-nuclear physics. It is a celebrated theory and one of its inventors, F. Wilczek, has termed it as '... our most perfect physical theory'. Part of this is related to the fact that QCD can be numerically simulated from first principles using the methods of lattice gauge theory. The computational demands of QCD are enormous and have not only played a role in the history of supercomputers but are also helping define their future. Here I will discuss the intimate relation of QCD and massively parallel supercomputers with focus on the Blue Gene supercomputer and QCD thermodynamics. I will present results on the performance of QCD on the Blue Gene as well as physics simulation results of QCD at temperatures high enough that sub-nuclear matter transitions to a plasma state of elementary particles, the quark gluon plasma. This state of matter is thought to have existed at around 10 microseconds after the big bang. Current heavy ion experiments are in the quest of reproducing it for the first time since then. And numerical simulations of QCD on the Blue Gene systems are calculating the theoretical values of fundamental parameters so that comparisons of experiment and theory can be made.
NASA Astrophysics Data System (ADS)
Larsen, Helena; Olsen, Jeppe; Hättig, Christof; Jørgensen, Poul; Christiansen, Ove; Gauss, Jürgen
1999-08-01
Static and frequency-dependent polarizabilities and first hyperpolarizabilities have been calculated for HF and Ne using full configuration interaction (FCI) and a hierarchy of coupled cluster models: coupled cluster singles (CCS), an approximate coupled cluster singles and doubles model (CC2), coupled cluster singles and doubles (CCSD), an approximate coupled cluster singles, doubles, and triples model (CC3), and coupled cluster singles, doubles, and triples (CCSDT). A previous study of BH concerning FCI benchmarking has been extended to include CC3 and static CCSDT values. Systematic improvements of the polarizabilities and the hyperpolarizabilities are found going from CCS to CCSD and from CCSD to CC3 or CCSDT. Little or no improvement of the polarizabilities and no improvement of the hyperpolarizabilities are seen when going from CCS to CC2. The CCSD results represent a significant improvement over CCS and CC2 but are again surpassed by the CC3 results which agree very well with the FCI values. The relative error for the static polarizability at the CC3 level is 0.11% for Ne and, respectively, 0.16% and 0.20% for αxx and αzz of HF. For βzzz and βzxx the errors are 0.50% and 1.7%, respectively. Only in the challenging case of BH does CCSDT improve the CC3 values. The dispersion for the polarizabilities and hyperpolarizabilities is predicted with increasing accuracy in the CCS-CC2-CCSD-CC3 sequence as expected from the increasing accuracy of the electronic excitation energies. For all molecules the effect of orbital relaxation has been investigated for the static properties. The inclusion of orbital relaxation gives results that are somewhat different from the unrelaxed results but are in general no improvement.
Towards the chiral limit in QCD
Shailesh Chandrasekharan
2006-02-28
Computing hadronic observables by solving QCD from first principles with realistic quark masses is an important challenge in fundamental nuclear and particle physics research. Although lattice QCD provides a rigorous framework for such calculations many difficulties arise. Firstly, there are no good algorithms to solve lattice QCD with realistically light quark masses. Secondly, due to critical slowing down, Monte Carlo algorithms are able to access only small lattice sizes on coarse lattices. Finally, due to sign problems it is almost impossible to study the physics of finite baryon density. Lattice QCD contains roughly three mass scales: the cutoff (or inverse lattice spacing) a{sup -1}, the confinement scale {Lambda}{sub QCD}, and the pion mass m{sub {pi}}. Most conventional Monte Carlo algorithms for QCD become inefficient in two regimes: when {Lambda}{sub QCD} becomes small compared to a{sup -1} and when m{sub {pi}} becomes small compared to {Lambda}{sub QCD}. The former can be largely controlled by perturbation theory thanks to asymptotic freedom. The latter is more difficult since chiral extrapolations are typically non-analytic and can be unreliable if the calculations are not done at sufficiently small quark masses. For this reason it has been difficult to compute quantities close to the chiral limit. The essential goal behind this proposal was to develop a new approach towards understanding QCD and QCD-like theories with sufficiently light quarks. The proposal was based on a novel cluster algorithm discovered in the strong coupling limit with staggered fermions [1]. This algorithm allowed us to explore the physics of exactly massless quarks and as well as light quarks. Thus, the hope was that this discovery would lead to the complete solution of at least a few strongly coupled QCD-like theories. The solution would be far better than those achievable through conventional methods and thus would be able to shed light on the chiral physics from a new
NASA Astrophysics Data System (ADS)
Stancil, Phillip
We propose to compute accurate collisional excitation rate coefficients for rovibrational transitions of CS, SiO, SO, NO, H_2O, and HCN due to H_2, He, or H impact. This extends our previous grant which focused on 3- and 4-atom systems to 4- and 5-atom collision complexes, with dynamics to be performed on 6-9 dimensional potential energy surfaces (PESs). This work, which uses fully quantum mechanical methods for inelastic scattering and incorporates full-dimensional PESs, pushes beyond the state-of-the-art for such calculations, as recently established by our group for rovibrational transitions in CO-H_2 in 6D. Many of the required PESs will be computed as part of this project using ab initio theory and basis sets of the highest level feasible and particular attention will be given to the long range form of the PESs. The completion of the project will result in 6 new global PESs and state-to-state rate coefficients for a large range of initial rovibrational levels for temperatures between 1 and 3000 K. The chosen collision systems correspond to cases where data are limited or lacking, are important coolants or diagnostics, and result in observable emission features in the infrared (IR). The final project results will be important for the analysis of a variety of interstellar and extragalactic environments in which the local conditions of gas density, radiation field, and/or shocks drive the level populations out of equilibrium. In such cases, collisional excitation data are critical to the accurate prediction and interpretation of observed molecular IR emission lines in protoplanetary disks, star-forming regions, planetary nebulae, embedded protostars, photodissociation regions, etc. The use of the proposed collisional excitation data will lead to deeper examination and understanding of the properties of many astrophysical environments, hence elevating the scientific return from the upcoming JWST, as well as from current (SOFIA, Herschel, HST) and past IR missions
Heavy Quarks, QCD, and Effective Field Theory
Thomas Mehen
2012-10-09
The research supported by this OJI award is in the area of heavy quark and quarkonium production, especially the application Soft-Collinear E ective Theory (SCET) to the hadronic production of quarkonia. SCET is an e ffective theory which allows one to derive factorization theorems and perform all order resummations for QCD processes. Factorization theorems allow one to separate the various scales entering a QCD process, and in particular, separate perturbative scales from nonperturbative scales. The perturbative physics can then be calculated using QCD perturbation theory. Universal functions with precise fi eld theoretic de nitions describe the nonperturbative physics. In addition, higher order perturbative QCD corrections that are enhanced by large logarithms can be resummed using the renormalization group equations of SCET. The applies SCET to the physics of heavy quarks, heavy quarkonium, and similar particles.
NASA Astrophysics Data System (ADS)
Wilczek, Frank
Introduction Symmetry and the Phenomena of QCD Apparent and Actual Symmetries Asymptotic Freedom Confinement Chiral Symmetry Breaking Chiral Anomalies and Instantons High Temperature QCD: Asymptotic Properties Significance of High Temperature QCD Numerical Indications for Quasi-Free Behavior Ideas About Quark-Gluon Plasma Screening Versus Confinement Models of Chiral Symmetry Breaking More Refined Numerical Experiments High-Temperature QCD: Phase Transitions Yoga of Phase Transitions and Order Parameters Application to Glue Theories Application to Chiral Transitions Close Up on Two Flavors A Genuine Critical Point! (?) High-Density QCD: Methods Hopes, Doubts, and Fruition Another Renormalization Group Pairing Theory Taming the Magnetic Singularity High-Density QCD: Color-Flavor Locking and Quark-Hadron Continuity Gauge Symmetry (Non)Breaking Symmetry Accounting Elementary Excitations A Modified Photon Quark-Hadron Continuity Remembrance of Things Past More Quarks Fewer Quarks and Reality
Urban, Federico R.; Zhitnitsky, Ariel R.
2010-08-30
We review two mechanisms rooted in the infrared sector of QCD which, by exploiting the properties of the QCD ghost, as introduced by Veneziano, provide new insight on the cosmological dark energy problem, first, in the form of a Casimir-like energy from quantising QCD in a box, and second, in the form of additional, time-dependent, vacuum energy density in an expanding universe. Based on [1, 2].
The CKM Matrix from Lattice QCD
Mackenzie, Paul B.; /Fermilab
2009-07-01
Lattice QCD plays an essential role in testing and determining the parameters of the CKM theory of flavor mixing and CP violation. Very high precisions are required for lattice calculations analyzing CKM data; I discuss the prospects for achieving them. Lattice calculations will also play a role in investigating flavor mixing and CP violation beyond the Standard Model.
Schvellinger, Martin
2008-07-28
We briefly review one of the current applications of the AdS/CFT correspondence known as AdS/QCD and discuss about the calculation of four-point quark-flavour current correlation functions and their applications to the calculation of observables related to neutral kaon decays and neutral kaon mixing processes.
QCD spectrum with three quark flavors
Bernard, Claude; Burch, Tom; Orginos, Kostas; Toussaint, Doug; DeGrand, Thomas A.; DeTar, Carleton; Datta, Saumen; Gottlieb, Steven; Heller, Urs M.; Sugar, Robert
2001-09-01
We present results from a lattice hadron spectrum calculation using three flavors of dynamical quarks -- two light and one strange -- and quenched simulations for comparison. These simulations were done using a one-loop Symanzik improved gauge action and an improved Kogut-Susskind quark action. The lattice spacings, and hence also the physical volumes, were tuned to be the same in all the runs to better expose differences due to flavor number. Lattice spacings were tuned using the static quark potential, so as a by-product we obtain updated results for the effect of sea quarks on the static quark potential. We find indications that the full QCD meson spectrum is in better agreement with experiment than the quenched spectrum. For the 0{sup ++} (a{sub 0}) meson we see a coupling to two pseudoscalar mesons, or a meson decay on the lattice.
Nucleon QCD sum rules in the instanton medium
Ryskin, M. G.; Drukarev, E. G. Sadovnikova, V. A.
2015-09-15
We try to find grounds for the standard nucleon QCD sum rules, based on a more detailed description of the QCD vacuum. We calculate the polarization operator of the nucleon current in the instanton medium. The medium (QCD vacuum) is assumed to be a composition of the small-size instantons and some long-wave gluon fluctuations. We solve the corresponding QCD sum rule equations and demonstrate that there is a solution with the value of the nucleon mass close to the physical one if the fraction of the small-size instantons contribution is w{sub s} ≈ 2/3.
LATTICE QCD AT FINITE DENSITY.
SCHMIDT, C.
2006-07-23
I discuss different approaches to finite density lattice QCD. In particular, I focus on the structure of the phase diagram and discuss attempts to determine the location of the critical end-point. Recent results on the transition line as function of the chemical potential (T{sub c}({mu}{sub q})) are reviewed. Along the transition line, hadronic fluctuations have been calculated; which can be used to characterize properties of the Quark Gluon plasma and eventually can also help to identify the location of the critical end-point in the QCD phase diagram on the lattice and in heavy ion experiments. Furthermore, I comment on the structure of the phase diagram at large {mu}{sub q}.
Nuclear Physics from Lattice QCD
William Detmold, Silas Beane, Konstantinos Orginos, Martin Savage
2011-01-01
We review recent progress toward establishing lattice Quantum Chromodynamics as a predictive calculational framework for nuclear physics. A survey of the current techniques that are used to extract low-energy hadronic scattering amplitudes and interactions is followed by a review of recent two-body and few-body calculations by the NPLQCD collaboration and others. An outline of the nuclear physics that is expected to be accomplished with Lattice QCD in the next decade, along with estimates of the required computational resources, is presented.
QCD and Top-Quark Results from the Tevatron
Zielinski, Marek; /Rochester U.
2006-10-01
Selected recent QCD and top-quark results from the Tevatron are reviewed, aiming to illustrate progression from basic studies of QCD processes to verification of perturbative calculations and Monte Carlo simulation tools, and to their applications in more novel and complex cases, like top-quark studies and searches for new physics.
The QCD evolution of TMD in the covariant approach
NASA Astrophysics Data System (ADS)
Efremov, A. V.; Teryaev, O. V.; Zavada, P.
2016-02-01
The procedure for calculation of the QCD evolution of transverse momentum dependent distributions within the covariant approach is suggested. The standard collinear QCD evolution together with the requirements of relativistic invariance and rotational symmetry of the nucleon in its rest frame represent the basic ingredients of our approach. The obtained results are compared with the predictions of some other approaches.
NASA Astrophysics Data System (ADS)
Ye, Lin-Hui
2015-09-01
Although the supercell method has been widely used for surface calculations, it only works well with short-ranged potentials, but meets difficulty when the potential decays very slowly into the vacuum. Unfortunately, the exact exchange-correlation potential of the density functional theory is asymptotically long ranged, and therefore is not easily handled by use of supercells. This paper illustrates that the authentic slab geometry, another technique for surface calculations, is not affected by this issue: It works equally well with both short- and long-ranged potentials, with the computational cost and the convergence speed being essentially the same. Using the asymptotically long-ranged Becke-Roussel'89 exchange potential as an example, we have calculated six surfaces of various types. We found that accurate potential values can be obtained even in extremely low density regions of more than 100 Å away from the surface. This high performance allows us to explore the asymptotic region, and prove with clean numerical evidence that the Becke-Roussel'89 potential satisfies the correct asymptotic behavior for slab surfaces, as it does for finite systems. Our finding further implies that the Slater component of the exact exchange optimized effective potential is responsible for the asymptotic behavior, not only for jellium slabs, but for slabs of any type. The Becke-Roussel'89 potential may therefore be used to build asymptotically correct model exchange potentials applicable to both finite systems and slab surfaces.
Some New/Old Approaches to QCD
DOE R&D Accomplishments Database
Gross, D. J.
1992-11-01
In this lecture I shall discuss some recent attempts to revive some old ideas to address the problem of solving QCD. I believe that it is timely to return to this problem which has been woefully neglected for the last decade. QCD is a permanent part of the theoretical landscape and eventually we will have to develop analytic tools for dealing with the theory in the infra-red. Lattice techniques are useful but they have not yet lived up to their promise. Even if one manages to derive the hadronic spectrum numerically, to an accuracy of 10% or even 1%, we will not be truly satisfied unless we have some analytic understanding of the results. Also, lattice Monte-Carlo methods can only be used to answer a small set of questions. Many issues of great conceptual and practical interest-in particular the calculation of scattering amplitudes, are thus far beyond lattice control. Any progress in controlling QCD in an explicit analytic, fashion would be of great conceptual value. It would also be of great practical aid to experimentalists, who must use rather ad-hoc and primitive models of QCD scattering amplitudes to estimate the backgrounds to interesting new physics. I will discuss an attempt to derive a string representation of QCD and a revival of the large N approach to QCD. Both of these ideas have a long history, many theorist-years have been devoted to their pursuit-so far with little success. I believe that it is time to try again. In part this is because of the progress in the last few years in string theory. Our increased understanding of string theory should make the attempt to discover a stringy representation of QCD easier, and the methods explored in matrix models might be employed to study the large N limit of QCD.
Some new/old approaches to QCD
Gross, D.J.
1992-11-01
In this lecture I shall discuss some recent attempts to revive some old ideas to address the problem of solving QCD. I believe that it is timely to return to this problem which has been woefully neglected for the last decade. QCD is a permanent part of the theoretical landscape and eventually we will have to develop analytic tools for dealing with the theory in the infra-red. Lattice techniques are useful but they have not yet lived up to their promise. Even if one manages to derive the hadronic spectrum numerically, to an accuracy of 10% or even 1%, we will not be truly satisfied unless we have some analytic understanding of the results. Also, lattice Monte-Carlo methods can only be used to answer a small set of questions. Many issues of great conceptual and practical interest-in particular the calculation of scattering amplitudes, are thus far beyond lattice control. Any progress in controlling QCD in an explicit analytic, fashion would be of great conceptual value. It would also be of great practical aid to experimentalists, who must use rather ad-hoc and primitive models of QCD scattering amplitudes to estimate the backgrounds to interesting new physics. I will discuss an attempt to derive a string representation of QCD and a revival of the large N approach to QCD. Both of these ideas have a long history, many theorist-years have been devoted to their pursuit-so far with little success. I believe that it is time to try again. In part this is because of the progress in the last few years in string theory. Our increased understanding of string theory should make the attempt to discover a stringy representation of QCD easier, and the methods explored in matrix models might be employed to study the large N limit of QCD.
Norniella, Olga; /Barcelona, IFAE
2005-01-01
Recent QCD measurements from the CDF collaboration at the Tevatron are presented, together with future prospects as the luminosity increases. The measured inclusive jet cross section is compared to pQCD NLO predictions. Precise measurements on jet shapes and hadronic energy flows are compared to different phenomenological models that describe gluon emissions and the underlying event in hadron-hadron interactions.
Transport coefficients in Yang-Mills theory and QCD.
Christiansen, Nicolai; Haas, Michael; Pawlowski, Jan M; Strodthoff, Nils
2015-09-11
We calculate the shear-viscosity-over-entropy-density ratio η/s in Yang-Mills theory from the Kubo formula using an exact diagrammatic representation in terms of full propagators and vertices using gluon spectral functions as external input. We provide an analytic fit formula for the temperature dependence of η/s over the whole temperature range from a glueball resonance gas at low temperatures, to a high-temperature regime consistent with perturbative results. Subsequently, we provide a first estimate for η/s in QCD. PMID:26406822
Exploring hyperons and hypernuclei with lattice QCD
Beane, S.R.; Bedaque, P.F.; Parreno, A.; Savage, M.J.
2003-01-01
In this work we outline a program for lattice QCD that wouldprovide a first step toward understanding the strong and weakinteractions of strange baryons. The study of hypernuclear physics hasprovided a significant amount of information regarding the structure andweak decays of light nuclei containing one or two Lambda's, and Sigma's.From a theoretical standpoint, little is known about the hyperon-nucleoninteraction, which is required input for systematic calculations ofhypernuclear structure. Furthermore, the long-standing discrepancies inthe P-wave amplitudes for nonleptonic hyperon decays remain to beunderstood, and their resolution is central to a better understanding ofthe weak decays of hypernuclei. We present a framework that utilizesLuscher's finite-volume techniques in lattice QCD to extract thescattering length and effective range for Lambda-N scattering in both QCDand partially-quenched QCD. The effective theory describing thenonleptonic decays of hyperons using isospin symmetry alone, appropriatefor lattice calculations, is constructed.
NASA Technical Reports Server (NTRS)
Bui, Trong T.
1993-01-01
New turbulence modeling options recently implemented for the 3-D version of Proteus, a Reynolds-averaged compressible Navier-Stokes code, are described. The implemented turbulence models include: the Baldwin-Lomax algebraic model, the Baldwin-Barth one-equation model, the Chien k-epsilon model, and the Launder-Sharma k-epsilon model. Features of this turbulence modeling package include: well documented and easy to use turbulence modeling options, uniform integration of turbulence models from different classes, automatic initialization of turbulence variables for calculations using one- or two-equation turbulence models, multiple solid boundaries treatment, and fully vectorized L-U solver for one- and two-equation models. Validation test cases include the incompressible and compressible flat plate turbulent boundary layers, turbulent developing S-duct flow, and glancing shock wave/turbulent boundary layer interaction. Good agreement is obtained between the computational results and experimental data. Sensitivity of the compressible turbulent solutions with the method of y(sup +) computation, the turbulent length scale correction, and some compressibility corrections are examined in detail. The test cases show that the highly optimized one-and two-equation turbulence models can be used in routine 3-D Navier-Stokes computations with no significant increase in CPU time as compared with the Baldwin-Lomax algebraic model.
NASA Technical Reports Server (NTRS)
Bui, Trong T.
1993-01-01
New turbulence modeling options recently implemented for the 3D version of Proteus, a Reynolds-averaged compressible Navier-Stokes code, are described. The implemented turbulence models include: the Baldwin-Lomax algebraic model, the Baldwin-Barth one-equation model, the Chien k-epsilon model, and the Launder-Sharma k-epsilon model. Features of this turbulence modeling package include: well documented and easy to use turbulence modeling options, uniform integration of turbulence models from different classes, automatic initialization of turbulence variables for calculations using one- or two-equation turbulence models, multiple solid boundaries treatment, and fully vectorized L-U solver for one- and two-equation models. Good agreements are obtained between the computational results and experimental data. Sensitivity of the compressible turbulent solutions with the method of y(+) computation, the turbulent length scale correction, and some compressibility corrections are examined in detail. Test cases show that the highly optimized one- and two-equation turbulence models can be used in routine 3D Navier-Stokes computations with no significant increase in CPU time as compared with the Baldwin-Lomax algebraic model.
NASA Astrophysics Data System (ADS)
Ono, Tomoya; Heide, Marcus; Atodiresei, Nicolae; Baumeister, Paul; Tsukamoto, Shigeru; Blügel, Stefan
2010-11-01
We have developed an efficient computational scheme utilizing the real-space finite-difference formalism and the projector augmented-wave (PAW) method to perform precise first-principles electronic-structure simulations based on the density-functional theory for systems containing transition metals with a modest computational effort. By combining the advantages of the time-saving double-grid technique and the Fourier-filtering procedure for the projectors of pseudopotentials, we can overcome the egg box effect in the computations even for first-row elements and transition metals, which is a problem of the real-space finite-difference formalism. In order to demonstrate the potential power in terms of precision and applicability of the present scheme, we have carried out simulations to examine several bulk properties and structural energy differences between different bulk phases of transition metals and have obtained excellent agreement with the results of other precise first-principles methods such as a plane-wave-based PAW method and an all-electron full-potential linearized augmented plane-wave (FLAPW) method.
Andersen, Jens O.; Leganger, Lars E.; Strickland, Michael; Su, Nan
2011-10-15
In this brief report we compare the predictions of a recent next-to-next-to-leading order hard-thermal-loop perturbation theory (HTLpt) calculation of the QCD trace anomaly to available lattice data. We focus on the trace anomaly scaled by T{sup 2} in two cases: N{sub f}=0 and N{sub f}=3. When using the canonical value of {mu}=2{pi}T for the renormalization scale, we find that for Yang-Mills theory (N{sub f}=0) agreement between HTLpt and lattice data for the T{sup 2}-scaled trace anomaly begins at temperatures on the order of 8T{sub c}, while treating the subtracted piece as an interaction term when including quarks (N{sub f}=3) agreement begins already at temperatures above 2T{sub c}. In both cases we find that at very high temperatures the T{sup 2}-scaled trace anomaly increases with temperature in accordance with the predictions of HTLpt.
Brodsky, Stanley J.; de Teramond, Guy F.; /Costa Rica U. /SLAC
2007-02-21
The AdS/CFT correspondence between string theory in AdS space and conformal .eld theories in physical spacetime leads to an analytic, semi-classical model for strongly-coupled QCD which has scale invariance and dimensional counting at short distances and color confinement at large distances. Although QCD is not conformally invariant, one can nevertheless use the mathematical representation of the conformal group in five-dimensional anti-de Sitter space to construct a first approximation to the theory. The AdS/CFT correspondence also provides insights into the inherently non-perturbative aspects of QCD, such as the orbital and radial spectra of hadrons and the form of hadronic wavefunctions. In particular, we show that there is an exact correspondence between the fifth-dimensional coordinate of AdS space z and a specific impact variable {zeta} which measures the separation of the quark and gluonic constituents within the hadron in ordinary space-time. This connection allows one to compute the analytic form of the frame-independent light-front wavefunctions, the fundamental entities which encode hadron properties and allow the computation of decay constants, form factors, and other exclusive scattering amplitudes. New relativistic lightfront equations in ordinary space-time are found which reproduce the results obtained using the 5-dimensional theory. The effective light-front equations possess remarkable algebraic structures and integrability properties. Since they are complete and orthonormal, the AdS/CFT model wavefunctions can also be used as a basis for the diagonalization of the full light-front QCD Hamiltonian, thus systematically improving the AdS/CFT approximation.
Nucleon Structure from Lattice QCD
Haegler, Philipp
2011-10-24
Hadron structure calculations in lattice QCD have seen substantial progress during recent years. We illustrate the achievements that have been made by discussing latest lattice results for a limited number of important observables related to nucleon form factors and generalized parton distributions. A particular focus is placed on the decomposition of the nucleon spin 1/2 in terms of quark spin and orbital angular momentum contributions. Results and limitations of the necessary chiral extrapolations based on ChPT will be briefly discussed.
Varelas, N.; D0 Collaboration
1997-10-01
We present recent results on jet production, dijet angular distributions, W+ Jets, and color coherence from p{anti p} collisions at {radical}s = 1.8 TeV at the Fermilab Tevatron Collider using the D0 detector. The data are compared to perturbative QCD calculations or to predictions of parton shower based Monte Carlo models.
The Chroma Software System for Lattice QCD
Robert Edwards; Balint Joo
2004-06-01
We describe aspects of the Chroma software system for lattice QCD calculations. Chroma is an open source C++ based software system developed using the software infrastructure of the US SciDAC initiative. Chroma interfaces with output from the BAGEL assembly generator for optimized lattice fermion kernels on some architectures. It can be run on workstations, clusters and the QCDOC supercomputer.
Effective charges and expansion parameters in QCD
Braaten, E.; Leveille, J.P.
1981-10-01
The momentum subtraction scheme MOM has been empirically successful in producing small QCD corrections to physical quantities at one loop order. By explicit calculations, we show that with a suitable shift in the renormalization scale, the minimal subtraction scheme coupling constant ..cap alpha../sub MS/ coincides with typical momentum scheme coupling constants at both one and two loop order.
Heavy quark masses from lattice QCD
NASA Astrophysics Data System (ADS)
Lytle, Andrew T.
2016-07-01
Progress in quark mass determinations from lattice QCD is reviewed, focusing on results for charm and bottom mass. These are of particular interest for precision Higgs studies. Recent determinations have achieved percent-level uncertainties with controlled systematics. Future prospects for these calculations are also discussed.
Hadronic and nuclear interactions in QCD
Not Available
1982-01-01
Despite the evidence that QCD - or something close to it - gives a correct description of the structure of hadrons and their interactions, it seems paradoxical that the theory has thus far had very little impact in nuclear physics. One reason for this is that the application of QCD to distances larger than 1 fm involves coherent, non-perturbative dynamics which is beyond present calculational techniques. For example, in QCD the nuclear force can evidently be ascribed to quark interchange and gluon exchange processes. These, however, are as complicated to analyze from a fundamental point of view as is the analogous covalent bond in molecular physics. Since a detailed description of quark-quark interactions and the structure of hadronic wavefunctions is not yet well-understood in QCD, it is evident that a quantitative first-principle description of the nuclear force will require a great deal of theoretical effort. Another reason for the limited impact of QCD in nuclear physics has been the conventional assumption that nuclear interactions can for the most part be analyzed in terms of an effective meson-nucleon field theory or potential model in isolation from the details of short distance quark and gluon structure of hadrons. These lectures, argue that this view is untenable: in fact, there is no correspondence principle which yields traditional nuclear physics as a rigorous large-distance or non-relativistic limit of QCD dynamics. On the other hand, the distinctions between standard nuclear physics dynamics and QCD at nuclear dimensions are extremely interesting and illuminating for both particle and nuclear physics.
QCD on the Massively Parallel Computer AP1000
NASA Astrophysics Data System (ADS)
Akemi, K.; Fujisaki, M.; Okuda, M.; Tago, Y.; Hashimoto, T.; Hioki, S.; Miyamura, O.; Takaishi, T.; Nakamura, A.; de Forcrand, Ph.; Hege, C.; Stamatescu, I. O.
We present the QCD-TARO program of calculations which uses the parallel computer AP1000 of Fujitsu. We discuss the results on scaling, correlation times and hadronic spectrum, some aspects of the implementation and the future prospects.
Recent progress in lattice QCD at finite temperature
Petreczky,P.
2009-02-01
I review recent progress in finite temperature lattice calculations,including the study of the nature of the deconfinement transition in QCD, equation of state, screening of static quarks and meson spectral functions.
Status of Average-x from Lattice QCD
Dru Renner
2011-09-01
As algorithms and computing power have advanced, lattice QCD has become a precision technique for many QCD observables. However, the calculation of nucleon matrix elements remains an open challenge. I summarize the status of the lattice effort by examining one observable that has come to represent this challenge, average-x: the fraction of the nucleon's momentum carried by its quark constituents. Recent results confirm a long standing tendency to overshoot the experimentally measured value. Understanding this puzzle is essential to not only the lattice calculation of nucleon properties but also the broader effort to determine hadron structure from QCD.
NASA Astrophysics Data System (ADS)
Lutz, Matthias F. M.; Lange, Jens Sören; Pennington, Michael; Bettoni, Diego; Brambilla, Nora; Crede, Volker; Eidelman, Simon; Gillitzer, Albrecht; Gradl, Wolfgang; Lang, Christian B.; Metag, Volker; Nakano, Takashi; Nieves, Juan; Neubert, Sebastian; Oka, Makoto; Olsen, Stephen L.; Pappagallo, Marco; Paul, Stephan; Pelizäus, Marc; Pilloni, Alessandro; Prencipe, Elisabetta; Ritman, Jim; Ryan, Sinead; Thoma, Ulrike; Uwer, Ulrich; Weise, Wolfram
2016-04-01
We report on the EMMI Rapid Reaction Task Force meeting 'Resonances in QCD', which took place at GSI October 12-14, 2015. A group of 26 people met to discuss the physics of resonances in QCD. The aim of the meeting was defined by the following three key questions: What is needed to understand the physics of resonances in QCD? Where does QCD lead us to expect resonances with exotic quantum numbers? What experimental efforts are required to arrive at a coherent picture? For light mesons and baryons only those with up, down and strange quark content were considered. For heavy-light and heavy-heavy meson systems, those with charm quarks were the focus. This document summarizes the discussions by the participants, which in turn led to the coherent conclusions we present here.
NASA Astrophysics Data System (ADS)
Shinotsuka, H.; Tanuma, S.; Powell, C. J.; Penn, D. R.
2012-01-01
We present mass collision electron stopping powers (SPs) for 41 elemental solids (Li, Be, graphite, diamond, glassy C, Na, Mg, Al, Si, K, Sc, Ti, V, Cr, Fe, Co, Ni, Cu, Ge, Y, Nb, Mo, Ru, Rh, Pd, Ag, In, Sn, Cs, Gd, Tb, Dy, Hf, Ta, W, Re, Os, Ir, Pt, Au, and Bi) that were calculated from experimental energy-loss-function data with the full Penn algorithm for electron energies between 50 eV and 30 keV. Improved sets of energy-loss functions were used for 19 solids. Comparisons were made of these SPs with SPs calculated with the single-pole approximation, previous SP calculations, and experimental SPs. Generally satisfactory agreement was found with SPs from the single-pole approximation for energies above 100 eV, with other calculated SPs, and with measured SPs.
Hadron scattering and resonances in QCD
NASA Astrophysics Data System (ADS)
Dudek, Jozef J.
2016-05-01
I describe how hadron-hadron scattering amplitudes are related to the eigenstates of QCD in a finite cubic volume. The discrete spectrum of such eigenstates can be determined from correlation functions computed using lattice QCD, and the corresponding scattering amplitudes extracted. I review results from the Hadron Spectrum Collaboration who have used these finite volume methods to study ππ elastic scattering, including the ρ resonance, as well as coupled-channel π >K, ηK scattering. Ongoing calculations are advertised and the outlook for finite volume approaches is presented.
QCD resummation for hadronic final states
NASA Astrophysics Data System (ADS)
Luisoni, Gionata; Marzani, Simone
2015-10-01
We review the basic concepts of all-order calculations in quantum chromodynamics (QCD) and their application to collider phenomenology. We start by discussing the factorization properties of QCD amplitudes and cross-sections in the soft and collinear limits and their resulting all-order exponentiation. We then discuss several applications of this formalism to observables which are of great interest at particle colliders. In this context, we describe the all-order resummation of event-shape distributions, as well as observables that probe the internal structure of hadronic jets.
NASA Astrophysics Data System (ADS)
Deur, Alexandre; Brodsky, Stanley J.; de Téramond, Guy F.
2016-09-01
We review the present theoretical and empirical knowledge for αs, the fundamental coupling underlying the interactions of quarks and gluons in Quantum Chromodynamics (QCD). The dependence of αs(Q2) on momentum transfer Q encodes the underlying dynamics of hadron physics-from color confinement in the infrared domain to asymptotic freedom at short distances. We review constraints on αs(Q2) at high Q2, as predicted by perturbative QCD, and its analytic behavior at small Q2, based on models of nonperturbative dynamics. In the introductory part of this review, we explain the phenomenological meaning of the coupling, the reason for its running, and the challenges facing a complete understanding of its analytic behavior in the infrared domain. In the second, more technical, part of the review, we discuss the behavior of αs(Q2) in the high momentum transfer domain of QCD. We review how αs is defined, including its renormalization scheme dependence, the definition of its renormalization scale, the utility of effective charges, as well as "Commensurate Scale Relations" which connect the various definitions of the QCD coupling without renormalization-scale ambiguity. We also report recent significant measurements and advanced theoretical analyses which have led to precise QCD predictions at high energy. As an example of an important optimization procedure, we discuss the "Principle of Maximum Conformality", which enhances QCD's predictive power by removing the dependence of the predictions for physical observables on the choice of theoretical conventions such as the renormalization scheme. In the last part of the review, we discuss the challenge of understanding the analytic behavior αs(Q2) in the low momentum transfer domain. We survey various theoretical models for the nonperturbative strongly coupled regime, such as the light-front holographic approach to QCD. This new framework predicts the form of the quark-confinement potential underlying hadron spectroscopy and
QCD sign problem for small chemical potential
Splittorff, K.; Verbaarschot, J. J. M.
2007-06-01
The expectation value of the complex phase factor of the fermion determinant is computed in the microscopic domain of QCD at nonzero chemical potential. We find that the average phase factor is nonvanishing below a critical value of the chemical potential equal to half the pion mass and vanishes exponentially in the volume for larger values of the chemical potential. This holds for QCD with dynamical quarks as well as for quenched and phase quenched QCD. The average phase factor has an essential singularity for zero chemical potential and cannot be obtained by analytic continuation from imaginary chemical potential or by means of a Taylor expansion. The leading order correction in the p-expansion of the chiral Lagrangian is calculated as well.
Mathematica and Fortran programs for various analytic QCD couplings
NASA Astrophysics Data System (ADS)
Ayala, César; Cvetič, Gorazd
2015-05-01
We outline here the motivation for the existence of analytic QCD models, i.e., QCD frameworks in which the running coupling A(Q2) has no Landau singularities. The analytic (holomorphic) coupling A(Q2) is the analog of the underlying pQCD coupling a(Q2) = αs(Q2)/π, and any such A(Q2) defines an analytic QCD model. We present the general construction procedure for the couplings Av (Q2) which are analytic analogs of the powers a(Q2)v. Three analytic QCD models are presented. Applications of our program (in Mathematica) for calculation of Av (Q2) in such models are presented. Programs in both Mathematica and Fortran can be downloaded from the web page: gcvetic.usm.cl.
Sharada, Shaama Mallikarjun; Bell, Alexis T. E-mail: bell@cchem.berkeley.edu; Head-Gordon, Martin E-mail: bell@cchem.berkeley.edu
2014-04-28
The cost of calculating nuclear hessians, either analytically or by finite difference methods, during the course of quantum chemical analyses can be prohibitive for systems containing hundreds of atoms. In many applications, though, only a few eigenvalues and eigenvectors, and not the full hessian, are required. For instance, the lowest one or two eigenvalues of the full hessian are sufficient to characterize a stationary point as a minimum or a transition state (TS), respectively. We describe here a method that can eliminate the need for hessian calculations for both the characterization of stationary points as well as searches for saddle points. A finite differences implementation of the Davidson method that uses only first derivatives of the energy to calculate the lowest eigenvalues and eigenvectors of the hessian is discussed. This method can be implemented in conjunction with geometry optimization methods such as partitioned-rational function optimization (P-RFO) to characterize stationary points on the potential energy surface. With equal ease, it can be combined with interpolation methods that determine TS guess structures, such as the freezing string method, to generate approximate hessian matrices in lieu of full hessians as input to P-RFO for TS optimization. This approach is shown to achieve significant cost savings relative to exact hessian calculation when applied to both stationary point characterization as well as TS optimization. The basic reason is that the present approach scales one power of system size lower since the rate of convergence is approximately independent of the size of the system. Therefore, the finite-difference Davidson method is a viable alternative to full hessian calculation for stationary point characterization and TS search particularly when analytical hessians are not available or require substantial computational effort.
Nonperturbative QCD corrections to electroweak observables
Dru B Renner, Xu Feng, Karl Jansen, Marcus Petschlies
2011-12-01
Nonperturbative QCD corrections are important to many low-energy electroweak observables, for example the muon magnetic moment. However, hadronic corrections also play a significant role at much higher energies due to their impact on the running of standard model parameters, such as the electromagnetic coupling. Currently, these hadronic contributions are accounted for by a combination of experimental measurements and phenomenological modeling but ideally should be calculated from first principles. Recent developments indicate that many of the most important hadronic corrections may be feasibly calculated using lattice QCD methods. To illustrate this, we will examine the lattice computation of the leading-order QCD corrections to the muon magnetic moment, paying particular attention to a recently developed method but also reviewing the results from other calculations. We will then continue with several examples that demonstrate the potential impact of the new approach: the leading-order corrections to the electron and tau magnetic moments, the running of the electromagnetic coupling, and a class of the next-to-leading-order corrections for the muon magnetic moment. Along the way, we will mention applications to the Adler function, the determination of the strong coupling constant and QCD corrections to muonic-hydrogen.
NLO QCD corrections to Zbb production with massive bottom quarks at the Fermilab Tevatron
Febres Cordero, F.; Reina, L.; Wackeroth, D.
2008-10-01
We calculate the next-to-leading order (NLO) QCD corrections to Zbb production in hadronic collisions including full bottom-quark mass effects. We present results for the total cross section and the invariant mass distribution of the bottom-quark jet pair at the Fermilab Tevatron pp collider. We perform a detailed comparison with a calculation that considers massless bottom quarks, as implemented in the Monte Carlo program MCFM. We find that neglecting bottom-quark mass effects overestimates the total NLO QCD cross section for Zbb production at the Tevatron by about 7%, independent of the choice of the renormalization and factorization scales. Moreover, bottom-quark mass effects can impact the shape of the bottom-quark pair invariant mass distribution, in particular, in the low invariant mass region.
Brodsky, Stanley J.; de Teramond, Guy F.; /Costa Rica U.
2012-02-16
The relation between the hadronic short-distance constituent quark and gluon particle limit and the long-range confining domain is yet one of the most challenging aspects of particle physics due to the strong coupling nature of Quantum Chromodynamics, the fundamental theory of the strong interactions. The central question is how one can compute hadronic properties from first principles; i.e., directly from the QCD Lagrangian. The most successful theoretical approach thus far has been to quantize QCD on discrete lattices in Euclidean space-time. Lattice numerical results follow from computation of frame-dependent moments of distributions in Euclidean space and dynamical observables in Minkowski spacetime, such as the time-like hadronic form factors, are not amenable to Euclidean lattice computations. The Dyson-Schwinger methods have led to many important insights, such as the infrared fixed point behavior of the strong coupling constant, but in practice, the analyses are limited to ladder approximation in Landau gauge. Baryon spectroscopy and the excitation dynamics of nucleon resonances encoded in the nucleon transition form factors can provide fundamental insight into the strong-coupling dynamics of QCD. New theoretical tools are thus of primary interest for the interpretation of the results expected at the new mass scale and kinematic regions accessible to the JLab 12 GeV Upgrade Project. The AdS/CFT correspondence between gravity or string theory on a higher-dimensional anti-de Sitter (AdS) space and conformal field theories in physical space-time has led to a semiclassical approximation for strongly-coupled QCD, which provides physical insights into its nonperturbative dynamics. The correspondence is holographic in the sense that it determines a duality between theories in different number of space-time dimensions. This geometric approach leads in fact to a simple analytical and phenomenologically compelling nonperturbative approximation to the full light
Generalized Parton Distributions from Lattice QCD
Orginos, Konstantinos
2007-10-01
I review recent results on moments of Generalized Parton Distribution functions (GPDs) from Lattice QCD. In particular, I discuss the methodology of lattice calculations, and how various systematic errors arising in these calculations are controlled. I conclude with an overview of the roadmap towards precision non-perturbative determination of moments of GPDs, and discuss the potential impact to the extraction of GPDs form experiment.
Harris, R.
1992-05-01
We present measurements of jet production and isolated prompt photon production in p{bar p} collisions at {radical}s = 1.8 TeV from the 1988--89 run of the Collider Detector at Fermilab (CDF). To test QCD with jets, the inclusive jet cross section (p{bar p} {yields} J + X) and two jet angular distributions (p{bar P} {yields} JJ + X) are compared to QCD predictions and are used to search for composite quarks. The ratio of the scaled jet cross sections at two Tevatron collision energies ({radical}s= 546 and 1800 GeV) is compared to QCD predictions for X{sub T} scaling violations. Also, we present the first evidence for QCD interference effects (color coherence) in third jet production (p{bar p} {yields} JJJ + X). To test QCD with photons, we present measurements of the transverse momentum spectrum of single isolated prompt photon production (p{bar p} {yields} {gamma} + X), double isolated prompt photon production (p{bar p} {yields} {gamma}{gamma} + X), and the angular distribution of photon-jet events (p{bar p} {yields} {gamma} J + X). We have also measured the isolated production ratio of {eta} and {pi}{sup 0} mesons (p{bar p} {yields} {eta} + X)/(p{bar p} {yields} {pi}{sup 0} + X) = 1.02 {plus minus} .15(stat) {plus minus} .23(sys).
Harris, R.; The CDF Collaboration
1992-05-01
We present measurements of jet production and isolated prompt photon production in p{bar p} collisions at {radical}s = 1.8 TeV from the 1988--89 run of the Collider Detector at Fermilab (CDF). To test QCD with jets, the inclusive jet cross section (p{bar p} {yields} J + X) and two jet angular distributions (p{bar P} {yields} JJ + X) are compared to QCD predictions and are used to search for composite quarks. The ratio of the scaled jet cross sections at two Tevatron collision energies ({radical}s= 546 and 1800 GeV) is compared to QCD predictions for X{sub T} scaling violations. Also, we present the first evidence for QCD interference effects (color coherence) in third jet production (p{bar p} {yields} JJJ + X). To test QCD with photons, we present measurements of the transverse momentum spectrum of single isolated prompt photon production (p{bar p} {yields} {gamma} + X), double isolated prompt photon production (p{bar p} {yields} {gamma}{gamma} + X), and the angular distribution of photon-jet events (p{bar p} {yields} {gamma} J + X). We have also measured the isolated production ratio of {eta} and {pi}{sup 0} mesons (p{bar p} {yields} {eta} + X)/(p{bar p} {yields} {pi}{sup 0} + X) = 1.02 {plus_minus} .15(stat) {plus_minus} .23(sys).
Toward lattice QCD simulation on AP1000
NASA Astrophysics Data System (ADS)
Ohta, Shigemi
AP1000 is Fujitsu Laboratory's experimental parallel computer consisting of up to 1024 microcomputers called cells. It is found that each AP1000 cell can sustain two to three megaflops computational speed for full QCD lattice numerical simulations in IEEE 64-bit floating point format
Exploring Hyperons and Hypernuclei with Lattice QCD
S.R. Beane; P.F. Bedaque; A. Parreno; M.J. Savage
2005-01-01
In this work we outline a program for lattice QCD that would provide a first step toward understanding the strong and weak interactions of strange baryons. The study of hypernuclear physics has provided a significant amount of information regarding the structure and weak decays of light nuclei containing one or two Lambda's, and Sigma's. From a theoretical standpoint, little is known about the hyperon-nucleon interaction, which is required input for systematic calculations of hypernuclear structure. Furthermore, the long-standing discrepancies in the P-wave amplitudes for nonleptonic hyperon decays remain to be understood, and their resolution is central to a better understanding of the weak decays of hypernuclei. We present a framework that utilizes Luscher's finite-volume techniques in lattice QCD to extract the scattering length and effective range for Lambda-N scattering in both QCD and partially-quenched QCD. The effective theory describing the nonleptonic decays of hyperons using isospin symmetry alone, appropriate for lattice calculations, is constructed.
Blazey, G.C.
1995-05-01
Selected recent Quantum Chromodynamics (QCD) results from the D0 and CDF experiments at the Fermilab Tevatron are presented and discussed. The inclusive jet and inclusive triple differential dijet cross sections are compared to next-to-leading order QCD calculations. The sensitivity of the dijet cross section to parton distribution functions (for hadron momentum fractions {approximately} 0.01 to {approximately} 0.4) will constrain the gluon distribution of the proton. Two analyses of dijet production at large rapidity separation are presented. The first analysis tests the contributions of higher order processes to dijet production and can be considered a test of BFKL or GLAP parton evolution. The second analysis yields a strong rapidity gap signal consistent with colorless exchange between the scattered partons. The prompt photon inclusive cross section is consistent with next-to-leading order QCD only at the highest transverse momenta. The discrepancy at lower momenta may be indicative of higher order processes impacting a transverse momentum or ``k{sub T}`` to the partonic interaction. The first measurement of the strong coupling constant from the Tevatron is also presented. The coupling constant can be determined from the ratio of W + 1jet to W + 0jet cross sections and a next-to-leading order QCD calculation.
Big bang nucleosynthesis and ΛQCD
NASA Astrophysics Data System (ADS)
Kneller, James P.; McLaughlin, Gail C.
2003-11-01
Big bang nucleosynthesis (BBN) has increasingly become the tool of choice for investigating the permitted variation of fundamental constants during the earliest epochs of the Universe. Here we present a BBN calculation that has been modified to permit changes in the QCD scale, ΛQCD. The primary effects of changing the QCD scale upon BBN are through the deuteron binding energy BD and the neutron-proton mass difference δmnp, which both play crucial roles in determining the primordial abundances. In this paper we show how a simplified BBN calculation allows us to restrict the nuclear data we need to just BD and δmnp yet still gives useful results so that any variation in ΛQCD may be constrained via the corresponding shifts in BD and δmnp by using the current estimates of the primordial deuterium abundance and helium mass fraction. The simplification predicts the helium-4 and deuterium abundances to within 1% and 50%, respectively, when compared with the results of a standard BBN code. But ΛQCD also affects much of the remaining required nuclear input so this method introduces a systematic error into the calculation and we find a degeneracy between BD and δmnp. We show how increased understanding of the relationship of the pion mass and/or BD to other nuclear parameters, such as the binding energy of tritium and the cross section of T+D→4He+n, would yield constraints upon any change in BD and δmnp at the 10% level.
NASA Astrophysics Data System (ADS)
These are the proceedings of the QCD Evolution 2015 Workshop which was held 26-30 May, 2015 at Jefferson Lab, Newport News, Virginia, USA. The workshop is a continuation of a series of workshops held during four consecutive years 2011, 2012, 2013 at Jefferson Lab, and in 2014 in Santa Fe, NM. With the rapid developments in our understanding of the evolution of parton distributions including low-x, TMDs, GPDs, higher-twist correlation functions, and the associated progress in perturbative QCD, lattice QCD and effective field theory techniques we look forward with great enthusiasm to the 2015 meeting. A special attention was also paid to participation of experimentalists as the topics discussed are of immediate importance for the JLab 12 experimental program and a future Electron Ion Collider.
Brodsky, Stanley J.; /SLAC
2007-07-06
I discuss a number of novel topics in QCD, including the use of the AdS/CFT correspondence between Anti-de Sitter space and conformal gauge theories to obtain an analytically tractable approximation to QCD in the regime where the QCD coupling is large and constant. In particular, there is an exact correspondence between the fifth-dimension coordinate z of AdS space and a specific impact variable {zeta} which measures the separation of the quark constituents within the hadron in ordinary space-time. This connection allows one to compute the analytic form of the frame-independent light-front wavefunctions of mesons and baryons, the fundamental entities which encode hadron properties and allow the computation of exclusive scattering amplitudes. I also discuss a number of novel phenomenological features of QCD. Initial- and final-state interactions from gluon-exchange, normally neglected in the parton model, have a profound effect in QCD hard-scattering reactions, leading to leading-twist single-spin asymmetries, diffractive deep inelastic scattering, diffractive hard hadronic reactions, the breakdown of the Lam Tung relation in Drell-Yan reactions, and nuclear shadowing and non-universal antishadowing--leading-twist physics not incorporated in the light-front wavefunctions of the target computed in isolation. I also discuss tests of hidden color in nuclear wavefunctions, the use of diffraction to materialize the Fock states of a hadronic projectile and test QCD color transparency, and anomalous heavy quark effects. The presence of direct higher-twist processes where a proton is produced in the hard subprocess can explain the large proton-to-pion ratio seen in high centrality heavy ion collisions.
Early Run 2 Hard QCD Results from the ATLAS Collaboration
NASA Astrophysics Data System (ADS)
Orlando, Nicola
2016-07-01
We provide an overview of hard QCD results based on data collected with the ATLAS detector in proton-proton collision at √s = 13 TeV at the Large Hadron Collider. The production of high transverse momentum jets, photons and photon-pairs were studied; the inclusive jet cross section is found to agree well with the prediction of perturbative QCD calculations performed at next-to-leading accuracy. The production cross sections for W and Z bosons in their e and μ decays was measured; in general, agreement is found with the expectation of next-to-next-to leading order QCD calculations and interesting sensitivities to the proton structure functions are already observed. The top production cross sections were measured in different top decay channels and found to agree with the state of the art QCD predictions.
Fister, T. T.; Nagel, K. P.; Vila, F. D.; Seidler, G. T.; Hamner, C.; Cross, J. O.; Rehr, J. J.; Univ. Washington; Washington State Univ.
2009-01-01
We report measurements of the nonresonant inelastic x-ray scattering (NRIXS) from the O 1s orbitals in ice Ih, and also report calculations of the corresponding spectra for ice Ih and several other phases of water ice. We find that the intermediate-energy fine structure may be calculated well using an ab initio real-space full multiple scattering approach and that it provides a strong fingerprint of the intermediate-range order for some ice phases. Both experiment and theory find that the intermediate-range fine structure, unlike the near-edge structure, is independent of momentum transfer (q) to very high q. These results have important consequences for future NRIXS measurements of high-pressure phases of ice.
Second-order QCD corrections to jet production at hadron colliders: the all-gluon contribution.
Gehrmann-De Ridder, A; Gehrmann, T; Glover, E W N; Pires, J
2013-04-19
We report the calculation of next-to-next-to-leading order QCD corrections in the purely gluonic channel to dijet production and related observables at hadron colliders. Our result represents the first next-to-next-to-leading order calculation of a massless jet observable at hadron colliders, and opens the path towards precision QCD phenomenology with the LHC. PMID:23679596
NASA Astrophysics Data System (ADS)
Hluchá, H.; Eberl, H.; Frisch, W.
2012-10-01
SFOLD (Sfermion Full One-Loop Decays) is a Fortran program package for calculating all sfermion two-body decay widths and the corresponding branching ratios at full one-loop level within the MSSM. The package adopts the SUSY Parameter Analysis convention and supports the SUSY Les Houches Accord input and output format. With the SFOLD package we found non-negligible electroweak corrections in bosonic decays of b˜,t˜ and τ˜. Program summaryProgram title: SFOLD Catalogue identifier: AEMZ_v1_0 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEMZ_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: GNU General Public License, version 3 No. of lines in distributed program, including test data, etc.: 263346 No. of bytes in distributed program, including test data, etc.: 1481697 Distribution format: tar.gz Programming language: Fortran 77. Computer: Workstation, PC. Operating system: Linux. RAM: approx. 500 Mbytes Classification: 11.1. External routines: LoopTools 2.6 (http://www.feynarts.de/looptools/), SLHALib 2.2 (http://www.feynarts.de/slha/) Nature of problem: If the MSSM is realized in nature, LHC will produce supersymmetric particles copiously. The best environment for a precise determination of the model parameters would be a high energy e+e- linear collider. Experimental accuracies are expected at the per-cent down to the per-mill level. These must be matched from the theoretical side. Therefore loop calculations are mandatory. Solution method: This program package calculates all sfermion two-body decay widths and the corresponding branching ratios at full one-loop level within the MSSM. The renormalization is done in the DR¯ scheme following the SUSY Parameter Analysis convention. The program supports the SUSY Les Houches Accord input and output format. Running time: The test provided just takes a few seconds to run.
Brodsky, Stanley J.; de Teramond, Guy F.; /SLAC /Southern Denmark U., CP3-Origins /Costa Rica U.
2011-01-10
AdS/QCD, the correspondence between theories in a dilaton-modified five-dimensional anti-de Sitter space and confining field theories in physical space-time, provides a remarkable semiclassical model for hadron physics. Light-front holography allows hadronic amplitudes in the AdS fifth dimension to be mapped to frame-independent light-front wavefunctions of hadrons in physical space-time. The result is a single-variable light-front Schroedinger equation which determines the eigenspectrum and the light-front wavefunctions of hadrons for general spin and orbital angular momentum. The coordinate z in AdS space is uniquely identified with a Lorentz-invariant coordinate {zeta} which measures the separation of the constituents within a hadron at equal light-front time and determines the off-shell dynamics of the bound state wavefunctions as a function of the invariant mass of the constituents. The hadron eigenstates generally have components with different orbital angular momentum; e.g., the proton eigenstate in AdS/QCD with massless quarks has L = 0 and L = 1 light-front Fock components with equal probability. Higher Fock states with extra quark-anti quark pairs also arise. The soft-wall model also predicts the form of the nonperturbative effective coupling and its {beta}-function. The AdS/QCD model can be systematically improved by using its complete orthonormal solutions to diagonalize the full QCD light-front Hamiltonian or by applying the Lippmann-Schwinger method to systematically include QCD interaction terms. Some novel features of QCD are discussed, including the consequences of confinement for quark and gluon condensates. A method for computing the hadronization of quark and gluon jets at the amplitude level is outlined.
Vector Meson Form Factors and Wave Functions from Holographic QCD
Hovhannes Grigoryan; Anatoly Radyushkin
2007-10-10
Based on the holographic dual model of QCD, we study 2- and 3-point functions of vector currents and derive form factors as well as wave functions for the vector mesons. As a result, generalized vector-meson dominance representation for form factors is obtained with a very specific VMD pattern. The calculated electric radius of the rho-meson is shown to be in a good agreement with predictions from lattice QCD.
Highly excited and exotic meson spectroscopy from lattice QCD
Christopher Thomas
2011-05-01
I will discuss recent progress in extracting highly excited and exotic meson spectra using lattice QCD. New results in the light meson sector will be presented, where a combination of techniques have enabled us to confidently identify the spin of extracted states. Highlights include many states with exotic quantum numbers and, for the first time in a lattice QCD calculation, spin-four states. I will conclude with comments on future prospects.
Charm and bottom hadronic form factors with QCD sum rules
Bracco, M. E.; Rodrigues, B. O.; Cerqueira, A. Jr.
2013-03-25
We present a brief review of some calculations of form factors and coupling constants in vertices with charm and bottom mesons in the framework of QCD sum rules. We first discuss the motivation for this work, describing possible applications of these form factors to charm and bottom decays processes. We first make a summarize of the QCD sum rules method. We give special attention to the uncertainties of the method introducing by the intrinsic variation of the parameters. Finally we conclude.
NASA Astrophysics Data System (ADS)
Zhao, Zhiqiang; Chen, Jun; Zhang, Zhaojun; Zhang, Dong H.; Lauvergnat, David; Gatti, Fabien
2016-05-01
Full quantum mechanical calculations of vibrational energies of methane and fluoromethane are carried out using a polyspherical description combining Radau and Jacobi coordinates. The Hamiltonian is built in a potential-optimized discrete variable representation, and vibrational energies are solved using an iterative eigensolver. This new approach can be applied to a large variety of molecules. In particular, we show that it is able to accurately and efficiently compute eigenstates for four different molecules : CH4, CHD3, CH2D2, and CH3F. Very good agreement is obtained with the results reported previously in the literature with different approaches and with experimental data.
Zhao, Zhiqiang; Chen, Jun; Zhang, Zhaojun; Zhang, Dong H; Lauvergnat, David; Gatti, Fabien
2016-05-28
Full quantum mechanical calculations of vibrational energies of methane and fluoromethane are carried out using a polyspherical description combining Radau and Jacobi coordinates. The Hamiltonian is built in a potential-optimized discrete variable representation, and vibrational energies are solved using an iterative eigensolver. This new approach can be applied to a large variety of molecules. In particular, we show that it is able to accurately and efficiently compute eigenstates for four different molecules : CH4, CHD3, CH2D2, and CH3F. Very good agreement is obtained with the results reported previously in the literature with different approaches and with experimental data. PMID:27250301
B ---> pi and B ---> K transitions from QCD sum rules on the light cone
Ball, P.
1998-09-01
I calculate the form factors describing semileptonic and penguin-induced decays of B mesons into light pseudoscalar mesons. The form factors are calculated from QCD sum rules on the light-cone including contributions up to twist 4, radiative corrections to the leading twist contribution and SU(3)-breaking effects. The theoretical uncertainty is estimated to be \\sim 15%. The heavy-quark-limit relations between semileptonic and penguin form factors are found to be valid in the full accessible range of momentum transfer.
Devlin, T.; CDF Collaboration
1996-10-01
The CDF collaboration is engaged in a broad program of QCD measurements at the Fermilab Tevatron Collider. I will discuss inclusive jet production at center-of-mass energies of 1800 GeV and 630 GeV, properties of events with very high total transverse energy and dijet angular distributions.
Radyushkin, Anatoly V.; Efremov, Anatoly Vasilievich; Ginzburg, Ilya F.
2013-04-01
We discuss some problems concerning the application of perturbative QCD to high energy soft processes. We show that summing the contributions of the lowest twist operators for non-singlet $t$-channel leads to a Regge-like amplitude. Singlet case is also discussed.
Lincoln, Don
2016-06-28
The strongest force in the universe is the strong nuclear force and it governs the behavior of quarks and gluons inside protons and neutrons. The name of the theory that governs this force is quantum chromodynamics, or QCD. In this video, Fermilab?s Dr. Don Lincoln explains the intricacies of this dominant component of the Standard Model.
Nathan Isgur
1997-03-01
The author presents an idiosyncratic view of baryons which calls for a marriage between quark-based and hadronic models of QCD. He advocates a treatment based on valence quark plus glue dominance of hadron structure, with the sea of q pairs (in the form of virtual hadron pairs) as important corrections.
THE TOP QUARK, QCD, AND NEW PHYSICS.
DAWSON,S.
2002-06-01
The role of the top quark in completing the Standard Model quark sector is reviewed, along with a discussion of production, decay, and theoretical restrictions on the top quark properties. Particular attention is paid to the top quark as a laboratory for perturbative QCD. As examples of the relevance of QCD corrections in the top quark sector, the calculation of e{sup +}e{sup -} + t{bar t} at next-to-leading-order QCD using the phase space slicing algorithm and the implications of a precision measurement of the top quark mass are discussed in detail. The associated production of a t{bar t} pair and a Higgs boson in either e{sup +}e{sup -} or hadronic collisions is presented at next-to-leading-order QCD and its importance for a measurement of the top quark Yulrawa coupling emphasized. Implications of the heavy top quark mass for model builders are briefly examined, with the minimal supersymmetric Standard Model and topcolor discussed as specific examples.
The Top Quark, QCD, And New Physics.
DOE R&D Accomplishments Database
Dawson, S.
2002-06-01
The role of the top quark in completing the Standard Model quark sector is reviewed, along with a discussion of production, decay, and theoretical restrictions on the top quark properties. Particular attention is paid to the top quark as a laboratory for perturbative QCD. As examples of the relevance of QCD corrections in the top quark sector, the calculation of e{sup+}e{sup -}+ t{bar t} at next-to-leading-order QCD using the phase space slicing algorithm and the implications of a precision measurement of the top quark mass are discussed in detail. The associated production of a t{bar t} pair and a Higgs boson in either e{sup+}e{sup -} or hadronic collisions is presented at next-to-leading-order QCD and its importance for a measurement of the top quark Yulrawa coupling emphasized. Implications of the heavy top quark mass for model builders are briefly examined, with the minimal supersymmetric Standard Model and topcolor discussed as specific examples.
Brodsky, Stanley J.; /SLAC /Southern Denmark U., CP3-Origins
2011-08-12
I review a number of topics where conventional wisdom in hadron physics has been challenged. For example, hadrons can be produced at large transverse momentum directly within a hard higher-twist QCD subprocess, rather than from jet fragmentation. Such 'direct' processes can explain the deviations from perturbative QCD predictions in measurements of inclusive hadron cross sections at fixed x{sub T} = 2p{sub T}/{radical}s, as well as the 'baryon anomaly', the anomalously large proton-to-pion ratio seen in high centrality heavy ion collisions. Initial-state and final-state interactions of the struck quark, the soft-gluon rescattering associated with its Wilson line, lead to Bjorken-scaling single-spin asymmetries, diffractive deep inelastic scattering, the breakdown of the Lam-Tung relation in Drell-Yan reactions, as well as nuclear shadowing and antishadowing. The Gribov-Glauber theory predicts that antishadowing of nuclear structure functions is not universal, but instead depends on the flavor quantum numbers of each quark and antiquark, thus explaining the anomalous nuclear dependence measured in deep-inelastic neutrino scattering. Since shadowing and antishadowing arise from the physics of leading-twist diffractive deep inelastic scattering, one cannot attribute such phenomena to the structure of the nucleus itself. It is thus important to distinguish 'static' structure functions, the probability distributions computed from the square of the target light-front wavefunctions, versus 'dynamical' structure functions which include the effects of the final-state rescattering of the struck quark. The importance of the J = 0 photon-quark QCD contact interaction in deeply virtual Compton scattering is also emphasized. The scheme-independent BLM method for setting the renormalization scale is discussed. Eliminating the renormalization scale ambiguity greatly improves the precision of QCD predictions and increases the sensitivity of searches for new physics at the LHC
Wang, Shuai; Luo, Xuenong; Wei, Wei; Zheng, Yadong; Dou, Yongxi; Cai, Xuepeng
2013-01-01
Individual genes or regions are still commonly used to estimate the phylogenetic relationships among viral isolates. The genomic regions that can faithfully provide assessments consistent with those predicted with full-length genome sequences would be preferable to serve as good candidates of the phylogenetic markers for molecular epidemiological studies of many viruses. Here we employed a statistical method to evaluate the evolutionary relationships between individual viral genes and full-length genomes without tree construction as a way to determine which gene can match the genome well in phylogenetic analyses. This method was performed by calculation of linear correlations between the genetic distance matrices of aligned individual gene sequences and aligned genome sequences. We applied this method to the phylogenetic analyses of porcine circovirus 2 (PCV2), measles virus (MV), hepatitis E virus (HEV) and Japanese encephalitis virus (JEV). Phylogenetic trees were constructed for comparisons and the possible factors affecting the method accuracy were also discussed in the calculations. The results revealed that this method could produce results consistent with those of previous studies about the proper consensus sequences that could be successfully used as phylogenetic markers. And our results also suggested that these evolutionary correlations could provide useful information for identifying genes that could be used effectively to infer the genetic relationships. PMID:24312527
Exploring Flavor Physics with Lattice QCD
NASA Astrophysics Data System (ADS)
Du, Daping; Fermilab/MILC Collaborations Collaboration
2016-03-01
The Standard Model has been a very good description of the subatomic particle physics. In the search for physics beyond the Standard Model in the context of flavor physics, it is important to sharpen our probes using some gold-plated processes (such as B rare decays), which requires the knowledge of the input parameters, such as the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements and other nonperturbative quantities, with sufficient precision. Lattice QCD is so far the only first-principle method which could compute these quantities with competitive and systematically improvable precision using the state of the art simulation techniques. I will discuss the recent progress of lattice QCD calculations on some of these nonpurturbative quantities and their applications in flavor physics. I will also discuss the implications and future perspectives of these calculations in flavor physics.
Exploring Three Nucleon Forces in Lattice QCD
Doi, Takumi
2011-10-21
We study the three nucleon force in N{sub f} = 2 dynamical clover fermion lattice QCD, utilizing the Nambu-Bethe-Salpeter wave function of the three nucleon system. Since parity-odd two nucleon potentials are not available in lattice QCD at this moment, we develop a new formulation to extract the genuine three nucleon force which requires only the information of parity-even two nucleon potentials. In order to handle the extremely expensive calculation cost, we consider a specific three-dimensional coordinate configuration for the three nucleons. We find that the linear setup is advantageous, where nucleons are aligned linearly with equal spacings. The lattice calculation is performed with 16{sup 3}x32 configurations at {beta} = 1.95, m{sub {pi}} = 1.13 GeV generated by CP-PACS Collaboration, and the result of the three nucleon force in triton channel is presented.
Lattice QCD Beyond Ground States
Huey-Wen Lin; Saul D. Cohen
2007-09-11
In this work, we apply black box methods (methods not requiring input) to find excited-state energies. A variety of such methods for lattice QCD were introduced at the 3rd iteration of the numerical workshop series. We first review a selection of approaches that have been used in lattice calculations to determine multiple energy states: multiple correlator fits, the variational method and Bayesian fitting. In the second half, we will focus on a black box method, the multi-effective mass. We demonstrate the approach on a toy model, as well as on real lattice data, extracting multiple states from single correlators. Without complicated operator construction or specialized fitting programs, the black box method shows good consistency with the traditional approaches.
Radiative Charmonium Physics from Lattice QCD
Jozef Dudek; Robert Edwards; David Richards; Nilmani Mathur
2007-08-09
Charmonium is an attractive system for the application of lattice QCD methods owing to the possibility of computing with the physical quark mass in a reasonable time. While the sub-threshold spectrum has been considered in some detail in previous works, it is only very recently that further properties such as radiative transitions and two-photon decays have come to be calculated; herein we discuss this recent progress.