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="f6f001ef-6786-4fe6-a27a-7e8981cbf412"><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="5f1eb70a-4940-4bba-94f8-42e95422b321"><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="9bbd16ab-fa88-4390-a471-941edc13b011"><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="54204a26-4d27-4abc-bc35-a6ce31a59dc5"><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="6296affc-6e11-4761-921d-88d8b2e3258d"><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="0a24bd81-0a3e-4a7a-8887-24f4e98cb554"><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="1b02df74-1967-4fed-837d-7d621a5f1929"><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() |