<?xml version="1.0" encoding="UTF-8"?>
<module xmlns="http://www.xml-cml.org/schema"
        xmlns:xi="http://www.w3.org/2001/XInclude"
        xmlns:nonsi="http://www.xml-cml.org/unit/nonSi/"
        xmlns:xsd="http://www.w3.org/2001/XMLSchema"
        xmlns:o="http://www.iochem-bd.org/dictionary/orca/"
        xmlns:convention="http://www.xml-cml.org/convention/"
        xmlns:cmlx="http://www.xml-cml.org/schema/cmlx"
        xmlns:cml="http://www.xml-cml.org/schema"
        xmlns:cc="http://www.xml-cml.org/dictionary/compchem/"
        xmlns:si="http://www.xml-cml.org/unit/si/"
        xmlns:nonsi2="http://www.iochem-bd.org/unit/nonSi2/"
        convention="convention:compchem"
        id="orca.log">
   <module dictRef="cc:jobList" id="jobList1">
      <module cmlx:templateRef="job" dictRef="cc:job" id="job">
         <module dictRef="cc:environment" id="environment">
            <parameterList>
               <parameter dictRef="cc:program">
                  <scalar dataType="xsd:string">Orca</scalar>
               </parameter>
               <parameter dictRef="cc:programVersion">
                  <scalar dataType="xsd:string" id="copy.0">3.0.1</scalar>
               </parameter>
               <parameter dictRef="cc:programSubversion">
                  <scalar dataType="xsd:string" id="copy.1">RELEASE</scalar>
               </parameter>
            </parameterList>
         </module>
         <module dictRef="cc:initialization" id="initialization">
            <parameterList>
               <module cmlx:templateRef="basis">
                  <list cmlx:templateRef="group">
                     <array dataType="xsd:integer" dictRef="o:group" size="3">1 2 3</array>
                     <array dataType="xsd:string" dictRef="o:primitive" size="3">17s11p7d1f 11s6p2d1f 5s1p</array>
                     <array dataType="xsd:string" dictRef="o:contraction" size="3">6s4p4d1f 5s3p2d1f 3s1p</array>
                  </list>
                  <list dictRef="atombasis">
                     <array dataType="xsd:integer" dictRef="cc:serial" size="4">0 1 2 3</array>
                     <array dataType="xsd:string" dictRef="cc:elementType" size="4">Fe O H H</array>
                     <array dataType="xsd:integer" dictRef="o:group" size="4">1 2 3 3</array>
                  </list>
               </module>
            </parameterList>
            <molecule id="molecule">
               <atomArray>
                  <atom elementType="Fe"
                        id="a1"
                        x3="-0.770308"
                        y3="-0.686274"
                        z3="0.0000"/>
                  <atom elementType="O"
                        id="a2"
                        x3="0.97924"
                        y3="-0.726037"
                        z3="0.0000"/>
                  <atom elementType="H"
                        id="a3"
                        x3="0.699794"
                        y3="0.543573"
                        z3="-0.0000"/>
                  <atom elementType="H"
                        id="a4"
                        x3="-0.310341"
                        y3="1.136589"
                        z3="0.0000"/>
               </atomArray>
               <bondArray>
                  <bond atomRefs2="a1 a2" order="S"/>
                  <bond atomRefs2="a2 a3" order="S"/>
               </bondArray>
               <formula concise="H 2 Fe 1 O 1">
                  <atomArray count="2 1 1" elementType="H Fe O"/>
               </formula>
               <property dictRef="cml:molmass">
                  <scalar units="unit:dalton">71.8444</scalar>
               </property>
            </molecule>
            <module dictRef="cc:userDefinedModule" id="otherComponents">
               <module cmlx:templateRef="scfsettings" dictRef="cc:userDefinedModule">
                  <module cmlx:templateRef="section" name="Hamiltonian">
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">Method</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">DFT(GTOs)</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">Exchange</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">PBE</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">XKappa</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">0.804000</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">Correlation</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">PBE</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">LDAOpt</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">PW91-LDA</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">PostSCFGGA</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">off</scalar>
                     </list>
                     <scalar dataType="xsd:string" dictRef="o:comment">Hybrid DFT is turned on</scalar>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">ScalHFX</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">0.250000</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">ScalDFX</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">0.750000</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">ScalDFC</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">1.000000</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">ScalLDAC</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">1.000000</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">Perturbative correction</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">0.000000</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">NL short-range parameter</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">6.900000</scalar>
                     </list>
                  </module>
                  <module cmlx:templateRef="section" name="General Settings">
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">IntName</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">run</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">HFTyp</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">UHF</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">Charge</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">1</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">Mult</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">6</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">NEL</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">35</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">Dim</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">88</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">ENuc</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">82.9701002563 Eh</scalar>
                     </list>
                  </module>
                  <module cmlx:templateRef="section" name="Convergence Acceleration">
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">CNVDIIS</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">on</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">DIISMaxIt</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">12</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">DIISStart</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">0.200000</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">DIISMaxEq</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">5</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">DIISBfac</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">1.050</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">DIISMaxC</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">10.000</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">CNVNR</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">off</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">CNVSOSCF</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">off</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">CNVShift</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">on</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">LevelShift</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">0.2500</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">ShiftErr</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">0.0010</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">CNVZerner</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">off</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">CNVDamp</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">on</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">DampFac</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">0.7000</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">DampMax</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">0.9800</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">DampMin</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">0.0000</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">DampErr</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">0.1000</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">CNVRico</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">off</scalar>
                     </list>
                  </module>
                  <module cmlx:templateRef="section" name="SCF Procedure">
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">MaxIter</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">125</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">SCFMode</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">Direct</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">Integral package</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">LIBINT</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">DirectResetFreq</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">20</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">Thresh</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">2.500e-11 Eh</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">TCut</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">2.500e-12 Eh</scalar>
                     </list>
                  </module>
                  <module cmlx:templateRef="section" name="Convergence Tolerance">
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">ConvCheckMode</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">Total+1el-Energy</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">TolE</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">1.000e-08 Eh</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">1-El. energy change</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">1.000e-05 Eh</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">TolErr</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">5.000e-07</scalar>
                     </list>
                  </module>
                  <module cmlx:templateRef="section" name="Diagonalization of the overlap matrix">
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">Smallest eigenvalue</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">4.177e-03</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">Time for diagonalization</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">0.013 sec</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">Threshold for overlap eigenvalues</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">1.000e-08</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">Number of eigenvalues below threshold</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">0</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">Time for construction of square roots</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">0.001 sec</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">Total time needed</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">0.013 sec</scalar>
                     </list>
                  </module>
               </module>
               <module cmlx:templateRef="scfsettings" dictRef="cc:userDefinedModule">
                  <module cmlx:templateRef="section" name="Hamiltonian">
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">Method</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">DFT(GTOs)</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">Exchange</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">PBE</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">XKappa</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">0.804000</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">Correlation</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">PBE</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">LDAOpt</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">PW91-LDA</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">PostSCFGGA</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">off</scalar>
                     </list>
                     <scalar dataType="xsd:string" dictRef="o:comment">Hybrid DFT is turned on</scalar>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">ScalHFX</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">0.250000</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">ScalDFX</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">0.750000</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">ScalDFC</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">1.000000</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">ScalLDAC</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">1.000000</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">Perturbative correction</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">0.000000</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">NL short-range parameter</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">6.900000</scalar>
                     </list>
                  </module>
                  <module cmlx:templateRef="section" name="General Settings">
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">IntName</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">run</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">HFTyp</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">UHF</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">Charge</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">1</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">Mult</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">6</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">NEL</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">35</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">Dim</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">88</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">ENuc</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">85.5024980232 Eh</scalar>
                     </list>
                  </module>
                  <module cmlx:templateRef="section" name="Convergence Acceleration">
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">CNVDIIS</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">on</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">DIISMaxIt</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">12</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">DIISStart</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">0.200000</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">DIISMaxEq</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">5</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">DIISBfac</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">1.050</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">DIISMaxC</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">10.000</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">CNVNR</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">off</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">CNVSOSCF</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">off</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">CNVShift</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">on</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">LevelShift</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">0.2500</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">ShiftErr</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">0.0010</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">CNVZerner</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">off</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">CNVDamp</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">on</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">DampFac</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">0.7000</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">DampMax</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">0.9800</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">DampMin</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">0.0000</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">DampErr</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">0.1000</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">CNVRico</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">off</scalar>
                     </list>
                  </module>
                  <module cmlx:templateRef="section" name="SCF Procedure">
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">MaxIter</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">125</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">SCFMode</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">Direct</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">Integral package</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">LIBINT</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">DirectResetFreq</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">20</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">Thresh</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">2.500e-11 Eh</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">TCut</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">2.500e-12 Eh</scalar>
                     </list>
                  </module>
                  <module cmlx:templateRef="section" name="Convergence Tolerance">
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">ConvCheckMode</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">Total+1el-Energy</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">TolE</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">1.000e-08 Eh</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">1-El. energy change</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">1.000e-05 Eh</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">TolErr</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">5.000e-07</scalar>
                     </list>
                  </module>
                  <module cmlx:templateRef="section" name="Diagonalization of the overlap matrix">
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">Smallest eigenvalue</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">3.795e-03</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">Time for diagonalization</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">0.003 sec</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">Threshold for overlap eigenvalues</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">1.000e-08</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">Number of eigenvalues below threshold</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">0</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">Time for construction of square roots</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">0.001 sec</scalar>
                     </list>
                     <list>
                        <scalar dataType="xsd:string" dictRef="cc:parameter">Total time needed</scalar>
                        <scalar dataType="xsd:string" dictRef="cc:value">0.004 sec</scalar>
                     </list>
                  </module>
               </module>
               <module cmlx:templateRef="input" id="initialization">
                  <module cmlx:templateRef="job">
                     <scalar dataType="xsd:integer" dictRef="o:charge">1</scalar>
                     <scalar dataType="xsd:integer" dictRef="cc:multiplicity">6</scalar>
                     <array dataType="xsd:string" dictRef="cc:keywords" size="5">PBE0 def2-TZVP TightSCF OptTS Freq</array>
                     <module cmlx:templateRef="block">
                        <scalar dataType="xsd:string" dictRef="o:type">PAL</scalar>
                        <scalar dataType="xsd:string" dictRef="o:parameters">NPROCS 8 END</scalar>
                     </module>
                     <module cmlx:templateRef="block">
                        <scalar dataType="xsd:string" dictRef="o:type">geom</scalar>
                        <scalar dataType="xsd:string" dictRef="o:parameters">Calc_Hess true</scalar>
                     </module>
                  </module>
               </module>
            </module>
         </module>
         <module dictRef="cc:calculation" id="calculation">
            <module dictRef="cc:userDefined" id="otherComponents">
               <module cmlx:templateRef="geometry">
                  <molecule id="molecule">
                     <atomArray>
                        <atom elementType="Fe"
                              id="a1"
                              x3="-0.738802"
                              y3="-0.647625"
                              z3="0.0000"/>
                        <atom elementType="O"
                              id="a2"
                              x3="0.972411"
                              y3="-0.705192"
                              z3="-0.0000"/>
                        <atom elementType="H"
                              id="a3"
                              x3="0.614203"
                              y3="0.540548"
                              z3="0.0000"/>
                        <atom elementType="H"
                              id="a4"
                              x3="-0.249425"
                              y3="1.08012"
                              z3="-0.0000"/>
                     </atomArray>
                     <bondArray>
                        <bond atomRefs2="a1 a3" order="S"/>
                        <bond atomRefs2="a1 a4" order="S"/>
                        <bond atomRefs2="a1 a2" order="S"/>
                        <bond atomRefs2="a2 a3" order="S"/>
                     </bondArray>
                     <formula concise="H 2 Fe 1 O 1">
                        <atomArray count="2 1 1" elementType="H Fe O"/>
                     </formula>
                     <property dictRef="cml:molmass">
                        <scalar units="unit:dalton">71.8444</scalar>
                     </property>
                  </molecule>
               </module>
               <module cmlx:templateRef="geometry">
                  <molecule id="molecule">
                     <atomArray>
                        <atom elementType="Fe"
                              id="a1"
                              x3="-0.700512"
                              y3="-0.587008"
                              z3="0.0000"/>
                        <atom elementType="O"
                              id="a2"
                              x3="0.999006"
                              y3="-0.75414"
                              z3="-0.0000"/>
                        <atom elementType="H"
                              id="a3"
                              x3="0.586743"
                              y3="0.529723"
                              z3="0.0000"/>
                        <atom elementType="H"
                              id="a4"
                              x3="-0.286851"
                              y3="1.079276"
                              z3="-0.0000"/>
                     </atomArray>
                     <bondArray>
                        <bond atomRefs2="a1 a4" order="S"/>
                        <bond atomRefs2="a1 a2" order="S"/>
                        <bond atomRefs2="a1 a3" order="S"/>
                     </bondArray>
                     <formula concise="H 2 Fe 1 O 1">
                        <atomArray count="2 1 1" elementType="H Fe O"/>
                     </formula>
                     <property dictRef="cml:molmass">
                        <scalar units="unit:dalton">71.8444</scalar>
                     </property>
                  </molecule>
               </module>
               <module cmlx:templateRef="geometry">
                  <molecule id="molecule">
                     <atomArray>
                        <atom elementType="Fe"
                              id="a1"
                              x3="-0.727797"
                              y3="-0.64234"
                              z3="0.0000"/>
                        <atom elementType="O"
                              id="a2"
                              x3="0.97919"
                              y3="-0.720129"
                              z3="-0.0000"/>
                        <atom elementType="H"
                              id="a3"
                              x3="0.596042"
                              y3="0.54124"
                              z3="0.0000"/>
                        <atom elementType="H"
                              id="a4"
                              x3="-0.249049"
                              y3="1.089079"
                              z3="-0.0000"/>
                     </atomArray>
                     <bondArray>
                        <bond atomRefs2="a1 a4" order="S"/>
                        <bond atomRefs2="a1 a3" order="S"/>
                        <bond atomRefs2="a1 a2" order="S"/>
                        <bond atomRefs2="a2 a3" order="S"/>
                     </bondArray>
                     <formula concise="H 2 Fe 1 O 1">
                        <atomArray count="2 1 1" elementType="H Fe O"/>
                     </formula>
                     <property dictRef="cml:molmass">
                        <scalar units="unit:dalton">71.8444</scalar>
                     </property>
                  </molecule>
               </module>
               <module cmlx:templateRef="geometry">
                  <molecule id="molecule">
                     <atomArray>
                        <atom elementType="Fe"
                              id="a1"
                              x3="-0.727572"
                              y3="-0.642564"
                              z3="0.0000"/>
                        <atom elementType="O"
                              id="a2"
                              x3="0.978725"
                              y3="-0.722633"
                              z3="-0.0000"/>
                        <atom elementType="H"
                              id="a3"
                              x3="0.595741"
                              y3="0.545483"
                              z3="0.0000"/>
                        <atom elementType="H"
                              id="a4"
                              x3="-0.248508"
                              y3="1.087564"
                              z3="-0.0000"/>
                     </atomArray>
                     <bondArray>
                        <bond atomRefs2="a1 a4" order="S"/>
                        <bond atomRefs2="a1 a3" order="S"/>
                        <bond atomRefs2="a1 a2" order="S"/>
                        <bond atomRefs2="a2 a3" order="S"/>
                     </bondArray>
                     <formula concise="H 2 Fe 1 O 1">
                        <atomArray count="2 1 1" elementType="H Fe O"/>
                     </formula>
                     <property dictRef="cml:molmass">
                        <scalar units="unit:dalton">71.8444</scalar>
                     </property>
                  </molecule>
               </module>
               <module cmlx:templateRef="geometry">
                  <molecule id="molecule">
                     <atomArray>
                        <atom elementType="Fe"
                              id="a1"
                              x3="-0.726951"
                              y3="-0.643166"
                              z3="0.0000"/>
                        <atom elementType="O"
                              id="a2"
                              x3="0.978742"
                              y3="-0.723645"
                              z3="-0.0000"/>
                        <atom elementType="H"
                              id="a3"
                              x3="0.595046"
                              y3="0.54575"
                              z3="0.0000"/>
                        <atom elementType="H"
                              id="a4"
                              x3="-0.248451"
                              y3="1.088912"
                              z3="-0.0000"/>
                     </atomArray>
                     <bondArray>
                        <bond atomRefs2="a1 a4" order="S"/>
                        <bond atomRefs2="a1 a3" order="S"/>
                        <bond atomRefs2="a1 a2" order="S"/>
                        <bond atomRefs2="a2 a3" order="S"/>
                     </bondArray>
                     <formula concise="H 2 Fe 1 O 1">
                        <atomArray count="2 1 1" elementType="H Fe O"/>
                     </formula>
                     <property dictRef="cml:molmass">
                        <scalar units="unit:dalton">71.8444</scalar>
                     </property>
                  </molecule>
               </module>
               <module cmlx:templateRef="finalspenergy">
                  <scalar dataType="xsd:double" dictRef="cc:energy">-1339.449443295649</scalar>
               </module>
               <module cmlx:templateRef="finalspenergy">
                  <scalar dataType="xsd:double" dictRef="cc:energy">-1339.452410037394</scalar>
               </module>
               <module cmlx:templateRef="finalspenergy">
                  <scalar dataType="xsd:double" dictRef="cc:energy">-1339.451140188418</scalar>
               </module>
               <module cmlx:templateRef="finalspenergy">
                  <scalar dataType="xsd:double" dictRef="cc:energy">-1339.452315332159</scalar>
               </module>
               <module cmlx:templateRef="finalspenergy">
                  <scalar dataType="xsd:double" dictRef="cc:energy">-1339.452310709526</scalar>
               </module>
               <module cmlx:templateRef="finalspenergy">
                  <scalar dataType="xsd:double" dictRef="cc:energy">-1339.452310871656</scalar>
               </module>
               <module cmlx:templateRef="loewdin">
                  <array dataType="xsd:integer" dictRef="cc:serial" size="4">0 1 2 3</array>
                  <array dataType="xsd:string" dictRef="cc:elementType" size="4">Fe O H H</array>
                  <array dataType="xsd:double" dictRef="x:charge" size="4">0.883806 -0.193583 0.249304 0.060473</array>
                  <array dataType="xsd:double" dictRef="x:spin" size="4">3.884912 0.867358 0.012062 0.235668</array>
               </module>
               <module cmlx:templateRef="mayer">
                  <module cmlx:templateRef="mayercharges">
                     <array dataType="xsd:integer" dictRef="cc:serial" size="4">0 1 2 3</array>
                     <array dataType="xsd:string" dictRef="cc:elementType" size="4">Fe O H H</array>
                     <array dataType="xsd:double" dictRef="o:na" size="4">25.0175 8.3092 0.6976 0.9757</array>
                     <array dataType="xsd:double" dictRef="o:za" size="4">26.0000 8.0000 1.0000 1.0000</array>
                     <array dataType="xsd:double" dictRef="o:qa" size="4">0.9825 -0.3092 0.3024 0.0243</array>
                     <array dataType="xsd:double" dictRef="o:va" size="4">5.7101 2.2577 0.9518 0.9970</array>
                     <array dataType="xsd:double" dictRef="o:bva" size="4">2.3090 1.9200 0.9495 0.9113</array>
                     <array dataType="xsd:double" dictRef="o:fa" size="4">3.4011 0.3377 0.0024 0.0856</array>
                  </module>
                  <module cmlx:templateRef="bonds">
                     <array dataType="xsd:double" dictRef="x:distance" size="5">1.4787 0.2194 0.6109 0.4355 0.2946</array>
                     <matrix cols="2" dataType="xsd:integer" dictRef="x:serial" rows="5">0 1 0 2 0 3 1 2 2 3</matrix>
                  </module>
               </module>
               <module cmlx:templateRef="mullikenpopulation">
                  <module cmlx:templateRef="atomicchargesspin">
                     <scalar dataType="xsd:double" dictRef="x:chargesum">1.0000000</scalar>
                     <scalar dataType="xsd:double" dictRef="x:spinsum">5.0000000</scalar>
                     <array dataType="xsd:integer" dictRef="cc:serial" size="4">0 1 2 3</array>
                     <array dataType="xsd:string" dictRef="cc:elementType" size="4">Fe O H H</array>
                     <array dataType="xsd:double" dictRef="x:charge" size="4">0.982477 -0.309176 0.302351 0.024348</array>
                     <array dataType="xsd:double" dictRef="x:spin" size="4">3.856024 0.879919 -0.027284 0.291341</array>
                  </module>
               </module>
               <module cmlx:templateRef="irspectrum">
                  <array dataType="xsd:integer" dictRef="cc:serial" size="5">7 8 9 10 11</array>
                  <array dataType="xsd:double" dictRef="cc:frequency" size="5">467.29 698.39 1062.04 1271.40 1555.97</array>
                  <array dataType="xsd:double" dictRef="o:t2" size="5">35.156402 2.893576 60.673627 15.461483 101.369210</array>
                  <matrix cols="3"
                          dataType="xsd:double"
                          dictRef="cc:displacement"
                          rows="5">-4.604941 3.735093 0.000000 1.330438 1.059958 0.000000 -0.000048 -0.000008 -7.789328 -2.702912 2.855828 0.000000 8.642844 5.164344 0.000000</matrix>
               </module>
               <module cmlx:templateRef="orbitalenergies">
                  <list cmlx:templateRef="orbital">
                     <array dataType="xsd:integer" dictRef="cc:serial" size="88">0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87</array>
                     <array dataType="xsd:double" dictRef="cc:occup" size="88">1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000</array>
                     <array dataType="xsd:double"
                            dictRef="cc:energy"
                            size="88"
                            units="nonsi:electronvolt">-6985.9849 -830.8537 -717.3014 -717.0894 -716.8890 -530.5648 -107.6949 -74.4147 -74.0036 -73.6508 -33.7770 -21.3064 -20.1045 -19.9024 -19.8252 -19.4917 -18.3679 -16.2359 -16.2047 -14.3188 -9.4361 -5.5674 -5.4519 -2.6021 -2.2343 -1.9622 -1.3559 -0.0578 0.0521 0.0882 2.7176 3.3737 5.2183 8.3365 9.1664 10.5177 15.3246 15.3708 15.9690 15.9968 16.5277 19.6530 26.5882 30.4203 31.4270 31.7766 35.1381 35.2880 38.0378 38.7898 39.3955 42.3424 45.9173 50.5112 52.4491 57.0564 60.8557 64.1645 64.6367 66.4956 71.0732 72.3954 82.0055 83.9996 84.3733 85.2064 86.5949 91.8800 99.3560 99.7241 99.9457 100.6268 100.9577 102.0572 107.8862 132.3732 133.2247 133.6188 134.3032 134.5002 140.8842 142.5709 162.9722 165.6314 166.0487 173.5696 179.7981 1173.2199</array>
                     <scalar dataType="xsd:string" dictRef="o:spin">UP</scalar>
                  </list>
                  <list cmlx:templateRef="orbital">
                     <array dataType="xsd:integer" dictRef="cc:serial" size="88">0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87</array>
                     <array dataType="xsd:double" dictRef="cc:occup" size="88">1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000</array>
                     <array dataType="xsd:double"
                            dictRef="cc:energy"
                            size="88"
                            units="nonsi:electronvolt">-6985.9398 -828.5926 -715.2909 -715.1936 -715.0818 -529.8673 -102.2085 -68.2474 -68.1298 -67.9795 -31.8894 -19.1709 -16.6308 -15.3587 -15.1857 -10.5380 -10.3258 -10.0038 -9.5220 -9.0060 -8.1783 -4.8060 -4.2244 -1.7331 -1.5672 -1.0303 -0.5868 0.4271 0.6482 0.6712 3.1907 3.7263 6.1068 8.9508 9.9406 10.9665 16.0350 17.0594 17.6785 17.8084 17.9354 21.1972 27.1216 31.3940 32.8831 33.2718 35.9977 36.7539 39.1154 40.6218 41.0642 43.7301 46.5330 51.9127 53.5096 58.2813 61.8363 65.2013 66.2479 67.7132 72.8500 73.8252 83.1966 87.4150 87.9320 88.2932 90.2247 94.5969 103.8963 104.0936 104.2006 104.4772 105.4612 105.7094 111.0212 133.4195 134.2152 134.9690 135.5545 135.7665 142.1566 143.8306 164.8469 167.3174 167.5486 175.7233 181.2019 1173.8470</array>
                     <scalar dataType="xsd:string" dictRef="o:spin">DOWN</scalar>
                  </list>
               </module>
               <module cmlx:templateRef="vibrations" dictRef="cc:vibrations">
                  <array dataType="xsd:double" dictRef="cc:frequency" size="12">0.00 0.00 0.00 0.00 0.00 0.00 -1996.02 467.29 698.39 1062.04 1271.40 1555.97</array>
                  <matrix cols="12"
                          dataType="xsd:double"
                          dictRef="cc:displacement"
                          rows="12">0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 -0.004608 0.002432 0.000000 0.003510 0.039876 0.000000 0.574157 -0.725709 0.000000 -0.374564 -0.041977 0.000000 -0.032464 -0.025023 0.000000 0.031739 0.071098 -0.000000 0.453123 -0.020928 0.000000 0.841854 0.278893 -0.000000 0.176242 -0.021293 0.000000 -0.651535 0.044028 -0.000000 0.225351 -0.114831 0.000000 0.350827 0.595819 -0.000000 0.000000 -0.000000 -0.005171 -0.000000 0.000000 0.034175 0.000000 -0.000000 -0.822951 0.000000 0.000000 0.567060 -0.004676 -0.006410 0.000000 0.038295 -0.037332 0.000000 -0.081275 -0.010848 0.000000 -0.267493 0.958541 0.000000 -0.003705 -0.004269 0.000000 -0.009945 -0.025270 0.000000 0.672336 0.671167 0.000000 -0.309195 -0.033579 0.000000</matrix>
               </module>
            </module>
         </module>
         <module dictRef="cc:finalization" id="finalization">
            <propertyList>
               <property dictRef="cc:frequencies">
                  <module cmlx:templateRef="vibrations" dictRef="cc:vibrations">
                     <array dataType="xsd:double" dictRef="cc:frequency" size="12">0.00 0.00 0.00 0.00 0.00 0.00 -1676.22 683.16 777.26 1022.33 1610.31 1784.31</array>
                     <matrix cols="12"
                             dataType="xsd:double"
                             dictRef="cc:displacement"
                             rows="12">0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 -0.007357 0.002678 0.000000 0.006056 0.047897 0.000000 0.573697 -0.760052 -0.000000 -0.262162 -0.148566 -0.000000 -0.032856 -0.024328 0.000000 0.033108 0.064515 -0.000000 0.440973 0.070288 -0.000000 0.854011 0.253650 -0.000000 0.187621 -0.026480 -0.000000 -0.703424 0.070878 0.000000 0.368895 -0.062671 -0.000000 0.400419 0.404871 -0.000000 -0.000000 -0.000000 -0.004236 -0.000000 0.000000 0.029473 0.000000 0.000000 -0.813675 0.000000 0.000000 0.580557 0.000852 -0.004176 0.000000 0.020682 -0.014306 0.000000 -0.316732 -0.399475 -0.000000 -0.058729 0.857905 -0.000000 -0.004988 -0.006456 -0.000000 0.008277 -0.035031 0.000000 0.604769 0.502031 0.000000 -0.459774 0.411662 0.000000</matrix>
                  </module>
               </property>
               <property dictRef="cc:irspectrum">
                  <module cmlx:templateRef="irspectrum">
                     <array dataType="xsd:integer" dictRef="cc:serial" size="5">7 8 9 10 11</array>
                     <array dataType="xsd:double" dictRef="cc:frequency" size="5">683.16 777.26 1022.33 1610.31 1784.31</array>
                     <array dataType="xsd:double" dictRef="o:t2" size="5">17.266521 6.127758 37.798653 18.139182 56.649226</array>
                     <matrix cols="3"
                             dataType="xsd:double"
                             dictRef="cc:displacement"
                             rows="5">-3.455389 2.307988 0.000000 2.474405 -0.071271 -0.000000 0.000002 -0.000021 -6.148061 -2.917240 -3.103045 -0.000000 6.928520 2.940211 -0.000000</matrix>
                  </module>
               </property>
            </propertyList>
            <molecule id="molecule">
               <atomArray>
                  <atom elementType="Fe"
                        id="a1"
                        x3="-0.727137"
                        y3="-0.643194"
                        z3="0.0000"/>
                  <atom elementType="O"
                        id="a2"
                        x3="0.978659"
                        y3="-0.723608"
                        z3="-0.0000"/>
                  <atom elementType="H"
                        id="a3"
                        x3="0.595196"
                        y3="0.546143"
                        z3="0.0000"/>
                  <atom elementType="H"
                        id="a4"
                        x3="-0.248332"
                        y3="1.088509"
                        z3="-0.0000"/>
               </atomArray>
               <bondArray>
                  <bond atomRefs2="a1 a4" order="S"/>
                  <bond atomRefs2="a1 a3" order="S"/>
                  <bond atomRefs2="a1 a2" order="S"/>
                  <bond atomRefs2="a2 a3" order="S"/>
               </bondArray>
               <formula concise="H 2 Fe 1 O 1">
                  <atomArray count="2 1 1" elementType="H Fe O"/>
               </formula>
               <property dictRef="cml:molmass">
                  <scalar units="unit:dalton">71.8444</scalar>
               </property>
            </molecule>
            <module dictRef="cc:userDefinedModule" id="otherComponents">
               <module cmlx:templateRef="totalenergy" dictRef="cc:userDefinedModule">
                  <scalar dataType="xsd:double" dictRef="cc:totalener" units="nonsi:hartree">-1339.44944330</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:nucrepener" units="nonsi:hartree">82.97010026</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:electener" units="nonsi:hartree">-1422.41954355</scalar>
                  <scalar dataType="xsd:double"
                          dictRef="cc:oneelecener"
                          units="nonsi:hartree">-2005.00596519</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:twoeener" units="nonsi:hartree">582.58642164</scalar>
                  <scalar dataType="xsd:double"
                          dictRef="cc:potentialEnergy"
                          units="nonsi:hartree">-2676.90315720</scalar>
                  <scalar dataType="xsd:double"
                          dictRef="cc:kineticenergy"
                          units="nonsi:hartree">1337.45371391</scalar>
                  <scalar dataType="xsd:double" dictRef="o:vircoeff">2.00149219</scalar>
                  <list id="dftcomponents">
                     <scalar dataType="xsd:double" dictRef="cc:alphae">20.000008427385</scalar>
                     <scalar dataType="xsd:double" dictRef="cc:betae">15.000009760159</scalar>
                     <scalar dataType="xsd:double" dictRef="cc:totale">35.000018187544</scalar>
                     <scalar dataType="xsd:double"
                             dictRef="o:exchangeener"
                             units="nonsi:hartree">-47.162577661347</scalar>
                     <scalar dataType="xsd:double"
                             dictRef="o:correlationener"
                             units="nonsi:hartree">-1.435466257453</scalar>
                     <scalar dataType="xsd:double" dictRef="o:xcener" units="nonsi:hartree">-48.598043918799</scalar>
                  </list>
               </module>
               <module cmlx:templateRef="spincontamination" dictRef="cc:userDefinedModule">
                  <scalar dataType="xsd:double" dictRef="cc:s2">8.772355</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:sideal">2.5</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:s2ideal">8.750000</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:s2deviation">0.022355</scalar>
               </module>
               <module cmlx:templateRef="thermochemistry" dictRef="cc:userDefinedModule">
                  <scalar dataType="xsd:double" dictRef="cc:temp" units="si:k">298.15</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:press" units="nonsi:atm">1.00</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:totalmass">73.86</scalar>
               </module>
               <module cmlx:templateRef="enthalpy" dictRef="cc:userDefinedModule">
                  <scalar dataType="xsd:double" dictRef="o:totalfree" units="nonsi:hartree">-1339.43468623</scalar>
                  <scalar dataType="xsd:double"
                          dictRef="o:thermalcorrenthalpy"
                          units="nonsi:hartree">0.00094421</scalar>
                  <scalar dataType="xsd:double"
                          dictRef="o:totalenthalpy"
                          units="nonsi:hartree">-1339.43374202</scalar>
               </module>
               <module cmlx:templateRef="gibbs" dictRef="cc:userDefinedModule">
                  <scalar dataType="xsd:double"
                          dictRef="o:totalenthalpy"
                          units="nonsi:hartree">-1339.43374202</scalar>
                  <scalar dataType="xsd:double"
                          dictRef="o:totalentropycorr"
                          units="nonsi:hartree">-0.02942170</scalar>
                  <scalar dataType="xsd:double"
                          dictRef="o:gibbsenthalpy"
                          units="nonsi:hartree">-1339.46316372</scalar>
                  <scalar dataType="xsd:double"
                          dictRef="o:gibbsminuselec"
                          units="nonsi:hartree">-0.01372043</scalar>
               </module>
               <module cmlx:templateRef="spincontamination" dictRef="cc:userDefinedModule">
                  <scalar dataType="xsd:double" dictRef="cc:s2">8.771977</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:sideal">2.5</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:s2ideal">8.750000</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:s2deviation">0.021977</scalar>
               </module>
               <module cmlx:templateRef="spincontamination" dictRef="cc:userDefinedModule">
                  <scalar dataType="xsd:double" dictRef="cc:s2">8.774087</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:sideal">2.5</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:s2ideal">8.750000</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:s2deviation">0.024087</scalar>
               </module>
               <module cmlx:templateRef="spincontamination" dictRef="cc:userDefinedModule">
                  <scalar dataType="xsd:double" dictRef="cc:s2">8.772835</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:sideal">2.5</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:s2ideal">8.750000</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:s2deviation">0.022835</scalar>
               </module>
               <module cmlx:templateRef="spincontamination" dictRef="cc:userDefinedModule">
                  <scalar dataType="xsd:double" dictRef="cc:s2">8.773013</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:sideal">2.5</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:s2ideal">8.750000</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:s2deviation">0.023013</scalar>
               </module>
               <module cmlx:templateRef="spincontamination" dictRef="cc:userDefinedModule">
                  <scalar dataType="xsd:double" dictRef="cc:s2">8.773053</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:sideal">2.5</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:s2ideal">8.750000</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:s2deviation">0.023053</scalar>
               </module>
               <module cmlx:templateRef="spincontamination" dictRef="cc:userDefinedModule">
                  <scalar dataType="xsd:double" dictRef="cc:s2">8.773061</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:sideal">2.5</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:s2ideal">8.750000</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:s2deviation">0.023061</scalar>
               </module>
               <module cmlx:templateRef="orbitalenergies" dictRef="cc:userDefinedModule">
                  <list cmlx:templateRef="orbital">
                     <array dataType="xsd:integer" dictRef="cc:serial" size="88">0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87</array>
                     <array dataType="xsd:double" dictRef="cc:occup" size="88">1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000</array>
                     <array dataType="xsd:double"
                            dictRef="cc:energy"
                            size="88"
                            units="nonsi:electronvolt">-6985.9860 -830.8720 -717.3195 -717.1500 -716.8611 -530.5619 -107.7204 -74.4301 -74.1225 -73.6106 -33.9719 -22.0070 -20.1651 -19.9924 -19.9125 -19.6835 -18.4652 -16.2381 -16.1780 -14.3830 -9.4413 -5.5854 -5.1886 -2.2875 -2.2235 -1.8096 -1.4153 -0.0281 0.0453 0.0841 2.8730 3.6313 5.6378 8.3916 9.2292 10.4674 15.1145 15.4412 16.0217 16.2963 16.7123 19.7292 26.5266 30.7067 31.7464 32.4845 34.5445 34.6318 36.7003 39.1587 41.3359 43.9252 47.2679 50.4153 53.2450 59.1030 61.3272 63.6393 64.6615 66.5023 71.8105 72.5195 82.8528 84.0023 84.1787 85.5283 86.1019 92.7746 99.2991 99.7379 100.0682 100.7317 102.1724 102.9853 109.0592 132.3799 133.1320 133.5215 134.4979 134.7290 141.3321 144.6169 162.6864 165.5276 166.0721 174.2757 180.8605 1174.1567</array>
                     <scalar dataType="xsd:string" dictRef="o:spin">UP</scalar>
                  </list>
                  <list cmlx:templateRef="orbital">
                     <array dataType="xsd:integer" dictRef="cc:serial" size="88">0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87</array>
                     <array dataType="xsd:double" dictRef="cc:occup" size="88">1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000</array>
                     <array dataType="xsd:double"
                            dictRef="cc:energy"
                            size="88"
                            units="nonsi:electronvolt">-6985.9392 -828.6158 -715.3140 -715.2286 -715.0853 -529.7979 -102.2602 -68.2879 -68.2224 -68.0417 -31.9055 -19.9995 -16.7006 -15.5674 -15.3565 -10.4897 -10.3687 -9.9762 -9.5453 -8.8458 -8.1756 -4.4480 -4.2193 -1.6900 -1.3509 -0.8959 -0.6405 0.4089 0.6374 0.6441 3.3387 4.0821 6.4546 9.0242 10.0404 10.9975 16.1797 16.7872 17.7804 17.9503 18.2605 21.2066 27.1550 32.0548 33.1017 33.5260 35.7844 36.1441 37.6997 41.0388 42.4407 45.2966 47.8785 51.9153 54.3010 60.3835 62.4118 64.7343 66.3770 67.4978 73.5959 74.2622 84.2311 87.1556 87.9144 88.4493 89.7176 95.5260 103.8618 104.0236 104.3877 104.4968 106.3484 106.6425 112.0074 133.4952 134.2093 135.2186 135.6866 136.1069 142.7939 146.0427 164.8145 167.0426 168.0550 176.5541 182.3846 1174.8444</array>
                     <scalar dataType="xsd:string" dictRef="o:spin">DOWN</scalar>
                  </list>
               </module>
               <module cmlx:templateRef="mullikenpopulation" dictRef="cc:userDefinedModule">
                  <module cmlx:templateRef="atomicchargesspin">
                     <scalar dataType="xsd:double" dictRef="x:chargesum">1.0000000</scalar>
                     <scalar dataType="xsd:double" dictRef="x:spinsum">5.0000000</scalar>
                     <array dataType="xsd:integer" dictRef="cc:serial" size="4">0 1 2 3</array>
                     <array dataType="xsd:string" dictRef="cc:elementType" size="4">Fe O H H</array>
                     <array dataType="xsd:double" dictRef="x:charge" size="4">0.997090 -0.297553 0.271402 0.029060</array>
                     <array dataType="xsd:double" dictRef="x:spin" size="4">3.829674 0.970482 -0.015681 0.215525</array>
                  </module>
               </module>
               <module cmlx:templateRef="loewdin" dictRef="cc:userDefinedModule">
                  <array dataType="xsd:integer" dictRef="cc:serial" size="4">0 1 2 3</array>
                  <array dataType="xsd:string" dictRef="cc:elementType" size="4">Fe O H H</array>
                  <array dataType="xsd:double" dictRef="x:charge" size="4">0.886829 -0.171063 0.211684 0.072550</array>
                  <array dataType="xsd:double" dictRef="x:spin" size="4">3.851977 0.956803 0.024555 0.166664</array>
               </module>
               <module cmlx:templateRef="mayer" dictRef="cc:userDefinedModule">
                  <module cmlx:templateRef="mayercharges">
                     <array dataType="xsd:integer" dictRef="cc:serial" size="4">0 1 2 3</array>
                     <array dataType="xsd:string" dictRef="cc:elementType" size="4">Fe O H H</array>
                     <array dataType="xsd:double" dictRef="o:na" size="4">25.0029 8.2976 0.7286 0.9709</array>
                     <array dataType="xsd:double" dictRef="o:za" size="4">26.0000 8.0000 1.0000 1.0000</array>
                     <array dataType="xsd:double" dictRef="o:qa" size="4">0.9971 -0.2976 0.2714 0.0291</array>
                     <array dataType="xsd:double" dictRef="o:va" size="4">5.7224 2.2706 0.9775 0.9999</array>
                     <array dataType="xsd:double" dictRef="o:bva" size="4">2.3838 1.8720 0.9754 0.9525</array>
                     <array dataType="xsd:double" dictRef="o:fa" size="4">3.3386 0.3986 0.0021 0.0474</array>
                  </module>
                  <module cmlx:templateRef="bonds">
                     <array dataType="xsd:double" dictRef="x:distance" size="5">1.5552 0.2667 0.5618 0.3174 0.3913</array>
                     <matrix cols="2" dataType="xsd:integer" dictRef="x:serial" rows="5">0 1 0 2 0 3 1 2 2 3</matrix>
                  </module>
               </module>
               <module cmlx:templateRef="finalspenergy" dictRef="cc:userDefinedModule">
                  <scalar dataType="xsd:double" dictRef="cc:energy">-1339.452311012693</scalar>
               </module>
               <module cmlx:templateRef="electricproperties" dictRef="cc:userDefinedModule">
                  <array dataType="xsd:double" dictRef="cc:dipole" size="9">2.38421 -4.01293 -1.62873 2.77738 -2.62705 0.15033 -0.00000 0.00000 0.00000</array>
                  <scalar dataType="xsd:double" dictRef="o:magnitude" units="nonsi2:au">1.63565</scalar>
                  <scalar dataType="xsd:double" dictRef="o:magnitude" units="nonsi2:debye">4.15749</scalar>
               </module>
               <module cmlx:templateRef="thermochemistry" dictRef="cc:userDefinedModule">
                  <scalar dataType="xsd:double" dictRef="cc:temp" units="si:k">298.15</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:press" units="nonsi:atm">1.00</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:totalmass">73.86</scalar>
               </module>
               <module cmlx:templateRef="innerenergy" dictRef="cc:userDefinedModule">
                  <module cmlx:templateRef="contributions">
                     <scalar dataType="xsd:double" dictRef="o:electronic" units="nonsi:hartree">-1339.45231101</scalar>
                     <scalar dataType="xsd:double" dictRef="o:zeropoint" units="nonsi:hartree">0.01338962</scalar>
                     <scalar dataType="xsd:double"
                             dictRef="o:thermalvibcorrection"
                             units="nonsi:hartree">0.00024322</scalar>
                     <scalar dataType="xsd:double"
                             dictRef="o:thermalrotcorrection"
                             units="nonsi:hartree">0.00141627</scalar>
                     <scalar dataType="xsd:double"
                             dictRef="o:thermaltrasncorrection"
                             units="nonsi:hartree">0.00141627</scalar>
                     <scalar dataType="xsd:double"
                             dictRef="o:thermaltotal"
                             units="nonsi:hartree">-1339.43584562</scalar>
                  </module>
                  <module cmlx:templateRef="corrections">
                     <scalar dataType="xsd:double" dictRef="o:thermalcorr" units="nonsi:hartree">0.00307577</scalar>
                     <scalar dataType="xsd:double"
                             dictRef="o:nonthermalcorr"
                             units="nonsi:hartree">0.01338962</scalar>
                     <scalar dataType="xsd:double" dictRef="o:totalcorr" units="nonsi:hartree">0.01646539</scalar>
                  </module>
               </module>
               <module cmlx:templateRef="enthalpy" dictRef="cc:userDefinedModule">
                  <scalar dataType="xsd:double" dictRef="o:totalfree" units="nonsi:hartree">-1339.43584562</scalar>
                  <scalar dataType="xsd:double"
                          dictRef="o:thermalcorrenthalpy"
                          units="nonsi:hartree">0.00094421</scalar>
                  <scalar dataType="xsd:double"
                          dictRef="o:totalenthalpy"
                          units="nonsi:hartree">-1339.43490141</scalar>
               </module>
            </module>
         </module>
      </module>
   </module>
</module>
