Wednesday, August 22, 2018

physical chemistry - When can a molecule be considered freely rotating at room temperature?


This question sparked from a long discussion in chat about the nature of HX2OX2 and whether that molecule can be considered to rotate around the OO axis (and hence display axial chirality) or not.


Considering two rather clear cases:
Ethane is considered to rotate freely around the CC bond. The activation energy for rotation (equivalent to the energy difference between the staggered and eclipsed conformations) is given as 12.5kJmol.


Ethene on the other hand is considered to not freely rotate around the C=C bond. The energy difference between the planar and perpendicular conformations (the rotation barrier) is given at 250kJmol. (I am unsure of this value which I found on Yahoo answers. If I understand correctly, this should correspond with the excitation from π to π which I found elsewhere as 743kJmol.)


Somewhere between these two extreme cases must lie some kind of barrier around which rotation is severely hindered to entirely inhibited.




Attempting to calculate this myself, I remembered the Boltzmann distribution: If two states differ in energy by ΔE, their relative populations F1 and F2 can be calculated as follows:


F2F1=eΔEkBT



(F2 being the higher energy state and kB being the Boltzmann constant.)


If I calculate relative populations from the energy differences 12.5,19,25 and 250kJmol for ethane, butane, HX2OX2 and ethene, respectively (which correspond to 0.13,0.20,0.26 and 2.6eVparticle) I get the following results at 300K:



  • Ethane: 6.55×103

  • Butane: 4.79×104

  • HX2OX2: 4.29×105

  • Ethene: 2.09×1044


Or logarithmic (ln) values of:




  • Ethane: 5.03

  • Butane: 7.64

  • HX2OX2: 10.1

  • Ethene: 101


Again, somewhere between butane and ethene must be a threshold-ish value, above which free rotation at 300K cannot be considered, but what range are we looking at; similar to what energy difference does this correspond to?


Asked the other way round, if I have an energy barrier of 25kJmol, what would I need to put on into the equation to calculate the temperature above which free rotation can be assumed?



Answer



I'm not a transition-state physical chemist, but I think a good approach to this problem is transition-state theory, specifically the Eyring equation:


k=κkbTheΔGRT



This equation tries to predict the rate constant k from an assumed pseudo-equilibrium between the transition-state and the starting material, which are assumed to differ in energy by ΔG. kb is Boltzmann's constant, h is Planck's constant, T is temperature, and κ is the "fudge factor" or transmission coefficient, the fraction of transition states that react productively instead of converting back to starting material.


In the conformational reactions you are talking about, I think a good starting place would be to assume the eclipsed conformations are the transition state, although I don't know if this has been shown to be 100% true experimentally or computationally. As a first approximation, let's assume there isn't a huge entropy term in ΔG, so that it is well-approximated by the energy level differences you found. If we also assume that κ is 1.0, then we can just multiply the eΔERT you've already calculated by kbTh=6.2×1012 s1 at 298 K to get (very) approximate rate constants.



  • Ethane: 6.55×1034.1×1010 s1

  • Butane: 4.79×1043.1×109 s1

  • HX2OX2: 4.29×1052.0×108 s1

  • Ethene: 2.09×10441.3×1031 s1


That shows that rotational equilibrium of everything except ethene should be very rapid, but that the half-life for ethene equilibration much greater than the age of the universe (4×1017 s).


Suppose we want to find an energy difference that should lead to a rate constant that is neither very fast nor very slow, say we'd like k1 hr1. Inverting our very very crude approximation gives:



ΔE=RTln(khκkbT)


Plugging in k=1 hr1=2.7×104s1 gives (if I did my math right) 93 kJ per mol. Some well-known atropoisomers do have estimated energy differences in that range.


No comments:

Post a Comment

periodic trends - Comparing radii in lithium, beryllium, magnesium, aluminium and sodium ions

Apparently the of last four, MgX2+ is closest in radius to LiX+. Is this true, and if so, why would a whole larger shell ($\ce{...