Hydrophobic Solutes Research Articles

Temperature Dependence of Hydrotropy.

The solubility of hydrophobic solutes increases dramatically with the temperature when hydrotropes are added to water. In this paper, the mechanism of this well-known observation will be explained via statistical thermodynamics through (i) enhanced enthalpy-hydrotrope number correlation locally (around the solute) that promotes the temperature dependence and (ii) hydrotrope self-association in the bulk solution that suppresses the temperature dependence. The contribution from (i), demonstrated to be dominant for urea as a hydrotrope, signifies the weakening of interaction energies around the solute (local) than in the bulk that accompanies incoming hydrotrope molecules. Thus, studying hydrotropic solubilization along the temperature and hydrotrope concentration provides complementary information on the local-bulk difference: the local accumulation of hydrotropes around the solute, driven by the enhanced local hydrotrope self-association, is also accompanied by the overall local weakening of energetic interactions, reflecting the fluctuational nature of hydrotrope association and the mediating role of water molecules.

Retention factor variation on wide range of mobile phase compositions in reversed-phase high-performance liquid chromatography; a short tutorial

Mobile phase composition remains the major experimental parameter influencing the separation process in reversed-phase high-performance liquid chromatography (RP-HPLC). Its influence on the retention of solutes is explained by the Hildebrand solubility parameter and solvophobic theory, which also considers molecular characteristics of the participants in a separation process. The main empirical dependences (linear and polynomial) between the retention factor and the composition of mobile phase (retention function) are discussed taking into consideration a wide range of mobile phase compositions. Two direct chromatographic indices are resulting from these equations: the extrapolated value of retention factor for zero organic content in mobile phase composition (log kw) and the solvent strength parameter (S), which for a linear retention function is its slope being constant over the entire mobile phase compositions. Another chromatographic index could result from the retention function for zero content of water in mobile phase, but this situation can be applied only for very hydrophobic solutes. For large domains of mobile phase composition, the retention function cannot be described by a single type of dependence. This explains the paradox of generating by extrapolation different logkwvalues for different organic components. For the mobile phase close to pure water the retention function is different from that at average water concentration.

The role of the water contact layer on hydration and transport at solid/liquid interfaces

Understanding the structure in the nanoscopic region of water that is in direct contact with solid surfaces, so-called contact layer, is key to quantifying macroscopic properties that are of interest to e.g. catalysis, ice nucleation, nanofluidics, gas adsorption, and sensing. We explore the structure of the water contact layer on various technologically relevant solid surfaces, namely graphene, MoS[Formula: see text], Au(111), Au(100), Pt(111), and Pt(100), which have been previously hampered by time and length scale limitations of ab initio approaches or force field inaccuracies, by means of molecular dynamics simulations based on ab initio machine learning potentials built using an active learning scheme. Our results reveal that the in-plane intermolecular correlations of the water contact layer vary greatly among different systems: Whereas the contact layer on graphene and on Au(111) is predominantly homogeneous and isotropic, it is inhomogeneous and anisotropic on MoS[Formula: see text], on Au(100), and on the Pt surfaces, where it additionally forms two distinct sublayers. We apply hydrodynamics and the theory of the hydrophobic effect, to relate the energy corrugation and the characteristic length-scales of the contact layer with wetting, slippage, the hydration of small hydrophobic solutes and diffusio-osmotic transport. Thus, this work provides a microscopic picture of the water contact layer and links it to macroscopic properties of liquid/solid interfaces that are measured experimentally and that are relevant to wetting, hydrophobic solvation, nanofluidics, and osmotic transport.

Slip Opacity and Fast Osmotic Transport of Hydrophobes at Aqueous Interfaces with Two-Dimensional Materials.

We investigate the interfacial transport of water and hydrophobic solutes on van der Waals bilayers and heterostructures formed by stacking graphene, hBN, and MoS2 using extensive ab initio molecular dynamics simulations. We compute water slippage and the diffusio-osmotic transport coefficient of hydrophobic particles at the interface by combining hydrodynamics and the theory of the hydrophobic effect. We find that slippage is dominated by the layer that is in direct contact with water and only marginally altered by the second layer, leading to a so-called "slip opacity". The screening of the lateral forces, where the liquid does not feel the forces coming from the second nearest layer, is one of the factors leading to the "slip opacity" in our systems. The diffusio-osmotic transport of small hydrophobes (with a radius below 2.5 Å) is also affected by the slip opacity, being dramatically enhanced by slippage. Furthermore, the direction of diffusio-osmotic flow is controlled by the solute size, with the flow in the opposite direction of the concentration gradient for smaller hydrophobes, and vice versa for larger ones. We connect our findings to the wetting properties of two-dimensional materials, and we propose that slippage and wetting can be controlled separately: whereas the slippage is mostly determined by the layer in closer proximity to water, wetting can be finely tuned by stacking different two-dimensional materials. Our study advances the computational design of two-dimensional materials and van der Waals heterostructures, enabling precise control over wetting and slippage properties for applications in coatings and water purification membranes.

Role of Deep Eutectic Solvent Precursors as Hydrotropes: Unveiling Synergism/Antagonism for Enhanced Kraft Lignin Dissolution.

Lignin holds significant potential as a feedstock for generating valuable aromatic compounds, fuels, and functional materials. However, achieving this potential requires the development of effective dissolution methods. Previous works have demonstrated the remarkable capability of hydrotropes to enhance the aqueous solubility of lignin, an amphiphilic macromolecule. Notably, deep eutectic solvents (DESs) have exhibited hydrotropic behavior, significantly increasing the aqueous solubility of hydrophobic solutes, making them attractive options for lignin dissolution. This study aimed at exploring the influence of hydrogen bond acceptors (HBAs) and hydrogen bond donors (HBDs) on the performance of DESs as hydrotropes for lignin dissolution, while possible dissolution mechanisms in different water/DES compositions were discussed. The capacity of six alcohols (glycerol, ethylene glycol, 1,3-propanediol, 1,4-butanediol, 1,5-pentanediol, and 1,6-hexanediol) and cholinium chloride to enhance the solubility of Kraft lignin in aqueous media was investigated. A correlation between solubility enhancement and the alkyl chain length of the alcohol was observed. This was rationalized upon the competition between hydrotrope-hydrotrope and solute-hydrotrope aggregates with the latter being maximized for 1,4-butanediol. Interestingly, the hydrotropic effect of DESs on lignin solubility is well represented by the independent sum of the dissolving contributions from the corresponding HBAs and HBDs in the diluted region. Conversely, in the concentrated region, the solubility of lignin for a certain hydrotrope concentration was always found to be higher for the pure hydrotropes rather than their combined HBA/HBD counterparts.

Interactions of adsorbing cosolutes with hydrophobic hydration shells.

The analysis of water density fluctuations in the hydration shell of nonpolar solutes provides insights into water-mediated interactions, especially hydrophobic interactions. These fluctuations are sensitive to small perturbations due to changes in thermodynamic conditions, such as temperature and pressure, but also to the presence of cosolutes, such as salts or small organic molecules. Herein, we investigate the effect of two classes of adsorbing cosolutes, using urea and methanol as representatives, on the fluctuations in energy and solvent density within the solvation shell of a model extended hydrophobic solute. We focus on the interactions of the cosolutes with the hydrophobic hydration shell, rather than with the solute itself, which though important remain largely unexplored. We calculate and analyze the interfacial partial molar energy of the cosolute, using a methodology based on the small system method. This approach provides correlated solvent density and energy fluctuations and allows us to decompose them into contributions due to interactions between the different components present in the solvation shell of the solute. The results show that adsorbed urea molecules interact more favorably with water than nonadsorbed urea molecules, which leads to the attenuation of interfacial density fluctuations and thus to the stabilization of the solvation shell. By contrast, the adsorbed methanol molecules interact preferably with other methanol molecules in the solvation shell, leading to a nano-phase segregated structure, which enhances interfacial fluctuations.

What is the best simulation approach for measuring local density fluctuations near solvo-/hydrophobes?

Measurements of local density fluctuations are crucial to characterizing the interfacial properties of equilibrium fluids. A specific case that has been well-explored involves the heightened compressibility of water near hydrophobic entities. Commonly, a spatial profile of local fluctuation strength is constructed from the measurements of the mean and variance of solvent particle number fluctuations in a set of contiguous subvolumes of the system adjacent to the solvo-/hydrophobe. An alternative measure proposed by Evans and Stewart [J. Phys.: Condens. Matter 27, 194111 (2015)] defines a local compressibility profile in terms of the chemical potential derivative of the spatial number density profile. Using Grand canonical Monte Carlo simulation, we compare and contrast the efficacy of these two approaches for a Lennard-Jones solvent at spherical and planar solvophobic interfaces and SPC/E water at a hydrophobic spherical solute. Our principal findings are as follows: (i) the local compressibility profile χ(r) of Evans and Stewart is considerably more sensitive to variations in the strength of local density fluctuations than the spatial fluctuation profile F(r) and can resolve much more detailed structure; and (ii) while the local compressibility profile is essentially independent of the choice of spatial discretization used to construct the profile, the spatial fluctuation profile exhibits a strong systematic dependence on the size of the subvolumes on which the profile is defined. We clarify the origin and nature of this finite-size effect.

S.O.S: Shape, orientation, and size tune solvation in electrocatalysis.

Current models to understand the reactivity of metal/aqueous interfaces in electrochemistry, e.g., volcano plots, are based on the adsorption free energies of reactants and products, which are often small hydrophobic molecules (such as in CO2 and N2 reduction). Calculations played a major role in the quantification and comprehension of these free energies in terms of the interactions that the reactive species form with the surface. However, solvation free energies also come into play in two ways: (i) by modulating the adsorption free energy together with solute-surface interactions, as the solute has to penetrate the water adlayer in contact with the surface and get partially desolvated (which costs free energy); (ii) by regulating transport across the interface, i.e., the free energy profile from the bulk to the interface, which is strongly non-monotonic due to the unique nature of metal/aqueous interfaces. Here, we use constant potential molecular dynamics to study the solvation contributions, and we uncover huge effects of the shape and orientation (on top of the already known size effect) of small hydrophobic and amphiphilic solutes on their adsorption free energy. We propose a minimal theoretical model, the S.O.S. model, that accounts for size, orientation, and shape effects. These novel aspects are rationalized by recasting the concepts at the base of the Lum-Chandler-Weeks theory of hydrophobic solvation (for small solutes in the so-called volume-dominated regime) into a layer-by-layer form, where the properties of each interfacial region close to the metal are explicitly taken into account.

Organic micropollutant removal and phosphate recovery by polyelectrolyte multilayer membranes: Impact of buildup interactions

Polyelectrolyte multilayer (PEM) deposition conditions can favorably or adversely affect the membrane filtration performance of various pollutants. Although pH and ionic strength have been proven to alter the characteristics of PEM, their role in determining the buildup interactions that control filtration efficacy has not yet been conclusively proved. A PEM constructed using electrostatic or non-electrostatic interactions from controlled deposition of a weak polyelectrolyte could retain both charged and uncharged pollutants from water. The fundamental relationship between polyelectrolyte charge density, PEM buildup interaction, and filtration performance was explored using a weak-strong electrolyte pair consisting of branching poly (ethyleneimine) and poly (styrene sulfonate) (PSS) across pH ranges of 4–10 and NaCl concentrations of 0 M–0.5 M. PEI/PSS multilayers at acidic pH were dominated by electrostatic interactions, which favored the selective removal of a charged solute, phosphate over chloride, while at alkaline pH, non-electrostatic interactions dominated, which favored the removal of oxybenzone (OXY), a neutral hydrophobic solute. The key factor determining these interactions was the charge density of PEI, which is controlled by pH and ionic strength of the deposition solutions. These findings indicate that the control of buildup interactions can largely influence the physico-chemical and transport characteristics of PEM membranes.

Permeability of Dermatological Solutes through the Short Periodicity Phase of Human Stratum Corneum Lipid Bilayers.

Determining the permeability of drug-like solutes through the densely packed and heterogeneous stratum corneum lipid layer presents a significant challenge. In this study, we employed umbrella sampling with a periodic weighing function applied to the center of mass of the lipid bilayers. Precise umbrella sampling was conducted with an interframe distance of 0.5 Å, spanning from the bilayer center to the water phase, and each frame was simulated for at least 20 ns. Autocorrelation functions, potential of mean force (PMF), and diffusivity profiles were analyzed for six solutes (testosterone, benzene, caffeine, ethanol, mannitol, and histidine). The results revealed that autocorrelations were dependent solely on the medium, whether water or lipid phase. Diffusivity and PMF profiles along the reaction coordinate were influenced by the hydrophilicity of the solute rather than its size. For hydrophobic solutes, the PMF curves exhibited a minimum at the bilayer center, while for hydrophilic solutes, the PMFs peaked at the bilayer center and lipid tails (where the lipid tails are not interacting with the cholesterol). Diffusivity curves were low at the bilayer center and water phase, with peaks observed at the headgroup or the boundary between fatty acid and cholesterol (1 nm from the bilayer center). The quantitative findings presented in this work hold significance for pharmacists and drug designers.

Water Fraction Dependence of the Aggregation Behavior of Hydrophobic Fluorescent Solutes in Water-Tetrahydrofuran.

This work investigates the water fraction dependence of the aggregation behavior of hydrophobic solutes in water-tetrahydrofuran (THF) and the elucidation of the role of THF using fluorescence microscopy, dynamic light scattering, neutron and X-ray scattering, and photoluminescence measurements. On the basis of the obtained results, the following model is proposed: hydrophobic molecules are molecularly dispersed in the low-water-content region (10-20 vol %), while they form mesoscopic particles upon increasing the water fraction to ∼30 vol %. This abrupt change is due to the composition fluctuation of the water-THF binary system to form hydrophobic areas in THF, followed by THF-rich droplets where hydrophobic solutes are incorporated and form loose aggregates. Further increasing the water content prompts the desolvation of THF, which decreases the particle size and generates tight aggregates of solute molecules. This model is consistent with the luminescence behavior of the solutes and will be helpful to control the aggregation state of hydrophobic solutes in various applications.

Diphenylsulfoxide exhibits an outstanding solubility in Methanol: Physicochemical and quantum chemical evidences. Comparison to acetonitrile

Diphenylsulfoxide (Ph2SO) represents a noteworthy organic substance composed of bulky amphiphilic molecules. The chemical interest in Ph2SO arises from a peculiar combination of the strongly polar and strongly nonpolar moieties which are uniquely packed. The solubility of Ph2SO in conventional organic solvents is a cornerstone of its prospective applications. In the present work, the full concentration ranges of the possible homogeneous solutions of diphenylsulfoxide (Ph2SO) in acetonitrile (ACN) and methanol (MEOH) have been investigated by densitometry, UV–vis, FTIR spectroscopy, and quantum-chemical calculations. For the first time, we unraveled an outstanding solubility of Ph2SO in MEOH. The homogeneous solutions form up to the molality of 0.9 mol kg−1. This observation is unexpected from the conventional chemical wisdom because Ph2SO contains twenty-two non-polar interaction sites vs. two polar interaction sites. Based on the volumetric measurements, we revealed that the partial molar volumes of Ph2SO depend somewhat on solvent (MEOH, ACN) peculiarities. The effects of solvent on the spectral characteristics of diphenylsulfoxide have been considered. The numerical simulations based on the systematic potential landscape exploration corroborate the existence of several types of supramolecular species involving Ph2SO and solvent molecules simultaneously. The peculiar anisotropic Ph2SO-MEOH coupling, which is a strong hydrogen bonding, along with the small size of the solvent molecules is responsible for the recorded outstanding solubility in the methanol media. In turn, the performance of ACN is somewhat less encouraging notwithstanding the way higher dipole moment of its molecules. As an amphiphilic solute, Ph2SO may boost the solubilities of hydrophobic solutes by interfacing its coupling with polar media. The reported results inspire further explorations of Ph2SO to design task-specific liquid systems.

Extraction of Polyphenols from Slovenian Hop (Humulus lupulus L.) Aurora Variety Using Deep Eutectic Solvents: Choice of the Extraction Method vs. Structure of the Solvent.

Polyphenols from Slovenian hops (Humulus lupulus L.) of the Aurora variety were extracted by different methods and using classical solvents and several deep eutectic solvents (DES) based on choline chloride as the hydrogen bond acceptor component. The obtained extract solutions were analyzed by HPLC for the content of extracted α- and β-acids and extracted xanthohumol. It was found that choline chloride:phenol DES concentrated aqueous solution had an extraction efficiency close to that of diethyl ether, which is considered one of the best classical extraction solvents for polyphenols from hops. The comparison of the extraction efficiency with other choline chloride-based DESs showed that the chemical similarity of the phenol ring in the solvent DES with the polyphenols in hops may be crucial for a highly efficient extraction with choline chloride:phenol DES. On the other hand, the choice of extraction method and the viscosity of the solvents tested seem to play only a minor role in this respect. As far as we know, this is the first study to attempt to relate extraction efficiency in the extraction of hydrophobic solutes to the compressibility of the DES extractants, the latter of which may be correlated with the extent of hydrophobic hydration around the DES components. In addition, using the heating and stirring method for the preparation of choline chloride-based DES concentrated aqueous solutions we found no support for the occurrence of water in two different roles (in the structural and in the dilution role) in these solvents.

A Heteronuclear NMR Study of Aqueous Lithium Salt Solutions of l-Alanine: Revealing Solute Hydrophobic Association through the NMR B' Coefficient.

In the present study, a set of heteronuclear NMR approaches has been adopted to investigate the solution behavior of a small hydrophobic solute l-alanine in the presence of lithium (Li) salts. The presence of salts plays a major role in determining the structure and solvation of biomolecules. It therefore evokes interest to understand the effect of Li salts on amino acids (alanine), the building block of biomolecules. The ionic solute dynamics in the present case has been probed using 1H, 7Li, and 13C nuclei available in the aqueous Li salt solution of l-alanine. Nuclear longitudinal spin relaxation of alanine protons was examined at a variable concentration range of three lithium salts, i.e., LiCl, Li2SO4 and LiClO4, to introduce the NMR B' coefficient for each salt defining ionic solute/solvent interaction in the solution. Analysis of the active relaxation mechanism of 7Li spin-lattice relaxation further revealed the presence of alanine in the solvation shell of Li ion depending on the anionic counterpart.

Aggregating fluids under centrifugal fields

We present a general theoretical framework designed to interpret the effect of centrifugation of ternary fluids that exhibit weak aggregation. The model includes a buoyancy term with thermodynamic description given by a Flory-Huggins model of ternary systems containing a hydrotrope. The system octanol-ethanol-water used as a numerical example is one of several ternary systems that exhibit dynamic aggregates even far from the critical point and undergo a spontaneous emulsification upon dilution with water (Ouzo effect). The model proposes an explanation from first principles of an increased responsiveness to the centrifugal fields. We illustrate the usage of this system in intensified extraction efficiency of hydrophobic solutes that is at the basis of the dispersive liquid-liquid extraction, and show that its capabilities can be extended from analytical to continuous preparative extraction.

Acidic Conditions Impact Hydrophobe Transfer across the Oil-Water Interface in Unusual Ways.

Molecular dynamics simulation and enhanced free energy sampling are used to study hydrophobic solute transfer across the water-oil interface with explicit consideration of the effect of different electrolytes: hydronium cation (hydrated excess proton) and sodium cation, both with chloride counterions (i.e., dissociated acid and salt, HCl and NaCl). With the Multistate Empirical Valence Bond (MS-EVB) methodology, we find that, surprisingly, hydronium can to a certain degree stabilize the hydrophobic solute, neopentane, in the aqueous phase and including at the oil-water interface. At the same time, the sodium cation tends to "salt out" the hydrophobic solute in the expected fashion. When it comes to the solvation structure of the hydrophobic solute in the acidic conditions, hydronium shows an affinity to the hydrophobic solute, as suggested by the radial distribution functions (RDFs). Upon consideration of this interfacial effect, we find that the solvation structure of the hydrophobic solute varies at different distances from the oil-liquid interface due to a competition between the bulk oil phase and the hydrophobic solute phase. Together with an observed orientational preference of the hydroniums and the lifetime of water molecules in the first solvation shell of neopentane, we conclude that hydronium stabilizes to a certain degree the dispersal of neopentane in the aqueous phase and eliminates any salting out effect in the acid solution; i.e., the hydronium acts like a surfactant. The present molecular dynamics study provides new insight into the hydrophobic solute transfer across the water-oil interface process, including for acid and salt solutions.

EPISOL: A software package with expanded functions to perform 3D-RISM calculations for the solvation of chemical and biological molecules.

Integral equation theory (IET) provides an effective solvation model for chemical and biological systems that balances computational efficiency and accuracy. We present a new software package, the expanded package for IET-based solvation (EPISOL), that performs 3D-reference interaction site model (3D-RISM) calculations to obtain the solvation structure and free energies of solute molecules in different solvents. In EPISOL, we have implemented 22 different closures, multiple free energy functionals, and new variations of 3D-RISM theory, including the recent hydrophobicity-induced density inhomogeneity (HI) theory for hydrophobic solutes and ion-dipole correction (IDC) theory for negatively charged solutes. To speed up the convergence and enhance the stability of the self-consistent iterations, we have introduced several numerical schemes in EPISOL, including a newly developed dynamic mixing approach. We show that these schemes have significantly reduced the failure rate of 3D-RISM calculations compared to AMBER-RISM software. EPISOL consists of both a user-friendly graphic interface and a kernel library that allows users to call its routines and adapt them to other programs. EPISOL is compatible with the force-field and coordinate files from both AMBER and GROMACS simulation packages. Moreover, EPISOL is equipped with an internal memory control to efficiently manage the use of physical memory, making it suitable for performing calculations on large biomolecules. We demonstrate that EPISOL can efficiently and accurately calculate solvation density distributions around various solute molecules (including a protein chaperone consisting of 120,715 atoms) and obtain solvent free energy for a wide range of organic compounds. We expect that EPISOL can be widely applied as a solvation model for chemical and biological systems. EPISOL is available at https://github.com/EPISOLrelease/EPISOL.

Understanding the physics of hydrophobic solvation.

Simulations of water near extended hydrophobic spherical solutes have revealed the presence of a region of depleted density and accompanying enhanced density fluctuations. The physical origin of both phenomena has remained somewhat obscure. We investigate these effects employing a mesoscopic binding potential analysis, classical density functional theory (DFT) calculations for a simple Lennard-Jones solvent, and Grand Canonical Monte Carlo (GCMC) simulations of a monatomic water (mw) model. We argue that the density depletion and enhanced fluctuations are near-critical phenomena. Specifically, we show that they can be viewed as remnants of the critical drying surface phase transition that occurs at bulk liquid-vapor coexistence in the macroscopic planar limit, i.e., as the solute radius Rs → ∞. Focusing on the radial density profile ρ(r) and a sensitive spatial measure of fluctuations, the local compressibility profile χ(r), our binding potential analysis provides explicit predictions for the manner in which the key features of ρ(r) and χ(r) scale with Rs, the strength of solute-water attraction ɛsf, and the deviation from liquid-vapor coexistence of the chemical potential, δμ. These scaling predictions are confirmed by our DFT calculations and GCMC simulations. As such, our theory provides a firm basis for understanding the physics of hydrophobic solvation.

Water model for hydrophobic cavities: structure and energy from quantum-chemical calculations.

This ab initio study aims to design a series of large water clusters having a hollow clathrate-like cage able to host hydrophobic solutes of various sizes. Starting from the (H2O)n (n = 18, 20, 24 and 28) hollow cages, water layers have been added in a stepwise manner in order to model the configuration of water molecules beyond the primary shell. The large (H2O)100, (H2O)120 and (H2O)140 clusters complete the hydrogen bonding network of the cage with optimal and regular tiling of the do-, tetra-decahedron and hexa-decahedron, respectively. This study is corroborated by an investigation of dense water clusters up to the (H2O)123 one, being highly consistent with experimental data on ice concerning the electronic and zero-point energies for aggregate formation at 0 K and enthalpy and entropy at 273 K. The cavity creation profoundly alters the orientation of water molecules compared with those found in dense clusters. Nevertheless, such a large reorganization is necessary to maximize the water-water attraction by making it similar to the one found in dense clusters. The cage formation is an endothermic process; however, the computed values are large compared with previous reports for hydrocarbon aqueous solutions. Larger clusters are required for a more fruitful comparison.

Shedding light on the formation and stability of mesostructures in ternary “Ouzo” mixtures

Hypothesis: Ternary systems made of water, a water-miscible solvent, and hydrophobic solutes spontaneously produce metastable particles by the “Ouzo effect” and thermodynamically stable “Surfactant-Free Micro Emulsions” (SFME). However, the use of different analyses has led to a variability in the criteria to determine the boundaries of the Ouzo domain. We hypothesized that this could be clarified by investigating the stability and the physical state of the particles.Experiments: We investigate four systems using both solid and liquid solutes and two different solvents, and achieved a careful investigation of their phase diagrams, using DLS, Nanoparticle Tracking Analysis, NMR, Multiple Light Scattering, electrophoretic mobility, and fluorescence analysis.Findings: Our results evidence that the transition from the monophasic to the Ouzo domains does not coincide with the cloudiness curve, and that compositions in the Ouzo domain can look fully transparent, in contrast to what is often considered. This transition is best determined by stability analysis. The cloudiness curve corresponds to the formation of particles with a large size dispersity. In the Ouzo domain, we observed an exchange of solute between the continuous phase and solute particles swollen with solvent. In addition, the particles are stabilized against coalescence by their high negative charge.