Can Polarity-Inverted Surfactants Self-Assemble in Nonpolar Solvents?

We investigate the self-assembly process of a surfactant with inverted polarity in water and cyclohexane using both all-atom and coarse-grained hybrid particle-field molecular dynamics simulations. Unlike conventional surfactants, the molecule under study, proposed in a recent experiment, is formed by a rigid and compact hydrophobic adamantane moiety, and a long and floppy triethylene glycol tail. In water, we report the formation of stable inverted micelles with the adamantane heads grouping together into a hydrophobic core and the tails forming hydrogen bonds with water. By contrast, microsecond simulations do not provide evidence of stable micelle formation in cyclohexane. Validating the computational results by comparison with experimental diffusion constant and small-angle X-ray scattering intensity, we show that at laboratory thermodynamic conditions the mixture resides in the supercritical region of the phase diagram, where aggregated and free surfactant states coexist in solution. Our simulations also provide indications as to how to escape this region to produce thermodynamically stable micellar aggregates.


INTRODUCTION
Life as we know it could not exist without water. In fact, living cells survive in environments mainly constituted by water. Cellular shape and functionality are determined by the presence of both the plasma and the cytoplasmic membrane, which define all of the necessary compartments for the organization of cellular matter and prevent the mixing of the cell with its external environment. To this aim, living organisms typically exploit biological lipids, amphiphile molecules comprising a strongly polar head group and one or more long hydrocarbon tails. 1 In aqueous solutions, these amphiphilic molecules tend to aggregate driven by "like-tolike" interactions that are usually referred to as the hydrophobic effect. 2 Within this general framework, water is unique because it forms hydrogen bonds with itself as well as with the polar moiety of the amphiphilic molecule. Hydrogen bonds play a special intermediate role as they have a strength of the order of 10−40 kJ/mol (corresponding to 5−10 k B T/bond at 298 K), much stronger than van der Waals interactions (≈1 kJ/mol) and considerably weaker than ionic or covalent bonds (≈100 kJ/mol or more). Also, hydrogen bonds have an intermediate orientation-dependence that is in between the strongly directional covalent and isotropic van der Waals interactions.
While this marvelous balance is the result of millions of years of evolution, it is possible to imagine that a similar outcome could be achieved in different biological environments under different conditions, such as those present in other planets of our universe. Although water has been detected in various thermodynamic states in our solar system, an alternative scenario suggests the possibility of using polarity-inverted membranes in nonpolar solvents, such as the hydrocarbons frequently found in earth-like systems (see the recent review by Sandstrom 3 ). Motivated by this idea, a number of related studies have recently been conducted. Pace and collaborators investigated protein stability in a nonaqueous solvent such as cyclohexane (C 6 H 12 ). 4 Hayashi et al. found that while proteins have well-defined unique structures in water, this is generally not the case in other nonpolar solvents. 5,6 The stabilities of single polar and hydrophobic amino acids in water and nonpolar solvents have also been studied (unpublished results). To close the triangle of life, the stability of B-DNA under nonaqueous conditions has also been recently assessed. 7 Notwithstanding the large number of studies that have been addressing the issue, the basic mechanism underlying a solvophobic effect in nonpolar solvents is still far from being fully understood. In its simplest terms, it could be stated as follows. If the polarity of amphiphilic molecules is inverted, so to have a hydrophobic (rather than polar) head and a polar (rather than hydrophobic) tail, would they self-assemble in nonpolar solvents, for example, C 6 H 12 ? If so, what would the driving force be?
Recently such an experiment has been performed on a newly synthesized molecule having exactly these features. 8 This molecule, henceforth referred to as ADOH, is formed by a rigid and compact adamantane moiety AD and a long and floppy tail that consists of a triethylene glycol (TEG) with a characteristic group O−CH 2 −CH 2 capped at the end by a hydroxyl group that is able to form hydrogen bonds ( Figure 1).
The self-assembly properties of ADOH were studied by measuring its diffusion coefficient in C 6 H 12 , using NMR spectroscopy at different concentrations. A monotonic decrease of the diffusion coefficient, which is a possible signal of the micellization process, was observed in the concentration range from 5 to 250 mM and this hypothesis was further supported by small-angle X-ray scattering (SAXS) measurements that appeared to indicate a critical micelle concentration (CMC) around 100 mM. While through the Stokes−Einstein equation this decrease in the diffusion coefficient can be ascribed to the appearance of large aggregates, it alone does not constitute proof of the existence of aggregates with a welldefined micellar shape, especially considering the fact that directional-dependent polar interactions, such as hydrogen bonds, are significantly stronger compared with nondirectional van der Waals interactions. Moreover, while SAXS experiments provide essential system information on properties such as micelle sizes, the interpretation of aggregates in terms of shapes requires postprocessing by modeling, which is prone to errors. Thus, a molecular picture of aggregation is often hard to obtain.
Molecular modeling can complement experiments by providing molecular resolution predictions on the spatial organization of the molecules. In general, surfactant aggregation is a challenging process to simulate as it is facilitated by slow diffusing molecules and typically occurs at very low concentrations requiring both very large system sizes and long simulation times. For this specific ADOH surfactant system, the expected drop in diffusion coefficient occurs at very high concentrations, making it within reach of standard allatom simulations. Nevertheless, all-atom simulations are computationally expensive and, even for this high-concentration regime, they may be affected by significant finite-size effects and may be limited in the range of accessible time scales.
Both of these drawbacks can be effectively addressed using coarse-grained modeling. 9,10 In coarse-grained models, a lower resolution representation of the molecular structure with effective potentials is used to lower the overall computational cost, thereby allowing for the study of larger and longer simulations. For soft matter systems in particular, such as ADOH, the coarse-grained methodology of hybrid particlefield molecular dynamics (hPF-MD) has already been proved to be particularly effective. In hPF-MD, the coarse-grained molecular resolution description is combined with density-field modeling of intermolecular interactions to yield a computationally efficient modeling of very large systems. Applications of hPF-MD have started from more conventional soft polymer mixtures and then have moved also to biological systems. 11−14 Examples from the literature include nanocomposites, nanoparticles, percolation phenomena in carbon nanotubes, 15−18 lamellar and nonlamellar phases of phospholipids, 19−22 and more recently polypeptides, 23 and polyelectrolytes. 24−28 These applications, in particular polyelectrolytic molecules and surfactants (Triton X-100), 29 give us sufficient ingredients for building a hPF-MD model for ADOH.
Thus, the aim of the present work is to provide a molecular understanding of the nature of the putative ADOH aggregates reported in ref 8, as well as the underlying physical driving forces. Using both all-atom and hPF-MD numerical simulations, the self-assembly properties of ADOH molecules will be studied both in C 6 H 12 and in water under the same conditions reported in the experiment. While for the former case this will provide a complementary description with respect to the experiment, the latter case in water represents a new prediction that could eventually lead to further experimental testing.
The rest of the article is organized as follows. Section 2 outlines the all-atom and hPF-MD methods used in the present article. Section 3 reports the results, while Section 3.3 provides a connection to the experimental findings. Finally, Section 4 includes the key messages of the present study as well as some perspectives for future work.

All-Atom
Simulations. NPT all-atom simulations were run using the all-atom optimized potentials for liquid simulations (OPLS-AA) force field 30 with a time step of 2 fs. The temperature was set to 300 K and the pressure to 1 atm. The coupling was ensured by applying the v-rescale thermostat, 31 with a relaxation time of 0.1 ps, and the Parrinello− Rahman barostat, 32 with coupling constant set to 3 ps and isothermal compressibility equal to 4.5 × 10 −5 bar −1 . Longrange electrostatics was calculated with the PME method using a fourth-order interpolation, a 0.16 nm Fourier spacing, and a 1.2 nm cutoff, which was the same also for the calculation of short-range van der Waals interactions. Bond lengths were constrained using the LINCS algorithm. 33 The water model used was the TIP4P, 34 while for C 6 H 12 we used the parameterization implemented in the OPLS-AA force field. 35 The all-atom simulations were run with the GROMACS 2018.4 software. 36 2.2. hPF-MD Approach. In hPF-MD, molecular dynamics is used to sample the phase space of a fully resolved molecular system composed of N mol molecules with Hamiltonian ] is the interaction energy that depends on the particle density ϕ.
To model nonbonded attraction and repulsion between particles, we employ the following interaction energy 37 where ϕ 0 is the total number density, χĩ j is the interaction term between species i and j, ϕ i (r) and ϕ j (r) are the number densities of the ith and jth species calculated at positions r, and κ is a compressibility term. The net effect of is an external potential V i acting on all particles of type i, which is obtained by The force acting on a particle of type i is obtained by gradient operation on V Calculations of the potentials and of the forces acting on the particles, and those used to integrate the equations of motion, are computed with a particle mesh approach. For more details, we refer to ref 38. Figure 1 shows the CG mapping chosen for ADOH, while in Table 1 we report the bead interaction matrix χĩ j used in the present work. The χṽ alues for C 6 H 12 and ADOH were selected from chemically similar moieties of Triton X-100 from ref 29.

hPF-MD Simulations.
CG simulations were run in the NVT ensemble using the OCCAM software. 38 The time step was set to 30 fs and the temperature was kept constant at 300 K by applying the Anderson thermostat 39 with a collision frequency of 7 ps −1 . The density field was updated every 20 time steps. In this model, each C 6 H 12 molecule is represented by a single bead, while a water bead comprises four molecules. For the simulations in C 6 H 12 , we added small alternating partial charges (+0.4 and −0.4, chosen equal as the all-atom charges assigned to ether oxygens) on ADOH tail beads. This was done to qualitatively mimic the weak electrostatic interactions that, from atomistic simulations of the system with and without partial charges on the TEG segment, were found crucial in describing the clustering of the surfactant in the apolar solvent. The electrostatic potential in this case follows the Poisson equation with a relative dielectric constant ϵ r = 5 and its computation is performed using the recently developed adopted PME approach as in refs 24, 25. The composition of all simulated systems is reported in Table 2.

RESULTS AND DISCUSSION
3.1. ADOH in Cyclohexane. We start our analysis by presenting simulation results for ADOH in C 6 H 12 , i.e., the same system studied experimentally in ref 8. In this experiment, results from SAXS and NMR suggested the onset of micelle formation, with a CMC of approximately 100 mM and an average radius estimated between 1.7 and 2.5 nm. As the contour length of ADOH is ≈1.5 nm, the emerging scenario was that of nearly spherical micelles with all of the TEG tails buried deep inside each micelle and the adamantane hydrocarbon heads in contact with the solvent. Hence, the analogue of a conventional spherical micelle with reversed polarity of both the surfactant and the solvent.
Spherical micelles can indeed be expected on the basis of the surfactant packing parameter 40 where V t and l t are the tail volume and length, respectively, and A c is the contact area between the head group and the tail. In the case of ADOH l 0.93 1.25 t = − nm, V t ≈ 0.603 nm 3 and A c ≈ 1.55 nm 2 , as obtained by approximating the head group to a sphere and assuming the contact area as half of the sphere surface. This leads to N s = 0.312 − 0.421, consistent with a surfactant forming spherical or cylindrical aggregates. We remark, however, that this argument assumes a tight packing in a straight conformation of the TEG tails into the core of the micelle, which contradicts the characteristic high flexibility associated with the chemical structure of the tail.
We first consider all-atom simulations to have full control over the driving forces at the microscopic level. The simulations were run at 200 mM, a concentration significantly higher than the putative experimental CMC at which an increase of the SAXS intensity was observed (≈100 mM). Figure 2A shows a snapshot of the final configuration after 450 ns, and no clear sign of any sort of aggregation is visible. Clearly, this may be an issue of the atomistic simulation time scale. Indeed, even in conventional surfactants, although early aggregation already occurs in the first few nanoseconds of simulation, the time scales for the stabilization of micelles usually extend to the microseconds regime and can be most suitably probed by coarse-grained models, such as hPF-MD. The absence of stable aggregates is however confirmed by hPF-MD, as depicted by the corresponding snapshot of Figure 2B obtained after 1800 ns. Zooming in on the all-atom snapshot of Figure 2A, it is possible to see a seven-molecule aggregate, two dimers kept together by tail−tail interactions of the hydroxyl groups, and two free monomers ( Figure 3A). This represents a typical transient molecular cluster that is frequently observed The Journal of Physical Chemistry B pubs.acs.org/JPCB Article during the simulation characterized by an irregular shape and a very short lifetime in the 10 1 −10 2 ps range. The absence of stable micelles at this concentration is confirmed by the practically flat time profile of the solvent-accessible surface area (SASA) obtained from all-atom simulations ( Figure 3B), indicating the absence of a core-collapse. The dynamical behavior of the system (see the movie provided in the SI) clearly indicates the presence of fast forming/disrupting dimers, trimers, and higher-order multimers. Nonetheless, aggregation does not seem to result in the formation of a   The Journal of Physical Chemistry B pubs.acs.org/JPCB Article typical core, ADOH diffusion appears dominated by the monomeric phase, and the nonpolar solvent does not show any preferential affinity for either the head or the tail group. In particular, the first peak of the radial distribution function (RDF) of C 6 H 12 against the adamantane head or the TEG tail atoms is found for both cases at ∼0.6 nm, while the profile shows the characteristic modulations typical of a good solvent ( Figure 4A,B). Moreover, the RDF plot does not change over time, indicating the absence of progressive aggregation during the simulation. The all-atom ADOH−ADOH center of mass RDF reported in Figure 3C shows a broad peak at ≈0.75 nm and tends to reach unity at a distance of r ≈ 2 nm, which is an indication of the presence of uncorrelated disorder beyond this range. Also shown in Figure 3C is the RDF from hPF-MD calculation that provides a consistent picture although with a peak at a slightly larger value, indicating the tendency for the ADOH molecules to settle at this distance on a longer time scale.
The absence of well-defined aggregates for ADOH in C 6 H 12 is finally confirmed by the cluster size distribution ( Figure 3D), which exhibits an exponential decay with cluster size N aggr . The average cluster size remains constant during the simulation around a value of ≈7, indicating that a steady state has been reached. This is further corroborated by hPF-MD simulations that do not show any formation of stable micelles even after 1800 ns (see Figure 2B) and that show the same exponential decay found for the all-atom case in the aggregation number analysis.

ADOH in Water.
To better understand the physics of ADOH in solution, it proves interesting to investigate the selfassembly of ADOH in water, a case that unfortunately was not covered in the experiment of ref 8. Figure 5A depicts the aggregation state for a 50 mM concentration of ADOH in water after only 50 ns of all-atom MD simulations. Even a very short simulation time window is sufficient to report evidence that ADOH self-assembles in regular spherical micelles, which is in agreement with the predictions given by its packing parameter of 0.312 and its simple molecular structure ( Figure  1). The bulky hydrophobic adamantane heads of the amphiphilic molecule promote fast aggregation to minimize the contact with water, grouping together into a well recognizable hydrophobic core, while the hydroxyl groups of the TEG tails stick out into the solvent. This is also well depicted by the time evolution of the RDF profile between the ADOH moieties and the solvent ( Figure 4C,D). In particular, at the end of the simulation, RDF profiles reach bulk values at much longer distances than in C 6 H 12 , a consequence of the micellar collapse that screens ADOH moieties from the solvent. There is also a clear differentiation between the RDF profiles of TEG oxygens and adamantane, indicating the expected clear preference of the solvent for the former. The opposite behavior instead does not occur in C 6 H 12 .
The formation of micelles at this concentration is confirmed by the analysis of the SASA that, unlike in C 6 H 12 , steadily decreases and levels off toward a stable value after 50 ns of MD simulations ( Figure 5B). This is further supported by the ADOH−ADOH center of mass RDF reported in Figure 5C. A very broad peak at ≈0.75 nm indicates the smallest ADOH− ADOH distance in a dimer (see Figure 5A), while a slow decaying signal that asymptotically reaches unity at 3 nm is indicative of the several additional characteristic ADOH− ADOH distances present in the micelle, all consistent signs of the presence of stable aggregates. The cluster size distribution analysis for the micelles ( Figure 5D; here, the bin size is five units) corroborates this picture with the visible presence of a multivalued distribution and a relatively large polydispersity in micelle size, which is mainly due to the relatively short simulation times. In fact, the plot of the average cluster size as a function of time is characterized by a gradual increase in the value, indicating that the final equilibrium state is not reached.
In any case, the all-atom results are sufficient to indicate that the simulated system is well above the critical micelle concentration (CMC). To obtain a rough estimate of the CMC of ADOH in water, we repeated simulations of ADOH at progressively lower concentrations using a hPF-MD approach. Thanks to the cheaper potentials and the intrinsically accelerated dynamics, hPF-MD yields a better convergence of the aggregation state of ADOH even at much lower concentrations than 50 mM.
The hPF-MD model was first validated against reference allatom data by running simulations at 50 mM. Figure 6 presents snapshots of the final configurations obtained after 100 ns of hPF-MD simulations at 10, 50, 100, and 200 mM concentrations. These snapshots provide clear indications of aggregation, also supported by the corresponding ADOH− ADOH RDF (not shown), all exhibiting a peak localized at 0.75 nm and a similar long-range decay (see Figure 5C). Note that the discrepancies on the height of the peak and higher distribution values at short range are a typical feature of hPF models and must be attributed to the soft nature of the hPF potential.
From the morphological viewpoint, the size of the aggregates increases with solute concentration and at 200 mM, micelles start fusing together to form tubular structures. This is consistent with the concentration-dependent smooth transition from spherical to tubular micelles that is commonly observed in more conventional amphiphilic surfactants and can be ascribed to the reaching of a critical packing of the hydrophobic moieties, above which they can no longer be accommodated into a compact spherical volume. 29,41,42 The absence of significant aggregated units at 10 mM suggests that, in water, the CMC for ADOH is between 10 and 50 mM. We repeated additional hPF-MD simulations at progressively higher ADOH dilutions (1, 2, 5, and 15 mM). By performing a linear fit of the concentration dependence of the free monomer fraction on the surfactant concentration at smaller and higher values than the CMC and determining their intersection point, we produce the best estimate for the CMC ≈ 13.5 mM (Figure 7). Here, we notice that the number of surfactant molecules present in the lowest concentration simulation is quite small, hence the value of the first point in Figure 7 could be lower if we accounted for possible finite-size errors. However, even changing the monomer fraction in the range 0.9−0.7, the error in the CMC is only ∼1−4%.
3.3. Reconciling with the Experiment. The interpretation of the experimental results of ADOH in C 6 H 12 in terms of formation of well-defined micelles above a CMC ≈ 100 mM 8 and the numerical simulations that found no evidence of a stable micellization process are in apparent contradiction.
In fact, formation of aggregates was originally suggested on the basis of two pieces of experimental evidence: (i) a significantly lower diffusion constant for ADOH compared to that of the pure solvent and (ii) a marked change in the SAXS intensity signal for concentrations above 100 mM.
Diffusion Coefficient. Diffusion coefficients for C 6 H 12 , water, and ADOH in the two solvents were estimated from the slope of the mean square deviation (MSD), σ 2 (t), calculated for all-atom simulations according to the Einstein diffusion equation, σ 2 (t) = 6Dt + A, and compared with experimental results obtained from 1 H two-dimensional diffusion-ordered spectroscopy (2D-DOSY) NMR measurements. Table 3 reports the comparison between experimental diffusion coefficients (first column) and the corresponding allatom estimates (third column) and provides evidence of a reasonably good agreement. The results of the second column will be discussed further below. On this basis, we now argue that the drop observed in the experiment might be due to other reasons and not related to the micellization process.
SAXS Spectrum. Figure 8 shows the experimental SAXS data for a 200 mM concentration of ADOH in C 6 H 12 taken from ref 8. Interestingly, the line is not characterized by a marked drop of the intensity as expected for regular spherical objects of a well-defined radius. Rather, at a higher scattering vector Q, it decays slowly and without showing any particular feature, a behavior compatible with the presence of irregular objects with no clear size. This is in very good qualitative agreement with our simulations, which show the absence of a well-defined organization of the surfactant. The molecular RDF, in particular, is characterized by a simple profile with a short-range peak and a fast decay ( Figure 3C). Assuming a similar situation in the experiment, with an exponentially decaying RDF, the spectral SAXS line would then have the following form where ξ is the Ornstein−Zernike correlation length. 44 As can be observed in Figure 8, such line shape fits well the experimental data. For comparison, we report the predicted SAXS spectrum obtained by Fourier transform of the RDF from the all-atom simulation, as well as its fitting by eq 6. Like for the diffusion coefficient, the agreement with the fit on the SAXS measurements is an indication that our computational model is well in line with the experimental findings.
Driving Forces of ADOH Interactions in C 6 H 12 . The standard understanding of micelle formation is based on the hydrophobic effect, which is an entropy-driven phenomenon. By contrast, in C 6 H 12 , the formation of ADOH clusters is likely to be penalized entropically due to the necessary confinement of TEG tails in the micellar core. Rather, ADOH aggregates might be stabilized by the enthalpy gain associated with the electrostatic interactions occurring among the polar TEG    To verify this hypothesis, we ran additional all-atom simulations of ADOH in both water and C 6 H 12 where all partial charges on the TEG segment of the amphiphile were set to zero and replaced by an effective dipole in the final hydroxyl group mimicking the total dipole of the tail. These artificial systems keep a coarse representation of the polarity of the TEG tail while at the same time removing the quadrupolar charge distributions characteristic of the glycol ether moieties. The presence of a finite dipole on the terminal OH also ensures a correct representation of the terminal H-bonding group, which was hypothesized to be crucial for ADOH aggregation in C 6 H 12 .
Simulations in water, discussed in Section 3.2, resulted in a faster formation of micelles compared to what was observed for the physical model (not shown). In C 6 H 12 , the absence of quadrupolar charges on the TEG segment reduces even further the aggregation of ADOH, resulting in a net 30% increase of its diffusion coefficient (see the second column in Table 3).
The final snapshot of this system in C 6 H 12 after 150 ns is shown in Figure 9A, while the SASA is reported in Figure 9B. Clear differences in the aggregation behavior are evident by inspecting the RDF between ADOH molecules ( Figure 9C). In the absence of partial charges on the TEG segment, the RDF does not show a clear first-neighbor peak. On the contrary, the RDF reaches unity immediately at r = 1 nm, which implies that each molecule is spatially uncorrelated from one another and that the distribution of surfactant in the sample is practically uniform. The distribution of the aggregation number ( Figure 9D) shows a marked difference with the distribution more skewed toward the monomer peak rather than larger aggregates and an average cluster size of N aggr = 4 significantly smaller than the original model with all of the partial charges (compare with the inset of Figure 3D).
Formation of hydrogen bonding between the hydroxyl groups at the tails of ADOH was also putatively indicated as a possible driving force for ADOH aggregation in C 6 H 12 . Figure  10 reports the number of hydrogen bonds per ADOH molecule (N H /N ADOH ) as a function of time in water and C 6 H 12 . As expected, the number of H-bonds between ADOH in water is marginal, as the strongly hydrophilic groups prefer to interact with the solvent (see Figure 10A). However, in C 6 H 12 , a significant fraction of ADOH is involved in H-bonds, with a very stable average of ≈0.25 H-bonds per molecule, as shown in Figure 10B. This number is anyway much smaller    Figure  3D), one would expect a value of ∼0.87 H-bonds per molecule, which is almost three times the one observed. This finding confirms that while H-bonds between ADOH tails can contribute to the binding interactions, they themselves are not strong enough to constitute a significant driving force for ADOH aggregation.
Overall, simulations of the artificial system without TEG charges confirm that cohesion of ADOH in C 6 H 12 is an enthalpy-dominated phenomenon mainly due to quadrupole− quadrupole interactions between the TEG polar regions. This is qualitatively different from the situation in water, where the collapse of the surfactant is determined by the entropy and where the same qualitative effect could be obtained by lumping all partial charges of TEG into an effective dipole on the hydroxyl group at the end of the tail. Our simulations in C 6 H 12 indicate the presence of small, labile ADOH clusters rapidly forming and dissolving in the sample. The presence of an appreciable density dishomogeneity in the absence of true stable aggregates is indicative of a mixture at supercritical conditions, confirmed by the shape of the SAXS profile, in qualitative agreement with those observed in other supercritical mixtures, for example, water/acetonitrile. 45

CONCLUSIONS
We performed a detailed computational analysis of the selfassembly process of ADOH surfactants composed of a hydrophobic adamantane head and a TEG tail terminating with an OH hydroxyl group. We have studied the mechanism driving the formation of aggregates in C 6 H 12 by comparing their sizes, distributions, interactions, as well as their diffusion coefficients with their experimental counterparts. 8 In addition, we have also studied the same system in water providing a qualitative prediction for all of the above quantities that could be experimentally tested. By a detailed analysis of the results under different concentrations, we have also been able to estimate a value of 13.5 mM for the CMC in water.
Our parallel analysis in water and C 6 H 12 has allowed us to underpin the different forces driving the self-assembly in the two solvents. In water, the conventional hydrophobic effect plays a major role in promoting the aggregation of the hydrophobic heads that are buried inside the micelles to avoid direct contact with water molecules. In C 6 H 12 , an analogous lipophobic effect does not occur, and formation of well-defined micelles with an inverted structure was not observed. The lack of such an effect may be attributed to concomitant factors, including the fact that TEG chains are not as lipophobic as adamantane is hydrophobic and that the conformational frustration of the flexible TEG moieties is insufficiently balanced by the solvent entropy release associated with the collapse of an aggregate.
Our investigation highlighted the fact that the aggregation of ADOH into low-weight oligomers is sufficient to explain the decrease of diffusion coefficient experimentally observed, as well as the previously reported SAXS data. Formation of these oligomers is determined by short-range attraction of the electrostatic quadrupoles distributed over the TEG segment. By contrast, in water, these interactions are favored by the low dielectric screening due to the weakly polarizable solvent.
However promising, ADOH in C 6 H 12 appears to linger in a supercritical phase without a clear collapse into well-defined self-assembled aggregates. Our study points to multiple possible routes that could produce significant steps toward such a goal. First, there is the need of decreasing the entropy of the free TEG chain to reduce the entropy loss upon confinement of the polar regions in the micellar core. This could be achieved by narrowing the accessible conformational space of the TEG chain already in the monomer, for example, by adding a second TEG unit to the ADOH structure. Additionally, it might be possible to chemically modify the TEG part or to use other polymers having an even stronger polar character. While these two options aim for raising the critical temperature of the systems, it is worth mentioning the even simpler idea of studying ADOH dissolved in low-freezingpoint solvents. This option should be considered in particular thinking of other thermodynamic regimes (high pressure, low temperature) that could be found in extraterrestrial environments. Synthetic effort in this direction may soon lead to the determination of surfactants having the ability of inverting their aggregation structure in solvents of radically different polarities.
We provide a short movie for the all-atom simulation of ADOH in C 6 H 12 to show the fast aggregation/ disgregation process (MPG)