“…2,[10][11][12][13][14] For the MnPc, the multiplet for ground state has been suggested either 4 E g or 4 B 2g state, 11,15,16 and for the CoPc, the 2 E g or 2 A 1g state. 11,13 For this, we previously proposed a constraint DFT approach, implemented in the full-potential linearized augmented plane-wave (FLAPW) method, 17 by imposing density matrix constraints on the d-orbitals occupation numbers, and our total energy calculations demonstrate that for the FePc, there are three stationary states of the 3 E g , 3 B 2g , and 3 A 2g symmetries and that the ground state is the 3 A 2g state with a planar magnetic anisotropy (MA). Here, we extend our DFT investigations for the MnPc and CoPc, and we demonstrate that for the MnPc (CoPc), the ground state is the 4 E g ( 2 A 1g ) state with the out-of-plane (planar) MA.…”