<div dir="ltr">Dear Vladimir,<div><br></div><div>Thanks for you reply.</div><div>Actually, I need to calculate the transition dipole moment <Psi_i| r |Psi_j> between two states, which can be evaluated using the MO_CUBE files.</div><div>However, as suggested by Prof. Hutter, using the MO coefficient should be more efficient. Since the MO coefficients are not normalized, the overlap</div><div>matrix is needed, so the dipole is then <C_i|S*r|C_j>. Unfortunately, I cannot use the diagonalization method. I will have to use ROKS with OT in my</div><div>applications. For OT, the ADDED_MOS are not allowed. So I cannot find a way around by now. Maybe I need turn back to MO_CUBE files.</div><div><br></div><div>Best,</div><div>Xiaoming<br><br>On Monday, October 15, 2018 at 4:38: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>I believe that for no application does one need to check CSC=1. If you do it for fun than you might as well try diagonalization once.</div><div><br></div><div>Yours,</div><div><br></div><div>Vladimir<br><br>пятница, 12 октября 2018 г., 17:47: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 suggestions. </div><div><br></div><div>Best,</div><div>Xiaoming</div><div><br></div><div>For my applications, however, I do need the OT method. </div><div>So I cannot use ADDED_MOS. <br><br>On Friday, October 12, 2018 at 10:04:18 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>you need to do the following:</div><div>1) use DIAGONALIZATION ins SCF</div><div>2) in SCF use ADDED_MOS #very large number</div><div><br></div><div>Yours,</div><div><br></div><div>Vladimir<br><br>пятница, 12 октября 2018 г., 14:53:59 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>I'm sorry I made a mistake. I mean could the overlap matrix be transformed to cartesian coordinates? Because I only have the MO of LUMO printed in cartesian coordinates.</div><div><br></div><div>Best,</div><div>Xiaoming <br><br>On Friday, October 12, 2018 at 8:48:45 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>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></blockquote></div></div></blockquote></div></div></blockquote></div></div></blockquote></div></div></blockquote></div></div>