<div dir="ltr"><div>Dear CP2K users,</div><div><br></div><div>I am trying to calculate the dipole moment and polarizability for <b>each OH bond</b> in a bulk water system (32 water in a cubic box with PBC).  The input file (H2O32_md.inp; pos.inc) for DFT MD simulation is attached here. </div><div><br></div><div>Here is my method to calculate the dipole moment for each OH bond.</div>I can use the Wannier centers to calculate the dipole moment for each water molecule (using $ \mu = \sum_i Z_i R_i -2 r^w_n$), and therefore the dipole moment of each OH bond in each water molecule.  This method works, but it is not very direct from CP2K code.<div><br></div><div>But for the polarizibility, even though I can obtain the polarizaibility (see PS) of the whole system,  I still do not know the method to extract the polarizability for neither each water molecule or OH bonds in them.</div><div><br></div><div>My question is: With <b>CP2K</b>, is it possbile to directly calculate <b>the dipole moment vector and polarizibility tensor</b> for <b>a OH bond</b> in a water molecule for the bulk water system?<br></div><div><div><div><br></div><div>With best regards,</div><div>jia</div><div><br></div><div><br></div><div><br></div><div>PS: the first 3 lines of the result for the polarizibility for the whole system:</div><div><br></div><div><div>jia@xil /home/jia/H2O32_md [16]$ head H2O32_md-1.polar</div><div>#  Â Step  Â Time [fs]  Â  Â  Â  Â  Â xx [a.u.]  Â  Â  Â  Â  Â xy [a.u.]  Â  Â  Â  Â  Â xz [a.u.]  Â  Â  Â  Â  Â yx [a.u.]  Â  Â  Â  Â  Â yy [a.u.]  Â  Â  Â  Â  Â yz [a.u.]  Â  Â  Â  Â  Â zx [a.u.]  Â  Â  Â  Â  Â zy [a.u.]  Â  Â  Â  Â  Â zz [a.u.]</div><div>  Â  Â  Â 0  Â  Â  Â 0.000  Â  Â  Â  268.31915474  Â  Â  Â  Â  1.84753325  Â  Â  Â  Â -7.77104321  Â  Â  Â  Â  1.80422537  Â  Â  Â  273.71170715  Â  Â  Â  Â  2.58088055  Â  Â  Â  Â -7.78446655  Â  Â  Â  Â  2.56412737  Â  Â  Â  271.68321643</div><div>  Â  Â  Â 1  Â  Â  Â 0.500  Â  Â  Â  268.28862178  Â  Â  Â  Â  1.67253194  Â  Â  Â  Â -7.78792499  Â  Â  Â  Â  1.63612816  Â  Â  Â  274.88797255  Â  Â  Â  Â  3.09590223  Â  Â  Â  Â -7.80366438  Â  Â  Â  Â  3.08627690  Â  Â  Â  273.27103200</div></div><div>  Â  Â  Â 2  Â  Â  Â 1.000  Â  Â  Â  269.95417495  Â  Â  Â  Â  1.47505442  Â  Â  Â  Â -7.87087433  Â  Â  Â  Â  1.44141475  Â  Â  Â  277.68566066  Â  Â  Â  Â  3.65694561  Â  Â  Â  Â -7.89349554  Â  Â  Â  Â  3.63530079  Â  Â  Â  276.41211218<br></div><div><br></div><div><br></div></div></div></div>