Development of the RF-MEAM Interatomic Potential for the Fe-C System to Study the Temperature-Dependent Elastic Properties.

One of the major impediments to the computational investigation and design of complex alloys such as steel is the lack of effective and versatile interatomic potentials to perform large-scale calculations. In this study, we developed an RF-MEAM potential for the iron-carbon (Fe-C) system to predict the elastic properties at elevated temperatures. Several potentials were produced by fitting potential parameters to the various datasets containing forces, energies, and stress tensor data generated using density functional theory (DFT) calculations. The potentials were then evaluated using a two-step filter process. In the first step, the optimized RSME error function of the potential fitting code, MEAMfit, was used as the selection criterion. In the second step, molecular dynamics (MD) calculations were employed to calculate ground-state elastic properties of structures present in the training set of the data fitting process. The calculated single crystal and poly-crystalline elastic constants for various Fe-C structures were compared with the DFT and experimental results. The resulting best potential accurately predicted the ground state elastic properties of B1, cementite, and orthorhombic-Fe7C3 (O-Fe7C3), and also calculated the phonon spectra in good agreement with the DFT-calculated ones for cementite and O-Fe7C3. Furthermore, the potential was used to successfully predict the elastic properties of interstitial Fe-C alloys (FeC-0.2% and FeC-0.4%) and O-Fe7C3 at elevated temperatures. The results were in good agreement with the published literature. The successful prediction of elevated temperature properties of structures not included in data fitting validated the potential's ability to model elevated-temperature elastic properties.

Model for the Simulation of the CnEm Nonionic Surfactant Family Derived from Recent Experimental Results.

Using a comprehensive set of recently published experimental results for training and validation, we have developed computational models appropriate for simulations of aqueous solutions of poly(ethylene oxide) alkyl ethers, an important class of micelle-forming nonionic surfactants, usually denoted CnEm. These models are suitable for use in simulations that employ a moderate amount of coarse graining and especially for dissipative particle dynamics (DPD), which we adopt in this work. The experimental data used for training and validation were reported earlier and produced in our laboratory using dynamic light scattering (DLS) measurements performed on 12 members of the CnEm compound family yielding micelle size distribution functions and mass-weighted mean aggregation numbers at each of several surfactant concentrations. The range of compounds and quality of the experimental results were designed to support the development of computational models. An essential feature of this work is that all simulation results were analyzed in a way that is consistent with the experimental data. Proper account is taken of the fact that a broad distribution of micelle sizes exists, so mass-weighted averages (rather than number-weighted averages) over this distribution are required for the proper comparison of simulation and experimental results. The resulting DPD force field reproduces several important trends seen in the experimental critical micelle concentrations and mass-averaged mean aggregation numbers with respect to surfactant characteristics and concentration. We feel it can be used to investigate a number of open questions regarding micelle sizes and shapes and their dependence on surfactant concentration for this important class of nonionic surfactants.

Challenge to Reconcile Experimental Micellar Properties of the CnEm Nonionic Surfactant Family.

We wished to compile a data set of results from the experimental literature to support the development and validation of accurate computational models (force fields) for an important class of micelle-forming nonionic surfactant compounds, the poly(ethylene oxide) alkyl ethers, usually denoted C nE m. However, careful examination of the experimental literature exposed a striking degree of variation in values reported for critical micelle concentrations (cmc) and mean aggregation numbers ( Nagg). This variation was so large that it masked important trends known to exist within this family of molecules, thereby rendering most of the literature data to be of limited utility for force field development. In this work, we describe some reasons for the wide variability in the experimental literature, and we present a set of cmc and aggregation number data for 12 C nE m compounds that we feel is appropriate to use for the construction of and validation of computational models. The cmc values we selected are from the existing experimental literature and represent a carefully chosen and consistent subset that conveys important trends seen by many of the experimental studies. However, for a corresponding and consistent set of weight-averaged aggregation numbers, we needed to perform new dynamic light scattering (DLS) experiments. The results of these experiments were carefully analyzed to obtain not just mean aggregation numbers but also the underlying micelle size distribution functions. Several trends observed in the cmc and Nagg observables are highlighted and serve as challenges for developers of force field and simulation methodology. The analysis of the DLS experiments accounts for the fact that a broad distribution of micelle sizes exists for many of these compounds and that one must be careful to use the appropriate weighted averages (e.g., mass-weighted vs number-weighted averages) in comparing results from different types of experiments and in comparing results from experiments with those from simulations.