  We investigate the atmospheric oxidation of mercury Hg(0) by halogens,
initiated by Br and I to yield Hg(I), and continued by I, Br, BrO, ClO, IO, NO2
and HO2 to yield Hg(II) or Hg(0) using computational methods with focus on
determining the impact of rising iodine levels. We calculate reaction
enthalpies and Gibbs free energies using the Coupled Cluster singlets,
doublets, and perturbative triplets method (CCSD(T)) with the ma-def2-TZVP
basis set and effective core potential to account for relativistic effects.
Additionally, we investigate the reaction kinetics using variational transition
state theory based on geometric scans of bond dissociations at the
CASPT2/ma-def2-TZVP level. We compare the results obtained from the CASPT2 and
CCSD(T) methods to help define the uncertainty. Our results provide insights
into the mechanisms of these reactions and their implications for mercury
depletion events and for the atmosphere as a whole. The reaction HgBr + Br ->
HgBr2 was found to be twice as fast as HgI + I -> HgI2, with reaction rate
coefficients of 8.8x10^-13 and 4.2x10^-13 cm^3 molecule^-1 s^-1 respectively.
The BrHg + BrO -> BrHgOBr reaction was about 7.2 times faster than the HgI + IO
-> IHgOI reaction with their rates being 3.3x10^-14 and 4.6x10^-15 cm^3
molecule^-1 s^-1 respectively. We investigate the HgXOY (X and Y=halogen)
complexes. We find that rising iodine levels will lead to shortened mercury
lifetime due to the impact of the HgI + I -> HgI2 reaction.
