Low Temperature Properties of Acetonitrile Confined in MCM-41

Probing Self-Diffusion of Guest Molecules in a Covalent Organic Framework: Simulation and Experiment

Covalent organic frameworks (COFs) are a class of porous materials whose sorption properties have so far been studied primarily by physisorption. Quantifying the self-diffusion of guest molecules inside their nanometer-sized pores allows for a better understanding of confinement effects or transport limitations and is thus essential for various applications ranging from molecular separation to catalysis. Using a combination of pulsed field gradient nuclear magnetic resonance measurements and molecular dynamics simulations, we have studied the self-diffusion of acetonitrile and chloroform in the 1D pore channels of two imine-linked COFs (PI-3-COF) with different levels of crystallinity and porosity. The higher crystallinity and porosity sample exhibited anisotropic diffusion for MeCN parallel to the pore direction, with a diffusion coefficient of Dpar = 6.1(3) × 10-10 m2 s-1 at 300 K, indicating 1D transport and a 7.4-fold reduction in self-diffusion compared to the bulk liquid. This finding aligns with molecular dynamics simulations predicting 5.4-fold reduction, assuming an offset-stacked COF layer arrangement. In the low-porosity sample, more frequent diffusion barriers result in isotropic, yet significantly reduced diffusivities (DB = 1.4(1) × 10-11 m2 s-1). Diffusion coefficients for chloroform at 300 K in the pores of the high- (Dpar = 1.1(2) × 10-10 m2 s-1) and low-porosity (DB = 4.5(1) × 10-12 m2 s-1) samples reproduce these trends. Our multimodal study thus highlights the significant influence of real structure effects such as stacking faults and grain boundaries on the long-range diffusivity of molecular guest species while suggesting efficient intracrystalline transport at short diffusion times.

Fluids and Electrolytes under Confinement in Single-Digit Nanopores

Confined fluids and electrolyte solutions in nanopores exhibit rich and surprising physics and chemistry that impact the mass transport and energy efficiency in many important natural systems and industrial applications. Existing theories often fail to predict the exotic effects observed in the narrowest of such pores, called single-digit nanopores (SDNs), which have diameters or conduit widths of less than 10 nm, and have only recently become accessible for experimental measurements. What SDNs reveal has been surprising, including a rapidly increasing number of examples such as extraordinarily fast water transport, distorted fluid-phase boundaries, strong ion-correlation and quantum effects, and dielectric anomalies that are not observed in larger pores. Exploiting these effects presents myriad opportunities in both basic and applied research that stand to impact a host of new technologies at the water-energy nexus, from new membranes for precise separations and water purification to new gas permeable materials for water electrolyzers and energy-storage devices. SDNs also present unique opportunities to achieve ultrasensitive and selective chemical sensing at the single-ion and single-molecule limit. In this review article, we summarize the progress on nanofluidics of SDNs, with a focus on the confinement effects that arise in these extremely narrow nanopores. The recent development of precision model systems, transformative experimental tools, and multiscale theories that have played enabling roles in advancing this frontier are reviewed. We also identify new knowledge gaps in our understanding of nanofluidic transport and provide an outlook for the future challenges and opportunities at this rapidly advancing frontier.

Simulations of the IR and Raman spectra of water confined in amorphous silica slit pores

Water in nano-scale confining environments is a key element in many biological, material, and geological systems. The structure and dynamics of the liquid can be dramatically modified under these conditions. Probing these changes can be challenging, but vibrational spectroscopy has emerged as a powerful tool for investigating their behavior. A critical, evolving component of this approach is a detailed understanding of the connection between spectroscopic features and molecular-level details. In this paper, this issue is addressed by using molecular dynamics simulations to simulate the linear infrared (IR) and Raman spectra for isotopically dilute HOD in D₂O confined in hydroxylated amorphous silica slit pores. The effect of slit-pore width and hydroxyl density on the silica surface on the vibrational spectra is also investigated. The primary effect of confinement is a blueshift in the frequency of OH groups donating a hydrogen bond to the silica surface. This appears as a slight shift in the total (measurable) spectra but is clearly seen in the distance-based IR and Raman spectra. Analysis indicates that these changes upon confinement are associated with the weaker hydrogen-bond accepting properties of silica oxygens compared to water molecules.

The joint effect of surface polarity and concentration on the structure and dynamics of acetonitrile solution: a molecular dynamics simulation study

The interfacial properties of the acetonitrile (ACN)-water-silica interface have great implications in both liquid chromatography and heterogeneous catalysis. We have performed molecular dynamics (MD) simulations of ACN and water binary solutions to give a comprehensive study of the collective effect of silica surface polarity and ACN concentration on interfacial structures and dynamics by tuning both surface charges and ACN concentration. MD simulation results indicate that many properties in the liquid-solid interface region undergo a monotonic change as the silica surface is tuned from polar to apolar due to the weakening of hydrogen bonding, while their dependence on ACN concentration is presumably governed by the preferential adsorption of water at the silica surface over ACN. However, at apolar surfaces, the interfacial structures of both water and ACN behave like the liquid-vapor interface, and this resemblance leads to an enrichment of ACN at the interface as well as accelerated dynamics, which is very different from that in the bulk solution. The organization of ACN molecules at both polar and apolar surfaces can be attributed to the amphiphilic nature of ACN, by which the micro-heterogeneity domain formed can persist both in the bulk and at the liquid-solid interface. Moreover, extending diffusion analysis to the second layer of the interface shows that the interfacial transport pathways at polar surfaces are likely very different from that of apolar surfaces. These simulation results give a full spectrum description of the ACN/water liquid-solid interface at the microscopic level and will be helpful for explaining related spectroscopic experiments and understanding the microscopic mechanisms of separation protocols in current chromatography applications.

Influence of High pH on the Organization of Acetonitrile at the Silica/Water Interface Studied by Sum Frequency Generation Spectroscopy

The acetonitrile-water mixture is one of the most commonly used solvents in hydrophilic interaction chromatography, which contains silica as the solid phase. As such, the silica/acetonitrile-water interface plays a large role in the separation of compounds. Varying the pH is one way to influence retention times, particularly of ionizable solutes, yet the influence of high pH is often unpredictable. To determine how the structure of this interface changes with pH, we utilized the surface specific technique sum frequency generation (SFG). Previous SFG studies at neutral pH have suggested the existence of acetonitrile bilayers at the aqueous silica interface even at low acetonitrile mole fractions. Here we find that the SFG signal from 2900 to 3040 cm^-1 at the silica/acetonitrile-water interface increased as we adjusted the aqueous pH from near neutral to high values. This increase in signal was attributed to a greater amount of aligned water which is consistent with an increase in silica surface charge at high pH. In contrast, complementary measurements of the silica/acetonitrile-deuterium oxide interface revealed that the acetonitrile methyl mode nearly vanished as the aqueous pH was increased. This loss of methyl mode signal is indicative of a decrease in the number density of acetonitrile molecules at the interface, as orientation analysis indicates no significant change in the net orientation of the outer leaflet of the acetonitrile bilayer over the pH range studied.

Nitriles at Silica Interfaces Resemble Supported Lipid Bilayers

Nitriles are important solvents not just for bulk reactions but also for interfacial processes such as separations, heterogeneous catalysis, and electrochemistry. Although nitriles have a polar end and a lipophilic end, the cyano group is not hydrophilic enough for these substances to be thought of as prototypical amphiphiles. This picture is now changing, as research is revealing that at a silica surface nitriles can organize into structures that, in many ways, resemble lipid bilayers. This unexpected organization may be a key component of unique interfacial behavior of nitriles that make them the solvents of choice for so many applications. The first hints of this lipid-bilayer-like (LBL) organization of nitriles at silica interfaces came from optical Kerr effect (OKE) experiments on liquid acetonitrile confined in the pores of sol-gel glasses. The orientational dynamics revealed by OKE spectroscopy suggested that the confined liquid is composed of a relatively immobile sublayer of molecules that accept hydrogen bonds from the surface silanol groups and an interdigitated, antiparallel layer that is capable of exchanging into the centers of the pores. This picture of acetonitrile has been borne out by molecular dynamics simulations and vibrational sum-frequency generation (VSFG) experiments. Remarkably, these simulations further indicate that the LBL organization is repeated with increasing disorder at least 20 Å into the liquid from a flat silica surface. Simulations and VSFG and OKE experiments indicate that extending the alkyl chain to an ethyl group leads to the formation of even more tightly packed LBL organization featuring entangled alkyl tails. When the alkyl portion of the molecule is a bulky t-butyl group, packing constraints prevent well-ordered LBL organization of the liquid. In each case, the surface-induced organization of the liquid is reflected in its interfacial dynamics. Acetonitrile/water mixtures are favored solvent systems for separations technologies such as hydrophilic interaction chromatography. Simulations had suggested that although a monolayer of water partitions to the silica surface in such mixtures, acetonitrile tends to associate with this monolayer. VSFG experiments reveal that, even at high water mole fractions, patches of well-ordered acetonitrile bilayers remain at the silica surface. Due to its ability to donate and accept hydrogen bonds, methanol also partitions to a silica surface in acetonitrile/methanol mixtures and can serve to take the place of acetonitrile in the sublayer closest to the surface. These studies reveal that liquid nitriles can exhibit an unexpected wealth of new organizational and dynamic behaviors at silica surfaces, and presumably at the surfaces of other chemically important materials as well. This behavior cannot be predicted from the bulk organization of these liquids. Our new understanding of the interfacial behavior of these liquids will have important implications for optimizing a wide range of chemical processes in nitrile solvents.

Non-Ergodic Electron Transfer in Mixed-Valence Charge-Transfer Complexes

Theories of activated transitions traditionally separate the dynamics and statistics of the thermal bath in the reaction rate into the preexponential frequency factor for the dynamics and a Boltzmann factor for the statistics. When the reaction rate is comparable to relaxation frequencies of the medium, the statistics loses ergodicity and the activation barrier becomes dependent on the medium dynamics. This scenario is realized for mixed-valence self-exchange electron transfer at temperatures near the point of solvent crystallization. These complexes, studied by Kubiak and coworkers, display anti-Arrhenius temperature dependence on lowering temperature when approaching crystallization; that is, the reaction rate increases nonlinearly in Arrhenius coordinates. Accordingly, the solvent relaxation slows down following a power temperature law. With this functional form for the relaxation time, nonergodic reaction kinetics accounts well for the observations.

Neutron spin echo measurements of monolayer and capillary condensed water in MCM-41 at low temperatures

Neutron spin echo measurements of monolayer and capillary condensed heavy water (D(2)O) confined in MCM-41 C10 (pore diameter 2.10 nm) were performed in a temperature range of 190-298 K. The intermediate scattering functions were analyzed by the Kohlrausch-Williams-Watts stretched exponential function. The relaxation times of confined D(2)O in the capillary condensed state follow remarkably well the Vogel-Fulcher-Tammann equation between 298 and 220 K, whereas below 220 K they show an Arrhenius type behavior. That is, the fragile-to-strong (FTS) dynamic crossover occurs, which has never been seen in experiments on bulk water. On the other hand, for monolayer D(2)O, the FTS dynamic crossover was not observed in the temperature range measured. The FTS dynamic crossover observed in capillary condensed water would take place in the central region of the pore, not near the pore surface. Because the tetrahedral-like water structure in the central region of the pore is more preserved than that near the pore surface, the FTS dynamic crossover would be concerned with the tetrahedral-like water structure.

Molecular-level mechanisms of vibrational frequency shifts in a polar liquid

A molecular-level analysis of the origins of the vibrational frequency shifts of the CN stretching mode in neat liquid acetonitrile is presented. The frequency shifts and infrared spectrum are calculated using a perturbation theory approach within a molecular dynamics simulation and are in good agreement with measured values reported in the literature. The resulting instantaneous frequency of each nitrile group is decomposed into the contributions from each molecule in the liquid and by interaction type. This provides a detailed picture of the mechanisms of frequency shifts, including the number of surrounding molecules that contribute to the shift, the relationship between their position and relative contribution, and the roles of electrostatic and van der Waals interactions. These results provide insight into what information is contained in infrared (IR) and Raman spectra about the environment of the probed vibrational mode.

Effect of confinement on the fluid properties of ammonia in mesopores of MCM-41 and SBA-15

The effect of pore size on capillary condensation and solid-liquid phase changes of ammonia in MCM-41 and SBA-15 was studied by adsorption and FTIR measurements of condensed phases at low temperatures. Adsorption isotherms are all typical type IV on the fully hydroxylated surfaces, without hysteresis loops in the smaller pores (d < 2.4 nm). In the larger pores, hysteresis loops appear at lower temperatures and disappear with increasing temperature, i.e., the capillary critical phenomenon was detected (hysteresis critical point). The criticality of the adsorption hysteresis loop is very similar to that for nonpolar nitrogen in mesopores of various shapes, suggesting that this is a universal phenomenon among fluids in mesopores. Freezing and melting of capillary-condensed ammonia were observed by FTIR spectroscopy. The melting temperature of capillary-condensed ammonia decreased with decreasing pore size, which is similar in the behavior of freezing. In the smaller pores (d < 2.4 nm); however, ammonia was not frozen. It is suggested that the capillary-condensed inner part, i.e., inside the ammonia monolayer, is affected too much by the pore wall and/or is too small in volume to crystallize. In the larger pores of SBA-15, crystallization is remarkably segregated from ammonia molecules strongly coordinated to surface hydroxyls.

Grand canonical Monte Carlo simulations of acetonitrile filling of silica pores of varying hydrophilicity/hydrophobicity

Grand canonical Monte Carlo simulations have been used to determine the equilibrium density of acetonitrile in model amorphous silica pores with varying radius and surface chemistry. Pores of diameter approximately 2-4 nm were considered with different ratios of surface -OH moieties to -OC(CH(3))(3) groups. The calculations found that the acetonitrile density in the interior of all the pores is essentially identical with that of the bulk liquid. On the other hand, a slightly elevated liquid density is observed near the pore surface for pores with only -OH surface moieties. Replacement of surface -OH groups with -OC(CH(3))(3) units lengthens the liquid/pore interfacial region as acetonitrile molecules can insert themselves between the -OC(CH(3))(3) units. The results indicate that the major effect of changing the surface functionality comes from the differences in excluded volume rather than hydrogen-bonding effects. Finally, the choice of the acetonitrile potential can qualitatively change the results.

Phase separation of acetonitrile-water mixtures and minimizing of ice crystallites from there in confinement of MCM-41

The effect of confinement of an acetonitrile-water mixture, whose correlation length was comparable to the pore size of the mesopores of MCM-41 (d=2.4-3.6 nm), on the phase changes was studied. Used techniques were low temperature differential scanning calorimetry and Fourier transform infrared spectroscopy, where the phase separation, lowering of the freezing and melting temperatures, and phase transitions of the acetonitrile were detected. The latter occurred in the mesopores at temperatures similar to that of the pure liquid, while the melting temperature of the water in the mesopores<3.1 nm decreased markedly at higher acetonitrile contents, suggesting a marked lowering of ice crystallite size.

Enthalpy and interfacial free energy changes of water capillary condensed in mesoporous silica, MCM-41 and SBA-15

The effect of confinement on the solid-liquid phase transitions of water was studied by using DSC and FT-IR measurements. Enthalpy changes upon melting of frozen water in MCM-41 and SBA-15 were determined as a function of pore size and found to decrease with decreasing pore size. The melting point also decreased almost monotonically with a decrease in pore size. Analysis of the Gibbs-Thomson relation on the basis of the thermodynamic data showed that there were two stages of interfacial free energy change after the constant region, i.e., below a pore size of 6.0 nm: a gradual decrease down to 3.4 nm and another decrease after a small jump upward. This fact demonstrates that the simple Gibbs-Thomson relation, i.e., a linear relation between the melting point change and the inverse pore size, is limited to the range not far from the melting point of bulk water. FT-IR measurements suggest that the decrease in enthalpy change and interfacial free energy change with decreasing pore size reflect the similarity of the structures of both liquid and solid phases of water in smaller pores at lower temperatures.