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="7d47adb3-2d83-4002-afdc-5f2971cd9d64"><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="60218bd0-e7c7-4f35-83b1-4888e8338209"><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="03f20658-7866-4fb6-85d6-05ab5c8d96e2"><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="504cb719-e1bf-41ff-aaff-5c04a688b745"><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="63493a46-c00f-4ec7-bec2-e8f331849452"><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="e48b162f-bbce-4da5-8fe5-db55ad30f232"><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="4ae6056f-7555-419b-b80d-7c7476c6876e"><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() |