Calculate dynamic correlation function for fermonic operators in ITensor

I find that in Proper way to construct fermion gates Mile’s reply mentioned that the apply function will automatically deal with the fermion string so maybe this should be the correct way.