Directly Modifying the Nonbonded Potential Based on the Standard Iterative Boltzmann Inversion Method for Coarse-Grained Force Fields

Derived Coarse-Grained Potentials for Semicrystalline Polymers with a Blended Multistate Iterative Boltzmann Inversion Method

In this article, we employ the multistate iterative Boltzmann inversion (MS-IBI) method to develop coarse-grained potentials capable of representing molecular structure in both the amorphous and crystalline phases of semicrystalline polymers with improved accuracy while allowing for tunable control over the dynamics governing the α-relaxation process. A unique feature of this method is that the potentials are blended using the product of the target structural distributions, for example, the radial density function, for each phase and a weighting factor. To demonstrate this approach, a family of potentials for polyethylene is developed where the weighting factor of the crystalline phase ranges is varied from zero, incorporating information only from the amorphous phase, to unity, where the model is trained from only the crystalline phase. The most accurate representation of structural distributions was obtained when the crystalline phases is weighted at 50%. However, we show that when the crystalline phase is weighted at 90%, the model more accurately represents dynamics of the α-relaxation process, with realistic predicted values of activation energy and diffusion rates, with relatively minor impact on accuracy in structure.

Permeability of antioxidants through a lipid bilayer model with coarse-grained simulations

Oxidative stress caused by pollution and lifestyle changes causes an excess of free radicals that react chemically with cell constituents leading to irreversible damage. There are molecules known as antioxidants that reduce the levels of free radicals. Some pigments of fruits and vegetables known as anthocyanins have antioxidant properties. Their interaction with the cell membrane becomes a crucial step in studying these substances. In this research, molecular dynamics simulations, particularly, coarse-grained molecular dynamics (CGMD) were used. This technique aims to replace functional groups with corresponding beads that represent their level of polarity and affinities to other chemical groups. Also, umbrella sampling was carried out to obtain free energy profiles that describe well the orientation and location of antioxidants in a membrane considering Trolox, Cyanidin, Gallic Acid, and Resveratrol molecules to study the structural effects they cause on it. It was concluded in this study that an antioxidant when crossing the membrane does not cause either damage to the structural properties or the loss of packing and stratification of phospholipids. it was also observed that the most reactive part of the molecules could easily approach area A prone to lipid oxidation, which can describe the antioxidant capacity of these molecules.Communicated by Ramaswamy H. Sarma.

Effectiveness of coarse graining degree and speedup on the dynamic properties of homopolymer

Evaluation of effective coarse graining (CG) degree and reasonable speedup relative to all-atomistic (AA) model was conducted to provide a basis for building appropriate larger-scale model. The reproducibility of atomistic conformation and temperature transferability both act as the analysis criteria to resolve the maximum acceptable CG degree. Taking short- and long time spans into account simultaneously in the estimation of computational speedup, a dynamic scaling factor is accessible in fitting mean squared displacement ratio of CG to AA as an exponential function. Computing loss in parallel running is an indispensable component in acceleration, which was also added in the evaluation. Subsequently, a quantified prediction of CG speedup arises as a multiplication of dynamic scaling factor, computing loss, time step, and the square of reduction in the number of degrees of freedom. Polyethylene oxide was adopted as a reference system to execute the direct Boltzmann inversion and iterative Boltzmann inversion. Bonded and non-bonded potentials were calculated in CG models with 1~4 monomers per bead. The effective CG degree was determined as two at the most with a speedup of four orders magnitude over AA in this study. Determination of effectiveness CG degree and the corresponding speedup prediction provide available tools in larger spatiotemporal-scale calculations.