A computational method to characterize framework aluminum in aluminosilicates.
Explore the source record for details and available documents.
Biomedical subjects
Publications and source records attributed to Berend Smit.
Explore the source record for details and available documents.
The recently proposed united atom force field by Dubbeldam et al. (Phys. Rev. Lett. 2004, 93, 088302) for the adsorption of alkanes in MFI-type zeolites was extended to other zeolites in this work. Its applicability to FER-type zeolites was evaluated in detail, for which the Henry coefficients, the isosteric heat of adsorption, the adsorption isotherms, as well as the locations of alkanes in the FER-type zeolites were computed and compared to experimental values. The results show that the new force field works well for FER zeolites. Furthermore, its applicability to MWW-, MTW-, CFI-, LTA-, and STF-type zeolites was investigated, and we found that the experimental isotherms could be accurately predicted except for STF-type zeolites. This work shows that the new united atom force field proposed by Dubbeldam et al. is applicable to most pure silica zeolites.
Disorder of intercalated surfactant molecules in clay minerals causes gradual swelling, rather than commonly assumed swelling in discrete steps.
The applicability of a recently proposed force field of Calero et al. (J. Am. Chem. Soc., 2004, 126, 11377) to Na-MOR zeolites is evaluated. The Henry law coefficients of ethane and C(5)-C(9) as well as the adsorption isotherms of ethane, propane, butane, and hexane in various Na-MOR zeolites are computed and compared with experimental values. These comparisons show that the new force field is suitable for Na-MOR zeolites. Furthermore, this force field is used to study the effects of sodium cations on the adsorption behavior of larger alkanes, such as C(4)-C(7), in MOR-type zeolites. These simulations give a better understanding of the underlying mechanisms of the cations' position and density influence on adsorption. In addition, a characteristic pressure named "reversal pressure" is introduced which characterizes the efficiency of the presence of cations in zeolites.
We perform grand-canonical molecular simulations to study the molecular mechanism of clay swelling hysteresis as a function of the relative humidity. In particular, we focus on the transition from the one- to the two-layer hydrate and the influence of three types of counterions (Li+, Na+, and K+). Our results cover the experimental relative humidity region where swelling and shrinking usually occur. We show that the thermodynamic origin of swelling hysteresis is a free-energy barrier separating the layered hydrates. This free-energy barrier is dominated by breaking and formation of hydrogen bonds between and within water layers. This network of water molecules is similar for all counterions, but the positions of these counterions depend upon their size. The relatively large K+ counterions show more affinity for clay surface adsorption, which increases the free-energy barrier and inhibits swelling. On the other hand, the relatively small Li+ counterions are quite well-accommodated in the water network, and thereby, they can form a new swelling state with a basal spacing of approximately 13.5 A. This new swelling state is an alternative explanation for the widely accepted simultaneous occurrence of two or more swelling phases.
A new force field for the simulation of the adsorption of cycloalkanes in nanoporous silica affords a significant improvement over any previously employed force field. The simulated isotherms reproduce the most salient features in the experimental isotherms extremely well. The study of cyclo-pentane, -hexane, and -heptane adsorption in MFI-type silica indicates an inflection for cyclopentane but not for cyclohexane at intermediate pressure. If corroborated by experiments, such an inflection point would afford an excellent calibration point for further force field developments. At low pressures, mixture isotherms of cyclohexane and n-hexane show a temperature dependence on the selectivity in accordance with recent results by J. P. Fox and S. P. Bates, J. Phys. Chem., 2004, 108, 17136. This dependence is caused by a difference in temperature dependence of the Henry coefficient for both molecules. At high pressures entropy effects due to packing always favor the sorption of n-hexane. Furthermore, the influence of the flexibility of the zeolite framework on the adsorption of these rather bulky molecules is investigated. It is found that this influence of the flexibility on the adsorption of cyclohexane is as small as with n-alkanes.
We investigated the phase behavior of double-tail lipids, as a function of temperature, headgroup interaction and tail length. At low values of the head-head repulsion parameter a(hh), the bilayer undergoes with increasing temperature the transitions from the subgel phase L(c) via the flat gel phase L(beta) to the fluid phase L(alpha). For higher values of a(hh), the transition from the L(c) to the L(alpha) phase occurs via the tilted gel phase L(beta)(') and the rippled phase P(beta)('). The occurrence of the L(beta)(') phase depends on tail length. We find that the rippled structure (P(beta)(')) occurs if the headgroups are sufficiently surrounded by water and that the ripple is a coexistence between the L(c) or L(beta)(') phase and the L(alpha) phase. The anomalous swelling, observed at the P(beta)(') --> L(alpha) transition, is not directly related to the rippled phase, but a consequence of conformational changes of the tails.
Biological membranes are complex and highly cooperative structures. To relate biomembrane structure to their biological function it is often necessary to consider simpler systems. Lipid bilayers composed of one or two lipid species, and with embedded proteins, provide a model system for biological membranes. Here we present a mesoscopic model for lipid bilayers with embedded proteins, which we have studied with the help of the dissipative particle dynamics simulation technique. Because hydrophobic matching is believed to be one of the main physical mechanisms regulating lipid-protein interactions in membranes, we considered proteins of different hydrophobic length (as well as different sizes). We studied the cooperative behavior of the lipid-protein system at mesoscopic time- and lengthscales. In particular, we correlated in a systematic way the protein-induced bilayer perturbation, and the lipid-induced protein tilt, with the hydrophobic mismatch (positive and negative) between the protein hydrophobic length and the pure lipid bilayer hydrophobic thickness. The protein-induced bilayer perturbation was quantified in terms of a coherence length, xi(P), of the lipid bilayer hydrophobic thickness profile around the protein. The dependence on temperature of xi(P), and the protein tilt-angle, were studied above the main-transition temperature of the pure system, i.e., in the fluid phase. We found that xi(P) depends on mismatch, i.e., the higher the mismatch is, the longer xi(P) becomes, at least for positive values of mismatch; a dependence on the protein size appears as well. In the case of large model proteins experiencing extreme mismatch conditions, in the region next to the so-called lipid annulus, there appears an undershooting (or overshooting) region where the bilayer hydrophobic thickness is locally lower (or higher) than in the unperturbed bilayer, depending on whether the protein hydrophobic length is longer (or shorter) than the pure lipid bilayer hydrophobic thickness. Proteins may tilt when embedded in a too-thin bilayer. Our simulation data suggest that, when the embedded protein has a small size, the main mechanism to compensate for a large hydrophobic mismatch is the tilt, whereas large proteins react to negative mismatch by causing an increase of the hydrophobic thickness of the nearby bilayer. Furthermore, for the case of small, peptidelike proteins, we found the same type of functional dependence of the protein tilt-angle on mismatch, as was recently detected by fluorescence spectroscopy measurements.
We use dissipative particle dynamics (DPD) and molecular models to simulate interacting oil/water/surfactant interfaces. The system comprises sections of two emulsion droplets separated by a film. The film is in equilibrium with a continuous phase, in analogy with the surface force apparatus. This is achieved by combining DPD with a Monte Carlo scheme to simulate a muVT ensemble. The setup enables the computation of surface forces as a function of the distance between the two interfaces, as well as the detection of film rupture. We studied monolayers of nonionic model surfactants at different densities and compared oil-water-oil and water-oil-water emulsion films. Between surfactant monolayers facing each other tails-on (water-oil-water films), we observed repulsive forces due to the steric interaction between overlapping hydrophobic tails. The repulsion increases with surfactant density. Conversely, no such repulsion is observed between surfactant monolayers facing each other heads-on. Instead, the film ruptures, the monolayers merge, and a channel forms between the two droplet phases. Film rupture can also be induced in the water-oil-water films by forcing the interfaces together. The separation at rupture increases for oil-water-oil films and decreases for water-oil-water films when the surfactant density increases. The results are in qualitative agreement with existing theories of emulsion stability in creams, in particular with the channel nucleation theory based on the natural curvature of surfactants.
We have developed a united atom force field able to accurately describe the adsorption properties of linear alkanes in the sodium form of FAU-type zeolites. This force field successfully reproduces experimental adsorption properties of n-alkanes over a wide range of sodium cation densities, temperatures, and pressures. The force field reproduces the sodium positions in dehydrated FAU-type zeolites known from crystallography, and it predicts how the sodium cations redistribute when n-alkanes adsorb. The cations in the sodalite cages are significantly more sensitive to the n-alkane loading than those in the supercages. We provide a simple expression that adequately describes the n-alkane Henry coefficient and adsorption enthalpy as a function of sodium density and temperature at low coverage. This expression affords an adequate substitute for complex configurational-bias Monte Carlo simulations. The applicability of the force field is by no means limited to low pressure and pure adsorbates, for it also successfully reproduces the adsorption from binary mixtures at high pressure.
Adsorption of alcohol molecules or other small amphiphilic molecules in the cell membrane can induce significant changes in the structure of the membrane. To understand the molecular mechanisms underlying these structural changes, we developed a mesoscopic membrane model. Molecular simulations on this model nicely reproduce the experimental phase diagrams. We find that alcohol can induce an interdigitated structure in which the normal bilayer structure changes into a monolayer in which the alcohol molecules screen the hydrophobic tails from the water phase. We compute the effect of the chain length of the alcohol on the phase behaviour of the membrane. At low concentrations of alcohol, the membrane has domains of the interdigitated phase that are in coexistence with the normal membrane phase. We use our model to clarify some of the experimental questions related to the structure of the interdigitated phase and put forward a simple model that explains the alcohol chain length dependence of the stability of this interdigitated phase.
The effect of the surfactant chain length n on the bending modulus kappa of surfactant monolayers is simulated with a mesoscopic oil-water-surfactant model. We confirm a power law, kappa is proportional to np, as predicted by mean-field theory and found experimentally, and find p approximately 1.5 at a constant surface density and p approximately 1.0 at a constant interfacial tension. This agrees quite well with both mean-field theory (p=2-3, assuming constant surface density) and experiments (at constant surface tension). Our results suggest that the previously reported agreement between theory and experiment may be fortuitous and caused by the difference in surfactant types.
We have used dissipative particle dynamics to simulate amphiphilic monolayers on the interface between oil and water. An ultralow interfacial tension is imposed by means of Monte Carlo to resemble the amphiphilic films that separate oil and water regions in microemulsions. We calculate the bending modulus by analyzing the undulation spectrum. By varying the surfactant chain length and topology we investigate the effect of surfactant structure and composition of the monolayer on the bending moduli. We find that increasing the thickness has a larger effect than increasing the density of the layer. This follows from the observations that at a given interfacial tension, the bending modulus increases with chain length and is larger for linear than branched surfactants. The increase with chain length is approximately linear, which is slower than the theoretical predictions at a fixed area. We also investigated a binary mixture of short and long surfactants compared to pure layers of the same average chain length. We find a roughly linear decrease in bending modulus with mole fraction of short surfactants. Furthermore, the mixed film has a lower bending modulus than the corresponding pure film for all mole fractions. Linking the bending moduli to the structure of the surfactants is an important step in predicting the stability of microemulsions.
In this study we introduce a mesoscopic lipid-water-alcohol model. Dissipative particle dynamics (DPD) simulations have been used to investigate the induced interdigitation of bilayers consisting of double-tail lipids by adding alcohol molecules to the bilayer. Our simulations nicely reproduce the experimental phase diagrams. We find that alcohol can induce an interdigitated structure where the common bilayer structure changes into monolayer in which the alcohol molecules screen the hydrophobic tails from the water phase. At low concentrations of alcohol the membrane has domains of the interdigitated phase that are in coexistence with the common membrane phase. We compute the effect of the chain length of the alcohol on the phase behavior of the membrane and show that the stability of the interdigitated phase depends on the length of the alcohol. We show that we can reproduce the experimental hydrophobic thickness of the bilayer for various combinations of lipids and alcohols. We use our model to clarify some of the experimental questions related to the structure of the interdigitated phase and put forward a simple model that explains the alcohol chain length dependence of the stability of this interdigitated phase.
Explore the source record for details and available documents.
Dissipative particle dynamics simulations are used to study the self-assembly of lipids into bilayers. With a simple mesoscopic lipid-water model, we observe the formation of the liquid crystalline phase L(alpha) and gel phases in which the tails are interdigitated L(betaI) or noninterdigitated L(beta). For double-tail lipids experiments show all three phases, while for single-tail lipids only L(beta) and L(alpha) are observed. We show that at sufficiently high head-head repulsion the L(betaI) is stable for single-tail lipids. This suggests that it might be possible to induce an L(beta)-->L(betaI) transition by adding chaotropic salts.
Explore the source record for details and available documents.
Recent developments in Configurational-Bias Monte Carlo (CBMC) techniques allow the accurate calculation of the sorption isotherms for alkanes, and their mixtures, in various zeolites. The CBMC simulations give new insights into subtle entropy effects affecting mixture adsorption. Three types of entropy effects can be distinguished. (1) Size entropy effects favour the component with the smaller number of C atoms because the smaller molecule finds it easier to fill in the 'gaps' within the zeolite matrix at high molecular loadings. (2) Configurational entropy effects come into play for mixtures of alkanes that differ in the degree of branching. For a mixture of linear and branched alkanes with the same number of C atoms, configurational entropy effects favour the linear isomer because such molecules 'pack' more efficiently within, say, the intersecting channel topology of MFI zeolite. (3) Length entropy effects comes into force for sorption of linear and branched alkanes within the cylindrical channels of say AFI and MOR zeolites; here the double branched alkane has the shortest length and can be packed more efficiently within the channels. We demonstrate that CBMC simulations allow the efficient screening of zeolite structures for a given separation duty and aid the development of novel separation processes exploiting entropy effects.