Engineering Research Center of Optical Instrument and System, Ministry of Education and Shanghai Key Lab of Modern Optical System, University of Shanghai for Science and Technology, No. 516 Jun Gong Road, Shanghai 200093, People’s Republic of China
2.
School of Electrical and Electronic Engineering, Nanyang Technological University, Singapore 639798, Singapore
Received:
30 October 2016
Accepted:
21 January 2017
Published:
25 January 2017
Au/graphene oxide/Ag sandwich structure thin film was fabricated. The effects of graphene oxide (GO) and bimetal on the structure and optical properties of metal silver films were investigated by X-ray diffraction (XRD), optical absorption, and Raman intensity measurements, respectively. Compared to silver thin film, Au/graphene oxide/Ag sandwich structure composite thin films were observed with wider optical absorption peak and enhanced absorption intensity. The Raman signal for Rhodamine B molecules based on the Au/graphene oxide/Ag sandwich nanostructure substrate were obviously enhanced due to the bimetal layer and GO layer with tunable absorption intensity and fluorescence quenching effects.
Citation: Ruijin Hong, Jialin Ji, Chunxian Tao, Daohua Zhang, Dawei Zhang. Fabrication of Au/graphene oxide/Ag sandwich structure thin film and its tunable energetics and tailorable optical properties[J]. AIMS Materials Science, 2017, 4(1): 223-230. doi: 10.3934/matersci.2017.1.223
Related Papers:
[1]
Ahmed Z. Abdullah, Adawiya J. Haider, Allaa A. Jabbar .
Pure TiO2/PSi and TiO2@Ag/PSi structures as controllable sensor for toxic gases. AIMS Materials Science, 2022, 9(4): 522-533.
doi: 10.3934/matersci.2022031
[2]
Angaraj Singh, Ajitanshu Vedrtnam, Kishor Kalauni, Aman Singh, Magdalena Wdowin .
Synthesis routes of zeolitic imidazolate framework-8 for CO2 capture: A review. AIMS Materials Science, 2025, 12(1): 118-164.
doi: 10.3934/matersci.2025009
[3]
Laís Weber Aguiar, Cleiser Thiago Pereira da Silva, Hugo Henrique Carline de Lima, Murilo Pereira Moises, Andrelson Wellington Rinaldi .
Evaluation of the synthetic methods for preparing metal organic frameworks with transition metals. AIMS Materials Science, 2018, 5(3): 467-478.
doi: 10.3934/matersci.2018.3.467
[4]
Jamal Alnofiay, Ahmed Al-Shahrie, Elsayed Shalaan .
Green synthesis of high-performance gallium oxide supercapacitor: A path to outstanding energy density. AIMS Materials Science, 2024, 11(6): 1065-1082.
doi: 10.3934/matersci.2024051
[5]
Akira Nishimura, Tadaki Inoue, Yoshito Sakakibara, Masafumi Hirota, Akira Koshio, Fumio Kokai, Eric Hu .
Optimum molar ratio of H2 and H2O to reduce CO2 using Pd/TiO2. AIMS Materials Science, 2019, 6(4): 464-483.
doi: 10.3934/matersci.2019.4.464
[6]
Akira Nishimura, Ryuki Toyoda, Daichi Tatematsu, Masafumi Hirota, Akira Koshio, Fumio Kokai, Eric Hu .
Optimum reductants ratio for CO2 reduction by overlapped Cu/TiO2. AIMS Materials Science, 2019, 6(2): 214-233.
doi: 10.3934/matersci.2019.2.214
[7]
Masatoshi Sakairi, Hirotaka Mizukami, Shuji Hashizume .
Effects of solution composition on corrosion behavior of 13 mass% Cr martensitic stainless steel in simulated oil and gas environments. AIMS Materials Science, 2019, 6(2): 288-300.
doi: 10.3934/matersci.2019.2.288
[8]
Pravina Kamini G., Kong Fah Tee, Jolius Gimbun, Siew Choo Chin .
Biochar in cementitious material—A review on physical, chemical, mechanical, and durability properties. AIMS Materials Science, 2023, 10(3): 405-425.
doi: 10.3934/matersci.2023022
[9]
Toni Varila, Henrik Romar, Tero Luukkonen, Ulla Lassi .
Physical activation and characterization of tannin-based foams enforced with boric acid and zinc chloride. AIMS Materials Science, 2019, 6(2): 301-314.
doi: 10.3934/matersci.2019.2.301
[10]
Endi Suhendi, Andini Eka Putri, Muhamad Taufik Ulhakim, Andhy Setiawan, Syarif Dani Gustaman .
Investigation of ZnO doping on LaFeO3/Fe2O3 prepared from yarosite mineral extraction for ethanol gas sensor applications. AIMS Materials Science, 2022, 9(1): 105-118.
doi: 10.3934/matersci.2022007
Abstract
Au/graphene oxide/Ag sandwich structure thin film was fabricated. The effects of graphene oxide (GO) and bimetal on the structure and optical properties of metal silver films were investigated by X-ray diffraction (XRD), optical absorption, and Raman intensity measurements, respectively. Compared to silver thin film, Au/graphene oxide/Ag sandwich structure composite thin films were observed with wider optical absorption peak and enhanced absorption intensity. The Raman signal for Rhodamine B molecules based on the Au/graphene oxide/Ag sandwich nanostructure substrate were obviously enhanced due to the bimetal layer and GO layer with tunable absorption intensity and fluorescence quenching effects.
1. Introduction
Metal--organic materials (MOMs) are a class of synthesized, often porous, and crystalline materials that have comprised the focus of a large amount of experimental and theoretical studies for the past few decades [1,2,3]. The application of MOMs is very diverse and has become essentially ubiquitous in scientific research, ranging from gas sorption [3,4,5,6,7,8,9,10,11,12,13,14,15,16,17], separation [4,18,19], catalysis [20,21,22,23], sensing [24,25], photoinduced electron transfer [26,27,28,29] and biological uses [30,31]. The major reason for this is that there are conceivably infinite variations of MOM structures that can be synthesized by altering their component metal ions, organic linkers, secondary building units (SBUs) [32,33], solvents, and thermodynamic or physical conditions.
This study focuses on explaining the sorption properties of CO2 and C2H2 in two isostructural porous materials. A method that separates these gases facilely from a mixture is highly sought for since CO2 is one of the main contaminants in C2H2. The latter is an important feedstock chemical that is used for myriad products, such as pharmaceuticals, plastics, and fuels. There is also a growing need for technology that can mitigate the effects that CO2 has on the environment. Indeed, the detrimental results of excess atmospheric CO2 are both already observed and predicted to propogate in human health and the environment if steps are not taken to combat it [34]. A material that can effectively sorb and/or separate these two gases is therefore desirable and an understanding of the mechanism of such properties in candidate materials can pave the way for better, safer, and cheaper porous materials.
Recent studies involved examining C2H2/CO2 separation in hybrid ultramicroporous materials (HUMs), a subclass of MOMs that combine ultramicropores with strong electrostatics from inorganic anions that line the pore surface [35,36]. The materials of interest in this particular study are classified as molecular porous materials (MPMs), a class of molecular solids that are held together through weak noncovalent interactions (e.g., hydrogen bonding) as opposed to the conventional fully covalently bonded interactions found in most MOMs. The combination of adenine and Cu+2 salts of Cl and Br forms self-assembled [Cu2(adenine)4Cl2]Cl2 and [Cu2(adenine)4Br2]Br2, known as MPM-1-Cl and MPM-1-Br, respectively [37]. Orthographic views of the 1 × 1 × 2 supercell of MPM-1-Br are shown in Figure 1.
Figure 1. (a) Orthographic a-axis view, (b) b-axis view, and (c) c-axis view of the 1 × 1 × 2 supercell of MPM-1-Br used for the simulations in this study. The crystal structure of MPM-1-Cl is very similar. Atom colors are C = cyan; N = blue; H = white; Cu = brown; Br = red.
The syntheses of MPM-1-Cl and MPM-1-Br are reported in references [38,39], respectively. Both MPMs feature an interesting hydrogen-bonding network in which (1) four hydrogen atoms from four different adenine linkers are hydrogen-bonded to a single halide ion and (2) an adenine linker from one [Cu2(adenine)4]4+ unit exhibits hydrogen-bonding contacts with an adenine linker from another unit. Unlike most MOMs that contain dinuclear copper paddlewheel clusters [3,40], these MPMs contain halogen atoms bonded to the Cu2+ ions of the copper paddlewheels at the axial position to create a saturated metal center. The two MPMs have a minor kink in the X--Cu--Cu--X (X = Cl or Br) motif such that the angle is not 180∘. Previous theoretical studies on MPMs [41] involved comparing the CO2 and H2 sorption mechanism in MPM-1-Cl and MPM-1-TIFSIX [42], with latter containing TiF2−6 ions in place of the halide ions. Overall, MPMs are of interest due to their thermal and water stability as well as their ease of synthesis compared to most MOMs. When these attributes are combined with their decent gas sorption and separation properties, MPMs could be good and interesting candidates for materials research.
The purpose of this study is to elucidate (with atomistic resolution) the sorption properties of CO2 and C2H2 in the isostructural MPM-1-Cl and MPM-1-Br analogues. Insights into the reason for why both MPMs prefer C2H2 over CO2 will be obtained. Another interest of this work is examining how the electronic environment (and thus MPM--sorbate interaction) is effected by the use of a different halogen in the material. It will be shown that the van der Waals radius of the halogen and its polarizability play a key role in determining the location of the favorable sorption site and the energetics of binding in the material. To our knowledge, this is the first study reporting simulations of C2H2 sorption in these materials as well as being the first theoretical study of MPM-1-Br.
2. Methods
The potential energy function for MPM-1-Cl was developed by our group in previous work [41] and utilized herein. The crystal structure for MPM-1-Br was obtained from reference [39]. For all simulations in both MPMs, the sorbent atoms were treated as rigid to accomodate a constant volume ensemble system. This approximation is especially valid when phononic effects are minor [43]. As with previous work on MPM-1-Cl [41], all atoms of MPM-1-Br were given Lennard-Jones 12--6 parameters, point partial charges, and scalar point polarizabilities to model repulsion/dispersion, stationary electrostatic, and explicit polarization, respectively. The Lennard-Jones parameters for all MPM atoms were taken from the Universal Force Field (UFF)[44], while the exponentially-damped polarizabilities for all atoms other than Cu were obtained from van Duijnen et al. [45]. The polarizability parameter for Cu2+ was determined in reference [46] and used herein. The point partial charges for the chemically unique atoms were determined by averaging atomic charges computed via a least squares fit [47] to the electrostatic potential surface of representative MPM fragments (see Supporting Information, Figure S1) generated by quantum mechanical calculations using NWChem [48].
The potentials used for CO2 were the nonpolarizable CO2-PHAST [49] and TraPPE [50] models and the polarizable CO2-PHAST* model [49]. C2H2-PHAST and C2H2-PHAST* [35,36] were used to simulate C2H2 sorption in both MPMs; the * denotes the inclusion of explicit polarization. Simulations of CO2 and C2H2 sorption in MPM-1-Cl and MPM-1-Br were performed using grand canonical Monte Carlo (GCMC) methods [51]. The sorption simulations, theoretical isosteric heat of adsorption (Qst) calculation, binding site analysis, etc. were carried out using the codes developed by our group [52,53]. Further details of performing the simulations and parametrizations are presented in the Supporting Information. All relevant simulation parameters for both MPMs and the gases simulated are also provided in the Supporting Information (Tables S1--S4).
Simulated annealing calculations were performed using the polarizable CO2 and C2H2 potentials in both MPMs to determine the classical binding energy for each sorbate in the respective materials. This is a computational method which involves executing a normal Monte Carlo procedure (using the MPM with a single sorbate molecule) that is tempered by reducing the temperature with each accepted perturbation [54]. As the temperature decreases, the guest molecule settles into an energy minima whose magnitude can be compared with the Qst near zero loading.
3. Summary of experimental data
3.1. C2H2 uptake
C2H2 uptake in MPM-1-Cl and MPM-1-Br was experimentally measured by Xie et al. [37] at 273 and 298 K and pressures up to 1 atm; the resulting isotherms (as estimated from [37]) are shown in Figure 2a. It can be observed that as the pressure increases at 273 K, the Cl variant exhibits mildly greater affinity for C2H2 than the Br analogue. This could be due to the Cl− ions having a smaller ionic radius than Br−, resulting in a slightly larger pore volume in MPM-1-Cl. This increased accessible volume in MPM-1-Cl relative to MPM-1-Br can allow for greater number of C2H2 molecules to occupy the pore space in the material at conditions approaching saturation.
Figure 2. Experimental gas sorption isotherms of (a) C2H2 and (b) CO2 for MPM-1-X (X = Br, circles or Cl, triangles). at two different temperatures (blue = 273 K; red = 298 K) and pressures up to 1 atm. All experimental data were estimated from [37].
Interestingly, the relative uptake trend is reversed at 298 K: MPM-1-Br shows slightly greater affinity for C2H2 compared to MPM-1-Cl for all pressures considered at this temperature. The observed phenomenon may be caused by the increased kinetic energy of the C2H2 molecules and the Br variant at higher temperatures. This permits the C2H2 molecules to diffuse through and collide with MPM-1-Br at a greater rate than at 273 K. Since Br is heavier than Cl, the Br− ions likely vibrate more slowly than the Cl− ions, thus making MPM-1-Br more rigid. When the C2H2 molecules are moving more quickly at higher temperatures, the favorable geometry of the MPM-1-Br framework is less flexible, whereas the higher kinetic energy of gas molecules may induce a slightly unfavorable conformational change for the Cl analogue. This is further emphasized from the fact that the primary binding site discovered for C2H2 in MPM-1-Br has a lower potential energy and closer distance to the Cu paddlewheel than the analogous site in MPM-1-Cl; this will be explained below in Section 4.3.
3.2. CO2 uptake
Xie et al. also measured CO2 sorption in MPM-1-Cl and MPM-1-Br under the same conditions as C2H2[37]; these isotherms are presented in Figure 2b. We note that Zaworotko's group also reported CO2 sorption data in the Cl analogue previously [41,42]. It can be seen that MPM-1-Cl sorbs a greater quantity of CO2 for nearly all state points considered. This phenomenon may be explained by Cl exhibiting a higher electronegativity than Br, which allows for greater electrostatic interactions between the halide ion and the CO2 molecule. Although CO2 has no permanent dipole moment in the bulk, it displays high induced dipoles when interacting with the MPM framework [41]. The interaction between the Cl− ion and the partial positive carbon atom of the CO2 molecule is fortified with a greater electronegativity compared to Br−. Our electronic structure calculations in both MPMs confirmed that the magnitude of the partial negative charges of the Cl− ions in MPM-1-Cl are greater than those for the Br− ions in MPM-1-Br (see Supporting Information, Table S1).
3.3. Isosteric heats of adsorption
Xie et al. [37] derived the experimental Qst values for C2H2 and CO2 in MPM-1-Cl and MPM-1-Br by applying the virial method [55] to the corresponding experimental isotherms. These Qst values are plotted as a function of gas loading in Figure 3. For both MPMs, the Qst for C2H2 is higher than that for CO2 at all loadings considered, thus indicating that these two MPMs display greater preference for C2H2 over CO2. Table 1 provides a summary of the zero-loading Qst values (Q0st) for C2H2 and CO2 in both MPMs.
Figure 3. Experimental isosteric heats of adsorption (Qst) for C2H2 (circles) and CO2 (triangles) for MPM-1-X (X = Br, blue or Cl, red). All experimental data were estimated from [37].
Table 1. Summary of data from key state points for gas sorption in MPM-1-Cl and MPM-1-Br. Experimental values [37] are compared with simulation using various potentials for each sorbate. Q0st is the Qst value for the lowest reported loading. Units are akJ mol−1, bmmol g−1.
With regards to C2H2, MPM-1-Cl exhibits a Q0st value of 28.57 kJ mol−1, while that for MPM-1-Br is 25.05 kJ mol−1. This suggests that MPM-1-Cl displays a more favorable initial loading site for C2H2 than MPM-1-Br according to experiment. The C2H2Qst plot for MPM-1-Cl exhibits an expected decrease in Qst with each successive C2H2 molecule sorbed in the MPM at low loading; however, after ~1.5 mmol g−1 loading is reached, the Qst values begin to increase and eventually surpass the Q0st value. The analogous plot for MPM-1-Br shows that the Qst continually increases as the loading increases. Overall, for both MPMs, the increase in Qst as a function of loading could be due to a combination of MPM--sorbate and cooperative sorbate--sorbate interactions [41,56]. Interestingly, while the experimentally reported C2H2Qst for MPM-1-Br is higher than that for MPM-1-Cl at high loadings (>1 mmol g−1), the reported C2H2 uptake is higher for MPM-1-Cl at 273 K and high pressures. We believe that this could be due to the empirical fitting method employed to obtain the Qst values from experimental isotherms. This is discussed in more detail in Section 4.1.3.
The Q0st value for CO2 in MPM-1-Cl and MPM-1-Br is 23.76 and 21.61 kJ mol−1, respectively. Further, the CO2Qst for MPM-1-Cl is greater than that for MPM-1-Br across the considered loading range. This implies that CO2 sorption is more favorable in MPM-1-Cl than MPM-1-Br, presumably due to the greater electronegativity and calculated partial negative charge of the Cl− ion as explained in the previous subsection. Interestingly, the CO2Qst plot for MPM-1-Cl shows the same behavior as the C2H2Qst plot for MPM-1-Br. Likewise, the shape of the Qst plot for CO2 in MPM-1-Br is similar to that for the C2H2Qst plot for MPM-1-Cl. As with C2H2, both MPMs exhibit increasing Qst with increasing loading in their CO2Qst plots, which indicates that CO2--CO2 interactions are favorable in these materials.
4. Simulation results and discussion
4.1. C2H2 sorption
4.1.1. MPM-1-Cl
The experimental and simulated C2H2 sorption isotherms in MPM-1-Cl at 273 and 298 K are displayed in Figure 4a. At 273 K, the calculated uptakes for the nonpolarizable C2H2-PHAST and polarizable C2H2-PHAST* models are in good agreement with experiment at pressures of 0.05 atm and lower. At higher pressures, however, the simulated uptakes for both potentials are significantly higher than experiment. At 298 K, both models generated sorption isotherms that show very good agreement with experiment for pressures up to 0.2 atm, while they notably oversorb experiment at greater pressures. The fact that the simulated uptakes for the two potentials are in close agreement with experiment at low pressures for both temperatures suggest proper modeling of the initial C2H2 binding sites in this material. The significant oversorption compared to experiment implies that the simulations produced inflated sorbate--sorbate interactions in the MPM relative to what is physically captured in experiment, which could be attributed to inadequacies in the present C2H2 potentials. It can be observed that the polarizable model produced uptakes that are somewhat closer to experiment at higher pressures.
Figure 4. Simulated gas sorption isotherms of (a) C2H2 and (b) CO2 for MPM-1-Cl compared with experiment at two different temperatures (blue = 273 K; red = 298 K) and pressures up to 1 atm. Experiment = pure solid line; C2H2/CO2-PHAST model = solid line with triangles; C2H2/CO2-PHAST* model = solid line with stars; TraPPE model = solid line with circles.
Figure 5a shows the simulated C2H2 sorption isotherms in MPM-1-Br at 273 and 298 K compared to experiment. As with MPM-1-Cl, simulations using the C2H2-PHAST and C2H2-PHAST* models in MPM-1-Br at both temperatures generally produced uptakes that are in good agreement with experiment at low pressures and those that are remarkably higher than experiment at higher pressures. Both models yielded comparable isotherms in MPM-1-Br at 273 K, with uptakes that are in closer agreement with experiment as the pressure approaches 1 atm. This could be a consequence of the material reaching C2H2 saturation under these conditions. At 298 K, the uptakes for the C2H2-PHAST* model is in better agreement with experiment at pressures lower than 0.1 atm.
Figure 5. Simulated gas sorption isotherms of (a) C2H2 and (b) CO2 for MPM-1-Br compared with experiment at two different temperatures (blue = 273 K; red = 298 K) and pressures up to 1 atm. Experiment = pure solid line; C2H2/CO2-PHAST model = solid line with triangles; C2H2/CO2-PHAST* model = solid line with stars; TraPPE model = solid line with circles.
As shown in Table 1, the experimental C2H2 uptakes for MPM-1-Br at both temperatures are higher than those for MPM-1-Cl at low pressures (0.05 atm), whereas such uptakes for the former are lower than those for the latter at higher pressures (1.0 atm). This trend is consistent with what was observed in the simulations when utilizing both C2H2 potentials. Therefore, MPM-1-Br exhibits greater interactions with the C2H2 molecules than MPM-1-Cl at low loadings, which could be explained by the Br− ions having a larger polarizability and van der Waals parameters than Cl− ions (see Supporting Information, Table S2). Nevertheless, at higher pressures, MPM-1-Br displays lower C2H2 uptake than the Cl analogue because less accessible space is available in the material under these conditions due to the larger ionic radius of the Br− ions. Indeed, the shape of the simulated C2H2 sorption isotherms for MPM-1-Br suggests that the material is near C2H2 saturation at ca. 3.5 mmol g−1 loading.
4.1.3. Isosteric heats of adsorption
The simulated C2H2Qst values for MPM-1-Cl and MPM-1-Br are compared with their corresponding experimental Qst plots in Figures 6a and 7a, respectively. It can be observed that simulations using the two models captured the apparent increase in Qst as the loading increases for both MPMs. This shape is consistent with the experimental C2H2Qst plot for MPM-1-Br. Overall, the trend in the simulated Qst values for MPM-1-Cl and MPM-1-Br support the experimental observation that both MPMs tend to facilitate C2H2 sorption such that initially loaded molecules tend to cooperatively interact with each other to lower the potential energy of subsequent guest molecules, thereby leading to enhanced sorbate--sorbate interactions at higher loadings. When considering all loadings, both models produced Qst values that are within the vicinity of experiment for both MPMs.
Figure 6. Simulated isosteric heat of sorption (Qst) for (a) C2H2 and (b) CO2 for MPM-1-Cl compared with experiment. Experiment = solid black line; Simulation = solid blue line with symbols (C2H2/CO2-PHAST = triangles; C2H2/CO2-PHAST* = stars; TraPPE model = circles).
Figure 7. Simulated isosteric heat of sorption (Qst) for (a) C2H2 and (b) CO2 for MPM-1-Br compared with experiment. Experiment = solid black line; Simulation = solid blue line with symbols (C2H2/CO2-PHAST = triangles; C2H2/CO2-PHAST* = stars; TraPPE model = circles).
The theoretical Q0st value for the C2H2-PHAST and C2H2-PHAST* models in MPM-1-Cl are 25.49 and 25.65 kJ mol−1, respectively, which are lower than the corresponding experimental value of 28.57 kJ mol−1 (Table 1). Simulations in MPM-1-Br yielded Q0st values of 30.61 and 27.55 kJ mol−1 for the nonpolarizable and polarizable models, respectively, which are higher than the experimental value for this MPM (25.05 kJ mol−1). In general, our simulations suggest that MPM-1-Br exhibits a higher Q0st value for C2H2 than MPM-1-Cl, whereas the opposite trend is observed according to experiment. From a molecular point of view, it is expected that MPM-1-Br would display a greater affinity for C2H2 due to the higher polarizability of the Br− ions. Such an effect is probably captured in our simulations that include both implicit and explicit polarization. This behavior is not reflected in the reported experimental C2H2Qst for both MPMs even though the experimental isotherms at both temperatures suggest that MPM-1-Br outperforms MPM-1-Cl for C2H2 uptake at low pressures (see Table 1). It is notable that the experimental Qst values are derived from an empirical fitting procedure on the measured sorption isotherms, whereas the theoretical values are obtained directly from GCMC simulation [57]. We could therefore attribute the discrepancy in the relative Q0st trend between experiment and simulation to the empirical fitting that was implemented to extract the experimental Qst values.
4.2. CO2 sorption
4.2.1. MPM-1-Cl
The simulated CO2 uptakes for three different potentials (CO2-PHAST, CO2-PHAST*, and TraPPE) in MPM-1-Cl at 273 and 298 K are shown in Figure 4b. The isotherms produced by the CO2-PHAST and TraPPE models oversorb experiment for nearly all state points considered. On the other hand, the isotherms for the CO2-PHAST* model are in excellent agreement with experiment across the considered pressure range at both temperatures. The strong agreement between the experimental and simulated isotherms for the CO2-PHAST* model demonstrates that this potential is highly accurate for simulations in heterogeneous media [49]. The TraPPE model yielded the highest calculated CO2 uptakes in MPM-1-Cl out of all models. Indeed, at 273 K, the calculated uptakes for this model are >1.5 mmol g−1 higher than experiment at pressures of 0.2 atm and beyond. Although the TraPPE model is widely used, it has been shown to generate isotherms that notably oversorb experiment in many porous materials [3,16,58,59,60].
The CO2-PHAST* model yielded the lowest CO2 uptakes in MPM-1-Cl out of all models for all considered state points even though it includes explicit many-body polarization. Generally, the theoretical gas uptake in MOMs increases with increasing terms in the potential energy function [46]. This was not the case here as the CO2-PHAST model gives uptakes that are higher than those for its polarizable counterpart. It is predicted that the inclusion of polarization produces a favorable sorbate ordering effect in this MPM, which causes the CO2 uptake to decrease compared to when only stationary electrostatic interactions are considered. A similar phenomenon was observed from theoretical studies in a HUM [56]. We note that CO2 sorption in MPM-1-Cl was investigated previously through GCMC simulations using only the CO2-PHAST* model [41]. Herein, we show the results that were obtained using different CO2 potentials in this material.
4.2.2. MPM-1-Br
A comparison of the experimental and simulated CO2 sorption isotherms for the three different potentials in MPM-1-Br at 273 and 298 K is presented in Figure 5b. The results for the three models in this MPM are consistent with that for MPM-1-Cl. Notably, the isotherms for the TraPPE model significantly oversorb experiment at both temperatures. The amount of oversorption for this model is generally proportional to the pressure. The CO2-PHAST model produces the next highest calculated CO2 uptakes for all pressures, followed by the CO2-PHAST* model. As observed in MPM-1-Cl, the CO2-PHAST* model yields uptakes that are in very good agreement with experiment for all pressures at both temperatures, thus further illustrating the accuracy of this potential for simulation CO2 sorption in porous materials. The isotherm for the nonpolarizable CO2-PHAST model is in good agreement with experiment at 298 K, but somewhat oversorbs experiment at 273 K. The increase in temperature may provide enough kinetic energy to decrease electrostatic interactions between the Br− ions and the CO2 molecule; this could explain why the CO2-PHAST model shows better agreement with experiment at the higher temperature.
As shown in Table 1, the experimental atmospheric CO2 uptakes in MPM-1-Br at 273 and 298 K are 2.54 and 1.56 mmol g−1, respectively. The corresponding values for MPM-1-Cl are 3.50 and 1.97 mmol g−1, respectively, which are higher than those for MPM-1-Br. Our simulations confirmed the experimental finding that the CO2 uptakes under these thermodynamic conditions are greater for MPM-1-Cl than MPM-1-Br; this is true regardless of the CO2 potential that was used. It is expected that CO2 sorption in these MPMs is governed by electrostatic interactions between the halide ion and the positively charged carbon atom of the CO2 molecule. MPM-1-Br displays weaker interactions with the CO2 molecules than MPM-1-Cl, probably due to the lower calculated partial negative charges for the Br− ions compared to the Cl− ions (see Supporting Information, Table S1). Lower magnitudes for the partial negative charge of the halide ion will result in decreased electrostatic interactions with the CO2 molecule.
Note, we also performed simulations of CO2 sorption in both MPMs using the EPM2 [61] model and the resulting simulated CO2 sorption isotherms are provided in the Supplementary Information, Figure S3. The calculated isotherms for this model are comparable to those for the CO2-PHAST model; both potentials generated uptakes that oversorb experiment at higher pressures. Overall, the polarizable CO2-PHAST* model produced isotherms that are in closest agreement with experiment. As a result, the binding sites that are discussed later (in Section 4.3) are based on simulations using this model.
4.2.3. Isosteric heats of adsorption
The simulated CO2 values for the CO2-PHAST, CO2-PHAST*, and TraPPE models for MPM-1-Cl and MPM-1-Br are displayed in Figures 6b and 7b, respectively. According to our simulations, the relative trend in the magnitude of the CO2Qst values in both MPMs is: TraPPE > CO2-PHAST > CO2-PHAST*. This is consistent with the pattern that was observed for the simulated uptakes for the three models in both MPMs. The simulated Qst values for the CO2-PHAST and CO2-PHAST* models in MPM-1-Cl are in good agreement with experiment as well as each other for all considered loadings. The TraPPE model Qst values for this MPM are somewhat higher than experiment across the loading range. The Qst plots for all models in MPM-1-Cl show increasing Qst with higher loadings, which is consistent with the shape for the experimental plot. Therefore, our simulations captured the expected cooperative CO2--CO2 interactions in this MPM.
In MPM-1-Br, the Qst values for all models are consistently higher than experiment for all loadings, with values produced by the polarizable CO2-PHAST* model being closer to experiment. The experimental Qst values in this MPM begins to increase when the material is partially saturated at ~1.0 mmol g−1 loading, indicating cooperativity between the CO2 molecules under these conditions. Unlike in MPM-1-Cl, the Qst plot for all three models in MPM-1-Br does not show much signs of sorbate--sorbate cooperativity in the material. Rather, the theoretical Qst values slightly decreases with increasing CO2 uptake, as does the experiment. In addition, the simulated Qst plots level off at a loading of ~0.4 mmol g−1, but there is neither a sustained or defined increasing trend in Qst as the uptakes continue to increase.
The calculated Q0st values for all three potentials in both MPMs are summarized in Table 1. The Q0st values produced by the CO2-PHAST, CO2-PHAST*, and TraPPE models in MPM-1-Cl are 24.33, 22.77 and 26.22 kJ mol−1, respectively, whereas the analogous quantities for the three potentials in MPM-1-Br are 25.57, 25.02 and 25.40 kJ mol−1, respectively. Simulations using the CO2-PHAST and CO2-PHAST* potentials suggest that the MPM--CO2 interaction is greater for MPM-1-Br than MPM-1-Cl. This is in contrast to what was observed for experiment as the empirical Q0st value for CO2 is higher for MPM-1-Cl than MPM-1-Br. A potential reason for why our simulations generated a higher initial loading CO2Qst value for MPM-1-Br is explained in next section where the binding sites are discussed.
4.3. Binding sites
Pham et al. [41] reported a binding site for CO2 in MPM-1-Cl in which the sorbate molecule aligns between two chlorine atoms, one that is bound to the copper paddlewheel complex and another which forms hydrogen-bonding interactions with four hydrogen atoms of two different adenine linkers (see Supporting Information, Figure S2). Interestingly, this site was minimally observed for simulations in MPM-1-Br. Indeed, the main binding site for CO2 in MPM-1-Br was actually between two Br atoms that are coordinated to the Cu2+ ions of two different [Cu2(adenine)4Br2]2+ complexes (Figure 8b). Notably, the positively charged carbon atom of the CO2 molecule can interact with two electronegative Br− ions simultaneously. We also observed this site in MPM-1-Cl for the simulations reported herein (Figure 8a). The larger van der Waals radius of the Br− ions causes this region to be more favorable for the CO2 molecules in MPM-1-Br than in MPM-1-Cl. This is because the larger size of the bromide ions will result in greater close-fitting interactions with the sorbate molecules. The CO2 molecule will exhibit a lower potential energy in the attractive well-depth potential region of both Br atoms as a result of the larger ionic radius. This could explain why the CO2-PHAST and CO2-PHAST* models produced CO2Q0st values that are greater for MPM-1-Br.
Figure 8. Molecular illustration of the primary binding site for CO2 observed in (a) MPM-1-Cl and (b) MPM-1-Br through simulation, which is nestled between two X (X = Cl, Br) atoms bound to two different [Cu2(adenine)4]4+ complexes. Atom colors are C = cyan; N = blue; H = white; Cu = brown; Cl = yellow; Br = red.
The electrostatic and electrodynamic (polarizable) effects serve to attract the positively charged carbon center of the CO2 molecule. The Lennard-Jones mixed σ parameter for Br--C used herein is ~3.6 Å, and the measured Br--CCO2 distance at the primary binding site is ~3.2 Å. The expected repulsion of the CO2 carbon atom is therefore counteracted by the attraction of the oxygen atoms on the molecule, which (being further away from the Br atoms) reside in the well-depth region of the dispersive potential. In addition, the calculated partial charge for this particular Br atom is −0.51 e−, which is slightly lower in magnitude than that for the analogous Cl atom in MPM-1-Cl (−0.57 e−) (see Supporting Information, Table S1). While the lower partial negative charge results in a weaker attraction between the halogen and the CO2 carbon atom, it also reduces the repulsion experienced by each electronegative oxygen atom on the sorbate. The net effect of all these interactions is a ``balancing act'' between the Br− ions in MPM-1-Br, which is a unique feature compared to the Cl analog.
MPM-1-Cl is able to sorb more CO2 overall than MPM-1-Br as demonstrated through both experimental measurements [37] and simulation (Figures 4 and 5). This could be due to the higher magnitude of the partial negative charge of the Cl− ion as mentioned above. Electrostatic interactions are especially important for CO2 sorption in both MPMs. It appears that the Coulombic attraction between the electropositive carbon atom of the CO2 molecule and the halide ion increases with increasing partial negative charge of the anion. Further, the Cl− ions tend to attract more electron density from the guest molecules via dispersive forces than the Br− ions due to having a higher electronegativity.
The primary binding site for C2H2 in both MPM-1-Cl and MPM-1-Br is located near the copper paddlewheel complex, where the sorbate molecule orients parallel to the X--Cu--Cu--X axis (Figure 9). At this site, the center-of-mass (COM) of the C2H2 molecule is positioned approximately 5.5 and 5.2 Å away from the geometric center of the Cu2+ ions in MPM-1-Cl and MPM-1-Br, respectively. This particular alignment of the C2H2 molecule about the copper paddlewheels allows the positively charged H atom on both sides of the sorbate to interact with the electronegative halide ions that are coordinated to the Cu2+ ions. The negatively charged COM of the C2H2 molecule can also interact somewhat with the Cu2+ ions of the copper paddlewheels.
Figure 9. Molecular illustration of the primary binding site for C2H2 observed in (a) MPM-1-Cl and (b) MPM-1-Br through simulation, which is aligned with the X--Cu--Cu--X (X = Cl, Br) axis of the dinuclear copper paddlewheel complex. Atom colors are C = cyan; N = blue; H = white; Cu = brown; Cl = yellow; Br = red.
Overall, there are more concurrent interactions between the C2H2 molecule and the framework about the primary binding site in both MPMs relative to that for CO2. Indeed, at the primary binding site, the C2H2 molecule exhibits two H⋅⋅⋅X− and two COM(C2H2)⋅⋅⋅Cu2+ interactions. On the other hand, there are only two C(CO2)⋅⋅⋅X− interactions for CO2 at the favorable binding site for this sorbate in both MPMs. The greater simultaneous interactions for C2H2 in MPM-1-X could explain why both MPMs display greater affinity and selectivity for C2H2 than CO2. We note that the longer COM(C2H2)⋅⋅⋅Cu2+ distance for MPM-1-Cl results in slightly weaker affinity at this site relative to MPM-1-Br, thus potentially explaining why our simulations produced lower initial C2H2Qst values for the Cl analogue.
The results for the classical binding energy calculations from simulated annealing are presented in Table 2. It is clear from these simulations that both materials favor C2H2 over CO2 (measured as a single molecule competition), which is consistent with experimental data. The energy minimum binding sites that were discovered through these simulations are the same as those shown in Figures 8 and 9. As observed in the results for the theoretical Qst values (Table 1), the calculated binding energies for MPM-1-Br are higher than those for MPM-1-Cl. This is likely due to the higher polarizability of Br used in these simulations.
Table 2. Classical binding energies (in kJ mol−1) for gases in MPM-1-Cl and MPM-1-Br as calculated by single-particle simulated annealing using the polarizable models for CO2 and C2H2. The starting temperature for all simulations was 500 K, to ensure free and full exploration of the phase space of the system.
This study aimed to elucidate the CO2 and C2H2 sorption properties and mechanisms in two isostructural MPMs, MPM-1-Cl and MPM-1-Br. Our simulations revealed that MPM-1-Cl exhibits a better overall sorption capacity for both gases than MPM-1-Br, which is consistent with the experimental measurements reported by Xie et al. [37]. This could be due to a combination of larger accessible space in the Cl analogue at higher pressures (a consequence of the smaller ionic radius for Cl−) and/or the fact that Cl is more electronegative than Br, which affords greater electrostatic interactions with the sorbate. At low loadings, however, the simulations suggest that MPM-1-Br is slightly more favorable toward both sorbates than MPM-1-Cl. We attribute this to the higher polarizablility and van der Waals parameters for the Br− ions. Notably, repulsion/dispersion interactions could dominate the sorption mechanism at such loadings.
It was discovered through our simulations that the primary binding site for C2H2 and CO2 in both MPMs are notably different between the two sorbates. Particularly, the C2H2 molecule prefers to align along the X--Cu--Cu--X axis of the copper paddlewheel to allow for favorable H⋅⋅⋅X− interactions on each side of the molecule, while the CO2 molecule seems to favor the region between the X− ions of two adjacent copper paddlewheel units. Both MPMs appear to be more selective toward C2H2 than CO2 because there are a greater number of concurrent interactions about the primary binding site for C2H2.
Herein, we demonstrated how substitution of the halide ion in two isostructural MPMs with the empirical formula [Cu2(adenine)4X2]X2 affects the gas sorption properties in the material through molecular simulations. We plan to investigate the CO2 and C2H2 sorption performance in MPM-1-TIFSIX [42] through theoretical studies in future work.
Acknowledgements
The authors acknowledge the National Science Foundation (Award No. DMR-1607989), including support from the Major Research Instrumentation Program (Award No. CHE-1531590). Computational resources were made available by a XSEDE Grant (No. TG-DMR090028) and by Research Computing at the University of South Florida. B.S. also acknowledges support from an American Chemical Society Petroleum Research Fund grant (ACS PRF 56673-ND6).
Conflict of interest
The authors declare no conflict of interest related to the content of this publication.
References
[1]
Fang N, Lee H, Sun C, et al. (2005) Sub-diffraction-limited optical imaging with a silver superlens. Science 308: 534-537. doi: 10.1126/science.1108759
[2]
Ozbay E (2006) Plasmonics: merging photonics and electronics at nanoscale dimensions. Science 311: 189-193. doi: 10.1126/science.1114849
[3]
Fleischmann M, Hendra PJ, Mcquillan AJ (1974) Raman spectra of pyridine adsorbed at a silver electrode. Chem Phys Lett 26: 163-166. doi: 10.1016/0009-2614(74)85388-1
[4]
Cheng F, Ellis AV, Voelcker NH (2012) Electrochemical synthesis of silver oxide nanowires, microplatelets and application as SERS substrate precursors. Electrochim Acta 59: 346-353. doi: 10.1016/j.electacta.2011.10.068
[5]
Braun G, Lee SJ, Dante M, et al. (2007) Surface-enhanced Raman spectroscopy for DNA detection by nanoparticle assembly onto smooth metal films. J Am Chem Soc 129: 6378-6379. doi: 10.1021/ja070514z
[6]
Grow AE, Wood LL, Claycomb JL, et al. (2003) New biochip technology for label-free detection of pathogens and their toxins. J Microbiol Meth 53: 221-233. doi: 10.1016/S0167-7012(03)00026-5
[7]
Traci RJ, Michelle DM, Christy LH, et al. (2000) Nanosphere lithography: tunable localized surface plasmon resonance spectra of silver nanoparticles. J Phys Chem B 104: 10549-10556. doi: 10.1021/jp002435e
[8]
Felidj N, Aubard J, Levi G, et al. (2003) Optimized surface-enhanced Raman scattering on gold nanoparticle arrays. Appl Phys Lett 82: 3095-3097. doi: 10.1063/1.1571979
[9]
Han Y, Lupitskyy R, Chou TM, et al. (2011) Effect of oxidation on surface-enhanced Raman scattering activity of silver nanoparticles: a quantitative correlation. Anal Chem 83: 5873-5880. doi: 10.1021/ac2005839
[10]
Gutés A, Maboudian R, Carraro C (2012) Gold-coated silver dendrites as SERS substrates with an improved lifetime. Langmuir 28: 17846-17850. doi: 10.1021/la303421s
[11]
Pande S, Ghosh SK, Praharaj S, et al. (2007) Synthesis of normal and inverted gold-silver core-shell architectures in β-Cyclodextrin and their applications in SERS. J Phys Chem C 111: 10806-10813. doi: 10.1021/jp0702393
[12]
Feng JJ, Gernert U, Sezer M, et al. (2009) Novel Au-Ag hybrid device for electrochemical SE(R)R spectroscopy in a wide potential and spectral range. Nano Lett 9: 298-303. doi: 10.1021/nl802934u
[13]
Liu F, Cao Z, Tang C, et al. (2010) Ultra-thin diamond-like carbon film coated silver nanoparticles-based substrates for surface-enhanced Raman spectroscopy. ACS Nano 4: 2643-2648. doi: 10.1021/nn100053s
[14]
Feng JJ, Gernert U, Hildebrandt P, et al. (2010) Induced SER-activity in nanostructured Ag-Silica-Au supports via long-range plasmon coupling. Adv Funct Mater 20: 1954-1961. doi: 10.1002/adfm.201000302
[15]
Mahurin SM, Bao L, Dai S (2004) Controlled layer-by-layer formation of ultrathin TiO2 on silver island films via a surface sol-gel method for Surface-Enhanced Raman Scattering. Anal Chem 76: 4531-4536. doi: 10.1021/ac049668c
[16]
Kim YK, Han SW, Min DH (2012) Graphene oxide sheath on Ag nanoparticle/graphene oxide hybrid Films as an antioxidative coating and enhancer of surface-enhanced Raman scattering. ACS Appl Mater Inter 4: 6545-6551. doi: 10.1021/am301658p
[17]
Sun C, Su KH, Valentine J, et al. (2011) Time Resolved Single-step Protease Activity Quantification Using Nanoplasmonic Resonator Sensors. ACS Nano 4: 978-984.
[18]
Kim KH, Baek YK, Jeon HJ, et al. (2012) Cylindrical posts of Ag/SiO2/Au multi-segment layer patterns for highly efficient surface enhanced Raman scattering. Nanotechnology 23: 315302. doi: 10.1088/0957-4484/23/31/315302
[19]
Su KH, Durant S, Steele JM, et al. (2006) Raman enhancement factor of a single tunable nanoplasmonic resonator. J Phys Chem B 110: 3964-3968. doi: 10.1021/jp055566u
[20]
Liu S, Tian JQ, Wang L, et al. (2011) A method for the production of reduced graphene oxide using benzylamine as a reducing and stabilizing agent and its subsequent decoration with Ag nanoparticles for enzymeless hydrogen peroxide detection. Carbon 49: 3158-3164. doi: 10.1016/j.carbon.2011.03.036
[21]
Liu S, Tian JQ, Wang L, et al. (2011) Aniline as a dispersing and stabilizing agent for reduced graphene oxide and its subsequent decoration with Ag nanoparticles for enzymeless hydrogen peroxide detection. J Colloid Interf Sci 363: 615-619. doi: 10.1016/j.jcis.2011.07.083
[22]
Hummers WS, Offeman RE (1958) Preparation of graphitic oxide. J Am Chem Soc 80: 1339-1339. doi: 10.1021/ja01539a017
[23]
Vosgröne T, Meixner AJ (2004) Surface and resonance enhanced micro-Raman spectroscopy of xanthene dyes at the single-molecule level. J Lumin 107: 13-20. doi: 10.1016/j.jlumin.2003.12.041
This article has been cited by:
1.
Douglas M. Franz, Jonathan L. Belof, Keith McLaughlin, Christian R. Cioce, Brant Tudor, Adam Hogan, Luciano Laratelli, Meagan Mulcair, Matthew Mostrom, Alejandro Navas, Abraham C. Stern, Katherine A. Forrest, Tony Pham, Brian Space,
MPMC and MCMD: Free High‐Performance Simulation Software for Atomistic Systems,
2019,
2,
2513-0390,
1900113,
10.1002/adts.201900113
Table 1. Summary of data from key state points for gas sorption in MPM-1-Cl and MPM-1-Br. Experimental values [37] are compared with simulation using various potentials for each sorbate. Q0st is the Qst value for the lowest reported loading. Units are akJ mol−1, bmmol g−1.
Table 2. Classical binding energies (in kJ mol−1) for gases in MPM-1-Cl and MPM-1-Br as calculated by single-particle simulated annealing using the polarizable models for CO2 and C2H2. The starting temperature for all simulations was 500 K, to ensure free and full exploration of the phase space of the system.