Genetic Parameterization of Interfacial Force Fields Based on Classical Bulk Force Fields and Ab Initio Data: Application to the Methanol-ZnO Interfaces

Understanding the Role of Solvent on the Growth of Zinc Oxide: Insight from Experiment and Molecular Dynamics Simulations

The controlled synthesis of nanoparticles with tailored shapes and morphologies has garnered significant attention, driven by the ever-growing demand for advanced materials with defined properties. In nanoparticle formation, various parameters influence the final product, and among these, the solvent plays a pivotal role, as it constitutes the major component of the reaction medium. In this work, the critical role of solvents in controlling the growth of zinc oxide (ZnO) nanoparticles was investigated, with a focus on simple primary alcoholic solvents as the reaction medium. A model reaction based on the direct solvolysis of anhydrous zinc acetylacetonate was employed to probe the influence of different primary alcohols, specifically methanol, ethanol, and their mixture. A substantial difference in the preferential growth direction of the ZnO nanocrystals in methanol and ethanol was observed through XRD and was further proven through TEM. Thereby, in ethanol, a preferential growth in the [001] direction was observed, resulting in short nanorods as primary particles, while this growth was inhibited in methanol, leading to platelet- or sheet-like primary particles. To unravel the underlying mechanisms responsible for the observed solvent-dependent variations, molecular dynamics (MD) simulations were employed using an optimized interface force field to model the ZnO-alcohol interaction. These simulations provide valuable insights into the preferential adsorption of the solvent molecules onto the polar (0001) and (0001̅) and nonpolar (101̅0) ZnO surfaces, shedding light on the fundamental interactions driving the shape control phenomenon. Essentially, the experimental observations on primary particle morphology could be explained well by the adsorption behavior determined by the MD simulations. Furthermore, this report provides an extensive comparison with various similar reaction systems for ZnO synthesis, deriving correlations with the findings from the model system. These insights contribute to a deeper understanding of the intricate interplay between solvent properties and nanoparticle growth, offering a valuable toolkit for designing and optimizing the synthesis of ZnO nanoparticles with specific shapes and functionalities.

Transferable Gaussian Attractive Potentials for Organic/Oxide Interfaces

Organic/oxide interfaces play an important role in many areas of chemistry and in particular for lubrication and corrosion. Molecular dynamics simulations are the method of choice for providing complementary insight to experiments. However, the force fields used to simulate the interaction between molecules and oxide surfaces tend to capture only weak physisorption interactions, discarding the stabilizing Lewis acid/base interactions. We here propose a simple complement to the straightforward molecular mechanics description based on "out-of-the-box" Lennard-Jones potentials and electrostatic interactions: the addition of an attractive Gaussian potential between reactive sites of the surface and heteroatoms of adsorbed organic molecules, leading to the Gaussian Lennard-Jones (GLJ) potential. The interactions of four oxygenated and four amine molecules with the typical and widespread hematite and γ-alumina surfaces are investigated. The root mean square deviation (RMSD) for all probed molecules is only 5.7 kcal/mol, which corresponds to an error of 23% over hematite. On γ-alumina, the RMSD is 11.2 kcal/mol using a single parameter for all five chemically inequivalent surface aluminum atoms. Applying GLJ to the simulation of organic films on oxide surfaces demonstrates that the mobility of the surfactants is overestimated by the simplistic LJ potential, while GLJ and other qualitatively correct potentials show a strong structuration and slow dynamics of the surface films, as could be expected from the first-principles adsorption energies for model head groups.

Efficient Parametrization of Force Field for the Quantitative Prediction of the Physical Properties of Ionic Liquid Electrolytes

The prediction of transport properties of room-temperature ionic liquids from nonpolarizable force field-based simulations has long been a challenge. The uniform charge scaling method has been widely used to improve the agreement with the experiment by incorporating the polarizability and charge transfer effects in an effective manner. While this method improves the performance of the force fields, this prescription is ad hoc in character; further, a quantitative prediction is still not guaranteed. In such cases, the nonbonded interaction parameters too need to be refined, which requires significant effort. In this work, we propose a three-step semiautomated refinement procedure based on (1) atomic site charges obtained from quantum calculations of the bulk condensed phase; (2) quenched Monte Carlo optimizer to shortlist suitable force field candidates, which are then tested using pilot simulations; and (3) manual refinement to further improve the accuracy of the force field. The strategy is designed in a sequential manner with each step improving the accuracy over the previous step, allowing the users to invest the effort commensurate with the desired accuracy of the refined force field. The refinement procedure is applied on N,N-diethyl-N-methyl-N-(2-methoxyethyl)ammonium bis(trifluoromethanesulfonyl)imide (DEME-TFSI), a front-runner as an electrolyte for electric double-layer capacitors and single-molecule-based devices. The transferability of the refined force field is tested on N,N-dimethyl-N-ethyl-N-methoxyethoxyethylammonium bis(trifluoromethanesulfonyl)imide (N_112,2O2O1-TFSI). The refined force field is found to be better at predicting both structural and transport properties compared to the uniform charge scaling procedure, which showed a discrepancy in the X-ray structure factor. The refined force field showed quantitative agreement with structural (density and X-ray structure factor) and transport properties-diffusion coefficients, ionic conductivity, and shear viscosity over a wide temperature range, building a case for the wide adoption of the procedure.