# Adding a new operator in HubbardSite Class

+1 vote

Hi, I recently have to compute the expectation value of pair-pair correlation function for s wave so I need the operator like C{i,up}C{i,down}. I add two new operators in the class HubbardSite as followings:

else
if(opname == "sWavedag")
{
Op.set(Em,UDP,-1);
}
else
if(opname == "sWave")
{
Op.set(UD,EmP,-1);
}
Can I get the correct IQMPO for the operator C{j,down}^{\dag}C{j,up}^{\dag}C{i,up}C{i,down} by using the AutoMPO like the following piece of the code:
ampo += 1 "sWavedag", j, "sWave",i ;
?