The kinetics of recombination of the P(+)H(A)(-) radical pair were compared in wild-type reaction centers from Rhodobacter sphaeroides and in seven mutants in which the free energy gap, ΔG, between the charge separated states P(+)B(A)(-) and P(+)H(A)(-) was either increased or decreased. Five of the mutant RCs had been described previously, and X-ray crystal structures of two newly constructed complexes were determined by X-ray crystallography. The charge recombination reaction was accelerated in all mutants with a smaller ΔG than in the wild-type, and was slowed in a mutant having a larger ΔG. The free energy difference between the state P(+)H(A)(-) and the PH(A) ground state was unaffected by most of these mutations. These observations were consistent with a model in which the P(+)H(A)(-) → PH(A) charge recombination is thermally activated and occurs via the intermediate state P(+)B(A)(-), with a mean rate related to the size of the ΔG between the states P(+)B(A)(-) and P(+)H(A)(-) and not the ΔG between P(+)H(A)(-) and the ground state. A more detailed analysis of charge recombination in the mutants showed that the kinetics of the reaction were multiexponential, and characterized by ~0.5, ~1-3, and 7-17 ns lifetimes, similar to those measured for wild-type reaction centers. The exact lifetimes and relative amplitudes of the three components were strongly modulated by the mutations. Two models were considered in order to explain the observed multiexponentiality and modulation, involving heterogeneity or relaxation of P(+)H(A)(-) states, with the latter model giving a better fit to the experimental results.