<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">
<HTML>
<HEAD>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=us-ascii">
<META NAME="Generator" CONTENT="MS Exchange Server version 08.00.0681.000">
<TITLE>CuO, LDA+U, Pseudopots, again</TITLE>
</HEAD>
<BODY>
<!-- Converted from text/rtf format -->

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">To All:</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">Please refer to a previous pw_forum thread &quot;CuO, LDA+U, Pseudopots,</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">etc.&quot;</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">I'm continuing to try to achieve a pwscf calculation which will yield a</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&quot;Hubbard gap&quot; in any compound containing Cu++ and O--.&nbsp; My efforts to</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">date have been unsucessful.&nbsp; My &quot;model&quot; has been a pseudo-rocksalt CuO</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">structure based on &quot;proxy&quot; lattice constants near those of FeO, NiO and</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">MgO.</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">I strongly suspect my problem involves use (or misuse or unavailability)</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">of a proper Cu pseudopotential.&nbsp; The issue of Cu++ bonded to O-- is</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">central to the physics of high temperature superconductivity, and a</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">failure of pwscf to properly resolve this difficulty would be</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">detrimential to the future deployment of quantum-espresso to this very</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">important and critical area of computational materials physics, not</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">just HTSC, but possibly GMR</FONT></SPAN><SPAN LANG="en-us"><FONT FACE="Calibri"> structures</FONT></SPAN><SPAN LANG="en-us"><FONT FACE="Calibri"> as well.</FONT></SPAN><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp; Note in what follows</FONT></SPAN><SPAN LANG="en-us"> <FONT FACE="Calibri">is that I changed occ_loc in tabd.f90 from 10 to 9.</FONT></SPAN><SPAN LANG="en-us"></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">I&nbsp; attach&nbsp; below&nbsp; 1)&nbsp; the&nbsp; list&nbsp; of&nbsp; PPs&nbsp; I've&nbsp; tried&nbsp; and their d-level</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">configurations.</FONT></SPAN><SPAN LANG="en-us">&nbsp;<FONT FACE="Calibri"> </FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">None&nbsp; really contain the&nbsp; generally assumed 3d9&nbsp; Hund's</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">Rule configuration for CuO complexes.&nbsp;&nbsp; Following this is the&nbsp; output 2)</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">of a &quot;hand written&nbsp; grep-er&quot; of the scf.out&nbsp; file for a typical&nbsp; example</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">(one of about 20) combinations of Cu PPs and the &quot;starting_ns_eigenvalue</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">suggestions&quot; to escape the &quot;minimum energy trap.&quot;&nbsp; I selected</FONT></SPAN><SPAN LANG="en-us"> <FONT FACE="Calibri">my</FONT></SPAN><SPAN LANG="en-us"><FONT FACE="Calibri"> &quot;best</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">attempt&quot; at this effort.&nbsp; The critical (IMHO) lines are pointed to by a</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&quot;--&gt;&quot; marker.&nbsp; BTW, I didn't try any Wannier projections since the</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">eigenvalue forcings didn't work, and my lesson from Example 25 is that</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">these have to be done</FONT></SPAN><SPAN LANG="en-us"> <FONT FACE="Calibri">successfully</FONT></SPAN><SPAN LANG="en-us"><FONT FACE="Calibri"> first.</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">Any of those in the pwscf community aware of a &quot;better&quot; PP for my</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">problem, please e-mail such to me.&nbsp; In the meantime, I'm experimenting</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&quot;rolling my own&quot; with the tools available at the Berlin-FHI site,</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">trying 3d9-4s2 as a first attempt--stay tuned.</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">Sorry for all my American slang.&nbsp; We didn't invent the language...the</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">Brits did</FONT></SPAN><SPAN LANG="en-us"><FONT FACE="Calibri">, and shoved it down our throats</FONT></SPAN><SPAN LANG="en-us"><FONT FACE="Calibri">.&nbsp; We just made English more user friendly.</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">-Paul</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">----------------------------------</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">1)</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">Cu.blyp-n-van_ak.UPF</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;nl pn&nbsp; l&nbsp;&nbsp; occ</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;3D&nbsp; 3&nbsp; 2&nbsp; 9.50</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;4S&nbsp; 4&nbsp; 0&nbsp; 1.00</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;4P&nbsp; 4&nbsp; 1&nbsp; 0.50</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">Cu.pbe-d-rrkjus.UPF</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;nl pn&nbsp; l&nbsp;&nbsp; occ</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;3D&nbsp; 3&nbsp; 2 10.00</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;3D&nbsp; 3&nbsp; 2&nbsp; 0.00</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;4P&nbsp; 2&nbsp; 1 -1.00</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;4S&nbsp; 1&nbsp; 0&nbsp; 1.00</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">Cu.pbe-n-van_ak.UPF</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;nl pn&nbsp; l&nbsp;&nbsp; occ</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;3D&nbsp; 3&nbsp; 2&nbsp; 9.50</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;4S&nbsp; 4&nbsp; 0&nbsp; 1.00</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;4P&nbsp; 4&nbsp; 1&nbsp; 0.50</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">Cu.pw91-n-van_ak.UPF</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;nl pn&nbsp; l&nbsp;&nbsp; occ</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;3D&nbsp; 3&nbsp; 2&nbsp; 9.50</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;4S&nbsp; 4&nbsp; 0&nbsp; 1.00</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;4P&nbsp; 4&nbsp; 1&nbsp; 0.50</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">Cu.pz-d-rrkjus.UPF</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;nl pn&nbsp; l&nbsp;&nbsp; occ</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;3D&nbsp; 3&nbsp; 2 10.00</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;3D&nbsp; 3&nbsp; 2&nbsp; 0.00</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;4P&nbsp; 2&nbsp; 1 -1.00</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;4S&nbsp; 1&nbsp; 0&nbsp; 1.00</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">Cu.pz-n-van_ak.UPF</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;nl pn&nbsp; l&nbsp;&nbsp; occ</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;3D&nbsp; 3&nbsp; 2&nbsp; 9.50</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;4S&nbsp; 4&nbsp; 0&nbsp; 1.00</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;4P&nbsp; 4&nbsp; 1&nbsp; 0.50</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">---------------------------------------</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">2)</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">Today is Wed Jan&nbsp; 9 19:07:17 PST 2008</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">Now is Wed Jan&nbsp; 9 19:10:07 PST 2008</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">SCF Run preparatory to post-processing</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">Script is ./CuO_3.scf.sh, Output is CuO_3.scf.data</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;&nbsp;&nbsp;&nbsp; kinetic-energy cutoff&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 30.0000&nbsp; Ry</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;&nbsp;&nbsp;&nbsp; charge density cutoff&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp; 240.0000&nbsp; Ry</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;&nbsp;&nbsp;&nbsp; convergence threshold&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.0E-06</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">--&gt;&nbsp;&nbsp;&nbsp; pwscf_exec&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = /home/pmpgrant/O-sesame/PW/pw.x</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;&nbsp;&nbsp; a_Ang&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 4.200</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;&nbsp;&nbsp; a_Bohr&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 7.93950</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;&nbsp;&nbsp; nbnd&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 20,</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;&nbsp;&nbsp; nelec&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = ,</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;&nbsp;&nbsp; occupations&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 'smearing',</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;&nbsp;&nbsp; smearing&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = gauss,</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;&nbsp;&nbsp; degauss&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.01,</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;&nbsp;&nbsp; la2F&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = .false.,</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;&nbsp;&nbsp; diagonalization&nbsp;&nbsp;&nbsp;&nbsp; = 'david',</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;&nbsp;&nbsp; starting_magnetization(1) = 0.0</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;&nbsp;&nbsp; starting_magnetization(2) = 0.5</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;&nbsp;&nbsp; starting_magnetization(3) = -0.5</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;&nbsp;&nbsp; nspin&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 2,</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;&nbsp;&nbsp; lda_plus_u&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = .true.,</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">--&gt;&nbsp;&nbsp;&nbsp; Hubbard_U(2)&nbsp;&nbsp;&nbsp; = 5.0</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">--&gt;&nbsp;&nbsp;&nbsp; Hubbard_U(3)&nbsp;&nbsp;&nbsp; = 5.0</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">--&gt;&nbsp;&nbsp;&nbsp; U_projection_type&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = norm-atomic</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">--&gt;&nbsp;&nbsp;&nbsp; starting_ns_eigenvalue(1,2,2) =&nbsp; 0.0</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;&nbsp;&nbsp; starting_ns_eigenvalue(2,2,2) =</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;&nbsp;&nbsp; starting_ns_eigenvalue(3,2,2) =</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;&nbsp;&nbsp; starting_ns_eigenvalue(4,2,2) =</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;&nbsp;&nbsp; starting_ns_eigenvalue(5,2,2) =</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">--&gt;&nbsp;&nbsp;&nbsp; starting_ns_eigenvalue(1,1,3) =&nbsp; 0.0</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;&nbsp;&nbsp; starting_ns_eigenvalue(2,1,3) =</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;&nbsp;&nbsp; starting_ns_eigenvalue(3,1,3) =</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;&nbsp;&nbsp; starting_ns_eigenvalue(4,1,3) =</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">&nbsp;&nbsp;&nbsp; starting_ns_eigenvalue(5,1,3) =</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">--&gt;&nbsp;&nbsp;&nbsp; TM_pseudopot&nbsp;&nbsp;&nbsp; = Cu.pz-d-rrkjus.UPF</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">--&gt;&nbsp;&nbsp;&nbsp; A_pseudopot&nbsp;&nbsp;&nbsp;&nbsp; = O.LDA.US.RRKJ3.UPF</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">Run_Type&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; npk&nbsp;&nbsp;&nbsp;&nbsp; Ef(eV)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Total E(Ry)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; CPU_Time&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Wall_Time</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">scf&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 8&nbsp;&nbsp;&nbsp; 11.1457&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -239.16751151&nbsp;&nbsp;&nbsp;&nbsp; 1m 8.04s&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2m50.40s</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">d-orbital occupations</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">138:eigenvalues:&nbsp; 1.0000000 1.0000000 1.0000000 1.0000000 1.0000000</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Wingdings">&#232;<FONT FACE="Courier New">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</FONT></FONT></SPAN><SPAN LANG="en-us"></SPAN><SPAN LANG="en-us"></SPAN><SPAN LANG="en-us"></SPAN><SPAN LANG="en-us"> <FONT FACE="Calibri">152:eigenvalues:&nbsp; 0.8000000 0.8000000 0.8000000 0.8000000 0.8000000</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Wingdings">&#232;<FONT FACE="Courier New">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</FONT></FONT></SPAN><SPAN LANG="en-us"> <FONT FACE="Calibri">167:eigenvalues:&nbsp; 0.8000000 0.8000000 0.8000000 0.8000000 0.8000000</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"></SPAN><SPAN LANG="en-us"><FONT FACE="Calibri">181:eigenvalues:&nbsp; 1.0000000 1.0000000 1.0000000 1.0000000 1.0000000</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">213:eigenvalues:&nbsp; 0.9976937 0.9976937 0.9984132 0.9984132 1.0008045</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">227:eigenvalues:&nbsp; 0.5164626 0.5164626 0.9680621 0.9823458 0.9823458</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">242:eigenvalues:&nbsp; 0.5164626 0.5164626 0.9680621 0.9823458 0.9823458</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">256:eigenvalues:&nbsp; 0.9976937 0.9976937 0.9984132 0.9984132 1.0008045</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">277:eigenvalues:&nbsp; 0.9976937 0.9976937 0.9984132 0.9984132 1.0008045</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">--&gt; 291:eigenvalues:&nbsp; 0.0000000 0.5164626 0.9680621 0.9823458 0.9823458</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">--&gt; 306:eigenvalues:&nbsp; 0.0000000 0.5164626 0.9680621 0.9823458 0.9823458</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">320:eigenvalues:&nbsp; 0.9976937 0.9976937 0.9984132 0.9984132 1.0008045</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">353:eigenvalues:&nbsp; 0.9936756 0.9936756 0.9968054 0.9968054 0.9991184</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">--&gt; 367:eigenvalues:&nbsp; 0.7870114 0.7870114 0.9938818 0.9938818 0.9969944</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">--&gt; 382:eigenvalues:&nbsp; 0.7867963 0.7867963 0.9938806 0.9938806 0.9969906</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">396:eigenvalues:&nbsp; 0.9936787 0.9936787 0.9968055 0.9968055 0.9991185</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">429:eigenvalues:&nbsp; 0.9914062 0.9914062 0.9971611 0.9971611 0.9994970</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">443:eigenvalues:&nbsp; 0.7968500 0.7968500 0.9944826 0.9944826 0.9975186</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">458:eigenvalues:&nbsp; 0.7971215 0.7971215 0.9944915 0.9944915 0.9975167</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">472:eigenvalues:&nbsp; 0.9914126 0.9914126 0.9971603 0.9971603 0.9994966</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">507:eigenvalues:&nbsp; 0.9667233 0.9667233 0.9969086 0.9969086 0.9992909</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">521:eigenvalues:&nbsp; 0.8723219 0.8723219 0.9959798 0.9959798 0.9984415</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">536:eigenvalues:&nbsp; 0.8713559 0.8713559 0.9959649 0.9959649 0.9984298</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">550:eigenvalues:&nbsp; 0.9663176 0.9663176 0.9968918 0.9968918 0.9992743</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">585:eigenvalues:&nbsp; 0.9525509 0.9525509 0.9963181 0.9963181 0.9986638</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">599:eigenvalues:&nbsp; 0.9127300 0.9127300 0.9961748 0.9961748 0.9985388</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">614:eigenvalues:&nbsp; 0.9146296 0.9146296 0.9962191 0.9962191 0.9985819</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">628:eigenvalues:&nbsp; 0.9537915 0.9537915 0.9963783 0.9963783 0.9987221</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">663:eigenvalues:&nbsp; 0.9458265 0.9458265 0.9962868 0.9962868 0.9986542</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">677:eigenvalues:&nbsp; 0.9104092 0.9104092 0.9959371 0.9959371 0.9983246</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">692:eigenvalues:&nbsp; 0.9047312 0.9047312 0.9957990 0.9957990 0.9982026</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">706:eigenvalues:&nbsp; 0.9417229 0.9417229 0.9961282 0.9961282 0.9984984</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">741:eigenvalues:&nbsp; 0.9354950 0.9354950 0.9960520 0.9960520 0.9984263</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">755:eigenvalues:&nbsp; 0.9174694 0.9174694 0.9958534 0.9958534 0.9982404</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">770:eigenvalues:&nbsp; 0.9179431 0.9179431 0.9958601 0.9958601 0.9982469</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">784:eigenvalues:&nbsp; 0.9358828 0.9358828 0.9960625 0.9960625 0.9984347</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">819:eigenvalues:&nbsp; 0.9341625 0.9341625 0.9960316 0.9960316 0.9984028</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">833:eigenvalues:&nbsp; 0.9222122 0.9222122 0.9959062 0.9959062 0.9982848</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">848:eigenvalues:&nbsp; 0.9216331 0.9216331 0.9958882 0.9958882 0.9982692</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">862:eigenvalues:&nbsp; 0.9336231 0.9336231 0.9960127 0.9960127 0.9983843</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">897:eigenvalues:&nbsp; 0.9327282 0.9327282 0.9960266 0.9960266 0.9983951</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">911:eigenvalues:&nbsp; 0.9284646 0.9284646 0.9959858 0.9959858 0.9983585</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">926:eigenvalues:&nbsp; 0.9283145 0.9283145 0.9959811 0.9959811 0.9983538</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">940:eigenvalues:&nbsp; 0.9325413 0.9325413 0.9960218 0.9960218 0.9983904</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">975:eigenvalues:&nbsp; 0.9318453 0.9318453 0.9960271 0.9960271 0.9983966</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">989:eigenvalues:&nbsp; 0.9297349 0.9297349 0.9960062 0.9960062 0.9983772</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">1004:eigenvalues:&nbsp; 0.9299729 0.9299729 0.9960136 0.9960136 0.9983841</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">1018:eigenvalues:&nbsp; 0.9320457 0.9320457 0.9960350 0.9960350 0.9984044</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">1053:eigenvalues:&nbsp; 0.9311035 0.9311035 0.9960273 0.9960273 0.9983963</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">1067:eigenvalues:&nbsp; 0.9316853 0.9316853 0.9960330 0.9960330 0.9984016</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">1082:eigenvalues:&nbsp; 0.9316968 0.9316968 0.9960336 0.9960336 0.9984023</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">1096:eigenvalues:&nbsp; 0.9310994 0.9310994 0.9960275 0.9960275 0.9983966</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">--&gt; 1131:eigenvalues:&nbsp; 0.9314650 0.9314650 0.9960338 0.9960338 0.9984027</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">--&gt; 1145:eigenvalues:&nbsp; 0.9312033 0.9312033 0.9960312 0.9960312 0.9984004</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">--&gt; 1160:eigenvalues:&nbsp; 0.9311989 0.9311989 0.9960309 0.9960309 0.9984001</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Calibri">--&gt; 1174:eigenvalues:&nbsp; 0.9314517 0.9314517 0.9960337 0.9960337 0.9984026</FONT></SPAN><SPAN LANG="en-us"></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><B><I></I></B></SPAN><SPAN LANG="en-us"><B><I></I></B></SPAN><SPAN LANG="en-us"><B><I></I></B></SPAN><B><I><SPAN LANG="en-us"><FONT SIZE=4 FACE="Comic Sans MS">Paul M. Grant, PhD</FONT></SPAN></I></B><SPAN LANG="en-us"></SPAN><SPAN LANG="en-us"></SPAN><SPAN LANG="en-us"></SPAN><SPAN LANG="en-us"></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><B></B></SPAN><SPAN LANG="en-us"><B></B></SPAN><SPAN LANG="en-us"><B></B></SPAN><B><SPAN LANG="en-us"><FONT SIZE=2 FACE="Arial">Principal, W2AGZ Technologies</FONT></SPAN></B><SPAN LANG="en-us"></SPAN><SPAN LANG="en-us"></SPAN><SPAN LANG="en-us"></SPAN><SPAN LANG="en-us"></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><B></B></SPAN><SPAN LANG="en-us"><B></B></SPAN><SPAN LANG="en-us"><B></B></SPAN><B><SPAN LANG="en-us"><FONT SIZE=2 FACE="Arial">Visiting Scholar, Applied Physics, Stanford University</FONT></SPAN></B><SPAN LANG="en-us"></SPAN><SPAN LANG="en-us"></SPAN><SPAN LANG="en-us"></SPAN><SPAN LANG="en-us"></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><B></B></SPAN><SPAN LANG="en-us"><B></B></SPAN><SPAN LANG="en-us"><B></B></SPAN><B><SPAN LANG="en-us"><FONT SIZE=2 FACE="Arial">EPRI Science Fellow (Retired)</FONT></SPAN></B><SPAN LANG="en-us"></SPAN><SPAN LANG="en-us"></SPAN><SPAN LANG="en-us"></SPAN><SPAN LANG="en-us"></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><B></B></SPAN><SPAN LANG="en-us"><B></B></SPAN><SPAN LANG="en-us"><B></B></SPAN><B><SPAN LANG="en-us"><FONT SIZE=2 FACE="Arial">IBM Research Staff Member Emeritus</FONT></SPAN></B><SPAN LANG="en-us"></SPAN><SPAN LANG="en-us"></SPAN><SPAN LANG="en-us"></SPAN><SPAN LANG="en-us"></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"></SPAN><A HREF="mailto:w2agz@pacbell.net"><SPAN LANG="en-us"><B></B></SPAN><SPAN LANG="en-us"><B></B></SPAN><SPAN LANG="en-us"><B></B></SPAN><B><SPAN LANG="en-us"><FONT SIZE=2 FACE="Arial">w2agz@pacbell.net</FONT></SPAN></B><SPAN LANG="en-us"></SPAN></A><SPAN LANG="en-us"></SPAN><SPAN LANG="en-us"></SPAN><SPAN LANG="en-us"></SPAN><SPAN LANG="en-us"></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"></SPAN><A HREF="http://www.w2agz.com"><SPAN LANG="en-us"><B></B></SPAN><SPAN LANG="en-us"><B></B></SPAN><SPAN LANG="en-us"><B></B></SPAN><B><SPAN LANG="en-us"><FONT SIZE=2 FACE="Arial">http://www.w2agz.com</FONT></SPAN></B><SPAN LANG="en-us"></SPAN></A><SPAN LANG="en-us"></SPAN><SPAN LANG="en-us"></SPAN><SPAN LANG="en-us"></SPAN><SPAN LANG="en-us"></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Times New Roman">&nbsp;</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"><FONT FACE="Times New Roman">&nbsp;</FONT></SPAN></P>

<P DIR=LTR><SPAN LANG="en-us"></SPAN></P>

</BODY>
</HTML>