<div dir="ltr">Hi Xiaoming,<div><br></div><div>overlap matrix is already in spherical Gaussians. The overlap is taken between basis functions (AO), not between MO(s). It is calculated before any MOs are obtained.</div><div><br></div><div>Yours,</div><div><br></div><div>Vladimir<br><br>пятница, 12 октября 2018 г., 14:34:43 UTC+2 пользователь Xiaoming Wang написал:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Hi Vladimir,<div><br></div><div>Thanks for your confirmation. Yes, with Cartesian off, the dimensions are consistent now. If cp2k does not print MO of vacant state, then is it possible to transform the overlap matrix to spherical coordinates?</div><div><br></div><div>Best,</div><div>Xiaoming<br><br>On Friday, October 12, 2018 at 8:18:30 AM UTC-4, Vladimir Rybkin wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear Xiaoming,<div><br></div><div>thank you for the data. I don't see any inconsistency now. You get 52 with Cartesian off and 54 with Cartesian on. Turning Cartesian off you can now multiply your C with your S. Isn't that what you want? Unfortunately, I can't comment on why orbital #5 is not printed. Perhaps, that is due to the fact that it is vacant.</div><div><br></div><div>Yours,</div><div><br></div><div>Vladimir<br><br>пятница, 12 октября 2018 г., 13:47:35 UTC+2 пользователь Xiaoming Wang написал:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Vladimir,<div><br></div><div>Thanks for you reply.</div><div><br></div><div>My input of the print part is:</div><div><p style="font-size:11px;line-height:normal;font-family:Menlo"><span>------------</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> &MO</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> EIGENVALUES</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> EIGENVECTORS</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> CARTESIAN</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> FILENAME mos</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> MO_INDEX_RANGE 4 5</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> OCCUPATION_NUMBERS</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> &EACH</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> QS_SCF 0</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> &END EACH</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> &END MO</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span>--------------</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span><br></span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span>With CARTESIAN on, I get</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span><br></span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span>----------------</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> MO EIGENVALUES, MO OCCUPATION NUMBERS, AND CARTESIAN MO EIGENVECTORS</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo;min-height:13px"><span></span><br></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 4 5</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> -0.260963 -0.034173</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo;min-height:13px"><span></span><br></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 2.000000 0.000000</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo;min-height:13px"><span></span><br></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 1 1 O 2s 0.000000 0.236782</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 2 1 O 3s 0.000000 -0.090781</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 3 1 O 4s 0.000000 -0.240258</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 4 1 O 3px -0.917807 -0.000000</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 5 1 O 3py -0.000000 -0.000000</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 6 1 O 3pz 0.000000 -0.238095</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 7 1 O 4px 0.020934 -0.000000</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 8 1 O 4py 0.000000 -0.000000</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 9 1 O 4pz -0.000000 0.013849</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 10 1 O 5px -0.003042 0.000000</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 11 1 O 5py 0.000000 0.000000</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 12 1 O 5pz -0.000000 0.099980</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 13 1 O 3dx2 -0.000000 0.000099</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 14 1 O 3dxy -0.000000 -0.000000</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 15 1 O 3dxz -0.005256 -0.000000</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 16 1 O 3dy2 -0.000000 0.000744</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 17 1 O 3dyz 0.000000 -0.000000</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 18 1 O 3dz2 0.000000 -0.000844</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 19 1 O 4dx2 -0.000000 0.000357</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 20 1 O 4dxy 0.000000 -0.000000</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 21 1 O 4dxz -0.022359 0.000000</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 22 1 O 4dy2 0.000000 0.004105</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 23 1 O 4dyz -0.000000 0.000000</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 24 1 O 4dz2 0.000000 -0.004463</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 25 1 O 5s -0.000000 -0.011663</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 26 1 O 5px 0.024750 -0.000000</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 27 1 O 5py -0.000000 0.000000</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 28 1 O 5pz 0.000000 -0.077018</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo;min-height:13px"><span></span><br></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 29 2 H 1s 0.000000 -0.116793</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 30 2 H 2s 0.000000 -0.003374</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 31 2 H 3s -0.000000 0.197029</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 32 2 H 2px 0.011575 -0.000000</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 33 2 H 2py 0.000000 0.008716</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 34 2 H 2pz -0.000000 -0.004748</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 35 2 H 3px 0.028901 -0.000000</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 36 2 H 3py 0.000000 0.013894</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 37 2 H 3pz -0.000000 -0.010138</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 38 2 H 4s -0.000000 0.421270</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 39 2 H 4px 0.026994 0.000000</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 40 2 H 4py 0.000000 0.041403</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 41 2 H 4pz -0.000000 -0.044975</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo;min-height:13px"><span></span><br></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 42 3 H 1s 0.000000 -0.116793</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 43 3 H 2s -0.000000 -0.003374</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 44 3 H 3s 0.000000 0.197029</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 45 3 H 2px 0.011575 -0.000000</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 46 3 H 2py 0.000000 -0.008716</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 47 3 H 2pz -0.000000 -0.004748</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 48 3 H 3px 0.028901 -0.000000</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 49 3 H 3py -0.000000 -0.013894</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 50 3 H 3pz -0.000000 -0.010138</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 51 3 H 4s 0.000000 0.421270</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 52 3 H 4px 0.026994 0.000000</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 53 3 H 4py 0.000000 -0.041403</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 54 3 H 4pz 0.000000 -0.044975</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo;min-height:13px"><span></span><br></p><p style="font-size:11px;line-height:normal;font-family:Menlo;min-height:13px"><span></span><br></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> Fermi energy: -0.260963</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo;min-height:13px"><span></span><br></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> </span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> HOMO-LUMO gap: 0.226790 = 6.17 eV</span></p><div><span>------------------------------<wbr>---------</span></div><div><span><br></span></div><div><span>However, with CARTESIAN OFF, I get only the orbital #4:</span></div><div><span><br></span></div><div><span>------------------------------<wbr>----------</span></div><div><span><p style="font-size:11px;line-height:normal;font-family:Menlo"><span> MO EIGENVALUES, MO OCCUPATION NUMBERS, AND SPHERICAL MO EIGENVECTORS</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo;min-height:13px"><span></span><br></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 4</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> -0.260963</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo;min-height:13px"><span></span><br></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 2.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo;min-height:13px"><span></span><br></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 1 1 O 2s 0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 2 1 O 3s 0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 3 1 O 4s 0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 4 1 O 3py -0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 5 1 O 3pz 0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 6 1 O 3px -0.917807</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 7 1 O 4py 0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 8 1 O 4pz -0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 9 1 O 4px 0.020934</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 10 1 O 5py 0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 11 1 O 5pz -0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 12 1 O 5px -0.003042</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 13 1 O 3d-2 -0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 14 1 O 3d-1 0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 15 1 O 3d0 0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 16 1 O 3d+1 -0.005256</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 17 1 O 3d+2 -0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 18 1 O 4d-2 0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 19 1 O 4d-1 -0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 20 1 O 4d0 0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 21 1 O 4d+1 -0.022359</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 22 1 O 4d+2 -0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 23 1 O 5s -0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 24 1 O 5py -0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 25 1 O 5pz 0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 26 1 O 5px 0.024750</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo;min-height:13px"><span></span><br></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 27 2 H 1s 0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 28 2 H 2s 0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 29 2 H 3s -0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 30 2 H 2py 0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 31 2 H 2pz -0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 32 2 H 2px 0.011575</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 33 2 H 3py 0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 34 2 H 3pz -0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 35 2 H 3px 0.028901</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 36 2 H 4s -0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 37 2 H 4py 0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 38 2 H 4pz -0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 39 2 H 4px 0.026994</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo;min-height:13px"><span></span><br></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 40 3 H 1s 0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 41 3 H 2s -0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 42 3 H 3s 0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 43 3 H 2py 0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 44 3 H 2pz -0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 45 3 H 2px 0.011575</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 46 3 H 3py -0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 47 3 H 3pz -0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 48 3 H 3px 0.028901</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 49 3 H 4s 0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 50 3 H 4py 0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 51 3 H 4pz 0.000000</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> 52 3 H 4px 0.026994</span></p> <p style="font-size:11px;line-height:normal;font-family:Menlo;min-height:13px"><span></span><br></p> <p style="font-size:11px;line-height:normal;font-family:Menlo;min-height:13px"><span></span><br></p> <p style="font-size:11px;line-height:normal;font-family:Menlo"><span> Fermi energy: -0.260963</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span>-----------------------------</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span><br></span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span>Why the MO #5 is not printed in spherical coordinate?</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span><br></span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span>Best,</span></p><p style="font-size:11px;line-height:normal;font-family:Menlo"><span>Xiaoming</span></p></span></div><br>On Friday, October 12, 2018 at 5:51:08 AM UTC-4, Vladimir Rybkin wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Xiaoming,<div><br></div><div>it is. I mean that by default it is not printed, only the basis is. So, please attach your input/output: this will help to stay on equal footing.</div><div><br></div><div>Yours,</div><div><br></div><div>Vladimir<br><br>четверг, 11 октября 2018 г., 16:26:57 UTC+2 пользователь Xiaoming Wang написал:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Vladimir,<div><br></div><div>Are you saying the option of printing MOs with eigenvector on is not printing the MO coefficients?</div><div><br></div><div>Best,</div><div>Xiaoming<br><br>On Thursday, October 11, 2018 at 9:36:49 AM UTC-4, Vladimir Rybkin wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear Xiaoming,<div><br></div><div>by default CP2K does not print MO coefficients. It prints the basis, which is Cartesian (e.i. 6 d orbitals instead of 5). So, input/output would be helpful.</div><div><br></div><div>Yours,</div><div><br></div><div>Vladimir<br><br>четверг, 11 октября 2018 г., 15:24:16 UTC+2 пользователь Xiaoming Wang написал:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi,<div><br></div><div>I just missed one thing. The MOs could also be printed with spherical coordinates. Then, it should be consistent.</div><div><br></div><div>Best,</div><div>Xiaoming<br><br>On Thursday, October 11, 2018 at 9:21:56 AM UTC-4, Xiaoming Wang wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Vladimir,<div><br></div><div>Thanks for your reply. I'm wondering why the eigenvectors and overlap are not using the same basis (or with same basis but different coordinates).</div><div><br></div><div>Best,</div><div>Xiaoming<br><br>On Thursday, October 11, 2018 at 4:28:32 AM UTC-4, Vladimir Rybkin wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear Xiaoming,<div><br></div><div>the 6 d orbitals correspond to the six components in the Cartesian coordinates, whereas 5 d orbitals are the same expressed is spherical coordinates. The coordinate transformation is involved in between. </div><div><br></div><div>Yours,</div><div><br></div><div>Vladimir<br><br>среда, 10 октября 2018 г., 16:16:56 UTC+2 пользователь Xiaoming Wang написал:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hello,<div><br></div><div>I'd like to check the equation C^T S C = 1, where C is MO coefficient and S the overlap matrix. </div><div>However, for my test calculation, C is 54*1 while S is 52*52. I checked the basis details, for the</div><div>MO coefficient, there are 6 d orbitals (dx2,dxy,dxz,dy2,dyz,dz2), while for the overlap matrix,</div><div>only 5 d orbitals (d-2,d-1,d0,d+1,d+2) are printed. Why is like that? Am I missing something?</div><div><br></div><div>Best,</div><div>Xiaoming </div></div></blockquote></div></div></blockquote></div></div></blockquote></div></div></blockquote></div></div></blockquote></div></div></blockquote></div></div></blockquote></div></div></blockquote></div></div></blockquote></div></div></blockquote></div></div>