On the coupling between ionic conduction and dipolar relaxation in deep eutectic solvents: Influence of hydration and glassy dynamics.

We have studied the ionic conductivity and the dipolar reorientational dynamics of aqueous solutions of a prototypical deep eutectic solvent (DES), ethaline, by dielectric spectroscopy in a broad range of frequencies (MHz-Hz) and for temperatures ranging from 128 to 283 K. The fraction of water in the DES was varied systematically to cover different regimes, starting from the pure DES and its water-in-DES mixtures to the diluted electrolyte solutions. Depending on these parameters, different physical states were examined, including low viscosity liquid, supercooled viscous liquid, amorphous solid, and freeze-concentrated solution. Both the ionic conductivity and the reorientational relaxation exhibited characteristic features of glassy dynamics that could be quantified from the deviation from the Arrhenius temperature dependence and non-exponential decay of the relaxation function. A transition occurred between the water-in-DES regime (<40 wt. %), where the dipolar relaxation and ionic conductivity remained inversely proportional to each other, and the DES-in-water regime (>40 wt. %), where a clear rotation-translation decoupling was observed. This suggests that for a low water content, on the timescale covered by this study (∼10-6 to 1 s), the rotational and transport properties of ethaline aqueous solutions obey classical hydrodynamic scaling despite these systems being presumably spatially microheterogeneous. A fractional scaling is observed in the DES-in-water regime due to the formation of a maximally freeze-concentrated DES aqueous solution coexisting with frozen water domains at sub-ambient temperature.

Dynamics of water confined in mesopores with variable surface interaction.

We have investigated the dynamics of liquid water confined in mesostructured porous silica (MCM-41) and periodic mesoporous organosilicas (PMOs) by incoherent quasielastic neutron scattering experiments. The effect of tuning the water/surface interaction from hydrophilic to more hydrophobic on the water mobility, while keeping the pore size in the range 3.5 nm-4.1 nm, was assessed from the comparative study of three PMOs comprising different organic bridging units and the purely siliceous MCM-41 case. An extended dynamical range was achieved by combining time-of-flight (IN5B) and backscattering (IN16B) quasielastic neutron spectrometers providing complementary energy resolutions. Liquid water was studied at regularly spaced temperatures ranging from 300 K to 243 K. In all systems, the molecular dynamics could be described consistently by the combination of two independent motions resulting from fast local motion around the average molecule position and the confined translational jump diffusion of its center of mass. All the molecules performed local relaxations, whereas the translational motion of a fraction of molecules was frozen on the experimental timescale. This study provides a comprehensive microscopic view on the dynamics of liquid water confined in mesopores, with distinct surface chemistries, in terms of non-mobile/mobile fraction, self-diffusion coefficient, residence time, confining radius, local relaxation time, and their temperature dependence. Importantly, it demonstrates that the strength of the water/surface interaction determines the long-time tail of the dynamics, which we attributed to the translational diffusion of interfacial molecules, while the water dynamics in the pore center is barely affected by the interface hydrophilicity.

More room for microphase separation: An extended study on binary liquids confined in SBA-15 cylindrical pores.

The confinement of liquid mixtures in porous channels provides new insight into fluid ordering at the nanoscale. In this study, we address a phenomenon of microphase separation, which appears as a novel fascinating confinement effect for fully miscible binary liquids. We investigate the structure of tert-butanol-toluene mixtures confined in the straight and mono-dispersed cylindrical nanochannels of SBA-15 mesoporous silicates (D = 8.3 nm). Small angle neutron scattering experiments on samples with carefully designed isotopic compositions are performed to systematically vary the scattering length density of the different compounds and assess the radial concentration profile of the confined phases. The resulting modulation of the Bragg reflections of SBA-15 is compared with the predictions from different core-shell models, highlighting a molecular-scale phase-separated tubular structure with the tert-butanol forming a layer at the pore surface, surrounding a toluene-rich core. The present structural study suggests that the microphase separation phenomenon in confinement, which so far had only been reported for a smaller pore size (D = 3.65 nm) and a unique mixture composition, must be considered as a general phenomenon. It also highlights the strength of neutron scattering method with isotopic substitution, which is a unique experimental approach to reveal this phenomenon.

Thermodynamics of binary gas adsorption in nanopores.

MCM-41 nanoporous silicas show a very high selectivity for monoalcohols over aprotic molecules during adsorption of a binary mixture in the gas phase. We present here an original use of gravimetric vapour sorption isotherms to characterize the role played by the alcohol hydrogen-bonding network in the adsorption process. Beyond simple selectivity, vapour sorption isotherms measured for various compositions help to completely unravel at the molecular level the step by step adsorption mechanism of the binary system in the nanoporous solid, from the first monolayers to the complete liquid condensation.

Thermotropic orientational order of discotic liquid crystals in nanochannels: an optical polarimetry study and a Landau-de Gennes analysis.

Optical polarimetry measurements of the orientational order of a discotic liquid crystal based on a pyrene derivative confined in parallelly aligned nanochannels of monolithic, mesoporous alumina, silica, and silicon as a function of temperature, channel radius (3-22 nm) and surface chemistry reveal a competition of radial and axial columnar orders. The evolution of the orientational order parameter of the confined systems is continuous, in contrast to the discontinuous transition in the bulk. For channel radii larger than 10 nm we suggest several, alternative defect structures, which are compatible both with the optical experiments on the collective molecular orientation presented here and with a translational, radial columnar order reported in previous diffraction studies. For smaller channel radii our observations can semi-quantitatively be described by a Landau-de Gennes model with a nematic shell of radially ordered columns (affected by elastic splay deformations) that coexists with an orientationally disordered, isotropic core. For these structures, the cylindrical phase boundaries are predicted to move from the channel walls to the channel centres upon cooling, and vice-versa upon heating, in accord with the pronounced cooling/heating hystereses observed and the scaling behavior of the transition temperatures with the channel diameter. The absence of experimental hints of a paranematic state is consistent with a biquadratic coupling of the splay deformations to the order parameter.

Discotic columnar liquid crystal studied in the bulk and nanoconfined states by molecular dynamics simulation.

A prototypical Gay Berne discotic liquid crystal was studied by means of molecular dynamics simulations both in the bulk state and under confinement in a nanoporous channel. The phase behavior of the confined system strongly differs from its bulk counterpart: the bulk isotropic-to-columnar transition is replaced by a continuous ordering from a paranematic to a columnar phase. Moreover, a new transition is observed at a lower temperature in the confined state, which corresponds to a reorganization of the intercolumnar order. It reflects the competing effects of pore surface interaction and genuine hexagonal packing of the columns. The translational molecular dynamics in the different phases has been thoroughly studied and discussed in terms of collective relaxation modes, non-Gaussian behavior, and hopping processes.

Crossover in structure and dynamics of a primary alcohol induced by hydrogen-bonds dilution.

Primary alcohols show a prominent Debye process in the dielectric relaxation located at a timescale longer than the main structural relaxation. Böhmer and co-workers studied dilution effects of the hydrogen bonding network of n-butanol (BuOH) with n-bromobutane (BuBr) to better understand the origin of this process. Interestingly, this work has evidenced a crossover in Debye relaxation time (τD) for a critical concentration in BuBr xc = 0.5. By using molecular dynamics simulations and NMR experiments we propose here to explore further dilution effects on the dipolar and translational dynamics. Moreover, we discuss the relation between structural and dynamical properties in the context of a detailed study of the microstructure and the H-bond network. The overall results are consistent with the existence of a topological change in the liquid structure occurring at about xc = 0.5 from a hypernetted percolating network to independent nanodomains of n-butanol molecules embedded in the n-bromobutane phase.

Configurational temperature and local properties of the anisotropic Gay-Berne liquid crystal model: applications to the isotropic liquid/vapor interface and isotropic/nematic transition.

Molecular simulations in the isothermal statistical ensembles require that the macroscopic thermal and mechanical equilibriums are respected and that the local values of these properties are constant at every point in the system. The thermal equilibrium in Monte Carlo simulations can be checked through the calculation of the configurational temperature, k(B)T(conf)=<|∇(r)U(r(N))|(2)>/<∇(r) (2)U(r(N))>, where ∇(r) is the nabla operator of position vector r. As far as we know, T(conf) was never calculated with the anisotropic Gay-Berne potential, whereas the calculation of T(conf) is much more widespread with more common potentials (Lennard Jones, electrostatic, ...). We establish here an operational expression of the macroscopic and local configurational temperatures, and we investigate locally the isotropic liquid phase, the liquid / vapor interface, and the isotropic-nematic transition by Monte Carlo simulations.

Molecular simulations of confined liquids: an alternative to the grand canonical Monte Carlo simulations.

Commonly, the confinement effects are studied from the grand canonical Monte Carlo (GCMC) simulations from the computation of the density of liquid in the confined phase. The GCMC modeling and chemical potential (µ) calculations are based on the insertion/deletion of the real and ghost particle, respectively. At high density, i.e., at high pressure or low temperature, the insertions fail from the Widom insertions while the performing methods as expanded method or perturbation approach are not efficient to treat the large and complex molecules. To overcome this problem we use a simple and efficient method to compute the liquid's density in the confined medium. This method does not require the precalculation of µ and is an alternative to the GCMC simulations. From the isothermal-isosurface-isobaric statistical ensemble we consider the explicit framework/liquid external interface to model an explicit liquid's reservoir. In this procedure only the liquid molecules undergo the volume changes while the volume of the framework is kept constant. Therefore, this method is described in the Np(n)AV(f)T statistical ensemble, where N is the number of particles, p(n) is the normal pressure, V(f) is the volume of framework, A is the surface of the solid/fluid interface, and T is the temperature. This approach is applied and validated from the computation of the density of the methanol and water confined in the mesoporous cylindrical silica nanopores and the MIL-53(Cr) metal organic framework type, respectively.

Extension and Limits of Cryoscopy for Nanoconfined Solutions.

This work investigates the phase behavior of aqueous solutions of glycerol confined in MCM-41 and SBA-15 nanoporous matrixes by calorimetry. Limitations due to overfilling and eutectic freezing are prevented by the absence of an external liquid reservoir and by the glass-forming property of glycerol. Consequently, the stability of nanoconfined ice in equilibrium with aqueous solutions is studied over a wide range of compositions. In confinement, a large temperature depression of the liquidus line is observed. A thermodynamic model accounting simultaneously for the cryoscopic and the Gibbs-Thomson effects gives a consistent view of the phase diagram for large pores (Rp = 4.15 nm). For smaller pores (Rp = 1.8 nm), it reveals that the water activity strongly deviates from the bulk solution with the same composition, indicating the possible role of concentration heterogeneities in determining the onset of ice freezing in strongly nanoconfined solutions.

Do Deep Eutectic Solvents Form Uniform Mixtures Beyond Molecular Microheterogeneities?

We have performed small-angle neutron scattering in a momentum transfer range (0.05 < Q < 0.5 Å-1) to study long-range order and concentration fluctuations in deep eutectic solvents (DESs) and their aqueous solutions. Ethaline (choline chloride/ethylene glycol), glycerol/lactic acid, and menthol/decanoic acid mixtures were selected to illustrate individually the case of ionic, nonionic, and hydrophobic mixtures. Carefully designed isotopic labeling was used to emphasize selectively the spatial correlations between the different solvent components. For ethaline DESs and their aqueous solutions, a weak low-Q peak observed only for certain compositions and some partial structure factors revealed the mesoscopic segregation of ethylene glycol molecules that do not participate in the solvation of ionic units, either because they are in excess with respect to the eutectic stoichiometry (1:4 neat ethaline) or substituted by water (4w-ethaline and higher aqueous dilutions). For the nonionic hydrophilic solutions, such a mesoscopic segregation was not observed. This indicates that the better balanced interactions between the three nonionic H-bonded components (water, lactic acid, and glycerol) favor homogeneous mixing. For the hydrophobic DESs, we observed an excess of coherent scattering intensity centered at Q = 0, which could be reproduced by a model of noninteracting spherical domains. Local concentration fluctuations are not excluded either. However, unlike liquid mixtures with a tendency to demix, we have found no evidence of expansion of domains with different compositions to a large scale.

Dynamic Heterogeneities in Liquid Mixtures Confined in Nanopores.

Binary liquid mixtures can exhibit nanosegregation, albeit being fully miscible and homogeneous at the macroscopic scale. This tendency can be amplified by geometrical nanoconfinement, leading to remarkable properties. This work investigates the molecular dynamics of tert-butanol (TBA)-toluene (TOL) mixtures confined in silica nanochannels by quasielastic neutron scattering and molecular dynamics simulation. It reveals a decoupling of the molecular motion of each constituent of the binary liquid, which can be followed independently by selective isotopic H/D labeling. We argue that this behavior is the signature of spatially segregated dynamic heterogeneities, which are due to the recently established core-shell nanophase separation induced by mesoporous confinement.

Structure and dynamics of a Gay-Berne liquid crystal confined in cylindrical nanopores.

Gay-Berne liquid crystals confined in two cylindrical nanopores with different pore sizes were studied by molecular dynamics simulation. Their structure and dynamics properties were obtained and compared with those of the bulk. Our data show that confinement changes the bulk isotropic-to-nematic transition to a continuous ordering from a paranematic to a nematic phase. Moreover, confinement strongly hinders the smectic translational order. The molecular dynamics is characterized by the translational diffusion coefficients and the first-rank reorientational correlation times. Very different characteristic times and temperature variations in the dynamics are observed in confinement. Spatially resolved quantities illustrate that confinement induces predominant structural and dynamical heterogeneities.

Characterization of molecular motions in biomolecular systems by elastic incoherent neutron scattering.

In the present work the role played by the instrumental resolution function in elastic incoherent neutron scattering (EINS) experiment is discussed. An important result consists in the definition of an equivalent time t(*), which depends both on the characteristic system time and on the resolution time, for which the spatial Fourier transform of EINS intensity profile and the self-distribution function (SDF) evaluated at t=t(*) are proportional. Then the equivalent time t(*) is introduced in the SDF procedure, an operational recipe for the mean square displacement determination. The new revised procedure is applied on data of myoglobin in trehalose dry environment and of hydrated homologous disaccharides (sucrose and trehalose).

Relation between static short-range order and dynamic heterogeneities in a nanoconfined liquid crystal.

We analyze the molecular dynamics heterogeneity of the liquid crystal 4-n-octyl-4'-cyanobiphenyl nanoconfined in porous silicon. We show that the temperature dependence of the dynamic correlation length xi_(wall) , which measures the distance over which a memory of the interfacial slowing down of the molecular dynamics persists, is closely related to the growth of the short-range static order arising from quenched random fields. More generally, this result may also shed some light on the connection between static and dynamic heterogeneities in a wide class of condensed and soft matter systems.

Phase diagram and glass transition of confined benzene.

We used differential scanning calorimetry, neutron scattering, and proton NMR to investigate the phase behavior, the structure, and the dynamics of benzene confined in a series of cylindrical mesoporous materials MCM-41 and SBA-15 with pore diameters, d, between 2.4 and 14 nm. With this multitechnique approach, it was possible to determine the structure and, for the first time to our knowledge, the density of confined benzene as a function of temperature and pore size. Under standard cooling rates, benzene partially crystallizes in SBA-15 matrixes (4.7

Evidence of anisotropic quenched disorder effects on a smectic liquid crystal confined in porous silicon.

We present a neutron scattering analysis of the structure of the smectic liquid crystal octylcyanobiphenyl (8CB) confined in one-dimensional nanopores of porous silicon films (PS). The smectic transition is completely suppressed, leading to the extension of a short-range ordered smectic phase aligned along the pore axis. It evolves reversibly over an extended temperature range, down to 50 K below the N-SmA transition in pure 8CB. This behavior strongly differs from previous observations of smectics in different one-dimensional porous materials. A coherent picture of this striking behavior requires that quenched disorder effects are invoked. The strongly disordered nature of the inner surface of PS acts as random fields coupling to the smectic order. The one-dimensionality of PS nanochannels offers perspectives on quenched disorder effects, of which observation has been restricted to homogeneous random porous materials so far.

High-resolution dielectric study reveals pore-size-dependent orientational order of a discotic liquid crystal confined in tubular nanopores.

We report a high-resolution dielectric study on a pyrene-based discotic liquid crystal (DLC) in the bulk state and confined in parallel tubular nanopores of monolithic silica and alumina membranes. The positive dielectric anisotropy of the DLC molecule at low frequencies (in the quasistatic case) allows us to explore the thermotropic collective orientational order. A face-on arrangement of the molecular discs on the pore walls and a corresponding radial arrangement of the molecules is found. In contrast to the bulk, the isotropic-to-columnar transition of the confined DLC is continuous, shifts with decreasing pore diameter to lower temperatures, and exhibits a pronounced hysteresis between cooling and heating. These findings corroborate conclusions from previous neutron and x-ray-scattering experiments as well as optical birefringence measurements. Our study also indicates that the relative simple dielectric technique presented here is a quite efficient method in order to study the thermotropic orientational order of DLC-based nanocomposites.

Collective molecular reorientation of a calamitic liquid crystal (12CB) confined in alumina nanochannels.

We study the smectic director structure of the rodlike liquid crystal 4-n-dodecyl-4'-cyanobiphenyl (12CB) confined in cylindrical cavities of 200 nm diameter in porous alumina templates by means of combined broadband dielectric spectroscopy, optical birefringence, and neutron scattering measurements. We show that the collective molecular orientation differs between entering the smectic A phase upon cooling from the isotropic state and entering the same phase upon heating while melting the confined crystal. We discuss this collective molecular realignment in terms of a competition between weak planar anchoring at the p-Al2O3/12CB interface and a preferred texture typical of the crystallization of rodlike molecules in nanochannels (Bridgman growth).

Criticality of an isotropic-to-smectic transition induced by anisotropic quenched disorder.

We report combined optical birefringence and neutron scattering measurements on the liquid crystal 12CB nanoconfined in mesoporous silicon layers. This liquid crystal exhibits strong nematic-smectic coupling responsible for a discontinuous isotropic-to-smectic phase transition in the bulk state. Confined in porous silicon, 12CB is subjected to strong anisotropic quenched disorder: a short-ranged smectic state evolves out of a paranematic phase. This transformation appears continuous, losing its bulk first-order character. This contrasts with previously reported observations on liquid crystals under isotropic quenched disorder. In the low temperature phase, both orientational and translational order parameters obey the same power law.