Prof. Slavka S. Tcholakova, Ph.D.
Head of the Laboratory for Active Formulations and Materials
Interests
- Formation, stability and rheology of foams and emulsions
- Mechanism of action of antifoams
- Surfactant aggregates and mesophases in bulk: structure, rheology, stability
- Biophysics & colloid science in fat digestion and drug delivery systems
- Defoamers and antifoams
- Physicochemical theoretical models for foamability, emulsification; foam rheology
- Natural based surfactants
Bio
Slavka Tcholakova received M.Sc. in Chemistry (1996), Ph.D. in Physical Chemistry (2004), Assistant Professor (2006), Associate Professor (2009), Professor (2013) in the Faculty of Chemistry and Pharmacy, Sofia University, Bulgaria. She is head of the Department of chemical and pharmaceutical engineering in Sofia University since January 2015. She has been a visiting researcher in the Research Center Paul Pascal, CNRS, Bordeaux, France (1997). Her research interests include: formation and stability of emulsions; rheology of foams and emulsions; protein adsorption in relation to emulsion stability; foam stability in the presence of antifoams; in-vitro models for digestion and bioavailability of hydrophobic molecules and the respective experimental methods. So far, she has published 129 research and review articles, cited over 4600 times in the scientific literature (h-index = 38). She has been leading 55 projects and participating in 30 projects with international companies (BASF, Unilever, Saint Gobain, Wacker, Lubrizol, PepsiCo, Altana, BYK Chemie, Productolysa, etc.). She has been a supervisor and co-supervisor of 14 completed PhD Theses, and 6 other Theses are under preparation. She was the recipient of the award “Best Young Scientist” for 2006 of the Sofia University Foundation “St. Kliment Ohridski”. She bears the National award “Pythagoras” (2018) for high scientific achievements in natural sciences.
Publications
Most recent publications
Emulsification in nearly Newtonian and non-Newtonian media of wormlike micelles
Emulsification experiments with four silicone oils, having viscosities ranging from 0.01 to 30 Pa.s, were conducted in two types of media: nearly Newtonian polyvinyl alcohol (PVA) solutions and non-Newtonian mixtures induced by worm-like micelles in solutions of sodium laureth sulfate and cocoamidopropyl betaine (BS) with NaCl. The increased viscosity of BS solutions upon the addition of NaCl did not significantly affect the drop size in the formed emulsions. In contrast, the increased viscosity of solutions with higher PVA concentrations significantly reduced the drop sizes for all silicone oils. A theoretical expression predicting the maximum drop size in both types of media (nearly Newtonian and non-Newtonian) was derived and validated against experimental data. The expression accounts for shear-thinning behavior in both the aqueous and oil phases. Interfacial stress dominates the breakage of less viscous oils, while viscous stress inside the breaking drop plays a leading role for more viscous oils. The formation of emulsions with similar sizes in non-Newtonian solutions of BS with different NaCl concentrations was explained by their strongly shear-thinning behavior, which leads to nearly similar viscosity at high shear rates, despite their zero-shear viscosities differing by more than two orders of magnitude.
Alkyl sucrose esters vs. Brijs: How chain length and temperature impact surface and foam properties
The primary objective of this study is to determine the similarities and differences in the surface, film, and foam properties of alkyl sucrose esters (SEs) with high monoester content (≥ 70 %) and polyoxyethylene alkyl ethers (Brijs) with a high number of ethoxy groups (≥ 20) in their head group. Experiments were conducted using surfactant molecules with alkyl chain lengths of 12, 16, and 18 carbon atoms at concentrations between 0.01 and 1 wt%, within a temperature range of 25 °C to 60 °C.
The lag time for surfactant adsorption increased with surfactant chain length, decreased with temperature, and significantly decreased with surfactant concentration for both types of studied surfactants. However, increasing the chain length from 12 to 18 carbon atoms led to a 10-fold increase in lag time for Brijs and more than a 600-fold increase for SEs. This effect rendered longer-chain SEs incapable of forming voluminous foam in Bartsch test and led to pronounced coalescence between the bubbles after their separation from the sparger in the foam rise method, resulting in foams with very large bubbles, which exhibited lower stability. The utilization of a Kenwood mixer for foam generation provided sufficient time for longer-chain SE molecules to adsorb on the bubble surfaces and to produce voluminous foams with small bubbles, which remained stable even at 60 °C. In contrast, foams generated from Brijs solutions are very unstable at 60 °C. The long-standing stability of SEs foam was attributed to the formation of mixed mono- and diesters adsorption layers on the bubble surfaces.
Mechanisms of dissolution and crystallization of amorphous glibenclamide
Amorphous solid dispersions enhance the dissolution and oral bioavailability of poorly water-soluble drugs. However, the link between polymer properties and formulation performance has not been fully clarified yet. We studied the effect of hydroxypropyl cellulose (HPC) polymers molecular weight (Mw) on the storage stability, dissolution kinetics and supersaturation stability of spray-dried amorphous glibenclamide (GLB) formulations. The solid-state stability of amorphous GLB during storage was significantly enhanced by both the 40 kDa (HPC-SSL) and 84 kDa (HPC-L) polymers, regardless of Mw differences. In contrast, HPC-SSL maintained significantly higher aqueous drug concentrations during dissolution, compared to HPC-L (its higher Mw analogue). Dedicated dissolution experiments, in situ optical microscopy and solid-state characterization revealed that aqueous drug concentrations were determined by the interplay between crystallization inhibition, drug ionization, wetting and solubilization effects: (1) HPC prevents surface nucleation, hence inhibiting crystallization, (2) intestinal colloids (bile salts and phospholipids) increase supersaturated drug concentrations via wetting and solubilization effects and (3) pH and drug ionization severely impact the degree of supersaturation. The better performance of the lower Mw HPC-SSL was due to its superior inhibition of surface crystallization during dissolution. These insights into the molecular mechanisms of dissolution and crystallization of amorphous solids provide foundation for rational formulation development.
Sucrose ester surfactants: Current understanding and emerging perspectives
Sucrose esters (SEs), derived from sucrose and fatty acids, are biodegradable and non-toxic surfactants increasingly favored as substitutes for petrochemically synthesized ones in food, cosmetics, and pharmaceuticals. SEs provide versatile hydrophilic–lipophilic properties, determined by the degree of sucrose esterification ranging from one to eight. The length of the fatty acid residues further influences the phase behavior of SEs, allowing creation of tailored formulations for specific applications. This review provides insights about our current understanding of the SEs phase behavior, their aggregation in aqueous and oily solutions, and its correlation with formulation outcomes. Furthermore, an overview of recent studies investigating SEs in various colloidal systems, including emulsions, foams, oleogels, and others, is provided. Novel concepts are discussed alongside future research directions, emphasizing the SEs potential as sustainable, functional ingredients.
Types of phases obtained by molecular dynamics simulations upon freezing of hexadecane-containing systems
Medium- to long-chain alkanes can form upon cooling intermediate phases between isotropic liquid and solid crystalline, called rotator phases, where relative freedom of the molecules to rotate about their long axis is combined with long range translational order. Rotator phases are well documented experimentally but the mechanism of their formation at the molecular level is still not fully explained. In a previous work, we have shown that molecular dynamics simulations can produce rotator phases upon cooling of hexadecane [S. Iliev et al., J. Col. Int. Sci., 2023, 638, 743]. The aim of the current work is to develop a procedure to identify the specific ordered phase obtained in the simulations. The influence of the cooling rate on the freezing process of hexadecane (bulk and surfactant-interfaced to water) is tested as well. Several parameters are combined to quantify the degree of ordering and the type of phase in the studied systems. These are the tilt angle of the molecules with respect to the crystallite plane, the radial distribution function of the centre of mass of the molecules in the crystallite, the percentage of the gauche torsion angles in the molecules, the angle of the second principal axis of each molecule with respect to the x axis of the coordinate system, and estimates from Voronoi analysis. The results show that the systems form a rotator phase, which transitions gradually towards the thermodynamically most stable triclinic crystal, and the transformation progresses to different extent depending on the system. The influence of the cooling rate is related only to the size of the largest crystallite formed, the other parameters of the freezing process remain unaffected. The work also presents a robust procedure for obtaining and identifying different types of ordered phases in alkane-containing systems with thoroughly tested computational protocol and a comprehensive set of structural analyses. Several key characteristics are advanced, compared to previous research [Ryckaert et al., Mol. Phys., 1989, 67, 957; Wentzel et al., J. Chem. Phys. 2011, 134, 224504], namely, a new methodology is proposed to compute the unit cell deformation parameter and azimuthal angle from MD simulation trajectories of the freezing process in alkane-containing systems. The suggested structural analysis, which is independent of the coordinate system, is applicable to any linear-chain system with polycrystalline structure.