Microphase separation of a miscible binary liquid mixture under confinement at the nanoscale

Recent experimental works suggested that the confinement into a cylindrical nanopore induced the microphase separation of a binary liquid, despite the miscible character of its bulk counterpart. A core–shell organization was evidenced such that one of the liquids was strongly anchored to the solid surface whereas the other was confined at the center of the pore. At the same time, a study based on atomistic simulations suggested a strong heterogeneity and the absence of a separation. In this work, by refining the solid–liquid interactions to qualitatively reproduce the experimental adsorption isotherms of both single liquids, the microphase separation and the core–shell structure are captured. By tuning the surface chemistry of the nanopore to mimic hydrophilic and hydrophobic confinement, we show that it is possible to control the structural characteristics of the core–shell structure.The molecular origin of the microphase separation is then ascribed to the strong hydrogen bonds and a commensurate arrangement between the confining material and both liquids.


INTRODUCTION
For three decades confinement effects at the nanoscale on the physics of fluids have been been intensively studied. [1][2][3][4][5][6][7][8][9][10][11][12][13][14][15] Many new properties have been then discovered such as the giant diffusion of liquids, [16][17][18] the apparition of new phases and new transitions, [19][20][21][22] giant dielectric properties, 11,12,[23][24][25] the increase of optical properties, 26 and the possible mixing of non-miscible mixtures. 27 These observations indicated that the classical understanding of the physics of liquids should be revisited in confined geometry. Whereas these effects have been largely investigated in the case of confined single components [1][2][3][4][5][6][7][8][9][10][11][12][13][14] and immiscible binary mixtures, [28][29][30] less works have been devoted to the confinement of miscible liquid mixtures. 15,[31][32][33] Recently, Muthulakshmi et al. 34 reported an experimental evidence of a partial phase separation of an ethanol-water mixture confined in mesoporous silica using positron annihilation lifetime spectroscopy. They showed that a small fraction of the ethanol molecules seemed to be anchored at the silica surface. 34 A similar result was also established by Guo et al. 14 who have studied confinement of the ethanol-water mixture between two planar silica walls. A partial and local separation between water and ethanol close to the silica surface was thus evidenced. 14 At the same time, Schmitz et al. 35 have exhibited that the glycol-water mixture could undergo an interfacial separation. This phenomenon was also numerically observed by You et al. 36 who exhibited a local demixing of binary hard-core Yukawa mixtures in a slitlike pore. More recently, Krycka et al. 32 displayed a separation between two confined apolar liquids. At same time, Harrach et al. 33 have studied the phase behavior of a mixture of isobutyric acid (iBA) and water confined in mesoporous SBA-15 silica material and have shown that the iBA-rich phase is close to the pore wall and the water-rich phase is in the center of the pores. Whereas these works only suggest a partial separation of two hydrogen bonds forming liquids or two apolar liquids near the solid surface, Morineau and co-workers 15 have recently provided a direct experimental structural evidence of the microphase separation of macroscopically miscible liquids consisting of hydrogen bonds forming liquid and an apolar one. Indeed, the structure of a mixture comprising toluene (TOL) and tert-butanol (TBA) molecules confined in a cylindrical silica nanopore (MCM-41) of radius 24 Å was explored by neutron scattering and compared with the miscible bulk one. 15 Using a core-shell (CS) model, the authors established, for the first time, a molecular-scale phase-separated tubular structure with the TBA molecules forming a layer at the pore surface (shell), surrounding a TOL-rich phase at the center of the pore (core). This observation was later extended down to low temperature in the glassy states and to the larger pores sizes of SBA-15. 37,38 A consistent picture, showing the highly selective segregation of polar molecules at the pore surface, was deduced from binary gas adsorption experiments. 31 Although these experiments revealed the presence of a microphase separation, the local organization and microscopic processes controlling the CS structure has to be clarified. In a recent numerical work, Essafri et al. 39 showed that the nanoconfinement of TBA/TOL mixture in a silica nanopore only induced strong heterogeneity. This result highlights then a scenario based on the local segregation rather than a core-shell organization of the confined TBA/TOL mixture. As suggested by the authors, the difference between experiment and simulation would be probably due to the hydrophilic character of the membrane that differs between both experimental and numerical studies. Therefore, the core-shell microphase separation of confined mixtures, which was experimentally evidenced by different studies, emerged as an open and challenging question for numerical methods. The current study paves the way to a thorough understanding of the phenomenon from carefully designed simulations. 15,31,37,38 Using atomistic simulations, we aim to capture the CS structure and to clarify the microscopic driving force ruling it. To do so, the solid-liquid interactions were refined to qualitatively reproduce the experimental adsorption isotherm of pure TBA and TOL components confined through the silica nanopore. 31 An illustration of three TBA, TOL, and silica framework is provided in Fig. 1a. This refinement allowed us to yield new insights into the molecular understanding of microphase separation. Moreover, by tuning the surface chemistry of the nanopore to mimic hydrophilic and hydrophobic porous materials, we investigated the role of the surface chemistry on the demixing.

RESULTS
Force field refinement To be in line with experiment, the confining medium corresponds to a cylindrical silica nanopore of radius of 24 Å. The models and the computational procedure are given in the Methods section. Adsorption isotherms were modeled using Monte Carlo simulation in the grand canonical statistical ensemble (GCMC). Liquid-liquid and liquid-solid interactions (U tot ) were taken into account by considering the electrostatic (U elec ) and the van der Waals interactions from the Lennard-Jones (LJ) model (U LJ ) such that where q i is the partial charge of i atom, N the number of particles, r ij the distance between i and j atoms, ε ij the depth of the potential well between i and j, and σ ij refers to the Van der Waals radius. To reproduce the adsorption isotherms of the pure components the solid-liquid interactions were refined. The OPLS (Optimized Potentials for Liquid Simulations) force field was used to model the TOL and TBA molecules. Indeed, it was shown that the thermodynamic properties of TBA/TOL mixtures were closely reproduced from the OPLS model. 39 Whereas the intramolecular contributions (bonds, bending, and dihedral angles) were conserved as original, the partial charges were calculated using ab-initio calculations. Silica material was modeled by considering the ClayFF force field. 40 All details are provided in the Methods section. To refine the LJ parameters ruling the solid-liquid interactions, we began to predict the adsorption isotherm by considering the initial force fields. Concerning the TBA adsorption through nanoporous silica, Fig. 2 shows a qualitative agreement between simulation and experiment that highlights the quality of ClayFF and OPLS force fields and their combining using the Lorentz-Berthelot mixing rules (ε ij ¼ ffiffiffiffiffiffiffi ffi ε ii ε jj p and σ ij = (σ ii + σ jj )/2). As shown in Fig. 2 the predicted adsorption isotherm of toluene strongly differs with experiment. To improve the silica-TOL interactions we adapted the procedure discussed in ref. 41 initially developed to derive the coarse-grained parameters from atomistic simulations. The LJ parameters have been optimized by minimizing the value of the function F ¼ 1 according to a procedure detailed in ref. 41 . In this relation, s i is the estimated statistical uncertainty, f sim i and f exp i are the values of the ith properties on n. In this work n = 2 and corresponds to the enthalpy of adsorption and adsorbed amount at low partial pressure. The minimum condition of F is that every partial derivative must be zero. 41 Optimized LJ parameters are provided in Table 1. Figure 2 shows that the predicted adsorption isotherms are in fair agreement with those obtained from experiments. 31 Interestingly, Fig. 2 exhibits that the capillary condensation pressure is qualitatively reproduced for both TOL and TBA components. These results allowed us to be confident in these refined interactions.  Fig. 1b, for both pure components, a layering structure was observed due to the confinement effect at the nanometric scale. 5 Whereas three layers are observed for both cases, the interfacial TBA layer seems closer to the silica surface, suggesting that the attractive interactions between the silica nanopore and TBA molecules are stronger. As shown in Fig. 3, the hydroxide groups (OH) of the TBA molecules point towards the silica surface. The driving force of this interfacial anchoring is probably the result of the strong hydrogen bonding (HB) network between porous silica and TBA molecules. This point will be discussed later. Interestingly, as shown in Fig. 3, this HB anchoring leads to a peculiar organization such that hydrophobic and hydrophilic (HB) domains are disposed in alternating succession. These results should probably impact the molecular structure of confined mixture. The density profile of the SiOH groups is also reported in Fig. 3. As shown in Fig. 3, the SiOH profile begins from 11.5 Å and presents a peak around 14 Å, highlighting the roughness of the surface. Calculation of the accessible volume by probing the porous volume by an atom of Argon allowed us to approximatively evaluate the pore radius at 12 Å.
From x TBA = 0.24 to 0.83, Fig. 4 (from Fig. 4a-d) shows that TBA molecules were adsorbed close to the interface constituting a first layer, whereas the toluene molecules are located right after. For x TBA = 0.83, the second and third layers correspond to a homogeneous mixture phase reminiscent of the bulk medium. Very interestingly, Fig. 4 highlights the progressive development of two phases as the TBA concentration decreases. Indeed, for x TBA = 0.71, the TBA molecules are almost only located close to the silica surface in the first adsorbed layer without TOL molecules. For x TBA = 0.49 that corresponds to the concentration where the core-shell organization was observed experimentally, the microphase separation is nearly complete. Interestingly, the distance between peaks of TOL and TBA seems to decrease as a function of the decrease in molar fraction in TBA. Indeed, whereas in Fig. 4a, b the distance between the interfacial TOL and TBA density (d TBA−TOL ) is~4 Å, the distance is greatly reduced to just~1 Å for x TBA = 0.49 and x TBA = 0.24. The difference in the d TBA−TOL as a function of the TBA molar fraction is due to the progressive increase in TOL. Indeed, for the TBA fraction lower than 0.5 corresponding to an unsaturated solid interface in TBA molecules, the SiOH surface can then adsorb the TOL molecules. The progressive adsorption in TOL molecules at the silica surface is then at the origin of the decrease of d TBA−TOL with the increase in TOL fraction and then is at the origin of the interpenetration of TOL molecules. That can be observed from the atomic density in We report in Fig. 6a, b the two-dimensional densities according to x and y directions for x TBA = 0.71, for TBA and TOL. As shown in Fig. 6, the demixing follows the cylindrical symmetry and a core-shell organization is highlighted with the TBA close to the interface (shell) and TOL at the center of the nanopore (core) that is in good agreement with the experiment. 15 To ensure the microphase separation from a thermodynamics point of view we computed the local surface tension (γ(r)) 42 along the radial distance. Local surface tension expresses the surface tension through the radial r-dependent local components of the pressure tensor. Local surface tension is a local interfacial property related to the difference between the normal (P n ) and tangential (P t ) pressure components such that γ(r) = P n (r) − P t (r) with  CH3 is the carbon of the methyl groups of toluene, CH the carbon of the CH groups of toluene, and C the carbon without hydrogen atoms in toluene The mechanical equilibrium involves that the normal pressure has to be constant even through the interface, whereas the tangential component is similar than the normal one in bulk phase and presents a negative peak at the interface highlighting a region under tension. Therefore, the difference P n (r) − P t (r) is null in bulk phase and presents a positive peak at the interface. Even if its definition is arbitrary (how the energy between a pairwise interaction is distributed between both particles), its calculation allows us to highlight the presence of an interface. Calculation of γ(r) is detailed in the Methods section. As shown in Fig. 6c, the local surface tension presents two peaks. Actually, a peak in surface tension exhibits an interface under tension and corresponds to an interface between two phases. From 0 to 8 Å no tension is observed and the local surface tension is close to 0 mN m −1 as expected for a homogeneous liquid phase. Figure 6c emphasizes then that two interfaces identified by the first peak in γ(r) located at 10.5 Å corresponds to the tension radius of surface of tension connected to the TBA-TOL liquid-liquid interface, and the another at 11.5 Å related to the liquid-solid interface. Figure 7a reports the profiles of the surface tension and the density of TBA, TOL, and the atoms of silica material. As shown in Fig. 7a, both peaks of tension are located close to the liquid-liquid and liquid-solid Gibbs dividing surfaces. The Gibbs model for the boundary between two phases consists of considering the finite interfacial region as a dividing surface which acts as a third phase of zero volume in which some magnitudes of physical properties such as density change abruptly. Local surface tension and density profile calculations then unambiguously highlight the phase separation of two confined TBA and TOL liquids at the nanoscale. The formation of a first TBA layer anchored at the silica surface is induced by the strong hydrogen bonds between the silanol groups SiOH of the silica material and the OH groups of the TBA molecules. Indeed, hydrogen bond numbers were computed by considering a geometrical criterium based on the quantum consideration 43 where two molecules are chosen as being hydrogen bonded only if their inter-oxygen distance is less than 3.5 Å, and simultaneously the oxygen-hydrogen is less than 2.5 Å. The profile of the hydrogen bond number (nHB) per molecule was then reported in Fig. 7b for x TBA = 0.24. As shown in Fig. 7b, a peak of nHB was observed close to the silica surface beyond the bulk value. As observed in Fig. 7b, this peak is the result of two contributions, the HBs between silica wall and the OH groups of TBA and HB between TBA molecules. This highlights a strong hydrogen bonding network between TBA and silica material that is at the origin of the core-shell organization. To ensure this hypothesis, MD simulations through a weakly hydrophilic silica nanopore (WH) was carried out. Figure 8a shows that the weak ability of WH to establish interfacial hydrogen bonds between the silanol groups and the TBA molecules explains the absence of microphase separation. This bears out the fact that the interfacial anchoring from HB is the driving force to observe a microphase segregation. Interestingly, Fig. 8b exhibits that the carbon atoms of TBA and TOL are preferentially adsorbed in comparison with the hydroxide groups of TBA that underlines a hydrophobic anchoring related to the interactions between silicon and carbon atoms. It seems then that the strong interfacial hydrogen bonding network between TBA molecules and the silica material rules the preferential adsorption of TBA at the interface. To strengthen this result, MD simulation of confined TBA/TOL mixture into a pure hydrophilic nanopore for a TBA concentration of x TBA = 0.24 was carried out. The ideal hydrophilic nanopore corresponds to a water nanotube (WNT) and details on pore building and force field are provided in the Methods section. As shown in Fig. 7c, TBA molecules are preferentially adsorbed at the silica surface given the strong hydrogen bonds between TBA and water molecules leading to a core-shell structure with a homogeneous TBA/TOL mixture at the center of pore. These results indicate that the highly hydrophilic surface induces a microphase separation of TBA/TOL mixture and could also suggest that a hydrophobic surface is not able to discriminate TBA and TOL molecules because TOL is fully hydrophobic. However, Fig. 7d shows that a carbon nanotube  (CNT) is also capable of separating TBA and TOL by forming a core-shell structure such that TOL and TBA molecules are respectively preferentially located at the CNT surface and at the center of pore. As shown in Fig. 9, this separation is the result of a commensurate organization due to stacking interactions between benzenic cycles of CNT and TOL molecules. To go a little further, the curvature effect was also investigated. Indeed, as often shown, the curvature induced by the cylindrical geometry can drastically impact the physical properties in relation to the slitlike confinement. 17,18 We then carried out MD simulation of a TBA/TOL mixture with a TBA concentration of x TBA = 0.5 confined through two silica slitlike pores such that the confinement in term of confined volume was conserved (two silica walls separated by a distance of 24 Å as illustrated in Fig. 10a). As observed in Fig. 10b, segregation between both liquids is observed highlighting that the separation is rather ruled by the interactions between the silica surface and TBA than the pore geometry. Eventually, pore size effect was also investigated and MD simulations of confined TBA/TOL mixture with a TBA concentration of x TBA = 0.5 confined in a silica nanopore of radius of 6 Å and 18 Å were carried out. As shown in Fig. 11, the microphase separation and the core-shell organization are still observed that corroborates that the driving force is connected to the interfacial anchoring by hydrogen bonds allowing that preferential adsorption of TBA molecules. Additionally, Morineau and colleagues 38 have recently shown, from Neutron scattering, that the core-shell organization was also observed with a pore radius of 43 Å. Therefore, the solid-liquid interactions and more especially the solid-liquid hydrogen bonds seem to overcome the bulk ones whatever the pore diameter leading to the microphase separation and the core-shell structure. Eventually, by comparing the structures obtained from the silica nanotube and from the CNT material (smoothness nanopore), we show that the core-shell organization is independent of the surface roughness. Although the core-shell structure is indirectly correlated to the capillary condensation, the main mechanism is related to the wetting on the silica surface and then to the interfacial interactions between TBA/TOL and the SiOH. Therefore, the microphase separation at the nanoscale cannot be related to any putative de-wetting phenomenon, nor it could be understood by considering solely the capillary condensation of the gas. Indeed, a vision of the intermolecular interactions and then the strength of wetting is needed to capture the origin of the microphase, what was performed in this work.

DISCUSSION
In this work, the TBA/TOL mixture that is homogeneous in the bulk liquid phase was confined at the nanoscale into cylindrical pores. We showed that the nano-confinement can induce, under a few conditions, microphase separation of binary mixture with a core-shell organization in line with Neutron scattering experiments. Phase separation was evidenced by highlighting both liquid-liquid and liquid-solid interfaces by means of the profiles of the density of center of mass and the local surface tension. By  Fig. 7 a Local surface tension (γ(r)) as a function of the radial distance (left axis) and radial densities of center of mass of tert-butanol (TBA) and toluene (TOL) molecules and radial density of atoms belonging to the silica material (right axis). b Profiles of the hydrogen bonds number per TBA molecules along the radial direction of the silica cylindrical nanopore. c Profile of the radial density of TBA and TOL liquids confined into a water nanotube (WNT) with a pore radius of 12 Å. d Profile of the radial density of TBA and TOL liquids confined into a carbon nanotube (CNT) with a pore radius of 12 Å computing the radial profile of the hydrogen bonds number and by tuning the hydrophilicity of the surface of the silica nanopore, we established that the phase separation was ruled at the molecular scale by the strong hydrogen bonds between TBA molecules and the silica material that favors the HB anchoring. Confinement through a carbon nanotube showed an inverse core-shell structure with TOL anchored at the solid surface highlighting an hydrophobic anchoring due to a commensurate organization implying the benzenic cycles stacking. These results shed light on the possibility to control the nanostructure of multicomponent fluids under confinement at the nanoscale by tuning the surface chemistry of nanopore. The nanoconfined phase separation seems to be independent to the pore size and is essentially connected to the strength of pore-fluid interactions. As shown in Fig. 8a, the microphase separation has not been established from a weakly hydrophilic pore, whereas the core-shell structure was also recovered from a hydrophobic matrix (CNT), highlighting that the phase separation is unambiguously connected to the liquid-solid interactions. Determination of a theoretical model predicting the phase separation as a function of the strength of the solid-liquid interactions could be possible by tuning progressively (thermodynamic integration) the pore-liquid interactions (Lennard-Jones parameters and partial charges). However, it will be difficult to connect this model with the realistic materials where the interactions are too complexes. We think that the so-developed strategy consisting of the refining of the pore-wall interactions by a comparison with experiment is probably the most relevant approach. Qualitatively, the microphase separation can be then predicted if one considers a mixture of an amphiphilic component and a hydrophobic one confined in pure hydrophilic and hydrophobic pores. Although the development of a theoretical model to predict the phase separation in a nanoconfined pore seems to be difficult to connect with realistic materials, it is possible to predict the optimal conditions to observe the core-shell organization that means the maximal number of TBA molecules anchored at the surface. Indeed, the pore volume can be calculated as V pore ¼ πR 2 pore L z where R pore is the pore radius of the nanopore while L z is the height of the nanotube, the volume of the core region is V core ¼ πR 2 core L z such that R core is the radius of the core region. The volume of the interfacial zone, i.e., the shell can be then evaluated as V shell ¼ πL z ðR 2 pore À R 2 core Þ such that R core = R pore − e shell with e shell corresponding to the thickness of the interfacial shell. e shell is computed from the density profile. From the density profile calculations, we found that e shell = 4 Å. If we consider the bulk density as the interfacial density, the number of TBA molecules interfacially anchored can be calculated as N TBA,s = V shell .ρ TBA,bulk . With R pore = 12 Å and e shell = 4 Å, an interfacial saturation of N TBA,s = 57 was predicted. Therefore, for the TBA molar fraction corresponding to a number of TBA molecules less than N TBA,s , the TBA will be fully anchored and a full core-shell organization will be observed (i.e., x TBA = 0.71, 0.49, and 0.24). Beyond N TBA,s , TBA molecules will also be located in the center of pore and in the intermediary layer. From GCMC simulations, we found that the solid interface is saturated for N TBA,s = 56 molecules that bears out our previous structural model. In Fig. 4a, the molar fraction x TBA = 0.83 corresponds to 79 TBA and 16 toluene molecules. Therefore, other TBA molecules (79-57) are distributed at the center of pore and in the intermediate layer. As shown in Fig. 1b, the first layer is located between 0 and 3 Å, the second between 3 and 7.5 Å, while the last is located between 7.5 and 12 Å. As the layer size corresponds to the molecule size and that TOL molecule is higher than TBA one, 3.5 Å vs 3.0 Å, it appears that TOL molecules have no place at the center of pore contrary to the intermediate layer where the layer size is 4.5 Å.
Eventually, the microphase separation raises exciting questions about its impact on the MD and the liquid flow in nanochannel, as suggested by the recent report on multiple glassy dynamics. 37 Further insight into this aspect could be attained from ongoing experimental spectroscopy and molecular simulation studies.

Models
Silica cylindrical nanopore with a hydrophilic surface was managed by applying the procedure proposed by Bródka and Zerda. 44 We generated a cylindrical cavity along the z axis of the cubic silica cell of 35.7 Å by removing the atoms within a cylinder of diameter (D) 24 Å. From their coordination numbers, we distinguished bridging oxygen (O b ) bonded to two silicon atoms from non bridging oxygens (Onb) bonded to only one silicon and bonded to one hydrogen atom (Hnb). An iterative procedure of atom (O and Si) removal was applied until only tetra-coordinated silicon atoms, bonded to a maximum of two Onbs, were present in the structure. Finally, non-bridging oxygens were saturated with hydrogen atoms to form surface hydroxyl groups. This procedure leads to a realistic description of the irregular inner surface of the porous silicate and of the interfacial interactions between the fluid and the matrix (see Fig. 2). The inner surface coverage of silanol groups was about 7.5 nm −2 , which corresponds to highly hydrated protonated silica pore. 44,45 Although the  Fig. 8 a Profile of the radial density of center of mass of tert-butanol (TBA) and toluene (TOL) molecules confined into the weakly hydrophilic silica nanopore. b Profile of the radial density of tertiary carbon (Ct), hydrogen atom of hydroxide group (HO) of TBA, and carbon of benzenic cycle (C) and methyl group (CH3) of TOL silica matrix was subsequently kept rigid, rotation around the Si-O bond of the hydroxyl groups was allowed from the SHAKE constraints algorithm, 10,41,46 where the distance between the oxygen and hydrogen atoms are kept fixed at 1.09 Å.
Intermolecular interactions are the sum of both electrostatic and dispersive-repulsive Lennard-Jones contributions. Silica framework was modeled using the ClayFF force field developed by Cygan et al. 40 Indeed, it has been shown that the potential developed by Bródka and Zerda 44 was too weak and underestimated the adsorbed amount at low pressure. 39,47 TBA and TOL were modeled by means of the non-polarizable flexible OPLS all atoms (AA) force field. 48 Indeed, the OPLS-AA models of TBA and TOL were found to well reproduce density and structure in liquid bulk phase. 48 Whereas the original Lennard-Jones parameters were initially conserved, the partial charges were calculated from ab-initio calculations. Calculation of the partial charges was carried out on the basis of 6-311G(d,p) Gaussiantype basis set. In the first time, geometry was optimized, whereas in the second stage the partial charges were calculated from the CHELPG (CHarges from ELectrostatic Potentials using a Grid-based) method. 49 These calculations were performed using the Gaussian code. 50 Force field parameters of silica nanopore, TOL, and TBA molecules are provided in the FIELD.txt file, the input file of DL_POLY software. 51 Crossed van der Waals interactions between TOL and TBA were calculated using the Lorentz-Berthelot mixing rules. The LJ parameters between the silica nanopore and TOL were optimized to qualitatively reproduce the experimental isotherm by following the procedure developed in ref. 41 . The so-optimized interactions between TOL molecules and silica material are provided in Table 1. Force field parameters of silica material, TOL, and TBA and crossed interactions can be found in the FIELD.txt file. Armchair  Fig. 11 Profiles of radial density of the center of mass of tert-butanol (TBA) and toluene (TOL) molecules for x TBA = 0.50 for three pore radii (R = 6, 12, and 18 Å) CNT of radius 12 Å with a pore length of 100 Å was modeled by considering the uncharged force field developed by Werder et al. 52 WNT was built by carving a cylindrical nanopore of radius 12 Å into an equilibrated cubic water box with a length box of 59 Å. Water molecules were modeled by considering the TIP4P/2005 model 53 and were considered as frozen. Weakly hydrophilic silica membrane was also modeled by combining the Universal force field (UFF) 54 for describing oxygen and hydrogen atoms, while silicon atoms were described using the DREIDING force field. 55 All details are provided in ref. 39 . Crossed Lennard-Jones parameters were calculated by considered the Lorentz-Berthelot mixing rules.
MD simulation MD simulations were performed using a time step of 0.002 ps to sample 10 ns (acquisition phase). The equilibration time corresponds to 10 ns. All MD simulations have been carried out with the DL_POLY package 51 using the combination of the velocity-Verlet algorithm and the Nose-Hoover thermostat. 56 Numbers of confined molecules of pure TBA and TOL and binary TBA/TOL were calculated using GCMC simulation at the saturation pressure vapor and 308 K. Compositions in TBA and TOL are provided in Table 2.

Monte Carlo simulation
Monte Carlo simulations were performed in the grand canonical statistical ensemble. Periodic boundary conditions were applied in the three directions. Silica nanopore was considered as rigid and only the hydrogen of silanol can rotate from an angular move implying SiOH groups. 44 Each cycle consisted of N randomly selected moves with fixed probabilities: translation of the center of mass of a randomly chosen, rotation of a randomly selected around its center of mass, change of the internal conformation using the configurational bias regrowth move 57 and insertion/deletion trial move. The frequencies of each type of move are 0.20 for translation, 0.20 for rotation, 0.20 for the change of the conformation, and 0.4 for insertion/deletion. GCMC simulations consisted of 700,0000 cycles.

Surface tension calculation
Surface tension was computed using the non-exponential test-area method TA2 58 based upon a thermodynamic route. 59 This method expresses the surface tension as a change in the free energy (F) for an infinitesimal change in the surface area performed in the constant-NVT ensemble. This infinitesimal change in the area is performed throughout a perturbation process for which the perturbed system (state A + ΔA and noted 1) is obtained from an infinitesimal change ΔA of the area A of the reference system (noted 0). The box dimensions ðL z =ð1 ± ξÞ where ξ → 0. The area of a cylindrical interface is A = 2πR e L z and ΔA = A [(1 ± ξ) −1/2 − 1], where R e is the radius equimolar dividing surface. These transformations conserve the volume of the box in the perturbed state. This means that it is possible to express the surface tension as a difference between the perturbed and reference states. U (0) (r N ) and U (1) (r′ N ) are the configurational energies of the systems with an area A and a configurational space r N , and an area A + ΔA and a configurational space r′ N , respectively. (1) where <···> 0 indicates that the average is carried out over the reference state. A local version of Eq. (1) can be obtained by assuming a decorrelation of the cylindrical slabs. 60 Rk ðr ij Þ À u Rk ðr ij ÞÞ ΔA ! * + 0 ; ( where k is the index of the cylindrical slab, R k the radius of the cylindrical shell, u Rk is the energy of the kth element, H(R ik ) is the Heaviside function with H(R ik ) = 1 for R i = R k and 0 otherwise and r ij is the distance between i and j molecules.

DATA AVAILABILITY
The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request. These numbers were taken from the grand canonical Monte Carlo (GCMC) simulations at the saturation pressure vapor and 308 K