Development of a coarse-grained model for simulations of tridecanoin liquid-solid phase transitions

Accessing the electronic structure of liquid crystalline semiconductors with bottom-up electronic coarse-graining

Understanding the relationship between multiscale morphology and electronic structure is a grand challenge for semiconducting soft materials. Computational studies aimed at characterizing these relationships require the complex integration of quantum-chemical (QC) calculations, all-atom and coarse-grained (CG) molecular dynamics simulations, and back-mapping approaches. However, these methods pose substantial computational challenges that limit their application to the requisite length scales of soft material morphologies. Here, we demonstrate the bottom-up electronic coarse-graining (ECG) of morphology-dependent electronic structure in the liquid-crystal-forming semiconductor, 2-(4-methoxyphenyl)-7-octyl-benzothienobenzothiophene (BTBT). ECG is applied to construct density functional theory (DFT)-accurate valence band Hamiltonians of the isotropic and smectic liquid crystal (LC) phases using only the CG representation of BTBT. By bypassing the atomistic resolution and its prohibitive computational costs, ECG enables the first calculations of the morphology dependence of the electronic structure of charge carriers across LC phases at the ∼20 nm length scale, with robust statistical sampling. Kinetic Monte Carlo (kMC) simulations reveal a strong morphology dependence on zero-field charge mobility among different LC phases as well as the presence of two-molecule charge carriers that act as traps and hinder charge transport. We leverage these results to further evaluate the feasibility of developing mesoscopic, field-based ECG models in future works. The fully CG approach to electronic property predictions in LC semiconductors opens a new computational direction for designing electronic processes in soft materials at their characteristic length scales.

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.

Polymorphic phase transitions in triglycerides and their mixtures studied by SAXS/WAXS techniques: In bulk and in emulsions

Triacylglycerols (TAGs) exhibit a monotropic polymorphism, forming three main polymorphic forms upon crystallization: α, β' and β. The distinct physicochemical properties of these polymorphs, such as melting temperature, subcell lattice structure, mass density, etc., significantly impact the appearance, texture, and long-term stability of a wide range products in the food and cosmetics industries. Additionally, TAGs are also of special interest in the field of controlled drug delivery and sustained release in pharmaceuticals, being a key material in the preparation of solid lipid nanoparticles. The present article outlines our current understanding of TAG phase behavior in both bulk and emulsified systems. While our primary focus are investigations involving monoacid TAGs and their mixtures, we also include illustrative examples with natural TAG oils, highlighting the knowledge transfer from simple to intricate systems. Special attention is given to recent discoveries via X-ray scattering techniques. The main factors influencing TAG polymorphism are discussed, revealing that a higher occurrence of structural defects in the TAG structure always accelerates the rate of the α → β polymorphic transformation. Diverse approaches can be employed based on the specific system: incorporating foreign molecules or solid particles into bulk TAGs, reducing drop size in dispersed systems, or using surfactants that remain fluid during TAG particle crystallization, ensuring the necessary molecular mobility for the polymorphic transformation. Furthermore, we showcase the role of TAG polymorphism on a recently discovered phenomenon: the creation of nanoparticles as small as 20 nm from initial coarse emulsions without any mechanical energy input. This analysis underscores how the broader understanding of the TAG polymorphism can be effectively applied to comprehend and control previously unexplored processes of notable practical importance.

Structural and thermodynamic properties of bulk triglycerides and triglyceride/water mixtures reproduced using a polarizable coarse-grained model

Triglycerides (TGs) play important roles in renewable energies, food production, medicine, and metabolism in organisms. Here, we developed a novel coarse-grained (CG) force field (FF) for triglycerides to reproduce both the structural and thermodynamic properties of bulk TGs, TG/air interfaces, and TG/water mixtures using molecular dynamics (MD) simulations. We rigorously optimized the bonded and nonbonded force parameters between the CG beads of TGs and nonbonded force parameters between TG beads and polarizable CG water beads by employing an efficient meta-multilinear interpolation parameterization algorithm recently developed by us. This CG FF performs very well in reproducing the percolating network of the TG bulk phase self-assembled in water and a variety of molecular conformations predicted by all-atom MD simulations. More importantly, it also correctly reproduces multiple experimentally measurable macroscopic thermodynamic properties, including the density and surface tensions of both the TG/air and TG/water interfaces. This paves the way for studying more complicated systems involving TGs on a large scale.

COGITO: A Coarse-Grained Force Field for the Simulation of Macroscopic Properties of Triacylglycerides

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.

Rapid Automated Quantification of Triacylglyceride Crystallinity in Molecular Dynamics Simulations

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.

Molecular Dynamics Simulation for Mechanism Elucidation of Food Processing and Safety: State of the Art

Molecular dynamics (MD) simulation is a useful technique to study the interaction between molecules and how they are affected by various processes and processing conditions. This review summarizes the application of MD simulations in food processing and safety, with an emphasis on the effects that emerging nonthermal technologies (for example, high hydrostatic pressure, pulsed electric field) have on the molecular and structural characteristics of foods and biomaterials. The advances and potential projection of MD simulations in the science and engineering aspects of food materials are discussed and focused on research work conducted to study the effects of emerging technologies on food components. It is expected by showing key case studies that it will stir novel developments as a valuable tool to study the effects of emerging food technologies on biomaterials. This review is useful to food researchers and the food industry, as well as researchers and practitioners working on flavor and nutraceutical encapsulations, dietary carbohydrate product developments, modified starches, protein engineering, and other novel food applications.

The triolein/aqueous interface and lipase activity studied by spectroscopic ellipsometry and coarse grained simulations

In spite of the importance of the triglyceride aqueous interface for processes like emulsification, surfactant interactions and lipase activity, relatively little is known about this interface compared to that between alkanes and water. Here, the contact between triolein and water was investigated in terms of water inclusion in the oil phase and orientation of the molecules at the interface. Coarse grained models of triglycerides in contact with water were constructed and correlated with experimental results of the changes in thickness and refractive index, obtained using spectroscopic ellipsometry of spin-coated triolein films. The topography of the layer was revealed by atomic force microscopy. Dry triolein and a triolein sample after equilibration with water were also compared structurally using small-angle X-ray scattering. Additionally, the kinetics of adsorption/activity of three different variants of the Thermomyces lanuginosus lipase (TLL) were investigated. The results show that uptake of water in the triolein phase leads to increase in thickness of the layer. The observed increase of thickness was further enhanced by an active lipase but reduced when an inactive mutant of the enzyme was applied.

Interdigitation between Triglycerides and Lipids Modulates Surface Properties of Lipid Droplets

Intracellular lipid droplets (LDs) are the main cellular site of metabolic energy storage. Their structure is unique inside the cell, with a core of esterified fatty acids and sterols, mainly triglycerides and sterol esters, surrounded by a single monolayer of phospholipids. Numerous peripheral proteins, including several that were previously associated with intracellular compartments surrounded by a lipid bilayer, have been recently shown to target the surface of LDs, but how they are able to selectively target this organelle remains largely unknown. Here, we use atomistic and coarse-grained molecular dynamics simulations to investigate the molecular properties of the LD surface and to characterize how it differs from that of a lipid bilayer. Our data suggest that although several surface properties are remarkably similar between the two structures, key differences originate from the interdigitation between surface phospholipids and core neutral lipids that occurs in LDs. This property is extremely sensitive to membrane undulations, unlike in lipid bilayers, and it strongly affects both lipid-packing defects and the lateral pressure profile. We observed a marked change in overall surface properties for surface tensions >10 mN/m, indicative of a bimodal behavior. Our simulations provide a comprehensive molecular characterization of the unique surface properties of LDs and suggest how the molecular properties of the surface lipid monolayer can be modulated by the underlying neutral lipids.

Transferability of the coarse-grained potentials for trans-1,4-polybutadiene

The CG torsion potential has a dual effect on the reproduction of crystallization of trans-1,4-polybutadiene in CG simulation.