<?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-47-02" />
  <metadata name="siesta:run_UUID"
   content="9215bbf0-9810-11ed-5545-eff3bff4984c" />
  <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="O" id="a1" ref="siesta:e001"
     x3="    0.00000000                               "
     y3="    0.00000000                               "
     z3="    0.00000000                               " />
    <atom elementType="H" id="a2" ref="siesta:e002"
     x3="    0.75700000                               "
     y3="    0.58600000                               "
     z3="    0.00000000                               " />
    <atom elementType="H" id="a3" ref="siesta:e002"
     x3="   -0.75700000                               "
     y3="    0.58600000                               "
     z3="    0.00000000                               " />
   </atomArray>
  </molecule>
  <lattice dictRef="siesta:ucell">
   <latticeVector units="siestaUnits:angstrom" dictRef="cml:latticeVector">
  1.376932944731E+01  0.000000000000E+00  0.000000000000E+00

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

   </latticeVector>
   <latticeVector units="siestaUnits:angstrom" dictRef="cml:latticeVector">
  0.000000000000E+00  0.000000000000E+00  1.062217830507E+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">Water molecule Broyden Optimization
   </scalar>
  </parameter>
  <parameter name="SystemLabel" dictRef="siesta:slabel">
   <scalar dataType="xsd:string">h2o_op_broyden
   </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">2
   </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">50.0000000000                                                                                       
   </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">3
   </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.250000000000                                                                                      
   </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">Broyden
   </scalar>
  </parameter>
  <parameter name="MD.NumCGSteps">
   <scalar dataType="xsd:integer" units="cmlUnits:countable">40
   </scalar>
  </parameter>
  <parameter name="MD.Steps">
   <scalar dataType="xsd:integer" units="cmlUnits:countable">40
   </scalar>
  </parameter>
  <parameter name="MD.MaxCGDispl">
   <scalar dataType="xsd:double" units="siestaUnits:Bohr">0.200000000000                                                                                      
   </scalar>
  </parameter>
  <parameter name="MD.MaxDispl">
   <scalar dataType="xsd:double" units="siestaUnits:Bohr">0.200000000000                                                                                      
   </scalar>
  </parameter>
  <parameter name="MD.MaxForceTol">
   <scalar dataType="xsd:double" units="siestaUnits:Ry_Bohr">0.388934876913E-03                                                                                  
   </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">2.81050622447                                                                                       
   </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="Geom. Optim" role="step" serial="1">
  <molecule>
   <atomArray>
    <atom elementType="O" id="a1" ref="siesta:e001"
     x3="    0.00000000                               "
     y3="    0.00000000                               "
     z3="    0.00000000                               " />
    <atom elementType="H" id="a2" ref="siesta:e002"
     x3="    0.75700000                               "
     y3="    0.58600000                               "
     z3="    0.00000000                               " />
    <atom elementType="H" id="a3" ref="siesta:e002"
     x3="   -0.75700000                               "
     y3="    0.58600000                               "
     z3="    0.00000000                               " />
   </atomArray>
  </molecule>
  <lattice dictRef="siesta:ucell">
   <latticeVector units="siestaUnits:Ang" dictRef="cml:latticeVector">
  7.286412448942E+00  0.000000000000E+00  0.000000000000E+00

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

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

   </latticeVector>
  </lattice>
  <crystal title="Lattice Parameters">
   <scalar title="a" dictRef="cml:a" units="units:angstrom                ">
   7.286412                                                                                            
   </scalar>
   <scalar title="b" dictRef="cml:b" units="units:angstrom                ">
   5.746952                                                                                            
   </scalar>
   <scalar title="c" dictRef="cml:c" units="units:angstrom                ">
   5.621012                                                                                            
   </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">23
    </scalar>
   </property>
   <property title="Number of non-zeros" dictRef="siesta:nnz">
    <scalar dataType="xsd:integer" units="cmlUnits:countable">529
    </scalar>
   </property>
  </propertyList>
  <propertyList>
   <property title="Mesh" dictRef="siesta:ntm">
    <array units="cmlUnits:countable" size="3" dataType="xsd:integer">
          32          30          24

    </array>
   </property>
   <property title="Requested Cut-Off" dictRef="siesta:g2max">
    <scalar dataType="xsd:double" units="siestaUnits:Ry">50.0000000000                                                                                       
    </scalar>
   </property>
   <property title="Actual Cut-Off" dictRef="siesta:g2mesh">
    <scalar dataType="xsd:double" units="siestaUnits:Ry">50.3842796991                                                                                       
    </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">-86.527792                                                                                          
     </scalar>
    </property>
    <property dictRef="siesta:Eions">
     <scalar dataType="xsd:double" units="siestaUnits:eV">815.854478                                                                                          
     </scalar>
    </property>
    <property dictRef="siesta:Ena">
     <scalar dataType="xsd:double" units="siestaUnits:eV">175.155695                                                                                          
     </scalar>
    </property>
    <property dictRef="siesta:Ekin">
     <scalar dataType="xsd:double" units="siestaUnits:eV">370.751214                                                                                          
     </scalar>
    </property>
    <property dictRef="siesta:Enl">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-67.181693                                                                                          
     </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.606153                                                                                          
     </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">1.705518                                                                                            
     </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">-116.213035                                                                                         
     </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">-467.643428                                                                                         
     </scalar>
    </property>
    <property dictRef="siesta:EtotK">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.242933                                                                                         
     </scalar>
    </property>
    <property dictRef="siesta:FreeEK">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.242933                                                                                         
     </scalar>
    </property>
   </propertyList>
   <propertyList title="SCF Cycle">
    <property dictRef="siesta:Eharrs">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-467.6434280                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">1.4388662                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">6.3040746                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.2429330                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.2429330                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-4.2552919                                                                                          
     </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">-466.1233531                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0466230                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">2.5360579                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.7623481                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.7623481                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-3.2051785                                                                                          
     </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">-465.8481731                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0250624                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.3466942                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8349358                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8349358                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.4515035                                                                                          
     </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">-465.8397682                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0320852                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.1541250                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8387434                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8387434                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.2209850                                                                                          
     </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">-465.8392185                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0015386                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.0358775                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8390972                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8390972                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.1700609                                                                                          
     </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">-465.8391631                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0002866                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.0208247                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8391318                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8391318                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.1668704                                                                                          
     </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">-465.8391617                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0003431                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.0023643                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8391500                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8391500                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.1616059                                                                                          
     </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">-465.8391504                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0003269                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.0004521                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8391503                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8391503                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.1598101                                                                                          
     </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">-465.8391503                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0000103                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.0001846                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8391503                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8391503                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.1596914                                                                                          
     </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">-465.839150                                                                                         
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.839150                                                                                         
     </scalar>
    </property>
    <property dictRef="siesta:Ebs">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-104.739997                                                                                         
     </scalar>
    </property>
    <property dictRef="siesta:E_Fermi">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.159691                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:E_KS_egg">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.839150                                                                                         
     </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="3" rows="3" dataType="xsd:double">
 -1.122340851816E-09 -4.593892076947E-01  1.805256629294E-10
  6.905180607010E-01  2.576834296892E-01 -2.750331681990E-11
 -6.905180595793E-01  2.576834304964E-01 -2.747120345050E-11

     </matrix>
    </property>
    <property dictRef="siesta:ftot">
     <array units="siestaUnits:evpa" size="3" dataType="xsd:double">
 -6.622490401887E-13  5.597765249087E-02  1.255511426590E-10

     </array>
    </property>
    <property dictRef="siesta:fmax">
     <scalar dataType="xsd:double" units="siestaUnits:evpa">0.690518060701                                                                                      
     </scalar>
    </property>
    <property dictRef="siesta:fres">
     <scalar dataType="xsd:double" units="siestaUnits:evpa">0.379688516666                                                                                      
     </scalar>
    </property>
    <property dictRef="siesta:cfmax">
     <scalar dataType="xsd:double" units="siestaUnits:evpa">0.690518060701                                                                                      
     </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">
 -9.137915787974E-04 -2.052165532581E-13 -2.926627067093E-18
  1.907970616985E-13 -1.609189899336E-04 -7.316567667731E-18
  2.194970300319E-18  9.757008813303E-15  1.028906143710E-03

    </matrix>
   </property>
  </module>
 </module>
 <module dictRef="Geom. Optim" role="step" serial="2">
  <molecule>
   <atomArray>
    <atom elementType="O" id="a1" ref="siesta:e001"
     x3="   -0.00000000                               "
     y3="   -0.00945494                               "
     z3="    0.00000000                               " />
    <atom elementType="H" id="a2" ref="siesta:e002"
     x3="    0.77121192                               "
     y3="    0.59130352                               "
     z3="   -0.00000000                               " />
    <atom elementType="H" id="a3" ref="siesta:e002"
     x3="   -0.77121192                               "
     y3="    0.59130352                               "
     z3="   -0.00000000                               " />
   </atomArray>
  </molecule>
  <lattice dictRef="siesta:ucell">
   <latticeVector units="siestaUnits:Ang" dictRef="cml:latticeVector">
  7.286412448942E+00  0.000000000000E+00  0.000000000000E+00

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

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

   </latticeVector>
  </lattice>
  <crystal title="Lattice Parameters">
   <scalar title="a" dictRef="cml:a" units="units:angstrom                ">
   7.286412                                                                                            
   </scalar>
   <scalar title="b" dictRef="cml:b" units="units:angstrom                ">
   5.746952                                                                                            
   </scalar>
   <scalar title="c" dictRef="cml:c" units="units:angstrom                ">
   5.621012                                                                                            
   </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">23
    </scalar>
   </property>
   <property title="Number of non-zeros" dictRef="siesta:nnz">
    <scalar dataType="xsd:integer" units="cmlUnits:countable">529
    </scalar>
   </property>
  </propertyList>
  <propertyList>
   <property title="Mesh" dictRef="siesta:ntm">
    <array units="cmlUnits:countable" size="3" dataType="xsd:integer">
          32          30          24

    </array>
   </property>
   <property title="Requested Cut-Off" dictRef="siesta:g2max">
    <scalar dataType="xsd:double" units="siestaUnits:Ry">50.0000000000                                                                                       
    </scalar>
   </property>
   <property title="Actual Cut-Off" dictRef="siesta:g2mesh">
    <scalar dataType="xsd:double" units="siestaUnits:Ry">50.3842796991                                                                                       
    </scalar>
   </property>
  </propertyList>
  <module dictRef="SCF" role="step" serial="1">
   <propertyList title="SCF Cycle">
    <property dictRef="siesta:Eharrs">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-466.0562337                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0388280                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.3083349                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8409046                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8409046                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.1031654                                                                                          
     </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">-465.8251379                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0133396                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.5068838                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8372200                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8372200                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.2980387                                                                                          
     </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">-465.8445871                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0080998                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.0134312                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8426062                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8426062                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.1767395                                                                                          
     </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">-465.8426217                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0012358                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.0014919                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8426148                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8426148                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.1827753                                                                                          
     </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">-465.8426149                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0000323                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.0010513                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8426149                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8426149                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.1832856                                                                                          
     </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">-465.8426149                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0000203                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.0005076                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8426149                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8426149                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.1834053                                                                                          
     </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">-465.842615                                                                                         
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.842615                                                                                         
     </scalar>
    </property>
    <property dictRef="siesta:Ebs">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-103.641612                                                                                         
     </scalar>
    </property>
    <property dictRef="siesta:E_Fermi">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.183405                                                                                           
     </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="3" rows="3" dataType="xsd:double">
  4.611004661340E-08  1.169434038819E+00  1.378846237860E-09
 -1.525860234489E-01 -3.197141998867E-01  1.067108760345E-10
  1.525860194161E-01 -3.197141935578E-01  1.170049231235E-10

     </matrix>
    </property>
    <property dictRef="siesta:ftot">
     <array units="siestaUnits:evpa" size="3" dataType="xsd:double">
  4.207717476339E-08  5.300056453748E-01  1.602562037018E-09

     </array>
    </property>
    <property dictRef="siesta:fmax">
     <scalar dataType="xsd:double" units="siestaUnits:evpa">1.16943403882                                                                                       
     </scalar>
    </property>
    <property dictRef="siesta:fres">
     <scalar dataType="xsd:double" units="siestaUnits:evpa">0.424077468577                                                                                      
     </scalar>
    </property>
    <property dictRef="siesta:cfmax">
     <scalar dataType="xsd:double" units="siestaUnits:evpa">1.16943403882                                                                                       
     </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">
 -5.318923464583E-04  4.056178875618E-13 -3.786718090326E-16
  1.288486249781E-12  5.892584745689E-05  2.776510371100E-14
 -1.372993355077E-15 -8.904801782906E-15  1.041333185195E-03

    </matrix>
   </property>
  </module>
 </module>
 <module dictRef="Geom. Optim" role="step" serial="3">
  <molecule>
   <atomArray>
    <atom elementType="O" id="a1" ref="siesta:e001"
     x3="    0.00000000                               "
     y3="    0.00176573                               "
     z3="    0.00000000                               " />
    <atom elementType="H" id="a2" ref="siesta:e002"
     x3="    0.76974787                               "
     y3="    0.58823588                               "
     z3="    0.00000000                               " />
    <atom elementType="H" id="a3" ref="siesta:e002"
     x3="   -0.76974787                               "
     y3="    0.58823588                               "
     z3="    0.00000000                               " />
   </atomArray>
  </molecule>
  <lattice dictRef="siesta:ucell">
   <latticeVector units="siestaUnits:Ang" dictRef="cml:latticeVector">
  7.286412448942E+00  0.000000000000E+00  0.000000000000E+00

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

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

   </latticeVector>
  </lattice>
  <crystal title="Lattice Parameters">
   <scalar title="a" dictRef="cml:a" units="units:angstrom                ">
   7.286412                                                                                            
   </scalar>
   <scalar title="b" dictRef="cml:b" units="units:angstrom                ">
   5.746952                                                                                            
   </scalar>
   <scalar title="c" dictRef="cml:c" units="units:angstrom                ">
   5.621012                                                                                            
   </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">23
    </scalar>
   </property>
   <property title="Number of non-zeros" dictRef="siesta:nnz">
    <scalar dataType="xsd:integer" units="cmlUnits:countable">529
    </scalar>
   </property>
  </propertyList>
  <propertyList>
   <property title="Mesh" dictRef="siesta:ntm">
    <array units="cmlUnits:countable" size="3" dataType="xsd:integer">
          32          30          24

    </array>
   </property>
   <property title="Requested Cut-Off" dictRef="siesta:g2max">
    <scalar dataType="xsd:double" units="siestaUnits:Ry">50.0000000000                                                                                       
    </scalar>
   </property>
   <property title="Actual Cut-Off" dictRef="siesta:g2mesh">
    <scalar dataType="xsd:double" units="siestaUnits:Ry">50.3842796991                                                                                       
    </scalar>
   </property>
  </propertyList>
  <module dictRef="SCF" role="step" serial="1">
   <propertyList title="SCF Cycle">
    <property dictRef="siesta:Eharrs">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.7551909                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0212232                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.2281342                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8491678                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8491678                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.1727129                                                                                          
     </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">-465.8414152                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0092241                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.3616816                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8473696                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8473696                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.0251177                                                                                          
     </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">-465.8510092                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0055498                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.0082378                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8500258                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8500258                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.1157447                                                                                          
     </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">-465.8500317                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0004872                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.0010693                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8500291                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8500291                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.1109304                                                                                          
     </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">-465.8500292                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0000119                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.0006774                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8500292                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8500292                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.1109440                                                                                          
     </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">-465.850029                                                                                         
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.850029                                                                                         
     </scalar>
    </property>
    <property dictRef="siesta:Ebs">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-104.099601                                                                                         
     </scalar>
    </property>
    <property dictRef="siesta:E_Fermi">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.110944                                                                                           
     </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="3" rows="3" dataType="xsd:double">
 -1.000997800705E-07 -3.546208596379E-02  6.058151172853E-09
  1.595438287010E-01 -8.691127493432E-03  4.550377770837E-10
 -1.595437776967E-01 -8.691201354603E-03  2.271821379191E-10

     </matrix>
    </property>
    <property dictRef="siesta:ftot">
     <array units="siestaUnits:evpa" size="3" dataType="xsd:double">
 -4.909551273185E-08 -5.284441481183E-02  6.740371087856E-09

     </array>
    </property>
    <property dictRef="siesta:fmax">
     <scalar dataType="xsd:double" units="siestaUnits:evpa">0.159543828701                                                                                      
     </scalar>
    </property>
    <property dictRef="siesta:fres">
     <scalar dataType="xsd:double" units="siestaUnits:evpa">0.762430925813E-01                                                                                  
     </scalar>
    </property>
    <property dictRef="siesta:cfmax">
     <scalar dataType="xsd:double" units="siestaUnits:evpa">0.159543828701                                                                                      
     </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">
 -6.817623758492E-04 -1.316287493384E-11 -9.916180446153E-14
 -2.037112822385E-11 -6.507663910668E-05  1.652177072170E-13
 -1.901019848024E-13  4.197049068246E-14  1.033783467593E-03

    </matrix>
   </property>
  </module>
 </module>
 <module dictRef="Geom. Optim" role="step" serial="4">
  <molecule>
   <atomArray>
    <atom elementType="O" id="a1" ref="siesta:e001"
     x3="   -0.00000000                               "
     y3="    0.00301512                               "
     z3="    0.00000000                               " />
    <atom elementType="H" id="a2" ref="siesta:e002"
     x3="    0.77408381                               "
     y3="    0.58736035                               "
     z3="    0.00000000                               " />
    <atom elementType="H" id="a3" ref="siesta:e002"
     x3="   -0.77408381                               "
     y3="    0.58736035                               "
     z3="    0.00000000                               " />
   </atomArray>
  </molecule>
  <lattice dictRef="siesta:ucell">
   <latticeVector units="siestaUnits:Ang" dictRef="cml:latticeVector">
  7.286412448942E+00  0.000000000000E+00  0.000000000000E+00

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

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

   </latticeVector>
  </lattice>
  <crystal title="Lattice Parameters">
   <scalar title="a" dictRef="cml:a" units="units:angstrom                ">
   7.286412                                                                                            
   </scalar>
   <scalar title="b" dictRef="cml:b" units="units:angstrom                ">
   5.746952                                                                                            
   </scalar>
   <scalar title="c" dictRef="cml:c" units="units:angstrom                ">
   5.621012                                                                                            
   </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">23
    </scalar>
   </property>
   <property title="Number of non-zeros" dictRef="siesta:nnz">
    <scalar dataType="xsd:integer" units="cmlUnits:countable">529
    </scalar>
   </property>
  </propertyList>
  <propertyList>
   <property title="Mesh" dictRef="siesta:ntm">
    <array units="cmlUnits:countable" size="3" dataType="xsd:integer">
          32          30          24

    </array>
   </property>
   <property title="Requested Cut-Off" dictRef="siesta:g2max">
    <scalar dataType="xsd:double" units="siestaUnits:Ry">50.0000000000                                                                                       
    </scalar>
   </property>
   <property title="Actual Cut-Off" dictRef="siesta:g2mesh">
    <scalar dataType="xsd:double" units="siestaUnits:Ry">50.3842796991                                                                                       
    </scalar>
   </property>
  </propertyList>
  <module dictRef="SCF" role="step" serial="1">
   <propertyList title="SCF Cycle">
    <property dictRef="siesta:Eharrs">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8774127                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0063425                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.0120690                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8507731                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8507731                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.0870786                                                                                          
     </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">-465.8508419                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0005024                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.0071478                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8508396                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8508396                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.0829964                                                                                          
     </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">-465.8508431                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0001340                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.0015937                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8508417                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8508417                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.0845749                                                                                          
     </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">-465.8508419                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0000797                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.0000738                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8508418                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8508418                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.0852175                                                                                          
     </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">-465.850842                                                                                         
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.850842                                                                                         
     </scalar>
    </property>
    <property dictRef="siesta:Ebs">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-103.946350                                                                                         
     </scalar>
    </property>
    <property dictRef="siesta:E_Fermi">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.085217                                                                                           
     </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="3" rows="3" dataType="xsd:double">
  1.788929974692E-06 -2.600150845445E-02  6.999647297723E-08
  2.898507235605E-02 -4.725559618418E-02  2.982494629469E-08
 -2.898555295023E-02 -4.725506540928E-02  2.797569390722E-08

     </matrix>
    </property>
    <property dictRef="siesta:ftot">
     <array units="siestaUnits:evpa" size="3" dataType="xsd:double">
  1.308335800911E-06 -1.205121700479E-01  1.277971131792E-07

     </array>
    </property>
    <property dictRef="siesta:fmax">
     <scalar dataType="xsd:double" units="siestaUnits:evpa">0.472555961842E-01                                                                                  
     </scalar>
    </property>
    <property dictRef="siesta:fres">
     <scalar dataType="xsd:double" units="siestaUnits:evpa">0.275328242931E-01                                                                                  
     </scalar>
    </property>
    <property dictRef="siesta:cfmax">
     <scalar dataType="xsd:double" units="siestaUnits:evpa">0.472555961842E-01                                                                                  
     </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">
 -6.247779410491E-04  6.994625885453E-11 -7.053039238014E-14
  1.066580443775E-10 -5.155818991373E-05  2.641941563201E-12
 -2.061567026573E-14 -6.743931428529E-12  1.034577479876E-03

    </matrix>
   </property>
  </module>
 </module>
 <module dictRef="Geom. Optim" role="step" serial="5">
  <molecule>
   <atomArray>
    <atom elementType="O" id="a1" ref="siesta:e001"
     x3="    0.00000018                               "
     y3="   -0.00124595                               "
     z3="    0.00000001                               " />
    <atom elementType="H" id="a2" ref="siesta:e002"
     x3="    0.77634692                               "
     y3="    0.58299268                               "
     z3="    0.00000000                               " />
    <atom elementType="H" id="a3" ref="siesta:e002"
     x3="   -0.77634697                               "
     y3="    0.58299273                               "
     z3="    0.00000000                               " />
   </atomArray>
  </molecule>
  <lattice dictRef="siesta:ucell">
   <latticeVector units="siestaUnits:Ang" dictRef="cml:latticeVector">
  7.286412448942E+00  0.000000000000E+00  0.000000000000E+00

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

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

   </latticeVector>
  </lattice>
  <crystal title="Lattice Parameters">
   <scalar title="a" dictRef="cml:a" units="units:angstrom                ">
   7.286412                                                                                            
   </scalar>
   <scalar title="b" dictRef="cml:b" units="units:angstrom                ">
   5.746952                                                                                            
   </scalar>
   <scalar title="c" dictRef="cml:c" units="units:angstrom                ">
   5.621012                                                                                            
   </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">23
    </scalar>
   </property>
   <property title="Number of non-zeros" dictRef="siesta:nnz">
    <scalar dataType="xsd:integer" units="cmlUnits:countable">529
    </scalar>
   </property>
  </propertyList>
  <propertyList>
   <property title="Mesh" dictRef="siesta:ntm">
    <array units="cmlUnits:countable" size="3" dataType="xsd:integer">
          32          30          24

    </array>
   </property>
   <property title="Requested Cut-Off" dictRef="siesta:g2max">
    <scalar dataType="xsd:double" units="siestaUnits:Ry">50.0000000000                                                                                       
    </scalar>
   </property>
   <property title="Actual Cut-Off" dictRef="siesta:g2mesh">
    <scalar dataType="xsd:double" units="siestaUnits:Ry">50.3842796991                                                                                       
    </scalar>
   </property>
  </propertyList>
  <module dictRef="SCF" role="step" serial="1">
   <propertyList title="SCF Cycle">
    <property dictRef="siesta:Eharrs">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8704057                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0030617                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.0172204                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8508966                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8508966                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.0732396                                                                                          
     </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">-465.8508659                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0006738                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.0230235                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8508900                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8508900                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.0828711                                                                                          
     </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">-465.8509049                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0003892                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.0014932                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8509009                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8509009                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.0770952                                                                                          
     </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">-465.8509011                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0001043                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.0001817                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8509010                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8509010                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.0776295                                                                                          
     </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">-465.8509010                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0000016                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.0001169                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8509010                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8509010                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.0776544                                                                                          
     </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">-465.850901                                                                                         
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.850901                                                                                         
     </scalar>
    </property>
    <property dictRef="siesta:Ebs">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-103.844396                                                                                         
     </scalar>
    </property>
    <property dictRef="siesta:E_Fermi">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.077654                                                                                           
     </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="3" rows="3" dataType="xsd:double">
 -3.506535806401E-05  2.813657348058E-01 -1.487176129193E-06
 -5.579955357708E-02 -8.699223530377E-02 -4.811837787681E-08
  5.581428699605E-02 -8.700330160213E-02 -1.389260760155E-07

     </matrix>
    </property>
    <property dictRef="siesta:ftot">
     <array units="siestaUnits:evpa" size="3" dataType="xsd:double">
 -2.033193909027E-05  1.073701978999E-01 -1.674220583085E-06

     </array>
    </property>
    <property dictRef="siesta:fmax">
     <scalar dataType="xsd:double" units="siestaUnits:evpa">0.281365734806                                                                                      
     </scalar>
    </property>
    <property dictRef="siesta:fres">
     <scalar dataType="xsd:double" units="siestaUnits:evpa">0.105689653788                                                                                      
     </scalar>
    </property>
    <property dictRef="siesta:cfmax">
     <scalar dataType="xsd:double" units="siestaUnits:evpa">0.281365734806                                                                                      
     </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">
 -5.871887675727E-04 -2.551377015511E-09 -6.604923611534E-12
 -2.641239544455E-09 -3.512529758803E-05 -5.026474060105E-11
 -4.619723341550E-11 -2.114059315959E-11  1.035262283869E-03

    </matrix>
   </property>
  </module>
 </module>
 <module dictRef="Geom. Optim" role="step" serial="6">
  <molecule>
   <atomArray>
    <atom elementType="O" id="a1" ref="siesta:e001"
     x3="   -0.00000073                               "
     y3="    0.00082702                               "
     z3="   -0.00000003                               " />
    <atom elementType="H" id="a2" ref="siesta:e002"
     x3="    0.77741348                               "
     y3="    0.57985666                               "
     z3="    0.00000000                               " />
    <atom elementType="H" id="a3" ref="siesta:e002"
     x3="   -0.77741313                               "
     y3="    0.57985642                               "
     z3="    0.00000000                               " />
   </atomArray>
  </molecule>
  <lattice dictRef="siesta:ucell">
   <latticeVector units="siestaUnits:Ang" dictRef="cml:latticeVector">
  7.286412448942E+00  0.000000000000E+00  0.000000000000E+00

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

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

   </latticeVector>
  </lattice>
  <crystal title="Lattice Parameters">
   <scalar title="a" dictRef="cml:a" units="units:angstrom                ">
   7.286412                                                                                            
   </scalar>
   <scalar title="b" dictRef="cml:b" units="units:angstrom                ">
   5.746952                                                                                            
   </scalar>
   <scalar title="c" dictRef="cml:c" units="units:angstrom                ">
   5.621012                                                                                            
   </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">23
    </scalar>
   </property>
   <property title="Number of non-zeros" dictRef="siesta:nnz">
    <scalar dataType="xsd:integer" units="cmlUnits:countable">529
    </scalar>
   </property>
  </propertyList>
  <propertyList>
   <property title="Mesh" dictRef="siesta:ntm">
    <array units="cmlUnits:countable" size="3" dataType="xsd:integer">
          32          30          24

    </array>
   </property>
   <property title="Requested Cut-Off" dictRef="siesta:g2max">
    <scalar dataType="xsd:double" units="siestaUnits:Ry">50.0000000000                                                                                       
    </scalar>
   </property>
   <property title="Actual Cut-Off" dictRef="siesta:g2mesh">
    <scalar dataType="xsd:double" units="siestaUnits:Ry">50.3842796991                                                                                       
    </scalar>
   </property>
  </propertyList>
  <module dictRef="SCF" role="step" serial="1">
   <propertyList title="SCF Cycle">
    <property dictRef="siesta:Eharrs">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8312758                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0054868                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.0702755                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8512824                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8512824                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.0692427                                                                                          
     </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">-465.8505955                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0027949                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.1087957                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8511250                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8511250                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.0242002                                                                                          
     </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">-465.8514495                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0016880                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.0027811                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8513617                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8513617                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.0518020                                                                                          
     </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">-465.8513628                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0001258                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.0004506                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8513623                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8513623                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.0504691                                                                                          
     </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">-465.8513623                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:dDmax">
     <scalar dataType="xsd:double" units="siestaUnits:none">0.0000055                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:dHmax">
     <scalar dataType="xsd:double" units="siestaUnits:eV">0.0003155                                                                                           
     </scalar>
    </property>
    <property dictRef="siesta:Etot">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8513623                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.8513623                                                                                        
     </scalar>
    </property>
    <property dictRef="siesta:Ef">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.0504661                                                                                          
     </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">-465.851362                                                                                         
     </scalar>
    </property>
    <property dictRef="siesta:FreeE">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-465.851362                                                                                         
     </scalar>
    </property>
    <property dictRef="siesta:Ebs">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-103.934213                                                                                         
     </scalar>
    </property>
    <property dictRef="siesta:E_Fermi">
     <scalar dataType="xsd:double" units="siestaUnits:eV">-2.050466                                                                                           
     </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="3" rows="3" dataType="xsd:double">
  1.487791177769E-04  7.865031526580E-04 -3.019878952808E-06
 -1.671000447425E-03 -4.598999839943E-03  4.327918431040E-08
  1.605961545551E-03 -4.552521451716E-03 -7.544769889599E-08

     </matrix>
    </property>
    <property dictRef="siesta:ftot">
     <array units="siestaUnits:evpa" size="3" dataType="xsd:double">
  8.374021590257E-05 -8.365018139000E-03 -3.052047467393E-06

     </array>
    </property>
    <property dictRef="siesta:fmax">
     <scalar dataType="xsd:double" units="siestaUnits:evpa">0.459899983994E-02                                                                                  
     </scalar>
    </property>
    <property dictRef="siesta:fres">
     <scalar dataType="xsd:double" units="siestaUnits:evpa">0.230671333381E-02                                                                                  
     </scalar>
    </property>
    <property dictRef="siesta:cfmax">
     <scalar dataType="xsd:double" units="siestaUnits:evpa">0.459899983994E-02                                                                                  
     </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">
 -6.150830227168E-04  1.139190967652E-08  4.134469278149E-11
  1.124927414378E-08 -6.548804600269E-05  3.580214462886E-11
  9.656810517712E-12  3.979614514605E-11  1.033311487667E-03

    </matrix>
   </property>
  </module>
 </module>
 <module title="Finalization">
  <molecule>
   <atomArray>
    <atom elementType="O" id="a1" ref="siesta:e001"
     x3="   -0.00000073                               "
     y3="    0.00082702                               "
     z3="   -0.00000003                               " />
    <atom elementType="H" id="a2" ref="siesta:e002"
     x3="    0.77741348                               "
     y3="    0.57985666                               "
     z3="    0.00000000                               " />
    <atom elementType="H" id="a3" ref="siesta:e002"
     x3="   -0.77741313                               "
     y3="    0.57985642                               "
     z3="    0.00000000                               " />
   </atomArray>
  </molecule>
  <lattice dictRef="siesta:ucell">
   <latticeVector units="siestaUnits:Ang" dictRef="cml:latticeVector">
  7.286412448942E+00  0.000000000000E+00  0.000000000000E+00

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

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

   </latticeVector>
  </lattice>
  <propertyList title="Eigenvalues">
   <property title="Fermi Energy" dictRef="siesta:E_Fermi">
    <scalar dataType="xsd:double" units="siestaUnits:ev">-2.05047                                                                                            
    </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="23" dataType="xsd:double">
 -2.454929711498E+01 -1.267054704318E+01 -8.576226453909E+00 -6.171143289642E+00
  1.562742483423E+00  3.954709430157E+00  9.753905528969E+00  1.148979718554E+01
  1.228771406199E+01  1.315015865062E+01  1.453099464469E+01  2.147656361024E+01
  2.392595640133E+01  2.731146479341E+01  2.873961732658E+01  3.178845236432E+01
  4.258014878569E+01  4.262606612933E+01  4.575220420130E+01  4.842015932271E+01
  4.981308288646E+01  6.995590362209E+01  9.589742325661E+01

     </array>
    </property>
   </propertyList>
  </propertyList>
  <propertyList title="Energy Decomposition">
   <property dictRef="siesta:Ebs">
    <scalar dataType="xsd:double" units="siestaUnits:eV">-103.934213                                                                                         
    </scalar>
   </property>
   <property dictRef="siesta:Eions">
    <scalar dataType="xsd:double" units="siestaUnits:eV">815.854478                                                                                          
    </scalar>
   </property>
   <property dictRef="siesta:Ena">
    <scalar dataType="xsd:double" units="siestaUnits:eV">175.070603                                                                                          
    </scalar>
   </property>
   <property dictRef="siesta:Ekin">
    <scalar dataType="xsd:double" units="siestaUnits:eV">349.719698                                                                                          
    </scalar>
   </property>
   <property dictRef="siesta:Enl">
    <scalar dataType="xsd:double" units="siestaUnits:eV">-61.529236                                                                                          
    </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">-1.354824                                                                                           
    </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">0.730222                                                                                            
    </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">-112.633347                                                                                         
    </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">-465.851362                                                                                         
    </scalar>
   </property>
   <property dictRef="siesta:EtotK">
    <scalar dataType="xsd:double" units="siestaUnits:eV">-465.851362                                                                                         
    </scalar>
   </property>
   <property dictRef="siesta:FreeEK">
    <scalar dataType="xsd:double" units="siestaUnits:eV">-465.851362                                                                                         
    </scalar>
   </property>
  </propertyList>
  <propertyList title="Final Energy">
   <property dictRef="siesta:Ebs">
    <scalar dataType="xsd:double" units="siestaUnits:eV">-103.9342129                                                                                        
    </scalar>
   </property>
   <property dictRef="siesta:Ekin">
    <scalar dataType="xsd:double" units="siestaUnits:eV">349.7196983                                                                                         
    </scalar>
   </property>
   <property dictRef="siesta:Uscf">
    <scalar dataType="xsd:double" units="siestaUnits:eV">380.8667154                                                                                         
    </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">-112.6333465                                                                                        
    </scalar>
   </property>
   <property dictRef="siesta:I-e">
    <scalar dataType="xsd:double" units="siestaUnits:eV">-1067.8143569                                                                                       
    </scalar>
   </property>
   <property dictRef="siesta:I-I">
    <scalar dataType="xsd:double" units="siestaUnits:eV">-15.9900726                                                                                         
    </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">-465.8513623                                                                                        
    </scalar>
   </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">
 -8.368696590480E-03  1.549960446767E-07  5.625276210464E-10
  1.530553741454E-07 -8.910172563034E-04  4.871168193913E-10
  1.313886325419E-10  5.414583916281E-10  1.405902943890E-02

   </matrix>
  </property>
  <propertyList title="Final Pressure">
   <property title="cell volume" dictRef="siesta:cellvol">
    <scalar dataType="xsd:double" units="siestaUnits:Ang__3">235.378011719                                                                                       
    </scalar>
   </property>
   <property title="Pressure of Solid" dictRef="siesta:pressSol">
    <scalar dataType="xsd:double" units="siestaUnits:kbar">-2.56314480503                                                                                      
    </scalar>
   </property>
   <property title="Pressure of Molecule" dictRef="siesta:pressMol">
    <scalar dataType="xsd:double" units="siestaUnits:kbar">-2.55266787072                                                                                      
    </scalar>
   </property>
  </propertyList>
  <property title="Electric dipole" dictRef="siesta:dipol">
   <array units="siestaUnits:Debye" size="3" dataType="xsd:double">
  2.336331319487E-06  1.398833671568E+00  1.856562178083E-07

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