StatisticalThermodynamics - Maple Help
For the best experience, we recommend viewing online help using Google Chrome or Microsoft Edge.

Online Help

All Products    Maple    MapleSim


Home : Support : Online Help : Toolboxes : Quantum Chemistry : Courses : Lessons : StatisticalThermodynamics

Statistical Thermodynamics: Calculating Thermodynamic Functions DH, DG, and DS

Copyright (c) RDMCHEM LLC 2019

 

 

Overview

Initialize

Specify Molecule

Calculate Thermodynamics

Calculate ΔH, ΔG, and ΔS (after completing above for all reactants and products)

Appendix

Overview

 

In this activity, you will apply statistical thermodynamics to calculate the enthalpy, free energy, and entropy of reaction for the combustion of ethene to form ethylene oxide at 298 K:

 

 

Rxn (1)

 

To simplify the calculations, you will treat each species as an ideal gas. The thermal entropy (S) (J/mol K) and the internal energy (E) (J/mol) are given by:

 

(1)

 

and

(2)

 

where q is the molecular partition function  , where the subscript t, e, r, and v refer to the translational, electronic, rotational, and vibrational contributions, respectively.  To calculate q, we further assume each species is in the ground electronic state (with degeneracy 2 spin + 1) with rovibrational energies given by harmonic oscillator and rigid rotor approximations. These assumptions lead to the following expressions for q, S, and E:

 

Translational:

(3)

(4)

(5)

where P = 1 atm = 100.325 kPa.

 

Electronic:

(6)

(7)

 

The thermal contribution to internal energy is 0 since the  does not depend on temperature.

 

Rotational:

Linear Molecule:

 

(8)

(9)

(10)

where  is the rotational symmetry number and  is the rotational temperature.  Each molecule we will consider in this activity is characterized by   

Nonlinear Polyatomic Molecule:

 

(11)

(12)

(13)

 

where  are the rotational temperatures corresponding to the principal moments of inertia.

 

Vibrational:

(14)

(15)

(16)

 

where s = 3N - 5 normal modes for a linear molecule and s = 3N - 6 for a nonlinear polyatomic.

 

Once the translational, electronic, rotational, and vibrational thermal contributions to the total entropy and internal energy have been calculated, total entropy (S), internal energy (E), enthalpy (H), and free-energy (G) can be calculated as follows:

 

(17)

(18)

(19)

(20)

 

Note we have corrected for the vibrational zero-point-energy for the internal energy.  We will calculate S, H, and G for each of the reactants and product in Rxn (1) and calculate ΔH, ΔG, and ΔS as

 

(21)

etc.

 

Initialize

 

We begin with initializing the QuantumChemistry package, along with the ScientificConstants and LinearAlgebra packages.

 

 

Specify Molecule

 

Here we specify the coordinates of the molecule for which we want to calculate thermodynamics as well as some related calculate parameters, such as symmetry number, charge, AObasis, etc.  Note that the geometries provided in .xyz data files are at the HartreeFock / STO-3G level.   If you wish to use a different method/basis combination, begin with the provided .xyz coordinates for each molecule and uncomment the geometry optimization line below.

 

(3.1)

Calculate Thermodynamics

 

If  not using Hartree-Fock and sto-3g, uncomment the following line:

 

 

Calculate electronic energy in J/mol using AOmethod, AObasis, spin, and charge defined above: 

 

(4.1)

 

 

To calculate the vibrational energies, perform a vibrational analysis:

 

(4.2)

 

 

Calculate principal moments of inertia:

 

(4.3)

 

 

Now we are ready to calculate partition functions!

 

Translational partition function:

 

(4.4)

 

Electronic partition function:

 

(4.5)

 

 

Rotational temperatures and partition function:

 

(4.6)

 

Vibrational temperatures and partition function:

 

(4.7)

 

ZPE correction in J/mol:

 

(4.8)

 

 

Now we can calculate total S, E, H, and G according to the equations provided in the Overview:

 

(4.9)

 

Repeat the above calculations for each reactant and product in Rxn (1). Don't forget to change the molec_index!

Calculate DH, DG, and DS (after completing above for all reactants and products)

 

 

(5.1)

(5.2)

(5.3)

 

 

NIST lists the following values for enthalpy of formation and free-energy of formation:

 

 

 (kJ/mol)

 (kJ/mol) (calc.)

S (1 bar) (J/molK)

ethylene

52.4

61.4

219.3

ethylene oxide

-52.6

-51.7

243.0

oxygen

0

0

205.15

   

Using Hess's Law, calculate ΔH, ΔG, and ΔS for Rxn (1) and compare to your calculated results.

 

 

(5.4)

(5.5)

(5.6)

(5.7)

 

Appendix

References

1) D. A. McQuarrie and J. D. Simon, Physical Chemistry: A Molecular Approach (University Science, New York, 1997).

2) JPCA scaling factors 2007

3) https://webbook.nist.gov/cgi/inchi/InChI=1S/C2H4/c1-2/h1-2H2

4) https://webbook.nist.gov/cgi/inchi/InChI=1S/C2H4O/c1-2-3-1/h1-2H2


Download Help Document