Skip to content

Bond Dissociation Energy

This workflow computes the homolytic bond dissociation energies of any hydrogen, flourine, chlorine, bromine, and iodine atoms in a molecule.

Steps

The workflow performs the following steps:

  1. Identify atoms for which the dissociation energy should be computed (e.g., H, F, Cl, Br, and I atoms).
  2. Optimize the geometry of the original molecule.
  3. For each dissociating atom, compute the geometrically relaxed energy of the molecule upon removing that atom.
  4. Compute the single-atom energies of any dissociating atoms.
  5. Evaluate the dissociation energies as the difference between the dissociated fragments and the original molecule.

Examples

The following example evaluates the bond dissociation energies of a halogenated methane molecule.

import sierra
from sierra.inputs import *

halogenated_methane = Molecule(
    data="""
    Cl   0.8836    0.0000    0.0000
    C   -0.8836    0.0000    0.0000
    F   -1.2405   -0.6985    0.7600
    H   -1.2405    1.0074    0.2249
    H   -1.2406   -0.3089   -0.9849
    """
)

inp = BondDissociationEnergyInput(
    molecule=halogenated_methane, energy_method=XTBMethod(model="GFN1")
)
result = sierra.run(inp)

print("Bond dissociation energies:")
#> Bond dissociation energies:
for atom_index, dissociation_energy in result.dissociation_energies.items():
    symbol = result.molecule.symbols[atom_index]
    print(f"Bond C-{symbol}: {dissociation_energy} Hartree")
    #> Bond C-Cl: 0.17975568409618425 Hartree
    #> Bond C-F: 0.23445643098841273 Hartree
    #> Bond C-H: 0.22056541588326617 Hartree
    #> Bond C-H: 0.22056541234047167 Hartree

BondDissociationEnergyInput

Model for describing input for the BDE workflow.

Fields

energy_method

The method used to compute the energies.

molecule

The molecule for which dissociation energies are to be computed.

workflow
  • Type: WorkflowIdentifier
  • Default: WorkflowIdentifier(name='bonddissociationenergy', image=None)

BondDissociationEnergyResult

Model for describing results from the BDE workflow.

Fields

All the fields in BondDissociationEnergyInput and the following:

computations
  • Type: Mapping[str, Any]
  • Default: {}
dissociation_energies
  • Type: Mapping[int, float]
  • Default: {}