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="c55391bf-9194-4d15-b1ae-fc8ce874c4c0"><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="22f7d275-7d40-451f-bff8-68b90b48b3cf"><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="8f72e24c-c6e1-4ab0-9d15-b3a6432e7fb6"><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="2d3d6206-61de-4498-84f8-54e0fe8ea21b"><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="c19681dc-c4f2-44d9-8af6-37e021820951"><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="2238eee7-bef0-4c9c-be35-b90762a4697a"><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="e223fa94-e5b2-40e5-b26d-ae1e5f68b9b4"><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() |