Energy Decomposition Analysis Research Articles

How reduced are nucleophilic gold complexes?

Nucleophilic formal gold(-I) and gold(I) complexes are investigated via Intrinsic Bond Orbital analysis and Energy Decomposition Analysis, based on density functional theory calculations. The results indicate gold(0) centres engaging in electron-sharing bonding with Al- and B- based ligands. Multiconfigurational (CASSCF) calculations corroborate the findings, highlighting the gap between the electonic structures and the oxidation state formalism.

Computational insights into CO2 binding reactions by intramolecular geminal group-IV+/phosphorus- and zirconium+/group-15-based frustrated Lewis pairs.

In order to reduce global warming, there is growing interest in the design of frustrated Lewis pair (FLP) molecules for CO2 capture. This research aims to investigate the influence of group IV (M) or group 15 (G15) elements on the reactivity of intramolecular geminal M+/G15-based frustrated Lewis pair (FLP) molecules in CO2 capture. Theoretical findings suggest that M+/P-FLP, Zr+/P-FLP, Zr+/As-FLP, and Zr+/Sb-FLP can readily undergo CO2 capture reactions without difficulty. Furthermore, Zr+/As-FLP and Zr+/Sb-FLP are predicted to undergo reversible CO2 binding reactions. Interestingly, our theoretical results suggest that the M-P bond length in isolated M+/P-FLP can serve as a criterion for assessing the free activation and free reaction energy of CO2 binding. To investigate the physical factors governing the reactivity trends for the capture of CO2 reactions by intramolecular geminal M+/G15-FLP, we employed frontier molecular orbital theory, energy decomposition analysis in conjunction with natural orbitals and chemical valence, and the activation strain model. Our theoretical information can assist experimental chemists in applying key factors in the design and synthesis of novel intramolecular geminal M+/G15-FLP molecules.

Exploring π-π interactions and electron transport in complexes involving a hexacationic host and PAH guest: a promising avenue for molecular devices.

Single isolated molecules and supramolecular host-guest systems, which consist of π-π stacking interactions, are emerging as promising building blocks for creating molecular electronic devices. In this article, we have investigated the noncovalent π-π interaction and intermolecular electron charge transport involved in a series of host-guest complexes formed between a cage-like host (H6+) and polycyclic aromatic hydrocarbon (PAH) guests (G1-G7) using different quantum chemical approaches. The host (H6+) consists of two triscationic π-electron-deficient trispyridiniumtriazine (TPZ3+) units that are bridged face-to-face by three ethylene-triazole-ethylene. Our theoretical calculations show that the perylene and naphthalene inclusion complexes G7⊂H and G1⊂H have the highest and lowest interaction energies, respectively. In addition, energy decomposition analysis (EDA) indicated that the dispersion interaction term, ΔEdisp, significantly contributes to the host-guest interaction and is correlated with the existence of π-π van der Waals interaction. Using the nonequilibrium Greens function (NEGF) method in combination with density functional theory (DFT), the current-voltage (I-V) curves of the complexes were estimated. The conductance values increased when the guests were embedded inside the host cavity. Notably, the complex G7⊂H has the maximum conductance value. Overall, this study provided the electron transport of the PAH inclusion host-guest complex through π-π interaction and provided a direction for the fabrication of future supramolecular molecular devices.

Hydrogen bond properties of Se in [ROH-Se(CH3)2] complexes (R = H, CH3, C2H5): matrix-isolation infrared spectroscopy and theoretical calculations.

Se is now considered as a potential centre for hydrogen bond interactions. The hydrogen bond acceptor ability of Se has been investigated in [ROH-Se(CH3)2] complexes (R = H, CH3, and C2H5) using matrix-isolation infrared spectroscopy and electronic structure calculations. The first impression of the IR spectra of the hydrogen bond complexes of [ROH-Se(CH3)2] in N2 and Ar matrices is presented here. Moreover, no spectroscopic data are available for the [HOH-Se(CH3)2] complex. Vibrational spectra in the OH stretching region indicate the formation of the [ROH-Se(CH3)2] complex under the matrix-isolation conditions. Comparison of the experimental spectra with the simulated vibrational frequencies at different levels of theory confirms the formation of the 1 : 1 cluster of [ROH-Se(CH3)2] stabilised by O-H⋯Se hydrogen bond interactions. Multiple conformers of the [CH3OH-Se(CH3)2] complex having marginally different stabilisation energies have been predicted from electronic structure calculations and signatures of the same have been observed under the cold conditions of matrix isolation. Conformer specific assignment of the 1 : 1 cluster of [C2H5OH-Se(CH3)2] (anti and gauche forms) has been carried out in both the matrices. Concentration dependent experiments indicate the formation of higher order clusters and/or mixed clusters along with the formation of a 1 : 1 cluster for CH3OH and C2H5OH. The nature of the selenium centred hydrogen bond has been delineated using AIM, NBO and energy decomposition analysis. A comparison of similar complexes of H2O, CH3OH, and C2H5OH with O, S and Se indicates that Se is not far away in hydrogen bond acceptor ability compared to O and S.

The iron group transition-metal (Fe, Ru, Os) coordination of Se-doped graphitic carbon (Se@g-C3N4) nanostructures for the smart therapeutic delivery of zidovudine (ZVD) as an antiretroviral drug: a theoretical calculation perspective.

This study employed density functional theory (DFT) computational techniques at the ωB97XD/def2svp level of theory to comprehensively explore the electronic behavior of Fe-group transition metal (Fe, Ru, Os) coordination of Se-doped graphitic carbon (Se@g-C3N4) nanosystems in the smart delivery of zidovudine (ZVD), an antiretroviral drug. The HOMO-LUMO results of the interactions show a general reduction in energy gap values across all complexes in the following order: ZVD_Se@C3N4 < ZVD_Ru_Se@C3N4 < ZVD_Fe_Se@C3N4 < ZVD_Os_Se@C3N4. ZVD_Se@C3N4 exhibits the smallest post-interaction band gap of 3.783 eV, while ZVD_Os_Se@C3N4 presents the highest energy band gap of 5.438 eV. Results from the corrected adsorption energy (BSSE) revealed that Os_Se@C3N4 and Ru_Se@C3N4 demonstrated more negative adsorption energies of -2.67 and -2.701 eV, respectively, pointing to a more favorable interaction between ZVD and these systems, thus potentially enhancing the drug delivery efficiency. The investigation into the drug release mechanism from the adsorbents involved a comprehensive examination of the dipole moment and the influence of pH, shedding light on the controlled release of ZVD. Additionally, investigating the energy decomposition analysis (EDA) revealed that ZVD_Ru_Se@C3N4 and ZVD_Fe_Se@C3N4 exhibited the same total energy of -787.7 kJ mol-1. This intriguing similarity in their total energy levels suggested that their stability was governed by factors beyond reactivity, possibly due to intricate orbital interactions. Furthermore, analyzing the bond dissociation energies showed that all systems exhibited negative enthalpy values, indicating that these systems were exothermic at both surface and interaction levels, thus suggesting that these processes emitted heat, contributing to the surrounding thermal energy.

A DFT study of the chemical bonding properties, aromaticity indexes and molecular docking study of some phenylureas herbicides

Herbicides have implied disastrous consequences towards the environment and human health. This practice urges scientists to investigate the physical, chemical and biological properties of these substances, hence avoiding the use of the most harmful pesticides. For this purpose, the molecular structure and chemical bonding properties of phenylurea herbicides namely: Fenuron (L1), Monuron (L2), Diuron (L3) and Chlorotoluron (L4), were calculated in water using density functional theory (DFT). The energy decomposition analysis (EDA) and the extended transition state natural orbitals for chemical Valence (ETS-NOCV) reveal the dominant ionic character in Carbon-Nitrogen bond between dimethylurea fragment and benzene ring. Besides, the interaction of these herbicides with the Human Serum Albumin (HSA) was undertaken by molecular modeling. The calculation of HOMA and FLU indexes indicate that the electronic delocalization is stronger in Diuron than the other compounds, mainly caused by the two chloro substituents effects on benzene. Good correlations are found between the calculated parameters such as structural parameters, Mulliken atomic charge, topological and bonding properties and aromaticity indexes. The Vinardo molecular docking results suggest that, the binding energies of the complexes formed between HSA target and investigated compounds have the following order: L3 (-27.57 kJ/mol) &lt; L2 (-25.56 kJ/mol) &lt; L4 (-24.94 kJ/mol) &lt; L1 (-24.10 kJ/mol), which confirmed that the Fenuron is the less harmful option between the studied herbicides especially against HSA.

Sulfur-arene interactions: the S⋯π and S-H⋯π interactions in the dimers of benzofuran⋯sulfur dioxide and benzofuran⋯hydrogen sulfide.

Non-covalent interactions between sulfur centers and aromatic rings play important roles in biological chemistry. We examined here the sulfur-arene interactions between the fused aromatic heterocycle benzofuran and two prototype sulfur divalent triatomics (sulfur dioxide and hydrogen sulfide). The weakly-bound adducts were generated in a supersonic jet expansion and characterized with broadband (chirped-pulsed) time-domain microwave spectroscopy. The rotational spectrum confirmed the detection of a single isomer for both heterodimers, consistent with the computational predictions for the global minima. The benzofuran⋯sulfur dioxide dimer exhibits a stacked structure with sulfur closer to benzofuran, while in benzofuran⋯hydrogen sulfide the two S-H bonds are oriented towards the bicycle. These binding topologies are similar to the corresponding benzene adducts, but offer increased interaction energies. The stabilizing interactions are described as S⋯π or S-H⋯π, respectively, using a combination of density-functional theory calculations (dispersion corrected B3LYP and B2PLYP), natural bond orbital theory, energy decomposition and electronic density analysis methods. The two heterodimers present a larger dispersion component, but nearly balanced by electrostatic contributions.

Glycine in a basket: protonated complexes of 1,1,n,n-tetramethyl[n](2,11)teropyrenophane (n = 7, 8, 9) with glycine in the gas-phase.

Protonated complexes composed of a basket-like host molecule 1,1,n,n-tetramethyl[n](2,11)teropyrenophanes (TMnTP) (n = 7, 8, 9) and glycine as a guest were studied in the gas phase by experimental and computational methods. Blackbody infrared radiative dissociation (BIRD) experiments of [(TMnTP)(Gly)]H+ not only provided the observed Arrhenius parameters (activation energies, Eobsa, and frequency factors, A) but also suggested the existence of two populations of isomeric complexes of [(TMnTP)(Gly)]H+, termed fast dissociating (FD) and slow dissociating (SD), due to their relative BIRD rate constants. Master equation modeling was conducted to obtain the threshold dissociation energies E0 of the host-guest complexes. The relative stabilities of the most stable of the n = 7, 8, or 9 [(TMnTP)(Gly)]H+ complexes followed the trend SD-[(TM7TP)(Gly)]H+ > SD-[(TM8TP)(Gly)]H+ > SD-[(TM9TP)(Gly)]H+ by both BIRD and energy resolved sustained off-resonance irradiation collision-induced dissociation experiments (ER-SORI-CID). Computed structures and energies of [(TMnTP)(Gly)]H+ were obtained using B3LYP-D3/6-31+G(d,p) and for all TMnTP molecules, the lowest-energy structures were ones where protonated glycine was within the cavity of the TMnTP, despite the TMnTP molecules having a proton affinity 100 kJ mol-1 higher than glycine. An independent gradient model based on the Hirshfeld partition (IGMH) and natural energy decomposition analysis (NEDA) were applied to visualize and reveal the nature of interactions between hosts and guest. The NEDA analysis suggested that the polarization (POL) component which described interactions between induced multipoles contributed the most to the [(TMnTP)(Gly)]H+ (n = 7, 8, 9) complexes.

Regiodivergent (3 + 2) annulation reactions of oxyallyl cations.

We report a new method for the regiodivergent dearomative (3 + 2) reaction between 3-substituted indoles and oxyallyl cations. Access to both regioisomeric products is possible and is contingent on the presence or absence of a bromine atom on the substituted oxyallyl cation. In this way, we are able to prepare molecules that contain highly-hindered, stereodefined, vicinal, quaternary centers. Detailed computational studies employing energy decomposition analysis (EDA) at the DFT level establishes that regiochemical control arises from either reactant distortion energy or orbital mixing and dispersive forces, depending on the oxyallyl cation. Examination of the Natural Orbitals for Chemical Valence (NOCV) confirms that indole acts as the nucleophilic partner in the annulation reaction.

Understanding the reactivity and selectivity of Diels-Alder reactions involving furans.

The reactivity and endo/exo selectivity of the Diels-Alder cycloaddition reactions involving furan and substituted furans as dienes have been computationally explored. In comparison to cyclopentadiene, it is found that furan is comparatively less reactive and also less endo-selective in the reaction with maleic anhydride as the dienophile. Despite that, both the reactivity and the selectivity can be successfully modified by the presence of substituents at either 2- or 3-positions of the heterocycle. In this sense, it is found that the presence of strong electron-donor groups significantly increases the reactivity of the system while the opposite is found in the presence of electron-withdrawing groups. The observed trends in both the reactivity and selectivity are analyzed quantitatively in detail by means of the activation strain model of reactivity in combination with the energy decomposition analysis methods.

Noncovalent interactions of aromatic heterocycles: rotational spectroscopy and theoretical calculations of the thiazole-CF4 and thiazole-SF6 complexes.

The complexes of thiazole with CF4 and SF6 have been investigated by Fourier transform microwave spectroscopy and quantum chemical calculations. One rotational spectrum was observed for the thiazole-CF4 complex. Experiments and theoretical computations confirmed that the observed structure of thiazole-CF4 is primarily formed due to N⋯CCF4 interaction with the C atom of CF4 located in the plane of the thiazole ring. The rotational transitions of thiazole-CF4 exhibit A/E torsional splitting induced by the internal rotation of the -CF3 top. The potential barrier of the -CF3 internal rotation is 0.2411(1) kJ mol-1, consistent with the calculated value (∼0.3 kJ mol-1). For the thiazole-SF6 complex, one conformer with SF6 located above the thiazole ring is detected. The observed structure of thiazole-SF6 is mainly stabilized by van der Waals interactions. The energy decomposition analysis reveals that the electrostatics and dispersion are the dominant attractive contributions to the formation of thiazole-CF4 and thiazole-SF6 dimers, whereas the weight of the dispersion term becomes more significant in the thiazole-SF6 complex compared to that of the thiazole-CF4 complex.

Bioinformatic analysis of highly consumed phytochemicals as P-gp binders to overcome drug-resistance.

P-glycoprotein (P-gp) is an adenosine triphosphate (ATP)-dependent membrane efflux pump for protecting cells against xenobiotic compounds. Unfortunately, overexpressed P-gp in neoplastic cells prevents cell entry of numerous chemotherapeutic agents leading to multidrug resistance (MDR). MDR cells may be re-sensitized to chemotherapeutic drugs via P-gp inhibition/modulation. Side effects of synthetic P-gp inhibitors encouraged the development of natural products. Molecular docking and density functional theory (DFT) calculations were used as fast and accurate computational methods to explore a structure binding relationship of some dietary phytochemicals inside distinctive P-gp binding sites (modulatory/inhibitory). For this purpose, top-scored docked conformations were subjected to per-residue energy decomposition analysis in the B3LYP level of theory with a 6-31g (d, p) basis set by Gaussian98 package. Consecutive application of computational techniques revealed binding modes/affinities of nutritive phytochemicals within dominant binding sites of P-gp. Blind docking scores for best-ranked compounds were superior to verapamil and rhodamine-123. Pairwise amino acid decomposition of superior docked conformations revealed Tyr303 as an important P-gp binding residue. DFT-based induced polarization analysis revealed major electrostatic fluctuations at the atomistic level and confirmed larger effects for amino acids with energy-favored binding interactions. Conformational analysis exhibited that auraptene and 7,4',7'',4'''-tetra-O-methylamentoflavone might not necessarily interact to P-gp binding sites through minimum energy conformations. Although there are still many hurdles to overcome, obtained results may propose a few nutritive phytochemicals as potential P-gp binding agents. Moreover; top-scored derivatives may have the chance to exhibit tumor chemo-sensitizing effects.

Shedding light on the physical nature of ion pair interactions involving carba-closo-dodecaborate anions. Insights from computation.

closo-Carboranes are structures that have been studied for many decades due to their vast applicability in medicine, catalysis, and energy storage. In order to go deeper into the physics behind the interactions of oppositely charged ions, which have potential applications in electrical energy storage and conversion devices, the present work aims to shed light on the physical nature of the interactions involving (R-CB11H11-, R = H, F, CH3, CF3) and Mq+ = Li+, Na+, Mg2+, Zn2+ ions. The bonding situations are evaluated in light of both canonical Kohn-Sham energy decomposition, EDA-NOCV, and local energy decomposition, LED, analyses. Electron and hole transports are also evaluated through charge transfer integrals. The findings reveal that such complexes present not only more significant electrostatic, but also non-negligible orbital contributions. Both energy decomposition analyses EDA-NOCV and DLPNO-LED confirm that the strength of ionic pair interactions (R-CB11H11- ↔ Mq+) is much more dependent on the nature of the cation employed than on the substituent R used. The results also reveal that charge transfers are not significant in such interacting systems.

In silico analysis of SOD1 aggregation inhibition modes of tertiary amine pyrazolone and pyrano coumarin ferulate as ALS drug candidates.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease, the familial form (fALS) of which is often cognate to mutations in the antioxidant enzyme Cu/Zn superoxide dismutase 1 (SOD1) leading to misfolding and aggregation. Two small molecules, a tertiary amine pyrazolone (TAP) and a pyrano coumarin ferulate (PCF) were suggested to be ALS drug candidates following experimental observation of their ability to inhibit SOD1 protein misfolding and aggregation. The present work aims at computational investigation of these experimentally proposed drug candidates to gain insight into their mechanism of SOD1 misfolding and aggregation inhibition. On the basis of molecular docking, molecular dynamics simulation, MM-PBSA and per-residue energy decomposition analysis, we examined the specific interactions of TAP and PCF with three probable binding sites of SOD1, namely, dimeric interface cavity, W32 and, UMP binding sites. Results suggest that the binding of TAP at W32 and at UMP sites are least probable due to absence of any favorable interaction. The binding of TAP to dimeric cavity is also unstable due to strong unfavorable interactions. In case of PCF, binding at the UMP site is least probable while binding at dimeric cavity is accompanied by unfavorable interactions. PCF, however, exhibits stable binding with the W32 binding site of SOD1 by stabilizing the solvent accessible hydrophobic residues, which otherwise would have acted as contact points for aggregation. Thus the results imply that compound PCF functions as an inhibitior of SOD1 misfolding/aggregation through direct interaction with the protein SOD1 at the W32 binding site. However, TAP is likely to act as an inhibitor through a different mechanism rather than direct interaction with the protein SOD1. These results apart from reinforcing previous experimental findings, shed light on the probable mechanism of action of the proposed drug candidates.

Discriminating and understanding molecular crystal polymorphism.

Polymorph discrimination for a molecular crystal has long been a challenging task, which, nonetheless, is a major concern in the pharmaceutical industry. In this work, we have investigated polymorph discrimination on three different molecular crystals, tetrolic acid, oxalic acid, and oxalyl dihydrazide, covering both packing polymorphism and conformational polymorphism. To gain more understanding, we have performed energy decomposition analysis based on many-body expansion, and have compared the results from the XO-PBC method, that is, the eXtended ONIOM method (XO) with the periodic boundary condition (PBC), with those from some commonly used dispersion corrected density functional theory (DFT-D) methods. It is shown here that, with the XYG3 doubly hybrid functional chosen as the target high level to capture the intra- and short-range intermolecular interactions, and the periodic PBE as the basic low level to take long range interactions into account, the XO-PBC(XYG3:PBE) method not only obtains the correct experimental stability orderings, but also predicts reasonable polymorph energy ranges for all three cases. Our results have demonstrated the usefulness of the present theoretical methods, in particular XO-PBC, while highlighted the importance of a better treatment of different kinds of interactions to be beneficial to polymorph control.

η6 -Metalated Aryl Iodides in Diels-Alder Cycloaddition Reactions: Mode of Activation and Catalysis.

The potential application of η6 -metalated aryl iodides as organocatalyst has been explored by means of computational methods. It is found that the enhanced halogen bonding donor ability of these species, in comparison with their demetalated counterparts, translates into a significant acceleration of the Diels-Alder cycloaddition reaction involving cyclohexadiene and methyl vinyl ketone. The factors behind this acceleration, the endo-exo selectivity of the process and the influence of the nature of the transition metal fragment in the activity of these species are quantitatively explored in detail by means of the combination of the Activation Strain Model of reaction and the Energy Decomposition Analysis methods.

On the Permeation of Polychlorinated Dibenzodioxins and Dibenzofurans through Lipid Membranes: Classical MD and Hybrid QM/MM-EDA Analysis.

The permeation of dioxin-like pollutants, namely, chlorinated dibenzodioxins and dibenzofurans, through lipid membranes has been simulated using classic molecular dynamics (CMD) combined with the umbrella sampling approach. The most toxic forms of chlorinated dibenzodioxin and dibenzofuran, 2,3,7,8-tetrachloro-p-dibenzodioxin (TCDD) and 2,3,7,8-tetrachlorodibenzofuran (TCDF), and a dioleyl-phosphatidylcholine (DOPC) lipid membrane of 50 Å wide have been chosen for our study. The free energy profile shows the penetration process is largely favoured thermodynamically (ΔG ≈ -12 kcal/mol), with a progressively decrease of the free energy until reaching the energy minima at distances of 8 Å and 9.5 Å from the centre of the membrane for, respectively, TCDD and TCDF. At the centre of the membrane, both molecules display subtle local maxima with free energy differences of 0.5 and 1 kcal/mol with respect to the energy minima for TCDD and TCDF, respectively. Furthermore, the intermolecular interactions between the molecules and the lipid membrane have been characterized at the minima and the local maxima using hybrid quantum mechanics/molecular mechanics energy decomposition analysis (QM/MM-EDA). Total interaction energies of -17.5 and -16.5 kcal/mol have been found at the energy minima for TCDD and TCDF, respectively. In both cases, the dispersion forces govern the molecule-membrane interactions, no significant changes have been found at the local maxima, in agreement with the classical free energy profile. The small differences found in the results obtained for TCDD and TCDF point out that the adsorption and diffusion processes through the cell membrane are not related to the different toxicity shown by these pollutants.

Redox Active cAAC-Fluorene/Indene Systems Displaying Solvatochromism, Green Luminescence and pH Sensing: Functionalization of Fluorenyl/Indenyl Rings with Radical Carbene.

Two new series of air stable compounds of cAACX = fluorene/indene (X = Me2 , Et2 , Cy) [cAAC = cyclic (alkyl) amino carbene] have been isolated and well characterized by X-ray single crystal diffraction, photoluminescence, cyclic voltammogram (CV) and electron paramagnetic resonance (EPR) studies. Fluorescence studies reveals green light emission of cAAC bonded fluorene, whereas free fluorene generally displays a violet emission. Interestingly, the sterically crowded cAAC-fluorene analogue display solvatochromism and CF3 CO2 H sensing in solution. CV of the these compounds show a quasi-reversible electron transfer process, indicating the functionalization of fluorene/indene with radical anionic form of carbene, confirmed by CV/EPR measurements. DFT/TDDFT calculations and energy decomposition analysis coupled with natural orbital for chemical valence (EDA-NOCV) have been carried out to study different aspects of bonding and electronic transitions. Such a class of redox active and thermally stable organic molecules may be suitable for molecule based spin memory devices in future.

C–O Bond Activation in Mononuclear Lanthanide Oxocarbonyl Complexes OLn(η2-CO) (Ln = La, Ce, Pr, and Nd)

Fundamental investigation of metal-CO interactions is of great importance for the development of high-performance catalysts to CO activation. Herein, a series of side-on bonded mononuclear lanthanide (Ln) oxocarbonyl complexes OLn(η2-CO) (Ln = La, Ce, Pr, and Nd) have been prepared and identified in solid argon matrices. The complexes exhibit uncommonly low C-O stretching bands near 1630 cm-1, indicating remarkable C-O bond activation in these Ln analogues. The η2-CO ligand in OLn(η2-CO) can be claimed as an anion on the basis of the experimental observations and quantum chemistry investigations, although the CO anion is commonly considered to be unstable with electron auto-detachment. The CO activation in OLn(η2-CO) is attributed to the photoinduced intramolecular charge transfer from LnO to CO rather than the generally accepted metal → CO π back-donation, which conforms to the traditional Dewar-Chatt-Duncanson motif. Energy decomposition analysis combined with natural orbitals for chemical valence calculations demonstrates that the bonding between LnO and η2-CO arises from the combination of dominant ionic forces (>76%) and normal Lewis "acid-base" interactions. The fundamental findings provide guidelines for the catalyst design of CO activation.

Computer-aided Screening of High Performance Covalent-Organic Frameworks for Removal of SO2 from Flue Gases

The combustion of fossil fuels can produce sulfur dioxide (SO2), which is a major source of the environmental problems such as acid rain, photochemical smog and ozone depletion. To solve these problems, it is urgent to develop new porous materials that are conducive to remove acid gases such as SO2. In this work, the adsorption and separation properties of SO2 from flue gases (SO2CO2N2) by covalent organic frameworks (COF) were systematically studied based on the grand canonical Monte Carlo simulations. The results of the structure-activity relationships show that the separation performance of SO2 is related to the pore volume (Vpore) and accessible surface area (ASA) of COF materials to some extent. Some promising candidates with high SO2 adsorption capacity and separation selectivity were identified by comprehensive evaluation metrics, and it can found that the existence of appropriate functional group and the small channels are more conducive to gas adsorption. In addition, the energy decomposition and IGM analysis indicate the nature of the adsorption of SO2 on CF3-TFP-TTA is mainly originated from the electrostatic attraction. These results can provide theoretical references for the design and synthesis of novel high-performance adsorptive separation material.