Surface coverages of bonded-phase ligands on silica: a computational study

Anionic Dye Removal with a Thin Cationic Polyaniline Coating on Cellulosic Biomaterial

This paper reports the development of novel adsorbent materials using polyaniline (PANI) grafted onto Posidonia (POS) fibers, aimed at efficiently removing phenol red (PSP), an anionic dye, from aqueous solutions. The synthesis involved the copolymerization of aniline grafted on the surface of the POS and aniline monomer in solution, resulting in a chemically bound thin PANI layer on the POS bioadsorbent. Structural characteristics and binding affinities of these adsorbents with PANI under its emeraldine salt (POS@PANI-ES) or emeraldine base (POS@PANI-EB) forms are reported. The rapid adsorption kinetics observed are attributed to enhanced accessibility to PANI adsorption sites on the POS surface. The binding percentages of PSP to POS@PANI-ES and POS@PANI-EB materials were found to be 97 and 50%, respectively, after 15 min of contact time. The Langmuir model for localized adsorption sites and the Volmer model for nonlocalized adsorption as a mobile layer were fitted to the experimental adsorption isotherms of PSP to POS@PANI-EB and POS@PANI-ES, yielding the thermodynamic parameters of adsorption. The adsorption capacities of PSP on POS@PANI-EB and POS@PANI-ES were 37.8 and 71.5 μmol g-1, respectively. The adsorption of PSP remained above 80% at moderate salt concentrations of around 0.1 mol L-1; however, higher concentrations of NaCl and CaCl2 in PSP solutions significantly reduced the adsorption on POS@PANI-ES.

Molecular Dynamics Simulations of Amylose- and Cellulose-Based Selectors and Related Enantioseparations in Liquid Phase Chromatography

In the last few decades, theoretical and technical advancements in computer facilities and computational techniques have made molecular modeling a useful tool in liquid-phase enantioseparation science for exploring enantioselective recognition mechanisms underlying enantioseparations and for identifying selector-analyte noncovalent interactions that contribute to binding and recognition. Because of the dynamic nature of the chromatographic process, molecular dynamics (MD) simulations are particularly versatile in the visualization of the three-dimensional structure of analytes and selectors and in the unravelling of mechanisms at molecular levels. In this context, MD was also used to explore enantioseparation processes promoted by amylose and cellulose-based selectors, the most popular chiral selectors for liquid-phase enantioselective chromatography. This review presents a systematic analysis of the literature published in this field, with the aim of providing the reader with a comprehensive picture about the state of the art and what is still missing for modeling cellulose benzoates and the phenylcarbamates of amylose and cellulose and related enantioseparations with MD. Furthermore, advancements and outlooks, as well as drawbacks and pitfalls still affecting the applicability of MD in this field, are also discussed. The importance of integrating theoretical and experimental approaches is highlighted as an essential strategy for profiling mechanisms and noncovalent interaction patterns.

Consequences of Cylindrical Pore Geometry for Interfacial Phenomena in Reversed-Phase Liquid Chromatography

The interfacial phenomena behind analyte separation in a reversed-phase liquid chromatography column take place nearly exclusively inside the silica mesopores. Their cylindrical geometry can be expected to shape the properties of the chromatographic interface with consequences for the analyte density distribution and diffusivity. To investigate this topic through molecular dynamics simulations, we introduce a cylindrical pore inside a slit pore configuration, where the inner curved and outer planar silica surface bear the same bonded phase. The present model replicates an average-sized (9 nm) mesopore in an endcapped C₁₈ column equilibrated with a mobile phase of 70/30 (v/v) water/acetonitrile. Simulations performed for ethylbenzene and acetophenone show that the surface curvature shifts the bonded phase and analyte density toward the pore center, decreases the solvent density in the bonded-phase region, increases the acetonitrile excess in the interfacial region, and considerably enhances the surface diffusivity of both analytes. Overall, the cylindrical pore provides a more hydrophobic environment than the slit pore. Ethylbenzene density is decidedly increased in the cylindrical pore, whereas acetophenone density is nearly equally distributed between the cylindrical and slit pore. The cylindrical pore geometry thus sharpens the discrimination between the apolar and moderately polar analytes while enhancing the mass transport of both.

Ultrathin Molecular-Layer-by-Layer Polyamide Membranes: Insights from Atomistic Molecular Simulations

In this study, we present an atomistic simulation study of several physicochemical properties of polyamide (PA) membranes formed from interfacial polymerization or from a molecular-layer-by-layer (mLbL) on a silicon wafer. These membranes are composed of meta-phenylenediamine (MPD) and benzene-1,3,5-tricarboxylic acid chloride (TMC) for potential reverse osmosis (RO) applications. The mLbL membrane generation procedure and the force field models were validated, by comparison with available experimental data, for hydrated density, membrane swelling, and pore size distributions of PA membranes formed by interfacial polymerization. Physicochemical properties such as density, free volume, thickness, the degree of cross-linking, atomic compositions, and average molecular orientation (which is relevant for the mLbL membranes) are compared for these different processes. The mLbL membranes are investigated systematically with respect to TMC monomer growth rate per substrate surface area, MPD/TMC ratio, and the number of mLbL deposition cycles. Atomistic simulations show that the mLbL deposition generates membranes with a constant film growth if both the TMC monomer growth rate and MPD/TMC monomer ratio are kept constant. The film growth rate increases with TMC monomer growth rate or MPD/TMC ratio. Furthermore, it was found on one hand that the mLbL membrane density and free volume varies significantly with respect to the TMC monomer growth rate, while on the other hand the degree of cross-linking and the atomic composition varies considerably with the MPD/TMC ratio. Additionally, it was found that both TMC and MPD orient at a tilted angle with respect to the substrate surface, where their angular distribution and average angle orientation depend on both the TMC growth rate and the number of deposition cycles. This study illustrates that molecular simulations can play a crucial role in the understanding of structural properties that can empower the design of the next generation RO membranes created from molecular-layer-by-layer (mLbL) on a silicon wafer.

Characterization of the structural morphology of chemically modified silica prepared by surface polymerization of a mixture of long and short alkyl chains using 13C and 29Si NMR spectroscopy

A series of bonded phases were prepared by the chemical modification of silica using the surface polymerization of trifunctional and difunctional ligands, and the structural morphology was characterized by solid-state nuclear magnetic resonance (NMR) spectroscopy using cross-polarization and magic angle spinning (CP/MAS). Mixed-phase surfaces were prepared using mixtures of trifunctional long-chain (C18) ligands with trifunctional and difunctional short-chain (C1) ligands, and these surfaces were compared to the corresponding single-phase surfaces consisting of only long- or short-chain ligands. For both types of mixed-phase surfaces, the incorporation of short chains increases the overall ligand density, the density of long chains, and the degree of cross-linking between ligands compared to that of the single-phase surface consisting exclusively of long chains. When the percentage of long-chain ligands in the mixture is high, a horizontally polymerized monolayer of chains is formed on the silica surface for both trifunctional and difunctional short chains. However, essentially all of the long chains adopt a trans conformation when trifunctional short chains are used, and a significant number of gauche defects are observed for the long chains when mixed with difunctional short chains. Furthermore, the ligands on the mixed-phase surface are more rigid when the short chains are trifunctional. When the percentage of trifunctional short chains is increased, some vertical polymerization occurs, caused by the molecular stacking of the highly reactive short chains near the surface. However, this does not preclude cross-linking between the ligands necessary to seal the surface, and the degree of cross-linking is quite high, suggesting that the short chains cross-link both vertically, away from the surface, and horizontally, across the surface. No such vertical polymerization is observed for the bulkier difunctional short chains. For both trifunctional and difunctional short chains, the surface chains are more mobile, with a greater number of gauche conformations among the long chains when the percentage of short-chain ligands in the reaction mixture is increased.

Computational study of enantioseparation by amylose tris(3,5-dimethylphenylcarbamate)-based chiral stationary phase

The mechanism of chiral separation on amylose tris(3,5-dimethylphenylcarbamate) is studied with docking simulations of enantiomers by molecular dynamics. All-atom models of amylose tris(3,5-dimethylphenylcarbamate) on the modified silica gel surface were constructed for the docking simulations of metalaxyl and benalaxyl. The elution orders and energetic differences were also predicted based on the intermolecular interactions, which were in agreement with the experimental results. The radial distribution function was employed to analyze the structural features of the enantiomer-chiral stationary phase complex and used to elucidate the mechanism of chiral separation. The separation of metalaxyl and benalaxyl is mainly controlled by the hydrogen bond. And the binding sites had slight differences for the pair of enantiomers, but obvious differences between different chemicals.

Retention Mechanism in Reversed-Phase Liquid Chromatography: A Molecular Perspective

A detailed, molecular-level understanding of the retention mechanism in reversed-phase liquid chromatography (RPLC) has eluded analytical chemists for decades. Through validated, particle-based Monte Carlo simulations of a model RPLC system consisting of dimethyloctadecylsilanes at a coverage of 2.9 micro mol/m2 on an explicit silica substrate with unprotected residual silanols in contact with a water/methanol mobile phase, we show that the molecular-level retention processes for nonpolar and polar analytes, such as alkanes and alcohols, are much more complex than what has been previously deduced from thermodynamic and theoretical arguments. In contrast to some previous assumptions, the simulations indicate that both partitioning and adsorption play a key role in the separation process and that the stationary phase in RPLC behaves substantially different from a bulk hydrocarbon phase. The retention of nonpolar methylene segments is dominated by lipophilic interactions with the retentive phase, while solvophilic interactions are more important for the retention of the polar hydroxyl group.

Chain conformation and solvent partitioning in reversed-phase liquid chromatography: Monte Carlo simulations for various water/methanol concentrations

Many structural models for the stationary phase in reversed-phase liquid chromatography (RPLC) systems have been suggested from thermodynamic and spectroscopic measurements and theoretical considerations. To provide a molecular picture of chain conformation and solvent partitioning in a typical RPLC system, a particle-based Monte Carlo simulation study is undertaken for a dimethyl octadecyl (C(18)) bonded stationary phase on a model siliceous substrate in contact with mobile phases having different methanol/water concentrations. Following upon previous simulations for gas-liquid chromatography and liquid-liquid phase equilibria, the simulations are conducted using the configurational-bias Monte Carlo method in the Gibbs ensemble and the transferable potentials for phase equilibria force field. The simulations are performed for a chain surface density of 2.9 micromol/m(2), which is a typical bonded-phase coverage for mono-functional alkyl silanes. The solvent concentrations used here are pure water, approximately 33 and 67% mole fraction of methanol and pure methanol. The simulations show that the chain conformation depends only weakly on the solvent composition. Most chains are conformationally disordered and tilt away from the substrate normal. The interfacial width increases with increasing methanol content and, for mixtures, the solvent shows an enhancement of the methanol concentration in a 10 Angstrom region outside the Gibbs dividing surface. Residual surface silanol groups are found to provide hydrogen bonding sites that lead to the formation of substrate bound water and methanol clusters, including bridging clusters that penetrate from the solvent/chain interfacial region all the way to the silica surface.

Molecular Dynamics Simulations of Alkylsilane Stationary-Phase Order and Disorder. 1. Effects of Surface Coverage and Bonding Chemistry

"Shape-selective" polymeric alkylsilane stationary phases are routinely employed over the more common monomeric phases in reversed-phase liquid chromatography (RPLC) to improve the separation of geometric isomers of shape-constrained solutes. We have investigated the molecular dynamics of chromatographic models that represent both monomeric and polymeric stationary phases with alkylsilane surface coverages and bonding chemistries typical of actual materials in an effort to elucidate the molecular-level structural features that control shape-selective separations. The structural characterization of these models is consistent with previous experimental observations of alkyl chain order and disorder: (1) alkyl chain order increases with increased surface coverage; and (2) monomeric and polymeric phases with similar surface coverages yield similar alkyl chain order (although subtle differences exist). In addition, a significant portion of the alkyl chain proximal to the silica surface is disordered (primarily gauche conformations) and the distal end is most ordered. Models that represent shape-selective RPLC phases possess a significant region of distal end chain order with primarily trans dihedral angle conformations. This is consistent with the view that the alkyl chains comprising polymeric stationary phases contain a series of well-defined and rigid voids in which shape-constrained solutes can penetrate and hence be selectively retained.

Molecular simulation study of the bonded-phase structure in reversed-phase liquid chromatography with neat aqueous solvent

The dramatic loss of retention in reversed-phase liquid chromatography when switching to 100% aqueous solvent and stopping flow (depressurizing) has long intrigued separation scientists. Recent experimental evidence suggests that the observed loss of retention is due to the loss of pore wetting with subsequent loss of solvent penetration in the porous matrix. One of the prevalent explanations of this phenomenon has been that the bonded phase chains, typically octadecyl silane bound to porous silica, would undergo significant conformational changes, viz. collapse, under pure aqueous conditions. As a definitive means toward elucidating the conformation of bonded-phase chains under pure aqueous conditions, configurational-bias Monte Carlo simulations in the Gibbs ensemble were carried out for a system of dimethyl octadecyl silane of intermediate coverage bound to the (111) face of beta-cristobalite and immersed in pure water. The results of two sets of simulations, which were started from two entirely different starting configurations as a validity check toward reaching the same equilibrium distribution of states, show that chains are neither clustering together nor laying on the surface but rather have a broad distribution of orientations and of conformational states. The interface between the bonded and solvent phases is rough on a molecular level, and clusters of water molecules are sometimes found to adsorb at the silica surface. This computational study lends further evidence that the driving force for the loss of retention when switching to pure aqueous conditions and depressurizing is not the collapse of bonded-phase chains.

Order and disorder in alkyl stationary phases

Covalently modified surfaces represent a unique state of matter that is not well described by liquid or solid phase models. The chemical bond in tethered alkanes imparts order to the surface in the form of anisotropic properties that are evident in chromatographic and spectroscopic studies. An understanding of the structure, conformation, and organization of alkyl-modified surfaces is requisite to the design of improved materials and the optimal utilization of existing materials. In recent years, the study of alkyl-modified surfaces has benefited from advances in modern analytical instrumentation. Aspects of alkyl chain conformation and motion have been investigated through the use of nuclear magnetic resonance spectroscopy, Fourier transform infrared spectroscopy, Raman spectroscopy, fluorescence spectroscopy, and neutron scattering studies. Chromatography provides complementary evidence of alkyl chain organization through interactions with solute probes. Computational simulations offer insights into the structure of covalently modified surfaces that may not be apparent through empirical observation. This manuscript reviews progress achieved in the study of the architecture of alkyl-modified surfaces.