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

Investigating hydrophobic environment in alkyl-group-functionalized silica particle with various chain lengths using absorption microspectroscopy

A well-known solvatochromic dye, Reichardt's dye (R-dye), was used to evaluate the hydrophobicity of alkyl-group-functionalized silica particles (ASPs) with different chain lengths. The absorption spectra of R-dye were measured in a single ASP in a mixed solution of water and an organic solvent (methanol (MeOH), ethanol (EtOH), acetonitrile (ACN), tetrahydrofuran (THF), or N,N-dimethylformamide (DMF)) using absorption microspectroscopy. The polarity parameter in the ASPs (ET), determined by the absorption maximum, was observed to be smaller than those in bulk solutions, indicating that R-dye was present in a more hydrophobic environment. In EtOH, THF, and DMF, R-dye was distributed within the alkyl chain layer including the organic solvent. An increase in the organic solvent content of the bulk solution led to a higher organic solvent concentration in the alkyl chain layer, resulting in a decrease in ET. In MeOH and ACN, the R-dye was distributed within the alkyl chain layer and concentrated phase. Moreover, with the increase in the organic molecule content, the distribution of R-dye in the concentrated phase became dominant in MeOH and ACN system, leading to an increase in the ET value. The findings presented in this paper are expected to attract the attention of a wide range of researchers in chromatography.

Separations by Shape: Molecular Shape Recognition in Liquid Chromatography

Molecular shape can provide a basis for chromatographic separations that is distinct from other interaction mechanisms often considered in liquid chromatography. Molecular shape recognition, or shape selectivity, is most evident for the separation of isomeric compounds that have constrained molecular structures, such as polycyclic aromatic hydrocarbons. A central feature of shape-selective columns is conformational order within the stationary phase; this aspect of stationary phase morphology is revealed through spectroscopic studies and molecular dynamics simulations, and is correlated with chromatographic performance.

A carbonylative coupling approach to alkyl stationary phases with variable embedded carbamate groups for high-performance liquid chromatography

A highly efficient carbonylative coupling method for the preparation of alkyl stationary phases with variable numbers of carbamate groups was established. The effectiveness of such method was verified through elemental analysis, X-ray photoelectron spectroscopy and solid-state ¹³C nuclear magnetic resonance spectroscopy of three as-synthesized stationary phases bearing different alkyl chains and different numbers of carbamate groups (octadecyl/one carbamate group, C18C; docosyl/three carbamate groups, C22C³; triacontyl/two carbamate groups, C30C²). The comparative evaluation of these stationary phases using a great variety of analytes, including three sets of isomers of alkylbenzenes, two sets of standard mixtures of polycyclic aromatic hydrocarbons (SRM 1647e and 869b), nine polychlorinated biphenyls, fiveteen N-substituted ureas, ten sulfonylureas, five xanthines and some other phytonutrients, revealed their remarkable applicability in reversed-phase liquid chromatography. Notably, the intercalated carbamate groups rendered the resultant stationary phases compatible with 100% aqueous mobile phase. The suppression of silanol activity was positively related to the number of polar groups embedded in the bonded selector, and the smallest peak tailing factor (1.14) for amitriptyline was obtained by C22C³. The molecular shape-related selectivity was found to be more closely related to the length of the selector's aliphatic chain, as supported by the lowest α_TBN/BaP value (0.31) by C30C². These carbamate-embedded alkyl stationary phases constituted another class of polar-embedded stationary phases possessing a single type of functional ligand.

Interfacial study of clathrates confined in reversed silica pores

Storing methane in clathrates is one of the most promising alternatives for transporting natural gas (NG) as it offers similar gas densities to liquefied and compressed NG while offering lower safety risks. However, the practical use of clathrates is limited given the extremely low temperatures and high pressures necessary to form these structures. Therefore, it has been suggested to confine clathrates in nanoporous materials, as this can facilitate clathrate's formation conditions while preserving its CH4 volumetric storage. Yet, the choice of nanoporous materials to be employed as the clathrate growing platform is still rather arbitrary. Herein, we tackle this challenge in a systematic way by computationally exploring the stability of clathrates confined in alkyl-grafted silica materials with different pore sizes, ligand densities and ligand types. Based on our findings, we are able to propose key design criteria for nanoporous materials favoring the stability of a neighbouring clathrate phase, namely large pore sizes, high ligand densities, and smooth pore walls. We hope that the atomistic insight provided in this work will guide and facilitate the development of new nanomaterials designed to promote the formation of clathrates.

Distinction of synthetic dl-α-tocopherol from natural vitamin E (d-α-tocopherol) by reversed-phase liquid chromatography. Enhanced selectivity of a polymeric C18 stationary phase at low temperature and/or at high pressure

Separation of diastereomers of dl-α-tocopherol was studied by reversed-phase liquid chromatography using three types of stationary phases, polymeric ODS, polymeric C30, and monomeric ODS. Polymeric ODS stationary phase (Inertsil ODS-P, 3mmID, 20cm) was effective for the separation of the isomers created by the presence of three chiral centers on the alkyl chain of synthetic dl-α-tocopherol. Considerable improvement of the separation of isomers was observed on ODS-P phase at high pressure and at low temperature. Complete separation of four pairs of diastereomers was achieved at 12.0°C, 536bar, while three peaks were observed when the separation was carried out either at 12.0°C at low pressure or at 20°C at 488bar. Higher temperature (30.0°C) with the ODS-P phase resulted in only partial separation of the diastereomers even at high pressure. Only slight resolution was observed for the mixture of diastereomers with the C30 stationary phase (Inertsil C30) at 12.0°C and 441bar, although the stationary phase afforded greater resolution for β- and γ-tocopherol than ODS-P. A monomeric C18 stationary phase did not show any separation at 12.0°C and 463bar. The results suggest that the binding site of the polymeric ODS-P phase is selective for flexible alkyl chains that provided the longest retention for the natural form, (R,R,R) form, and the enantiomer, (S,S,S) form, of dl-α-tocopherol.

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.

Molecular simulation studies of reversed-phase liquid chromatography

Over the past 20 years, molecular simulation methods have been applied to the modeling of reversed-phase liquid chromatography (RPLC). The purpose of these simulations was to provide a molecular-level understanding of: (i) the structure and dynamics of the bonded phase and its interface with the mobile phase, (ii) the interactions of analytes with the bonded phase, and (iii) the retention mechanism for different analytes. However, the investigation of chromatographic systems poses significant challenges for simulations with respect to the accuracy of the molecular mechanics force fields and the efficiency of the sampling algorithms. This review discusses a number of aspects concerning molecular simulation studies of RPLC systems including the historical development of the subject, the background needed to understand the two prevalent techniques, molecular dynamics (MD) and Monte Carlo (MC) methods, and the wealth of insight provided by these simulations. Examples from the literature employing MD approaches and from the authors' laboratory using MC methods are discussed. The former can provide information on chain dynamics and transport properties, whereas the latter techniques are uniquely suited for the investigation of phase and sorption equilibria that underly RPLC retention, and both can be used to elucidate the bonded-chain conformations and solvent distributions.

Ionic polymeric amphiphiles with cholesterol mesogen: adsorption and organization characteristics at the air/water interface from Langmuir film balance studies

Ionic polymeric amphiphiles consisting of cholesterol mesogen were investigated for the interfacial adsorption characteristics at the air/water interface using a Langmuir film balance with an aim to understand the influence of ionic segment from 2-acrylamido-2-methyl-1-propane sulfonic acid (AMPS) on the packing behavior of cholesterol at the interface. From surface pressure (pi)-area (A) isotherm characteristics, it is demonstrated that the homopolymer and the copolymer C consisting of 0.15 mol fraction CAB segments exhibit the most expanded structures contributing to surface area of about 84 A(2)/molecule. It is shown that the copolymer B with 0.1 mol fraction CAB provides optimum hydrophilic liphophilic balance to form the most compact structures contributing to a surface area of 35.75 A(2)/molecule. The high surface pressure, >40 mN/m, in contrast to that of PAMPS demonstrates significant adsorption of the copolymers at the interface. An interesting correlation among interfacial packing characteristics, thermal behavior, and solution structures is demonstrated. From molecular models developed for CAB, it is shown that the horizontal orientation of the linker group with respect to cholesterol chain in CAB underlies the expanded structures observed in PCAB and copolymer C.

Effect of the surface coverage of endcapped C18-silica on the excess adsorption isotherms of commonly used organic solvents from water in reversed phase liquid chromatography

The excess adsorption isotherms of methanol, ethanol, 2-propanol, acetonitrile, and tetrahydrofuran from water were measured on five different silica-based packing materials by the minor disturbance method. These materials were prepared with the same lot of 5-microm particles (average pore size 90 A), all endcapped with trimethylchlorosilane (TMS), and bonded to octadecyl chains with different surface coverages (0, 0.42, 1.01, 2.03, and 3.15 micromol/m2). The relative adsorption of one eluent by respect to a second one informs on the heterogeneity of the material (alkyl-bonded and bare silica regions) and on the accessibility of the unreacted silanol groups to the mobile phase. It is shown that the total surface area of the adsorbent decreases with increasing degree of surface coverage with octadecyl chains and that the relative surface area of the regions occupied by accessible silanol groups to the regions occupied by alkyl-bonded groups decreases. For the five columns, an average of 10% of the adsorbent surface area is covered of bare silica accessible to the liquid phase, with a minimum of 5% with tetrahydrofuran and a maximum of 12% with ethanol or 2-propanol. Increasing the surface coverage by the C18 chains causes a significant increase of the attraction potential of the hydrophobic surface toward the organic solvent. This result is confirmed by the increase of the number of adsorbate monolayers with increasing bonding density of the octadecyl chains. This number is twice larger for the 315C18 column than for the C1 column.

Adsorption of Acridine Orange at a C8,18/Water/Acetonitrile Interface

Fully atomistic simulations are used to characterize the molecular dynamics (MD) of acridine orange (3,6-dimethylaminoacridine) at a chromatographic interface. Multiple 1 ns MD simulations were performed for acridine orange at the interface between three different acetonitrile/water mixtures (0/100, 20/80, and 50/50) with C8 and C18 alkyl chains. The diffusion coefficient, D, of acridine orange in pure solvent was found to be 4 times smaller at the water/C18 interface (D = 0.022 x 10(-4) cm2/s) than in bulk water (D = 0.087 x 10(-4) cm2/s), in qualitative agreement with experiment. Rotational reorientation times were 20 and 700 ps, which also agree favorably with the measured time scales of 130 and 740 ps. Contrary to experiment, the simulations found that for increasing surface coverage, the diffusion coefficient for acridine decreased. Detailed analysis of the solvent structure showed that the transport properties of acridine were primarily governed by the solvent distribution above the functionalized surface. The solvent structure, in turn, was largely determined by the surface consisting of the silica layer, the alkyl chains, and their functionalization.

Stationary phase characterization and method development

The properties of stationary phases and their characterization methods are reviewed. New and significant developments have occurred in the last few years, and new methods for stationary phase characterization have become available. The characterization methods are discussed, and the differences between the different methods are pointed out. In addition, method development approaches are reviewed, with special emphasis on recent developments that employ multiple parameters in parallel. Also, the renewed interest of temperature as a tool in method development is surveyed.

Atomistic simulations to compute surface properties of poly(N-vinyl-2-pyrrolidone) (PVP) and blends of PVP/chitosan

Atomistic simulations were performed on poly(N-vinyl-2-pyrrolidone) (PVP) and its blends with chitosan (CS) in different ratios using molecular mechanics (MM) and molecular dynamics (MD) simulations in three-dimensionally periodic and effective two-dimensionally periodic condensed phases. Four independent microstructures were generated to analyze their surface properties. The calculated surface-energy values for PVP compared quite well with the experimental data reported in the literature. The density profile was analyzed, and the structure of the films showed an interior region of the bulk density. Various components of the energetic interactions (torsional, van der Waals, etc.) were examined to gain deeper insight into the nature of regular and anomalous interactions between the bulk and the surface films. Surface energies of PVP/CS blends were computed by MD simulations using the bulk pressure-volume-temperature (PVT) parameters. Bulk properties such as the cohesive energy density (CED) and solubility parameter (delta) were calculated using MM and MD simulations in the NVT ensemble under periodic boundary conditions. The Flory equation of state was used to compute the thermal expansion coefficient as well as PVT parameters. These surface-energy values agreed well with the surface-energy data calculated using the Zisman equation, which were also in accordance with the experimental observations. The results from this study suggest that computer simulations would provide valuable information on polymers and polymer-blend surfaces.

Experimental evidence of the influence of the surface chemistry of the packing material on the column pressure drop in reverse-phase liquid chromatography

The permeabilities of six columns packed with different packing materials (neat silica, C(1) endcapped silica at 3.92 micro mol/m(2), C(18) bonded and endcapped silica with 0.42, 1.01, 2.03, and 3.15 micro mol/m(2) of C(18) bonded chains) were measured. All these materials derive from the same batch of spherical particles, 5 micro m in diameter. The columns have the same tube inner diameter (phi=0.460+/-0.003 cm) and length (L=15.000+/-0.003 cm). The experimental conditions were the same, flow-rate (F(v)=1.000+/-0.003 mL/min) and temperature (295 K). Nevertheless, it was found that the column permeability decreases significantly, by about 25%, from the neat silica column to the one packed with the highest density of C(18)-bonded silica (3.15 micro mol/m(2)). The results measured on two duplicate columns were very reproducible. Accurate (+/-0.5 %) measurements of the hold-volumes with concentrated and dilute solutions of NO(3)(-) showed that the columns had all nearly the same external porosity. The result cannot be explained by the error made on the volume of the column tube either as it was measured accurately for all the columns. The residual explanation is that the interstitial velocity distribution between the packed particles depends on the chemical nature of the external surface of these particles.

Selectivity in reversed-phase separations

The selectivity difference between 15 different stationary phases was measured using a large number of analytes at 2 or 3 different pH values (3, 7 and 10) with acetonitrile and methanol as the mobile phase modifiers. The packings discussed include standard C(8) and C(18) packings, packings with embedded polar groups, a phenyl packing, a pentafluoro-phenyl packing, an adamantylethyl packing and others. The major selectivity differences observed are discussed in detail. Specific effects such as pi-pi interactions on phenyl packings or hydrogen-bond interactions on phases with embedded polar groups are confirmed.

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.

Identification of isolated cavity features within molecular dynamics simulated chromatographic surfaces

Highly ordered morphological features were characterized for molecular dynamics simulated alkyl-modified silica models that represent chromatographic materials with enhanced shape recognition capability. Deep cavities (8-10A wide) within the alkyl chains were identified for C18 polymeric models corresponding to shape-selective RPLC stationary phases. The all-trans conformational distal-end segments of these isolated cavities averaged over a 100 ps simulation time interval were observed to increase (up to 15 A) in models with an increase in both surface coverage and corresponding shape selectivity. Similar-structure cavities with significant alkyl chain ordered regions (>11A) were isolated from two independent C18 models (differing in bonding chemistry, density and temperature) that represent highly shape-selective materials. The size and depth of these ordered regions increased (up to 28 A) for the extended-length C30 alkyl phase models. These initial results offer a physical representation of alkyl-modified surfaces that may facilitate the identification of potential molecular features that may be involved in the shape-selective retentive processes, as well as illustrating the potential for such computational techniques to predict the molecular recognition capabilities of novel analyte-specific sorbents.

Heterogeneity of the Adsorption Mechanism of Low Molecular Weight Compounds in Reversed-Phase Liquid Chromatography

The retention mechanism in RPLC mode was investigated based on the acquisition of adsorption isotherm data by frontal analysis measurements and their modeling. This work is a review of the results of four years of adsorption data measurements. The data were acquired on a wide variety of brands of C18-silica columns (from Akzo Nobel, Bishoff, Hypersil, Merck, Phenomenex, Supelco, Vydac, and Waters) with several low molecular weight compounds such as phenol (94 g/mol), caffeine (194 g/mol), tryptophan (204 g/mol), sodium 2-naphthalenesulfonate (235 g/mol), and propranololium chloride (295 g/mol). The mobile phase was a mixture of methanol and water at variable composition. The adsorption isotherms were all convex upward (langmuirian), and the degree of heterogeneity of the adsorption system was determined from the calculation of the adsorption energy distribution using the expectation-maximization method. The adsorption isotherm parameters (number of types of adsorption sites, surface concentration of each type of site, and difference between the adsorption energies E(i) - E(j) on sites i and j), obtained from the mathematical fit of the adsorption data to the appropriate multi-Langmuir adsorption isotherm model, were analyzed and compared. The results allow the drawing of general conclusions regarding the relationships between the size of the analyte and the adsorption properties (saturation capacities, adsorption energies) characterizing the retention mechanism in RPLC mode for neutral, anionic, and cationic compounds.

Adsorption mechanism in reversed-phase liquid chromatography

The effect of the bonding density of the octadecyl chains onto the same silica on the adsorption and retention properties of low molecular weight compounds (phenol, caffeine, and sodium 2-naphthalene sulfonate) was investigated. The same mobile phase (methanol:water, 20:80, v/v) and temperature (T = 298 K) were applied and two duplicate columns (A and B) from each batch of packing material (neat silica, simply endcapped or C1 phase, 0.42, 1.01, 2.03, and 3.15 micromol/m2 of C18 alkyl chains) were tested. Adsorption data of the three compounds were acquired by frontal analysis (FA) and the adsorption energy distributions (AEDs) were calculated using the expectation-maximization method. Results confirmed earlier findings in linear chromatography of a retention maximum at an intermediate bonding density. From a general point of view, the saturation capacity of the adsorbent tends to decrease with increasing bonding density, due to the vanishing space intercalated between the C18 bonded chains and to the decrease of the specific surface area of the stationary phase. The equilibrium constants are maximum for an intermediary bonding density (approximately 2 micromol/m2). An enthalpy-entropy compensation was found for the thermodynamic parameters of the isotherm data. Weak equilibrium constants (small deltaH) and high saturation capacities (large deltaS) were observed at low bonding densities, higher equilibrium constants and lower saturation capacities at high bonding densities, the combinations leading to similar apparent retention in RPLC. The use of a low surface coverage column is recommended for preparative purposes.