[julia] how to create local identity operator

Question

Dear all

How to create local identity operator at site j? Indirectly, one can use the identity "S^zjS^zj" or "projUp+projDn", which have built-in functions like op("Z",sites,j) and op("projDn",sites,j). Is there a more standard way?

Moreover, how about high-spin sites, fermion/electron sites?

Best,
Junsen

Answer 1

Hi Junsen,
You can use the operator name "Id" to do this. So like op("Id",sites,j) or alternatively op("Id",ind) where ind is a Index object.

This particular operator ("Id") is actually defined for all indices, not just ones carrying special tags.

Hope that is the answer you need -

Miles

Comments

thank you Miles, that's what I need :)