I am confused about how to define a harmonic trap in ALPS, shall I add it directly to the SITETERM as follows?
<SITETERM type="0" site="i">
<PARAMETER name="x" value="i"/>
-mu*n(i)+V*n(i)*(x-L/2)*(x-L/2)
</SITETERM>
where parameters mu and V, operator n(i) are defined in my_models.xml
but, L the length of lattice is defined in a separate independent file my_lattices.xml.
Is this kind of definition of harmonic trap OK in ALPS?
Bo-Nan
--
Stay foolish,Stay hungry.