RMG-Java icon indicating copy to clipboard operation
RMG-Java copied to clipboard

Rotor counting inaccurate for certain species with triple bonds

Open gmagoon opened this issue 14 years ago • 1 comments

RMG will incorrectly count rotors for certain species with triple bonds, which will affect the accuracy of Cp,Inf used in Wilhoit polynomial fitting. Examples include:

  1. species where "internal" bond rotation corresponds to rotation of the entire molecule: e.g. CH3CCH (should be 0 rotors, but RMG will count the H3C-C rotor)
  2. species with "degenerate" internal rotations: e.g. dimethylacetylene: CH3CCCH3 (should be 1 rotor, but RMG will count both H3C-C rotors)

Relevant section of code is ChemGraph's calculateInternalRotor().

gmagoon avatar Apr 19 '10 03:04 gmagoon

A related potential issue involves resonant radicals like allyl (H2C=CH-CH2.). RMG will count one rotor, but one could make an argument for 0 or 2 rotors (depending on the temperature range of interest...I'm not sure of the barrier off the top of my head).

gmagoon avatar Apr 19 '10 16:04 gmagoon