Density Functional Theory Study on the Interactions of Metal Ions with Long Chain Deprotonated Carboxylic Acids.

In this work, interactions between carboxylate ions and calcium or sodium ions are investigated via density functional theory (DFT). Despite the ubiquitous presence of these interactions in natural and industrial chemical processes, few DFT studies on these systems exist in the literature. Special focus has been placed on determining the influence of the multibody interactions (with up to 4 carboxylates and one metal ion) on an effective pair-interaction potential, such as those used in molecular mechanics (MM). Specifically, DFT calculations are employed to quantify an effective pair-potential that implicitly includes multibody interactions to construct potential energy curves for carboxylate-metal ion pairs. The DFT calculated potential curves are compared to a widely used molecular mechanics force field (OPLS-AA). The calculations indicate that multibody effects do influence the energetic behavior of these ionic pairs and the extent of this influence is determined by a balance between (a) charge transfer from the carboxylate to the metal ions which stabilizes the complex and (b) repulsion between carboxylates, which destabilizes the complex. Additionally, the potential curves of the complexes with 1 and 2 carboxylates and one counterion have been examined to higher separation distance (20 Å) by the use of relaxed scan optimization and constrained density functional theory (CDFT). The results from the relaxed scan optimization indicate that near the equilibrium distance, the charge transfer between the metal ion and the deprotonated carboxylic acid group is significant and leads to non-negligible differences between the DFT and MM potential curves, especially for calcium. However, at longer separation distances the MM calculated interaction potential functions converge to those calculated with CDFT, effectively indicating the approximate domain of the separation distance coordinate where charge transfer between the ions is occurring.

Ab Initio Molecular Dynamics Study on the Interactions between Carboxylate Ions and Metal Ions in Water.

The interaction between a carboxylate anion (deprotonated propanoic acid) and the divalent Mg(2+), Ca(2+), Sr(2+), Ba(2+) metal ions is studied via ab initio molecular dynamics. The main focus of the study is the selectivity of the carboxylate-metal ion interaction in aqueous solution. The interaction is modeled by explicitly accounting for the solvent molecules on a DFT level. The hydration energies of the metal ions along with their diffusion and mobility coefficients are determined and a trend correlated with their ionic radius is found. Subsequently, a series of 16 constrained molecular dynamics simulations for every ion is performed, and the interaction free energy is obtained from thermodynamic integration of the forces between the metal ion and the carboxylate ion. The results indicate that the magnesium ion interacts most strongly with the carboxylate, followed by calcium, strontium, and barium. Because the interaction free energy is not enough to explain the selectivity of the reaction observed experimentally, more detailed analysis is performed on the simulation trajectories to understand the steric changes in the reaction complex during dissociation. The solvent dynamics appear to play an important role during the dissociation of the complex and also in the observed selectivity behavior of the divalent ions.

Heavy crude oils/particle stabilized emulsions.

Fluid characterization is a key technology for success in process design for crude oil mixtures in the future offshore. In the present article modern methods have been developed and optimized for crude oil applications. The focus is on destabilization processes in w/o emulsions, such as creaming/sedimentation and flocculation/coalescence. In our work, the separation technology was based on improvement of current devices to promote coalescence of the emulsified systems. Stabilizing properties based on particles was given special attention. A variety of particles like silica nanoparticles (AEROSIL®), asphalthenes, wax (paraffin) were used. The behavior of these particles and corresponding emulsion systems was determined by use of modern analytical equipment, such as SARA fractionation, NIR, electro-coalescers (determine critical electric field), Langmuir technique, pedant drop technique, TG-QCM, AFM.

Pore structural characteristics, size exclusion properties and column performance of two mesoporous amorphous silicas and their pseudomorphically transformed MCM-41 type derivatives.

Highly ordered mesoporous silicas such as, mobile composition of matter, MCM-41, MCM-48, and the SBA-types of materials have helped to a large extent to understand the formation mechanisms of the pore structure of adsorbents and to improve the methods of pore structural characterization. It still remains an open question whether the high order, the regularity of the pore system, and the narrow pore size distribution of the materials will lead to a substantial benefit when these materials are employed in liquid phase separation processes. MCM-41 type 10 microm beads are synthesized following the route of pseudomorphic transformation of highly porous amorphous silicas. Highly porous silicas and the pseudomorphically transformed derivatives are characterized by nitrogen sorption at 77 K and by inverse size-exclusion chromatography (ISEC) employing polystyrene standards. Applying the network model developed by Grimes, we calculated the pore connectivity n(T) of the materials. The value of n(T) varies between the percolation threshold of the lattice and values of n(T) > 10, the latter being the limiting value above which the material can be considered to be almost infinitely connected such that the ISEC behavior of the material calculated with the pore network model is the same when calculated with a parallel pore model which assumes an infinite connectivity. One should expect that the pore connectivity is reflected in the column performance, when these native and unmodified materials are packed into columns and tested with low molecular weight analytes in the Normal Phase LC mode. As found in a previous study on monolithic silicas and highly porous silicas, the slope of the plate height (HETP) - linear velocity (u) curve decreased significantly with enhanced pore connectivity of the materials. First results on the pseudomorphically transformed MCM-41 type silicas and their highly porous amorphous precursors showed that (i) the transformation did not change the pore connectivity (within the limits detectable by ISEC) from the starting material to the final product and (ii) the slope of the HETP versus u curve for dibutylphtalate did not change significantly after the pseudomorphic transformation.

The effect of the pore structure and zeta potential of porous polymer monoliths on separation performance in ion-exchange mode.

Most often, in bioseparations involving charged macromolecules, the chromatographic systems have low Reynolds and high Peclet numbers. For such systems, an expression is developed and presented in this work for evaluating the throughput in polymeric monoliths where ion-exchange adsorption occurs, as a function of (i) the pressure drop along the length of the monolith, (ii) the functional form and width of the throughpore-size distribution of the monolith, and (iii) the magnitude of the zeta potential on the surface of the throughpores of the monolith. Gaussian and log-normal throughpore-size distributions whose mean throughpore-size and standard deviation values are based on experimentally measured throughpore-size distribution data by mercury porosimetry employed on polymeric monoliths are used in this work, and their effect on the throughput relative to that obtained from a polymeric monolith having a uniform throughpore-size distribution is studied for different values of the ratio of the standard deviation to the mean throughpore-size. The results indicate that relatively modest increases in the throughput, when compared with the throughput that could be achieved in a polymeric monolith having a uniform throughpore-size distribution, could be obtained in polymeric monoliths having disperse throughpore-size distributions, and the magnitude of the increase becomes larger when the disperse distribution is skewed to larger throughpore sizes. Furthermore, the results of this work indicate that, under certain conditions, relatively modest increases in the throughput of a charged analyte could also be achieved by altering the value of the zeta potential on the surface of the throughpores of the monolith. Due to the difficulties inherent in controlling the functional form and width of the throughpore-size distribution during the synthesis of polymeric monoliths, it would appear to be more practical to increase the value of the throughput of a charged analyte by altering the value of the zeta potential through prudent selection of the ion-exchange surface functional groups and fine-tuned with the pH of the mobile phase. Thus, for ion-exchange chromatography systems, the zeta potential could be considered an important parameter for column designers and operators to manipulate, since its alteration could increase the through-put of a charged analyte in polymeric monoliths or in columns packed with charged particles.

Impact of pore structural parameters on column performance and resolution of reversed-phase monolithic silica columns for peptides and proteins.

In this work, monolithic silica columns with the C4, C8, and C18 chemistry and having various macropore diameters and two different mesopore diameters are studied to access the differences in the column efficiency under isocratic elution conditions and the resolution of selected peptide pairs under reversed-phase gradient elution conditions for the separation of peptides and proteins. The columns with the pore structural characteristics that provided the most efficient separations are then employed to optimize the conditions of a gradient separation of a model mixture of peptides and proteins based on surface chemistry, gradient time, volumetric flow rate, and acetonitrile concentration. Both the mesopore and macropore diameters of the monolithic column are decisive for the column efficiency. As the diameter of the through-pores decreases, the column efficiency increases. The large set of mesopores studied with a nominal diameter of approximately 25 nm provided the most efficient column performance. The efficiency of the monolithic silica columns increase with decreasing n-alkyl chain length in the sequence of C18

Pore structural characterization of monolithic silica columns by inverse size-exclusion chromatography.

In this work, a parallel pore model (PPM) and a pore network model (PNM) are developed to provide a state-of-art method for the calculation of several characteristic pore structural parameters from inverse size-exclusion chromatography (ISEC) experiments. The proposed PPM and PNM could be applicable to both monoliths and columns packed with porous particles. The PPM and PNM proposed in this work are able to predict the existence of the second inflection point in the experimental exclusion curve that has been observed for monolithic materials by accounting for volume partitioning of the polymer standards in the macropores of the column. The appearance and prominence of the second inflection point in the exclusion curve is determined to depend strongly on the void fraction of the macropores (flow-through pores), (b) the nominal diameter of the macropores, and (c) the radius of gyration of the largest polymer standard employed in the determination of the experimental ISEC exclusion curve. The conditions that dictate the appearance and prominence of the second inflection point in the exclusion curve are presented. The proposed models are applied to experimentally measured ISEC exclusion curves of six silica monoliths having different macropore and mesopore diameters. The PPM and PNM proposed in this work are able to determine the void fractions of the macropores and silica skeleton, the pore connectivity of the mesopores, as well as the pore number distribution (PND) and pore volume distribution (PVD) of the mesopores. The results indicate that the mesoporous structure of all materials studied is well connected as evidenced by the similarities between the PVDs calculated with the PPM and the PNM, and by the high pore connectivity values obtained from the PNM. Due to the fact that the proposed models can predict the existence of the second inflection point in the exclusion curves, the proposed models could be more applicable than other models for ISEC characterization of chromatographic columns with small diameter macropores (interstitial pores) and/or large macropore (interstitial pore) void fractions. It should be noted that the PNM can always be applied without the use of the PPM, since the PPM is an idealization that considers an infinitely connected porous medium and for materials having a low (<6) pore connectivity the PPM would force the PVD to a lower average diameter and larger distribution width as opposed to properly accounting for the network effects present in the real porous medium.