Next Article in Journal
An Investigation into the Potential of Targeting Escherichia coli rne mRNA with Locked Nucleic Acid (LNA) Gapmers as an Antibacterial Strategy
Next Article in Special Issue
Modeling Chemical Reactivity at the Interfaces of Emulsions: Effects of Partitioning and Temperature
Previous Article in Journal
Characterizing the Volatile and Sensory Profiles, and Sugar Content of Beeswax, Beebread, Bee Pollen, and Honey
Previous Article in Special Issue
The Isoelectric Point of an Exotic Oxide: Tellurium (IV) Oxide
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Methyl-Cyclohexane Methanol (MCHM) Isomer-Dependent Binding on Amorphous Carbon Surfaces

by
William A. Alexander
Department of Chemistry, The University of Memphis, Memphis, TN 38152, USA
Molecules 2021, 26(11), 3411; https://doi.org/10.3390/molecules26113411
Submission received: 8 May 2021 / Revised: 25 May 2021 / Accepted: 28 May 2021 / Published: 4 June 2021
(This article belongs to the Special Issue Interfacial Chemistry)

Abstract

:
In January 2014, over 10,000 gallons of methyl-cyclohexane methanol (MCHM) leaked into the Elk River in West Virginia, in a chemical spill incident that contaminated a large portion of the state’s water supply and left over 300,000 residents without clean water for many days and weeks. Initial efforts to remove MCHM at the treatment plant centered on the use of granulated activated carbon (GAC), which removed some of the chemical from the water, but MCHM levels were not lowered to a “non-detect” status until well after the chemical plume had moved downstream of the intake. Months later, MCHM was again detected at the outflow (but not the inflow) at the water treatment facility, necessitating the full and costly replacement of all GAC in the facility. The purpose of this study is to investigate the hypothesis that preferential absorbance of one of the two MCHM isomers, coupled with seasonal variations in water temperature, explain this contrary observation. Calculated intermolecular potentials between ovalene (a large planar polycyclic aromatic hydrocarbon) and the MCHM isomers were compared to physisorption potentials of MCHM onto an amorphous carbon model. While a molecular mechanics (MM) force field predicts no difference in the average interaction potentials between the cis- and trans-MCHM with the planar ovalene structure, MM predicts that the trans isomer binds stronger than the cis isomer to the amorphous carbon surface. Semi-empirical and density functional theory also predict stronger binding of trans-MCHM on both the planar and amorphous surfaces. The differences in the isomer binding strengths on amorphous carbon imply preferential absorbance of the trans isomer onto activated charcoal filter media. Considering seasonal water temperatures, simple Arrhenius kinetics arguments based on these predicted binding energies help explain the environmental observations of MCHM leeching from the GAC filters months after the spill. Overall, this work shows the important implications that can arise from detailed interfacial chemistry investigations.

Graphical Abstract

1. Introduction

Many aspects of interfacial chemistry are subtle and detailed [1,2,3,4], but they can have important implications for human health and environmental conservation [4,5,6], ranging from the vital transport of breathing gases across the lung surfactant layer [4] to the potentially deleterious effects of gas reactions on the surface of smoggy aerosols [5,6,7,8,9] and dusty grains [10,11] in the atmosphere. We likely associate interfacial chemistry with being particularly important in heterogenous catalysis processes, which have been the center of industrial advancement over many decades [12,13]. Herein, we investigate a seemingly small facet of an interfacial problem that was at the nexus of concerns regarding industrial malfeasance, human health, and environmental impacts following what a former US National Science Foundation division program director characterized as “one of the largest human-made environmental disasters in this century” [14].
On 9 January 2014, residents in WV’s Kanawha Valley awoke anticipating a reprieve from a bitter cold polar vortex that had plunged temperatures from an unseasonably warm 10 to −20 °C within the same day [15,16]. After about 48 h of the chill, as fast-frozen pipes began to thaw, some residents fearing their plumbing would break allowed their faucets to drip and leak what would soon be contaminated water. Little did they know that the frost jack effect was also working its destructive forces on the already weakened floor of a chemical storage tank—number 396—situated on the banks of the Elk River [17], just 1.5 miles upstream from the single drinking water intake [18] that serviced all of their homes and businesses (see Figure 1a). Tank 396 contained a mixture of compounds used as coal frothing agents—surfactants designed to foam and hold onto small coal particles at the surface of coal floatation vessels [19,20]; this allows dirt and rock to easily separate from the high-value coal fines. Two holes in the bottom of Tank 396 allowed passage of over 10,000 gallons of the coal washing mixture to escape within ~8 h—mostly during the early morning hours of 9 January [17,21]. Residents smelled and complained of a heavy licorice odor during their breakfast plans and morning commutes [17]. The main component in Tank 396 was 4-methyl cyclohexane methanol, or MCHM (see Figure 1b). MCHM is an oily compound with a higher boiling point (~200 °C) than water [22] and a sweet licorice odor that our human noses can detect at about 100-fold lower concentrations than the best analytical methods that were available at the time [23,24,25,26]. Beyond the ruptured tank bottom, multiple secondary levels of containment failed at the poorly maintained facility, allowing a nearly unimpeded flow into the Elk River [17].
Because of the warning smell, personnel from the water treatment facility were able to assess the spill in time to make quick decisions about the contaminated water drifting their way [17]. Unfortunately, they were hindered in their decision-making process by the lack of data available regarding the physicochemical properties, toxicity, and mitigation strategies for the MCHM mixture [17,27,28]. While the MCHM mixture is used regionally in quite high quantities, it was still classified as a low volume production chemical, and, like thousands of chemicals in commerce today, MCHM is grandfathered into the Toxic Substances Control Act of 1976 [28,29,30,31,32]. Based on limited guidance form the chemical facility, the Kanawha Valley Treatment Plant’s water intake was never closed in fear of depressurizing the system and losing firefighting capabilities [17]. Contaminated water entered the facility, where operators believed turning up the potassium permanganate and using copious granulated activated carbon (GAC) would sufficiently remove the spilled components from the drinking water [17]. These mitigation efforts were insufficient, and contaminated water was distributed across the 2200 miles of plumbing in the single distribution network that serves over 300,000 people [27,28]. Within a few hours, and over the next few days, residents were turning up in local emergency rooms with exposure-related symptoms [28]. A “Do Not Use” order was enacted by the water company and the governor [33,34]. Residents and businesses across nine counties shut off their taps and waited for a solution [17,28].
The treatment strategies were largely ineffective [17]. Sampling data from 5:00 p.m. on 9 January showed an MCHM concentration at the plant intake of 13.7 ppm, which the treatment facility was only able to decrease to 4.6 ppm in samples of the finished water [17]. A day later, samples showed that MCHM levels of 1.04–3.35 ppm at the intake were only being reduced to 1.02–2.4 ppm [17]. Eventually, sample levels were able to consistently maintain <2 ppb levels of MCHM in the drinking water systems, however, issues regularly manifested out in the distribution system. School children returning from long weekends would smell licorice from the taps, weeks after flushing guidance, prompting school cancellations [35,36,37]. One reason for these pop-up incidents is that MCHM had previously seeped into the polymer plumbing during the long contact time while the “Do Not Use” order was in effect. Once the lines were flushed and filled with clean water, the now-reversed concentration gradient across the polymer surface made the MCHM leech back into the water. This was particularly bad in schools where personnel and students would be absent for a long weekend, allowing the stagnant water to build up these MCHM concentrations for many days [38]. As improved spring weather arrived, and surface temperatures increased, however, it was announced on 24 March that MCHM was again being detected (in samples collected on 21 and 22 January) at the finished water output of the treatment facility, but not at the intake [39]. It seemed the GAC was releasing adsorbed MCHM into the water system, necessitating the full-scale removal and replacement of all GAC filters, at a large cost [40].
A simple question motivates this study: Why did the GAC not work better to adsorb the MCHM and keep it out of people’s premise plumbing? The answers may lie in the varying behavior of MCHM cis and trans isomers and, importantly, how these isomers participate in interfacial interactions with environmentally relevant interfaces. In this study, we examine the isomeric dependence on the physisorption potential between MCHM and carbonaceous materials using computational molecular modelling approaches. The physisorption potential is predicted using computation at the molecular mechanic force field, semi-empirical, and density functional levels of theory. We consider MCHM in contact with both planar and amorphous surfaces, with surface-site sampling performed using molecular dynamics techniques. We then use the predicted differences in surface binding in kinetic analysis to determine whether seasonal surface water temperatures explain the observed MCHM concentration trends at the treatment facility.

2. Results

2.1. Computed MCHM Physisorption Potentials with Carbon Surfaces

2.1.1. cis- and trans-MCHM Plus Ovalene PAH Structures

First, we investigated how cis- and trans-MCHM physisorb on a model of an ideal flat carbon surface, such as graphene. We employed the polycyclic aromatic hydrocarbon (PAH) ovalene (C32H16) for this model surface. cis- and trans-MCHM were allowed to optimize on the ovalene surface to form MCHM:PAH complexes, and a predicted physisorption potential was computed using the MMFF94x force field. Conformational averaged physisorption potentials of 6.2 and 6.1 kcal mol−1 were predicted for the cis-MCHM:PAH and trans-MCHM:PAH complexes, respectively. The energies obtained from these force field computations are not generally considered to be quantitatively accurate, but comparisons can be made between the energies with more confidence. In this case, with the difference between the two complexes of only ~0.1 kcal mol−1, the molecular mechanics force field would indicate that on average, the cis- and trans-MCHM show no difference in physisorption potential on our planar carbon model.
The MCHM:PAH isomer complexes were evaluated at higher levels of theory to confirm whether the MMFF94x predictions would hold. The force field-predicted geometries were optimized with the semi-empirical PM7 Hamiltonian. In addition to PM7, energies were also calculated using the popular density functional theory (DFT) method B3LYP with a modest basis set and inclusion of empirical dispersion (labelled B3LYP-D3 in Table 1) (these methods are further detailed in Section 4.1.1). The intermolecular attraction between MCHM and ovalene is predicted to be over twice as strong with PM7 and DFT (13.9 and 13.5 kcal mol−1, respectively, for the cis-MCHM:PAH complex) than with MMFF94x (6.2 kcal mol−1). Comparing the relative cis- and trans-MCHM:PAH results indicate that there is in fact a small but significant difference between the physisorption potentials of the cis- and trans-MCHM isomers with ovalene. The PM7 results indicate that the trans isomer should physisorb about 1.9 kcal mol−1 stronger than the cis. DFT calculations at the PM7 geometries indicate a similar but slightly smaller value, all of which are summarized in Table 1. It seems that some aspects of the attractive potential interaction vary between the cis- and trans-MCHM isomers interacting with the planar surface, and these subtle effects are not captured by the simple molecular mechanics model.

2.1.2. cis- and trans-MCHM Plus AC Structures

While the MCHM isomers show no preferential absorbance to the planar surface, GAC and other carbonaceous filter media are not perfect planar structures. We next examine how the cis- and trans-MCHM physisorb on a more realistic model of these amorphous surfaces. We used as a starting structure, an amorphous carbon model specifically designed with realism in mind, rather than a random model [41]. After equilibration and sampling of this AC slab model, MCHM was allowed to optimize on the surface to form an adsorbed MCHM:AC species. Energy optimization was performed using two different force fields, MMFF94x, and Amber14:EHT within MOE2019. The discrepancy in predicting quantitative energies with MD is apparent in our results. The particular parameters within Amber14:EHT result in predicted physisorption potentials ~4 kcal mol−1 stronger than those predicted by MMFF94x. As with the PAH surface, though, we are interested and more confident in the relative energy differences. Both force fields predict similar results: the trans-MCHM binds to the AC surface about 2 kcal mol−1 stronger. trans-MCHM:AC is predicted to be 2.4 (Amber14:EHT) and 1.3 (MMFF94x) kcal mol−1 greater than the cis-MCHM:AC on average. The true difference in physisorption potential may be even greater, as inspection of individual values showed the Amber14:EHT potentials consistently achieving 3–4 kcal mol−1 stronger binding with the trans isomer than the cis. To further investigate this possibility, the PM7 semi-empirical method was used to calculate the physisorption potential at the Amber14:EHT geometries. The PM7 method predicts even stronger binding to the AC (>10 kcal mol−1 stronger than MMFF94x predictions), with the relative energy difference indicating that the trans species adsorbs stronger that the cis species by 4.7 kcal mol−1. Overall, these data imply that the true difference in binding between the two isomer species on amorphous carbon surfaces is likely to be 4 kcal mol−1 or greater.

2.2. Kinetic Analysis of MCHM Isomer Desorption from Carbon Surfaces

With the prediction that the trans-MCHM binds ~4 kcal mol−1 stronger (or even more) to AC than cis-MCHM, we can examine how this would manifest in the relative desorption kinetics of these two species. Our best estimate for water temperatures is 0–2 °C and ~10 °C for 19 January and 22 March, respectively; we therefore use 273 and 283 K for the relevant kinetic temperatures. Using the Arrhenius equation, MCHM isomers with two different binding energies will desorb at different rates at a particular temperature given by the ratio,
r c i s r t r a n s = e Δ ε k b T ,
where Δ ε is the difference in physisorption potentials (the MCHM-AC binding energy) between the two isomers. This ratio varies with temperature, and as temperature increases the desorption rate of the more tightly bound species “catches” up with the desorption rate of the more loosely bound species. How quickly the rates catch up with each other depends on the magnitude of Δ ε . In Table 2, we include the results for Δ ε values between 1.3 and 8 kcal mol−1 and from 0 to 30 °C (273–303 K).

3. Discussion

3.1. Isomeric Differences in MCHM Properties

In the wake of the MCHM spill, much work was carried out to characterize the properties and behavior of MCHM, and to a lesser extent, the other components of the crude MCHM mixture in Tank 396. Various computational models are regularly used in incident response situations as they can predict—almost instantaneously—physicochemical, fate and transport, and biological activity properties [28]. However, it was quickly pointed out by researchers that most quick field-level estimations did not account for differences in isomers [25,27,29,31,42,43]. While there was some variation between groups due to methods and sampling, it was determined that MCHM made up approximately 85% of the spilled liquid, with a cis:trans ratio of about 1:1.75 [43]. Of particular relevance to examine in this study are experimental and computed estimates of partitioning coefficients, aqueous solubility, and carbon sorption behavior. The octanol–water partitioning coefficient, log KOW, can be estimated using these programs. EPA’s EPISuite [44] returns a log KOW value of 2.55, and ARChem’s SPARC [45] predicts the same value to be 2.43 [46]. Other programs may report these log KOW values as the more specific log P or log D, which vary depending on whether the solute ionizes in solution. As MCHM does not ionize in water to any appreciable extent unless the pH is very low (MCHM’s pKa <2), log P and log D should not significantly differ at neutral pH. At a pH of 7.4, ACD/Labs®’ Percepta® [47] and ChemAxon’s Chemicalize [48] platforms predict log P (log D) values 2.36 and 1.88 (2.24 and 1.88), respectively. None of these models consider the isomeric nature of the MCHM molecule. Experimentally, log KOW was determined using reverse-phase HPLC to be 2.35 and 2.46 for the cis- and trans-MCHM isomers, respectively [43].
While these estimates are not horrible in comparison with the experiment, larger discrepancies are seen in estimates of water solubility. EPISuite predicts a water solubility value of 2024 mg L−1, while Chemicalize predicts 969 mg L−1. Within SPARC and COSMO-RS [49], a particular temperature can be specified. SPARC predicts solubility values of 2900 and 1340 mg L−1 at 21 and 5 °C, respectively, while COSMO-RS gives 2575 and 2660 mg L−1 at the same temperatures [46]. The two methods predict contrary temperature trends. Experimentally, aqueous solubilities at 23 °C were determined using the slow-stir method to 2600 and 2020 mg L−1 for the cis- and trans-MCHM, respectively [43]. At 4 °C, the values are 2780 and 2240 mg L−1 [43]. These aqueous solubilities reported in Ref. [43] are based on the experimentally determined value for a 50:50 mole fraction mixture, and as such are one-half of the value expected for pure compounds of cis- or trans-MCHM. To facilitate direct comparison with the other predictions in the text, we have extrapolated this result directly, so there may be small discrepancies to the true solubility of the pure isomeric compounds. The experimental results clearly show large solubility differences between the isomers, which are not able to be predicted from the algorithmic models. Much of the difference between solubility and partitioning behavior can be explained by considering the specific dipole character of the individual isomers. Of two similar sized and charged compounds, increased aqueous solubility and increased partitioning out of the octanol phase and into the water phase should be seen for the compound with the greater dipole moment. A detailed quantum chemical study revealed that to fully appreciate this difference in the MCHM isomers, both Boltzmann averaging of the low-lying conformers, as well as solvation effects had to be treated [42]. Computed dipole moments in the gas phase predicted the cis isomer to be more polar; when solvated, the trans species is predicted to be more polar, which is consistent with the experimental results [43].

3.2. MCHM Sorption and Desorption

Experimental groups have also examined MCHM sorption and desorption in the lab on polymers [38], coal [19,20,50], activated carbons [31,50], and riverine sediment [31]. Studies of polymer plumbing were mentioned earlier in context of the school closings. They showed that MCHM readily sorbed and desorbed from polyethylene and epoxy pipe linings, and even cause polyurethane to swell and deform [38]. Taken together, the coal and activated carbon experimental results indicate that while GAC readily sorbs MCHM, it also allows significant desorption [31,38,50]; on the other hand, while coal does not sorb as much MCHM, the fraction that does sorb is held tightly and does not readily desorb [19,20,31,50]. Only one of these studies examined the specific sorption differences in light of the specific isomers [31]. Weidhaas and co-workers observed an increase in the MCHM cis/trans ratio in solution after allowing sorption onto powdered activated carbon [31]. An even larger increase in the ratio was seen in laboratory work on adsorption to sediments from the Elk River [31]; major components of these sediments are carbonate- and silica-based silts and clays, inferring that the isomeric differences in sorption and desorption could be even greater in other filter media, e.g., sand filters.
Our calculations support Weidhaas’s observations. Because the trans-MCHM species adsorbs to AC with a substantially larger physisorption potential (by ~4 kcal mol−1 or more) than the cis, the trans will preferentially adsorb on the surface of amorphous carbons. The molecular mechanics predictions indicate that the MCHM binds with a strength comparable to a strong hydrogen bond, with predictions of about 7.2 and 11.1 kcal mol−1 for the MMFF94x and Amber14:EHT force fields, respectively; this is not entirely surprising based on the presence of an alcohol group. What is more surprising is that the trans species binds with a significantly greater higher physisorption potential (8.5 and 13.4 kcal mol−1, for the same force fields) than the cis. Even though it is on the opposite end of the molecule and separated from the hydroxyl group with a nearly-free rotatable methylene unit, the specific orientation of the cis-methyl group must still sterically interfere with the MCHM’s optimal orientation in a way that it does not when it is in the trans position. In ongoing work in our lab, we are pursuing higher level theory to confirm these energetics using a cluster model of the AC slab. In Figure 2b, the green mesh indicates a 12 Å interaction sphere from the MCHM center of mass that is scooped out from the surface and capped with hydrogens. This smaller size model (circa 600 atoms vs. 6000 in the slab) can be approached with higher levels theory. These results will be reported in different context in the future, but here we can report results of semi-empirical calculations using the PM7 Hamiltonian in Gaussian16 [51]. Using the cluster model, single-point PM7 energies computed at the Amber14:EHT geometries predict that the binding of MCHM may actually be much stronger than the molecular mechanics models predict, in the range of 15–23 kcal mol−1 (average of 18.6 kcal mol−1) for the cis-MCHM species and a range of 20–28 kcal mol−1 (average of 23.2 kcal mol−1) for the trans-MCHM species. Additional preliminary PM7 energy computations at the PM7-optimized cluster geometry indicate that MCHM would physisorb even more strongly (~25% deeper, e.g., 19–33 kcal mol−1 rather than 15–28 kcal mol−1). Next, we will discuss how these energy differences manifest in predicted desorption kinetic rates.

3.3. Kinetic Analysis of MCHM Isomer Desorption from Carbon Surfaces

We have three pieces of information that can be used to infer the relative rate of desorption for different isomers of MCHM: (1) the MCHM-AC physisorption potential, E a ; (2) differences in physisorption potential between the cis- and trans-MCHM, Δ ε ; and (3) a change in seasonal water temperatures, ∆T. There remain uncertainties about the exact quantitative values of the physisorption potentials, but we now have good estimations of the reasonable ranges. Table 2 displays the relative rates of cis:trans desorption kinetics at specific temperatures as a function of the difference between their physisorption potentials. The actual desorption rates are also proportional to the actual physisorption potential (as opposed to the difference between the isomers). In this way, the physisorption potential acts as the activation energy barrier for desorption, and the rate is proportional to the population of sorbed molecules that have a larger kinetic energy than the physisorption potential. This high-energy fraction,   frac h i ε , of the population is obtained from integrating the Maxwell–Boltzmann distribution of energies,
frac h i ε = E a 2 π ( π k b T ) 3 2 ε 1 2 e ε k b T d ε ,
where E a here is the physisorption potential and the integral is taken over the possible kinetic energies, ε . At a particular temperature, the fraction of molecules in the distribution with sufficient energy to overcome the physisorption barrier decreases as E a increases. Conversely, when a particular E a is set, the fraction of sufficiently high-energy molecules increases as T is increased. We can combine these fractions with the relative rate data in Table 2 to arrive at an enhancement factor value, EF, that indicates how the total desorption rate is enhanced when the temperature is increased, and as a function of the activation energy to desorption. We plot in Figure 3 the EF values for temperature changes of 10 and 30 K over a range of hypothetical cis-MCHM physisorption potential values for Δ ε = 2.4 kcal mol−1. This energy difference corresponds to that predicted by the Amber14:EHT force field; the trend is similar for any value of Δ ε . All these values are reported relative to the cis-MCHM desorption rate at 0 °C. More details of this derivation are in the Appendix A.
Inspection of Figure 3 reveals that an increase of only 10 degrees in the water temperature can have a substantial impact on the desorption rates. Taking the MMFF94x predicted cis-MCHM-AC binding energy of 7.2 kcal mol−1, the total desorption rate at 10 °C is 59% higher than at 0 °C. This EF continues to increase as the physisorption potential increases, such that with an E a value of 16 kcal mol−1, there is approximately 3 times the desorption rate at 10 °C as would be seen at 0 °C. Furthermore, also plotted are comparable EF values for a 30 °C temperature, indicative of peak summer water temperatures on the Elk River. The EF is quite impressive, with expected ~3.7-, 7.6-, and 20-fold desorption rate increases compared to that at 0 °C, based on the 7.2, 11, and 16 kcal mol−1  E a values in the ranges predicted by the MMFF94x, Amber14:EHT, and PM7 computations, respectively.
Digging slightly deeper into the EF values in Figure 3, let us consider the maximal plotted values at each temperature and think about what these values really mean. (For simplicity, we assume a 1:1 cis:trans ratio throughout this analysis.) For the ∆T = +10K EF of ~3 (with a representative value for E a = 16 kcal mol−1), this means for every 1000 MCHM molecules that would desorb at 0 °C, 3000 will desorb at 10 °C in the same timeframe, but the total MCHM is not the whole story. Of those 1000 MCHM molecules desorbing at 0 °C, 991 would be the cis isomer and only 9 would be the trans isomer. Once we move to 10 °C, three times as many molecules would desorb, but the individual isomer ratios change. The 3000 desorbing MCHM molecules consist of 2967 cis and 33 trans species. The “isomer specific” EF would then be 2.99 for cis-, but 3.65 for trans-MCHM. The same analysis gives 20,000 MCHM molecules at 30 °C relative to the 1000 at 0 °C. Of these, 19,706 are cis and 294 are trans, yielding isomer specific EF values of 19.9 for the cis, but 32.7 for the trans.

3.4. Odor Implications in Drinking Water

Thinking about the isomer-specific EF values is especially important for the case of MCHM within the overall context of this investigation. The MCHM spill contaminated residential drinking water supplies with what was a highly odorous compound. In odor threshold testing, it was determined that not only is there a slightly different perceived smell between the cis and trans isomers, but also different odor detection limits. “Licorice odor” was always the headline and the warning sign during and after the MCHM spill. However, it is only the trans-MCHM species that is strongly associated with the infamous telltale sweet licorice scent [25]. The cis-MCHM has more subtle scent characteristics, with perceptions ranging from fruit to mint to nail polish [25]. It is also the trans species which the human nose is uniquely sensitive to. The odor threshold concentrations (OTC) for the trans-MCHM isomer was determined to be 0.060 ppb-v in air, while the cis-MCHM threshold concentration was 120 ppb-v in air [25], meaning the human nose is expected to be ~2000 times more sensitive to the trans-MCHM isomer. An even more extreme difference was seen based on participant age, with data suggesting that subjects <50 years old (min age 20) were over 100,000 times more sensitive to detection of the trans species compared to the cis [25]. From our kinetic analysis, we know that the trans species is the minority component of MCHM desorption from carbonaceous sorbent material, but with it being perceivable at a 2000× lower concentration, the increase in odor as a function of temperature should be striking. Much smaller changes in the trans-MCHM concentrations could be enough to be detected through smell than with the cis-MCHM concentrations. Taste and odor problems in water supplies can drive fear and engender outrage, prompting the public to demand solutions to what is thought to be a risky situation [52]. Even if the contaminant concentration is safe from a health perspective, odor and taste perception are important consumer issues for drinking water operations. Our analysis herein shows that partitioning behavior at the carbon filter interface may have played a significant role in the on-going drinking water quality issues following the MCHM spill in 2014.

4. Materials and Methods

4.1. Molecular Modelling

This section provides more details of the models used in our molecular modelling approach. Structures mentioned in the proceeding sections are included in Figure 4 for visual reference.

4.1.1. cis- and trans-MCHM Plus Ovalene PAH Structures

Ovalene is the smallest symmetrical planar polycyclic aromatic hydrocarbon (PAH) that is large enough in both planar dimensions to allow an MCHM molecule to “lay flat” on the center and still have about one benzene unit in either direction before the edge. Starting structures for cis- and trans-MCHM interacting with the ovalene (C32H14) were arranged by hand to ensure relatively close contact with the center of the ovalene molecule (see Figure 2a), minimizing edge effects. Previous work showed that both MCHM isomers can adopt a plethora of low-energy local minima [42], so it is important to sample the various conformations of MCHM as well as the relative orientations between the MCHM and the ovalene “surface.” Both aspects were accounted for using the LowModeMD search program in the MOE software package [53]. The LowModeMD search algorithm in MOE generates conformations using short (~1 ps) molecular dynamics (MD) simulations at constant temperature, followed by all-atom geometry optimizations [54]. By default, the energy minimization technique in MOE utilizes the MMFF94x force field [55] and incorporates a distance-dependent solvation effect with a dielectric constant of 80 for the “exterior” or solvent-exposed parts of the system. MMFF94x is a general all-atom force field parameterized for small organic molecules, which incorporates improvements over the initial MMFF94 (Merck molecular force field) (for instance, MMFF94x predicts planar conjugated nitrogens, while MMFF94 predicts tetrahedral nitrogen conjugation) [55]. MCHM:PAH complexes with relative energies higher than 7.0 kcal mol−1 from the lowest energy conformer were discarded. These searches provided 120 conformers for the cis-MCHM:PAH complex from 800 poses, and 138 conformers for the trans-MCHM:PAH from 809 poses. Intermolecular energies were determined by taking the MMFF94x energy computed for the MCHM:PAH complexes and subtracting that computed for a lone ovalene. Then, the average of the five lowest MCHM conformers (reported previously, see [42,43]) was subtracted to obtain an apparent physisorption potential that is relative to the first few low-lying conformers.
The MCHM:PAH isomer complexes were also evaluated at higher levels of theory. The MMFF94x force field-predicted geometries were further optimized with the semi-empirical PM7 Hamiltonian in Gaussian16. Of the initial 120 cis- and 138 trans-MCHM:PAH complexes, 114 and 129 structures, respectively, converged to unique geometries. In addition to PM7, energies were calculated using the popular density functional theory (DFT) method B3LYP with a modest basis set, 6–31G(d), as implemented in Gaussian16. Empirical dispersion was included using the D3 version of Grimme’s dispersion method with Becke–Johnson damping [56]. The DFT energies were calculated at the PM7-optimized geometries. For both PM7 and DFT computations, the physisorption potential was taken as the difference in the computed heats of formation between the MCHM:PAH complex and the sum of the lone MCHM molecule and the lone PAH, calculated with the same geometry as was in the optimized complex.

4.1.2. Amorphous Carbon (AC) Model

We used a previously thermally equilibrated periodic slab model of amorphous carbon as a starting structure [41]. The initial structure was originally developed using a semi-empirical tight binding approach implemented in a canonical Monte Carlo algorithm. An advantage of the model is that it gives a particularly good description of the sp, sp2, and sp3 hybrid bonds; the resulting structure was also shown to be more accurate to describe Cu, Ag, and Au atom adsorption than a model constructed from randomly generated carbon positions [41]. To prepare a slab model of an amorphous carbon surface, the z-dimension of the 5 × 5 × 3 nm periodic box was extended to 25 nm, while the x- and y-dimensions were maintained with periodic boundary conditions, creating two solid-vacuum interfaces. Dangling bonds were saturated with hydrogens in MOE, and the surface was allowed to equilibrate through a molecular dynamics run in MOE 2013. Molecular dynamics were carried out with NAMD 2.9 [57] at 300 K. An initial 100 ps equilibration run was followed by a 500 ps production run. Molecular geometry and velocity snapshots of the AC slab taken along the production run were saved for later use. The resulting AC surface was then used to study adsorption of MCHM as a surrogate model for GAC and similar materials. It is noted, however, that our AC model is of pure carbon capped with hydrogen. Real GAC and other materials will have some level of oxidized functional groups on the surface as well, and this simple model does not capture effects from those groups.

4.1.3. cis- and trans-MCHM Plus AC Structures

AC Slab Model. Adsorption sites on the central ~3 × 3 nm area of one face of the AC were sampled for each isomer by roughly aligning a subset of the MCHM:PAH conformations with the AC surface. The PAH was aligned with the AC plane, and the PAH removed, leaving the MCHM in a favorable starting orientation. Randomized orientations were selected, and the three lowest MCHM:PAH conformers for each isomer were intentionally included in the sample set. The MCHM isomer was then allowed to minimize to the rigid AC surface using one of two force fields, MMFF94x or Amber14:EHT [58,59], as implemented in MOE2019 [53]. Predicted physisorption potentials for each MCHM:AC were taken as the intermolecular van der Waals interaction energy computed between the MCHM and the AC.
AC Cluster Model. Subsequently, the PM7 semi-empirical method was used to calculate the physisorption potential at the Amber14:EHT geometries using a cluster model of the surface. To form the cluster model, all atoms within a 12 Å interaction sphere as measured from the MCHM center of mass were selected (see Figure 2b). Exterior atoms were discarded and dangling bonds were automatically capped with hydrogen within MOE2019; structures that resulted in excess charge or non-singlet multiplicities were discarded. About two-thirds of the Amber14:EHT MCHM:AC geometries remained and were used as input for PM7 computations in Gaussian16. The physisorption potential was taken as the difference in the computed heats of formation between the cluster complex and the sum of the lone MCHM molecule and the lone cluster, calculated at the same geometry as was in the cluster complex. Because several structures were discarded automatically, it is especially important to ensure that the differences between the cis (M = 18.6, SD = 2.4) and trans (M = 23.2, SD = 2.3) isomers are significant and not biased by this sampling. Statistical analysis indicates the differences in the means is significant with the independent two-tailed t-test yielding p < .001, with a size effect of 2.006 (Cohen’s d value).

4.2. Surface Water Temperature Estimations

4.2.1. Historical Average Water Temperatures

Ideally, water temperature data from the water treatment plant, and specifically of the GAC filters, would be used in the kinetic analysis herein. Those data were not readily available. Water temperatures after the spill were reported as 3.5 ± 0.1 °C and 3.2 ± 0.7 °C on 16 January and 20 January, and 24.6 °C on 18 August [31], a few days and then months after the relevant days. However, historical daily water temperature measurements were found to be available from a USGS gauging station (USGS 03198000) on the Kanawha River in South Charleston, WV, (see Figure 1) only ~4.4 river miles downstream from the water facility intake [60]. These data span the period from 1 April 1950—9 April 1985; other than a gap in coverage from 7 August 1970—1 October 1971, 93% of days have recorded minimum and maximum daily water temperatures. Average historical daily water temperatures were calculated from measurements taken within one week of the two significant dates in the MCHM spill story—9 January and 22 March. Average historical water temperature highs and lows for January 2–16 were 5.58 (407 temps, σ = 2.45, range 0.5–15.0 °C) and 4.21 °C (407 temps, σ = 2.29, range 0.5–10.6 °C), respectively. For March 15–29, average highs and lows were 8.96 °C (420 temps, σ = 1.76, range 4.0–14.5 °C) and 8.46 °C (420 temps, σ = 1.69, range 3.9–14.0 °C), respectively. This shows that historically there should be an expected ~4 °C difference in river water temperatures in the area from 9 January to 22 March. This historic average does not however account for incidental variation in water temperatures due to local weather and microclimate events.

4.2.2. Accounting for Effects of Ambient Air Temperature (i.e., Weather)

Rather than rely on historical water temperature data, which at best was 23 years old relative to the spill incident, it would be convenient to use ambient atmospheric weather data (i.e., regular weather reports), as these data are available hourly from the nearby Yeager Airport located only 0.6 and 2.2 miles from the spill site and water intake, respectively (see Figure 1) [61]. Attempts to correlate historical water data with the weather reports to produce a simple predictive model were not fruitful. It should perhaps be obvious that the high heat capacity of water largely resists the sometimes-quick changes in weather. This is easily seen in the historical water temperature data, which show average daily temperature variations of only 1.08 (2–16 January, σ = 1.14), 0.49 (15–29 March, σ = 0.46), and 1.05 (year round, σ = 0.91). For a specific example of how contrary these measures can be, on 9 January 1966, reported temperatures were 7.8/7.8 °C (hi/lo, water) and 4.4/−3.9/−12.2 °C (hi/avg/lo, air), while on 9 January 1984 those same measures were 2.0/2.0 °C (hi/lo, water) and 15/4.4/−3.3 °C (hi/avg/lo, air). It should also be noted that precipitation events likely shift water temperatures more quickly than air temperature changes, but we make no attempt to correlate/correct for the effects in this study.
While weather may not change water temperatures as fast as air temperatures, our predictions based on historical averages can be hedged up or down based on the unique and specific weather conditions within the days leading up to the spill incident (~9 January) and re-appearance of MCHM at the plant (~21 March). Historical weather averages report air temps of ~1.1 °C for Jan, in line with the 9 January 2014 average temperature of 0.97 °C; however, leading up to this date, temperatures were well below average, with average daily temperatures of −6.7, −16.0, and −8.0, for 8, 7, and 6 January, respectively, with a new all-time record low recorded on 7 January of −19.4 °C. These sustained lows would certainly drive down the water temperature, leading us to an estimate of 0–2 °C used in the kinetic analysis. Historical weather averages report air temps of ~8.9 °C for March, and the 22 March 2014 average temperature was close at 8.6 °C. However, while 22 March was at the average temperatures, the days previous were unseasonably hot. Daily average temperatures were reported as 12.7, 14.5, and 3.5 °C, for 21, 20, and 19 March, respectively, with highs of 20.6 and 20 °C on 21 and 19 March, respectively. The sustained higher temperatures would increase the water temperature, lending confidence to estimating water temps above 10 °C used in the kinetic analysis. Finally, we note that the maximum water temperature recorded in the historical data was 35 °C, while temperatures at the gauging station recorded in mid–late summer consistently reached a degree or two above 30 °C. Coupled with the reported ~25 °C water temp for 18 August [31], this gives validation to perform a kinetic analysis at these higher temperatures.

5. Environmental Implications and Concluding Remarks

Overall, our findings indicate that MCHM may physisorb to planar carbon surfaces with attractive potentials of around 15 kcal mol−1 and may bind even stronger to amorphous/rough carbon surfaces (circa 20–25 kcal mol−1). Importantly, however, physisorption is not equal with both isomers. The trans-MCHM isomer is predicted to bind to amorphous carbon surfaces ~4 kcal mol−1 or more strongly than the cis-MCHM isomer. The preferential binding effect is more subtle on planar carbon, such that higher levels of theory than molecular mechanics force fields are needed to sufficiently predict this difference. A kinematic analysis of relative desorption rates taking these differential adsorption strengths of the MCHM isomers into account indicates that an overall enhancement factor of trans desorption relative to cis could be 20-fold or higher based on a shift from winter-to-summer seasonal water temperatures.
With regard to the aftermath of the 14 January MCHM spill, our data and analysis help to explain the contrary results of MCHM detection in the finished water months after it was undetectable in the source water. Of course, our model is simple and uses ideal, clean surfaces in our computations. Our model does not account for complicating effects found in the field, such as dissolved biomatter, competition with other would-be adsorbates, etc., but the basic concepts still hold. Nevertheless, even this simplistic level of analysis indicates that (1) binding of MCHM to GAC and other carbonaceous filter media is not strong enough for permanent removal from the water system, (2) differences in the isomer physisorption potential would lead to differential sorption and desorption behaviors, and (3) desorption from filter media into the water system should be expected as seasonal temperatures increase. Such insights could have influenced decision making to promulgate mitigation and removal of the contaminated GAC filter media prior to spring temperatures allowing further release of the compound into the drinking water supply. Our data indicate that the Elk River water temperature changes from Jan (~0 °C) to Mar (~10 °C) could have resulted in increasing the total rate of MCHM desorption by as much as a factor of three; had the GAC filters been allowed to continue to desorb MCHM into the heat of the summer (where water temperatures routinely exceed 30 °C), MCHM desorption rates could have increased to over 20 times the initial desorption rates in Jan. Because of the different isomeric desorption rates, as water temperature increases, odor perception of the increased MCHM desorption could see even more drastic scaling than the desorption rates indicate. Beyond the Elk River incident, ongoing development of computational libraries of amorphous surfaces in the author’s lab can be combined with higher level QM/MM calculations to refine the predictions herein, and to make similar predictions for other substances. The general methodology detailed in this study can be applied for future chemicals whose under-characterized or orphan status necessitate estimates of surface binding to many types of amorphous interfaces relevant for filter media, environmental interfaces, premise plumbing, etc.

Funding

This research was supported by the National Science Foundation through grant numbers CBET-1435289 and CHE-1905207. The APC was funded by CHE-1905207.

Data Availability Statement

Data presented in this study may be made available on request from the corresponding author.

Acknowledgments

The author is thankful to Andrew Whelton, Nathan DeYonker, and Katherine A. Charbonnet for helpful discussions regarding aspects of this work. The High Performance Computing Center and the Computational Research on Materials Institute at The University of Memphis (CROMIUM) also provided generous resources for this research.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Samples of the compounds are not available from the author.

Appendix A

Total MCHM Desorption Rate Enhancement Factor. The amount of cis-MCHM that desorbs is proportional to the fraction of cis molecules, frac h i ε , c i s , that have enough energy to overcome a given physisorption potential, E a , at a given temperature, T; this fraction is obtained by integrating the Maxwell–Boltzmann distribution (Equation (2)) over the E a limits. The fraction of trans-MCHM molecules that desorb relative to frac h i ε , c i s is obtained using the relative rates from the Arrhenius equation (Equation (1)) for a particular temperature, T, and given the difference between the isomers’ physisorption potentials, Δ ε , as
frac h i ε ,   t r a n s = frac h i ε ,   c i s r c i s r t r a n s .
The total amount of MCHM desorption at a given set of T, E a , and Δ ε values is
total   MCHM ( T ,   E a ,   Δ ε ) = amount   of   cis - MCHM + amount   of   trans - MCHM   ,
which is proportional to the computed fractions,
total   MCHM ( T ,   E a ,   Δ ε )   frac h i ε ,   c i s + frac h i ε ,   t r a n s
  frac h i ε ,   c i s + frac h i ε ,   c i s r c i s r t r a n s .
The enhancement factor plotted in Figure 3 is given by the ratio
EF = total   MCHM ( T 2 ,   E a ,   Δ ε ) total   MCHM ( T 1 ,   E a ,   Δ ε )
where all EF values are defined in this paper relative to the total MCHM proportions with T1 = 273K (0 °C) for each set of E a ,   Δ ε values. In other words, we set the denominator equal to 1 for each set of E a ,   Δ ε values at 0 °C. The EF value is then always given in this paper as the enhancement of the total MCHM desorption rate at different temperatures, relative to that at 0 °C. As an example, if we obtain, using Equation (A5), an EF value of 10 at T2, the total rate of MCHM is 10-fold more at T2 than what it would be at 0 °C.

References

  1. Alexander, W.A. Particle Beam Scattering from the Vacuum-Liquid Interface. In Physical Chemistry of Gas-Liquid Interfaces; Faust, J.A., House, J.E., Eds.; Elsevier: Amsterdam, The Netherlands, 2008. [Google Scholar] [CrossRef]
  2. Nathanson, G.M.; Davidovits, P.; Worsnop, D.R.; Kolb, C.E. Dynamics and Kinetics at the Gas−Liquid Interface. J. Phys. Chem. 1996, 100, 13007–13020. [Google Scholar] [CrossRef]
  3. Alexander, W.A.; Zhang, J.; Murray, V.J.; Nathanson, G.M.; Minton, T.K. Kinematics and dynamics of atomic-beam scattering on liquid and self-assembled monolayer surfaces. Faraday Discuss. 2012, 157, 355–374. [Google Scholar] [CrossRef]
  4. Davidovits, P.; Kolb, C.E.; Williams, L.R.; Jayne, J.T.; Worsnop, D.R. Mass Accommodation and Chemical Reactions at Gas−Liquid Interfaces. Chem. Rev. 2006, 106, 1323–1354. [Google Scholar] [CrossRef]
  5. Goldstein, A.H.; Galbally, I.E. Known and Unexplored Organic Constituents in the Earth’s Atmosphere. Environ. Sci. Technol. 2007, 41, 1514–1521. [Google Scholar] [CrossRef] [Green Version]
  6. Mohr, C.; Huffman, J.A.; Cubison, M.J.; Aiken, A.C.; Docherty, K.S.; Kimmel, J.R.; Ulbrich, I.M.; Hannigan, M.; Jimenez, J.L. Characterization of Primary Organic Aerosol Emissions from Meat Cooking, Trash Burning, and Motor Vehicles with High-Resolution Aerosol Mass Spectrometry and Comparison with Ambient and Chamber Observations. Environ. Sci. Technol. 2009, 43, 2443–2449. [Google Scholar] [CrossRef]
  7. Donaldson, D.J.; Vaida, V. The Influence of Organic Films at the Air−Aqueous Boundary on Atmospheric Processes. Chem. Rev. 2006, 106, 1445–1461. [Google Scholar] [CrossRef] [PubMed]
  8. Lu, J.W.; Flores, J.M.; Lavi, A.; Abo-Riziq, A.; Rudich, Y. Changes in the Optical Properties of Benzo [a] pyrene-Coated Aerosols upon Heterogeneous Reactions with NO2 and NO3. Phys. Chem. Chem. Phys. 2011, 13, 6484–6492. [Google Scholar] [CrossRef] [PubMed]
  9. Freidlander, S.K. Smoke, Dust, and Haze: Fundamentals of Aerosol Dynamics, 2nd ed.; Oxford University Press: New York, NY, USA, 2000. [Google Scholar]
  10. Hatch, C.D.; Grassian, V.H. 10th Anniversary Review: Applications of analytical techniques in laboratory studies of the chemical and climatic impacts of mineral dust aerosol in the Earth’s atmosphere. J. Environ. Monit. 2008, 10, 919–934. [Google Scholar] [CrossRef]
  11. Prather, K.A.; Hatch, C.D.; Grassian, V.H. Analysis of atmospheric aerosols. Annu. Rev. Anal. Chem. 2008, 1, 485–514. [Google Scholar] [CrossRef] [PubMed]
  12. Ertl, G.; Knözinger, H.; Weitkamp, J. Handbook of Heterogeneous Catalysis; Wiley-VCH: Weinheim, Germany, 1997. [Google Scholar]
  13. Ertl, G. Reactions at Surfaces: From Atoms to Complexity (Nobel Lecture). Angew. Chem. Int. Ed. 2008, 47, 3524–3535. [Google Scholar] [CrossRef]
  14. NSF News Release 14-018. NSF Awards Rapid Response Grants to Study West Virginia Chemical Spill. National Science Foundtion. [Online]. 30 January 2014. Available online: https://www.nsf.gov/news/news_summ.jsp?cntn_id=130304 (accessed on 30 April 2021).
  15. Charleston WV Weather History (6 January 2014), Yeager Airport Station. Available online: https://www.wunderground.com/history/daily/us/wv/charleston/KCRW/date/2014-1-6 (accessed on 29 April 2021).
  16. West Virginia Department of Health and Human Resources. 2014 Winter Storm Information for West Virginia. 7 January 2014. Available online: https://dhhr.wv.gov/healthprep/eventsandannouncements/Pages/2014-Winter-Storm-Information----WV.aspx (accessed on 4 May 2021).
  17. U.S. Chemical Safety and Hazard Investigation Board. Investigation Report No. 2014-01-I-WV. “Chemical Spill Contaminates Public Water Supply in Charleston, West Virginia”. Available online: https://www.csb.gov/assets/1/20/final_freedom_industries_investigation_report_(5-11-2017).pdf?15829 (accessed on 28 April 2021).
  18. Markham, P.; Gianato, J.; Hoyer, J. “After Action Review Emergency Response To 9 January 2014 Freedom Industries Chemical Leak” Charleston, WV, USA. 2014. Available online: https://emd.wv.gov/West%20Virginia%20Public%20Water%20Supply%20Study%20Commission/Documents/After%20Action%20Review.PDF (accessed on 4 May 2021).
  19. He, Y.T.; Noble, A.A.; Ziemkiewicz, P. Investigation of MCHM transport mechanisms and fate: Implications for coal beneficiation. Chemosphere 2015, 127, 158–163. [Google Scholar] [CrossRef] [PubMed]
  20. Noble, A.; He, Y.T.; Ziemkiewicz, P. Partitioning Behavior of 4-Methyl Cyclohexane Methanol in Two Appalachian Coal Preparation Plants. Int. J. Coal Prep. Util. 2016, 36, 69–81. [Google Scholar] [CrossRef]
  21. West Virginia Attorney General, “Elk River Chemical Spill Incident Report” 8 January 2015. Charleston, WV, USA. Available online: http://ago.wv.gov/Documents/010815-ElkRiverChemicalSpill.PDF (accessed on 4 May 2021).
  22. Luttrell, W.E. 4-Methylcyclohexane methanol (MCHM). J. Chem. Health Saf. 2015, 22, 39–41. [Google Scholar] [CrossRef]
  23. McGuire, M.J. Odor Study Technical Memorandum; WVTAP: Santa Monica, CA, USA, 2014. [Google Scholar]
  24. Ward, K., Jr. Virginia Tech Experts Pinpoint MCHM Odors West Virginia Gazette. [Online]. 26 March 2014. Available online: https://www.wvgazettemail.com/news/special_reports/virginia-tech-experts-pinpoint-mchm-odors/article_f928c3ee-96b6-5ce0-b853-bd93482a837f.html (accessed on 5 May 2021).
  25. Gallagher, D.L.; Phetxumphou, K.; Smiley, E.; Dietrich, A.M. Tale of two isomers: Complexities of human odor perception for cis- and trans-4-methylcyclohexane methanol from the chemical spill in West Virginia. Environ. Sci. Technol. 2015, 49, 1319–1327. [Google Scholar] [CrossRef]
  26. Sain, A.; Dietrich, A.M.; Smiley, E.; Gallagher, D.L. Assessing human exposure and odor detection during showering with crude 4-(methylcyclohexyl) methanol (MCHM) contaminated drinking water. Sci. Total Environ. 2015, 538, 298–305. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Whelton, A.J.; McMillan, L.; Novy, C.L.-R.; White, K.D.; Huang, X. Case study: The crude MCHM chemical spill investigation and recovery in West Virginia USA. Environ. Sci. Water Res. Technol. 2017, 3, 312–332. [Google Scholar] [CrossRef] [Green Version]
  28. Weidhaas, J.L.; Dietrich, A.M.; DeYonker, N.J.; Dupont, R.R.; Foreman, W.T.; Gallagher, D.; Gallagher, J.E.G.; Whelton, A.J.; Alexander, W.A. Enabling Science Support for Better Decision-Making when Responding to Chemical Spills. J. Environ. Qual. 2016, 45, 1490–1500. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Cooper, W.J. Responding to crisis: The West Virginia chemical spill. Environ. Sci. Technol. 2014, 48, 3095. [Google Scholar] [CrossRef]
  30. Schnoor, J.L. Re-emergence of emerging contaminants. Environ. Sci. Technol. 2014, 48, 11019–11020. [Google Scholar] [CrossRef]
  31. Weidhaas, J.; Lin, L.-S.; Buzby, K. A case study for orphaned chemicals: 4-methylcyclohexanemethanol (MCHM) and propylene glycol phenyl ether (PPH) in riverine sediment and water treatment processes. Sci. Total Environ. 2017, 574, 1396–1404. [Google Scholar] [CrossRef]
  32. USEPA. The Use of Structure-Activity Relationships (SAR) in the High Production Volume Chemicals Challenge Program. Available online: https://nepis.epa.gov/Exe/ZyPURL.cgi?Dockey=P1008MA9.TXT (accessed on 4 May 2021).
  33. West Virginia, Office of the Governor, “Gov. Tomblin’s Update on State of Emergency” West Virginia Press. [Online]. 10 January 2014. Available online: https://wvpress.org/copydesk/insight/gov-tomblins-update-state-emergency/ (accessed on 5 May 2021).
  34. West Virginia American Water, American Water Company Customer Advisory ** Alert **. [Online]. 9 January 2014. Available online: https://dhhr.wv.gov/healthprep/eventsandannouncements/Pages/American-Water-Company-Customer-Advisory.aspx (accessed on 5 May 2021).
  35. Associated Press. School Closes after Licorice Odor Present. Herald Dispatch. [Online]. 17 February 2014. Available online: https://www.herald-dispatch.com/news/recent_news/school-closes-after-licorice-odor-present/article_5b79ae09-f1bf-5b27-8045-8ab9db7a7baa.html (accessed on 4 May 2021).
  36. Jenkins, J. Kanawha Closes More Schools over Licorice Odor. West Virginia MetroNews. [Online]. 6 February 2014. Available online: https://wvmetronews.com/2014/02/06/kanawha-closes-more-schools-over-licorice-odor/ (accessed on 4 May 2021).
  37. Mays, M. Water Response Team Called to Four Schools Monday. Charleston Gazette-Mail. [Online]. 17 February 2014. Available online: https://www.wvgazettemail.com/news/special_reports/water-response-team-called-to-four-schools-monday/article_872ca250-83c7-57a2-9a80-02a55f8ed269.html (accessed on 4 May 2021).
  38. Ahart, M.; Gallagher, D.L.; Scardina, P.; Dietrich, A.M. Industrial Spills and Water Distribution: Crude MCHM Sorption and Desorption in Polymer Pipes and Linings. J. Environ. Eng. 2016, 142, 04016045. [Google Scholar] [CrossRef]
  39. Mistich, D. Update: Testing at WVAM Indicates Filtering Process Adds Trace Amounts of MCHM to Water Supply. WVPB. [Online]. 25 March 2014. Available online: https://www.wvpublic.org/news/2014-03-25/update-testing-at-wvam-indicates-filtering-process-adds-trace-amounts-of-mchm-to-water-supply (accessed on 29 April 2021).
  40. Associated Press. WVAW Changing Carbon Filters at Charleston Plant. WVPB. [Online]. 1 April 2014. Available online: https://www.wvpublic.org/news/2014-04-01/wvaw-changing-carbon-filters-at-charleston-plant (accessed on 5 May 2021).
  41. Ricolleau, C.; Le Bouar, Y.; Amara, H.; Landon-Cardinal, O.; Alloyeau, D. Random vs. realistic amorphous carbon models for high resolution microscopy and electron diffraction. J. Appl. Phys. 2013, 114, 213504. [Google Scholar] [CrossRef]
  42. De Yonker, N.J.; Charbonnet, K.A.; Alexander, W.A. Dipole moments of trans- and cis-(4-methylcyclohexyl)methanol (4-MCHM): Obtaining the right conformer for the right reason. Phys. Chem. Chem. Phys. 2016, 18, 17856–17867. [Google Scholar] [CrossRef]
  43. Dietrich, A.M.; Thomas, A.; Zhao, Y.; Smiley, E.; Shanaiah, N.; Ahart, M.; Charbonnet, K.A.; DeYonker, N.J.; Alexander, W.A.; Gallagher, D.L. Partitioning, Aqueous Solubility, and Dipole Moment Data for cis- and trans-(4-Methylcyclohexyl)methanol, Principal Contaminants of the West Virginia Chemical Spill. Environ. Sci. Technol. Lett. 2015, 2, 123–127. [Google Scholar] [CrossRef] [Green Version]
  44. US EPA. Estimation Programs Interface Suite™ for Microsoft® Windows, v 4.11; United States Environmental Protection Agency: Washington, DC, USA, 2017.
  45. SPARC Developed by ARChem, Inc. Danielsville, GA, USA.
  46. Whelton, A.J.; McMillan, L.; Connell, M.; Kelley, K.M.; Gill, J.P.; White, K.D.; Gupta, R.; Dey, R.; Novy, C. Residential Tap Water Contamination Following the Freedom Industries Chemical Spill: Perceptions, Water Quality, and Health Impacts. Environ. Sci. Technol. 2015, 49, 813–823. [Google Scholar] [CrossRef] [Green Version]
  47. ACD/Labs Percepta Platform-Phys Chem Module, Advanced Chemistry Development. Available online: http://www.chemspider.com/Chemical-Structure.105625.html (accessed on 2 May 2021).
  48. ChemAxon Chemicalize Platform. Available online: Chemicalize.com (accessed on 2 May 2021).
  49. COSMOS-RS, SCM, Theoretical Chemistry, Vrije Universiteit, Amsterdam, The Netherlands. Available online: https://www.scm.com/COSMO-RS/ (accessed on 2 May 2021).
  50. Jeter, T.S.; Sarvera, E.A.; McNair, H.M.; Rezaeea, M. 4-MCHM sorption to and desorption from granular activated carbon and raw coal. Chemosphere 2016, 157, 106–165. [Google Scholar] [CrossRef] [Green Version]
  51. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian 16, Revision B.01; Gaussian, Inc.: Wallingford, CT, USA, 2016. [Google Scholar]
  52. Ward, K., Jr. MCHM Odors Show Chemical not Gone, Experts Say. Charleston Gazette-Mail. [Online]. 18 March 2014. Available online: https://www.wvgazettemail.com/news/special_reports/mchm-odors-show-chemical-not-gone-experts-say/article_8f5da982-a634-58b8-8418-872c7c9fbe6d.html (accessed on 4 May 2021).
  53. Chemical Computing Group. Molecular Operating Environment (MOE); Chemical Computing Group Inc.: Montreal, QC, Canada, 2019. [Google Scholar]
  54. Labute, P. Low Mode MD-Implicit Low Mode Velocity Filtering Applied to Conformational Search of Macrocycles and Protein Loops. J. Chem. Inf. Model. 2010, 50, 792–800. [Google Scholar] [CrossRef] [PubMed]
  55. Halgren, T.A. The Merck molecular force field. I. Basis, form, scope, parameterization, and performance of MMFF94. J. Comp. Chem. 1996, 17, 490–641. [Google Scholar] [CrossRef]
  56. Grimme, S.; Ehrlich, S.; Goerigk, L. Effect of the damping function in dispersion corrected density functional theory. J. Comp. Chem. 2011, 32, 1456–1465. [Google Scholar] [CrossRef]
  57. Phillips, J.C.; Hardy, D.J.; Maia, J.D.C.; Stone, J.E.; Ribeiro, J.V.; Bernardi, R.C.; Buch, R.; Giacomo Fiorin, G.; Henin, J.; Jiang, W.; et al. Scalable molecular dynamics on CPU and GPU architectures with NAMD. J. Chem. Phys. 2020, 153, 044130. [Google Scholar] [CrossRef] [PubMed]
  58. Maier, J.A.; Martinez, C.; Kasavajhala, K.; Wickstrom, L.; Hauser, K.E.; Simmerling, C. ff14SB: Improving the accuracy of protein side chain and backbone parameters from ff99SB. J. Chem. Theory Comput. 2015, 11, 3696–3713. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Gerber, P.R.; Müller, K. MAB, a generally applicable molecular force field for structure modelling in medicinal chemistry. J. Comput. Aided Mol. Des. 1995, 9, 251–268. [Google Scholar] [CrossRef] [PubMed]
  60. US Geological Survey. National Water Information System Data Available on the World Wide Web (USGS Water Data for the Nation). Available online: https://waterdata.usgs.gov/nwis/inventory/?site_no=03198000&agency_cd=USGS (accessed on 10 April 2021).
  61. Charleston, WV Weather History. Available online: https://www.wunderground.com/history/daily/us/wv/charleston/KCRW (accessed on 30 April 2021).
Figure 1. (a) Locations in the Charleston, WV area relevant to this study. Base map from Google Maps, inset images from Reference [17]; (b) Isomers of 4-methyl-cyclohexane methanol, the major component of the fluid mixture spilled into the Elk River.
Figure 1. (a) Locations in the Charleston, WV area relevant to this study. Base map from Google Maps, inset images from Reference [17]; (b) Isomers of 4-methyl-cyclohexane methanol, the major component of the fluid mixture spilled into the Elk River.
Molecules 26 03411 g001
Figure 2. Planar and amorphous carbon surface models used in this study. (a) Starting structure of cis-MCHM on ovalene PAH. In this view, it can be seen that the ovalene is just large enough to accommodate the maximum footprint of the MCHM molecule. (b) Minimum structure of cis-MCHM adsorption site on an amorphous carbon surface. The green mesh indicates the 12 Å boundary for building the cluster model. Atoms (colored) within the mesh are subject to higher level computation using the PM7 semi-empirical method (see details in Section 4.1.3 of the text).
Figure 2. Planar and amorphous carbon surface models used in this study. (a) Starting structure of cis-MCHM on ovalene PAH. In this view, it can be seen that the ovalene is just large enough to accommodate the maximum footprint of the MCHM molecule. (b) Minimum structure of cis-MCHM adsorption site on an amorphous carbon surface. The green mesh indicates the 12 Å boundary for building the cluster model. Atoms (colored) within the mesh are subject to higher level computation using the PM7 semi-empirical method (see details in Section 4.1.3 of the text).
Molecules 26 03411 g002
Figure 3. Enhancement in the predicted total desorption rate for MCHM from amorphous carbon as a function of hypothetical cis-MCHM physisorption potentials for temperature increases of 10 and 30 K. Enhancement is relative to the rate at 273 K. Results are plotted for the case where the difference in physisorption potentials for cis- and trans-MCHM is 2.4 kcal mol−1; the trend is similar for any value of Δ ε . The triangles highlight relevant values predicted by the molecular mechanics and semi-empirical predictions detailed in the text.
Figure 3. Enhancement in the predicted total desorption rate for MCHM from amorphous carbon as a function of hypothetical cis-MCHM physisorption potentials for temperature increases of 10 and 30 K. Enhancement is relative to the rate at 273 K. Results are plotted for the case where the difference in physisorption potentials for cis- and trans-MCHM is 2.4 kcal mol−1; the trend is similar for any value of Δ ε . The triangles highlight relevant values predicted by the molecular mechanics and semi-empirical predictions detailed in the text.
Molecules 26 03411 g003
Figure 4. Structures used in this study.
Figure 4. Structures used in this study.
Molecules 26 03411 g004
Table 1. Differences in physisorption potentials between MCHM isomers and carbonaceous surfaces.
Table 1. Differences in physisorption potentials between MCHM isomers and carbonaceous surfaces.
Difference in Isomers, Etrans—Ecis 1
SystemMMFF94xAmber14:EHTPM7B3LYP-D3
MCHM:PAH0.1-1.91.4 2
MCHM:AC1.32.44.7 3-
1 Difference in predicted physisorption potential in kcal mol−1. 2 B3LYP/6-31G(d) energies, with D3-BJ dispersion, at the PM7 optimized geometries. 3 Calculations were carried out on the cluster model with Amber14:EHT geometries (see Figure 2).
Table 2. Predicted relative desorption rates of MCHM isomers from amorphous carbon at different temperatures.
Table 2. Predicted relative desorption rates of MCHM isomers from amorphous carbon at different temperatures.
r c i s / r t r a n s
Δ ε   1   0 °C10 °C20 °C30 °C
1.31.29 × 1011.18 × 1011.08 × 1011.00 × 101
2.41.06 × 1029.02 × 1017.74 × 1016.70 × 101
4.02.40 × 1031.82 × 1031.41 × 1031.11 × 103
8.05.53 × 1063.20 × 1061.92 × 1061.19 × 106
1 Difference in physisorption potential in kcal mol−1.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Alexander, W.A. Methyl-Cyclohexane Methanol (MCHM) Isomer-Dependent Binding on Amorphous Carbon Surfaces. Molecules 2021, 26, 3411. https://doi.org/10.3390/molecules26113411

AMA Style

Alexander WA. Methyl-Cyclohexane Methanol (MCHM) Isomer-Dependent Binding on Amorphous Carbon Surfaces. Molecules. 2021; 26(11):3411. https://doi.org/10.3390/molecules26113411

Chicago/Turabian Style

Alexander, William A. 2021. "Methyl-Cyclohexane Methanol (MCHM) Isomer-Dependent Binding on Amorphous Carbon Surfaces" Molecules 26, no. 11: 3411. https://doi.org/10.3390/molecules26113411

APA Style

Alexander, W. A. (2021). Methyl-Cyclohexane Methanol (MCHM) Isomer-Dependent Binding on Amorphous Carbon Surfaces. Molecules, 26(11), 3411. https://doi.org/10.3390/molecules26113411

Article Metrics

Back to TopTop