Programs in Physics & Physical Chemistry
|[Licence| Download | New Version Template] aeja_v1_0.tar.gz(616 Kbytes)|
|Manuscript Title: MathQCDSR: a Mathematica package for QCD sum rules calculations|
|Authors: Lai Wang, Frank X. Lee|
|Program title: MathQCDSR|
|Catalogue identifier: AEJA_v1_0|
Distribution format: tar.gz
|Journal reference: Comput. Phys. Commun. 182(2011)1721|
|Programming language: Mathematica.|
|Computer: PCs and Workstations.|
|Operating system: Any OS that supports Mathematica.The package has been tested under Windows XP, Macintosh OS X, and Linux.|
|Keywords: Quantum chromodynamics, QCD sum rule method, baryon masses and magnetic moments, Mathematica, FeynCalc.|
|PACS: 12.38.-t, 11.55.Hx, 13.40.Em, 14.20.Gk, 14.20.Jn.|
External routines: FeynCalc (http://www.feyncalc.org/). It is a freely available Mathematica package for high-energy physics calculations. Here it is used primarily to handle gamma-matrix algebra.
Nature of problem:
The QCD sum rule method is a nonperturbative approach to solving quantum chromodynamics (QCD), the fundamental theory of the strong force. The approach establishes a direct link between hadron phenomenology and the QCD vacuum structure via a few QCD parameters called vacuum condensates and susceptibilities. It has been widely applied in nuclear and particle physics to gain insight into various aspects of strong-interaction physics.
First, QCD sum rules are constructed by evaluating correlation functions from two perspectives. On the quark level, it leads to a function of QCD parameters and the Borel mass parameter M. On the hadronic level, it leads to a function of phenomenological parameters and the same M. By numerically matching the two sides over a range in M, the phenomenological parameters can be extracted. The construction involves a large amount of gamma-matrix algebra, Fourier transform, and Borel transform. The matching usually involves searching for minimum χ2. We employ a Monte-Carlo-based procedure to perform the analysis which allows for realistic error estimates
The package deals with only standard (SVZ) QCD sum rule calculations. It can be easily adapted to handle other variants of the method (like finite-energy sum rules). Due to the use of FeynCalc, two of the notebooks (qcdsr2pt-construction.nb and qcdsr3pt-construction.nb) only run on version 6.0 of Mathematica. The other two can run on any version.
The package consists of the following 4 notebooks.
For mass calculations, qcdsr2ptconstruction.nb and qcdsr2pt-analysis.nb take about a minute each to run on a laptop. For magnetic moment calculations, qcdsr3pt-construction.nb can take up to 10 minutes for a given particle, and qcdsr3pt-analysis.nb typically a few minutes, depending on the number of Monte-Carlo samples
|Disclaimer | ScienceDirect | CPC Journal | CPC | QUB|