We have employed two-to-one mapping scheme to develop three coarse-grained (CG) water models, namely, 1-, 2-, and 3-site CG models. Here, for the first time, particle swarm optimization (PSO) and gradient descent methods were coupled to optimize the force-field parameters of the CG models to reproduce the density, self-diffusion coefficient, and dielectric constant of real water at 300 K. The CG MD simulations of these new models conducted with various timesteps, for different system sizes, and at a range of different temperatures are able to predict the density, self-diffusion coefficient, dielectric constant, surface tension, heat of vaporization, hydration free energy, and isothermal compressibility of real water with excellent accuracy. The 1-site model is ∼3 and ∼4.5 times computationally more efficient than 2- and 3-site models, respectively. To utilize the speed of 1-site model and electrostatic interactions offered by 2- and 3-site models, CG MD simulations of 1:1 combination of 1- and 2-/3-site models were performed at 300 K. These mixture simulations could also predict the properties of real water with good accuracy. Two new CG models of benzene, consisting of beads with and without partial charges, were developed. All three water models showed good capacity to solvate these benzene models.

Optimizing force-field (FF) parameters to perform molecular dynamics (MD) simulations is a challenging and time-consuming process. We present a novel FF optimization framework that integrates MD simulations with particle swarm optimization (PSO) algorithm and artificial neural network (ANN). This new ANN-assisted PSO framework was used to develop transferable coarse-grained (CG) models for DO and DMF as a proof of concept. The PSO algorithm was used to generate the set of input FF parameters for the MD simulations of the CG models of these solvents, which were optimized to reproduce their experimental properties. Herein, for the first time, a reverse approach was employed for on-the-fly training of the ANN model, where results (solvent properties) obtained from the MD simulations and their corresponding FF parameters were used as inputs and outputs, respectively. The ANN model was then required to predict a set of new FF parameters, which were tested for their ability to predict the desired experimental properties. This new framework can be extended to integrate any optimization algorithm with ANN and MD simulations to accelerate the FF development.

We have utilized an approach that integrates molecular dynamics (MD) simulations with particle swarm optimization (PSO) to accelerate the development of coarse-grained (CG) models of hydrocarbons. Specifically, we have developed new transferable CG beads, which can be used to model the hydrocarbons (C5 to C17) and reproduce their experimental properties with good accuracy. First, the PSO method was used to develop the CG beads of the decane model represented with a 2:1 (2-2-2-2-2) mapping scheme. This was followed by the development of the nonane model described with hybrid 2-2-3-2 and 3:1 (3-3-3) mapping schemes. The force-field parameters for these three CG models were optimized to reproduce four experimentally observed properties including density, enthalpy of vaporization, surface tension, and self-diffusion coefficient at 300 K. The CG MD simulations conducted with these new CG models of decane and nonane, at different timesteps, for various system sizes, and at a range of different temperatures, were able to predict their density, enthalpy of vaporization, surface tension, self-diffusion coefficient, expansibility, and isothermal compressibility with good accuracy. Moreover, a comparison of structural features obtained from the CG MD simulations and the CG beads of mapped all-atom trajectories of decane and nonane showed very good agreement. To test the chemical transferability of these models, we have constructed the models for hydrocarbons ranging from pentane to heptadecane, by using different combinations of the CG beads of decane and nonane. The properties of pentane to heptadecane predicted by these new CG models showed excellent agreement with the experimental data.

Four different machine learning (ML) regression models: artificial neural network, k-nearest neighbors, gaussian process regression and random forest were built to backmap coarse-grained models to all-atom models. The ML models...

We present a computational framework that integrates coarse-grained (CG) molecular dynamics (MD) simulations and a data-driven machine-learning (ML) method to gain insights into the conformations of polymers in solutions. We employ this framework to study conformational transition of a model thermosensitive polymer, poly­(N-isopropyl­acrylamide) (PNIPAM). Here, we have developed the first of its kind, a temperature-independent CG model of PNIPAM that can accurately predict its experimental lower critical solution temperature (LCST) while retaining the tacticity in the presence of an explicit water model. The CG model was extensively validated by performing CG MD simulations with different initial conformations, varying the radius of gyration of chain, the chain length, and the angle between the adjacent monomers of the initial configuration of PNIPAM (total simulation time = 90 μs). Moreover, for the first time, we utilize the nonmetric multidimensional scaling (NMDS) method, a data-driven ML approach, to gain further insights into the mechanisms and pathways of this coil-to-globule transition by analyzing CG MD simulation trajectories. NMDS analysis provides entirely new insights and shows multiple metastable states of PNIPAM during its coil-to-globule transition above the LCST.