Computational Study of Adsorption and Separation of CO(2), CH(4), and N(2) by an rht-Type Metal-Organic Framework

In this work, a computational study is performed to evaluate the adsorption-based separation of CO(2) from flue gas (mixtures of CO(2) and N(2)) and natural gas (mixtures of CO(2) and CH(4)) using microporous metal organic framework Cu-TDPAT as a sorbent material. The results show that electrostatic interactions can greatly enhance the separation efficiency of this MOF for gas mixtures of different components. Furthermore, the study also suggests that Cu-TDPAT is a promising material for the separation of CO(2) from N(2) and CH(4), and its macroscopic separation behavior can be elucidated on a molecular level to give insight into the underlying mechanisms. On the basis of the single-component CO(2), N(2), and CH(4) isotherms, binary mixture adsorption (CO(2)/N(2) and CO(2)/CH(4)) and ternary mixture adsorption (CO(2)/N(2)/CH(4)) were predicted using the ideal adsorbed solution theory (IAST). The effect of H(2)O vapor on the CO(2) adsorption selectivity and capacity was also examined. The applicability of IAST to this system was validated by performing GCMC simulations for both single-component and mixture adsorption processes.     

Probing the dynamics of CO2 and CH4 within the porous zirconium terephthalate UiO-66(Zr): a synergic combination of neutron scattering measurements and molecular simulations

Quasi-elastic neutron scattering (QENS) measurements combined with molecular dynamics (MD) simulations were conducted to deeply understand the concentration dependence of the self- and transport diffusivities of CH(4) and CO(2), respectively, in the humidity-resistant metal-organic framework UiO-66(Zr). The QENS measurements show that the self-diffusivity profile for CH(4) exhibits a maximum, while the transport diffusivity for CO(2) increases continuously at the loadings explored in this study. Our MD simulations can reproduce fairly well both the magnitude and the concentration dependence of each measured diffusivity. The flexibility of the framework implemented by deriving a new forcefield for UiO-66(Zr) has a significant impact on the diffusivity of the two species. Methane diffuses faster than CO(2) over a broad range of loading, and this is in contrast to zeolites with narrow windows, for which opposite trends were observed. Further analysis of the MD trajectories indicates that the global microscopic diffusion mechanism involves a combination of intracage motions and jump sequences between tetrahedral and octahedral cages.     

Rate‐Based Modeling of CO2 Absorption into Piperazine‐Activated Aqueous N‐Methyldiethanolamine Solution: Kinetic and Mass Transfer Analysis

A comprehensive two‐dimensional mathematical model based on surface renewal theory has been developed to analyze the CO₂ absorption into piperazine (PZ)‐activated aqueous N‐methyldiethanolamine (MDEA) solvent by taking into account the structured packed bed column hydraulics, mass transfer resistances, and chemical reactions. The modeling results have been validated with the experimental data reported in the literature, and they have been found to be in good agreement with the experimental results. The effects of amine concentration, liquid temperature, initial CO₂ partial pressure, liquid flow rate, and CO₂ loading on the mass transfer performance have been evaluated in terms of overall mass transfer coefficient (K _Ga_v). The overall mass transfer coefficient and absorption flux of CO₂ into aqueous MDEA+PZ blended solution have been calculated over the CO₂ partial pressure range of 4–16 kPa, temperature range of 298–333 K, and solvent concentration of 1–3 M. To evaluate the performance of different solvents on separation process, some common industrial chemical absorbents including monoethanolamine (MEA), diethanolamine (DEA), triethylamine (TEA), MDEA and PZ were compared with a MDEA+PZ blended solution. The results indicate that CO₂ absorption reaction with PZ is faster than that with MDEA, but also adding small amounts of PZ as a promoter to MDEA solvents improves significantly the absorption rate. The results show that CO₂ absorption reaction with the MDEA+PZ blended solution is faster than that with TEA and MDEA, also comparable with DEA, but slower than those with MEA and PZ. The modeling results illustrate that the K _Ga_v enhances with increasing the solvent concentration, liquid temperature, and liquid flow rate, but reduces with increasing the CO₂ loading and initial CO₂ partial pressure. In addition, the reaction kinetics in terms of enhancement factor was found to decrease as the CO₂ loading enhances and increase as the operating temperature rises.     

Modeling of CO2 solubility and carbamate concentration in DEA,MDEA and their mixtures using the Deshmukh–Mather model

The equilibrium of CO₂ and carbamate concentration data for the absorption of CO₂ in aqueous solutions of single and mixed amines was analyzed using the Deshmukh–Mather model. Data on CO₂ loading in aqueous solutions of DEA and MDEA and their mixtures at various temperature (303–323 K) and CO₂ partial pressure (0.09–100 kPa) together with carbamate concentrations in case of DEA and its mixtures with MDEA were fitted simultaneously to generate the different interaction parameters required to calculate the activity coefficients in the model. Using the generated interaction parameters, the model was applied to correlate the CO₂ loading in solutions of DEA and MDEA and their mixtures reported in the literature as well as those obtained in our laboratory and was found to be able to give a good estimation of CO₂ loading and carbamate concentration over a wide range of operating conditions in both single and mixed amine solutions.     

The capture of carbonyl sulfide by N-methyldiethanolamine: A systematic density functional theory investigation

In this contribution, the mechanism of carbonyl sulfide (COS) absorption by N-methyldiethanolamine (MDEA) aqueous solution was explored via theoretical computations. Detailed reaction mechanisms were analyzed using density functional theory (DFT) calculations at the B3LYP-D3 level of theory. In total, four different pathways for COS absorption by MDEA have been considered. The most favorable pathway for the removal of COS is a three-step mechanism including the hydrolysis, proton transfer, and dissociation of CO₂, and hydrolysis is the rate-determining step. The mechanisms of the COS absorption by different amines were investigated, and the calculated results suggest that the total energy barrier for the COS absorption by MDEA is comparable to that by monoethanolamine (MEA), diethanolamine (DEA), and diisopropylamine (DIPA), indicating the COS absorption by all the four amines are feasible, while MDEA gives a better performance in terms of thermodynamics.     

二氧化碳: 一种同时提供亲核中心保护和亲电中心活化的试剂 XI:2-取代四氢噻唑的新合成方法

四氢噻唑和正丁基锂反应生成的N-锂代四氢噻唑(2)与二氧化碳作用生成N-羟酸锂盐3; 3然后在叔丁基锂作用下失去2-位质子形成相应的2-位碳阴离子4; 该碳阴离子和一系列亲电试剂反应, 并在酸性条件下除去保护基团得产物-2-取代四氢噻唑7. 在此过程中, 二氧化碳既是氨基的保护基团, 同时又活化了2-位亲电中心.此方法不需分离各步中间体, 产率中等.     

Mechanisms and kinetics for sorption of CO2 on bicontinuous mesoporous silica modified with n-propylamine

We studied equilibrium adsorption and uptake kinetics and identified molecular species that formed during sorption of carbon dioxide on amine-modified silica. Bicontinuous silicas (AMS-6 and MCM-48) were postsynthetically modified with (3-aminopropyl)triethoxysilane or (3-aminopropyl)methyldiethoxysilane, and amine-modified AMS-6 adsorbed more CO(2) than did amine-modified MCM-48. By in situ FTIR spectroscopy, we showed that the amine groups reacted with CO(2) and formed ammonium carbamate ion pairs as well as carbamic acids under both dry and moist conditions. The carbamic acid was stabilized by hydrogen bonds, and ammonium carbamate ion pairs formed preferably on sorbents with high densities of amine groups. Under dry conditions, silylpropylcarbamate formed, slowly, by condensing carbamic acid and silanol groups. The ratio of ammonium carbamate ion pairs to silylpropylcarbamate was higher for samples with high amine contents than samples with low amine contents. Bicarbonates or carbonates did not form under dry or moist conditions. The uptake of CO(2) was enhanced in the presence of water, which was rationalized by the observed release of additional amine groups under these conditions and related formation of ammonium carbamate ion pairs. Distinct evidence for a fourth and irreversibly formed moiety was observed under sorption of CO(2) under dry conditions. Significant amounts of physisorbed, linear CO(2) were detected at relatively high partial pressures of CO(2), such that they could adsorb only after the reactive amine groups were consumed.     

Effects of solute structure on local solvation and solvent interactions: results from UV/vis spectroscopy and molecular dynamics simulations

Solvation of heterocyclic amines in CO(2)-expanded methanol (MeOH) has been explored with UV/vis spectroscopy and molecular dynamics (MD) simulations. A synergistic study of experiments and simulations allows exploration of solute and solvent effects on solvation and the molecular interactions that affect absorption. MeOH-nitrogen hydrogen bonds hinder the n-pi* transition; however, CO(2) addition causes a blue shift relative to MeOH because of Lewis acid/base interactions with nitrogen. Effects of solute structure are considered, and very different absorption spectra are obtained as nitrogen positions change. MD simulations provide detailed solvent clustering behavior around the solute molecules and show that the local solvent environment and ultimately the spectra are sensitive to the solute structure. This work demonstrates the importance of atomic-level information in determining the structure-property relationships between solute structure, local salvation, and solvatochromism.     

Selective n,n-dibenzylation of primary aliphatic amines with dibenzyl carbonate in the presence of phosphonium salts

In the presence of catalytic amounts of tetraalkylphosphonium salts and under solventless conditions, primary aliphatic amines (RNH(2): R = PhCH(2), Ph(CH(2))(2), n-decyl, and 1-naphthylmethyl) are efficiently N-benzylated to the corresponding RN(CH(2)Ph)(2), using dibenzyl carbonate as the benzylating reagent. Compared to the reaction run without salt, where the competitive formation of the benzyl carbamate is favored, the phosphonium salt promotes high selectivity toward the benzylated amine and an increase of the reaction rate as well. However, in a single case explored for an amino acidic compound, namely 4-(aminomethyl)benzoic acid [4-(NH(2)CH(2))C(6)H(4)CO(2)H], both N,N-dibenzylation and esterification of the acid group were observed. Analysis of the IR vibrational modes of benzylamine in the presence of tetrabutylphosphonium bromide supports the hypothesis that this enhanced selectivity may be due to an acid-base interaction between the salt and the amine, which increases the steric bulk of the amine and favors attack of the nucleophile on the less hindered alkyl terminus of dibenzyl carbonate.     

Intramolecular trapping of an intermediate in the reduction of imines by a hydroxycyclopentadienyl ruthenium hydride: support for a concerted outer sphere mechanism

Reduction of imines by [2,5-Ph2-3,4-Tol(2)(eta(5)-C(4)COH)]Ru(CO)2H (1) produces kinetically stable ruthenium amine complexes. Reduction of an imine possessing an intramolecular amine was studied to distinguish between inner sphere and outer sphere mechanisms. 1,4-Bn(15)NH(c-C(6)H(10))=NBn (12) was reduced by 1 in toluene-d8 to give 85% of [2,5-Ph2-3,4-Tol(2)(eta(4)-C(4)CO)](CO)(2)RuNHBn(c-C(6)H(10))(15)NHBn (16-RuN,15N), resulting from coordination of the newly formed amine to the ruthenium center, and 15% of trapping product [2,5-Ph2-3,4-Tol(2)(eta(4)-C(4)CO)](CO)(2)Ru(15)NHBn(c-C(6)H(10))NHBn (16-Ru(15)N,N), resulting from coordination of the intramolecular trapping amine. These results provide support for an outer sphere transfer of hydrogen to the imine to generate a coordinatively unsaturated intermediate, which can be trapped by the intramolecular amine. An opposing mechanism, requiring coordination of the imine nitrogen to ruthenium prior to hydrogen transfer, cannot readily explain the observation of the trapping product 16-Ru(15)N,N.     

In silico screening of metal-organic frameworks in separation applications

Porous materials such as metal-organic frameworks (MOFs) and zeolitic imidazolate frameworks (ZIFs) offer considerable potential for separating a variety of mixtures such as those relevant for CO(2) capture (CO(2)/H(2), CO(2)/CH(4), CO(2)/N(2)), CH(4)/H(2), alkanes/alkenes, and hydrocarbon isomers. There are basically two different separation technologies that can be employed: (1) a pressure swing adsorption (PSA) unit with a fixed bed of adsorbent particles, and (2) a membrane device, wherein the mixture is allowed to permeate through a micro-porous crystalline layer. In view of the vast number of MOFs, and ZIFs that have been synthesized there is a need for a systematic screening of potential candidates for any given separation task. Also of importance is to investigate how MOFs and ZIFs stack up against the more traditional zeolites such as NaX and NaY with regard to their separation characteristics. This perspective highlights the potency of molecular simulations in determining the choice of the best MOF or ZIF for a given separation task. A variety of metrics that quantify the separation performance, such as adsorption selectivity, working capacity, diffusion selectivity, and membrane permeability, are determined from a combination of Configurational-Bias Monte Carlo (CBMC) and Molecular Dynamics (MD) simulations. The practical utility of the suggested screening methodology is demonstrated by comparison with available experimental data.     

Thermodynamics of aqueous potassium carbonate,piperazine, and carbon dioxide

CO₂ solubility was measured in a wetted-wall column in 0.6–3.6 molal (m) piperazine (PZ) and 2.5–6.2 m potassium ion (K⁺) at 40–110 °C. Piperazine speciation was determined using ¹H NMR for 0.6–3.6 m piperazine (PZ) and 3.6–6.2 m potassium ion (K⁺) at 25–70 °C. The capacity of CO₂ in solution increases as total solute concentration increases and compares favorably with estimates for 7 m (30 wt.%) monoethanolamine (MEA). The presence of potassium in solution increases the concentration of CO₃²⁻/HCO₃⁻ in solution, buffering the solution. The buffer reduces protonation of the free amine, but increases the amount of carbamate species. These competing effects yield a maximum fraction of reactive species at a potassium to piperazine ratio of 2:1.A rigorous thermodynamic model was developed, based on the electrolyte nonrandom two-liquid (ENRTL) theory, to describe the equilibrium behavior of the solvent. Modeling work established that the carbamate stability of piperazine and piperazine carbamate resembles primary amines and gives approximately equal values for the heats of reaction, ΔH_rxn (18.3 and 16.5 kJ/mol). The pK_a of piperazine carbamate is twice that of piperazine, but the ΔH_rxn values are equivalent (∼−45 kJ/mol). Overall, the heat of CO₂ absorption is lowered by the formation of significant quantities of HCO₃⁻ in the mixed solvent and strongly depends on the relative concentrations of K⁺ and PZ, ranging from −40 to −75 kJ/mol.     

Nanometer-scale ion aggregates in aqueous electrolyte solutions: guanidinium carbonate

Neutron diffraction with isotopic substitution (NDIS) experiments and molecular dynamics (MD) simulations have been used to characterize the structure of aqueous guanidinium carbonate (Gdm2CO3) solutions. The MD simulations found very strong hetero-ion pairing in Gdm2CO3 solution and were used to determine the best structural experiment to demonstrate this ion pairing. The NDIS experiments confirm the most significant feature of the MD simulation, which is the existence of strong hetero-ion pairing between the Gdm+ and CO3(2-) ions. The neutron structural data also support the most interesting feature of the MD simulation, that the hetero-ion pairing is sufficiently strong as to lead to nanometer-scale aggregation of the ions. The presence of such clustering on the nanometer length scale was then confirmed using small-angle neutron scattering experiments. Taken together, the experiment and simulation suggest a molecular-level explanation for the contrasting denaturant properties of guanidinium salts in solution.     

Insights into uranyl chemistry from molecular dynamics simulations

First-principles and purely classical molecular dynamics (MD) simulations for complexes of the uranyl ion (UO(2)(2+)) are reviewed. Validation of Car-Parrinello MD simulations for small uranyl complexes in aqueous solution is discussed. Special attention is called to the mechanism of ligand-exchange reactions at the uranyl centre and to effects of solvation and hydration on coordination and structural properties. Large-scale classical MD simulations are surveyed in the context of liquid-liquid extraction, with uranyl complexes ranging from simple hydrates to calixarenes, and nonaqueous phases from simple organic solvents and supercritical CO(2) to ionic liquids.     

铑催化合成气转化为乙醇反应中甲酰基中间体的化学捕获

本文采用化学捕获法对铑基催化剂上合成气转化反应中的甲酰基中间体进行了化学捕获,在CO+2D_2反应后,用CH_3I进行的化学捕获反应中生成了CH_3CHO、CH_3CDO两种形式的乙醛;补充的Ar吹扫实验显示DCO的甲基化反应对生成的CH_3CDO有重要贡献。因此,甲酰基的确是合成气反应中的C_1含氧中间体。根据这一结果,初步探讨了合成气反应中CH_x物种的生成途径。     

Molecular dynamics simulations for the growth of CH_4-CO_2 mixed hydrate

Molecular dynamics simulations are performed to study the growth mechanism of CH4-CO2 mixed hydrate in xCO2= 75%, xCO2= 50%, and xCO2= 25% systems at T = 250 K, 255 K and 260 K, respectively. Our simulation results show that the growth rate of CH4-CO2 mixed hydrate increases as the CO2 concentration in the initial solution phase increases and the temperature decreases. Via hydrate formation, the composition of CO2 in hydrate phase is higher than that in initial solution phase and the encaging capacity of CO2 in hydrates increases with the decrease in temperature. By analysis of the cage occupancy ratio of CH4 molecules and CO2 molecules in large cages to small cages, we find that CO2 molecules are preferably encaged into the large cages of the hydrate crystal as compared with CH4 molecules. Interestingly, CH4 molecules and CO2 molecules frequently replace with each other in some particular cage sites adjacent to hydrate/solution interface during the crystal growth process. These two species of guest molecules eventually act to stabilize the newly formed hydrates, with CO2 molecules occupying large cages and CH4 molecules occupying small cages in hydrate.     

A novel ionic liquids-based scrubbing process for efficient CO2 capture

A novel alkanolamine-based ionic liquid, N-methyl-diethanolammonium tetrafluoroborate ([MDEA][BF₄]), was synthesized in our laboratory. The ionic liquid-based composite solution consisting of N-methyl-diethanolamine (MDEA), [MDEA][BF₄], piperazine (PZ) and H₂O was investigated for CO₂ capture. The optimal performance for CO₂ capture was found at 45 °C, 1.50 MPa, probably due to a synergistic action of the reaction and the transport. No apparent corrosion was found on stainless and carbon steel with the above composite solution. This finding is very significant to the promotion of its engineering application.     

Molecular dynamics simulations for the growth of CH4-CO2 mixed hydrate

Molecular dynamics simulations are performed to study the growth mechanism of CH4-CO2 mixed hydrate in xco2 = 75%, xco2 = 50%, and zco2 = 25% systems at T = 250 K, 255 K and 260 K, respectively. Our simulation results show that the growth rate of CH4-CO2 mixed hydrate increases as the CO2 concentration in the initial solution phase increases and the temperature decreases. Via hydrate formation, the composition of CO2 in hydrate phase is higher than that in initial solution phase and the encaging capacity of CO2 in hydrates increases with the decrease in temperature. By analysis of the cage occupancy ratio of CH4 molecules and CO2 molecules in large cages to small cages, we find that CO2 molecules are preferably encaged into the large cages of the hydrate crystal as compared with CH4 molecules. Interestingly, CH4 molecules and CO2 molecules frequently replace with each other in some particular cage sites adjacent to hydrate/solution interface during the crystal growth process. These two species of guest molecules eventually act to stabilize the newly formed hydrates, with CO2 molecules occupying large cages and CH4 molecules occupying small cages in hydrate.     

Self-diffusion and transport diffusion of light gases in metal-organic framework materials assessed using molecular dynamics simulations

Metal-organic framework (MOF) materials pose an interesting alternative to more traditional nanoporous materials for a variety of separation processes. Separation processes involving nanoporous materials can be controlled by either adsorption equilibrium, diffusive transport rates, or a combination of these factors. Adsorption equilibrium has been studied for a variety of gases in MOFs, but almost nothing is currently known about molecular diffusion rates in MOFs. We have used equilibrium molecular dynamics (MD) to probe the self-diffusion and transport diffusion of a number of small gas species in several MOFs as a function of pore loading at room temperature. Specifically, we have studied Ar, CH4, CO2, N2, and H2 diffusion in MOF-5. The diffusion of Ar in MOF-2, MOF-3, and Cu-BTC has been assessed in a similar manner. Our results greatly expand the range of MOFs for which data describing molecular diffusion is available. We discuss the prospects for exploiting molecular transport properties in MOFs in practical separation processes and the future role of MD simulations in screening families of MOFs for these processes.     

Solvent effects on isomerization and spectral properties of photochromic-switching diarythene derivatives in polar and apolar solutions

The photocyclization behavior and dynamic conformational transition of photochromic switches of diarythene derivatives in solutions are investigated by using the density functional theory (DFT) and molecular dynamics (MD) simulations. Three possible conformations, antiparallel (anti), parallel (para), and twist, for the open-ring isomers of 1,2-bis(2-methylbenzothiophene-3-yl)maleic anhydride are located. Both PCM-B3LYP/6-31G* calculations and MD simulations demonstrate that anti and twist open-ring isomers can interconvert freely in n-hexane and acetonitrile solutions at room temperature. The statistical ratio of twist to anti isomers from MD simulations is 2.09 in n-hexane and 1.07 in CH(3)CN, in qualitative agreement with those (1.18 in n-hexane and 1.05 in CH(3)CN) estimated from Arrhenius analysis of DFT activation energies. The solvent polarity has little influence on the isomerization of open-ring isomers in the ground state. Due to the evident charge transfer upon excitations, the solvent effects on the electronic structures and absorption spectra of low-lying excited states (S(1) and S(2)) are more significant. For such charge-transfer excited states, the long-range corrected functional CAM-B3LYP gives better agreement with the experimental spectra than B3LYP. The solvent polarity and polarization of the charge-transfer excited states are crucial for fabricating the novel functionalized photochromic molecular switches.