<?xml version="1.0" encoding="UTF-8" ?>
<cml convention="CMLComp" xmlns="http://www.xml-cml.org/schema"
 xmlns:siesta="http://www.uam.es/siesta/namespace"
 xmlns:siestaUnits="http://www.uam.es/siesta/namespace/units"
 xmlns:xsd="http://www.w3.org/2001/XMLSchema"
 xmlns:fpx="http://www.uszla.me.uk/fpx"
 xmlns:dc="http://purl.org/dc/elements/1.1/"
 xmlns:units="http://www.uszla.me.uk/FoX/units"
 xmlns:cmlUnits="http://www.xml-cml.org/units/units"
 xmlns:siUnits="http://www.xml-cml.org/units/siUnits"
 xmlns:atomicUnits="http://www.xml-cml.org/units/atomic">
 <metadataList>
  <metadata name="siesta:Program" content="Siesta" />
  <metadata name="siesta:Version" content="v4.1-b4" />
  <metadata name="siesta:Arch" content="unknown" />
  <metadata name="siesta:Flags" content="mpif90 -O2 -fPIC -ftree-vectorize" />
  <metadata name="siesta:PPFlags"
   content="-DFC_HAVE_ABORT -DMPI -DCDF -DNCDF -DNCDF_4 -DSIESTA__DIAG_2STAGE" />
  <metadata name="siesta:StartTime" content="2023-01-19T14-03-52" />
  <metadata name="siesta:run_UUID"
   content="8a55d810-980a-11ed-4dd5-b42addee316e" />
  <metadata name="siesta:Mode" content="Parallel" />
  <metadata name="siesta:Nodes" content="        12" />
  <metadata name="siesta:NetCDF" content="true" />
 </metadataList>
 <module title="Initial System">
  <molecule>
   <atomArray>
    <atom elementType="H" id="a1" ref="siesta:e002"
     x3="    0.54909300                               "
     y3="    2.15544000                               "
     z3="   -0.33830000                               " />
    <atom elementType="P" id="a2" ref="siesta:e003"
     x3="    0.00004600                               "
     y3="    0.00000300                               "
     z3="   -0.00004600                               " />
    <atom elementType="O" id="a3" ref="siesta:e001"
     x3="   -1.52650300                               "
     y3="   -0.44109200                               "
     z3="   -0.28837100                               " />
    <atom elementType="H" id="a4" ref="siesta:e002"
     x3="   -1.72075800                               "
     y3="   -1.26327900                               "
     z3="    0.22747600                               " />
    <atom elementType="O" id="a5" ref="siesta:e001"
     x3="    0.63508400                               "
     y3="   -0.27454000                               "
     z3="    1.29608100                               " />
    <atom elementType="O" id="a6" ref="siesta:e001"
     x3="    0.84483000                               "
     y3="   -0.58810600                               "
     z3="   -1.24864600                               " />
    <atom elementType="O" id="a7" ref="siesta:e001"
     x3="   -0.25213100                               "
     y3="    1.55933500                               "
     z3="   -0.39017200                               " />
    <atom elementType="H" id="a8" ref="siesta:e002"
     x3="    0.37803100                               "
     y3="   -0.63081200                               "
     z3="   -2.13200000                               " />
   </atomArray>
  </molecule>
  <lattice dictRef="siesta:ucell">
   <latticeVector units="siestaUnits:angstrom" dictRef="cml:latticeVector">
  2.194225476387E+01  0.000000000000E+00  0.000000000000E+00

   </latticeVector>
   <latticeVector units="siestaUnits:angstrom" dictRef="cml:latticeVector">
  0.000000000000E+00  2.147685056073E+01  0.000000000000E+00

   </latticeVector>
   <latticeVector units="siestaUnits:angstrom" dictRef="cml:latticeVector">
  0.000000000000E+00  0.000000000000E+00  2.230268879633E+01

   </latticeVector>
  </lattice>
  <property dictRef="siesta:shape">
   <scalar>molecule
   </scalar>
  </property>
 </module>
 <parameterList title="Input Parameters">
  <parameter name="SystemName" dictRef="siesta:sname">
   <scalar dataType="xsd:string">H3PO4 molecule
   </scalar>
  </parameter>
  <parameter name="SystemLabel" dictRef="siesta:slabel">
   <scalar dataType="xsd:string">h3po4_3
   </scalar>
  </parameter>
  <parameter name="LongOutput" dictRef="siesta:verbosity">
   <scalar dataType="xsd:boolean">false
   </scalar>
  </parameter>
  <parameter title="NumberOfSpecies" dictRef="siesta:ns">
   <scalar dataType="xsd:integer" units="cmlUnits:countable">3
   </scalar>
  </parameter>
  <parameter name="WriteDenChar">
   <scalar dataType="xsd:boolean">false
   </scalar>
  </parameter>
  <parameter name="WriteMullikenPop">
   <scalar dataType="xsd:integer" units="cmlUnits:dimensionless">0
   </scalar>
  </parameter>
  <parameter name="MatelNRTAB" dictRef="siesta:matel_nrtab">
   <scalar dataType="xsd:integer" units="cmlUnits:countable">1024
   </scalar>
  </parameter>
  <parameter name="MeshCutOff" dictRef="siesta:g2max">
   <scalar dataType="xsd:double" units="siestaUnits:Ry">300.000000000                                                                                       
   </scalar>
  </parameter>
  <parameter name="NetCharge" dictRef="siesta:NetCharge">
   <scalar dataType="xsd:double" units="siestaUnits:e__">0.00000000000                                                                                       
   </scalar>
  </parameter>
  <parameter name="MaxSCFIterations" dictRef="siesta:maxscf">
   <scalar dataType="xsd:integer" units="cmlUnits:countable">1000
   </scalar>
  </parameter>
  <parameter name="MinSCFIterations" dictRef="siesta:minscf">
   <scalar dataType="xsd:integer" units="cmlUnits:countable">0
   </scalar>
  </parameter>
  <parameter name="DM.NumberPulay" dictRef="siesta:maxsav">
   <scalar dataType="xsd:integer" units="cmlUnits:countable">7
   </scalar>
  </parameter>
  <parameter name="DM.NumberBroyden" dictRef="siesta:broyden_maxit">
   <scalar dataType="xsd:integer" units="cmlUnits:countable">0
   </scalar>
  </parameter>
  <parameter name="DM.MixSCF1" dictRef="siesta:mix">
   <scalar dataType="xsd:boolean">true
   </scalar>
  </parameter>
  <parameter name="DM.PulayOnFile" dictRef="siesta:pulfile">
   <scalar dataType="xsd:boolean">false
   </scalar>
  </parameter>
  <parameter name="DM.MixingWeight" dictRef="siesta:wmix">
   <scalar dataType="xsd:double" units="cmlUnits:dimensionless">0.300000000000                                                                                      
   </scalar>
  </parameter>
  <parameter name="DM.OccupancyTolerance" dictRef="siesta:occtol">
   <scalar dataType="xsd:double" units="cmlUnits:dimensionless">
   0.100000000000E-11                                                                                  
   </scalar>
  </parameter>
  <parameter name="DM.NumberKick" dictRef="siesta:nkick">
   <scalar dataType="xsd:integer" units="cmlUnits:countable">0
   </scalar>
  </parameter>
  <parameter name="DM.KickMixingWeight" dictRef="siesta:wmixkick">
   <scalar dataType="xsd:double" units="cmlUnits:dimensionless">0.500000000000                                                                                      
   </scalar>
  </parameter>
  <parameter name="SCF.Harris.Converge" dictRef="siesta:ReqHarrisConv">
   <scalar dataType="xsd:boolean">false
   </scalar>
  </parameter>
  <parameter name="SCF.Harris.Tolerance" dictRef="siesta:Harris_tolerance">
   <scalar dataType="xsd:double" units="siestaUnits:eV">0.100000000000E-03                                                                                  
   </scalar>
  </parameter>
  <parameter name="SCF.DM.Converge" dictRef="siesta:ReqDMConv">
   <scalar dataType="xsd:boolean">true
   </scalar>
  </parameter>
  <parameter name="SCF.DM.Tolerance" dictRef="siesta:dDtol">
   <scalar dataType="xsd:double" units="siestaUnits:eAng_3">0.100000000000E-03                                                                                  
   </scalar>
  </parameter>
  <parameter name="SCF.EDM.Converge" dictRef="siesta:ReqEDMConv">
   <scalar dataType="xsd:boolean">false
   </scalar>
  </parameter>
  <parameter name="SCF.EDM.Tolerance" dictRef="siesta:EDM_tolerance">
   <scalar dataType="xsd:double" units="siestaUnits:eVeAng_3">
   0.100000000000E-02                                                                                  
   </scalar>
  </parameter>
  <parameter name="SCF.H.Converge" dictRef="siesta:ReqHConv">
   <scalar dataType="xsd:boolean">true
   </scalar>
  </parameter>
  <parameter name="SCF.H.Tolerance" dictRef="siesta:dHtol">
   <scalar dataType="xsd:double" units="siestaUnits:eV">0.100000000000E-02                                                                                  
   </scalar>
  </parameter>
  <parameter name="SCF.FreeE.Converge" dictRef="siesta:ReqEnergyConv">
   <scalar dataType="xsd:boolean">false
   </scalar>
  </parameter>
  <parameter name="SCF.FreeE.Tolerance" dictRef="siesta:dEtol">
   <scalar dataType="xsd:double" units="siestaUnits:eV">0.100000000000E-03                                                                                  
   </scalar>
  </parameter>
  <parameter name="DM.UseSaveDM" dictRef="siesta:usesavedm">
   <scalar dataType="xsd:boolean">false
   </scalar>
  </parameter>
  <parameter name="NeglNonOverlapInt" dictRef="siesta:negl">
   <scalar dataType="xsd:boolean">false
   </scalar>
  </parameter>
  <parameter name="SolutionMethod" dictRef="siesta:SCFmethod">
   <scalar dataType="xsd:string">diagon
   </scalar>
  </parameter>
  <parameter name="ElectronicTemperature" dictRef="siesta:etemp">
   <scalar dataType="xsd:double" units="siestaUnits:Ry">0.190000000000E-02                                                                                  
   </scalar>
  </parameter>
  <parameter name="FixSpin" dictRef="siesta:fixspin">
   <scalar dataType="xsd:boolean">false
   </scalar>
  </parameter>
  <parameter name="TotalSpin" dictRef="siesta:totalspin">
   <scalar dataType="xsd:double" units="siestaUnits:eSpin">0.00000000000                                                                                       
   </scalar>
  </parameter>
  <parameter name="MD.TypeOfRun">
   <scalar dataType="xsd:string">Single-Point
   </scalar>
  </parameter>
  <parameter name="MD.BulkModulus">
   <scalar dataType="xsd:double" units="siestaUnits:Ry_Bohr__3">
   0.679773000000E-02                                                                                  
   </scalar>
  </parameter>
 </parameterList>
 <propertyList title="k-points" dictRef="siesta:kpoints">
  <property dictRef="siesta:nkpnt">
   <scalar dataType="xsd:integer" units="cmlUnits:countable">1
   </scalar>
  </property>
  <kpoint coords="  0.00000000  0.00000000  0.00000000" weight="  1.00000000" />
  <property dictRef="siesta:kcutof">
   <scalar dataType="xsd:double" units="siestaUnits:angstrom">5.68252767459                                                                                       
   </scalar>
  </property>
 </propertyList>
 <property dictRef="siesta:kscell">
 <!--In matrix, row (first) index is fastest-->
  <matrix units="siestaUnits:Ang" columns="3" rows="3" dataType="xsd:integer">
           1           0           0
           0           1           0
           0           0           1

  </matrix>
 </property>
 <property dictRef="siesta:kdispl">
  <array units="siestaUnits:Ang" size="3" dataType="xsd:double">
  0.000000000000E+00  0.000000000000E+00  0.000000000000E+00

  </array>
 </property>
 <module dictRef="Single-Point" role="step" serial="1">
  <molecule>
   <atomArray>
    <atom elementType="H" id="a1" ref="siesta:e002"
     x3="    0.54909300                               "
     y3="    2.15544000                               "
     z3="   -0.33830000                               " />
    <atom elementType="P" id="a2" ref="siesta:e003"
     x3="    0.00004600                               "
     y3="    0.00000300                               "
     z3="   -0.00004600                               " />
    <atom elementType="O" id="a3" ref="siesta:e001"
     x3="   -1.52650300                               "
     y3="   -0.44109200                               "
     z3="   -0.28837100                               " />
    <atom elementType="H" id="a4" ref="siesta:e002"
     x3="   -1.72075800                               "
     y3="   -1.26327900                               "
     z3="    0.22747600                               " />
    <atom elementType="O" id="a5" ref="siesta:e001"
     x3="    0.63508400                               "
     y3="   -0.27454000                               "
     z3="    1.29608100                               " />
    <atom elementType="O" id="a6" ref="siesta:e001"
     x3="    0.84483000                               "
     y3="   -0.58810600                               "
     z3="   -1.24864600                               " />
    <atom elementType="O" id="a7" ref="siesta:e001"
     x3="   -0.25213100                               "
     y3="    1.55933500                               "
     z3="   -0.39017200                               " />
    <atom elementType="H" id="a8" ref="siesta:e002"
     x3="    0.37803100                               "
     y3="   -0.63081200                               "
     z3="   -2.13200000                               " />
   </atomArray>
  </molecule>
  <lattice dictRef="siesta:ucell">
   <latticeVector units="siestaUnits:Ang" dictRef="cml:latticeVector">
  1.161133654918E+01  0.000000000000E+00  0.000000000000E+00

   </latticeVector>
   <latticeVector units="siestaUnits:Ang" dictRef="cml:latticeVector">
  0.000000000000E+00  1.136505534918E+01  0.000000000000E+00

   </latticeVector>
   <latticeVector units="siestaUnits:Ang" dictRef="cml:latticeVector">
  0.000000000000E+00  0.000000000000E+00  1.180206994918E+01

   </latticeVector>
  </lattice>
  <crystal title="Lattice Parameters">
   <scalar title="a" dictRef="cml:a" units="units:angstrom                ">
   11.611337                                                                                           
   </scalar>
   <scalar title="b" dictRef="cml:b" units="units:angstrom                ">
   11.365055                                                                                           
   </scalar>
   <scalar title="c" dictRef="cml:c" units="units:angstrom                ">
   11.802070                                                                                           
   </scalar>
   <scalar title="alpha" dictRef="cml:alpha"
    units="units:degree                  ">90.000000                                                                                           
   </scalar>
   <scalar title="beta" dictRef="cml:beta"
    units="units:degree                  ">90.000000                                                                                           
   </scalar>
   <scalar title="gamma" dictRef="cml:gamma"
    units="units:degree                  ">90.000000                                                                                           
   </scalar>
  </crystal>
  <propertyList title="Orbital info">
   <property title="Number of orbitals in unit cell" dictRef="siesta:no_u">
    <scalar dataType="xsd:integer" units="cmlUnits:countable">76
    </scalar>
   </property>
   <property title="Number of non-zeros" dictRef="siesta:nnz">
    <scalar dataType="xsd:integer" units="cmlUnits:countable">5698
    </scalar>
   </property>
  </propertyList>
  <propertyList>
   <property title="Mesh" dictRef="siesta:ntm">
    <array units="cmlUnits:countable" size="3" dataType="xsd:integer">
         128         120         128

    </array>
   </property>
   <property title="Requested Cut-Off" dictRef="siesta:g2max">
    <scalar dataType="xsd:double" units="siestaUnits:Ry">300.000000000                                                                                       
    </scalar>
   </property>
   <property title="Actual Cut-Off" dictRef="siesta:g2mesh">
    <scalar dataType="xsd:double" units="siestaUnits:Ry">308.120821401                                                                                       
    </scalar>
   </property>
  </propertyList>
  <module dictRef="SCF" role="step" serial="1">
   <propertyList title="Energy Decomposition">
    <property dictRef="siesta:Ebs">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-364.425086                                                                                         
     </scalar>
    </property>
    <property dictRef="siesta:Eions">
     <scalar dataType="xsd:double" units="siestaUnits:eV">3463.417331                                                                                         
     </scalar>
    </property>
    <property dictRef="siesta:Ena">
     <scalar dataType="xsd:double" units="siestaUnits:eV">761.570620                                                                                          
     </scalar>
    </property>
    <property dictRef="siesta:Ekin">
     <scalar dataType="xsd:double" units="siestaUnits:eV">1553.234838                                                                                         
     </scalar>
    </property>
    <property dictRef="siesta:Enl">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-258.470844                                                                                         
     </scalar>
    </property>
    <property dictRef="siesta:Eldau">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.000000                                                                                            
     </scalar>
    </property>
    <property dictRef="siesta:DEna">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-100.726560                                                                                         
     </scalar>
    </property>
    <property dictRef="siesta:Eso">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.000000                                                                                            
     </scalar>
    </property>
    <property dictRef="siesta:DUscf">
     <scalar dataType="xsd:double" units="siestaUnits:eV">27.828765                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:DUext">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.000000                                                                                            
     </scalar>
    </property>
    <property dictRef="siesta:Enegf">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.000000                                                                                            
     </scalar>
    </property>
    <property dictRef="siesta:Exc">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-474.381342                                                                                         
     </scalar>
    </property>
    <property dictRef="siesta:Ecorrec">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.000000                                                                                            
     </scalar>
    </property>
    <property dictRef="siesta:Emad">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.000000                                                                                            
     </scalar>
    </property>
    <property dictRef="siesta:Emeta">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.000000                                                                                            
     </scalar>
    </property>
    <property dictRef="siesta:Emm">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.000000                                                                                            
     </scalar>
    </property>
    <property dictRef="siesta:Ekinion">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.000000                                                                                            
     </scalar>
    </property>
    <property dictRef="siesta:EharrsK">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1922.131498                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:EtotK">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1954.361854                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:FreeEK">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1954.361854                                                                                        
     </scalar>
    </property>
   </propertyList>
   <propertyList title="SCF Cycle">
    <property dictRef="siesta:Eharrs">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1922.1314980                                                                                       
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">1.7810599                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">13.7968903                                                                                          
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1954.3618543                                                                                       
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1954.3618543                                                                                       
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.0744253                                                                                           
     </scalar>
    </property>
   </propertyList>
  </module>
  <module dictRef="SCF" role="step" serial="2">
   <propertyList title="SCF Cycle">
    <property dictRef="siesta:Eharrs">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1973.5631643                                                                                       
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.6214498                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">5.2599081                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1966.3070651                                                                                       
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1966.3070651                                                                                       
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-0.3548012                                                                                          
     </scalar>
    </property>
   </propertyList>
  </module>
  <module dictRef="SCF" role="step" serial="3">
   <propertyList title="SCF Cycle">
    <property dictRef="siesta:Eharrs">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1967.5225988                                                                                       
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.1945461                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">1.2054452                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1967.1416013                                                                                       
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1967.1416013                                                                                       
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-0.8246960                                                                                          
     </scalar>
    </property>
   </propertyList>
  </module>
  <module dictRef="SCF" role="step" serial="4">
   <propertyList title="SCF Cycle">
    <property dictRef="siesta:Eharrs">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1967.2116592                                                                                       
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0493069                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.3523307                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1967.1832046                                                                                       
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1967.1832046                                                                                       
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1.0162871                                                                                          
     </scalar>
    </property>
   </propertyList>
  </module>
  <module dictRef="SCF" role="step" serial="5">
   <propertyList title="SCF Cycle">
    <property dictRef="siesta:Eharrs">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1967.1878455                                                                                       
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0293362                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.2150505                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1967.1882999                                                                                       
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1967.1882999                                                                                       
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1.0703379                                                                                          
     </scalar>
    </property>
   </propertyList>
  </module>
  <module dictRef="SCF" role="step" serial="6">
   <propertyList title="SCF Cycle">
    <property dictRef="siesta:Eharrs">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1967.1909322                                                                                       
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0078320                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.1604120                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1967.1897013                                                                                       
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1967.1897013                                                                                       
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1.0673109                                                                                          
     </scalar>
    </property>
   </propertyList>
  </module>
  <module dictRef="SCF" role="step" serial="7">
   <propertyList title="SCF Cycle">
    <property dictRef="siesta:Eharrs">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1967.1898911                                                                                       
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0135138                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.0864476                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1967.1905961                                                                                       
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1967.1905961                                                                                       
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1.0196931                                                                                          
     </scalar>
    </property>
   </propertyList>
  </module>
  <module dictRef="SCF" role="step" serial="8">
   <propertyList title="SCF Cycle">
    <property dictRef="siesta:Eharrs">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1967.1909499                                                                                       
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0027992                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.0209199                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1967.1908200                                                                                       
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1967.1908200                                                                                       
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1.0268404                                                                                          
     </scalar>
    </property>
   </propertyList>
  </module>
  <module dictRef="SCF" role="step" serial="9">
   <propertyList title="SCF Cycle">
    <property dictRef="siesta:Eharrs">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1967.1908411                                                                                       
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0008653                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.0008181                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1967.1908359                                                                                       
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1967.1908359                                                                                       
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1.0280532                                                                                          
     </scalar>
    </property>
   </propertyList>
  </module>
  <module dictRef="SCF" role="step" serial="10">
   <propertyList title="SCF Cycle">
    <property dictRef="siesta:Eharrs">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1967.1908361                                                                                       
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0000412                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.0002592                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1967.1908360                                                                                       
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1967.1908360                                                                                       
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1.0276833                                                                                          
     </scalar>
    </property>
   </propertyList>
  </module>
  <module title="SCF Finalization">
   <propertyList title="Energies and spin">
    <property dictRef="siesta:E_KS">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1967.190836                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1967.190836                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:Ebs">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-462.173982                                                                                         
     </scalar>
    </property>
    <property dictRef="siesta:E_Fermi">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1.027683                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:E_KS_egg">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-1967.190836                                                                                        
     </scalar>
    </property>
   </propertyList>
   <propertyList title="Forces">
   <!--Output: matrix fa(1:3,1:na_u)-->
    <property dictRef="siesta:forces">
    <!--In matrix, row (first) index is fastest-->
     <matrix units="siestaUnits:evpa" columns="8" rows="3" dataType="xsd:double">
 -6.197355209508E-01 -8.429398962120E-01  6.435049422438E-02
 -3.864484040841E-01  1.453251741970E+00 -2.024776678133E+00
 -1.005664472089E+00 -9.308548778506E-01  2.628911965965E-01
 -5.511234922795E-02  2.917646843209E-01 -1.919077862671E-01
  1.060274763758E+00 -8.014668316046E-01  2.418642537562E+00
  4.588450017855E-01 -7.019179044139E-01 -1.311349479598E+00
  3.815107595949E-01  1.326788463031E+00 -1.680143136781E-01
  1.590633543660E-01  2.132694335046E-01  9.530990777172E-01

     </matrix>
    </property>
    <property dictRef="siesta:ftot">
     <array units="siestaUnits:evpa" size="3" dataType="xsd:double">
 -7.266866848267E-03  7.894812745607E-03  2.935048424611E-03

     </array>
    </property>
    <property dictRef="siesta:fmax">
     <scalar dataType="xsd:double" units="siestaUnits:evpa">2.41864253756                                                                                       
     </scalar>
    </property>
    <property dictRef="siesta:fres">
     <scalar dataType="xsd:double" units="siestaUnits:evpa">0.967907507748                                                                                      
     </scalar>
    </property>
    <property dictRef="siesta:cfmax">
     <scalar dataType="xsd:double" units="siestaUnits:evpa">2.41864253756                                                                                       
     </scalar>
    </property>
   </propertyList>
   <property title="Stress" dictRef="siesta:stress">
   <!--In matrix, row (first) index is fastest-->
    <matrix units="siestaUnits:evpa3" columns="3" rows="3" dataType="xsd:double">
 -1.086147286472E-04  4.232501517033E-05 -3.772480418396E-05
  4.232617963461E-05 -3.632228950255E-05  2.371491685230E-05
 -3.772094410655E-05  2.371501445493E-05 -1.239436682304E-04

    </matrix>
   </property>
  </module>
 </module>
 <module title="Finalization">
  <molecule>
   <atomArray>
    <atom elementType="H" id="a1" ref="siesta:e002"
     x3="    0.54909300                               "
     y3="    2.15544000                               "
     z3="   -0.33830000                               " />
    <atom elementType="P" id="a2" ref="siesta:e003"
     x3="    0.00004600                               "
     y3="    0.00000300                               "
     z3="   -0.00004600                               " />
    <atom elementType="O" id="a3" ref="siesta:e001"
     x3="   -1.52650300                               "
     y3="   -0.44109200                               "
     z3="   -0.28837100                               " />
    <atom elementType="H" id="a4" ref="siesta:e002"
     x3="   -1.72075800                               "
     y3="   -1.26327900                               "
     z3="    0.22747600                               " />
    <atom elementType="O" id="a5" ref="siesta:e001"
     x3="    0.63508400                               "
     y3="   -0.27454000                               "
     z3="    1.29608100                               " />
    <atom elementType="O" id="a6" ref="siesta:e001"
     x3="    0.84483000                               "
     y3="   -0.58810600                               "
     z3="   -1.24864600                               " />
    <atom elementType="O" id="a7" ref="siesta:e001"
     x3="   -0.25213100                               "
     y3="    1.55933500                               "
     z3="   -0.39017200                               " />
    <atom elementType="H" id="a8" ref="siesta:e002"
     x3="    0.37803100                               "
     y3="   -0.63081200                               "
     z3="   -2.13200000                               " />
   </atomArray>
  </molecule>
  <lattice dictRef="siesta:ucell">
   <latticeVector units="siestaUnits:Ang" dictRef="cml:latticeVector">
  1.161133654918E+01  0.000000000000E+00  0.000000000000E+00

   </latticeVector>
   <latticeVector units="siestaUnits:Ang" dictRef="cml:latticeVector">
  0.000000000000E+00  1.136505534918E+01  0.000000000000E+00

   </latticeVector>
   <latticeVector units="siestaUnits:Ang" dictRef="cml:latticeVector">
  0.000000000000E+00  0.000000000000E+00  1.180206994918E+01

   </latticeVector>
  </lattice>
  <propertyList title="Eigenvalues">
   <property title="Fermi Energy" dictRef="siesta:E_Fermi">
    <scalar dataType="xsd:double" units="siestaUnits:ev">-1.02768                                                                                            
    </scalar>
   </property>
   <property title="Number of k-points" dictRef="siesta:nkpoints">
    <scalar dataType="xsd:integer" units="cmlUnits:countable">1
    </scalar>
   </property>
   <propertyList dictRef="siesta:kpt_band">
    <kpoint coords="  0.00000000  0.00000000  0.00000000" weight="  1.00000000" />
    <property dictRef="siesta:eigenenergies">
     <array units="siestaUnits:ev" size="76" dataType="xsd:double">
 -2.777462140461E+01 -2.591110919115E+01 -2.559744323857E+01 -2.433585471645E+01
 -1.661174588683E+01 -1.483040395184E+01 -1.332885262942E+01 -1.151401642132E+01
 -1.074829508212E+01 -1.067669862614E+01 -9.500819120638E+00 -9.034030511039E+00
 -8.674439390167E+00 -8.042713204679E+00 -7.377971820068E+00 -7.129300212161E+00
  7.224624006649E-01  1.626894967404E+00  2.078852254979E+00  2.842156354742E+00
  3.494378467279E+00  3.862987472967E+00  5.194303232580E+00  1.104832725744E+01
  1.478000580064E+01  1.529582090972E+01  1.717440023939E+01  1.726129104253E+01
  1.735571698469E+01  1.778009093590E+01  1.800244003241E+01  1.841124067812E+01
  1.987672062553E+01  2.053270957090E+01  2.072399166379E+01  2.089295505579E+01
  2.173016221040E+01  2.183452098193E+01  2.240250011807E+01  2.273964139257E+01
  2.309040468050E+01  2.358766905098E+01  2.379917739079E+01  2.432558415022E+01
  2.454382855781E+01  2.495671429778E+01  2.615085195679E+01  2.811598632027E+01
  2.865293381965E+01  2.932730969402E+01  2.954513625048E+01  3.176067928586E+01
  3.191889210902E+01  3.341590250585E+01  3.462841584067E+01  3.533177806681E+01
  3.581847058461E+01  3.741763066269E+01  3.774306503987E+01  3.825140693088E+01
  3.943074182312E+01  4.000378266064E+01  4.169328099672E+01  4.392496892931E+01
  4.562307978374E+01  4.604453350680E+01  4.712554717436E+01  4.883956202000E+01
  4.998373571798E+01  5.178296941807E+01  5.348666716424E+01  5.907949703413E+01
  6.937758289032E+01  7.154627929805E+01  7.692831774968E+01  7.948490385959E+01

     </array>
    </property>
   </propertyList>
  </propertyList>
  <propertyList title="Energy Decomposition">
   <property dictRef="siesta:Ebs">
    <scalar dataType="xsd:double" units="siestaUnits:eV">-462.173982                                                                                         
    </scalar>
   </property>
   <property dictRef="siesta:Eions">
    <scalar dataType="xsd:double" units="siestaUnits:eV">3463.417331                                                                                         
    </scalar>
   </property>
   <property dictRef="siesta:Ena">
    <scalar dataType="xsd:double" units="siestaUnits:eV">761.570620                                                                                          
    </scalar>
   </property>
   <property dictRef="siesta:Ekin">
    <scalar dataType="xsd:double" units="siestaUnits:eV">1396.970955                                                                                         
    </scalar>
   </property>
   <property dictRef="siesta:Enl">
    <scalar dataType="xsd:double" units="siestaUnits:eV">-203.813078                                                                                         
    </scalar>
   </property>
   <property dictRef="siesta:Eldau">
    <scalar dataType="xsd:double" units="siestaUnits:eV">0.000000                                                                                            
    </scalar>
   </property>
   <property dictRef="siesta:DEna">
    <scalar dataType="xsd:double" units="siestaUnits:eV">-13.917366                                                                                          
    </scalar>
   </property>
   <property dictRef="siesta:Eso">
    <scalar dataType="xsd:double" units="siestaUnits:eV">0.000000                                                                                            
    </scalar>
   </property>
   <property dictRef="siesta:DUscf">
    <scalar dataType="xsd:double" units="siestaUnits:eV">6.297906                                                                                            
    </scalar>
   </property>
   <property dictRef="siesta:DUext">
    <scalar dataType="xsd:double" units="siestaUnits:eV">0.000000                                                                                            
    </scalar>
   </property>
   <property dictRef="siesta:Enegf">
    <scalar dataType="xsd:double" units="siestaUnits:eV">0.000000                                                                                            
    </scalar>
   </property>
   <property dictRef="siesta:Exc">
    <scalar dataType="xsd:double" units="siestaUnits:eV">-450.882541                                                                                         
    </scalar>
   </property>
   <property dictRef="siesta:Ecorrec">
    <scalar dataType="xsd:double" units="siestaUnits:eV">0.000000                                                                                            
    </scalar>
   </property>
   <property dictRef="siesta:Emad">
    <scalar dataType="xsd:double" units="siestaUnits:eV">0.000000                                                                                            
    </scalar>
   </property>
   <property dictRef="siesta:Emeta">
    <scalar dataType="xsd:double" units="siestaUnits:eV">0.000000                                                                                            
    </scalar>
   </property>
   <property dictRef="siesta:Emm">
    <scalar dataType="xsd:double" units="siestaUnits:eV">0.000000                                                                                            
    </scalar>
   </property>
   <property dictRef="siesta:Ekinion">
    <scalar dataType="xsd:double" units="siestaUnits:eV">0.000000                                                                                            
    </scalar>
   </property>
   <property dictRef="siesta:EharrsK">
    <scalar dataType="xsd:double" units="siestaUnits:eV">-1967.190836                                                                                        
    </scalar>
   </property>
   <property dictRef="siesta:EtotK">
    <scalar dataType="xsd:double" units="siestaUnits:eV">-1967.190836                                                                                        
    </scalar>
   </property>
   <property dictRef="siesta:FreeEK">
    <scalar dataType="xsd:double" units="siestaUnits:eV">-1967.190836                                                                                        
    </scalar>
   </property>
  </propertyList>
  <propertyList title="Final Energy">
   <property dictRef="siesta:Ebs">
    <scalar dataType="xsd:double" units="siestaUnits:eV">-462.1739817                                                                                        
    </scalar>
   </property>
   <property dictRef="siesta:Ekin">
    <scalar dataType="xsd:double" units="siestaUnits:eV">1396.9709553                                                                                        
    </scalar>
   </property>
   <property dictRef="siesta:Uscf">
    <scalar dataType="xsd:double" units="siestaUnits:eV">2836.3395191                                                                                        
    </scalar>
   </property>
   <property dictRef="siesta:Eldau">
    <scalar dataType="xsd:double" units="siestaUnits:eV">0.0000000                                                                                           
    </scalar>
   </property>
   <property dictRef="siesta:Eso">
    <scalar dataType="xsd:double" units="siestaUnits:eV">0.0000000                                                                                           
    </scalar>
   </property>
   <property dictRef="siesta:DUext">
    <scalar dataType="xsd:double" units="siestaUnits:eV">0.0000000                                                                                           
    </scalar>
   </property>
   <property dictRef="siesta:Enegf">
    <scalar dataType="xsd:double" units="siestaUnits:eV">0.0000000                                                                                           
    </scalar>
   </property>
   <property dictRef="siesta:Exc">
    <scalar dataType="xsd:double" units="siestaUnits:eV">-450.8825414                                                                                        
    </scalar>
   </property>
   <property dictRef="siesta:I-e">
    <scalar dataType="xsd:double" units="siestaUnits:eV">-6941.1048237                                                                                       
    </scalar>
   </property>
   <property dictRef="siesta:I-I">
    <scalar dataType="xsd:double" units="siestaUnits:eV">1191.4860548                                                                                        
    </scalar>
   </property>
   <property dictRef="siesta:Ekinion">
    <scalar dataType="xsd:double" units="siestaUnits:eV">0.0000000                                                                                           
    </scalar>
   </property>
   <property dictRef="siesta:Etot">
    <scalar dataType="xsd:double" units="siestaUnits:eV">-1967.1908360                                                                                       
    </scalar>
   </property>
  </propertyList>
  <propertyList title="Force Summary">
  <!--Output: matrix fa(1:3,1:na_u)-->
   <property dictRef="siesta:forces">
   <!--In matrix, row (first) index is fastest-->
    <matrix units="siestaUnits:evpa" columns="8" rows="3" dataType="xsd:double">
 -6.197355209508E-01 -8.429398962120E-01  6.435049422438E-02
 -3.864484040841E-01  1.453251741970E+00 -2.024776678133E+00
 -1.005664472089E+00 -9.308548778506E-01  2.628911965965E-01
 -5.511234922795E-02  2.917646843209E-01 -1.919077862671E-01
  1.060274763758E+00 -8.014668316046E-01  2.418642537562E+00
  4.588450017855E-01 -7.019179044139E-01 -1.311349479598E+00
  3.815107595949E-01  1.326788463031E+00 -1.680143136781E-01
  1.590633543660E-01  2.132694335046E-01  9.530990777172E-01

    </matrix>
   </property>
   <property dictRef="siesta:ftot">
    <array units="siestaUnits:evpa" size="3" dataType="xsd:double">
 -7.266866848267E-03  7.894812745607E-03  2.935048424611E-03

    </array>
   </property>
  </propertyList>
  <property dictRef="siesta:stress">
  <!--In matrix, row (first) index is fastest-->
   <matrix units="siestaUnits:eV_Ang__3" columns="3" rows="3"
    dataType="xsd:double">
 -1.477790275028E-03  5.758656914044E-04 -5.132761407661E-04
  5.758815348726E-04 -4.941938065137E-04  3.226604157090E-04
 -5.132236213249E-04  3.226617436709E-04 -1.686352761209E-03

   </matrix>
  </property>
  <propertyList title="Final Pressure">
   <property title="cell volume" dictRef="siesta:cellvol">
    <scalar dataType="xsd:double" units="siestaUnits:Ang__3">1557.44225190                                                                                       
    </scalar>
   </property>
   <property title="Pressure of Solid" dictRef="siesta:pressSol">
    <scalar dataType="xsd:double" units="siestaUnits:kbar">1.95378838786                                                                                       
    </scalar>
   </property>
   <property title="Pressure of Molecule" dictRef="siesta:pressMol">
    <scalar dataType="xsd:double" units="siestaUnits:kbar">-0.251779325201E-01                                                                                 
    </scalar>
   </property>
  </propertyList>
  <property title="Electric dipole" dictRef="siesta:dipol">
   <array units="siestaUnits:Debye" size="3" dataType="xsd:double">
 -9.060557207789E-01 -1.246428318271E-01 -2.656023517781E+00

   </array>
  </property>
 </module>
 <metadata name="siesta:EndTime" content="2023-01-19T14-03-58" />
 <metadata name="dc:contributor" content="Siesta-CML" />
</cml>
