Mathematical and computational modeling of fats and triacylglycerides

Fats and oils are found in many food products; however, their macroscopic properties are difficult to predict, especially when blending different fats or oils together. With difficulties in sourcing specific fats or oils, whether due to availability or pricing, food companies may be required to find alternative sources for these ingredients, with possible differences in ingredient performance. Mathematical and computational modeling of these ingredients can provide a quick way to predict their properties, avoiding costly trials or manufacturing problems, while, most importantly, keeping the consumers happy. This review covers a range of mathematical models for triacylglycerides (TAGs) and fats, namely, models for the prediction of melting point, solid fat content, and crystallization temperature and composition. There are a number of models that have been designed for both TAGs and fats and which have been shown to agree very well with empirical measurements, using both kinetic and thermodynamic approaches, with models for TAGs being used to, in turn, predict fat properties. The last section describes computational models to simulate the behavior of TAGs using molecular dynamics (MD). Simulation of TAGs using MD, however, is still at an early stage, although the most recent papers on this topic are bringing this area up to speed.

Comparative study of natural wax-based W/O emulsion gels: Microstructure and macroscopic properties

In this paper, six kinds of natural wax, including sunflower wax (SFX), rice bran wax (RBX), carnauba Brazilian wax (CBX), beeswax (BWX), candelilla wax (CDX), and sugarcane wax (SGX) were used to prepare water-in-oil (W/O) emulsion gels. Microstructures and rheological properties of all emulsion gels were investigated by microscopy, confocal laser scanning microscope (CLSM), scanning electron microscopy (SEM), and rheometer, respectively. By comparing polarized light images of wax-based emulsion gels and corresponding wax-based oleogels, it could be found that dispersed water droplets greatly affected the crystal distribution and hindered crystal growth. Polarized light microscopy and CLSM images proved that natural wax could perform a dual-stabilization mechanism by interfacial crystallization and crystal networks. SEM images illustrated all waxes except SGX were platelets and formed networks by stacking on top of each other, while flocs-like SGX was easier to adsorb on the interface and formed a "crystalline shell". The surface area and pore formed by different wax varied wildly, which accounted for their differences in the gelation ability, oil binding capacity, and strength of the crystal network. The rheological study showed that all wax had solid-like properties and wax-based oleogels with denser crystal networks correspond to emulsion gels with higher modules. The dense crystal network and interfacial crystallization could improve the stability of W/O emulsion gels proved by recovery rates and critical strain. All the above proved that natural wax-based emulsion gels can be used as stable, low-fat, and thermal-sensitive fat mimics.

Reproduction of macroscopic properties of unsaturated triacylglycerides using a modified NERD force field

Unsaturated triacylglycerides are found in many commonly consumed foods, such as cooking oils, nuts and chocolate. There are however very few publications on Molecular Dynamics simulations of such molecules, and, to the best of our knowledge, no such published research on crystalline mono-unsaturated triacylglycerides. The work described in this paper is an evaluation of different force fields (GROMOS96 and NERD) to determine the best force field parameters to reproduce the crystalline and melted macroscopic properties of such molecules accurately. The best results were obtained by modifying the NERD force field, through which we were able to reproduce the crystalline and melted density as well as crystal dimensions of mono-unsaturated triacylglycerides.

Effects of wax concentration and carbon chain length on the structural modification of fat crystals

Natural waxes are cost-effective and potential fat crystallization modifiers; however, there is limited information about their implementation in solid triacylglycerol (TAG) oil. Herein, we investigated the effects of two natural waxes, namely, candelilla wax (CLW) and rice bran wax (RW), with different concentrations (2, 4, 6, and 8 wt%) and carbon chains on the crystal growth and structure of palm kernel stearin (PKS85). CLW significantly accelerated the PKS85 crystallization process. Both waxes could induce a new hydrocarbon chain with the lengths of 3.70 and 4.15 Å during the TAG crystallization, respectively. Particularly, X-ray diffraction (XRD) indicated that PKS85 combined with CLW showed a similar lamellar thickness (d₀₀₁) and crystal domain size (ξ) with pure PKS85, whereas that of PKS85 containing RW was 1.7-1.8 and 1.5-1.8 fold higher, respectively. This result corresponded to the carbon chain length of CLW and RW, which was double and quadruple that of PKS85, respectively. Further, these variations were reflected in the crystal microstructures of PKS85 with CLW and RW, where the former showed small homogeneous crystals, while the latter displayed large rod-like layered crystals. In addition, the firmness significantly increased when CLW and RW were added, which is possibly attributed to the fact that the waxes became the backbone of the crystal "fence". Our findings give clear insight into the interaction between TAGs and wax molecules in the crystallization process, which can help guide the utilization of natural waxes in the modification of fats.

Ultra small angle x-ray scattering in complex mixtures of triacylglycerols

Ultra-small angle x-ray scattering (USAXS) has been used to elucidate, in situ, the aggregation structure of unsheared model edible oils. Each system comprised one or two solid lipids and a combination of liquid lipids. The 3D nano- to micro-structures of each system were characterized. The length scale investigated, using the Bonse-Hart camera at beamline ID-15D at the Advanced Photon Source, ANL, ranged from 300 Å-10 µm. Using the Unified Fit model, level-1 analysis showed that the scatterers were 2D objects with either a smooth, a rough, or a diffuse surface. These 2D objects had an average radius of gyration Rg1 between 200-1500 Å. Level-2 analysis displayed a slope between -1 and -2. Use of the Guinier-Porod model gave s ≈ 1 thus showing that it was cylinders (TAGwoods) aggregating with fractal dimension 1 ≤ D2 ≤ 2. D2 = 1 is consistent with 1D structures formed from TAGwoods, while D2 = 2 implies that the TAGwoods had formed structures characteristic of diffusion or reaction limited cluster-cluster aggregation (DLCA/RLCA). These aggregates exhibited radii of gyration, Rg2, between 2500 and 6500 Å. Level-3 analyses showed diffuse surfaces, for most of the systems. These interpretations are in accord with theoretical models which studied crystalline nano-platelets (CNPs) coated with nano-scale layers arising from phase separation at the CNP surfaces. These layers could be due to either liquid-liquid phase separation with the CNPs coated, uniformly or non-uniformly, by a diffuse layer of TAGs, or solid-liquid phase separation with the CNPs coated by a rough layer of crystallites.A fundamental understanding of the self-organizing structures arising in these systems helps advance the characterization of fat crystal networks from nanometres to micrometres. This research can be used to design novel fat structures that use healthier fats via nano- and meso-scale structural engineering.