Full Text Available

Note: Clicking the button above will open the full text document at the original institutional repository in a new window.

Development of a group contribution method to predict the enthalpy of formation of energetic azoles

Energetic materials (EMs) are a class of high-nitrogen content material that stores a large amount of chemical energy released upon external stimuli (e.g., friction, thermal shock, or electrostatic discharge). Azoles, which are five-membered heterocyclic aromatic com- pounds containing a nitrogen at...

Full description

Saved in:
Bibliographic Details
Main Author: Coetzee, Megan
Other Authors: Venter, Gerhard
Format: Thesis
Language:English
English
Published: Department of Chemistry 2026
Subjects:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1869483664486694912
access_status_str Open Access
author Coetzee, Megan
author2 Venter, Gerhard
author_browse Coetzee, Megan
Venter, Gerhard
author_facet Venter, Gerhard
Coetzee, Megan
author_sort Coetzee, Megan
collection Thesis
description Energetic materials (EMs) are a class of high-nitrogen content material that stores a large amount of chemical energy released upon external stimuli (e.g., friction, thermal shock, or electrostatic discharge). Azoles, which are five-membered heterocyclic aromatic com- pounds containing a nitrogen atom and at least one other heteroatom, are an ideal source of EMs. The enthalpy of formation (∆fH) is a critical thermodynamic quantity that is required to estimate the detonation performance of EMs (e.g., detonation pressure and velocity). Because physical measurements of ∆fH require time and resources, group con- tribution methods (GCMs) provide an alternative that is quick and cost-effective. A GCM is built on the principle that a property of a molecule can be estimated using the sum of individual contributions associated with its smaller structural units, or groups. The enthalpy of formation of a cyclic compound can be obtained by summing the corre- sponding group additive values (GAV) determined from the acyclic reference molecules and adding a correction accounting for the ring strain (RS). This ring-strain correction (RSC) is generally positive and is calculated from the difference between the known value of ∆fH of the cyclic compound and the sum of the acyclic GAVs. However, if ∆fH of an aromatic compound is calculated in the same way, the resulting correction must account for both the RS and the aromatic stabilisation energy (ASE), where this sum is typically negative. Alternatively, an aromatic compound may also be built directly as a sum over GAVs that have been determined using aromatic reference molecules. The former model has been applied to pyrrole derivatives only, and to the best of our knowledge, although the latter model has been applied to azines, it has not been attempted for azoles. In this work, GCMs based on both approaches were developed using ab initio quantum mechanical (QM) calculations due to the limited availability of physical measurements of ∆fH of azoles. The uncertainty arising primarily from the systematic error or bias associated with the G4(MP2)-6X composite method was first quantified using 47 neutral CHN-containing molecules. Reference ∆fH values for these molecules were obtained from Active Thermochemical Tables (ATcT). A correction of +1.11 kJ mol−1 resulted, while the standard uncertainty associated with this correction is 3.04 kJ mol−1, leading to a 95 % confidence interval of 6.08 kJ mol−1 for the calculated values. Subsequently, a database of 72 acyclic molecules and multiple linear regression (MLR) was used to determine 42 Benson-type GAVs, which were then applied to estimate ∆fH of ten unsubstituted azoles, including pyrrole, two diazoles, four triazoles, two tetrazoles, and pentazole. The sum of RS and ASE was determined as the difference between the calculated ∆fH values and the sum of the acyclic GAVs, providing a GCM based on the strain-centred approach. To extend the work to energetic azoles, the aromatic group approach was then applied and a set of 33 GAVs was obtained using MLR of a database of the calculated ∆fH values of 60 energetic azoles. In addition to unsubstituted azoles, these azoles were also functionalised with a methyl group and explosophoric functional groups such as amine, azide, nitro, and nitramide. The resulting GCM showed a mean absolute error (MAE) of 3.22 kJ mol−1 and maximum error of 10.95 kJ mol−1 for the enthalpy of formation across all functionalised and unfunctionalised azoles.
format Thesis
id oai:open.uct.ac.za:11427/43377
institution University of Cape Town (South Africa)
language English
eng
last_indexed 2026-07-01T04:02:35.762Z
license_str Not specified — see source repository
provenance_str_mv Harvested via OAI-PMH from UCTD — University of Cape Town Open Access Repository
publishDate 2026
publishDateRange 2026
publishDateSort 2026
publisher Department of Chemistry
publisherStr Department of Chemistry
record_format dspace
source_str UCTD — University of Cape Town Open Access Repository
spelling oai:open.uct.ac.za:11427/43377 Development of a group contribution method to predict the enthalpy of formation of energetic azoles Coetzee, Megan Venter, Gerhard Energetic materials group additive values Energetic materials (EMs) are a class of high-nitrogen content material that stores a large amount of chemical energy released upon external stimuli (e.g., friction, thermal shock, or electrostatic discharge). Azoles, which are five-membered heterocyclic aromatic com- pounds containing a nitrogen atom and at least one other heteroatom, are an ideal source of EMs. The enthalpy of formation (∆fH) is a critical thermodynamic quantity that is required to estimate the detonation performance of EMs (e.g., detonation pressure and velocity). Because physical measurements of ∆fH require time and resources, group con- tribution methods (GCMs) provide an alternative that is quick and cost-effective. A GCM is built on the principle that a property of a molecule can be estimated using the sum of individual contributions associated with its smaller structural units, or groups. The enthalpy of formation of a cyclic compound can be obtained by summing the corre- sponding group additive values (GAV) determined from the acyclic reference molecules and adding a correction accounting for the ring strain (RS). This ring-strain correction (RSC) is generally positive and is calculated from the difference between the known value of ∆fH of the cyclic compound and the sum of the acyclic GAVs. However, if ∆fH of an aromatic compound is calculated in the same way, the resulting correction must account for both the RS and the aromatic stabilisation energy (ASE), where this sum is typically negative. Alternatively, an aromatic compound may also be built directly as a sum over GAVs that have been determined using aromatic reference molecules. The former model has been applied to pyrrole derivatives only, and to the best of our knowledge, although the latter model has been applied to azines, it has not been attempted for azoles. In this work, GCMs based on both approaches were developed using ab initio quantum mechanical (QM) calculations due to the limited availability of physical measurements of ∆fH of azoles. The uncertainty arising primarily from the systematic error or bias associated with the G4(MP2)-6X composite method was first quantified using 47 neutral CHN-containing molecules. Reference ∆fH values for these molecules were obtained from Active Thermochemical Tables (ATcT). A correction of +1.11 kJ mol−1 resulted, while the standard uncertainty associated with this correction is 3.04 kJ mol−1, leading to a 95 % confidence interval of 6.08 kJ mol−1 for the calculated values. Subsequently, a database of 72 acyclic molecules and multiple linear regression (MLR) was used to determine 42 Benson-type GAVs, which were then applied to estimate ∆fH of ten unsubstituted azoles, including pyrrole, two diazoles, four triazoles, two tetrazoles, and pentazole. The sum of RS and ASE was determined as the difference between the calculated ∆fH values and the sum of the acyclic GAVs, providing a GCM based on the strain-centred approach. To extend the work to energetic azoles, the aromatic group approach was then applied and a set of 33 GAVs was obtained using MLR of a database of the calculated ∆fH values of 60 energetic azoles. In addition to unsubstituted azoles, these azoles were also functionalised with a methyl group and explosophoric functional groups such as amine, azide, nitro, and nitramide. The resulting GCM showed a mean absolute error (MAE) of 3.22 kJ mol−1 and maximum error of 10.95 kJ mol−1 for the enthalpy of formation across all functionalised and unfunctionalised azoles. 2026-06-25T08:18:45Z 2026-06-25T08:18:45Z 2026 2026-06-25T07:47:03Z Thesis / Dissertation Masters MSc http://hdl.handle.net/11427/43377 en eng application/pdf Department of Chemistry Faculty of Science University of Cape Town
spellingShingle Energetic materials
group additive values
Coetzee, Megan
Development of a group contribution method to predict the enthalpy of formation of energetic azoles
thesis_degree_str Master's
title Development of a group contribution method to predict the enthalpy of formation of energetic azoles
title_full Development of a group contribution method to predict the enthalpy of formation of energetic azoles
title_fullStr Development of a group contribution method to predict the enthalpy of formation of energetic azoles
title_full_unstemmed Development of a group contribution method to predict the enthalpy of formation of energetic azoles
title_short Development of a group contribution method to predict the enthalpy of formation of energetic azoles
title_sort development of a group contribution method to predict the enthalpy of formation of energetic azoles
topic Energetic materials
group additive values
url http://hdl.handle.net/11427/43377
work_keys_str_mv AT coetzeemegan developmentofagroupcontributionmethodtopredicttheenthalpyofformationofenergeticazoles