Two-flavor full QCD configurations by CP-PACS

The two-flavor full QCD configurations in JLDG are generated by the CP-PACS Collaboration on the CP-PACS parallel computer at the Center for Computational Sciences, University of Tsukuba.



Simulation parameters

These gauge configurations are generated with a renormalization-group improved gauge action and a mean field improved clover quark action at three values of beta = 6/g2, corresponding to lattice spacings of a \approx 0.22, 0.16 and 0.11 fm, and four sea quark masses corresponding to mPS/mV \approx 0.8, 0.75, 0.7 and 0.6. The sizes of lattice are chosen to be 123x24, 163x32 and 243x48 so that the physical spatial size is kept constant at La \approx 2.5 fm. The scale a is fixed by Mrho = 768.4 MeV.

beta lattice cSW a [fm] La [fm] mPS/mV for sea quarks
#traj (thermalized)
1.80 123x24 1.60 0.2150(22) 2.580(26) 0.807(1)
6250
0.753(1)
5000
0.694(2)
7000
0.547(4)
5250
1.95 163x32 1.53 0.1555(17) 2.489(27) 0.804(1)
7000
0.752(1)
7000
0.690(1)
7000
0.582(3)
5000
2.10 243x48 1.47 0.1076(13) 2.583(31) 0.806(1)
4000
0.755(2)
4000
0.691(3)
4000
0.576(3)
4000
2.20 243x48 1.44 0.0865(33) 2.076(79) 0.799(3)
2000
0.753(4)
2000
0.705(5)
2000
0.632(7)
2000

Configurations are generated with the Hybrid Monte Carlo (HMC) algorithm. At the main coupling constants beta=1.8-2.1, runs are made with a length of 4000-7000 HMC unit-trajectories per sea quark mass. The additional runs at beta=2.2 are stopped at 2000 HMC trajectories per sea quark mass because the physical lattice size La turned out to be significantly smaller than the other three lattices. (Configurations at beta=2.2 are not stored in JLDG) For the inversion of the quark matrix during the HMC update we use the even/odd preconditioned BiCGStab algorithm.

For the clover coefficient cSW we adopt a mean field improved choice defined by cSW = (W1x1)-3/4 = (1-0.8412/beta)-3/4 where for the plaquette W1x1 the value calculated in one-loop perturbation theory is substituted. This choice is based on our observation that the one-loop calculation of W1x1 reproduces the measured values well, with a difference of at most 8% in our simulations. We also note that the one-loop value of cSW = 1+0.678/beta differs from our choice cSW = 1+0.631/beta+... only by a few per cent. The leading scaling violation with our choice of the clover coefficient is O(g2a).

For more details, see


Configuration ensembles

The table below gives a list of ensembles of configurations. A consecutive set of HMC trajectories for a given value of beta, kappa and cSW is called an ensemble. A series name (A, B, C, etc) is given to distinguish multiple HMC runs. The step size dt is sometimes changed in the course of HMC evolution as indicated in the table below.

ensemble # beta kappa cSW series lattice dt initial
traj
final
traj
#config size/config
[MB]
total
[GB]
1 1.80 0.1409 1.60 A 123x24 0.0100
0.0330
210
710
700
6750
47
595
22.815
2 1.80 0.1430 1.60 A 123x24 0.0080
0.0250
210
700
680
5450
47
472
22.812
3 1.80 0.1445 1.60 A 123x24 0.0065
0.0200
0.0167
210
510
580
500
570
3830
30
7
322
22.89
1.80 0.1445 1.60 B 123x24 0.0065 2103880 362 22.89
4 1.80 0.1464 1.60 A 123x24 0.0033 2601740 148 22.84
1.80 0.1464 1.60 B 123x24 0.0033 5901630 104 22.83
1.80 0.1464 1.60 C 123x24 0.0033 6001620 103 22.83
1.80 0.1464 1.60 D 123x24 0.0066 8901420 54 22.82
ensemble # beta kappa cSW series lattice dt initial
traj
final
traj
#config size/config
[MB]
total
[GB]
5 1.95 0.1375 1.53 A 163x32 0.0313 1107100 599 72.044
6 1.95 0.1390 1.53 A 163x32 0.0200
0.0250
110
150
140
7130
4
682
72.050
7 1.95 0.1400 1.53 A 163x32 0.0125
0.0156
0.0185
10
70
420
60
410
3420
6
35
294
72.025
1.95 0.1400 1.53 B 163x32 0.0125
0.0156
0.0185
10
70
480
60
470
3570
6
40
305
72.026
8 1.95 0.1410 1.53 A 163x32 0.0100
0.0080
10
250
240
2500
24
223
72.018
1.95 0.1410 1.53 B 163x32 0.0100
0.0080
10
240
230
2500
22
226
72.018
ensemble # beta kappa cSW series lattice dt initial
traj
final
traj
#config size/config
[MB]
total
[GB]
9 2.10 0.1357 1.47 A 243x48 0.0200 2054205 798 364.5291
10 2.10 0.1367 1.47 A 243x48 0.0160 2154235 788 364.5288
11 2.10 0.1374 1.47 A 243x48 0.0143 2054200 767 364.5280
12 2.10 0.1382 1.47 A 243x48 0.0075 2051505 253 364.593
2.10 0.1382 1.47 B 243x48 0.0075 2051835 320 364.5117
2.10 0.1382 1.47 C 243x48 0.0075 2051265 213 364.578

Naming of Ensembles and Configurations (for ILDG users)

An ensemble is identified by markovChainURI named e.g. as A configuration is labeled by dataLFN named e.g. as Basename of the dataLFN is formmated as
RgClover spatial size x temporal size - \beta\kappa\c_sw - series - trajectory number

To get full list of markovChainURIs and LFNs, please use QCDml faceted navigation or follow this link.

Location of Configurations (for JLDG users)

Configurations are placed in the gfarm directory
  /gfarm/public/ILDG/JLDG/CP-PACS/RCNF2/
and classified in terms of β and κ. For example, the directory
  /gfarm/public/ILDG/JLDG/CP-PACS/RCNF2/RC12x24-B1800/K014090
contains configurations at β=1.800 and κ=0.14090. Filename is the same as basename of dataLFN. All ensemble and configuration XML files can be found in
  /gfarm/public/ILDG/JLDG/CP-PACS/RCNF2/QCDml


Last update: 03 Sept. 2009, T.Yoshie