Automated Generation of DKH Property Operators

The evaluation of DKH property operators up to any predefined order in the external potential requires a computational scheme, which cannot rely on a direct noniterative numerical method since odd operators have no representation in a truly two-component formalism (which includes the one-component approximation). Only the product of two odd operators may be evaluated. 

Each order of the upper left block of the transformed property operator 

---