Upon cooling, alkanes can form intermediate phases between liquid and crystal. They are called "rotator" or "plastic" phases and have long-range positional order with rotational freedom around the long molecular axis which gives them non-trivial and useful visco-plastic properties. We expect that the formation and structure of rotator phases formed in freezing alkanes can be understood much deeper by tracking the process at molecular level with atomistic molecular dynamics. We defined an appropriate CHARMM36-based computational protocol for simulating the freezing of hexadecane, which contained a sufficiently long (500 ns) equilibrium sampling of the frozen states. We employed it to simulate successfully the freezing of bulk and interface-contacting hexadecane and to provide a pioneering clarification of the effect of surfactant on the crystallization mechanism and on the type of intermolecular ordering in the crystallites. The devised computational protocol was able to reproduce the experimentally observed polycrystalline structure formed upon cooling. However, different crystallization mechanisms were established for the two types of models. Crystallites nucleate at random locations in the bulk and start growing rapidly within tens of nanoseconds. In contrast, the surfactants freeze first during the fast cooling (<1 ns), followed by rapid hexadecane freezing, with nucleation starting along the entire surfactant adsorption layer. Thereby, the hexadecane molecules form rotator phases which transition into a more stable ordered phase. This collective transition is first-time visualized directly. The developed robust computational protocol creates a foundation for future in-depth modelling and analysis of solid-state alkane-containing, incl. lipid, structures.

TiO₂-based materials are of great practical interest in several technological areas. Both the size and the morphology of the TiO₂ particles are of critical importance for their applications. The current study explores the effect of several factors on the outcome of the TiO₂ particle synthesis via the so-called ‘two-emulsion method’. In this technique, two water-in-oil emulsions—each of them containing different reactant in the dispersed water drops—are mixed under well controlled conditions. Upon such mixing, partial coalescence of the water drops from the two emulsions leads to mixing of the drop content, with chemical reaction occurring within the drops, and to synthesis of Ti(OH)₄ particles. Afterwards, the latter are transformed by emulsion heating into TiO₂ particles and aggregates of predominantly anatase structure. Our results show that—depending on the precursor and surfactant concentrations, oil viscosity, emulsification time, and mixing speed—the obtained nanoparticles could aggregate either on the drop surface, forming capsules with a very smooth surface, or inside the water droplets, thus leading to hierarchically structured aggregates of micrometer size. The spherical smooth capsules are constructed of very small monodisperse TiO₂ nanoparticles with size below 5 nm. The hierarchical bulk aggregates, on the other hand, are formed from bigger primary particles of sub-micrometer size. The obtained results show that one can obtain various TiO₂ structures by controlling the conditions during the emulsion preparation and mixing

Lipid-based formulations (LBF) enhance oral drug absorption by promoting drug solubilization and supersaturation. The aim of the study was to determine the effect of the lipid carrier type, drop size and surfactant concentration on the rate of fenofibrate release in a bicarbonate-based in vitro digestion model. The effect of the lipid carrier was studied by preparing type I LBF with drop size ≈ 2 µm, based on medium-chain triglycerides (MCT), sunflower oil (SFO), coconut oil (CNO) and cocoa butter (CB). The drop size and surfactant concentration effects were assessed by studying MCT and SFO-based formulations with a drop size between 400 nm and 14 µm and surfactant concentrations of 1 or 10%. A filtration through a 200 nm filter followed by HPLC analysis was used to determine the aqueous fenofibrate, whereas lipid digestion was followed by gas chromatography. Shorter-chain triglycerides were key in promoting a faster drug release. The fenofibrate release from long-chain triglyceride formulations (SFO, CNO and CB) was governed by solubilization and was enhanced at a smaller droplet size and higher surfactant concentration. In contrast, supersaturation was observed after the digestion of MCT emulsions. In this case, a smaller drop size and higher surfactant had negative effects: lower peak fenofibrate concentrations and a faster onset of precipitation were observed. The study provides new mechanistic insights on drug solubilization and supersaturation after LBF digestion, and may support the development of new in silico prediction models.

The saponin escin, extracted from horse chestnut seeds, forms adsorption layers with high viscoelasticity and low gas permeability. Upon deformation, escin adsorption layers often feature surface wrinkles with characteristic wavelength. In previous studies, we investigated the origin of this behavior and found that the substantial surface elasticity of escin layers may be related to a specific combination of short-, medium-, and long-range attractive forces, leading to tight molecular packing in the layers. In the current study, we performed atomistic molecular dynamics simulations of 441 escin molecules in a dense adsorption layer with an area per molecule of 0.49 nm². We found that the surfactant molecules are less submerged in water and adopt a more upright position when compared to the characteristics determined in our previous simulations with much smaller molecular models. The number of neighbouring molecules and their local orientation, however, remain similar in the different-size models. To maintain their preferred mutual orientation, the escin molecules segregate into well-ordered domains and spontaneously form wrinkled layers. The same specific interactions (H-bonds, dipole–dipole attraction, and intermediate strong attraction) define the complex internal structure and the undulations of the layers. The analysis of the layer properties reveals a characteristic wrinkle wavelength related to the surface lateral dimensions, in qualitative agreement with the phenomenological description of thin elastic sheets.

In the current work, a comprehensive procedure for structural analysis of quasilinear organic molecules arranged in a polycrystalline sample generated by molecular dynamics is developed. A linear alkane, hexadecane, is used as a test case because of its interesting behavior upon cooling. Instead of a direct transition from isotropic liquid to the solid crystalline phase, this compound forms first a short-lived intermediate state known as a “rotator phase”. The rotator phase and the crystalline one are distinguished by a set of structural parameters. We propose a robust methodology to evaluate the type of ordered phase obtained after a liquid-to-solid phase transition in a polycrystalline assembly. The analysis starts with the identification and separation of the individual crystallites. Then, the eigenplane of each of them is fit and the tilt angle of the molecules relative to it is computed. The average area per molecule and the distance to the nearest neighbors are estimated by a 2D Voronoi tessellation. The orientation of the molecules with respect to each other is quantified by visualization of the second molecular principal axis. The suggested procedure may be applied to different quasilinear organic compounds in the solid state and to various data compiled in a trajectory.