Transition State Calculations Research Articles

Permeability of gas molecules through sub-nanometer nitrogen-terminated porous graphene membranes: A DFT study

We present a systematic theoretical investigation of gas permeability through single-layer nanoporous graphene membranes with N-terminated (pyridinic and mixed pyridinic/pyrrolic) pores for He, H2, N2, O2, CO, CO2, H2O, and CH4. Our study is based on density functional theory transition-state calculations and the kinetic theory of gasses. We systematically evaluate pores of varying sizes and shapes up to 5.7 Å in diameter. According to our findings, membranes with pores in the size regime (4.5 Å - 5.0 Å) may exhibit industrially acceptable permeance to several gas molecules and advanced selectivity. Notably, we found that, among a few others, a pore with the same topological features at the pore boundary as those of the characteristic pore of g-C3N4 carbon nitride, is particularly promising. In addition, membranes with larger pores, ∼5.5 Å - 5.7 Å, can effectively separate CH4 from the other molecules. Our findings support the advanced potential of single-layer graphene membranes with nitrogen-terminated sub-nanometer pores, for gas separation applications.

Quantum Chemical Characterization of Rotamerism in Thio-Michael Additions for Targeted Covalent Inhibitors.

Myotonic dystrophy type I (DM1) is the most common form of adult muscular dystrophy and is a severe condition with no treatment currently available. Recently, small-molecule ligands have been developed as targeted covalent inhibitors that have some selectivity for and covalently inhibit cyclin-dependent kinase 12 (CDK12). CDK12 is involved in the transcription of elongated RNA sections that results in the DM1 condition. The covalent bond is achieved after nucleophilic addition to a Michael acceptor warhead. Previous studies of the conformational preferences of thio-Michael additions have focused on characterizing the reaction profile based on the distance between the sulfur and β-carbon atoms. Rotamerism, however, has not been investigated extensively. Here, we use high-level quantum chemistry calculations, up to coupled cluster with single, double, and perturbative triple excitations [CCSD(T)], to characterize the nucleophilic addition of an archetypal nucleophile, methanethiolate, to various nitrogen-containing Michael acceptors which are representative of the small-molecule covalent inhibitors. By investigating the structural, energetic, and electronic properties of the resulting enolates, as well as their reaction profiles, we show that synclinal additions are generally energetically favored over other additions due to the greater magnitude of attractive noncovalent interactions permitted by the conformation. The calculated transition states associated with the addition process indicate that synclinal addition proceeds via lower energetic barriers than antiperiplanar addition and is the preferred reaction pathway.

Efficiency and Mechanism of Catalytic Siloxane Exchange in Vitrimer Polymers: Modeling and Density Functional Theory Investigations.

Of late, siloxane-containing vitrimers have gained significant interest due to their fast dynamic characteristics over a reasonable temperature range (180-220 °C), making them well-suited for diverse applications. The exchange reaction pathway in the siloxane vitrimers is accountable for the covalent adaptive network, with the reaction's effectiveness being regulated by either organic or organometallic catalysts. However, directly studying the exchange reaction pathway in the bulk phase using experimental approaches is challenging because of the intricate and interconnected structure of these vitrimers. Here, we perform comprehensive density functional theory (DFT) and experimental investigations to discover the detailed catalytic efficacy of siloxane exchange and provide direction for the reaction process using a 1,5,7-triazabicyclo[4.4.0]dec-5-ene (TBD) catalyst. The calculated transition barrier energy and catalytic efficiency of hexamethyldisiloxane and dihydroxy-dimethylsilane exchange derived from the nudged elastic band with transition-state calculations strongly agree with the experimental findings. In addition, Fukui indices, along with partial charges, are employed to evaluate the nucleophilic and electrophilic behaviors of silanol and siloxane molecules. Our analysis revealed that by utilizing the Fukui indices of both the acid and the base, we can make an approximate estimation of the respective kinetics of the SN2 process in the siloxane exchange reaction mechanism. These findings establish a foundation for comprehending a crucial aspect of the exchange mechanism in siloxane vitrimer systems and could aid in the development of novel catalysts.

Disentangling the Conformational Space and Structural Preferences of Tetrahydrofurfuryl Alcohol Using Rotational Spectroscopy and Computational Chemistry.

The influence of the hydroxymethyl (CH2OH) group on the tetrahydrofuran (THF) ring structure was investigated by disentangling the gas phase conformational landscape of the sugar analogue tetrahydrofurfuryl alcohol (THFA). By combining rotational spectroscopy (6-20 GHz) and quantum chemical calculations, transitions corresponding to two stable conformers of THFA and their 13C isotopologues were observed and assigned in the rotational spectrum. The positions of the C atoms were precisely determined to unambiguously distinguish between nearly isoenergetic pairs of conformers that differ in their ring configurations: envelope (E) versus twist (T). The rotational spectrum confirms that the E ring geometry is favoured when the CH2OH fragment lies gauche (-) to the THF backbone (OCCO ~-60°) whereas the T form is more stable for the gauche (+) alignment of the substituent (OCCO ~+60°). The observed spectral intensities suggest that conformational relaxation of the THF geometry (E↔T) to the more stable form readily occurs within the pairs of g- and g+ conformers which is consistent with the low barriers (1.5-1.7 kJ mol-1) for conversion determined via transition state calculations. Insights into the intramolecular hydrogen bonding and other weak interactions stabilizing the lowest energy structures of THFA were derived and rationalized using non-covalent interaction analyses.

Computational Insights into the Mechanism of Boyland-Sims Oxidation Reaction

The Boyland-Sims oxidation is the reaction between peroxydisulfate ions and arylamines to form the aryl amine ortho sulfate with minor para isomer. The long-standing mechanism involves a nucleophilic attack by the amine on peroxide oxygen to form arylhydroxylamine-O-sulfonate intermediate, which subsequently rearranges to the arylamine o-sulfate. Marjanović et al. have challenged this long-standing mechanism, proposing a nitrenium ion intermediate as the reactive species instead of the uncharged amine. The nitrenium ion proposal was based on relatively low-level calculations. We previously computed the energetics of reaction intermediates in both mechanisms for aniline, 2,4-dinitroaniline, and N,N-dimethylaniline using high-level of density functional theory (B3LYP/6-311++G**) calculations. We present here new computations using density functional theory (B3LYP/6-311++G**) to model both arylhydroxylamine-O-sulfonate and nitrenium ion pathways in a series of aromatic amines. The transition state calculations revealed two possible energetically feasible pathways for the rearrangement of arylhydroxylamine-O-sulfonate to arylamine o-sulfate. Additionally, we computationally modeled the feasibility of the reaction of aniline and tetrathionate to yield 2-aminothiophenol product via Boyland-Sims Oxidation reaction. The results indicate tetrathionate ion is less reactive with aniline than the peroxydisulfate ion.

Exploring the conformational landscape, hydrogen bonding, and internal dynamics in the diallyl ether and diallyl sulfide monohydrates.

The conformational spaces of the diallyl ether (DAE) and diallyl sulfide (DAS) monohydrates were explored using rotational spectroscopy from 6 to 19GHz. Calculations at the B3LYP-D3(BJ)/aug-cc-pVTZ level suggested significant differences in their conformational behavior, with DAE-w exhibiting 22 unique conformers and DAS-w featuring three stable structures within 6 kJ mol-1. However, only transitions from the lowest energy conformer of each were experimentally observed. Spectral analysis confirmed that binding with water does not alter the conformational preference for the lowest energy structure of the monomers, but it does influence the relative stabilities of all other conformers, particularly in the case of DAE. Non-covalent interaction and quantum theory of atoms in molecules analyses showed that the observed conformer for each complex is stabilized by two intermolecular hydrogen bonds (HBs), where water primarily interacts with the central oxygen or sulfur atom of the diallyl compounds, along with secondary interactions involving the allyl groups. The nature of these interactions was further elucidated using symmetry-adapted perturbation theory, which suggests that the primary HB interaction with S in DAS is weaker and more dispersive in nature compared to the primary HB in DAE. This supports the experimental observation of a tunneling splitting exclusively in the rotational spectrum of DAS-w, as the weaker contact allows water to undergo internal motions within the complex, as shown based on calculated transition state structures for possible tunneling pathways.

Pericyclic Cascade Reactions Affording Thienyl- and Benzo[b]thienyl-Fused Architectures: A Synthetic and Density Functional Theory Analysis of Asynchronous Diels-Alder Reactions of Thioarylmaleimides.

Herein, we report the synthesis and characterization of a series of thioarylmaleimides and their varied propensity toward highly selective domino Diels-Alder (D-A)/rearrangement, D-A/ene/elimination, and D-A/oxidation reactions to give three types of thienyl-fused architectures. Stereochemical assignment was achieved using a combination of nuclear magnetic resonance (NMR) studies, gauge independent atomic orbital (GIAO) NMR chemical shift calculations, and DP4+ analysis. Transition-state calculations support an asynchronous concerted mechanism and provide support for rationalizing the observed regio- and stereoselectivity.

Mechanistic study of the atomic layer deposition of cobalt: a combined mass spectrometric and computational approach.

Cobaltcarbonyl-tert-butylacetylene (CCTBA) is a conventional precursor for the selective atomic layer deposition of Co onto silicon surfaces. However, a limited understanding of the deposition mechanism of such cobalt precursors curbs rational improvements on their design for increased efficiency and tuneable selectivity. The impact of using a less reactive internal alkyne instead of a terminal alkyne was investigated using experimental and computational methods. Using electrospray-ionization mass spectrometry, the formation of CCTBA analogs and their gas phase decomposition pathways were studied. Decomposition experiments show very similar decomposition pathways between the two complexes. The internal alkyne dissociates from the Co complex at slightly lower energies than the terminal alkyne, suggesting that an internal alkynyl ligand may be more suited to low temperature ALD. In addition, transition state calculations using the nudged elastic band method confirm an increased reaction barrier between the internal alkyne and the Si-H surface bonds on Si(111). These results suggests that using a less reactive internal alkyne will result in fewer embedded carbon impurities during deposition onto Si wafers. DFT calculations using the PBE functional and periodic boundary conditions also predict increased surface binding with the metal centers of the internal alkynyl complex.

High energy barrier hydroxyl radical dissociation mechanism of a low shock sensitivity dihydroxylammonium 5,5'-bistetrazole-1,1'-diolate (TKX-50) explosive.

As a representative of the new generation of high-energy explosives, TKX-50 has attracted widespread attention due to its remarkably low sensitivity toward shock. However, the reported decomposition barriers of TKX-50 (∼37 kcal mol-1) are comparable to those of commonly used explosives. The mechanism of its low shock sensitivity remains unclear. In this study, using an ab initio molecular dynamics method combined with a multiscale shock simulation technique and transition state calculations (at the B2PLYP-D3/Def2TZVP level), we discovered an unconventional reaction pathway of TKX-50 under shock, and its rate-controlling step is the dissociation of the hydroxyl radical (OH) from the anion ring after proton transfer, followed by ring rupture and the production of H2O and N2. The barrier for this OH dissociation reaction is as high as 51.9 kcal mol-1. In contrast, under thermal stimuli, TKX-50 prefers to open rings directly after proton transfer without losing the OH. The corresponding barrier is 35.4 kcal mol-1, which is in good agreement with previous studies. The reason for the unconventional reaction pathway of TKX-50 under shock may be the suppression of anion ring opening in thermal decomposition by steric hindrance upon shock compression. In addition, the dominant N2 generation pathway under shock releases less energy than pyrolysis which further explains the low shock sensitivity of TKX-50. This study comprehensively elucidates the different reaction mechanisms of TKX-50 under thermal and shock conditions and proposes a crucial reaction pathway leading to its low shock sensitivity. These findings will contribute to the understanding and application of tetrazole anionic energetic salts.

Modular synthesis of aryl amines from 3-alkynyl-2-pyrones.

The synthesis of aryl amines from 3-alkynyl-2-pyrones and various amines is described. Mechanistically, the aryl amines are proposed to arise from the 3-alkynyl-2-pyrone substrates through their selective opening in a 1,6-fashion by secondary amines followed by decarboxylation and an unexpected rearrangement. The proposed mechanism is supported by quantum chemical transition-state calculations, which are consistent with the regiochemical outcome. The scope of this transformation spans a variety of 3-alkynyl-2-pyrones and a range of secondary amines. The influence of the secondary amine coupling partners on reaction efficiency was elucidated through data-driven modeling as well as scope exploration. These latter studies revealed that the steric bulk of the secondary amine coupling partner under the reaction conditions serves as a strong indicator of overall reaction efficiency.

Molecular field analysis in half-titanocene complexes: computational study towards data-driven in silico optimization of single-site olefin polymerization catalysts

The control of multiple reaction outcomes in olefin polymerization has been demonstrated through data-driven catalyst optimization based on molecular field analysis, using data collected from DFT-based transition-state calculations.

Small molecule activation by sila/germa boryne species.

The inability of p-block elements to participate in π-backbonding restricts them from activating small molecules like CO, H2 , and so forth. However, the development of the main group metallomimetics became a new pathway, where the main-group elements like boron can bind and activate small molecules like CO and H2 . The concept of the frustrated Lewis pair, Boron-Boron multiple bonds, and borylene are previously illustrated. Some of these reported classes of boron species can mimic the jobs of the metal complexes. Hence, we have theoretically studied the binding of CO/N2 molecules at B-center of elusive species like sila/germa boryne stabilized by donor base ligands (cAAC)BE(Me)(L), where E  Si, L  cAACMe , NHCMe , PMe3 , E  Ge, L  cAACMe and (NHCMe )BE(Me)(cAACMe )). The substitutional analogues of (cAACR )BSiR1 (cAAC) and E  P, L  cAACMe ) have been studied by density functional theory (DFT), natural bond orbital, QTAIM calculations and energy decomposition analysis (EDA) coupled with natural orbital for chemical valence (NOCV) analyses. The computed bond dissociation energy and inner stability analyses by the EDA-NOCV method showed that the CO molecule can bind at the B-center of the above-mentioned species due to stronger σ-donor ability while binding of N2 has been theoretically predicted to be weak. The energy barrier for the CO binding is estimated to be 13-14 kcal/mol by transition state calculation. The change of partial triple bond character to single bond nature of the BSi bond and the bending of CBSi bond angle of sila-boryne species are the reason for the activation energy. Our study reveals the ability of such species to bind and activate the CO molecule to mimic the transition metal-containing complexes. We have additionally shown that binding of Fe(CO)4 and Ni(CO)3 is feasible at Si-center after binding of CO at the B-center.

OH Roaming during the Ozonolysis of α-Pinene: A New Route to Highly Oxygenated Molecules?

The formation of low-volatility organic compounds in the ozonolysis of α-pinene, the dominant atmospheric monoterpene, provides an important route to aerosol formation. In this work, we consider a previously unexplored set of pathways for the formation of highly oxygenated molecules in α-pinene ozonolysis. Pioneering, direct experimental observations of Lester and co-workers have demonstrated a significant production of hydroxycarbonyl products in the dissociation of Criegee intermediates. Theoretical analyses indicate that this production arises from OH roaming-induced pathways during the OO fission of the vinylhydroperoxides (VHPs), which in turn come from internal H transfers in the Criegee intermediates. Ab initio kinetics computations are used here to explore the OH roaming-induced channels that arise from the ozonolysis of α-pinene. For computational reasons, the calculations consider a surrogate for α-pinene, where two spectator methyl groups are replaced with H atoms. Multireference electronic structure calculations are used to illustrate a variety of energetically accessible OH roaming pathways for the four VHPs arising from the ozonolysis of this α-pinene surrogate. Ab initio transition-state theory-based master equation calculations indicate that for the dissociation of stabilized VHPs, these OH roaming pathways are kinetically significant with a branching that generally increases from ∼20% at room temperature up to ∼70% at lower temperatures representative of the troposphere. For one of the VHPs, this branching already exceeds 60% at room temperature. For the overall ozonolysis process, these branching ratios would be greatly reduced by a limited branching to the stabilized VHP, although there would also be some modest roaming fraction for the nonthermal VHP dissociation process. The strong exothermicities of the roaming-induced isomerizations/additions and abstractions suggest new routes to fission of the cyclobutane rings. Such ring fissions would facilitate further autoxidation reactions, thereby providing a new route for producing highly oxygenated nonvolatile precursors to aerosol formation.

High-temperature effects for transition state calculations in solids.

Transition state calculation is a critical technique to understand and predict versatile dynamical phenomena in solids. However, the transition state results obtained at 0K are often utilized for the prediction or interpretation of dynamical processes at high temperatures, and the error bars of such an approximation are largely unknown. In this benchmark study, all the major temperature effects, including lattice expansion, lattice vibration, electron excitation, and band-edge shift, are evaluated with first-principles calculations for defect diffusion in solids. With the inclusion of these temperature effects, the notable discrepancies between theoretical predictions at 0K and the experimental diffusivities at high temperatures are dramatically reduced. In particular, we find that lattice expansion and lattice vibration are the dominant factors lowering the defect formation energies and hopping barriers at high temperatures, but the electron excitation exhibits minor effects. In sharp contrast to typical assumptions, the attempt frequency with lattice expansion and vibration varies significantly with materials: several THz for aluminum bulk but surprisingly over 500THz for 4H-SiC. For defects in semiconductors, the band-edge shift is also significant at high temperatures and plays a vital role in defect diffusion. We expect that this study would help accurately predict the dynamical processes at high temperatures.

Fully Automatized Optimization of Ring-Opening Reactions in Lactone Derivatives via Two-Step Machine Learning.

Cyclization and cycloreversion of organic compounds are fundamental kinetic processes in the design of functional molecules, molecular machines, nanoscale sensors, and switches in the field of molecular and nanoelectronics. We present a fully automatic computational platform for the design of a class of five- and six-membered ring lactones by optimizing the ring-opening reaction rate. Starting from a minimal initial parent set, our algorithm generates iteratively cascades of pools of candidate lactone derivatives where optimization and down-selection are performed without human supervision. We employ the density functional theory combined with the transition state theory to elucidate the exact mechanism leading to the lactone ring-opening reaction. On the basis of the analysis of the reaction pathway and the frontier molecular orbitals, we identify a simple descriptor that can easily correlate with the reaction rate. Consequently, we can omit computationally expensive transition state calculations and deduce the reaction rate from simple ground-state and ionic calculations. To accelerate the platform, we use a data set of the order of 800 molecules to train machine learning models for the prediction of targeted chemical properties, reducing the computational time by a 90% factor. We developed an evolutionary algorithm capable of generating data sets 3 orders of magnitude larger than the initial parent set. Thus, we can explore a large domain of chemical space using minimal computational effort. Our entire platform is modular, and our current implementation for lactone can be further generalized to more complex systems via substitution of the quantum chemical and fingerprinting modules.

Dissociative Electron Attachment to 5-Iodo-4-thio-2'-deoxyuridine: A Potential Radiosensitizer of Hypoxic Cells.

In the search for effective radiosensitizers for tumor cells, halogenated uracils have attracted more attention due to their large cross section for dissociation upon the attachment of low-energy electrons. In this study, we investigated dissociative electron attachment (DEA) to 5-iodo-4-thio-2'-deoxyuridine, a potential radiosensitizer using a crossed electron-molecule beam experiment coupled with quadrupole mass spectrometry. The experimental results were supported by calculations on the threshold energies of formed anions and transition state calculations. We show that low-energy electrons with kinetic energies near 0 eV may effectively decompose the molecule upon DEA. The by far most abundant anion observed corresponds to the iodine anion (I-). Due to the associated bond cleavage, a radical site is formed at the C5 position, which may initiate strand break formation if the molecule is incorporated into a DNA strand. Our results reflect the conclusion from previous radiolysis studies with the title compound, suggesting its potential as a radiosensitizer.

Molecular Insights into Antioxidant Efficiency of Melanin: A Sustainable Antioxidant for Natural Rubber Formulations.

N-(1,3-Dimethyl butyl)-N'-phenyl-p-phenylenediamine (6-PPD) is a worldwide antioxidant commonly added to delay the thermo-oxidative degradation of tire rubbers. Unfortunately, 6PPD and its transformation product 6PPD-quinone are toxic to aquatic organisms (e.g., coho salmon). Herein, we explore the free radical scavenging activity and protective mechanism of melanin (MLN) on natural rubber's (NR's) oxidative resistance using molecular dynamics (MD) and quantum mechanical (QM) calculations. The relationship between the molecular structure and the chemical nature of the antioxidant molecules via transition state calculations is explored to unravel the reaction mechanisms of antioxidants interacting with peroxy radicals (ROO·) of NR with the estimation of reaction barriers. Following this, the radical scavenging activity of antioxidants was quantified via a hydrogen atom transfer mechanism and bond dissociation energy calculations. Parallel MD simulations were considered to study the interfacial interactions of antioxidant molecules with polymer chains and fillers with a quantifiable structure-property correlation. Given these results, the nanocomposite (NR-MLN-SiO2) with natural antioxidant melanin manifested outstanding antioxidant properties by preferentially bagging the ROO· radicals, thus improving NR's thermal-oxidative aging relative to 6-PPD. The MD results revealed that the intermolecular interactions at the NR/antioxidant interface benefited the antioxidant MLN to bind tightly to the NR in NR-MLN-SiO2 composite, thus exhibiting improved dispersion, O2 barrier properties, and thermo-oxidative stability, which could extend the service life of NR products (e.g., tires). In addition, as a sustainable antioxidant, MLN could replace toxic antioxidants like 6-PPD. More importantly, the QM/MD simulations provided a fundamental understanding of the mechanistic pathways of antioxidant molecules in NR composites, which are conducive to designing high-performance and sustainable green elastomers.

Excellent Anode Performance of N-, P-, and As-Doped Graphdiynes for Lithium-Ion Batteries.

Recently, graphdiyne (GDY) as a two-dimensional planar carbon allotrope has received significant research attention in the fields of rechargeable batteries, catalysis, biomedicine, and so forth. However, the theoretical capacity of a perfect GDY anode is only 744 mA h/g in the configuration of LiC3, encouraging further efforts to increase the capacity. In this study, we explore the anode performance of N-, P-, and As-doped GDYs by using first-principles calculations. Ab initio molecular dynamics simulations show that the doped GDYs can remain stable at 1000 K, indicating good thermal stability. With the loss of part acetylenic linkages, the rhomboid-like pores produce more Li sites, and the theoretical capacities reach 2209, 2031, and 1681 mA h/g for the N-, P-, and As-doped GDYs, respectively. In addition, the transition-state calculations indicate that the Li diffusion barriers of the three doped GDYs are similar to the perfect GDY. This study demonstrates that doping is an effective strategy to improve the anode performance of GDY.

Theoretical investigation of hydroxylated analogues of valinomycin as potassium transporter

Valinomycin is a potent ionophore known for its ability to transport potassium ions across biological membranes. The study focuses on the hydroxylated analogues of valinomycin (HyVLMs) and compares their energy profiles and capabilities for transporting potassium ions across phospholipid membranes. Using metadynamics, we investigated the energy profiles of wildtype valinomycin (VLM_1) and its three hydroxylated analogues (VLM_2, VLM_3, and VLM_4). We observed that all analogues exhibited energy maxima in the centre of the membrane and preferred positions below the phospholipid heads. Furthermore, the entry barriers for membrane penetration were similar among the analogues, suggesting that the hydroxyl group did not significantly affect their passage through the membrane. Transition state calculations provided insights into the ability of valinomycin analogues to capture potassium ions, with VLM_4 showing the lowest activation energy and VLM_2 displaying the highest. Our findings contribute to understanding the mechanisms of potassium transport by valinomycin analogues and highlight their potential as ionophores. The presence of the hydroxyl group is of particular importance because it paves the way for subsequent chemical modifications and the synthesis of new antiviral agents with reduced intrinsic toxicity.

DFT studies on exposure of sulfur impregnated and sulfur functionalized activated carbon to Hg0 vapors

For removal of poisonous vapor emissions such as mercury, it is necessary to select suitable materials based on an understanding of their properties and interactions with the vapor. As mercury has a high affinity towards sulfur, it’s adsorption using sulfur-impregnated activated carbon was explored in this study. The impregnation of sulfur on activated carbon followed by the adsorption of Hg0 vapors was computationally investigated using DFT simulations. Sulfur adsorption was investigated on activated carbon with armchair edge, zigzag edge, and graphene surface. Sulfur adsorption was investigated on activated carbon with edge functional groups such as hydroxyl and carboxylic acid. Activated carbon with edge functional groups such as sulfonic acid, sulfenic acid, and sulfinic acid was further investigated for the adsorption of Hg0 vapors. Among the edge functional groups on the activated carbon, the hydroxyl group was most favored for sulfur adsorption and, subsequently, Hg0 vapors. This was quantified in terms of shortest bond lengths, strongest binding energies, and maximum charge transfer. Among the sulfur-containing functional groups on activated carbon, sulfenic acid was the most favored for the adsorption of Hg0 vapors. Transition state calculations were carried out, and a reaction pathway was proposed for the adsorption of Hg0 on sulfur-impregnated activated carbon.