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="6b6eaa1f-30c0-483c-bfc5-2d171ddf49f1"><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="0a513828-4889-41d9-8169-d4e6803b806d"><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="48dbc071-57ab-4ca1-be68-ce6d788e2858"><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="d175eb1f-8364-4305-9559-ed7e2dfa856a"><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="8f496d9b-e375-467a-9d92-f7b037c22c71"><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="dc5008e3-e922-4753-9ccf-e9824e23e121"><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="3bb35fdd-e441-4adc-a9e1-aef6e511f3b4"><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() |