13

Detailed balance is an important property of many classes of physical systems. It can be written as $$ \frac{p_{i \to j}}{p_{j \to i}} = e^{\frac{\Delta G}{k_B T}},\tag{1} $$ where $i$ and $j$ represent microscopic states of the entire system; $p$ represents the probability for the system to transition from one state to another in a particular finite time period; and $\Delta G$ represents the difference in free energy between the two states. (Whether to use the Gibbs or Helmholtz free energy or some other potential depends on the ensemble.)

Many systems obey detailed balance, but not all do. The Earth cannot fluctuate backwards in its orbit around the Sun, because this would violate conservation of angular momentum. An RLC circuit does not obey detailed balance because the fluctuations have a distinctive ringing time. For these systems the correct formula is $$ \frac{p_{i \to j}}{p_{j' \to i'}} = e^{\frac{\Delta G}{k_B T}},\tag{2} $$ where $i'$ and $j'$ represent states identical to $i$ and $j$, except that all velocities and magnetic fields have been reversed. (In quantum mechanics, they represent something like the complex conjugates of states $i$ and $j$.)

Both of these formulae guarantee that the system will obey the second law (on average), but $(1)$ is substantially stronger, because it guarantees that not only will the system tend toward equilibrium in the thermodynamic limit, but that it will not oscillate as it approaches the equilibrium. (It can still oscillate far away from equilibrium, however.)

My question is about chemical kinetics. Here we universally assume equation $(1)$ and not $(2)$. This puts strong constraints on the reaction rates, and leads to the well-known result that near-equilibrium oscillations are impossible in chemical systems. I've recently been discussing this topic with a very experienced researcher in nonlinear dynamics, and I found myself unable to convince him that $(1)$ rather than $(2)$ is a good assumption in the case of chemistry.

So I thought I'd ask here and see if anyone can help me out: what is the argument that leads us to assume the 'strong' form of detailed balance in chemical systems, rather than the weaker form in equation $(2)$?

N. Virgo
  • 35,274

1 Answers1

2

For physicists, the term "detailed balance" is used in one of the following contexts:

1) A closed system A is observed at a coarse-grained level (e.g. A is a vessel with a macroscopic amount of chemicals and one only tracks/observes the global concentrations of the chemicals).

2) An observable system A (e.g. a mixture of chemicals) is weakly in contact with an external "bath" B. The latter is in thermodynamic equilibrium.

Microscopically then, the contents of A resp. A+B must obey deterministic, time-reversible equations of motion. Careful: when we hit the rewind-button for the motion of a particle in time, its instantaneous position $x$ is not changed, but its velocity $v$ of course flips to $-v$. The former variables are called "odd" and the latter are "even" under time-exchange. Now, when the global system ($A$ resp. $A+B$) is in thermal equilibrium, every microscopic trajectory has the same likelihood as the time-reversed trajectory. So microscopically the probability to see a certain transition-event from a micro-state $\alpha$ to $\beta$ is equal to the probability to see a transition from $\beta$ to $\alpha$. For our macroscopic observations (concerning $P_{i \to j}$, it remains simply to count the number $W(i)$ of micro-states $\alpha$ that correspond to a certain macrostate $i$. The $\log$ of this microstate-multiplicity of the macrostate $i$ is precisely the entropy/ free energy (depending on the context) $G(i)$. Also there's the multiplication of the log by $k_B T$. We then obtain $$1=\frac{W(i)P_{i \to j}}{W(j')P_{j' \to i'}}=\frac{W(i)P_{i \to j}}{W(j)P_{j' \to i'}}.$$ So $$\frac{P_{i \to j}}{P_{j' \to i'}}=\frac{W(j)}{W(i)}=e^{\frac{G(j)-G(i)}{k_BT}}.$$ To proceed finally to your question why in the context of chemical kinetics, the time-reversed state $i'$ of $i$ is simply $i$: You presumably keep track of the macroscopic concentrations of your chemicals and those data are the content of your states {i}. Now, since the microscopic positions of particles are even under time-exchange, also concentrations of those particles in a fixed region of space are even variables. So indeed, upon applying time-reversal, an array of concentrations $i$ remains the same array of concentrations $i$.

Hecatonchires
  • 633
  • 3
  • 12