Negative pH, efflorescent mineralogy, and consequences for environmental restoration at the Iron Mountain Superfund site, California

The Richmond Mine of the Iron Mountain copper deposit contains some of the most acid mine waters ever reported. Values of pH have been measured as low as −3.6, combined metal concentrations as high as 200 g/liter, and sulfate concentrations as high as 760 g/liter. Copious quantities of soluble metal sulfate salts such as melanterite, chalcanthite, coquimbite, rhomboclase, voltaite, copiapite, and halotrichite have been identified, and some of these are forming from negative-pH mine waters. Geochemical calculations show that, under a mine-plugging remediation scenario, these salts would dissolve and the resultant 600,000-m3 mine pool would have a pH of 1 or less and contain several grams of dissolved metals per liter, much like the current portal effluent water. In the absence of plugging or other at-source control, current weathering rates indicate that the portal effluent will continue for approximately 3,000 years. Other remedial actions have greatly reduced metal loads into downstream drainages and the Sacramento River, primarily by capturing the major acidic discharges and routing them to a lime neutralization plant. Incorporation of geochemical modeling and mineralogical expertise into the decision-making process for remediation can save time, save money, and reduce the likelihood of deleterious consequences.

Dynamical hierarchy in transition states: Why and how does a system climb over the mountain?

How a reacting system climbs through a transition state during the course of a reaction has been an intriguing subject for decades. Here we present and quantify a technique to identify and characterize local invariances about the transition state of an N-particle Hamiltonian system, using Lie canonical perturbation theory combined with microcanonical molecular dynamics simulation. We show that at least three distinct energy regimes of dynamical behavior occur in the region of the transition state, distinguished by the extent of their local dynamical invariance and regularity. Isomerization of a six-atom Lennard–Jones cluster illustrates this: up to energies high enough to make the system manifestly chaotic, approximate invariants of motion associated with a reaction coordinate in phase space imply a many-body dividing hypersurface in phase space that is free of recrossings even in a sea of chaos. The method makes it possible to visualize the stable and unstable invariant manifolds leading to and from the transition state, i.e., the reaction path in phase space, and how this regularity turns to chaos with increasing total energy of the system. This, in turn, illuminates a new type of phase space bottleneck in the region of a transition state that emerges as the total energy and mode coupling increase, which keeps a reacting system increasingly trapped in that region.