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="625a0f40-3d56-4cf1-81b1-c28e14ade6fc"><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="7536be4b-006a-43c2-a3ec-d0863b5ac5f8"><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="fec142c1-8cbb-4b61-8b39-701d7de4461d"><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="aa0069e1-b0b6-4145-8718-1b1ad777606c"><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="a0b6fde0-d8da-4b4b-8d87-bd0d4619b0ee"><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="542dc11e-26c4-48eb-bd71-09b571e6d1e8"><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="304485d2-6244-4052-8cdb-703370e0bd61"><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() |