<div dir="ltr">Hi Thomas,<div><br></div><div>Thanks for your comments. I tried different setups for the SCF loop and also tried the </div><div>diagonalization method but with no success. I have a concern about the XC potential</div><div>because with US or SPZ if the XC potentials were modified slightly different from the</div><div>default ones I got the similar convergence issue. However, with the default correct XC</div><div>potentials the US and SPZ scheme could get convergence fairly easily. Below is my</div><div>implementation which is coded in qs_vxc.F.</div><div><br></div><div>---</div><div><div> ! The Sio-Giustino scheme </div><div> ! E(rho_alpha,rho_beta)-b/2*(E(rho_alpha,rho_beta) + E(rho_beta-m,rho_beta) -2*E(rho_beta,rho_beta))</div><div> IF (dft_control%sic_method_id .EQ. sic_mauri_giustino .AND. .NOT. sic_scaling_b_zero) THEN</div><div> </div><div> rho_r(1)%pw => rho_struct_r(2)%pw</div><div> rho_r(2)%pw => rho_struct_r(2)%pw</div><div> </div><div> </div><div> IF (rho_g_valid) THEN</div><div> rho_g(1)%pw => rho_struct_g(2)%pw</div><div> rho_g(2)%pw => rho_struct_g(2)%pw</div><div> ENDIF</div><div><br></div><div> IF (my_just_energy) THEN</div><div> exc_dd = xc_exc_calc(rho_r=rho_r, tau=tau, &</div><div> rho_g=rho_g, xc_section=xc_section, &</div><div> pw_pool=xc_pw_pool)</div><div> ELSE</div><div> ! virial untested</div><div> CPASSERT(.NOT. compute_virial)</div><div> CALL xc_vxc_pw_create1(vxc_rho=my_vxc_rho, vxc_tau=my_vxc_tau, rho_r=rho_r, &</div><div> rho_g=rho_g, tau=tau, exc=exc_dd, &</div><div> xc_section=xc_section, &</div><div> pw_pool=xc_pw_pool, &</div><div> compute_virial=.FALSE., &</div><div> virial_xc=virial_xc_tmp)</div><div> END IF</div><div> </div><div> ! and take care of the potential</div><div> IF (.NOT. my_just_energy) THEN</div><div> ! </div><div> vxc_rho(2)%pw%cr3d = vxc_rho(2)%pw%cr3d/2.0_dp + 2.0_dp*my_vxc_rho(1)%pw%cr3d</div><div> CALL pw_release(my_vxc_rho(1)%pw)</div><div> CALL pw_release(my_vxc_rho(2)%pw)</div><div> DEALLOCATE (my_vxc_rho)</div><div> ENDIF</div><div> </div><div> ALLOCATE (rho_m_rspace(2), rho_m_gspace(2))</div><div> CALL pw_pool_create_pw(xc_pw_pool, rho_m_gspace(1)%pw, &</div><div> use_data=COMPLEXDATA1D, &</div><div> in_space=RECIPROCALSPACE)</div><div> CALL pw_pool_create_pw(xc_pw_pool, rho_m_rspace(1)%pw, &</div><div> use_data=REALDATA3D, &</div><div> in_space=REALSPACE)</div><div> CALL pw_copy(rho_struct_r(1)%pw, rho_m_rspace(1)%pw)</div><div> CALL pw_axpy(rho_struct_r(2)%pw, rho_m_rspace(1)%pw, alpha=-1._dp)</div><div> CALL pw_copy(rho_struct_g(1)%pw, rho_m_gspace(1)%pw)</div><div> CALL pw_axpy(rho_struct_g(2)%pw, rho_m_gspace(1)%pw, alpha=-1._dp)</div><div> </div><div> CALL pw_pool_create_pw(xc_pw_pool, rho_m_gspace(2)%pw, &</div><div> use_data=COMPLEXDATA1D, &</div><div> in_space=RECIPROCALSPACE)</div><div> CALL pw_pool_create_pw(xc_pw_pool, rho_m_rspace(2)%pw, &</div><div> use_data=REALDATA3D, &</div><div> in_space=REALSPACE)</div><div> CALL pw_copy(rho_struct_r(2)%pw, rho_m_rspace(2)%pw)</div><div> CALL pw_copy(rho_struct_g(2)%pw, rho_m_gspace(2)%pw)</div><div> </div><div> CALL pw_axpy(rho_m_rspace(1)%pw, rho_m_rspace(2)%pw, alpha=-1._dp)</div><div> CALL pw_axpy(rho_m_gspace(1)%pw, rho_m_gspace(2)%pw, alpha=-1._dp)</div><div> </div><div> rho_r(1)%pw => rho_m_rspace(2)%pw</div><div> rho_r(2)%pw => rho_struct_r(2)%pw</div><div> </div><div> IF (rho_g_valid) THEN</div><div> rho_g(1)%pw => rho_m_gspace(2)%pw</div><div> rho_g(2)%pw => rho_struct_g(2)%pw</div><div> ENDIF</div><div> </div><div> IF (my_just_energy) THEN</div><div> exc_dm = xc_exc_calc(rho_r=rho_r, tau=tau, &</div><div> rho_g=rho_g, xc_section=xc_section, &</div><div> pw_pool=xc_pw_pool)</div><div> ELSE</div><div> ! virial untested</div><div> CPASSERT(.NOT. compute_virial)</div><div> CALL xc_vxc_pw_create1(vxc_rho=my_vxc_rho, vxc_tau=my_vxc_tau, rho_r=rho_r, &</div><div> rho_g=rho_g, tau=tau, exc=exc_dm, &</div><div> xc_section=xc_section, &</div><div> pw_pool=xc_pw_pool, &</div><div> compute_virial=.FALSE., &</div><div> virial_xc=virial_xc_tmp)</div><div> END IF</div><div> </div><div> IF (.NOT. my_just_energy) THEN</div><div> vxc_rho(1)%pw%cr3d = vxc_rho(1)%pw%cr3d/2.0_dp + my_vxc_rho(1)%pw%cr3d/2.0_dp</div><div> vxc_rho(2)%pw%cr3d = vxc_rho(2)%pw%cr3d - my_vxc_rho(2)%pw%cr3d/2.0_dp - my_vxc_rho(1)%pw%cr3d</div><div> CALL pw_release(my_vxc_rho(1)%pw)</div><div> CALL pw_release(my_vxc_rho(2)%pw)</div><div> DEALLOCATE (my_vxc_rho)</div><div> ENDIF</div><div> </div><div> exc = exc/2.0_dp + exc_dd - exc_dm/2.0_dp</div><div> </div><div> </div><div> DO ispin = 1, 2</div><div> CALL pw_pool_give_back_pw(xc_pw_pool, rho_m_rspace(ispin)%pw)</div><div> CALL pw_pool_give_back_pw(xc_pw_pool, rho_m_gspace(ispin)%pw)</div><div> ENDDO</div><div> DEALLOCATE (rho_m_rspace)</div><div> DEALLOCATE (rho_m_gspace)</div><div> ENDIF</div><div>--- </div></div><div><br></div><div>Best,</div><div>Xiaoming<br><br>On Wednesday, January 1, 2020 at 5:02:49 AM UTC-5, tkuehne wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div style="word-wrap:break-word">Dear Xiaoming, <div><br></div><div>happy new year as well! Though I can’t say anything definite, but having a quick glance </div><div>into Feliciano’s PRB paper I don’t see any apparent mistake. Regarding your convergence </div><div>issue (which shouldn’t be taken as a test to prove correctness or the opposite!) I want to </div><div>comment that the excessive usage of the OUTER_LOOP and hence too short inner loop </div><div>is not the smartest thing to do and may also lead to the behavior you observe. This is </div><div>particularly true when using the CG minimizer in which case you are only computing 10 </div><div>gradients before restarting the inner loop with the default STEPSIZE of 0.15, which in my </div><div>experience is basically always too large except for trivially cases such as water … </div><div><br></div><div>Cheers, </div><div>Thomas</div><div><br><div><blockquote type="cite"><div>Am 01.01.2020 um 06:02 schrieb Xiaoming Wang <<a href="javascript:" target="_blank" gdf-obfuscated-mailto="_vebk44eBgAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">wx...@gmail.com</a>>:</div><br><div><div dir="ltr" style="font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">Happy new year!<div><br></div><div>Really appreciate if anyone could confirm my derivation of the XC potential.</div><div><br></div><div>Best,</div><div>Xiaoming<br><br>On Friday, December 27, 2019 at 11:50:58 AM UTC-5, Xiaoming Wang wrote:<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hello all,<div><br></div><div>I'm trying to implement the SIC scheme proposed in PRB 99, 235139 (2019) by Giustino et al. </div><div><br></div><div>The new SIC scheme is a minor modification on the one proposed by Mauri [PRB 71, 205210 (2005)]. </div><div><br></div><div>Since the Mauri SIC is implemented in CP2K, incorporation of the Giustino scheme should not need </div><div><br></div><div>much coding work. The Hartree self-interaction part is similar for all the schemes, so I'm only interested</div><div><br></div><div>in the XC part. The implementation of the XC part of the Mauri scheme is coded in qs_vxc.F. </div><div><br></div><div>There are two types of SIC functionals for the Mauri scheme: Mauri_SPZ and Mauri_US. </div><div><br></div><div>The corresponding SIC corrrected XC functionals are:</div><div><br></div><div> Mauri_SPZ: Exc = Exc [ alpha, beta ] - Exc [ alpha - beta, 0 ]</div><div><br></div><div> Mauri_US: Exc = Exc [ alpha, beta ] - Exc [ alpha, beta ] + Exc [ beta, beta ]</div><div><br></div><div>And the XC potentials are:</div><div><br></div><div> Mauri_SPZ: Vxc_up = Vxc_up [ alpha, beta ] - Vxc_up [ alpha - beta, 0 ]</div><div> Vxc_dn = Vxc_dn [ alpha, beta ] + Vxc_up [ alpha - beta, 0 ]</div><div><br></div><div> Mauri_US: Vxc_up = 0</div><div> <span> </span>Vxc_dn = Vxc_up [ beta, beta ] + Vxc_dn [ beta, beta ]</div><div><br></div><div>For the Giustino scheme, the XC functional is:</div><div><br></div><div> <span> </span>Exc = 0.5*Exc [ alpha, beta ] + Exc [ beta, beta ] - 0.5* Exc [ beta - m, beta ]</div><div><br></div><div>where m = alpha - beta. Based on my understanding, the XC potentials are:</div><div><br></div><div> Vxc_up = 0.5*Vxc_up [ alpha, beta ] + 0.5*Vxc_up [ beta - m, beta ]</div><div><br></div><div> Vxc_dn = 0.5*Vxc_dn [ alpha, beta ] + Vxc_up [ beta, beta ] + Vxc_dn [ beta, beta ]</div><div> </div><div> <span> </span>-0.5*Vxc_dn [ beta - m, beta ] - Vxc_up [ beta - m, beta ]</div><div><br></div><div>Please can anyone correct me if my understanding above is wrong. </div><div><br></div><div>I implemented the Giustino SIC based on the above equations by slightly modifying the US scheme </div><div><br></div><div>as already in CP2K. However, I never get convergence for the SCF calculations (see part of the log</div><div><br></div><div>file below). For the same system, the SPZ and US schemes both can easily get convergence.</div><div><br></div><div>So I'm wondering if anyone could give me any comments on this?</div><div><br></div><div><br></div><div><br></div><div>Best,</div><div><br></div><div>Xiaoming</div><div><br></div><div><br></div><div><br></div><div>--</div><div><div> 1 OT CG 0.15E+00 14.2 0.00014039 -8153.9988425820 -8.15E+03</div><div> 2 OT LS 0.10E+00 5.2 -8154.0035499679</div><div> 3 OT CG 0.10E+00 10.6 0.00020232 -8154.0032891957 -4.45E-03</div><div> 4 OT LS 0.46E-01 6.0 -8154.0015689085</div><div> 5 OT CG 0.46E-01 10.8 0.00009816 -8154.0041296811 -8.40E-04</div><div> 6 OT LS 0.31E-01 5.1 -8154.0048751897</div><div> 7 OT CG 0.31E-01 13.1 0.00007670 -8154.0046628243 -5.33E-04</div><div> 8 OT LS 0.26E-01 5.2 -8154.0051281540</div><div> 9 OT CG 0.26E-01 12.0 0.00007395 -8154.0050568034 -3.94E-04</div><div> <span> </span>10 OT LS 0.21E-01 6.2 -8154.0053917653</div><div> <span> </span>11 OT CG 0.21E-01 11.2 0.00007210 -8154.0053264707 -2.70E-04</div><div> <span> </span>12 OT LS 0.16E-01 5.3 -8154.0055646935</div><div> <span> </span>13 OT CG 0.16E-01 10.4 0.00005042 -8154.0055112362 -1.85E-04</div><div> <span> </span>14 OT LS 0.15E-01 5.3 -8154.0056327656</div><div> <span> </span>15 OT CG 0.15E-01 11.4 0.00004942 -8154.0056270301 -1.16E-04</div><div> <span> </span>16 OT LS 0.14E-01 5.3 -8154.0057343493</div><div> <span> </span>17 OT CG 0.14E-01 11.5 0.00004857 -8154.0057262676 -9.92E-05</div><div> <span> </span>18 OT LS 0.13E-01 5.4 -8154.0058189369</div><div> <span> </span>19 OT CG 0.13E-01 10.7 0.00004785 -8154.0058097707 -8.35E-05</div><div> <span> </span>20 OT LS 0.11E-01 5.0 -8154.0058883901</div><div> <span> </span>----------------------------<wbr>------- OT ------------------------------<wbr>---------</div><div> <span> </span>----------------------------<wbr>------- OT ------------------------------<wbr>---------</div><div> 1 OT CG 0.15E+00 21.8 0.00004727 -8154.0058790748 -6.93E-05</div><div> 2 OT LS 0.11E+00 4.7 -8154.0065750586</div><div> 3 OT CG 0.11E+00 10.9 0.00004074 -8154.0064374937 -5.58E-04</div><div> 4 OT LS 0.13E+00 5.4 -8154.0072008089</div><div> 5 OT CG 0.13E+00 11.3 0.00002783 -8154.0072885701 -8.51E-04</div><div> 6 OT LS 0.51E+00 5.5 -8154.0078327851</div><div> 7 OT CG 0.51E+00 11.1 0.00007874 -8154.0089341020 -1.65E-03</div><div> 8 OT LS 0.21E+00 5.2 -8154.0051598661</div><div> 9 OT CG 0.21E+00 11.3 0.00017344 -8154.0086760074 2.58E-04</div><div> <span> </span>10 OT LS 0.98E-01 5.3 -8154.0049564999</div><div> <span> </span>11 OT CG 0.98E-01 11.2 0.00006510 -8154.0082092588 4.67E-04</div><div> <span> </span>12 OT LS 0.58E-01 5.7 -8154.0086320114</div><div> <span> </span>13 OT CG 0.58E-01 10.8 0.00003585 -8154.0085028738 -2.94E-04</div><div> <span> </span>14 OT LS 0.10E+00 5.3 -8154.0088394530</div><div> <span> </span>15 OT CG 0.10E+00 12.0 0.00007085 -8154.0090772735 -5.74E-04</div><div> <span> </span>16 OT LS 0.44E-01 5.4 -8154.0084763907</div><div> <span> </span>17 OT CG 0.44E-01 11.1 0.00003408 -8154.0088708613 2.06E-04</div><div> <span> </span>18 OT LS 0.72E-01 5.3 -8154.0090917387</div><div> <span> </span>19 OT CG 0.72E-01 11.4 0.00005874 -8154.0092230991 -3.52E-04</div><div> <span> </span>20 OT LS 0.30E-01 5.5 -8154.0089057509</div><div> <span> </span>----------------------------<wbr>------- OT ------------------------------<wbr>---------</div><div> <span> </span>----------------------------<wbr>------- OT ------------------------------<wbr>---------</div><div> 1 OT CG 0.15E+00 22.0 0.00001421 -8154.0091066057 1.16E-04</div><div> 2 OT LS 0.18E+00 5.9 -8154.0092180076</div><div> 3 OT CG 0.18E+00 10.7 0.00005904 -8154.0092398262 -1.33E-04</div><div> 4 OT LS 0.96E-01 5.5 -8154.0094065127</div><div> 5 OT CG 0.96E-01 11.2 0.00007661 -8154.0094075743 -1.68E-04</div><div> 6 OT LS 0.42E-01 5.6 -8154.0088928557</div><div> 7 OT CG 0.42E-01 10.9 0.00001444 -8154.0092323658 1.75E-04</div><div> 8 OT LS 0.57E-01 5.1 -8154.0092666992</div><div> 9 OT CG 0.57E-01 11.2 0.00001390 -8154.0092787816 -4.64E-05</div><div> <span> </span>10 OT LS 0.68E-01 5.2 -8154.0093186446</div><div> <span> </span>11 OT CG 0.68E-01 10.5 0.00002940 -8154.0093259578 -4.72E-05</div><div> <span> </span>12 OT LS 0.26E-01 5.0 -8154.0092084126</div><div> <span> </span>13 OT CG 0.26E-01 11.0 0.00001386 -8154.0092843540 4.16E-05</div><div> <span> </span>14 OT LS 0.13E-01 5.4 -8154.0092839420</div><div> <span> </span>15 OT CG 0.13E-01 11.0 0.00001387 -8154.0092841905 1.64E-07</div><div> <span> </span>16 OT LS 0.16E-01 5.4 -8154.0092934733</div><div> <span> </span>17 OT CG 0.16E-01 11.0 0.00001375 -8154.0092962332 -1.20E-05</div><div> <span> </span>18 OT LS 0.20E-01 5.4 -8154.0093077912</div><div> <span> </span>19 OT CG 0.20E-01 10.4 0.00001362 -8154.0093106170 -1.44E-05</div><div> <span> </span>20 OT LS 0.24E-01 5.3 -8154.0093242773</div><div> <span> </span>----------------------------<wbr>------- OT ------------------------------<wbr>---------</div><div> <span> </span>----------------------------<wbr>------- OT ------------------------------<wbr>---------</div><div> 1 OT CG 0.15E+00 22.4 0.00001343 -8154.0093267803 -1.62E-05</div><div> 2 OT LS 0.15E+00 5.4 -8154.0094121812</div><div> 3 OT CG 0.15E+00 11.6 0.00007460 -8154.0094132278 -8.64E-05</div><div> 4 OT LS 0.74E-01 5.3 -8154.0093027278</div><div> 5 OT CG 0.74E-01 10.9 0.00007707 -8154.0094484130 -3.52E-05</div><div> 6 OT LS 0.33E-01 5.6 -8154.0090697156</div><div> 7 OT CG 0.33E-01 14.3 0.00001494 -8154.0093119312 1.36E-04</div><div> 8 OT LS 0.46E-01 6.5 -8154.0093411677</div></div><div><br></div><div>--</div><div><br></div><div><br></div><div><br></div></div></blockquote></div></div><div style="font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><br></div><span style="font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;float:none;display:inline!important">--<span> </span></span><br style="font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><span style="font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;float:none;display:inline!important">You received this message because you are subscribed to the Google Groups "cp2k" group.</span><br style="font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><span style="font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;float:none;display:inline!important">To unsubscribe from this group and stop receiving emails from it, send an email to<span> </span></span><a href="javascript:" style="font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px" target="_blank" gdf-obfuscated-mailto="_vebk44eBgAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@<wbr>googlegroups.com</a><span style="font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;float:none;display:inline!important">.</span><br style="font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><span style="font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;float:none;display:inline!important">To view this discussion on the web visit<span> </span></span><a href="https://groups.google.com/d/msgid/cp2k/a2bed8c8-339a-493f-8c48-27809a74d1a6%40googlegroups.com?utm_medium=email&utm_source=footer" style="font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/msgid/cp2k/a2bed8c8-339a-493f-8c48-27809a74d1a6%40googlegroups.com?utm_medium\x3demail\x26utm_source\x3dfooter';return true;" onclick="this.href='https://groups.google.com/d/msgid/cp2k/a2bed8c8-339a-493f-8c48-27809a74d1a6%40googlegroups.com?utm_medium\x3demail\x26utm_source\x3dfooter';return true;">https://groups.google.<wbr>com/d/msgid/cp2k/a2bed8c8-<wbr>339a-493f-8c48-27809a74d1a6%<wbr>40googlegroups.com</a><span style="font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;float:none;display:inline!important">.</span></div></blockquote></div><br></div><br><br><div>
<div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">==============================</div><div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">Thomas D. Kühne</div><div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">Dynamics of Condensed Matter</div><div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">Chair of Theoretical Chemistry</div><div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">University of Paderborn</div><div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">Warburger Str. 100</div><div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">D-33098 Paderborn</div><div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">Germany</div><div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><a href="javascript:" target="_blank" gdf-obfuscated-mailto="_vebk44eBgAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">td...@mail.upb.de</a></div><div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">+49/(0)5251/60-5726</div></div></div></div>
</div>
<br></div></blockquote></div></div>