Empirical Potential Structure Refinement Technique Research Articles

Modeling the Structure of Ketoprofen-Poly(vinylpyrrolidone) Amorphous Solid Dispersions with Empirical Potential Structure Refinements of X-ray Scattering Data.

The metastability of amorphous formulations poses barriers to their safe and widespread commercialization. The propensity of amorphous solid dispersions (ASDs) to crystallize is directly linked to their molecular structure. Amorphous structures are inherently complex and thus difficult to fully characterize by experiments, which makes structural simulations an attractive route for investigating which structural characteristics correlate with ASD stability. In this study, we use empirical potential structure refinement (EPSR) to create molecular models of ketoprofen-poly(vinylpyrrolidone) (KTP/PVP) ASDs with 0-75 wt % drug loading. The EPSR technique uses X-ray total scattering measurements as constraints, yielding models that are consistent with the X-ray data. We perform several simulations to assess the sensitivity of the EPSR approach to input parameters such as intramolecular bond rotations, PVP molecule length, and PVP tacticity. Even at low drug loading (25 wt %), ∼40% of KTP molecules participate in KTP-KTP hydrogen bonding. The extent of KTP-PVP hydrogen bonding does not decrease significantly at higher drug loadings. However, the models' relative uncertainties are too large to conclude whether ASDs' lower stabilities at high drug loadings are due to changes in drug-excipient hydrogen bonding or a decrease in steric hindrance of KTP molecules. This study illustrates how EPSR, combined with total scattering measurements, can be a powerful tool for investigating structural characteristics in amorphous formulations and developing ASDs with improved stability.

Hydrophobic hydration of the hydrocarbon adamantane in amorphous ice.

Hydrophobic molecules are by definition difficult to hydrate. Previous studies in the area of hydrophobic hydration have therefore often relied on using amphiphilic molecules where the hydrophilic part of a molecule enabled the solubility in liquid water. Here, we show that the hydrophobic adamantane (C10H16) molecule can be fully hydrated through vapour codeposition with water onto a cryogenic substrate at 80 K resulting in the matrix isolation of adamantane in amorphous ice. Using neutron diffraction in combination with the isotopic substitution method and the empirical potential structure refinement technique, we find that the first hydration shell of adamantane is well structured consisting of a hydrogen-bonded cage of 28 water molecules that is also found in cubic structure II clathrate hydrates. The four hexagonal faces of the 51264 cage are situated above the four methine (CH) groups of adamantane whereas the methylene (CH2) groups are positioned below the edges of two adjoining pentagonal faces. The oxygen atoms of the 28 water molecules can be categorised on the basis of symmetry equivalences as twelve A, twelve B and four C oxygens. The water molecules of the first hydration shell display orientations consistent with those expected for a clathrate-hydrate-type cage, but also unfavourable ones with respect to the hydrogen bonding between the water molecules. Annealing the samples at 140 K, which is just below the crystallisation temperature of the matrix, removes the unfavourable orientations and leads to a slight increase in the structural order of the first hydration shell. The very closest water molecules display a tendency for their dipole moments to point towards the adamantane which is attributed to steric effects. Other than this, no significant polarisation effects are observed which is consistent with weak interactions between adamantane and the amorphous ice matrix. FT-IR spectroscopy shows that the incorporation of adamantane into amorphous ice leads to a weakening of the hydrogen bonds. In summary, the matrix-isolation of the highly symmetric adamantane in amorphous ice provides an interesting test case for hydrophobic hydration. Studying the structure and spectroscopic properties of water at the interface with hydrophobic hydrocarbons is also relevant for astrophysical environments, such as comets or the interstellar medium, where amorphous ice and hydrocarbons have been shown to coexist in large quantities.

Local structure and orientational ordering in liquid bromoform

The neutron diffraction data of liquid bromoform (CHBr3) at 25°C was analysed using the Empirical Potential Structure Refinement technique in combination with H/D isotopic substitution. Compared to liquid chloroform (CHCl3), CHBr3 displays more spatially defined intermolecular contacts. A preference for polar stacking with collinear alignment of dipole moments is observed for the most closely approaching CHBr3 molecules, although to a lesser extent than in chloroform. Consistent with this, and in line with dielectric spectroscopy, the Kirkwood correlation factor from the structural model of CHBr3 is smaller than that of CHCl3. The net antiparallel alignment of dipole moments in CHBr3, as suggested by dielectric spectroscopy, must be due to weak but persistent long-range orientation correlations in CHBr3, which counteract the local polar stacking.

The Double-Faced Nature of Hydrogen Bonding in Hydroxy-Functionalized Ionic Liquids Shown by Neutron Diffraction and Molecular Dynamics Simulations.

We characterize the double-faced nature of hydrogen bonding in hydroxy-functionalized ionic liquids by means of neutron diffraction with isotopic substitution (NDIS), molecular dynamics (MD) simulations, and quantum chemical calculations. NDIS data are fit using the empirical potential structure refinement technique (EPSR) to elucidate the nearest neighbor H⋅⋅⋅O and O⋅⋅⋅O pair distribution functions for hydrogen bonds between ions of opposite charge and the same charge. Despite the presence of repulsive Coulomb forces, the cation-cation interaction is stronger than the cation-anion interaction. We compare the hydrogen-bond geometries of both "doubly charged hydrogen bonds" with those reported for molecular liquids, such as water and alcohols. In combination, the NDIS measurements and MD simulations reveal the subtle balance between the two types of hydrogen bonds: The small transition enthalpy suggests that the elusive like-charge attraction is almost competitive with conventional ion-pair formation.

Application of the Empirical Potential Structure Refinement Technique to a Borosilicate Glass of Nuclear Interest

We present how the Empirical Potential Structure Refinement (EPSR) technique developed by A. Soper and co-workers can be applied to borosilicate glasses of nuclear interest to build realistic atomistic models. In particular, we show how the atomistic configuration for a six-oxide glass is modified by the introduction of the empirical potential during the simulation. EPSR strengths and weaknesses are then brought out by comparing simulation results with experimental data.

Preference for Isolated Water Molecules in a Concentrated Glycerol–Water Mixture

Neutron diffraction coupled with hydrogen/deuterium isotopic substitution has been used to investigate the structure of a concentrated glycerol water (4:1 mole fraction) solution. The neutron diffraction data were used to constrain a three-dimensional computational model that is experimentally relevant using the empirical potential structure refinement technique. From interrogation of this model, we find that glycerol-glycerol hydrogen bonding is largely unperturbed by the presence of water in the solution. We find that glycerol-water hydrogen bonding is prevalent, suggesting that water molecules effectively take the place of glycerol molecules in this concentrated solution. In contrast, we find that water-water hydrogen bonding is significantly perturbed. While the first coordination shell of water in the concentrated solution remains similar to that of pure water, water-water hydrogen bonding is greatly diminished beyond the first neighbor distance. Interestingly, the majority of water molecules exist as single monomers in the concentrated glycerol solution. The preference of isolated water molecules results in a solution that is well mixed with optimal glycerol-water hydrogen bonding. These results highlight the importance of preferential hydrogen bonding in aqueous solutions and suggest a mechanism for cryoprotection by which glycerol effectively hydrogen bonds with water, resulting in a disrupted hydrogen-bonded water network.

Influence of Concentration and Anion Size on Hydration of H+ Ions and Water Structure

Neutron diffraction experiments with hydrogen isotope substitution on aqueous solutions of HCl and HBr have been performed at concentrations ranging from 1:17 to 1:83 solute per water molecules, at ambient conditions. Data are analyzed using the empirical potential structure refinement technique in order to extract information on both the ion hydration shells and the microscopic structure of the solvent. It is found that the influence of these solutes on the water structure is less concentration dependent than that of salts or hydroxides. Moreover protons readily form a strong H-bond with a water molecule upon solvation, at all proportions. The majority of them is also bonded via a longer bond to another water molecule, giving a prepeak in the g(OwOw). At high solute concentration, the second water molecule may be substituted by the counterion. In particular at solute concentrations of the order of 1:17 or higher, all protons have an anion within a distance of 4.5 A.

Hydration of Sodium, Potassium, and Chloride Ions in Solution and the Concept of Structure Maker/Breaker

Neutron diffraction data with hydrogen isotope substitution on aqueous solutions of NaCl and KCl at concentrations ranging from high dilution to near-saturation are analyzed using the Empirical Potential Structure Refinement technique. Information on both the ion hydration shells and the microscopic structure of the solvent is extracted. Apart from obvious effects due to the different radii of the three ions investigated, it is found that water molecules in the hydration shell of K+ are orientationally more disordered than those hydrating a Na+ ion and are inclined to orient their dipole moments tangentially to the hydration sphere. Cl- ions form instead hydrogen-bonded bridges with water molecules and are readily accommodated into the H-bond network of water. The results are used to show that concepts such as structure maker/breaker, largely based on thermodynamic data, are not helpful in understanding how these ions interact with water at the molecular level.

Structure and interactions in simple solutions

Neutron scattering with hydrogen/deuterium isotopic substitution techniques has been used to investigate the full range of structural interactions in a dilute 0.02 mol fraction solution of tertiary butanol in water, both in the absence and in the presence of a small amount of sodium chloride. Emphasis is given to the detailed pictures of the intermolecular interactions that have been derived using the empirical potential structure refinement technique. Analysis has been performed to the level of the spatial density distribution functions that illustrate the orientational dependence of the intermolecular interactions between all combinations of molecular and ionic components. The results show the key structural motifs involved in the interactions between the various components in a complex aqueous system. They underline the structural versatility of the water molecule in accommodating a range of different kinds of interactions while retaining its characteristic first-neighbour interaction geometry. Within this framework, the results highlight the complex interplay between the polar, non-polar and charged molecular interactions that exist in the system.

The structure of a concentrated aqueous solution of tertiary butanol: Water pockets and resulting perturbations

The structure of a 0.86 mole fraction solution of tertiary butanol in water at room temperature has been investigated by neutron diffraction with hydrogen/deuterium isotope substitution. The intermolecular correlations between the alcohol and water molecules have been extracted using the empirical potential structure refinement technique. The results highlight the creation of water pockets within the solution structure. These pockets mediate the interactions that occur between the alcohol hydroxyl groups, reducing the level of direct alcohol–alcohol hydrogen bonding. The results show that only a small amount of water is required to drive this solution in the direction of hydrophobic behavior more commonly associated with the water-rich solution compositions.

Determination of the orientational pair correlation function of a molecular liquid from diffraction data

Abstract The structure of a molecular liquid is defined by means of the orientational pair correlation function (OPCF), g( r ,ω 1 ,ω 2 ) , which determines the likelihood of a molecule being found at position r, with orientation ω2, given a molecule with orientation ω1 at the origin. Diffraction experiments and most computer simulations of liquids however present only the site-site radial distribution functions, which by definition contain less information than the OPCF. By expanding the OPCF as a series of generalised spherical harmonic functions (SHARM), and by showing that this series, even when truncated, contains important orientational information, it is proposed that the coefficients to this series should provide the essential database for the structure of a molecular liquid as derived from a given set of diffraction data. The most reliable method of estimating these coefficients is to run an auxilliary computer simulation which incorporates all known physical constraints on molecular positions as well as reproducing the measured radial distribution functions or diffraction data as accurately as possible. The Empirical Potential Structure Refinement (EPSR) method of performing the simulation is a robust method of obtaining these structural coefficients and the results appear to be relatively insensitive to the details of the assumed reference intermolecular potential. The possibility of using the EPSR technique as a route to establishing an accurate intermolecular potential is also discussed. Some examples of applying the method to liquid water are presented.