I would like to know if the Constrained DFT is implemented within the plane wave basis in NWCHEM? If not, do you have a plan to implement it in the future? I think it will become attractive for people who want to study a charged system on surface or something like that. Thanks in advance.