Skip to content

thermochemistry ​

Implementation level ​

TypeStatus
CML extraction template
HTML5 representation

Template attributes ​

AttributeValue
sourceADF log
idthermochemistry
pattern\s*Statistical\sThermal\sAnalysis.*
endPattern\s*\*{10}+\s*
endPattern2\s*={10,}+\s*$\s*\S+.*$\s*={10,}+
endPattern3~
endOffset0
xml:basefrequencyanalysis/thermochemistry.xml

Input ​

 ============================
 Statistical Thermal Analysis  ***  ideal gas assumed  ***
 ============================
  
 Pressure:                  1.000000 atm.
 Temperature:             300.000000 K



 Moments of Inertia (and direction vectors)
 ==========================================

       41342.8675      42584.1405      47497.7524
 ------------------------------------------------
           0.9979          0.0640          0.0012
           0.0640         -0.9973         -0.0372
           0.0012         -0.0372          0.9993


 The rotational contribution to the molecular entropy includes
 a term, dependent on the symmetry number sigma. The results 
 reported below were computed using sigma = 1, determined
 from the point group symmetry of the input geometry (NOSYM).
 If this is not the correct symmetry, please contact SCM to 
 report a bug.


     Temp                                                       Transl     Rotat    Vibrat     Total
     ----                                                       ------     -----    ------     -----

     300.00   Entropy (cal/mole-K):                             47.312    39.674   110.986   197.972
              Internal Energy (Kcal/mole):                       0.894     0.894   335.070   336.859
              Constant Volume Heat Capacity (cal/mole-K):        2.981     2.981   179.493   185.455
 
 ************************************************************************************************

Output text ​

xml
<comment class="example.output" id="thermochemistry"> 
        <module cmlx:templateRef="thermochemistry">       
            <scalar dataType="xsd:double" dictRef="cc:press" units="nonsi:atm">1.0</scalar>
            <scalar dataType="xsd:double" dictRef="cc:temp" units="si:k">300.0</scalar>
            <array dataType="xsd:double" dictRef="cc:moi" size="3">41342.8675 42584.1405 47497.7524</array>         
            <scalar dataType="xsd:integer" dictRef="cc:symmnumber">1</scalar>
            <scalar dataType="xsd:string" dictRef="cc:pointgroup">NOSYM</scalar>
            <module cmlx:lineCount="6" cmlx:templateRef="energies">
               <scalar dataType="xsd:double" dictRef="cc:temp" units="si:k">300.0</scalar>
               <list cmlx:templateRef="entropy">               
                  <scalar dataType="xsd:double" dictRef="cc:transl" units="nonsi2:cal.mol-1.K-1">47.312</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:rotat" units="nonsi2:cal.mol-1.K-1">39.674</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:vibrat" units="nonsi2:cal.mol-1.K-1">110.986</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:total" units="nonsi2:cal.mol-1.K-1">197.972</scalar>
               </list>
               <list cmlx:templateRef="internalEnergy">
                  <scalar dataType="xsd:double" dictRef="cc:transl" units="nonsi2:kcal.mol-1">0.894</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:rotat" units="nonsi2:kcal.mol-1">0.894</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:vibrat" units="nonsi2:kcal.mol-1">335.07</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:total" units="nonsi2:kcal.mol-1">336.859</scalar>
               </list>
               <list cmlx:templateRef="heat">
                  <scalar dataType="xsd:double" dictRef="cc:transl" units="nonsi2:cal.mol-1.K-1">2.981</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:rotat" units="nonsi2:cal.mol-1.K-1">2.981</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:vibrat" units="nonsi2:cal.mol-1.K-1">179.493</scalar>
                  <scalar dataType="xsd:double" dictRef="cc:total" units="nonsi2:cal.mol-1.K-1">185.455</scalar>
               </list>
            </module>
        </module> 
    </comment>

Template definition ​

xml
<templateList>  <template id="temppressure" pattern="\s*Pressure:.*" endPattern="\s*Temperature:.*" endOffset="1">    <record id="pressure">\s*Pressure:&#123;F,cc:press&#125;.*</record>    <record id="temperature">\s*Temperature:&#123;F,cc:temp&#125;.*</record>    <transform process="pullup" xpath=".//cml:scalar" />
        </template>  <template id="inertia" pattern="\s*Moments\sof\sInertia.*" endPattern=".*\d\s*$\s*" endOffset="2">    <record repeat="3" />    <record>&#123;3F,cc:moi&#125;</record>                   
        </template>  <template id="symmetry" pattern="\s*The\srotational\scontribution.*" endPattern="\s*" endPattern2="~">    <record repeat="2">.*</record>    <record>.*sigma\s=\s&#123;I,cc:symmnumber&#125;.*</record>    <record>.*point\sgroup\ssymmetry\sof\sthe\sinput\sgeometry\s\(&#123;X,cc:pointgroup&#125;\).*</record>    <transform process="pullup" xpath=".//cml:scalar" />
        </template>  <template id="energies" pattern="\s*Temp\s*Transl.*" endPattern="\s*Constant\sVolume\sHeat.*" endOffset="1" repeat="*">    <record repeat="3" />    <record id="entropy">&#123;F,cc:temp&#125;Entropy\s*\(cal/mole-K\):&#123;F,cc:transl&#125;&#123;F,cc:rotat&#125;&#123;F,cc:vibrat&#125;&#123;F,cc:total&#125;</record>    <record id="internalEnergy">\s*Internal\sEnergy\s\(Kcal/mole\):&#123;F,cc:transl&#125;&#123;F,cc:rotat&#125;&#123;F,cc:vibrat&#125;&#123;F,cc:total&#125;</record>    <record id="heat">\s*Constant\sVolume\sHeat\sCapacity\s\(cal/mole-K\):&#123;F,cc:transl&#125;&#123;F,cc:rotat&#125;&#123;F,cc:vibrat&#125;&#123;F,cc:total&#125;</record>
        </template>   
    </templateList>
<transform process="pullup" xpath=".//cml:scalar" />
<transform process="pullup" xpath=".//cml:array" repeat="2" />
<transform process="addUnits" xpath=".//cml:list[@cmlx:templateRef='entropy' or @cmlx:templateRef='heat']/cml:scalar" value="nonsi2:cal.mol-1.K-1" />
<transform process="addUnits" xpath=".//cml:list[@cmlx:templateRef='internalEnergy']/cml:scalar" value="nonsi2:kcal.mol-1" />
<transform process="addUnits" xpath="./cml:scalar[@dictRef='cc:press']" value="nonsi:atm" />
<transform process="addUnits" xpath=".//cml:scalar[@dictRef='cc:temp']" value="si:k" />
<transform process="pullup" xpath=".//cml:list[@cmlx:templateRef='entropy']/cml:scalar[@dictRef='cc:temp']" />
<transform process="delete" xpath=".//cml:list[count(*)=0]" />
<transform process="delete" xpath=".//cml:list[count(*)=0]" />
<transform process="delete" xpath=".//cml:module[count(*)=0]" />