Abstract:The MARTINI force field has become increasingly popular in material science by virtue of its versatility and its building block approach. However, force field parametrization is primarily based on reproducing structural and thermodynamic properties of the building blocks in solution, which may not be the most suitable strategy for regular, extended surfaces at solid-liquid interfaces, e.g., the planar facets of nanoparticles. In this study, we propose a novel strategy for parametrizing MARTINI planar surfaces in contact with liquids. While the choice of bead types in MARTINI models for molecules and ions in solution is often dictated by their chemistry, we suggest that surface-liquid interactions can be refined through a top-down approach. Specifically, we propose tuning the density of particles in the model surface to reproduce available experimental data on contact angles. In addition, we suggest several possible solutions to mitigate an artifact of the force field, that is, the excessive ordering of the liquid in contact with regular surfaces. Our findings demonstrate that water ordering can be reduced by tuning surface roughness and the level of coarse-graining of water. Eventually, we validate our approach on silica surfaces functionalized by alkyl chains.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the challenges encountered in modeling the solid - liquid interface in the MARTINI 3 force field. Specifically: 1. **Excessive ordering of liquid on the solid surface**: When simulating the solid - liquid interface, water molecules using the coarse - grained (CG) model will form too many ordered layers when in contact with the solid surface, which is inconsistent with the experimentally observed phenomena. This excessive ordering phenomenon is particularly evident on very regular surfaces (such as the planar crystal faces of nanoparticles). 2. **Accurate modeling of surface - liquid interactions**: Traditional MARTINI force - field parameterizations are mainly based on the structure and thermodynamic properties of building blocks in solution, which may not be applicable to regularly extended surfaces on the solid - liquid interface. Therefore, new strategies need to be developed to optimize the interactions between these surfaces and the liquid. 3. **Selection of surface geometry and bead types**: Selecting appropriate bead types and surface geometries is crucial for accurately describing the solid - liquid interface. However, the definition of the actual surface at the atomic scale is often unclear, solvent - dependent reconstruction may occur, and there may be defects or functionalized molecules. ### Solutions To address the above challenges, the authors propose the following strategies: 1. **Introducing surface disorder**: - **Static disorder**: Introduce surface roughness by randomly displacing the positions of surface beads. - **Dynamic disorder**: Allow surface beads to vibrate near their equilibrium positions. 2. **Using a mixture of water beads of different sizes**: - Mix water beads of different sizes (W, SW, TW) to reduce the excessive ordering of water molecules when in contact with the surface. 3. **Adjusting the surface bead density**: - Optimize the interaction between the surface and the solvent by adjusting the density of the surface beads so that the simulated contact angle is consistent with the experimental data. ### Experimental verification The authors verified the effectiveness of these strategies through a series of molecular dynamics (MD) simulations. For example, they used different types of surface models (p - surface, h - surface, r - surface) and mixtures of water beads of different sizes for simulations, and evaluated the ordering degree of water molecules by analyzing the number of water molecule layers and density profiles. ### Conclusion Through these strategies, the authors successfully reduced the excessive ordering phenomenon of water molecules when in contact with the solid surface, and were able to match the contact angles measured experimentally by adjusting the surface bead density. These methods provide important guidance for more accurate simulation of the solid - liquid interface in the MARTINI 3 force field. Hope this information is helpful to you! If you have any other questions, feel free to continue asking.

Modeling Solid-Liquid Interfaces in MARTINI 3

The surface tension of Martini 3 water mixtures

Coarse-Graining the Liquid-Liquid Interfaces with the MARTINI Force Field: How Is the Interfacial Tension Reproduced?

Martini 3: a general purpose force field for coarse-grained molecular dynamics

The Martini Model in Materials Science

Martini 3 Coarse-Grained Force Field for Carbohydrates

Many-Body Dissipative Particle Dynamics with the MARTINI "Lego" approach

Simulations of Interfacial Tension of Liquid-Liquid Ternary Mixtures Using Optimized Parametrization for Coarse-Grained Models

Coarse Grained modeling of Zeolitic Imidazolate Framework-8 using MARTINI Force Fields

Building complex membranes with Martini 3

Martinize2 and Vermouth: Unified Framework for Topology Generation

Bridging the gap between mesoscopic and molecular models of solid/liquid interfaces out-of-equilibrium

Martini 3 coarse-grained models for carbon nanomaterials

First-principles molecular dynamics simulations at solid-liquid interfaces with a continuum solvent

Automated Coarse-Grained Mapping Algorithm for the Martini Force Field and Benchmarks for Membrane–Water Partitioning

Recalibration of MARTINI-3 Parameters for Improved Interactions between Peripheral Proteins and Lipid Bilayers

Solid-Fluid Interaction with Surface-Tension-Dominant Contact

Martignac: Computational workflows for reproducible, traceable, and composable coarse-grained Martini simulations

Mesoscale modeling of near-contact interactions for complex flowing interfaces

Reparameterization of Polarizable Force Fields for Studying Ion Transfer across Liquid-Liquid Interfaces