Search

With the advent of new experimental techniques, measurements of individual, per-position, vapor pressure isotope effects (VPIEs) became possible. Frequently, they are in opposite directions (larger and smaller than unity), leading to the cancellation when only bulk values are determined. This progress has not been yet paralleled by the theoretical description of phase change processes that would allow for computational prediction of the values of these isotope effects. Herein, we present the first computational protocol that allowed us to predict carbon VPIEs for ethanol—the molecule of great importance in authentication protocols that rely on the precise information about position-specific isotopic composition. Only the model comprising explicit treatment of the surrounding first-shell molecules provided good agreement with the measured values of isotope effects. Additionally, we find that the internal vibrations of molecules of the model to predict isotope effects work better than the entire set of normal modes of the system.

We have performed calculations of possible oxidation pathways of 4-nitrotoluene (4NT) by permanganate anion and evaluated relative contributions of oxidation of the methyl group and aromatic ring. Although a few theory levels matched the experimental results obtained by compound specific isotope analysis (CSIA) for 4NT, they failed in reproducing results for other nitrotoluene derivatives studied previously [Wijker, R.S.; Adamczyk, P.; Bolotin, J.; Paneth, P.; Hofstetter, T.B. Environ. Sci. Technol., 2013, 47, 13459-13468]. This discrepancy prompted us to reevaluate the experimental isotopic fractionation of carbon and hydrogen for 4NT on which the relative contributions of the oxidation channels has been based. Using position specific isotope analysis (PSIA) for hydrogen isotopic fractionation we have found that the previously determined value was incorrect. Reexamination of theory levels that are in agreement with these new findings indicated that while a better agreement for this particular case can be reached, overall, the previously used B3LYP functional expressed in the 6-31+G(d,p) basis set with inclusion of the polarized continuum model (PCM) of aqueous solution remains the theoretical level of choice in modeling oxidation of nitrotoluene derivatives by permanganate.