Abstract:The spontaneous formation of micelles in aqueous solutions is governed by the amphipathic nature of surfactants and is practically interesting due to the regular use of micelles as membrane mimics, for the characterization of protein structure, and for drug design and delivery. We performed a systematic characterization of the finite-size effect observed in single-component dodecylphosphocholine (DPC) micelles with the coarse-grained MARTINI model. Of multiple coarse-grained solvent models investigated using large system sizes, the non-polarizable solvent model was found to most-accurately reproduce SANS spectra of 100 mM DPC in aqueous solution. We systematically investigated the finite-size effect at constant 100 mM concentration in 23 systems of sizes 40 to 150 DPC, confirming the finite-size effect to manifest as an oscillation in the mean micelle aggregation number about the thermodynamic aggregation number as the system size increases, mostly diminishing once the system supports formation of three micelles. The accuracy of employing a multiscale simulation approach to avoid finite-size effects in the micelle size distribution and SANS spectra using MARTINI and CHARMM36 was explored using multiple long timescale 500-DPC coarse-grained simulations which were back-mapped to CHARMM36 all-atom systems. It was found that the MARTINI model generally occupies more volume than the all-atom model, leading to the formation of micelles that are of a reasonable radius of gyration, but are smaller in aggregation number. The systematic characterization of the finite-size effect and exploration of multiscale modeling presented in this work provides guidance for the accurate modeling of micelles in simulations.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the influence of system size on the micelle aggregation number (AN) distribution, namely the finite - size effect, during the simulation of surfactant micelle self - assembly. Specifically, by systematically characterizing dodecylphosphocholine (DPC) micelles in systems of different sizes, the author explored how to select an appropriate system size to avoid or reduce the influence of the finite - size effect on the simulation results when using the coarse - grained (CG) MARTINI model for simulation. In addition, the study also investigated the effectiveness of the multi - scale simulation method, that is, reverse - mapping from the CG model to the all - atom (AA) model, to more accurately simulate the properties of micelles, such as radius distribution, SANS spectra, etc. ### Main research questions: 1. **Influence of finite - size effect on micelle aggregation number distribution**: Study the change of the average aggregation number of micelles with system size under different system sizes (40 to 150 DPC molecules), and whether this change conforms to theoretical predictions. 2. **Selection of the optimal system size**: Determine the minimum system size that can produce an aggregation number distribution close to the thermodynamic limit when simulating the self - assembly of DPC micelles. 3. **Effectiveness of the multi - scale simulation method**: Evaluate the accuracy of the method of reverse - mapping from the CG model to the AA model in simulating the properties of micelles, especially the radius distribution and SANS spectra of micelles. ### Research methods: - **Selection of the CG model**: Four different MARTINI solvent models (non - polarized water, polarized water, improved polarized water, improved polarized water with polarized ions) were compared, and the non - polarized water model that best matched the experimental SANS spectra was selected. - **Characterization of the finite - size effect**: Systems containing 40 to 150 DPC molecules were constructed, with one system for every 5 DPC molecules. Large - scale simulations were carried out and the aggregation number distribution of micelles was analyzed. - **Multi - scale simulation**: A CG model of 500 DPC molecules was used for simulation, and then reverse - mapped to the CHARMM36 all - atom model. The differences in micelle properties between the two models were compared. ### Key findings: - **Characterization of the finite - size effect**: As the system size increases, the average aggregation number of micelles shows an oscillation phenomenon, and when the system can form three micelles, the finite - size effect basically disappears. - **Optimal system size**: It was determined that system sizes of 45, 90, and 135 DPC molecules can better avoid the finite - size effect, and the micelle aggregation numbers under these system sizes are close to the thermodynamic limit values. - **Effectiveness of multi - scale simulation**: After reverse - mapping to the AA model, the radius distribution and SANS spectra of micelles are closer to the experimental results, but there are still some deviations, indicating that the CG model may overestimate the volume of micelles in some aspects. Through these studies, the author provides guidance for future simulations of micelle self - assembly, especially in terms of system size and model selection.

Finite-size effects and optimal system sizes in simulations of surfactant micelle self-assembly

Non-equilibrium effects of micelle formation as studied by a minimum particle-based model

Coarse-grained simulations of concentration-dependent molecular self-assembly of Polysorbate 80 in water

Influence of hydrophobic micelle structure on crystallization of the photosynthetic RC-LH1-PufX complex from Rhodobacter blasticus

Simulating micelle self-assembly to assess potential for viscosity build in surfactant formulations

DPD simulations of anionic surfactant micelles: a critical role for polarisable water models

Molecular dynamics simulation of self-assembly of n-decyltrimethylammonium bromide micelles

Analytical modeling of micelle growth. 4. Molecular thermodynamics of wormlike micelles from ionic surfactants: theory vs. experiment

Physical Insight into the Conditions Required in the Solid-Phase Molecular Self-Assembly of SDS Revealed by Coarse-Grained Molecular Dynamics Simulation

Monte Carlo Simulation for the Self-Aggregation of Nonionic Surfactants

Finite size effects and optimization of the calculation of the surface tension in surfactant mixtures at liquid/vapour interfaces

Micellization of Telechelic Associative Polymers: Self-Consistent Field Modeling and Comparison with Scaling Concepts.

Experimental and simulation study of reverse micelles formed by aerosol-OT and water in non-polar solvents

Quantifying the Critical Micelle Concentration of Nonionic and Ionic Surfactants by Self-Consistent Field Theory

A molecular dynamics study of free energy of micelle formation for sodium dodecyl sulfate in water and its size distribution.

Molecular Simulation of a Dichain Surfactant/Water/Carbon Dioxide System. 1. Structural Properties of Aggregates

Molecular thermodynamic modeling of specific ion effects on micellization of ionic surfactants

Simple and complex micelles in amphiphilic mixtures: a coarse-grained mean-field study

Efficient Simulation of Tunable Lipid Assemblies Across Scales and Resolutions

Experimental and Molecular Dynamics Studies