Comparison of model potentials for molecular-dynamics simulations of silica

Structural, thermomechanical, and dynamic properties of pure silica SiO2 are calculated with three different model potentials, namely, the potential suggested by van Beest, Kramer, and van Santen (BKS) [Phys. Rev. Lett. 64, 1955 (1990)], the fluctuating-charge potential with a Morse stretch term for the short-range interactions proposed by Demiralp, Cagin, and Goddard (DCG)[Phys. Rev. Lett. 82, 1708 (1999)], and a polarizable force field proposed by Tangney and Scandolo (TS) [J. Chem. Phys. 117, 8898 (2002)]. The DCG potential had to be modified due to flaws in the original treatment. While BKS reproduces many thermomechanical properties of different polymorphs rather accurately, it also shows qualitatively wrong trends concerning the phononic density of states, an absence of the experimentally observed anomaly in the c/a ratio at the quartz alpha-beta transition, pathological instabilities in the beta-cristobalite phase, and a vastly overestimated transition pressure for the stishovite I --> II transition. These shortcomings are only partially remedied by the modified DCG potential but greatly improved by the TS potential. DCG and TS both reproduce a pressure-induced transition from alpha-quartz to quartz II, predicted theoretically based on the BKS potential.     

Special wetting behavior of a graphitic nanofiber-modified epoxy generalized for rough textured fabric surfaces

Wetting behavior of a polymer resin used as matrix on fabric surfaces is one of the key attributes for making high quality structural composites. Though incorporation of various functionalized nanoparticles can stimulate improvements to many properties of epoxy resins, there has not been any report on wettability of any nano-modified epoxy on rough inclined fabric surfaces. In this research work, wetting behavior of a previously developed nano-epoxy resin modified by a type of reactive graphitic nanofibers (r-GNFs) was investigated. The observation results revealed that a unique wetting behavior was discovered from the nano-epoxy on rough fabric surfaces due to the contribution of the r-GNFs. Based on this dramatically improved wettability of the epoxy, a concept of dry–wet contact model was proposed to interpret the different wetting phenomenon observed from the nano-epoxy and that of the pure epoxy. The improved wetting characteristics of the nano-epoxy system will be essential for enabling future energy efficient infusion processing for manufacturing high quality and high-performance structural composite applications.     

Mixed micellization and surface properties of non-ionic/cationic surfactants

We investigate the surface properties of aqueous binary mixtures of our cationic surfactant O-dodecyl-N,N′-diisopropylisourea hydrochloride (ISO-DIC C₁₂) with commercially available nonionic surfactant polyoxyethylene p-(1,1,3,3-tetramethylbutyl)phenyl ether (TritonX-100) at different temperatures (288 to 303?K). The micellization behavior of the binary systems is studied by determining the surface tension and other important physicochemical parameters, such as the critical micelle concentration (CMC), surface tension at the CMC(γ_cmc), Krafft Temperature (T_K), maximum excess concentration (Γ_max), minimum surface area per molecule (A_min), surface pressure at the CMC (П_cmc), and the adsorption efficiency (pC20) at the air/water interface. The study has additionally covered the calculation of thermodynamic parameters of micellization, including the standard Gibbs free energy, the standard enthalpy, the standard entropy, the free energy, and the Gibbs free energy of adsorption at air/water interface. The CMC values of the binary systems determined by experimental data are used to evaluate the micellar composition in the mixed micelle, the interaction parameter β and the activity coefficients f₁(ISO-DIC C₁₂) and f₂ (polyoxyethylene p-(1,1,3,3-tetramethylbutyl)phenyl ether) using the theoretical treatment proposed by Clint and Rubingh. Our results reveal that the proposed binary systems possess enhanced surface activity compared to those of the individual surfactants.     

A radiometric method for the determination of non-ionic surfactants

A radiometric method for the determination of non-ionic surfactants of ethylene oxide type has been elaborated. The principle of the method is the formation of an associate¹³³Ba-non-ionic surfactant — dicarbolide and its extraction into a mixture nitrobenzene-chloroform /41/. The extracted activity of¹³³Ba is proportional to the non-ionic surfactant concentration in the sample. The factors influencing the determination of non-ionic surfactants in the real waste solution after washing have been studied. The possibility of the determination of non-ionic surfactants by the calibration graph method and by the method of standard addition was verified.     

Syntheses and properties of novel non-ionic fluorinated multichains “star-like” surfactants

A series of star-like nonionic surfactants (with two hydrophobic and two hydrophilic chains) with different lengths of hydrophilic and hydrophobic arms were synthesised on the basis of pyromellitic acid dianhydride. The hydrophilic arms were formed by polyoxyethylene and hydrophobic ones either by perfluoro- or by alkyl chains. The adsorption monolayers (Gibbs monolayers) were studied by surface pressure (π) measurements as a function of time for different surfactant concentrations. For the spread monolayers (Langmuir monolayers), the measurements of the surface pressure (π) versus the molecular area (A) as well as the relaxation measurements of the area (A) as a function of time at constant surface pressure were performed. The comparison between the characteristic parameters of two types of monolayers was made in order to understand the effect of the preparation conditions on the structure of these monolayers.It was found that decreasing the fluoroalkyl chain length induced a systematical decrease in the stability of Langmuir monolayers, which is manifested as the Marangoni-Gibbs viscoelasticity of the monolayers. For the surfactants, which have a large number of oxyethylene groups, adsorption at the air/water interface from the bulk solution required extremely long times to reach equilibrium due to the diffusion from the solution and to the conformational rearrangements at the interface. The observation of a hysteresis in the compression/decompression curves for these compounds is explained by the presence of the residual organic solvent molecules absorbed by oxyethylenic chains. A novel model describing the kinetics of desorption or rearrangement of molecules during the lateral compression was suggested, allowing the estimation of both characteristic time of this process and areas per molecule at the equilibrium from the relaxation curves A(t).     

Silicon‐modified surfactants and wetting: IV. Spreading behaviour of trisiloxane surfactants on energetically different solid surfaces

The spreading behaviour of defined trisiloxane surfactants of general formula [(CH₃)₃SiO]₂ CH₃Si(CH₂)₃(OCH₂CH₂) _nOCH₃ (n = 3–9) on five different solid surfaces has been investigated. Maximum spreading areas and rates are found on non‐polar or slightly polar surfaces of 30 to 40 mN m⁻¹ surface energy. Extremely low or high surface energies substantially reduce the spreading rates. On non‐polar surfaces rapid spreading is observed for 1 wt % solutions of the relatively short‐chained penta‐ and hexa‐ethylene glycol derivatives. On slightly polar surfaces dilute 0.1 wt % solutions of longer‐chained derivatives spread faster. This spreading pattern shift coincides with a change of the phase behaviour. Solutions of Silwet L77 do not prefer one specific surface, since 1 wt % solutions abruptly stop spreading after a few seconds and the maximum spreading rates are found for 0.1 wt % solutions. Therefore, Silwet L77 essentially belongs among the long‐chained derivatives. Copyright © 2000 John Wiley & Sons, Ltd.     

Synthesis and properties of new glucocationic surfactants: model structures for marking cationic surfactants with carbohydrates

[reaction: see text] In this work, we report the synthesis of a new series of glucocationic surfactants, a class of surfactants we introduced very recently. The preparation of the surfactants is based on the synthesis of the 2-bromoethyl-2,3,4,6-tetra-O-acetyl-beta-D-glucopyranoside, whose preparation was studied in order to improve yields and stereoselectivity of this key intermediate. These glucocationic amphiphiles were prepared and studied as a model of cationic surfactants marked with a carbohydrate moiety. The use of carbohydrates as markers on cationic lipids was recently introduced to induce recognition by specific receptors, present on the surface of cell membranes. The chemicophysical characterization of these model structures can give more insight on the aggregation behavior. Conductivity and surface tension measurements were performed in order to characterize the compounds from the amphiphilic point of view. The results showed a different effect of the glucosidic moiety on the cmc value with respect to the glucopyridinium cationic surfactants. The surfactants also showed the tendency to form premicellar aggregates in solution when the hydrophobicity is raised.     

Silicon-modified carbohydrate surfactants. IV. The impact of substructures on the wetting behaviour of siloxanyl-modified carbohydrate surfactants on low-energy surfaces

The siloxanyl-modified carbohydrate surfactants investigated consist of the four structural elements: (1) siloxanyl moiety; (2) spacer; (3) carbohydrate unit; and (4) modifying element. By static surface tension (γ_sν — γ_sβ α) measurements the contact angles of the aqueous surfactant solutions above the critical micelle formation concentration (cmc) on nonpolar perfluorinated surfaces (FEP^® plate) were determined. Although the siloxanyl units were found to have a high capacity to level out the interfacial properties, both surface tension and wetting tension react independently to defined changes in the chemical structure of the surfactant molecules. The results of spreading experiments on polyproylene show good correlation with the dependences found by wetting meaurements. © 1997 John Wiley & Sons, Ltd.     

Evolution and potential of the BiAS procedure for the determination of non-ionic surfactants

The evolution of the BiAS procedure is reviewed and its standard recommended version was compared with the modified version combined with the indirect tensammetric method (BiAS-ITM). New applications of the use of BiAS-ITM for the determination of nonionic surfactants (in the presence of hydrocarbons or adsorbed on particles) and polyethylene glycols were discussed.     

An extraction—spectrophotometmc method for the determination of non-ionic surfactants

The surfactant is extracted into 1,2-dichlorobenzene as a neutral adduct with potassium tetrathiocyanatozincate(II), and zinc(II) in the extract is determined spectrophotometrically after addition of l-(2-pyridylazo)-2-naphthol and triethanolamine. With a 150-ml water sample, the limit of detection is 15 μg l^-1 (as Triton X-100). The method requires only one extraction and is applicable, without modification, to fresh, estuarine and sea-water samples.     

An atomic absorption spectrometric method for the determination of non-ionic surfactants

A method is described for the determination of non-ionic surfactants in the concentration range 0.05–2 mg l^-1.Surfactant molecules are extracted into 1,2-dichlorobenzene as a neutral adduct with potassium tetrathiocyanatozincate(II) and the determination is completed by atomic absorption spectrometry. With a 150-ml water sample, the limit of detection is 0.03 mg l^-1(as Triton X-100).The method requires a single phase separation step, is applicable, without modification, to fresh, estuarine and sea-water samples and is relatively free from interference by anionic surfactants; the presence of up to 5 mg l^-1 of anionic surfactant (as LAS) introduces an error of no more than 0.07 mg l^-1 (as Triton X-100) in the apparent non-ionic surfactant concentration.     

Synchronization of trajectories in canonical molecular-dynamics simulations: observation, explanation, and exploitation

For two methods commonly used to achieve canonical-ensemble sampling in a molecular-dynamics simulation, the Langevin thermostat and the Andersen [H. C. Andersen, J. Chem. Phys. 72, 2384 (1980)] thermostat, we observe, as have others, synchronization of initially independent trajectories in the same potential basin when the same random number sequence is employed. For the first time, we derive the time dependence of this synchronization for a harmonic well and show that the rate of synchronization is proportional to the thermostat coupling strength at weak coupling and inversely proportional at strong coupling with a peak in between. Explanations for the synchronization and the coupling dependence are given for both thermostats. Observation of the effect for a realistic 97-atom system indicates that this phenomenon is quite general. We discuss some of the implications of this effect and propose that it can be exploited to develop new simulation techniques. We give three examples: efficient thermalization (a concept which was also noted by Fahy and Hamann [S. Fahy and D. R. Hamann, Phys. Rev. Lett. 69, 761 (1992)]), time-parallelization of a trajectory in an infrequent-event system, and detecting transitions in an infrequent-event system.     

Silicon‐modified surfactants and wetting: V. The spreading behaviour of trimethylsilane surfactants on energetically different solid surfaces

The spreading behaviour of defined trimethylsilane‐based surfactants of general formula (CH₃)₃Si(CH₂)₆(OCH₂CH₂) _nOCH₃, n = 2–6, on five different solid surfaces at 21 °C has been investigated. Compounds bearing short diethylene and triethylene glycol hydrophiles do not spread. For the longer‐chained tetraethylene to hexaethylene glycol derivatives, the ability to spread depends on the surface energy. Rapid spreading is restricted to the slightly polar surface of 40 mN m⁻¹ surface energy. Lower or higher surface energies considerably reduce the spreading rates. The phase behaviour of the solutions substantially influences the spreading process. The dispersed systems of the tetraethylene glycol derivative spread constantly over long time intervals. The dispersions of the pentaethylene glycol analogue are very close to the temperature for a transition into the one‐phase state. A retardation of the spreading process occurs after a few seconds. Micellar solutions of the hexaethylene glycol derivative either spread very slowly or stop spreading after a few seconds. The largest spreading areas and highest initial spreading rates were found for the 0.1 wt% solutions. Copyright © 1999 John Wiley & Sons, Ltd.     

Tests of the homogeneous nucleation theory with molecular-dynamics simulations. I. Lennard-Jones molecules

Two kinds of the homogeneous nucleation theory exist at the present: the classical nucleation theory and the semiphenomenological model. To test them, we performed molecular-dynamics (MD) simulations of nucleation from vapor to liquid with 5000-20,000 Lennard-Jones-type molecules. Simulations were done for various values of supersaturation ratios (from 2 to 10) and temperatures (from 80 to 120 K). We compared the size distribution of clusters in MD simulations with those in the theoretical models because the number density of critical clusters governs the nucleation rate. We found that the semiphenomenological model achieves excellent agreements in size distributions of the clusters with all MD simulations we done. The classical theory underestimates the number density of the clusters in the temperature range of 80-100 K, but overestimates in 100-120 K. The semiphenomenological model also predicts well the nucleation rate in MD simulations, while the classical nucleation theory does not. Our results confirmed the validity of the semiphenomenological model for Lennard-Jones-type molecules.     

Monte Carlo simulation methods for computing the wetting and drying properties of model systems

We introduce general Monte Carlo simulation methods for determining the wetting and drying properties of model systems. We employ an interface-potential-based approach in which the interfacial properties of a system are related to the surface excess free energy of a thin fluid film in contact with a surface. Two versions of this approach are explored: a "spreading" method focused on the growth of a thin liquid film from a surface in a mother vapor and a "drying" method focused on the growth of a thin vapor film from a surface in a mother liquid. The former provides a direct measure of the spreading coefficient while the latter provides an analogous drying coefficient. When coupled with an independent measure of the liquid-vapor surface tension, these coefficients enable one to compute the contact angle. We also show how one can combine information gathered from application of the spreading and drying methods at a common state point to obtain direct measures of the contact angle and liquid-vapor surface tension. The computational strategies introduced here are applied to two model systems. One includes a monatomic Lennard-Jones fluid that interacts with a structureless substrate via a long-ranged substrate potential. The second model contains a monatomic Lennard-Jones fluid that interacts with an atomistically detailed substrate via a short-ranged potential. Expanded ensemble techniques are coupled with the interface potential approach to compile the temperature- and substrate strength-dependence of various interfacial properties for these systems. Overall, we find that the approach pursued here provides an efficient and precise means to calculate the wetting and drying properties of model systems.     

Structure, thermodynamics, and liquid-vapor equilibrium of ethanol from molecular-dynamics simulations using nonadditive interactions

We present a molecular-dynamics simulation study of the bulk and liquid-vapor interfacial properties of ethanol using a polarizable force field based on the fluctuating charge (FQ) formalism, as well as the nonpolarizable CHARMM22 force field. Both models are competitive with respect to the prediction of ambient liquid properties such as liquid density, enthalpy of vaporization, dielectric constant, and self-diffusion constants. The polarizable model predicts an average condensed-phase dipole moment of 2.2 D associated with an induced liquid-phase dipole moment of 0.6 D; though qualitatively in agreement with earlier nonadditive models as well as recent Car-Parinello calculations, the current FQ model underestimates the condensed-phase dipole moment. In terms of liquid structure, both models are in agreement with recent neutron-diffraction results of liquid ethanol structure, although the polarizable model predicts the hydroxyl-hydrogen-hydroxyl-hydrogen structure factor in closer agreement with the experimental data. In terms of interfacial properties, both models predict ambient surface tension to within 4% of the experimental value of 22.8 dyncm, while overestimating the surface excess entropy by almost a factor of 2. Both models display the characteristic preferential orientation of interfacial molecules. The polarizable model allows for a monotonic variation of the average molecular dipole moment from the bulk value to that of the vapor phase. Consequently, there is a dramatic difference in the surface potential predicted by the polarizable and nonpolarizable models. The polarizable model estimates a surface potential of -209+/-3 mV, while the nonpolarizable model yields a value of -944+/-10 mV. Finally, based on the vapor-liquid equilibrium simulation data from several temperatures, we estimate the critical properties of both models. As observed with other FQ models for associating fluids (such as water and methanol), and counter to what one would anticipate by modeling more physically the electrostatic response to local environment, the current FQ model underestimates the critical temperature and overestimates the critical density of ethanol; moreover, the FQ model is, in this respect, equivalent to the underlying fixed-charge model. These results further suggest the need to revisit polarizable models in terms of quantitative vapor-liquid equilibrium prediction.