News

Jul 21, 2023

Experimental and theoretical study of the effect of different functionalities of graphene oxide/polymer composites on selective CO2 capture

Scientific Reports volume 12, Article number: 15992 (2022) Cite this article There is a constant need for versatile technologies to reduce the continuously increasing concentration of CO2 in the

Scientific Reports volume 12, Article number: 15992 (2022) Cite this article

There is a constant need for versatile technologies to reduce the continuously increasing concentration of CO2 in the atmosphere, able to provide effective solutions under different conditions (temperature, pressure) and composition of the flue gas. In this work, a combination of graphene oxide (GO) and functionalized waterborne polymer particles was investigated, as versatile and promising candidates for CO2 capture application, with the aim to develop an easily scalable, inexpensive, and environmentally friendly CO2 capture technology. There are huge possibilities of different functional monomers that can be selected to functionalize the polymer particles and to provide CO2-philicity to the composite nanostructures. Density functional theory (DFT) was employed to gain a deeper understanding of the interactions of these complex composite materials with CO2 and N2 molecules, and to build a basis for efficient screening for functional monomers. Estimation of the binding energy between CO2 and a set of GO/polymer composites, comprising copolymers of methyl methacrylate, n-butyl acrylate, and different functional monomers, shows that it depends strongly on the polymer functionalities. In some cases, there is a lack of cooperative effect of GO. It is explained by a remarkably strong GO-polymer binding, which induced less effective CO2-polymer interactions. When compared with experimental results, in the cases when the nanocomposite structures presented similar textural properties, the same trends for selective CO2 capture over N2 were attained. Besides novel functional materials for CO2 capture and a deeper understanding of the interactions between CO2 molecules with various materials, this study additionally demonstrates that DFT calculations can be a shorter route toward the efficient selection of the best functionalization of the composite materials for selective CO2 capture.

Significant and continuous augmentation of the concentration of greenhouse gases in the atmosphere has turned into one of the most fundamental and persistent problems nowadays because fossil fuel reserves are still affordable and developing countries are in the process of economic growth. Among various greenhouse gases, carbon dioxide (CO2) is a key player in the global-warming scenario1. Even though the global CO2 capture capacity has reached 40 million tons by 2020, gigatonnes per year of CO2 have to be captured to have a significant impact on climate change2. Since green energy technologies are far away from the point to replace fossil fuel energy sources, the reduction of CO2 emissions and, thus, the reduction of global warming, is one of the most challenging environmental issues nowadays. Therefore, the development of efficient, selective, and low-cost carbon-capture technologies is crucial2. Strategies such as chemical/physical adsorption3, enzymatic conversion4, and membrane separation5 have emerged as potential solutions.

Various adsorbents have been proposed for CO2 capture6,7,8, namely, porous polymers, ion-exchange resins, covalent- and metal-organic frameworks, zeolites, silica- and alumina-based materials, metal oxides, etc. However, most of them suffer from low adsorption capacities (or a long time is required for saturation), lack of good chemical/thermal stability and/or selectivity relative to other gases, or they have reduced activity in the presence of moisture, such as zeolite-based adsorbents9, 10.

Carbon-based adsorbents are arising as a promising alternative to overcome most of the mentioned drawbacks, owing to one of the highest adsorption capacities and relatively low energy requirements for regeneration6,7,8. Furthermore, features such as large surface area, stability in cycle operations, porous structure that can be easily functionalized, and fast adsorption kinetics endorse them as one of the most promising adsorbents. Among these materials, due to lower production costs, graphene and its derivatives have been considered for commercial use11. Aimed to further improve the adsorption capacity and separation ability, functionalization of graphene surface with heteroatoms (N, S, O, etc.) has been widely investigated, as well as the production of composites with polymers such as polypyrrole12, polyaniline13, polyindole14, polythiophene15, mono-, di-, and triethylene-triamine16, tetraethylenepentamine17, poly(diallyldimethylammonium chloride)/polystyrene sulfonate18, poly(dimethylsiloxane)19, polyether block amide20, polyethylene-imine21, and also with metal–organic frameworks22.

The addition of polymer nanoparticles onto graphene platelets, together with the functionalization of the surface, improves the physicochemical properties and provides easier handling. Such composite platelets are more durable and present improved stability in cycle operations. The polymer particles are usually produced by emulsion polymerization of different combinations of (meth)acrylic monomers, functionalized onto the surface by using a minor amount of functional monomers during the synthesis. The selection of highly CO2-philic functional monomers may be decisive towards the development of efficient composite platelets for CO2 selective adsorption and, thus, the screening and evaluation of different monomers is a fundamental step23, 24.

Computational chemistry has been revealed as a powerful tool in materials science25, 26. In particular, density functional theory (DFT) calculations may provide a detailed insight into the interactions between CO2 and different molecules, giving useful information for the pursuit of potential new adsorbents. The first computational studies, regarding the interaction of CO2 and graphene, were performed even before the successful isolation of graphene. Cinke et al.27 used benzene and coronene as a model of graphene to estimate its interaction with CO2. Allouche and Ferro28 performed mixed quantum mechanics/molecular mechanics ONIOM (B3LYP:UFF) calculations on a moderate-sized cluster model of graphene. According to their results, CO2 can not be adsorbed. Xu et al.29 did more precise ONIOM calculations (B3LYP:DFTB-D) on a smaller cluster. They obtained small interaction only after correction on zero-point vibrational energies. Cabrera-Sanfelix30 found that CO2 can be better physisorbed on top of vacancies in the graphene structure. Liu and Wilcox31 performed more detailed calculations and stated that the CO2 binding energy on a monovacancy is 4-fold higher than that on the defectless graphene. The same authors also showed that hydroxyl and carbonyl groups in graphene establish stronger interactions with CO2 due to higher electron densities32. Sun et al.33 showed improved adsorption near nitrogen atoms in N-doped graphene and, by means of energy decomposition analyses, demonstrated the relevance of dispersion interactions in this process. Wang et al.34 performed a similar study on N and O co-doped graphene and established that the heteroatoms bind CO2 more strongly. They also highlighted the relevant role of dispersion interactions. Dasgupta et al.35 used higher levels of theory and small molecular models to estimate the binding energy between graphene functionalized with different groups and CO2. Seema et al.15 synthesized S-doped microporous carbon materials by chemical activation of reduced graphene oxide/polythiophene (rGO/PTh) and showed, using DFT calculations, that the interaction energy of CO2 with thiophene is higher than that with pyrrole, in agreement with the observed higher adsorption capacity of rGO/PTh. These authors suggested that the reason for this stronger attraction is due to a larger negative charge of S in thiophene than that of N in pyrrole. Finally, polymeric systems have also been considered. Patel et al.36 calculated the binding energies for CO2 adsorption on several azo-bridged covalent organic polymers. Although simple models were used, accurate trends in binding energy were obtained. In summary, quantum chemical calculations confirm the well-known fact that carbon materials doped with heteroatoms show enhanced CO2 adsorption. Das et al.37 investigated the role of surface OH groups in triazine-based N-rich porous organic polymers for enhancing CO2 capture. Using the ωB97XD density functional method, capable of a reliable description of both short-range and long-range interactions, with, a relatively small 3-21G* basis set, they found that CO2 mainly interacts through the N–N⋯O and O–H⋯O hydrogen bonds, but also through the weak N⋯C, N–O⋯C, and C–H⋯O interactions. Ullah et al.38 performed a combined experimental and theoretical study of CO2 adsorption by amine and amide porous polymers. Using the B3LYP-D3/6-311++G** level of theory they obtained a moderate agreement with experiments due to the different morphologies of samples.

As has been previously mentioned, GO/polymer composites are very promising candidates for commercial applications, and the selection of optimal functional monomer can significantly increase adsorption capacity. However, comprehensive investigation, both experimental and theoretical, of the effects that different functionalities of GO/polymer composites have on selective CO2 capture are sparse. Even more, to the best of our knowledge, there is no data in the literature on the theoretical estimation of binding energy for the CO2 adsorption on these composites. Thus, in this work, we have performed quantum chemical calculations to estimate the binding energy of CO2 with a set of GO/polymer composites. More precisely, in the pursuit of new materials for CO2 capture, we investigated copolymers of methyl methacrylate (MMA), n-butyl acrylate (BA), and different functional monomers. Furthermore, the functional monomers, which, according to the theoretical study provided the strongest interactions with CO2, have been used for the synthesis of functionalized MMA/BA particles and combined with graphene oxide platelets. The synthesis of composites is specially designed to give a reliable platform for comparison with the theoretical results. Finally, the synthesized composites have been evaluated experimentally for CO2 capture. Besides, a similar study has been performed for N2 to evaluate the CO2/N2 selectivity of these materials. The observed experimental trends of CO2 adsorption capacities are in accordance with the trends of computed binding energies, demonstrating that the DFT calculations are a useful tool for the development of new materials for CO2 capture. Additionally, it is convenient to say that our procedure for synthesis has advantages over many similar reports since it is done by low energy synthesis procedure (using ascorbic acid as redox initiators under mild conditions, i.e. without any solvents, high temperature, and C-footprint treatments, nor the use of aggressive acid–base treatments). Having that in mind, it is aligned with the green chemistry trends and has the prospect to scale up well.

All geometry optimizations were performed within density functional theory (DFT) by using the long-range corrected ωB97XD functional39 together with the 6–31+G(d) basis set. For the calculations including the GO model, the smaller 6–31G(d) basis set was used to relieve the computational effort. The zero-point vibrational energies (ZPVE) were evaluated within the harmonic oscillator approximation at the same level of theory. All optimized structures showed no imaginary frequencies, meaning that all structures are minima on the potential energy surface. The CO2 and N2 molecules were placed in different positions near every functional group of adsorbents and the structure with the lowest energy (highest binding energy) is presented. Also, when the structure of composites was optimized, the polymer was placed so that as many interactions between functional groups of GO and polymer are achieved. No less than six rearrangements were checked. The electronic energies were then refined by single-point calculations performed on the optimized structures using the 6–311++G(2df,2p) basis set. All calculations were carried out with the Gaussian 16 package40.

As a compromise between computational cost and accuracy, the model of GO is constructed to be the one that can guarantee that the polymer remains within the line which connects centroids of outer benzene rings. Graphene is randomly functionalized with two carboxyls, two epoxies, and one hydroxyl group in order to approximate the average functionalization of the real GO, that is, to achieve an adequate ratio between carbon and other atoms, as well as the uniformity of the distribution. Graphene edges were fully terminated by hydrogen atoms. Lengths of the same types of bonds in the interior of the graphene plane and near the edges were similar, indicating that edge effects are not significantly pronounced. In particular, the molecular formula for the model of the GO platelet was C69H22O7 (Supplementary Fig. S1). This is the smallest molecular model that would allow including the main interactions between the GO platelet and both the functional monomer and CO2/N2.

The pursuit of the global energy minimum for the GO/monomer/CO2 and GO/monomer/N2 clusters would demand the exploration of a huge conformational space, which is out of our present computational capabilities. Thus, in order to find the most stable configuration, only several conformations were considered, where the CO2 and N2 molecules were located in different positions, near every functional group. Besides, the functional monomer was arranged in such a way that as many interactions as possible were established with the GO platelet. The energy differences between all the conformations explored were very small, and the structure with the lowest energy is used for the discussions.

An aqueous dispersion of graphene oxide sheets of 4 mg mL−1 (Graphenea) was used as supplied. The monolayer content in the dispersion was > 95% and in a pH range between 2.2 and 2.5. The elemental analysis of graphene oxide layers was provided in the technical data sheet from Graphenea: C (49–56%), H (0–1%), N (0–1%), S (2–4%), and O (41–50%). Technical monomers, methyl methacrylate (MMA, Quimidroga) and butyl acrylate (BA, Quimidroga), were used as supplied without any further purification. Sodium 4-vinylbenzenesulfonate (NaSS, Sigma-Aldrich), glycidyl methacrylate (GMA, Acros Organics), 2-hydroxyethyl methacrylate (HEMA, Sigma-Aldrich), and 2-aminoethyl methacrylate hydrochloride (AEMH, Sigma-Aldrich) were used as functional monomers. Tert-butyl hydroperoxide solution (TBHP, Sigma-Aldrich) and l-ascorbic acid (AsA, Sigma-Aldrich) were employed as redox initiators. Furthermore, sodium dodecyl sulfate (SDS, Sigma-Aldrich) and hexadecyltrimethyl ammonium chloride (HAC, Sigma-Aldrich) were employed as emulsifiers. Sodium bicarbonate (NaHCO3, Sigma-Aldrich) was used as a buffer. Deionized water was used throughout the experimental work.

The batch emulsion polymerization process was used for the synthesis of functionalized polymer nanoparticles in aqueous dispersion (polymer latex). As main monomers, MMA and BA in 50/50 weight ratio were used, to which 3 wt% of functional monomer was added, i.e. NaSS, GMA, HEMA, and AEMH. In all cases, the same formulation was employed as described in Table 1, with a solids content of the final aqueous dispersions of 20 wt%. Two different types of surfactants were used, SDS in GMA, HEMA, and NaSS systems and HAC in AEMH system.

The reactions were performed in a 250 mL jacketed glass reactor equipped with a reflux condenser, temperature probe, nitrogen and feeding inlet, and stainless-steel stirrer rotating at 200 rpm. Pre-emulsion (monomers, surfactant, buffer, and water) was charged in the reactor, and then the temperature was raised to 70 °C, after which an aqueous solution of TBHP initiator was added as a shot, whereas AsA was fed for 120 min into the reactor. The reaction mixture was then kept at 70 °C for an additional 30 min before cooling to room temperature.

Composite platelets were prepared by simple blending of aqueous dispersions of GO and polymer particles (1:0.5 weight ratio of GO:polymer), which were left agitated overnight at room temperature. Afterward, the hybrid dispersion was subjected to a freeze-drying process, in Telstar LyoQuest 55 at − 49 °C and 0.2 mbar for 3 days. To check how much polymer was incorporated onto the GO platelets, the aqueous phase of dispersion was analyzed gravimetrically. The agitation of GO and polymer mixture was stopped at different time periods. When the agitation was stopped, the platelets precipitated and the residual aqueous phase was analyzed. 2 mL of the dispersion extracted at different times was dried in the oven. The polymer quantity present in the aqueous phase was calculated from the difference in weight of wet samples. It was found that after 3 h agitation of the dispersion mixture of GO and polymer, there was no polymer left in the dispersion, or in other words, all added polymer was completely incorporated onto the GO platelets.

This method of synthesis was selected in order to prevent the development of complex, hierarchical porous morphology, typical for these systems, as we reported previously23, 24, expecting to produce composite platelets that differ just in the functionalization, introduced by the selected functional monomers. In this way, a solid platform for comparison with the theoretical models used for the calculations was provided.

The z-average particle size (dz) of the polymer particles was measured by Dynamic Light Scattering Spectroscopy (DLS), using the Malvern Zetasizer Nano ZS. Before the measurement, a fraction of the latex was sufficiently diluted with deionized water in order to avoid multiple scattering. Analyses were carried out at a temperature of 25 °C. The reported particle sizes are the average of three repeated measurements per sample.

The gel contents of the polymer particles (fraction of polymer insoluble in THF due to the presence of crosslinked and branched polymer chains) were determined by Soxhlet extraction. A few drops of latex were placed on glass fiber square pads and dried overnight at 60 °C. Then, the filter, together with the dried polymer, was subjected to a continuous extraction with THF under reflux in the Soxhlet for 24 h. After that, the filter was dried overnight at 60 °C in order to weigh the gel polymer fraction.

The molecular weights corresponding to the soluble fraction of the polymers were determined by gel permeation chromatography (GPC). The soluble part from the Soxhlet extraction was first dried, redissolved in GPC grade THF at a concentration of 2 mg mL−1, and finally, the solution was filtered (polyamide Ф = 45 μm) before injection into the GPC instrument via an autosampler (Waters 717). The GPC consisted of a pump (LC-20A, Shimadzu), a differential refractometer (Waters 2410), and three columns in series (Styragel HR2, HR4, and HR6, with pores sizes ranging from 102 to 106 Å). The chromatograms were obtained at 35 °C using a THF flow rate of 1 mL min−1. The equipment was calibrated using narrow polystyrene standards, thus, the presented molecular weights are relative to this standard.

In terms of graphene-polymer composite materials, thermal stability and the amount of oxygen-containing functional groups presented in the monolithic structure were studied by the TGA500 apparatus (TA instrument). Samples of around 2 mg were heated under a nitrogen atmosphere (90 mL min−1) from 25 to 800 °C, at a rate of 10 °C min−1.

The surface morphology of the composites was analyzed by scanning electron microscopy (SEM) using Hitachi TM3030 scanning electron microscope at 15 kV after the samples were coated with a thin gold layer, whereas the structure of the composites was analyzed by transmission electron microscopy (TEM) using Tecnai TM G2 20 Twin device at 200 kV (FEI Electron Microscopes). Prior to analysis, the samples were embedded in epoxy resin, from which ultra-thin sections (80 nm) were cut with a diamond knife on Leica EMFC6 ultramicrotome device and placed on a 200 mesh copper grid.

The textural properties of the monoliths were characterized by N2 adsorption–desorption, performed at − 196 °C in a Micromeritics ASAP 2020. Prior to the measurements, the samples were degassed at 110 °C for 8 h under vacuum. From N2 adsorption–desorption isotherms, the specific surface area and the adsorption average pore width (4 V/A) were calculated from the Brunauer–Emmett–Teller (BET) equation. Moreover, the t-plot method was used to estimate the micropore volume. Finally, the pore volume was calculated using the method proposed by Barrett–Joyner–Halenda. The CO2 adsorption capacities of the 2D graphene-polymer composites were determined using a TGA analyzer. Prior to the adsorption measurements, the samples were heated to 100 °C in an N2 atmosphere at a flow of 50 mL min−1 and held at that temperature for 30 min. The samples were then allowed to cool to 25 °C. Once the temperature reached 25 °C, the gas was changed to pure CO2 at a flow rate of 50 mL·min−1 until a constant weight was observed. The weight change of the sample was interpreted as the CO2 adsorption capacity. The same measurements were performed for N2 adsorption. These values were corrected by taking into account the buoyancy effect of gas change during the measurements.

In this work, a set of functional monomers, namely, acrylamide (Am), 2-aminoethyl methacrylate hydrochloride (AEMH), 4-bromostyrene (BS), 2-hydroxyethyl methacrylate (HEMA), glycidyl methacrylate (GMA), methyl 2-chloro acrylate (MClA), and sodium 4-vinylbenzenesulfonate (NaSS), see Fig. 1, have been computationally studied in the pursuit of new materials to improve the CO2 adsorption capacity of GO/polymer composites. To this end, the binding energy of CO2 in these materials has been estimated using the theoretical methods of quantum chemistry and the most promising composites have been evaluated experimentally.

Molecular structures of the monomers: acrylamide (Am), 2-aminoethyl methacrylate hydrochloride (AEMH), methyl 2-chloro acrylate (MClA), hydroxyethyl methacrylate (HEMA), glycidyl methacrylate (GMA), sodium 4-vinylbenzenesulfonate (NaSS), 4-bromostyrene (BS), butyl acrylate (BA), and methyl methacrylate (MMA).

The results are organized as follows: (1) the computational studies comprise the analysis of the binding of CO2 with the bare functional monomers, then with MMA/BA/monomer copolymers, and finally with GO/copolymer composites. These composites have also been used to analyze the binding energy for N2. (2) The experimentally evaluated CO2 and N2 capture abilities of selected composites.

The theoretical investigation was divided into four parts. In the first part, we calculated the binding energies between a set of functional monomers and CO2. In the second one, it was examined how the introduction of MMA and BA monomers affects the values of binding energy and the trends between them. In the third part, for some of the functional monomers, we estimated the binding energy between CO2 and GO/MMA/BA/monomer composites, as well as the interaction between the GO platelet and the copolymer. Finally, in the last part, the previously selected composites were used to study the binding energy with N2.

In this subsection, the CO2 affinity of the functional monomers represented in Fig. 1 was analyzed in terms of binding energies and geometrical features. In Table 2 are collected binding energies and selected geometrical parameters of the complexes formed by CO2 and the functional monomers. The optimized geometries are represented in Fig. 2. The binding energy was estimated, using the ZPVE-corrected energies (E0), as:

Optimized geometries of the CO2-monomer complexes.

These energies include the correction of the basis set superposition error (BSSE) by means of the counterpoise method41, 42.

Inspecting Table 2, it is observed that the weakest interactions with CO2 correspond to BS and MClA monomers (− 5.5 kJ·mol−1, for both of them). As a consequence, the C=O bond lengths of CO2 remain almost unchanged compared to the free molecule (1.165 Å) and the O=C=O angle deviates slightly from linearity (179.0° in both cases). Also, the CO2-monomer distances (D) are longer than for other monomers.

The highest binding energy corresponds to NaSS (− 29.1 kJ·mol−1), which is reflected in a larger deformation of CO2 molecule (R1 = 1.174 Å, R2 = 1.156 Å and α = 176.2°). Also, the distance between S=O groups and CO2 is relatively small considering the presence of the sodium cation. Furthermore, the relative position of this cation with respect to the S=O groups is modified after the introduction of CO2.

The remaining four functional monomers show intermediate binding energies. Again, for larger binding energies, larger deformations of CO2 and, in general, smaller distances between this molecule and the monomer are observed. The only clear exception is AEMH, which shows one of the lowest binding energies (− 7.2 kJ·mol−1), but a rather small CO2-monomer distance (D = 2.934 Å) and remarkable changes in CO2 geometry (R1 = 1.172 Å, R2 = 1.158 Å and α = 176.0°). This monomer comprises a molecular complex between the HCl and the NH2 moiety. After the adsorption of CO2, the molecular geometry of AEMH is significantly changed; the NH2-HCl interaction is remarkably affected, in such a way that the acidic proton is almost transferred to the amino group (the N–H distance is shortened from 1.612 to 1.460 Å, while the H–Cl distance is enlarged from 1.377 to 1.717 Å). Thus, the capacity of the nitrogen atom to interact with CO2 is decreased, which is reflected in the low binding energy. Finally, AEMH and NaSS are the only functional monomers for which geometry is changed with the adsorption of CO2, due to the presence of weakly bonded molecules or ions.

In this subsection, the CO2 affinity of copolymers composed of methyl methacrylate (MMA), butyl acrylate (BA), and the functional monomers are evaluated. The introduction of MMA/BA in the system notably increases the binding energy for CO2 in all cases (between 4 and 12 kJ·mol−1, depending on the functional monomer), except for Am, see Table 3. This may be explained by the presence of the carbonyl groups of the MMA/BA copolymer, which allows additional interactions with CO2. The molecular structures of the copolymer-CO2 complexes are in the Supplementary Information (Supplementary Figs. S2–S8). It is convenient to say that values of calculated energies are in the range of those which can be found in the literature for CO2 adsorption by polymers7, 36, 37, which indicates that the size of the model and level of theory are adequate.

The decrease observed in Am-containing copolymer with respect to the free Am (from − 16.1 to − 13.4 kJ·mol−1) can be ascribed to the interaction of the NH2 of Am with the C=O of BA (see Supplementary Fig. S2). Now, the amino group is less available to interact with CO2 and, consequently, CO2 is displaced closer to the carbonyl group of Am. For AEMH, the presence of MMA/BA allows CO2 to locate in a position closer to the amino group, and several bonds in AEMH change significantly (see Supplementary Fig. S3). For BS, the introduction of MMA and BA also changes the positioning of CO2, which is now located in a parallel plane over the benzene ring of BS, allowing for π–π interactions that increase the binding energy (see Supplementary Fig. S4). In the case of MClA, CO2 is now located relatively far from the Cl atom and closer to the carbonyl group, enhancing notably the affinity (see Supplementary Fig. S5). For the GMA system, CO2 is located between the C=O group of BA and the epoxide group of GMA (see Supplementary Fig. S6). The CO2-BA distance is smaller than the CO2-GMA one in the free monomer. As a consequence, the binding energy is increased. Similar to the previously discussed cases, in copolymer with HEMA, CO2 is adsorbed both by BA and HEMA. In the free monomer, CO2 mainly interacts with the carbonyl group; however, after including the MMA and BA, CO2 is located close to both the OH group of HEMA and the carbonyl group of BA (see Supplementary Fig. S7a) and, thus, the largest increment of the binding energy (11.5 kJ·mol−1) is observed. The introduction of CO2 notably alters the molecular structure of the copolymer (see Supplementary Fig. S7b). Finally, in the copolymer including NaSS, the presence of MMA/BA allows the CO2 to locate in a position where π–π interactions with the benzene ring of functional monomer (i.e. NaSS) are favored, as in the copolymer with BS (see Supplementary Fig. S8).

In this section, the performance of the GO/copolymer composites in CO2 capture is analyzed for the four copolymers that showed the highest affinity with CO2, concretely, those which include AEMH, GMA, HEMA, and NaSS functional monomers. In this case, the interaction energy between the GO platelet and the copolymer was estimated, as in the previous section, by using the ZPVE-corrected energies (E0):

The results were collected in Table 4. As in the case of the values from Table 3, calculated energies are similar to those measured and calculated for GO and its composite with polymers15, 33, 34, which indicates the reliability of the results. It is remarkable how the binding energy is increased for the composite including GMA (from − 18.3 to − 54.0 kJ·mol−1), while a softer increase is found for AEMH (from − 14.7 to − 19.3 kJ·mol−1). However, the binding energy is decreased for the composites including NaSS (from − 33.4 to − 23.1 kJ·mol−1) and HEMA (from − 24.2 to − 17.0 kJ mol−1). These last two cases are striking since it is well known that GO is a good CO2 adsorbent11 and, therefore, a cooperative effect is expected, enhancing the capacity of the copolymer.

The composite with GMA presents the highest binding affinity for CO2 (− 54.0 kJ·mol−1). Without CO2, the copolymer is oriented in such a way that the carbonyl group of MMA is close to one of the carboxyl groups of GO. The interaction energy of GO and the copolymer is − 137.7 kJ·mol−1. However, the introduction of CO2 significantly changes the positioning of the copolymer and, thus, after the adsorption, the epoxy group of GMA is now able to interact with the acid group of GO, with an Oepox–HCOOH distance of only 1.721 Å. The CO2 is placed between the carboxyl group of GO and the carbonyl group of MMA and, thus, a cooperative effect between GO and the copolymer is observed (see Supplementary Fig. S9).

When the AEMH-containing copolymer is placed in the GO platelet, the binding energy is increased to a lower extent than in the case of GMA. AEMH is placed in such a way that the Cl‒ ion interacts with the OH group of the GO surface, while the amino group is close to the carboxyl group. This makes the copolymer-GO interaction energy to be notably higher (− 172.2 kJ·mol−1) than for the GMA copolymer. The molecule of CO2 interacts mainly with the amino group of AEMH, as in the copolymer without GO, weakening the GO-copolymer interaction (the Cl−⋯HO distance increases from 2.259 to 2.284 Å). Thus, in the presence of CO2, the copolymer is attached to GO mainly by the interaction between the carbonyl group of AEMH and the nearby carboxyl group. The small increase observed in the binding energy can be due to a rather small model of GO. In Supplementary Fig. S10 can be observed how the copolymer is covering the whole GO platelet and, therefore, the CO2 is able to interact mainly only with the copolymer, in such a way that the effect of GO is almost absent. A larger model probably would increase this binding energy.

In the case of NaSS, there is a strong interaction between the copolymer and GO. Na+ is placed near the epoxy group and one of the oxygens from the sulfonate group interacts with the carboxyl group of GO. Also, BA moiety is oriented so that the carbonyl group is in the vicinity of another epoxy group of the GO surface. The distance between Na+ and the epoxy group is increased after CO2 adsorption, i.e. Na+ moves towards the carboxyl group. Besides, the sulfonate anion moves further from the COOH and C=O groups of GO. In the composite, CO2 is placed closer to the sulfonate group, while in the isolated copolymer there is an interaction with the carbonyl of BA and the aromatic ring of NaSS. Thus, the increase in the interaction energy provided by the GO platelet is compensated by a weaker interaction with the copolymer, and the net effect is a small decrease in binding energy (see Supplementary Fig. S11).

Finally, the HEMA composite presents the lowest binding affinity toward CO2. Also, of all four investigated composites, the one with HEMA has the lowest interaction energy between the GO platelet and the copolymer. The copolymer is placed in such a way that OH groups of HEMA and GO form hydrogen bonds, while oxygen atom from ether moiety of the carboxyl group of MMA forms another with the carboxyl group of GO. CO2 is positioned close to the hydroxyl group of HEMA and in relatively close vicinity of the carbonyl group of BA (see Supplementary Fig. S12). After the adsorption, the copolymer moves further from GO. For instance, two mentioned hydrogen bonds change the length from 2.050 to 2.124 Å and from 1.831 to 1.843 Å, respectively. As in the case of NaSS, CO2 is positioned far from the GO and interacts mainly with the copolymer. Therefore, as a result, binding energy decreases.

In this subsection CO2/N2, selectivity of investigated materials was estimated. More precisely, energies by which four composites selected in the previous subsection bind N2 were calculated. The results were collected in Table 5. The composite with AEMH binds N2 with the lowest energy. The nitrogen molecule is positioned close to the butyl group of BA and the carbonyl group of MMA (see Supplementary Fig. S13). Similarly, in the case of NaSS, N2 is placed between the butyl group of BA and one of the oxygens from the sulfonate group, to which it is close (see Supplementary Fig. S14). Therefore, the binding energy is only slightly higher than in the previously discussed composite. In the composite with HEMA, bonding occurs through the ether group of HEMA and carbonyl group of BA (see Supplementary Fig. S15), and thus, since N2 is placed close to the two heteroatoms (i.e. oxygens), the energy of adsorption is higher than in the case of the AEMH and the NaSS. Lastly, the composite with GMA has the highest binding energy, and nitrogen is placed close to the oxygens of the ether group of GMA and carbonyl group of BA (see Supplementary Fig. S16). For any of the composites, the introduction of N2 does not remarkably change the relative position of the copolymer in the composite. The position at which CO2 is primarily bound differs from that where N2 is adsorbed. Moreover, in the case of GMA and AEMH, these two molecules bind CO2 at significantly different positions.

For the experimental work, GMA, AEMH, HEMA, and NaSS were selected as functional monomers added in a small amount (3%) to the main monomer mixture made of MMA/BA, to produce functionalized polymer nanoparticles. Table 6 presents the characteristics of the polymer particles. Despite that the syntheses were performed under identical conditions, the polymer dispersions have distinct features, which means that the type of functional monomer affected the polymerization process and the polymer properties. NaSS-, GMA- and HEMA-functionalized particles have a similar average particle size of about 80 nm, whereas AEMH functionalized particles are much larger and polydispersed. Obviously, the last system was less colloidally stable during emulsion polymerization, likely due to the cationic surfactant employed for colloidal stabilization in order to avoid possible ionic interactions between the surfactant and AEMH. As a consequence of the larger average size of AEMH functionalized particles, shorter polymer chains were produced due to the lower number of radicals per particle, reducing the possibility of bimolecular termination. The THF insoluble polymer fraction (gel fraction) that indicates the presence of cross-linked and branched structures was very low. The gel fraction of NaSS-containing polymer was slightly increased, which is probably the result of the incorporation of sulfonate moieties within MMA/BA chains, lowering their solubility in THF, especially high molar mass chains. GMA- and HEMA-containing polymers presented an increased amount of insoluble fraction, likely corresponding to crosslinked and branched structures.

The nanocomposite structures were created by mixing the functionalized polymer particles dispersion with GO nanoplatelets dispersion, during which process the polymer particles were adsorbed onto the platelet surface. By gravimetrical analysis of resultant water, after the composite platelets were recuperated, it was shown that the whole amount of polymer from the dispersion was incorporated within the composite structure. Therefore, the resulting composites contain GO and polymer in a 1:0.5 weight ratio. In Fig. 3, the morphology of the nanocomposites determined by SEM is presented under two magnifications. Porous materials were obtained in all cases, probably due to the hydrophobic interaction of the composite rGO-polymer platelets placed in aqueous dispersion. Namely polymer nanoparticles interacted with the oxygen functional groups at the GO surface, creating H bondings (H donor groups such as COOH and OH in GO interacted with H accepting groups in polymers, such as carbonyl, sulfonate, and amino). As the oxygen functional groups on GO provide amphiphilic character to GO, by decreasing the number of these functionalities, GO platelets became more hydrophobic. This induced aggregation, wrinkling, and crumpling of the platelets, forming the porous structure, which processes were even promoted during drying. However, the presence of different functional monomers, even in as low amounts as 3 wt% with respect to polymer, influences the structure and morphology. Placed onto the polymer particles surface, the different functionalities affected polymer—GO interactions. While nanocomposites functionalized with NaSS, GMA, and HEMA (Fig. 3a,c,g, respectively) have fluffy structures with very well-developed pores of about 5–10 μm, the nanocomposite functionalized with AEHM is less porous (Fig. 3e). The cationic nature of the last probably induced ionic complexing with the numerous anionic oxygen-containing functional groups of GO, giving rise to more compact composite structures. Under higher magnification (Fig. 3b,d,f,h), no important differences between the four nanocomposites may be noticed, which creates a stable platform for comparison of the interaction forces between CO2 and the respective functionalities in the nanocomposites.

Morphology of the GO/polymer composites with different functional monomers under different magnification: (a) and (b) NaSS; (c) and (d) GMA; (e) and (f) AEMH; (g) and (h) HEMA.

The structure of the nanocomposites was evaluated by TEM imaging (Fig. 4), where the cross-section of the composite platelets may be observed. The black areas represent the GO platelets, whereas the white areas correspond to the polymer. A peculiar combination of these two phases may be observed, in which the platelets wrap the single polymer particle or aggregates of a few of them, creating composite honeycomb-like structures. The presence of GO platelets probably prevented the full particle coalescence and formation of large polymer areas. The thickness of the composite platelets depends on the functional monomer, therefore, the composite platelets functionalized with NaSS present a thickness of about 200 nm, those functionalized with GMA and HEMA have a thickness of 250–300 nm, whereas the AEMH functionalized platelets, with a thickness of 500–1000 nm, are the thickest. Two possible causes can affect the formation of thicker AEMH-containing composite platelets. On one hand, the size of AEMH functionalized polymer particles is double on average in comparison to other functionalized particles, and on the other, the ionic interactions between cationic polymer particles and anionic GO increased the likelihood of platelets aggregation.

TEM images of (a) NaSS; (b) GMA; (c) AEMH; and (d) HEMA functionalized composites. Scale bars in all images are 500 nm.

TGA curves of the nanocomposites, shown in Supplementary Fig. S17, were used to determine the content of oxygen functionalities in each composite. Thermal degradation occurred in three steps, assigned as follows. The humidity is lost until 100 °C, the weight loss between 100 and 260 °C corresponds to a loss of oxygen functionalities distributed onto GO, whereas the polymer is degraded between 300 and 400 °C23. The advantage that the oxygen functionalities over GO are lost in the distinct region than the polymer itself (including the functionalities containing oxygen within the polymer chains), provided the possibility to calculate their relative contents. The content of oxygen functionalities (originating from GO) within the composites and their textural properties are presented In Table 7. The quantity of the oxygen functional groups is similar in all composites.

Relatively modest specific BET surface area was observed for all nanocomposites (Table 7), which is not surprising, as the synthesis of the nanocomposites was altered in order to limit the development of the porous structures and to provide a base to investigate the effect of functionalities on the CO2 selective capture. The observed porosity in SEM images is surface morphology on the micron level, while the textural properties from Table 7 demonstrate that no deep meso- and micropores were developed. NaSS functionalized composite presents the lowest BET surface area and total pore volume, indicating that this composite is less porous and more compact than others, which might be due to the aromatic ring from NaSS functional monomer that interacts more tightly with GO than other polymers. Despite this, only NaSS functionalized composite is characterized by microporosity, although in small quantity. It might provide compensation in the case of CO2 adsorption, as it is known that microporosity is one of determining characteristics of the CO2 adsorption capacity23. The textural properties of GMA and AEMH composites are similar.

The CO2 and N2 adsorption capacity of nanocomposites is given in Supplementary Table S1 and in Fig. 5. For comparison, the CO2 adsorption of neat polymers is also shown. Taking into account that the gas adsorption was studied at atmospheric conditions, CO2 adsorption is within the range typical for carbon-based nanomaterials (0.5–1 mmol g−1), but it is still much lower than those we previously achieved by similar materials (3.7 mmol g−1)24. Surprisingly, contrary to the previous study24, N2 adsorption was very similar to that of CO2 adsorption, indicating that nanocomposites have very low selectivity. The modest capacity for selective CO2 capture of the present nanocomposites might be due to a few reasons. The most important for this study is that the composite synthesis was directed to obtain possible similar materials and thus to provide a good basis for the comparison using theoretical studies.

Adsorption of CO2 and N2 of the functionalized GO/polymer composites.

The amount of CO2 adsorbed by neat polymer materials is approximately one order of magnitude lower than that adsorbed by nanocomposites. (Supplementary Table S1). Considering that the polymers are not porous, this result is not surprising. Nevertheless, it is clear that a combination of polymers with GO might be a useful way to increase their CO2 adsorption performance.

The comparison of the experimental results with the theoretical prediction of the binding energy between nanocomposites and CO2 is not straightforward. Even with the intention to eliminate the effect of the textural characteristics of the synthesized nanocomposites, only GMA and AEMH functionalized nanocomposites presented similar porous morphology, providing a basis for comparison. According to Table 4, GMA functionalized composite presents the highest binding energy, much higher than that of AEMH (54 kJ mol−1 versus 19 kJ mol−1). According to Supplementary Table S1, GMA composite adsorbs 0.7 mmol g−1, whereas AEMH composite adsorbs 0.5 mmol g−1, thus, the same trend is followed. On the other hand, Table 5 presents that GMA-based composite has N2 binding energy of almost 10 kJ·mol−1, whereas AEMH composite only 7 kJ·mol−1. The same trend is followed by experimental results, thus GMA composites adsorbed 0.61 mmol g−1 of N2, and AEMH adsorbed 0.34 mmol·g−1. Therefore if, by appropriate design of synthesis procedure, the effect of morphology over the gas absorption capacity of the material can be eliminated, a comparison with the theoretical prediction of the adsorption can be performed. In such a case, the presented theoretical study seems to be an excellent tool to predict the interaction of functionalized composite structures with CO2 and N2, which can be useful for the selection of functionalization of composites for application in gas adsorption.

Moreover, the comparison of the theoretical prediction and experimental results for HEMA composites provides evidence of the importance of the porous structures’ morphology for the selective CO2 capture. According to Table 4, HEMA functionalized composite presented three-time lower binding energy for CO2 and slightly lower N2 binding energy, as shown in Table 5. Experimental results show that this composite presents the highest CO2 adsorption (1 mmol·g−1) and the highest selectivity over N2 (Fig. 5). Considering the textural properties of HEMA composites presented in Table 7, i.e. the highest BET surface area and total pore volume, this composite is clearly the most porous when compared to the other studied nanocomposites in this work, and accordingly, it achieved the highest performance for selective CO2 capture.

We have performed a computational and experimental study in the pursuit of new GO/copolymer composites that may be promising candidates for CO2 capture materials. In particular, the interactions between CO2 and seven functional monomers, namely, acrylamide (Am), 2-aminoethyl methacrylate hydrochloride (AEMH), methyl 2-chloro acrylate (MClA), hydroxyethyl methacrylate (HEMA), glycidyl methacrylate (GMA), sodium 4-vinylbenzenesulfonate (NaSS), and 4-bromostyrene (BS), have been studied by means of density functional theory (DFT) calculations. The binding energies with CO2 were analyzed first for the isolated functional monomers and then for the copolymers with MMA/BA, in terms of the interactions with the different moieties and structural rearrangements. The incorporation of MMA/BA enhances the CO2 affinity due to the presence of the carbonyl groups of the MMA/BA copolymer, allowing for new interactions with CO2. The only exception is Am, in which case CO2 does not interact with carbonyl groups. From this study, four copolymers including AEMH, GMA, HEMA, and NaSS have been selected for analysis of adsorption of the GO/polymer composites.

The calculations show that the binding energy for CO2 is increased for the composites including AEMH and GMA, especially for the latter, while for NaSS and HEMA is decreased. This striking result, where the cooperative effect of the GO seems to be absent, may be explained by a remarkably strong interaction between the copolymer and GO, in such a way that the CO2, while having limited interaction with GO, now shows a less effective interaction with the copolymer.

These four composites have also been synthesized experimentally and the CO2 selective adsorption over N2 has been evaluated. Even with a synthesis procedure designed to decrease the development of highly porous structures, only GMA and AEHM functionalized composites presented very similar textural structures allowing us to study the effect of different functionalities on the selective CO2 capture. The trends observed experimentally were exactly predicted by the theoretical study. Moreover, the difference between the theoretical and experimental results in the case of HEMA functionalized composite confirmed the importance of the effect of morphology. In summary, this work shows that DFT is a useful tool for screening functional monomers for the design and synthesis of new GO-based materials for CO2 capture.

Some of the datasets used and/or analyzed during the current study are presented in the Supplementary Information file. The datasets used and/or analyzed during the current study, which are not shown in the Supplementary Information file, are available from the corresponding author upon reasonable request.

IPCC, 2018: Global Warming of 1.5°C. An IPCC Special Report on the impacts of global warming of 1.5°C above pre-industrial levels and related global greenhouse gas emission pathways, in the context of strengthening the global response to the threat of climate change, sustainable development, and efforts to eradicate poverty (eds. Masson-Delmotte, V. et al).

Smit, B., Reimer, J. R., Oldenburg, C. M. & Bourg, I. C. Introduction to Carbon Capture and Sequestration (Imperial College Press, 2014).

Book Google Scholar

Yu, C.-H., Huang, C.-H. & Tan, C.-S. A review of CO2 capture by absorption and adsorption. Aerosol Air Qual. Res. 12, 745–769 (2012).

Article CAS Google Scholar

Shi, J. et al. Enzymatic conversion of carbon dioxide. Chem. Soc. Rev. 44, 5981–6000 (2015).

Article CAS PubMed Google Scholar

Dong, G.-X., Li, H.-Y. & Chen, V. Challenges and opportunities for mixed-matrix membranes for gas separation. J. Mater. Chem. A 1, 4610–4630 (2013).

Article CAS Google Scholar

Lu, A.-H. & Hao, G.-P. Porous materials for carbon dioxide capture. Annu. Rep. Prog. Chem. Sect. A: Inorg. Chem. 109, 484–503 (2013).

Article CAS Google Scholar

Wang, J. et al. Recent advances in solid sorbents for CO2 capture and new development trends. Energy Environ. Sci. 7, 3478–3518 (2014).

Article CAS Google Scholar

Lee, S.-Y. & Park, S.-J. A review on solid adsorbents for carbon dioxide capture. J. Ind. Eng. Chem. 23, 1–11 (2015).

Article Google Scholar

Chew, T. L., Ahmad, A. L. & Bhatia, S. Ordered mesoporous silica (OMS) as an adsorbent and membrane for separation of carbon dioxide (CO2). Adv. Colloid Interface Sci. 153, 43–57 (2010).

Article CAS PubMed Google Scholar

Li, J.-R. et al. Carbon dioxide capture-related gas adsorption and separation in metal-organic frameworks. Coord. Chem. Rev. 255, 1791–1823 (2011).

Article CAS Google Scholar

Najafabadi, A. T. Emerging applications of graphene and its derivatives in carbon capture and conversion: Current status and future prospects. Renew. Sust. Energ. Rev. 41, 1515–1545 (2015).

Article Google Scholar

Chandra, V. et al. Highly selective CO2 capture on N-doped carbon produced by chemical activation of polypyrrole functionalized graphene sheets. Chem. Commun. 48, 735–737 (2012).

Article CAS Google Scholar

Kemp, K. C., Chandra, V., Saleh, M. & Kim, K. S. Reversible CO2 adsorption by an activated nitrogen doped graphene/polyaniline material. Nanotechnology 24, 235703 (2013).

Article ADS PubMed Google Scholar

Saleh, M., Chandra, V., Kemp, K. C. & Kim, K. S. Synthesis of N-doped microporous carbon via chemical activation of polyindole-modified graphene oxide sheets for selective carbon dioxide adsorption. Nanotechnology 24, 255702 (2013).

Article ADS PubMed Google Scholar

Seema, H. et al. Highly selective CO2 capture by S-doped microporous carbon materials. Carbon 66, 320–326 (2014).

Article CAS Google Scholar

Zhao, Y., Ding, H. & Zhong, Q. Preparation and characterization of aminated graphite oxide for CO2 capture. Appl. Surf. Sci. 258, 4301–4307 (2012).

Article ADS CAS Google Scholar

Gadipelli, S., Lu, Y., Skipper, N. T., Yildirim, T. & Guo, Z. Design of hyperporous graphene networks and their application in solid-amine based carbon capture systems. J. Mater. Chem. A 5, 17833–17840 (2017).

Article CAS Google Scholar

Heo, J. et al. Highly permeable graphene Oxide/polyelectrolytes hybrid thin films for enhanced CO2/N2 separation performance. Sci. Rep. 7, 456 (2016).

Article ADS Google Scholar

Ha, H. et al. Gas permeation and selectivity of poly(dimethylsiloxane)/graphene oxide composite elastomer membranes. J. Membr. Sci. 518, 131–140 (2016).

Article CAS Google Scholar

Shen, J. et al. Membranes with fast and selective gas-transport channels of laminar graphene oxide for efficient CO2 capture. Angew. Chem. Int. Ed. Engl. 127, 588–592 (2015).

Article Google Scholar

Liu, F.-Q. et al. Covalent grafting of polyethyleneimine on hydroxylated three-dimensional graphene for superior CO2 capture. J. Mater. Chem. A 3, 12252–12258 (2015).

Article CAS Google Scholar

Kumar, R., Raut, D., Ramamurty, U. & Rao, C. N. Remarkable improvement in the mechanical properties and CO2 uptake of MOFs brought about by covalent linking to graphene. Angew. Chem. Int. Ed. 55, 7857–7861 (2016).

Article CAS Google Scholar

Politakos, N. et al. Graphene-based monolithic nanostructures for CO2 capture. Ind. Eng. Chem. Res. 59, 8612–8621 (2020).

Article CAS Google Scholar

Politakos, N. et al. Reduced graphene oxide/polymer monolithic materials for selective CO2 capture. Polymers 12, 936 (2020).

Article CAS PubMed Central Google Scholar

Elliott, J. A. Novel approaches to multiscale modelling in materials science. Int. Mater. Rev. 56, 207–225 (2011).

Article CAS Google Scholar

Ruipérez, F. Application of quantum chemical methods in polymer chemistry. Int. Rev. Phys. Chem. 38, 343–403 (2019).

Article Google Scholar

Cinke, M., Li, J., Bauschlicher, C. W. Jr., Ricca, A. & Meyyappa, M. CO2 adsorption in single-walled carbon nanotubes. Chem. Phys. Lett. 376, 761–766 (2003).

Article ADS CAS Google Scholar

Allouche, A. & Ferro, Y. Dissociative adsorption of atmospheric molecules at vacancies on the graphite (0 0 0 1) surface of samples exposed to plasma. J. Nucl Mater. 363, 117–121 (2007).

Article ADS Google Scholar

Xu, S. C., Irle, S., Musaev, D. G. & Lin, M. C. Quantum chemical prediction of reaction pathways and rate constants for dissociative adsorption of COx and NOx on the graphite (0001) surface. J. Phys. Chem. B. 110, 21135–21144 (2006).

Article CAS PubMed Google Scholar

Cabrera-Sanfelix, P. Adsorption and reactivity of CO2 on defective graphene sheets. J. Phys. Chem. A 113, 493–498 (2009).

Article CAS PubMed Google Scholar

Liu, Y. & Wilcox, J. CO2 adsorption on carbon models of organic constituents of gas shale and coal. Environ. Sci. Technol. 45, 809–814 (2011).

Article ADS CAS PubMed Google Scholar

Liu, Y. & Wilcox, J. Effects of surface heterogeneity on the adsorption of CO2 in microporous carbons. Environ. Sci. Technol. 46, 1940–1947 (2012).

Article ADS CAS PubMed Google Scholar

Sun, F. et al. Highlighting the role of nitrogen doping in enhancing CO2 uptake onto carbon surface: A combined experimental and computational analysis. J. Mater. Chem. A 4, 18248–18252 (2016).

Article CAS Google Scholar

Wang, Y. et al. Nitrogen and oxygen codoped porous carbon with superior CO2 adsorption performance: A combined experimental and DFT calculation study. Ind. Eng. Chem. Res. 58, 13390–13400 (2019).

Article CAS Google Scholar

Dasgupta, T., Punnathanam, S. N. & Ayappa, K. G. Effect of functional groups on separating carbon dioxide from CO2/N2 gas mixtures using edge functionalized graphene nanoribbons. Chem. Eng. Sci. 121, 279–291 (2015).

Article CAS Google Scholar

Patel, H. A. et al. Directing the structural features of N2-phobic nanoporous covalent organic polymers for CO2 capture and separation. Chem. Eur. J. 20, 772–780 (2014).

Article CAS PubMed Google Scholar

Das, S. K., Bhanja, P., Kundu, S. K., Mondal, S. & Bhaumik, A. Role of surface phenol-OH groups in N-rich porous organic polymers for enchancing the CO2 uptake and CO2/N2 selectivity: Experimental and computational study. ACS Appl. Mater. Interfaces 10, 23812–23824 (2018).

Article Google Scholar

Ullah, R., Patel, H., Aparicio, S., Yavuz, C. T. & Atilhan, M. A combined experimental and theoretical study on gas adsorbtion performance of amine and amide porous polymers. Micropor. Mesopor. Mat. 279, 61–72 (2021).

Article Google Scholar

Chai, J.-D. & Head-Gordon, M. Long-range corrected hybrid density functionals with damped atom–atom dispersion corrections. Phys. Chem. Chem. Phys. 10, 6615–6620 (2008).

Article CAS PubMed Google Scholar

Gaussian 16, Revision C.01. Frisch, M. J. et al. Gaussian, Inc., Wallingford CT, (2016).

Boys, S. F. & Bernardi, F. The calculation of small molecular interactions by the difference of separate total energies. Some procedures with reduced errors. Mol. Phys. 19, 553–566 (1970).

Article ADS CAS Google Scholar

Simon, S., Duran, M. & Dannenberg, J. J. How does basis set superposition error change the potential surfaces for hydrogen bonded dimers?. J. Chem. Phys. 105, 11024–11031 (1996).

Article ADS CAS Google Scholar

Download references

Spanish Government (BES-2017-080221) is gratefully acknowledged for its financial support. The authors would like to acknowledge the contribution of the COST Action CA 15107. The authors thank SGIker (UPV/EHU, ERDF, EU) for technical and human support.

POLYMAT and Departamento de Química Aplicada, Facultad de Ciencias Químicas, University of the Basque Country UPV/EHU, Joxe Mari Korta Zentroa, Tolosa Hiribidea, 72, 20018, Donostia-San Sebastián, Spain

Branislav Stankovic, Iranzu Barbarin, Oihane Sanz & Radmila Tomovska

Faculty of Physical Chemistry, University of Belgrade, Studentski Trg 12-16, Belgrade, 11050, Republic of Serbia

Branislav Stankovic

IKERBASQUE, Basque Foundation for Science, María Díaz de Haro 3, 48013, Bilbao, Spain

Radmila Tomovska

POLYMAT and Physical Chemistry Department, Faculty of Pharmacy, University of the Basque Country, 01006, Vitoria-Gasteiz, Spain

Fernando Ruipérez

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

B.S. performed the density functional theory (DFT) calculations and wrote the first version of the manuscript; I.B. performed the experimental study, synthesis of all the materials and wrote the comments on the figures and tables placed in the manuscript; O.S. performed the adsoprtion-desorption measurements, prepared the figures and discussions related to them; R.T. supervised the experimental work and corrected the final version of the manuscript; F.R supervised the theoretical part of the work and corrected the final version of the manuscript. All authors reviewed the manuscript.

Correspondence to Radmila Tomovska or Fernando Ruipérez.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Stankovic, B., Barbarin, I., Sanz, O. et al. Experimental and theoretical study of the effect of different functionalities of graphene oxide/polymer composites on selective CO2 capture. Sci Rep 12, 15992 (2022). https://doi.org/10.1038/s41598-022-20189-5

Download citation

Received: 04 April 2022

Accepted: 09 September 2022

Published: 26 September 2022

DOI: https://doi.org/10.1038/s41598-022-20189-5

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

Scientific Reports (2023)

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.