A coarse-grained model for polylactide: glass transition temperature and conformational properties

A Temperature-Transferable Coarse-Grained Model for Poly(lactic Acid) Melts

We present a temperature-transferable coarse-grained (CG) model for poly(lactic acid) (PLA), specifically designed to replicate its volumetric properties and solubility parameter in the molten state. The CG-bonded potentials were derived by using the iterative Boltzmann inversion (IBI) optimization method to match structural properties from detailed atomistic models. A parametrization workflow was employed to determine nonbonded interaction parameters with temperature-dependent corrections that provide agreement with the target properties across the melting temperature range. The CG model successfully replicates key features of the PLA melt. It satisfactory reproduces the density and solubility parameter, maintains the dependence of chain conformation on molecular weight, and captures the dynamic behavior through agreement in scaled mean squared displacement and diffusion coefficients with the atomistic model. Additionally, the CG model offers much faster equilibration compared with the atomistic model. The proposed model is expected to be particularly useful for investigating the miscibility characteristics of PLA in various blends and composites that remain challenging to explore using fully atomistic simulations or experiments.

Physics-Informed Deep Learning Approach for Reintroducing Atomic Detail in Coarse-Grained Configurations of Multiple Poly(lactic acid) Stereoisomers

Multiscale modeling of complex molecular systems, such as macromolecules, encompasses methods that combine information from fine and coarse representations of molecules to capture material properties over a wide range of spatiotemporal scales. Being able to exchange information between different levels of resolution is essential for the effective transfer of this information. The inverse problem of reintroducing atomistic degrees of freedom in coarse-grained (CG) molecular configurations is particularly challenging as, from a mathematical point of view, it is an ill-posed problem; the forward mapping from the atomistic to the CG description is typically defined via a deterministic operator ("one-to-one" problem), whereas the reversed mapping from the CG to the atomistic model refers to creating one representative configuration out of many possible ones ("one-to-many" problem). Most of the backmapping methods proposed so far balance accuracy, efficiency, and general applicability. This is particularly important for macromolecular systems with different types of isomers, i.e., molecules that have the same molecular formula and sequence of bonded atoms (constitution) but differ in the three-dimensional configurations of their atoms in space. Here, we introduce a versatile deep learning approach for backmapping multicomponent CG macromolecules with chiral centers, trained to learn structural correlations between polymer configurations at the atomistic level and their corresponding CG descriptions. This method is intended to be simple and flexible while presenting a generic solution for resolution transformation. In addition, the method is aimed to respect the structural features of the molecule, such as local packing, capturing therefore the physical properties of the material. As an illustrative example, we apply the model on linear poly(lactic acid) (PLA) in melt, which is one of the most popular biodegradable polymers. The framework is tested on a number of model systems starting from homopolymer stereoisomers of PLA to copolymers with randomly placed chiral centers. The results demonstrate the efficiency and efficacy of the new approach.

Temperature-Transferable Coarse-Grained Models for Volumetric Properties of Poly(lactic Acid)

A new coarse-grained (CG) model, for which each monomer is mapped as one bead at its center of mass, was developed for simulating the volumetric properties of the polylactide (PLA) bulk. The three bonded CG potentials are first parametrized against the strain energies of the dimer, trimer, and tetramer, and the nonbonded CG potentials are then optimized to match the melt densities of the decamer. With the derived CG potentials, molecular dynamics (MD) simulations are found to reproduce thermal expansion and glass transition. By rescaling the dihedral and nonbonded potentials with temperature-independent factors, the glass transition temperature (Tg) is also satisfactorily restored with little modifications on the volumetric expansive coefficients at both rubbery and glassy states. Therefore, the finally optimized CG potentials exhibit excellent temperature transferability, as rationalized by the Simha-Boyer relation. Furthermore, it is confirmed that the dihedral torsions and nonbonded interactions play key roles in glass transition. Also, the simulated bulk moduli and conformational properties in a wide temperature range compare well with the referenced data. The proposed multiscale scheme has great potential in simulating thermo-mechanical properties of PLA and other polymers.

A Study on the Effect of Doping Metallic Nanoparticles on Fracture Properties of Polylactic Acid Nanofibres via Molecular Dynamics Simulation

All-atom molecular dynamics simulations are conducted to elucidate the fracture mechanism of polylactic acid nanofibres doped with metallic nanoparticles. Extensional deformation is applied on polymer nanofibres decorated with spherical silver nanoparticles on the surface layer. In the obtained stress-strain curve, the elastic, yield, strain softening and fracture regions are recognized, where mechanical parameters are evaluated by tracking the stress, strain energy and geometrical evolutions. The energy release rate during crack propagation, which is a crucial factor in fracture mechanics, is calculated. The results show that the presence of doping nanoparticles improves the fracture properties of the polymer nanofibre consistently with experimental observation. The nanoparticles bind together polymer chains on the surface layer, which hinders crack initiation and propagation. The effect of the distribution of nanoparticles is studied through different doping decorations. Additionally, a discussion on the variation of internal energy components during uniaxial tensile loading is provided to unravel the deformation mechanism of nanoparticle-doped nanofibres.

Machine-Learning Enabled New Insights into the Coil-to-Globule Transition of Thermosensitive Polymers Using a Coarse-Grained Model

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-isopropylacrylamide) (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.