Anisotropic relaxation of nuclear spins dipolar energy of water molecules in two-dimensional nanopores - A single crystal NMR study.

Energy transfer from Zeeman to dipolar order discovered by Jeener et al. is usually observed in solids with a strong dipole-dipole interaction of nuclear spins. It is not observed in liquids, where fast molecular motion completely averages this interaction. The intermediate case, when the dipole-dipole interaction of nuclear spins is only partially averaged, has been poorly studied. We report on the first measurement of an angular-dependent proton spin relaxation of a dipolar reservoir in mobile water molecules confined in the interlayer pores of a vermiculite single crystal. In this layered crystal, the intramolecular dipole-dipole interactions of nuclear spins are only partially averaged due to the restricted anisotropic molecular motion in nanopores. We show that this allows the formation of dipolar echo. We measured the spin-lattice relaxation times of the dipolar order T1D at different angles between the normal to the crystal surface and the applied magnetic field and obtained a distinct angular dependence of T1D. The minimum relaxation rate R1D was found around the magic angle of 54.74°.

New insights into the protein stabilizing effects of trehalose by comparing with sucrose.

Disaccharides are well known to be efficient stabilizers of proteins, for example in the case of lyophilization or cryopreservation. However, although all disaccharides seem to exhibit bioprotective and stabilizing properties, it is clear that trehalose is generally superior compared to other disaccharides. The aim of this study was to understand this by comparing how the structural and dynamical properties of aqueous trehalose and sucrose solutions influence the protein myoglobin (Mb). The structural studies were based on neutron and X-ray diffraction in combination with empirical potential structure refinement (EPSR) modeling, whereas the dynamical studies were based on quasielastic neutron scattering (QENS) and molecular dynamics (MD) simulations. The results show that the overall differences in the structure and dynamics of the two systems are small, but nevertheless there are some important differences which may explain the superior stabilizing effects of trehalose. It was found that in both systems the protein is preferentially hydrated by water, but that this effect is more pronounced for trehalose, i.e. trehalose forms less hydrogen bonds to the protein surface than sucrose. Furthermore, the rotational motion around dihedrals between the two glucose rings of trehalose is slower than in the case of the dihedrals between the glucose and fructose rings of sucrose. This leads to a less perturbed protein structure in the case of trehalose. The observations indicate that an aqueous environment closest to the protein molecules is beneficial for an efficient bioprotective solution.

Atomistic molecular dynamics simulations of tubulin heterodimers explain the motion of a microtubule.

Microtubules are essential parts of the cytoskeleton that are built by polymerization of tubulin heterodimers into a hollow tube. Regardless that their structures and functions have been comprehensively investigated in a modern soft matter, it is unclear how properties of tubulin heterodimer influence and promote the self-assembly. A detailed knowledge of such structural mechanisms would be helpful in drug design against neurodegenerative diseases, cancer, diabetes etc. In this work atomistic molecular dynamics simulations were used to investigate the fundamental dynamics of tubulin heterodimers in a sheet and a short microtubule utilizing well-equilibrated structures. The breathing motions of the tubulin heterodimers during assembly show that the movement at the lateral interface between heterodimers (wobbling) dominates in the lattice. The simulations of the protofilament curvature agrees well with recently published experimental data, showing curved protofilaments at polymerization of the microtubule plus end. The tubulin heterodimers exposed at the microtubule minus end were less curved and displayed altered interactions at the site of sheet closure around the outmost heterodimers, which may slow heterodimer binding and polymerization, providing a potential explanation for the limited dynamics observed at the minus end.

Size and Refractive Index Determination of Subwavelength Particles and Air Bubbles by Holographic Nanoparticle Tracking Analysis.

Determination of size and refractive index (RI) of dispersed unlabeled subwavelength particles is of growing interest in several fields, including biotechnology, wastewater monitoring, and nanobubble preparations. Conventionally, the size distribution of such samples is determined via the Brownian motion of the particles, but simultaneous determination of their RI remains challenging. This work demonstrates nanoparticle tracking analysis (NTA) in an off-axis digital holographic microscope (DHM) enabling determination of both particle size and RI of individual subwavelength particles from the combined information about size and optical phase shift. The potential of the method to separate particle populations is demonstrated by analyzing a mixture of three types of dielectric particles within a narrow size range, where conventional NTA methods based on Brownian motion alone would fail. Using this approach, the phase shift allowed individual populations of dielectric beads overlapping in either size or RI to be clearly distinguished and quantified with respect to these properties. The method was furthermore applied for analysis of surfactant-stabilized micro- and nanobubbles, with RI lower than that of water. Since bubbles induce a phase shift of opposite sign to that of solid particles, they were easily distinguished from similarly sized solid particles made up of undissolved surfactant. Surprisingly, the dependence of the phase shift on bubble size indicates that only those with 0.15-0.20 µm radius were individual bubbles, whereas larger bubbles were actually clusters of bubbles. This label-free means to quantify multiple parameters of suspended individual submicrometer particles offers a crucial complement to current characterization strategies, suggesting broad applicability for a wide range of nanoparticle systems.

DOPC versus DOPE as a helper lipid for gene-therapies: molecular dynamics simulations with DLin-MC3-DMA.

Ionizable lipids are important compounds of modern therapeutic lipid nano-particles (LNPs). One of the most promising ionizable lipids (or amine lipids) is DLin-MC3-DMA. Depending on their pharmaceutical application these LNPs can also contain various helper lipids, such as phospho- and pegylated lipids, cholesterol and nucleic acids as a cargo. Due to their complex compositions the structures of these therapeutics have not been refined properly. Therefore, the role of each lipid in the pharmacological properties of LNPs has not been determined. In this work an atomistic model for the neutral form of DLin-MC3-DMA was derived and all-atom molecular dynamics (MD) simulations were carried out in order to investigate the effect of the phospholipid headgroup on the possible properties of the shell-membranes of LNPs. Bilayers containing either DOPC or DOPE lipids at two different ratios of DLin-MC3-DMA (5 mol% and 15 mol%) were constructed and simulated at neutral pH 7.4. The results from the analysis of MD trajectories revealed that DOPE lipid headgroups associated strongly with lipid tails and carbonyl oxygens of DLin-MC3-DMA, while for DOPC lipid headgroups no significant associations were observed. Furthermore, the strong associations between DOPE and DLin-MC3-DMA result in the positioning of DLin-MC3-DMA at the surface of the membrane. Such an interplay between the lipids slows down the lateral diffusion of all simulated bilayers, where a more dramatic decrease of the diffusion rate is observed in membranes with DOPE. This can explain the low water penetration of lipid bilayers with phosphatidylethanolamines and, probably, can relate to the bad transfection properties of LNPs with DOPE and DLin-MC3-DMA.

Stabilization of proteins embedded in sugars and water as studied by dielectric spectroscopy.

In many products proteins have become an important component, and the long-term properties of these products are directly dependent on the stability of their proteins. To enhance this stability it has become common to add disaccharides in general, and trehalose in particular. However, the mechanisms by which disaccharides stabilize proteins and other biological materials are still not fully understood, and therefore we have here used broadband dielectric spectroscopy to investigate the stabilizing effect of the disaccharides trehalose and sucrose on myoglobin, with the aim to enhance this understanding in general and to obtain specific insights into why trehalose exhibits extraordinary stabilizing properties. The results show the existence of three or four clearly observed relaxation processes, where the three common relaxations are the local (ß) water relaxation below the glass transition temperature (Tg), the structural α-relaxation of the solvent, observed above Tg, and an even slower protein relaxation due to large-scale conformational protein motions. For the trehalose containing samples with less than 50 wt% myoglobin a fourth relaxation process was observed due to a ß-relaxation of trehalose below Tg. This latter process, which was assigned to intramolecular rotations of the monosaccharide rings in trehalose, could not be detected for high protein concentrations or for the sucrose containing samples. Since sucrose has previously been found to form more intramolecular hydrogen bonds at the present hydration levels, it is likely that this rotation becomes too slow to be observed in the case of sucrose. However, this sugar relaxation has probably less influence on the protein stability below Tg, where the better stabilizing effect of trehalose on proteins can be explained by our observation that trehalose slows down the water relaxation more than sucrose does. Finally, we show that the α-relaxation of the solvent and the large-scale protein motions exhibit similar temperature dependences, which suggests that these protein motions are slaved by the α-relaxation. Furthermore, the α-relaxation of the trehalose solution is slower than for the corresponding sucrose solution, and thereby also the protein motions become slower in the trehalose solution, which explains the more efficient stabilizing effect of trehalose on proteins above Tg.

Water dynamics in the hydration shells of biological and non-biological polymers.

The dynamics of water at supercooled temperatures in aqueous solutions of different types of solutes has been deeply analyzed in the literature. In these previous works and in most of the cases, a single relaxation of water molecules is observed. In this work, we analyze the dynamics of water in solutions for which a dual relaxation of water molecules is experimentally measured. We discuss the criteria for observing these two water relaxations in these specific solutions and their most likely origins. We also discuss how these two water relaxations relate to the relaxation behavior of bulk water and how the slower one is coupled to the solute dynamics and is essential for the dynamics and functional properties of proteins.

Motions of water and solutes-Slaving versus plasticization phenomena.

It is well-accepted that hydration water is crucial for the structure, dynamics, and function of proteins. However, the exact role of water for the motions and functions of proteins is still debated. Experiments have shown that protein and water dynamics are strongly coupled but with water motions occurring on a considerably faster time scale (the so-called slaving behavior). On the other hand, water also reduces the conformational entropy of proteins and thereby acts as a plasticizer of them. In this work, we analyze the dynamics (using broadband dielectric spectroscopy) of some specific non-biological water solutions in a broad concentration range to elucidate the role of water in the dynamics of the solutes. Our results demonstrate that at low water concentrations (less than 5 wt. %), the plasticization phenomenon prevails for all the materials analyzed. However, at higher water concentrations, two different scenarios can be observed: the slaving phenomenon or plasticization, depending on the solute analyzed. These results generalize the slaving phenomenon to some, but not all, non-biological solutions and allow us to analyze the key factors for observing the slaving behavior in protein solutions as well as to reshaping the slaving concept.

Stable Air Nanobubbles in Water: the Importance of Organic Contaminants.

Nanobubbles, surprisingly stable submicrometer gas bubbles in water, appear to be common in water and biological fluids and are of great interest in technical applications ranging from ultrasound contrast agents to flotation in the mining industry. Nanobubbles on surfaces have been more researched than freely floating bulk nanobubbles, and the reason for their stability appears to be better explained. The stability of bulk nanobubbles is less well explained, several theories exist, and even their existence is sometimes questioned. In the present study, an attempt was made to generate nanobubbles through hydrodynamic cavitation as well as through vigorous shaking in test tubes, and it was found that none of these methods generated a detectable concentration of possible bulk nanobubbles if pure water was used, with or without a small addition of NaCl, the equipment was cleaned properly, and certain plastic materials were avoided. These results indicate that trace organic contaminants are necessary for nanobubble stabilization. Experiments were also made with the dissolution of a high concentration of inorganic salts, which generated bubbles by creating air supersaturation. Light scattering submicron particles were found in all solutions and appeared to be actual gas bubbles in at least one case. However, in many cases, these light scattering particles were unaffected by vacuum and pressure and appear, therefore, to be something else other than air bubbles. It is concluded that, in future research on nanobubble stability, it is very important to avoid contamination, as well as to ascertain that light scattering objects really are bubbles and not oil droplets or particles.

Confined Water as Model of Supercooled Water.

Water in confined geometries has obvious relevance in biology, geology, and other areas where the material properties are strongly dependent on the amount and behavior of water in these types of materials. Another reason to restrict the size of water domains by different types of geometrical confinements has been the possibility to study the structural and dynamical behavior of water in the deeply supercooled regime (e.g., 150-230 K at ambient pressure), where bulk water immediately crystallizes to ice. In this paper we give a short review of studies with this particular goal. However, from these studies it is also clear that the interpretations of the experimental data are far from evident. Therefore, we present three main interpretations to explain the experimental data, and we discuss their advantages and disadvantages. Unfortunately, none of the proposed scenarios is able to predict all the observations for supercooled and glassy bulk water, indicating that either the structural and dynamical alterations of confined water are too severe to make predictions for bulk water or the differences in how the studied water has been prepared (applied cooling rate, resulting density of the water, etc.) are too large for direct and quantitative comparisons.

Possible relations between supercooled and glassy confined water and amorphous bulk ice.

In this paper we discuss apparent contradictions in the literature between dynamical results on supercooled confined water obtained by different experimental methods. The reason for the lack of a clear glass transition of confined water is also discussed. Dielectric relaxation data and results from differential scanning calorimetry measurements provide a consistent picture, but it is still unclear why the glass transition related structural (α) relaxation disappears before the normal time-scale of a calorimetric glass transition (i.e. about 100 s) is reached. From recent results on amorphous bulk ice we propose that this anomalous phenomenon may not be an effect of confinement, but an intrinsic property of water when it transforms to a crystal-like glassy state, probably around 225 K. Thus, the results from the studies of confined water in the so-called no man's land (the temperature range 150-235 K) where bulk water rapidly crystallizes may be of more relevance for supercooled and glassy bulk water than previously thought. Furthermore, the structural difference between glassy water (or amorphous ice) and crystalline ice is likely to be rather small, due to the large degree of disorder in crystalline ice.

Dynamics of aqueous binary glass-formers confined in MCM-41.

Dielectric permittivity measurements were performed on water solutions of propylene glycol (PG) and propylene glycol monomethyl ether (PGME) confined in 21 Å pores of the silica matrix MCM-41 C10 in wide frequency (10(-2)-10(6) Hz) and temperature (130-250 K) ranges. The aim was to elucidate how the formation of large hydrogen bonded structural entities, found in bulk solutions of PGME, was affected by the confined geometry, and to make comparisons with the dynamic behavior of the PG-water system. For all solutions the measurements revealed four almost concentration independent relaxation processes. The intensity of the fastest process is low compared to the other relaxation processes and might be caused by both hydroxyl groups of the pore surfaces and by local motions of water and solute molecules. The second fastest process contains contributions from both the main water relaxation as well as the intrinsic ß-relaxation of the solute molecules. The third fastest process is the viscosity related α-relaxation. Its concentration independency is very different compared to the findings for the corresponding bulk systems, particularly for the PGME-water system. The experimental data suggests that the surface interactions induce a micro-phase separation of the two liquids, resulting in a full molecular layer of water molecules coordinating to the hydrophilic hydroxyl groups on the surfaces of the silica pores. This, in turn, increases the geometrical confinement effect for the remaining solution even more and prevents the building up of the same type of larger structural entities in the PGME-water system as in the corresponding bulk solutions. The slowest process is mainly hidden in the high conductivity contribution at low frequencies, but its temperature dependence can be extracted for the PGME-water system. However, its origin is not fully clear, as will be discussed.

Dynamics of supercooled water in a biological model system of the amino acid L-lysine.

The dynamics of supercooled water in aqueous solutions of the single amino acid L-lysine has been studied by broadband dielectric spectroscopy. The chosen biological system is unique in the sense that the water content is high enough to fully dissolve the amino acid, but low enough to avoid crystallisation to ice at any temperature. This is not possible to achieve for proteins or other larger biomolecules, where either hydrated samples without ice or solutions with large quantities of ice, or a cryoprotectant sugar, have to be studied at low temperatures. Thus, it is a key finding to be able to study water and biomolecular dynamics in a non-crystallized and biologically realistic solution at supercooled temperatures. Here, we focus on the water dynamics in this unique biological solution of L-lysine and water. We show that this unique system also gives rise to unique water dynamics, since, for the first time, a continuation of a cooperative (α-like) water relaxation is observed after a crossover to a more local ß-like water relaxation has occurred with decreasing temperature. This implies that the supercooled water in the biological solution shows a twofold relaxation behaviour, with one relaxation identical to the main relaxation of water in hard confinements and one relaxation almost identical to the main water relaxation in ordinary aqueous solutions.

Anomalous dynamics of aqueous solutions of di-propylene glycol methylether confined in MCM-41 by quasielastic neutron scattering.

The molecular dynamics of solutions of di-propylene glycol methylether (2PGME) and H2O (or D2O) confined in 28 Å pores of MCM-41 have been studied by quasielastic neutron scattering and differential scanning calorimetry over the concentration range 0-90 wt.% water. This system is of particular interest due to its pronounced non-monotonic concentration dependent dynamics of 2PGME in the corresponding bulk system, showing the important role of hydrogen bonding for the dynamics. In this study we have elucidated how this non-monotonic concentration dependence is affected by the confined geometry. The results show that this behaviour is maintained in the confinement, but the slowest diffusive dynamics of 2PGME is now observed at a considerably higher water concentration; at 75 wt.% water in MCM-41 compared to 30 wt.% water in the corresponding bulk system. This difference can be explained by an improper mixing of the two confined liquids. The results suggest that water up to a concentration of about 20 wt.% is used to hydrate the hydrophilic hydroxyl surface groups of the silica pores, and that it is only at higher water contents the water becomes partly mixed with 2PGME. Hence, due to this partial micro-phase separation of the two liquids larger, and thereby slower relaxing, structural entities of hydrogen bonded water and 2PGME molecules can only be formed at higher water contents than in the bulk system. However, the Q-dependence is unchanged with confinement, showing that the nature of the molecular motions is preserved. Thus, there is no indication of localization of the dynamics at length scales of less than 20 Å. The dynamics of both water and 2PGME is strongly dominated by translational diffusion at a temperature of 280 K.

Glass transition and relaxation dynamics of propylene glycol-water solutions confined in clay.

The molecular dynamics of aqueous solutions of propylene glycol (PG) and propylene glycol methylether (PGME) confined in a two-dimensional layer-structured Na-vermiculite clay has been studied by broadband dielectric spectroscopy and differential scanning calorimetry. As typical for liquids in confined geometries the intensity of the cooperative α-relaxation becomes considerably more suppressed than the more local ß-like relaxation processes. In fact, at high water contents the calorimetric glass transition and related structural α-relaxation cannot even be observed, due to the confinement. Thus, the intensity of the viscosity related α-relaxation is dramatically reduced, but its time scale as well as the related glass transition temperature Tg are for both systems only weakly influenced by the confinement. In the case of the PGME-water solutions it is an important finding since in the corresponding bulk system a pronounced non-monotonic concentration dependence of the glass transition related dynamics has been observed due to the growth of hydrogen bonded relaxing entities of water bridging between PGME molecules [J. Sjöström, J. Mattsson, R. Bergman, and J. Swenson, Phys. Chem. B 115, 10013 (2011)]. The present results suggest that the same type of structural entities are formed in the quasi-two-dimensional space between the clay platelets. It is also observed that the main water relaxation cannot be distinguished from the ß-relaxation of PG or PGME in the concentration range up to intermediate water contents. This suggests that these two processes are coupled and that the water molecules affect the time scale of the ß-relaxation. However, this is most likely true also for the corresponding bulk solutions, which exhibit similar time scales of this combined relaxation process below Tg. Finally, it is found that at higher water contents the water relaxation does not merge with, or follow, the α-relaxation above Tg, but instead crosses the α-relaxation, indicating that the two relaxation processes are independent of each other. This can only occur if the two processes do not occur in the same parts of the confined solutions. Most likely the hydration shell of the interlayer Na(+) ions is causing this water relaxation, which does not participate in the α-relaxation at any temperature.

The inhibition of fibril formation of lysozyme by sucrose and trehalose.

The two disaccharides, trehalose and sucrose, have been compared in many studies due to their structural similarity. Both possess the ability to stabilise and reduce aggregation of proteins. Trehalose has also been shown to inhibit the formation of highly structured protein aggregates called amyloid fibrils. This study aims to compare how the thermal stability of the protein lysozyme at low pH (2.0 and 3.5) is affected by the presence of the two disaccharides. We also address the anti-aggregating properties of the disaccharides and their inhibitory effects on fibril formation. Differential scanning calorimetry confirms that the thermal stability of lysozyme is increased by the presence of trehalose or sucrose. The effect is slightly larger for sucrose. The inhibiting effects on protein aggregation are investigated using small-angle X-ray scattering which shows that the two-component system consisting of lysozyme and water (Lys/H2O) at pH 2.0 contains larger aggregates than the corresponding system at pH 3.5 as well as the sugar containing systems. In addition, the results show that the particle-to-particle distance in the sugar containing systems (Lys/Tre/H2O and Lys/Suc/H2O) at pH 2.0 is longer than at pH 3.5, suggesting larger protein aggregates in the former. Finally, the characteristic distance separating ß-strands in amyloid fibrils is observed for the Lys/H2O system at pH 2.0, using wide-angle X-ray scattering, while it is not clearly observed for the sugar containing systems. This study further shows that the two disaccharides stabilise the native fold of lysozyme by increasing the denaturation temperature. However, other factors, such as a weakening of hydrophobic interactions and hydrogen bonding between proteins, might also play a role in their inhibitory effect on amyloid fibril formation.

Comparison of Sucrose and Trehalose for Protein Stabilization Using Differential Scanning Calorimetry.

The disaccharide trehalose is generally acknowledged as a superior stabilizer of proteins and other biomolecules in aqueous environments. Despite many theories aiming to explain this, the stabilization mechanism is still far from being fully understood. This study compares the stabilizing properties of trehalose with those of the structurally similar disaccharide sucrose. The stability has been evaluated for the two proteins, lysozyme and myoglobin, at both low and high temperatures by determining the glass transition temperature, Tg, and the denaturation temperature, Tden. The results show that the sucrose-containing samples exhibit higher Tden than the corresponding trehalose-containing samples, particularly at low water contents. The better stabilizing effect of sucrose at high temperatures may be explained by the fact that sucrose, to a greater extent, binds directly to the protein surface compared to trehalose. Both sugars show Tden elevation with an increasing sugar-to-protein ratio, which allows for a more complete sugar shell around the protein molecules. Finally, no synergistic effects were found by combining trehalose and sucrose. Conclusively, the exact mechanism of protein stabilization may vary with the temperature, as influenced by temperature-dependent interactions between the protein, sugar, and water. This variability can make trehalose to a superior stabilizer under some conditions and sucrose under others.

On the interactions between RNA and titrateable lipid layers: implications for RNA delivery with lipid nanoparticles.

Characterising the interaction between cationic ionisable lipids (CIL) and nucleic acids (NAs) is key to understanding the process of RNA lipid nanoparticle (LNP) formation and release of NAs from LNPs. Here, we have used different surface techniques to reveal the effect of pH and NA type on the interaction with a model system of DOPC and the CIL DLin-MC3-DMA (MC3). At only 5% MC3, differences in the structure and dynamics of the lipid layer were observed. Both pH and %MC3 were shown to affect the absorption behaviour of erythropoietin mRNA, polyadenylic acid (polyA) and polyuridylic acid (polyU). The adsorbed amount of all studied NAs was found to increase with decreasing pH and increasing %MC3 but with different effects on the lipid layer, which could be linked to the NA secondary structure. For polyA at pH 6, adsorption to the surface of the layer was observed, whereas for other conditions and NAs, penetration of the NA into the layer resulted in the formation of a multilayer structure. By comparison to simulations excluding the secondary structure, differences in adsorption behaviours between polyA and polyU could be observed, indicating that the NA's secondary structure also affected the MC3-NA interactions.

Different behavior of water in confined solutions of high and low solute concentrations.

Water-glycerol solutions confined in 21 Å pores of the silica matrix MCM-41 C10 have been studied using differential scanning calorimetry (DSC) and broadband dielectric spectroscopy (BDS). The results suggest a micro-phase separation caused by the confinement. Likely the water molecules coordinate to the hydroxyl surface groups of the pores, leaving most of the glycerol molecules in the centre of the pores. This makes the dynamics of glycerol almost concentration independent up to water concentrations of about 85 wt%. However, at higher water concentrations no substantial clustering of glycerol molecules should occur and the glass transition related dynamics exhibit an anomalous behaviour. Instead of a common plasticization effect of water, as for the corresponding bulk solutions (when no ice is formed), it is evident that water acts as an anti-plasticizer in the confinement at high water concentrations. We propose that the increased water concentration slows down the glass transition related dynamics in the deeply supercooled regime due to that a rigid hydrogen bonded network structure of water molecules is formed at low temperatures and low glycerol concentrations. This is in contrast to the situation in a homogenously mixed bulk solution of a high solute concentration where the water molecules will be less hydrogen bonded, and therefore are typically more mobile than the surrounding solute molecules. An almost complete hydrogen bonded network of water molecules may, even in confinements, be sufficiently rigid to slow down the relaxation of embedded solute molecules. It can also be expressed the other way around, i.e. small amounts of glycerol act as a plasticizer for water, due to its breaking up of the nearly tetrahedral network structure. From the here observed concentration dependent behaviour of the deeply supercooled bulk and confined solutions it seems, furthermore, evident that the Tg value of bulk water cannot be estimated from extrapolations of aqueous solutions.

The temperature dependent structure of liquid 1-propanol as studied by neutron diffraction and EPSR simulations.

The structure of liquid 1-propanol is investigated as a function of temperature using neutron diffraction together with Empirical Potential Structure Refinement modelling. The combined diffraction and computer modelling analysis demonstrates that propanol molecules form hydrogen bonded clusters with a relatively wide size distribution, which broadens at lower temperatures. We find that the cluster size distribution is well described by a recently proposed statistical model for branched H-bonded networks [P. Sillrén, J. Bielecki, J. Mattsson, L. Börjesson, and A. Matic, J. Chem. Phys. 136, 094514 (2012)]. The average cluster size increases from ~3 to 7 molecules, whilst the standard deviation of the size distribution increases from 3.3 to 8.5 as the temperature is decreased from 293 to 155 K. The clusters are slightly branched, with a higher degree of branching towards lower temperatures. An analysis of the cluster gyration tensor (R(mn)) reveals an average elongated ellipsoidal shape with axes having proportions 1:1.4:1.9. We find that the average radius of gyration has a cluster size dependence consistent with that of fractal clusters, R(g) ∝ n(1/D), with a fractal dimension D ≈ 2.20, which is close to D = 2.00 expected for an ideal random walk or D = 2.11 expected for reaction limited aggregation. The characteristic angles between the H-bonded OH-groups that constitute the clusters show only a weak temperature dependence with O-H···O angles becoming more narrowly distributed around 180° at lower temperatures.