This study is carried out using the COGITO force field to determine whether the thermodynamic melting point of pure triacylglyceride crystals can be predicted using molecular dynamics simulations. The triacylglycerides used in this study are both saturated and unsaturated, as well as symmetrical and asymmetrical, to test the robustness of both the force field and the direct heating methodology described in this paper. Given the nonequilibrium nature of a melting system, a larger number of simulations are required to ensure that the results are sufficiently converged, that is, with little fluctuation and a small confidence interval. The study also highlights the importance of the presence of defects, in this case as voids, to lower the melting nucleation energy barrier of the crystals and avoid superheating of the systems being tested. The size of these defects is much larger than what would be found in a physical crystal, however, the simple and robust procedure that was developed allows the accurate prediction of melting points of the different triacylglycerides.

No abstract

The use of molecular dynamics simulations is becoming ever more widespread; however, the application of this to pure triacylglyceride (TAG) systems is not. In this study, we are presenting the development, and validation, of a new force field (FF), which we have called the COarse-Grained Interchangeable Triacylglyceride-Optimized FF. The FF has been developed using both a bottom-up and top-down approach for different parameters, with the non-bonded parameters being optimized using a Bayesian optimization method. While the FF was developed using monounsaturated TAGs, results show that it is also suitable for fully saturated TAGs. Description of molecules which were not used during the development of the FF is carried out simply by interchanging the bead in the molecule topologies. Results show that the FF can reproduce the macroscopic properties (density and lattice parameters) of pure TAGs as both crystals and melt with high accuracy, as well as reproduce the differences in enthalpies.

The relative stability of crystalline polymorphs and the transition between crystalline and melt phases are key parameters in determining the physical properties of triacylglycerides used in food. However, while the determination of properties experimentally is well-defined, the ability to predict the onset of melting and discriminate between polymorphs is less well-defined within a molecular dynamics simulation environment. In this work, we present metrics for measuring the crystallinity, including a new metric, the near-neighbor occupancy time, giving a rapid determination of how many, and which, molecules are found in a crystal over a simulation trajectory, and the polymorphic determination of triacylglycerides over a simulation trajectory.