Evidence for the Direct Observation of Molecular Exchange of a Liquid at the Solid/Liquid Interface

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.

Orientational Pair Correlations and Local Structure of Benzonitrile from Molecular Dynamics Simulations with Comparisons to Experiments

We present an experimentally parametrized molecular dynamics study of single-molecule and collective orientational relaxation in neat benzonitrile through the analysis of the reorientational anisotropy and polarizability anisotropy time correlation function (PA-TCF). The simulations show that the PA-TCF is dominated by collective reorientation after 20 ps. Collective reorientation is found to be slower than single-molecule reorientation by a factor of 1.67, consistent with recent experiments. The simulations provide direct evidence of local antiparallel benzonitrile configurations. These structures, which have been the center of some debate, are responsible for the slower rate of collective versus single-molecule reorientation in the liquid. Further structural analysis indicates that significant Coulombic interactions between the nitrile group and hydrogen atoms on adjacent molecules play a role in the formation of the antiparallel structures. The single-molecule dynamics reflected in the anisotropy are complex and consist of a ballistic regime, restricted angular diffusion, and spatially anisotropic free diffusion. The principal components of the rotational diffusion tensor are independently obtained and shown to reproduce the free diffusion regime of the anisotropy for each principal axis according to the predictions of a previous theory.

The role of solute polarity on methanol-silica interfacial solvation: a molecular dynamics study

The solvation structure and dynamics of small organic molecules at the methanol-silica interface are important for understanding the reaction dynamics in heterogeneous catalysis as well as the transport mechanisms in liquid chromatography. The role of solute polarity in interfacial solvation at the methanol-silica interface has been investigated via umbrella sampling molecular dynamics (MD) simulations and 1,3-propanediol and n-pentane were selected as representative species of polar and apolar solutes. Free energy calculations reveal that it took a similar free energy cost to transfer both solute molecules from the interface to the bulk, despite the huge difference in their polarities. The 1,3-propendiol molecule can penetrate the adsorbed methanol layer and form hydrogen bonds with the silica surface with its backbone perpendicular to the silica surface, resulting in a significant slowdown of translational and rotational dynamics. Further analysis of solvent density and solute orientations suggest that at the minimum of the adsorption free energy curve, the 1,3-propanediol molecule is in a desolvated state, while n-pentane is in a solvated state. The collective effect of solute concentration has also been studied by unbiased MD simulations, and the free energy barriers of transferring the solute molecule from the interface to bulk, as well as the parallel diffusion coefficients at the interface, show a non-monotonic dependence on solute concentration, which can be related to the crowded environment in the interfacial layers.

Molecular dynamics simulations of amino acid adsorption and transport at the acetonitrile–water–silica interface: the role of side chains

The solvation and transport of amino acid residues at liquid–solid interfaces have great importance for understanding the mechanism of separation of biomolecules in liquid chromatography. This study uses umbrella sampling molecular dynamics simulations to study the adsorption and transport of three amino acid molecules with different side chains (phenylalanine (Phe), leucine (Leu) and glutamine (Gln)) at the silica–water–acetonitrile interface in liquid chromatography. Free energy analysis shows that the Gln molecule has stronger binding affinity than the other two molecules, indicating the side chain polarity may play a primary role in adsorption at the liquid–solid interface. The Phe molecule with a phenyl side chain exhibits stronger adsorption free energy than Leu with a non-polar side chain, which can be ascribed to the better solvated configuration of Phe. Further analysis of molecular orientations found that the amino acid molecules with apolar side chains (Phe and Leu) have ‘standing up’ configurations at their stable adsorption state, where the polar functional groups are close to the interface and the side chain is far from the interface, whereas the amino acid molecule with a polar side chain (Gln) chooses the ‘lying’ configuration, and undergoes a sharp orientation transition when the molecule moves away from the silica surface. Extending our simulation studies to systems with different solute concentrations reveals that there is a decrease in the adsorption free energy as well as surface diffusion as the solute concentration increases, which is related to the crowding in the interfacial layers. This simulation study gives a detailed microscopic description of amino acid molecule solvation and transport at the acetonitrile–water–silica interface in liquid chromatography and will be helpful for understanding the retention mechanism for amino acid separation.

The solvation and transport of amino acid residues at liquid–solid interfaces have great importance for understanding the mechanism of separation of biomolecules in liquid chromatography.

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.

The effect of surface polarity on the structure and collective dynamics of liquid ethanol

Typical configurations of ethanol during polarity modulation.

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.

Nanostructures and Dynamics of Isochorically Confined Amorphous Drug Mediated by Cooling Rate, Interfacial, and Intermolecular Interactions

The production and stabilization of amorphous drugs by the nanoconfinement effect has recently become a research hotspot in pharmaceutical sciences. Herein, two guest/host systems, indomethacin (IMC) and griseofulvin (GSF) confined in anodic aluminum oxide (AAO) templates with different pore diameters (25-250 nm) are investigated by differential scanning calorimetry (DSC) and broadband dielectric spectroscopy (BDS). The crystallization of the confined drugs is suppressed, and their glass transition temperatures show an evident pore-size dependency. Moreover, a combination of dielectric and calorimetric results demonstrates that the significant change in the temperature dependence of the structural relaxation time during the cooling process is attributed to the vitrification of the interfacial molecules and the local density heterogeneity under isochoric confinement. Interestingly, compared with the case of IMC/AAO, which can be described by a typical two-layer model, GSF/AAO presents an rare scenario of three glass transition temperatures under fast cooling (40-10 K/min), indicating that there exists a thermodynamic nonequilibrium interlayer between the bulk-like core and interfacial layer. In contrast, the slow cooling process (0.5 K/min) would lead confined GSF into the stable core-shell nanostructure. Using surface modification, the interfacial effect is confirmed to be an important reason for the different phenomena between these two guest/host systems, and intermolecular hydrogen bonding is also suggested to be emphasized considering the long-range effect of interfacial interactions. Our results not only provide insight into the glass transition behavior of geometrically confined supercooled liquids, but also offer a means of adjusting and stabilizing the nanostructure of amorphous drugs under two-dimensional confinement.

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.

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.

Ultrafast dynamics in complex fluids observed through the ultrafast optically-heterodyne-detected optical-Kerr-effect (OHD-OKE)

The ultrafast molecular dynamics of complex fluids have been recorded using the optically-heterodyne-detected optical-Kerr-effect (OHD-OKE). The OHD-OKE method is reviewed and some recent refinements to the method are described. Applications to a range of complex fluids, including microemulsions, polymer melts and solutions, liquid crystal and ionic liquids are surveyed. The level of detail attainable with the OHD-OKE method in these complex fluids is discussed. The prospects for future experiments are discussed.

Ultrafast Orientational Dynamics of Nanoconfined Benzene

Ultrafast optical Kerr effect spectroscopy has been used to study the orientational dynamics of benzene and benzene-d(6) confined in nanoporous sol-gel glass monoliths with a range of average pore sizes. All of the observed orientational diffusion of confined benzene is found to occur on a slower time scale than in the bulk, even in pores with diameters that are significantly larger than a benzene molecule. The orientational dynamics of benzene-d(6) are found to be inhibited to a lesser extent than those of benzene, which is attributed to the differences in wetting properties of the two liquids on silica. The decays are fit well by a sum of two exponentials, the faster of which depends on pore size. Similar results are found in pores that have been modified with trimethylsilyl groups, although the relaxation is faster than in unmodified pores. Comparison to Raman line width data for confined benzene-d(6) suggests that the liquid exhibits significant structuring at the pore walls, with the benzene molecules lying flat on the surfaces of unmodified pores.

Ultrafast Dynamics of Liquid Poly(ethylene glycol)s and Crown Ethers Studied by Femtosecond Raman-Induced Kerr Effect Spectroscopy

Ultrafast molecular dynamics of liquid poly(ethylene glycol)s, tetra(ethylene glycol), penta(ethylene glycol), and poly(ethylene glycol) with the molecular weight of 600, and crown ethers, 12-crown-4 and 15-crown-5, have been investigated by means of femtosecond optical heterodyne-detected Raman-induced Kerr effect spectroscopy. Picosecond Kerr transients of poly(ethylene glycol)s and crown ethers are characterized by a biexponential function with the time constants of about 2 and 20 ps. Both the faster and slower time constants do not vary much among the five oligo(ethylene oxide)s. Femtosecond dynamics is discussed based on the Kerr (depolarized Raman) spectra obtained by Fourier transform deconvolution analysis of the high time resolution Kerr transients. The broad low-frequency band (0-200 cm(-1)) in the Kerr spectrum is analyzed by two Brownian oscillators. The spectral shapes of linear poly(ethylene glycol) and cyclic crown ether are very different. Both the low- and high-frequency Brownian oscillators for crown ethers show lower frequency and broader spectral features than those for poly(ethylene glycol)s. The comparison of the low-frequency spectra of poly(ethylene glycol)s and crown ethers shows that the low-frequency spectrum of 15-crown-5 is closer to that of poly(ethylene glycol)s than that of 12-crown-4 is. The difference of the low-frequency spectra between poly(ethylene glycol) and crown ether is discussed with the concepts of molecular conformation and liquid density. The features of the observed intramolecular vibrational bands are also correlated with the molecular conformations.

Ultrafast molecular dynamics of liquid aromatic molecules and the mixtures with CCl4

The ultrafast molecular dynamics of liquid aromatic molecules, benzene, toluene, ethylbenzene, cumene, and 1,3-diphenylpropane, and the mixtures with CCl(4) have been investigated by means of femtosecond optical heterodyne-detected Raman-induced Kerr effect spectroscopy. The picosecond Kerr transients of benzene, toluene, ethylbenzene, and cumene and the mixtures with CCl(4) show a biexponential feature. 1,3-Diphenylpropane and the mixtures with CCl(4) show triexponential picosecond Kerr transients. The slow relaxation time constants of the aromatic molecules and the mixtures with CCl(4) are qualitatively described by the Stoke-Einstein-Debye hydrodynamic model. The ultrafast dynamics have been discussed based on the Kerr spectra in the frequency range of 0-800 cm(-1) obtained by the Fourier transform analysis of the Kerr transients. The line shapes of the low-frequency intermolecular spectra located at 0-180 cm(-1) frequency range have been analyzed by two Brownian oscillators ( approximately 11 cm(-1) and approximately 45 cm(-1) peaks) and an antisymmetric Gaussian function ( approximately 65 cm(-1) peak). The spectrum shape of 1,3-diphenylpropane is quite different from the spectrum shapes of the other aromatic molecules for the low magnitude of the low-frequency mode of 1,3-diphenylpropane and/or an intramolecular vibration. Although the concentration dependences of the low- and intermediate-frequency intermolecular modes (Brownian oscillators) do not show a significant trend, the width of high-frequency intermolecular mode (antisymmetric Gaussian) becomes narrower with the higher CCl(4) concentration for all the aromatics mixtures with CCl(4). The result indicates that the inhomogeneity of the intermolecular vibrational mode in aromatics/CCl(4) mixtures is decreasing with the lower concentration of aromatics. The intramolecular vibrational modes of the aromatic molecules observed in the Kerr spectra are also shown with the calculation results based on the density functional theory.

Solvation dynamics of formamide and N,N-dimethylformamide in aerosol OT reverse micelles

The solvation dynamics of formamide and N,N-dimethylformamide in Aerosol OT reverse micelles has been investigated in this work. The solvation dynamics of formamide and N,N-dimethylformamide in the reverse micelles is more than 100 times slower than that of the pure solvents. The solvation dynamics of formamide in the reverse micelle solution depends strongly on the molar ratio between formamide and Aerosol OT (w = [polar solvent]/[Aerosol OT]), but that of N,N-dimethylformamide in the reverse micelle solution shows a tiny w dependence. We have estimated the interaction energies of the geometry-optimized clusters of a simple model of the Aerosol OT polar headgroup (CH3SO3-) and formamide or N,N-dimethylformamide by ab initio calculations (the second-order Møller-Plesset perturbation theory) to find their interactions. The interaction energies of the mimic clusters estimated by the ab initio calculations and the features of the slow solvation dynamics and w dependence in formamide and N,N-dimethylformamide reverse micelles are discussed.

Ultrafast dynamics in aqueous polyacrylamide solutions

We have investigated the ultrafast dynamics of aqueous polyacrylamide ([-CH(2)CH(CONH(2))-](n), or PAAm) solutions using femtosecond optical heterodyne-detected Raman-induced Kerr effect spectroscopy (OHD-RIKES). The observed aqueous PAAm dynamics are nearly identical for both M(w) = 1500 and 10 000. Aqueous propionamide (CH(3)CH(2)CONH(2), or PrAm) solutions were also studied, because PrAm is an exact model for the PAAm constitutional repeat unit (CRU). The longest time scale dynamics observed for both aqueous PAAm and PrAm solutions occur in the 4-10 ps range. Over the range of concentrations from 0 to 40 wt %, the picosecond reorientation time constants for the aqueous PAAm and PrAm solutions scale linearly with the solution concentration, despite the fact that the solution shear viscosities vary exponentially from 1 to 264 cP. For a given value of solution concentration in weight percent, constant ratios of measured reorientation time constants for PAAm to PrAm are obtained. This ratio of PAAm to PrAm reorientation time constants is equal to the ratio of the volume for the PAAm constitutional repeat unit (-CH(2)CHCONH(2)-) to the molecular volume of PrAm. For these reasons, we assign the polymer reorientation dynamics to motions of the entire constitutional repeat unit, not only side group motions. Simple molecular dynamics simulations of H[-CH(2)CH(CONH(2))-](7)H in a periodic box with 180 water molecules support this assignment. Amide-amide and amide-water hydrogen-bonding interactions lead to strongly oscillatory femtosecond dynamics in the Kerr transients, peaking at 80, 410, and 750 fs.