Towards nuclear reactions from lattice QCD
"In this talk I will motivate the evaluation of nuclear reactions cross sections from Lattice Quantum Chromodynamics (LQCD) and discuss challenges associated with such calculations. In particular, I will explore the connection between the energy spectrum of a three-body system in a finite volume and infinite volume scattering matrix elements using an effective field theoretical approach. The implication of this formalism for studying systems composed of a particle and a bound-state below the bound-state break-up will be discussed. I will briefly discuss the similarities of the three-body problem and that of two-body coupled-channels systems and will comment on extensions of the formalism above the break-up threshold."
#!/usr/common/usg/python/2.7.1/bin/python
import h5py
import os
hh=h5py.File('cl21_32_256.hdf5','w')
flist=os.listdir('GLUONIC')
<ac:structured-macro ac:name="unmigrated-wiki-markup" ac:schema-version="1" ac:macro-id="c9a60f16-584c-48ec-a42d-5029ecc6a8c3"><ac:plain-text-body><![CDATA[clist=[]]]></ac:plain-text-body></ac:structured-macro>
<ac:structured-macro ac:name="unmigrated-wiki-markup" ac:schema-version="1" ac:macro-id="f0a9058b-b7a5-48ae-a38b-3dbeed525aed"><ac:plain-text-body><![CDATA[slist=[]]]></ac:plain-text-body></ac:structured-macro>
for f in flist:
print f
c=f.replace('.hdf5','')
clist.append(c)
ii=h5py.File('GLUONIC/'+f,'r')
ll=list(ii)
ll.remove('info')
<ac:structured-macro ac:name="unmigrated-wiki-markup" ac:schema-version="1" ac:macro-id="010d5a31-f454-4f59-9662-de25b9e97875"><ac:plain-text-body><![CDATA[ ns=ii['info/N_src'][()]]]></ac:plain-text-body></ac:structured-macro>
slist.append(ns)
<ac:structured-macro ac:name="unmigrated-wiki-markup" ac:schema-version="1" ac:macro-id="9c1dc73b-ce78-43c3-8102-38e4948dd451"><ac:plain-text-body><![CDATA[ if f == flist[0]:]]></ac:plain-text-body></ac:structured-macro>
<ac:structured-macro ac:name="unmigrated-wiki-markup" ac:schema-version="1" ac:macro-id="fd86148d-cd8e-405f-9eb6-6ebfce9a43f2"><ac:plain-text-body><![CDATA[ nt=ii['info/Nt'][()]]]></ac:plain-text-body></ac:structured-macro>
<ac:structured-macro ac:name="unmigrated-wiki-markup" ac:schema-version="1" ac:macro-id="86d836f1-528d-42a6-8c5d-e7a3e9e6cb77"><ac:plain-text-body><![CDATA[ nx=ii['info/Nx'][()]]]></ac:plain-text-body></ac:structured-macro>
for l in ll:
hh.create_group(l)
for l in ll:
<ac:structured-macro ac:name="unmigrated-wiki-markup" ac:schema-version="1" ac:macro-id="f808822e-f52c-490f-87f6-2dc1cf2f3b32"><ac:plain-text-body><![CDATA[ hh[l].create_dataset(c,data=ii[l])]]></ac:plain-text-body></ac:structured-macro>
hh.create_dataset('config_list',data=clist)
hh.create_dataset('nsrc_list',data=slist)
hh.close() |