<div dir="ltr"><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000"><div class="gmail_default" style="font-family:verdana,sans-serif">Dear Viktor,</div><div class="gmail_default" style="font-family:verdana,sans-serif"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif">As Fabien says, two different :MMT indicates that you have two minima. The mixer only finds local minima, in general not the global although you can be lucky if they are close. You probably have to look at the f density matrix which I suspect is slightly different for the two cases. Take the density matrix into something else and look at the eigenvectors [1] -- they are probably different (different parity or something).</div><div class="gmail_default" style="font-family:verdana,sans-serif"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif">[1] I have some code "somewhere" where I added this to lapw2, and I think it is a good general diagnostic particularly for f electrons.</div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Jan 22, 2019 at 1:47 PM <a href="mailto:tran@theochem.tuwien.ac.at">tran@theochem.tuwien.ac.at</a> <<a href="mailto:tran@theochem.tuwien.ac.at">tran@theochem.tuwien.ac.at</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi,<br>
<br>
Maybe the two different energies with 30000 k-points correspond<br>
to different local minima.<br>
<br>
If not already done, then do the following:<br>
1) restore the calculation with 20000 k-points (restore -f ...)<br>
2) execute "x kgen" to generate the mesh with 30000 k-points<br>
3) start the calculation (runsp_lapw -ec 0.0001 -cc 0.0001)<br>
<br>
FT<br>
<br>
<br>
On Tuesday 2019-01-22 20:34, Victor Zenou wrote:<br>
<br>
>Date: Tue, 22 Jan 2019 20:34:09<br>
>From: Victor Zenou <<a href="mailto:zanov@post.bgu.ac.il" target="_blank">zanov@post.bgu.ac.il</a>><br>
>Reply-To: A Mailing list for WIEN2k users <<a href="mailto:wien@zeus.theochem.tuwien.ac.at" target="_blank">wien@zeus.theochem.tuwien.ac.at</a>><br>
>To: A Mailing list for WIEN2k users <<a href="mailto:wien@zeus.theochem.tuwien.ac.at" target="_blank">wien@zeus.theochem.tuwien.ac.at</a>><br>
>Subject: Re: [Wien] proper k-points for Nd<br>
><br>
>Thanks Peter and Laurie for your advises. I hope we all stay friends after my finding :)<br>
>The reason I didn't progress with higher k-points (I guess it doesn't work backwards, higher k-points to lower) is that I actually tried that few times and it<br>
>worked (at other phase),but sometimes I got error message after the first cycle. In order to be systematic I tried to avoid that.<br>
>According to Peter's advice I tried again and got an error message (I can't recall if it was the same as previous) at lapw2 -up of the first cycle.<br>
>uplapw2.error file showed: " 'LAPW2" - semicore band-ranges too large, ghostbands ? ". I 've deleted that file and re- run SCF while unmarking the "iterative<br>
>diag" box and for some reason it worked. It worked another time, as I run 2 parallel sessions/cases with different k-points (Nd-SP and Nd-SP2).<br>
>I have good news and bad ones. The good news <br>
>20,000 k-points gave the same energy as 15,000 (ec 0.0001, cc 0.0001): 77040.469768(3) and 77040.469768(2), with same MMTOT (16.645, 16.647). :CHARGE<br>
>convergence -.000216 and -.000422, for 15000 and 20000 k-points!<br>
><br>
>Now for the troubling news:  <br>
>30,000 k-points gave something different -77040.506092(3), MMTOT=16.091; CHARGE convergence -.000070.<br>
>Running the same at different case, as I said I run 2 parallel cases using the same case.struc file) I got for the same 30,000 k-points -77040.4874669(3),<br>
>MMTOT=16.682;  CHARGE convergence -.000070.<br>
>It seems that using same k-points gave me two different energies. As I understand repeating calculation should give the same results. And if these two<br>
>calculations are the same, unless I did something stuoid,  should give the same result. Also, the different MMTOT may suggest something wrong with the spin<br>
>polarized calculations. Maybe I should do the k-test without spin polarization, as I usually do.<br>
>I used both to try 60,000 points, and later I will follow Laurie's and Peter HDLOs.<br>
>Victor<br>
><br>
><br>
><br>
>‫בתאריך יום א׳, 20 בינו׳ 2019 ב-20:01 מאת ‪Peter Blaha‬ <‪<a href="mailto:pblaha@theochem.tuwien.ac.at" target="_blank">pblaha@theochem.tuwien.ac.at</a>‬>:‬<br>
>      I disagree with Lauries point: Eventually I want to do even an<br>
>      "unphysical" calculation, of course it might be you have to pay a price<br>
>      for it, because this unphysical model as eg. a spike in the DOS at EF.<br>
>      And mixer should lead to convergence, but sometimes has problems<br>
>      (Lauries mixer is still one of the big advantages of WIEN2k. I know<br>
>      other codes, which would not be able to converge a 5f system at all).<br>
><br>
>      In any case, put a HDLO (5f) and increase LVNS=6 first.<br>
><br>
>      Most importantly: why are you starting from scratch when doing the<br>
>      k-mesh tests. This is very error prone and in particular burns cpu-hours<br>
>      completely unnecessarily.<br>
><br>
>      One of my main "advises" is always: do an init_lapw only ONCE ! (except<br>
>      you have to change spheres).<br>
><br>
>      You will get much better internal convergence when you continue using<br>
>      the previous calculation:<br>
><br>
>      loop:<br>
>      runsp ....<br>
>      save XXX-k-points<br>
>      x kgen (increase the mesh always at least by a factor of 2 (maybe even 4).<br>
>      goto loop<br>
><br>
>      I'm not so much surprised about the k-mesh for a magnetic 5d-element<br>
>      calculation. I start with 10000 k-points for every metallic spd-element<br>
>      and test it with 50000 k. You may need 100 000 or even more. (optics for<br>
>      fcc Al is not fully converged with 1 000 000 k-points !). It does not<br>
>      cost much cpu, since you should start with a nearly converged calculation.<br>
><br>
>      Your :DIS is still quite large, try to bring it below 0.0001<br>
><br>
>      You can try TEMPS 0.002 as Laurence advised. This will most likely<br>
>      improve convergence, but most likely also reduce MMT a little bit, which<br>
>      is ok when you compare with room temp experiments, but not for 4k<br>
>      measurements ...<br>
><br>
><br>
><br>
>      Am 20.01.2019 um 17:50 schrieb Laurence Marks:<br>
>      > Sorry Fabien, I am dubious that k-points is the issue, I suspect a bad<br>
>      > model.<br>
>      ><br>
>      > Let me explain, repeating some things I have said before or are in the<br>
>      > mixer notes albeit perhaps buried.<br>
>      ><br>
>      > An "easy" assumption that many people make is that dft problems are well<br>
>      > posed, so have "simple" fixed-point solutions. In fact this is not true.<br>
>      > In reality dft problems have stacks of numerical issues which leads to<br>
>      > "noise". You can think about this in terms of finding the minimum of<br>
>      > "F(x) = x*x + Constant*Noise(x)" where Noise(x) is some randomly<br>
>      > generated noise. The larger Constant is, the larger the deviation is<br>
>      > from x=0 of the apparent minimum. Averaging over more k-points may<br>
>      > reduce Constant, but often in my experience does relatively little.<br>
>      ><br>
>      > A symptom of noise is poor convergence -- which is clearly what is<br>
>      > occurring here.<br>
>      ><br>
>      > If you improve the convergence rate you will reduce the noise. Some<br>
>      > suggestions:<br>
>      ><br>
>      > 1) Use HDLO and reduce RMT<br>
>      > 2) Reduce RMT and use HDLO<br>
>      > 3) Larger LM (e.g. 8)<br>
>      > 4) Larger gmax in case.in2 (e.g. 16)<br>
>      > 5) Use TEMPS 0.0018 (room temperature)<br>
>      > 6) Peter's suggestions<br>
>      ><br>
>      > On Sun, Jan 20, 2019 at 10:27 AM <a href="mailto:tran@theochem.tuwien.ac.at" target="_blank">tran@theochem.tuwien.ac.at</a><br>
>      > <mailto:<a href="mailto:tran@theochem.tuwien.ac.at" target="_blank">tran@theochem.tuwien.ac.at</a>> <<a href="mailto:tran@theochem.tuwien.ac.at" target="_blank">tran@theochem.tuwien.ac.at</a><br>
>      > <mailto:<a href="mailto:tran@theochem.tuwien.ac.at" target="_blank">tran@theochem.tuwien.ac.at</a>>> wrote:<br>
>      ><br>
>      >     Hi,<br>
>      ><br>
>      >     If not too expensive, you could try more k-points (20000 and 25000)<br>
>      >     to see if the fluctuations persist.<br>
>      ><br>
>      >     Questions: how many atoms in the unit cell? Are you optimizing the<br>
>      >     position of atoms?<br>
>      ><br>
>      >     F. Tran<br>
>      ><br>
>      >     On Sunday 2019-01-20 17:17, Victor Zenou wrote:<br>
>      ><br>
>      >      >Date: Sun, 20 Jan 2019 17:17:55<br>
>      >      >From: Victor Zenou <<a href="mailto:zanov@post.bgu.ac.il" target="_blank">zanov@post.bgu.ac.il</a><br>
>      >     <mailto:<a href="mailto:zanov@post.bgu.ac.il" target="_blank">zanov@post.bgu.ac.il</a>>><br>
>      >      >Reply-To: A Mailing list for WIEN2k users<br>
>      >     <<a href="mailto:wien@zeus.theochem.tuwien.ac.at" target="_blank">wien@zeus.theochem.tuwien.ac.at</a><br>
>      >     <mailto:<a href="mailto:wien@zeus.theochem.tuwien.ac.at" target="_blank">wien@zeus.theochem.tuwien.ac.at</a>>><br>
>      >      >To: A Mailing list for WIEN2k users<br>
>      >     <<a href="mailto:wien@zeus.theochem.tuwien.ac.at" target="_blank">wien@zeus.theochem.tuwien.ac.at</a><br>
>      >     <mailto:<a href="mailto:wien@zeus.theochem.tuwien.ac.at" target="_blank">wien@zeus.theochem.tuwien.ac.at</a>>><br>
>      >      >Subject: Re: [Wien] proper k-points for Nd<br>
>      >      ><br>
>      >      >Dear Peter, LDM and wien2k users<br>
>      >      >I did some more calculation with energy and CHARGE convergence<br>
>      >     criteria<br>
>      >      >equal to 10^-4 Ryd and 10^-3, respectively. The energy usually<br>
>      >     converged to<br>
>      >      >1*10-5 or even 1*10^-6. Still I get energy fluctuation up to<br>
>      >     4*10-2 Ryd<br>
>      >      >(please see below)<br>
>      >      >I didn't check AFM yet, but this is not the point. The point is to get<br>
>      >      >"stable" calculation for optimal k-points.<br>
>      >      >I have to say that for each calculation I totally aresed all file<br>
>      >     except<br>
>      >      >case.struc.<br>
>      >      >Also I used all defaults as well as gmax=12, also<br>
>      >     rmt(min)*kmax=8.5, for all<br>
>      >      >calculations.<br>
>      >      >I can repeat all calculation with smaller RMT, such as 2.5, or to<br>
>      >     check<br>
>      >      >HDLO.<br>
>      >      >I think I will start with smaller RMT. Any comments / suggestions?<br>
>      >      >Victor<br>
>      >      ><br>
>      >      >K-pointsIBZ E [Ryd], 10-4 CHARGE convergence CHARGE DISTANCE  MMTOT<br>
>      >      >3000 208 -77040.50516761 -.000504 0.000825   16.00318<br>
>      >      >5000 280 -77040.48748499 -.000544 0.000758   16.68523<br>
>      >      >7500 455 -77040.51668342 -.000862 0.000183   16.62547<br>
>      >      >10000 540 -77040.48747210 -.000817 0.000499   16.67544<br>
>      >      >15000 840 -77040.46976815 -.000216 0.001362   16.64536<br>
>      >      ><br>
>      >      ><br>
>      >      ><br>
>      ><br>
>      ><br>
>      ><br>
>      > --<br>
>      > Professor Laurence Marks<br>
>      > "Research is to see what everybody else has seen, and to think what<br>
>      > nobody else has thought", Albert Szent-Gyorgi<br>
>      > <a href="http://www.numis.northwestern.edu" rel="noreferrer" target="_blank">http://www.numis.northwestern.edu</a> <<a href="http://www.numis.northwestern.edu" rel="noreferrer" target="_blank">http://www.numis.northwestern.edu</a>> ;<br>
>      > Corrosion in 4D: <a href="http://MURI4D.numis.northwestern.edu" rel="noreferrer" target="_blank">MURI4D.numis.northwestern.edu</a><br>
>      > <<a href="http://MURI4D.numis.northwestern.edu" rel="noreferrer" target="_blank">http://MURI4D.numis.northwestern.edu</a>><br>
>      > Partner of the CFW 100% program for gender equity,<br>
>      > <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__www.cfw.org_100-2Dpercent&d=DwIDaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=l1d2pIyr6QG4zudFD9mIvqSRlEOa2QsiTeaYL3aqbw4&s=1e2MOt8EgguN2ubmBL7trr8gJAXljkp2s01vt13BpLc&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__www.cfw.org_100-2Dpercent&d=DwIDaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=l1d2pIyr6QG4zudFD9mIvqSRlEOa2QsiTeaYL3aqbw4&s=1e2MOt8EgguN2ubmBL7trr8gJAXljkp2s01vt13BpLc&e=</a> <<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__www.cfw.org_100-2Dpercent&d=DwIDaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=l1d2pIyr6QG4zudFD9mIvqSRlEOa2QsiTeaYL3aqbw4&s=1e2MOt8EgguN2ubmBL7trr8gJAXljkp2s01vt13BpLc&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__www.cfw.org_100-2Dpercent&d=DwIDaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=l1d2pIyr6QG4zudFD9mIvqSRlEOa2QsiTeaYL3aqbw4&s=1e2MOt8EgguN2ubmBL7trr8gJAXljkp2s01vt13BpLc&e=</a>><br>
>      > Co-Editor, Acta Cryst A<br>
>      ><br>
>      > _______________________________________________<br>
>      > Wien mailing list<br>
>      > <a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
>      > <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__zeus.theochem.tuwien.ac.at_mailman_listinfo_wien&d=DwIDaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=l1d2pIyr6QG4zudFD9mIvqSRlEOa2QsiTeaYL3aqbw4&s=HqjVLD2WpfrSwRKFbqf7V88ADhDECm1QctI1U57Y7xA&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__zeus.theochem.tuwien.ac.at_mailman_listinfo_wien&d=DwIDaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=l1d2pIyr6QG4zudFD9mIvqSRlEOa2QsiTeaYL3aqbw4&s=HqjVLD2WpfrSwRKFbqf7V88ADhDECm1QctI1U57Y7xA&e=</a><br>
>      > SEARCH the MAILING-LIST at:  <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__www.mail-2Darchive.com_wien-40zeus.theochem.tuwien.ac.at_index.html&d=DwIDaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=l1d2pIyr6QG4zudFD9mIvqSRlEOa2QsiTeaYL3aqbw4&s=_FHptNqbf_hSe5ztGYQwgom849ky-4GtwrFrc8MOZhQ&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__www.mail-2Darchive.com_wien-40zeus.theochem.tuwien.ac.at_index.html&d=DwIDaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=l1d2pIyr6QG4zudFD9mIvqSRlEOa2QsiTeaYL3aqbw4&s=_FHptNqbf_hSe5ztGYQwgom849ky-4GtwrFrc8MOZhQ&e=</a><br>
>      ><br>
><br>
>      --<br>
>      --------------------------------------------------------------------------<br>
>      Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna<br>
>      Phone: +43-1-58801-165300             FAX: +43-1-58801-165982<br>
>      Email: <a href="mailto:blaha@theochem.tuwien.ac.at" target="_blank">blaha@theochem.tuwien.ac.at</a>    WIEN2k: <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__www.wien2k.at&d=DwIDaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=l1d2pIyr6QG4zudFD9mIvqSRlEOa2QsiTeaYL3aqbw4&s=AoRHoP7PeTtPCMY9sk_zfa0dbR4Ia5Q5efrMBa2qMq4&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__www.wien2k.at&d=DwIDaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=l1d2pIyr6QG4zudFD9mIvqSRlEOa2QsiTeaYL3aqbw4&s=AoRHoP7PeTtPCMY9sk_zfa0dbR4Ia5Q5efrMBa2qMq4&e=</a><br>
>      WWW:<br>
>      <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__www.imc.tuwien.ac.at_tc-5Fblaha-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D&d=DwIDaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=l1d2pIyr6QG4zudFD9mIvqSRlEOa2QsiTeaYL3aqbw4&s=DSEHUVbTJi_6QfalAunhz2X7sOekjMlmFnCBfffZyEI&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__www.imc.tuwien.ac.at_tc-5Fblaha-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D&d=DwIDaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=l1d2pIyr6QG4zudFD9mIvqSRlEOa2QsiTeaYL3aqbw4&s=DSEHUVbTJi_6QfalAunhz2X7sOekjMlmFnCBfffZyEI&e=</a><br>
><br>
>      _______________________________________________<br>
>      Wien mailing list<br>
>      <a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
>      <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__zeus.theochem.tuwien.ac.at_mailman_listinfo_wien&d=DwIDaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=l1d2pIyr6QG4zudFD9mIvqSRlEOa2QsiTeaYL3aqbw4&s=HqjVLD2WpfrSwRKFbqf7V88ADhDECm1QctI1U57Y7xA&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__zeus.theochem.tuwien.ac.at_mailman_listinfo_wien&d=DwIDaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=l1d2pIyr6QG4zudFD9mIvqSRlEOa2QsiTeaYL3aqbw4&s=HqjVLD2WpfrSwRKFbqf7V88ADhDECm1QctI1U57Y7xA&e=</a><br>
>      SEARCH the MAILING-LIST at:  <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__www.mail-2Darchive.com_wien-40zeus.theochem.tuwien.ac.at_index.html&d=DwIDaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=l1d2pIyr6QG4zudFD9mIvqSRlEOa2QsiTeaYL3aqbw4&s=_FHptNqbf_hSe5ztGYQwgom849ky-4GtwrFrc8MOZhQ&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__www.mail-2Darchive.com_wien-40zeus.theochem.tuwien.ac.at_index.html&d=DwIDaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=l1d2pIyr6QG4zudFD9mIvqSRlEOa2QsiTeaYL3aqbw4&s=_FHptNqbf_hSe5ztGYQwgom849ky-4GtwrFrc8MOZhQ&e=</a><br>
><br>
><br>
></blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><span style="font-size:12.8px">Professor Laurence Marks</span><br></div><div dir="ltr"><span style="font-size:12.8px">"Research is to see what everybody else has seen, and to think what nobody else has thought", </span><span style="font-size:12.8px">Albert Szent-Gyorgi</span><br><a href="http://www.numis.northwestern.edu" target="_blank">www.numis.northwestern.edu</a> ; <span style="font-size:12.8px">Corrosion in 4D: </span><a href="http://MURI4D.numis.northwestern.edu" style="font-size:12.8px" target="_blank">MURI4D.numis.northwestern.edu</a><div><span style="font-size:12.8px">Partner of the CFW 100% program for gender equity, </span><a href="http://www.cfw.org/100-percent" style="font-size:12.8px" target="_blank">www.cfw.org/100-percent</a></div><div>Co-Editor, Acta Cryst A</div></div></div></div></div></div></div></div></div></div></div>