<div dir="ltr">Hi Nuri,<div>Thanks for your clarification.</div><div><br></div><div>Best,</div><div>Xiaoming<br><br>On Friday, October 12, 2018 at 7:47:40 AM UTC-4, Nuri Yazdani 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><div>You are right, the MOs are not normalised. However, and I am not 100% sure about this, I think that they shouldnt necessarily be normalized, since the basis functions (i.e. the gaussian orbital on each atom) are not orthogonal to one another. This is where the overlap matrix comes in.. </div><div>So, I think, if one wants to construct for example the electron density from the MOs and the basis functions, one should use the MO coefficients as is (i.e. do not renormalise the MOs yourself). </div><div>I think, and again I am not totally confident about this, that if you print out the overlap matrix S, you can get the dipole moment (missing some prefactors, and for example the x component) by</div><div><br></div><div>A_ij,x = <mo_i| S x |mo_j></div><div><br></div><div>The problem that I have is that I have already run really large MD calculations, and I only was printing out the MOs. I think the best I can do then is to use</div><div><br></div><div>A_ij,x ~ <mo_i| x |mo_j></div><div><br></div><div>which from my understanding, ignores any contribution to A_ij,x from orbitals on two different atoms. For my system, A_ij,x will be dominated by s -> p transitions on one particular atomic species, so this approximation should be OK, but it may fail terribly for systems where the transition goes from an orbital on one atom type to an orbital on a different atom type..</div><div><br>On Wednesday, October 10, 2018 at 2:49:34 PM UTC+2, 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 Nuri,<div><br></div><div>Are you using the MOS? The mo coefficients are not normalized, how do you solve this problem?</div><div><br></div><div>Best,</div><div>Xiaoming<br><br>On Tuesday, October 9, 2018 at 7:31:19 AM UTC-4, Nuri Yazdani wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Prof. Hutter,<div><br></div><div>I am trying to do such calculations. I have printed out the MOS from my calculations to do this, however, I am quite confused by the spatial representation and ordering of the d-orbitals... </div><div><br></div><div>I am using the MOLOPT-DZVP basis set. For example, for Cs DZVP-MOLOPT-SR-GTH 1 2 0 2 6 3 2 1, I have 3x S, 6x P, and 5x D orbitals (this counting is also consistent with the length of my MOSs).</div><div><br></div><div>The S and P orbitals I can construct from the basis set, and (if CGS are a linear combination of gaussians) the p orbitals are px = x*CGS, py = y*CGS, pz =z*CGS, however, what is the ordering and form of the d orbitals?? i.e. what is r^2 equal to in each of the 5 last lines of this page: <a href="https://www.cp2k.org/basis_sets" rel="nofollow" target="_blank" onmousedown="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fwww.cp2k.org%2Fbasis_sets\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNFpZrgEKEYE6kdbDW5k7sfvb8siug';return true;" onclick="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fwww.cp2k.org%2Fbasis_sets\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNFpZrgEKEYE6kdbDW5k7sfvb8siug';return true;">https://www.cp2k.org/basis_<wbr>sets</a>? </div><div><br></div><div>Cheers,</div><div>Nuri</div><div><br></div><div><br></div><div><br><br>On Thursday, May 31, 2018 at 11:28:04 AM UTC+2, jgh wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex">Hi <br> <br>this property is on our TO DO list. However, I cannot say when <br>it will become available. <br> <br>If you want to calculate it yourself from the MO cube files <br>you need to use the Berry phase algorithm, meaning you need <br> <br>mu = IMAG LOG <phi_a | exp(i*k*r) |phi_b> <br> <br>However, this can be calculated much more efficiently from <br>the atomic basis functions. <br> <br>regards <br> <br>Juerg Hutter <br>------------------------------<wbr>------------------------------<wbr>-- <br>Juerg Hutter Phone : ++41 44 635 4491 <br>Institut für Chemie C FAX : ++41 44 635 6838 <br>Universität Zürich E-mail: <a rel="nofollow">hut...@chem.uzh.ch</a> <br>Winterthurerstrasse 190 <br>CH-8057 Zürich, Switzerland <br>------------------------------<wbr>------------------------------<wbr>--- <br> <br>-----<a rel="nofollow">cp...@googlegroups.com</a> wrote: ----- <br>To: cp2k <<a rel="nofollow">cp...@googlegroups.com</a>> <br>From: Xiaoming Wang <br>Sent by: <a rel="nofollow">cp...@googlegroups.com</a> <br>Date: 05/24/2018 05:07PM <br>Subject: [CP2K:10343] transition dipole moment <br> <br>Hi, <br> <br>Is it possible for CP2K to output the transition dipole moment <phi_A | r | phi_B> between two KS states, say HOMO to LUMO transition? <br>I tried to calculate the integral by printing the MO_CUBE files which are the phi_A and phi_B and then do the integration. This is fine for <br>molecules or in other word non-periodic systems. But for periodic systems, the result was not as expected, maybe due to the position <br>operator is ill-defined in that case. So how to evaluate the transition dipole integral for periodic systems with the MO_CUBE files known? <br> <br>Best, <br>Xiaoming <br> -- <br> You received this message because you are subscribed to the Google Groups "cp2k" group. <br> To unsubscribe from this group and stop receiving emails from it, send an email to <a rel="nofollow">cp2k+...@googlegroups.com</a>. <br> To post to this group, send email to <a rel="nofollow">cp...@googlegroups.com</a>. <br> Visit this group at <a href="https://groups.google.com/group/cp2k" rel="nofollow" target="_blank" onmousedown="this.href='https://groups.google.com/group/cp2k';return true;" onclick="this.href='https://groups.google.com/group/cp2k';return true;">https://groups.google.com/<wbr>group/cp2k</a>. <br> For more options, visit <a href="https://groups.google.com/d/optout" rel="nofollow" target="_blank" onmousedown="this.href='https://groups.google.com/d/optout';return true;" onclick="this.href='https://groups.google.com/d/optout';return true;">https://groups.google.com/d/<wbr>optout</a>. <br> <br></blockquote></div></div></blockquote></div></div></blockquote></div></div></div></blockquote></div></div>