Interfacial properties of the ionic liquid [bmim][triflate] over a wide range of temperatures

The potential of supramolecular synthon to develop coamorphous systems with tailored physical stability: Mechanistic insights integrating kinetics and thermodynamics

Coamorphous drug delivery systems have received increasing interest owing to their potential to improve the solubility, dissolution and bioavailability of poorly water-soluble drugs. However, the crystallization risk is one of major limitations in their application. It has been widely recognized that the coformer plays a vital role in physical stability of coamorphous formulation. Unfortunately, the screen of optimal coformer still adopts a trial-and-error method, which is time-consuming and expensive. Herein, a supramolecular synthon approach based on the interaction between functional groups, was exploited to design coamorphous systems (CMs) consisting of lurasidone hydrochloride (LH) and three coformers, saccharin (SAC), L-tryptophan (TRP), and L-cysteine hydrochloride (CYS). X-ray powder diffraction suggested the order of physical stability of the coamorphous systems was ranked as LH-CYS CM > LH-TRP CM > LH-SAC CM. The charge-assisted hydrogen bond between LH and coformer was confirmed by infrared spectroscopy and solid-state 13C NMR. Moreover, structural, electronic, and molecular interaction information, especially hydrogen bonding interactions, were quantified by theoretical calculations, including miscibility calculations, molecular dynamics simulations and quantum chemical calculations. It was revealed that LH-CYS CM exhibited the best miscibility, strongest binding energy and strongest H-bond with partially covalent character, demonstrating the significant role of supramolecular synthon in stabilizing coamorphous formulations. Interestingly, LH-TRP CM, not LH-CYS CM, exhibited the lowest molecular mobility among three coamorphous systems, which was inconsistent with their physical stability. But from thermodynamic perspective, the order of configurational entropy and physical stability of coamorphous systems was completely consistent. We shed light on the comprehensive effects of molecular mobility and configurational entropy on physical stability of coamorphous systems. Importantly, the relationship between supramolecular synthon and kinetic/thermodynamic mechanisms was also discussed, and the positive correlation between configurational entropy and intermolecular interactions was proposed in this paper. Our findings demonstrated the great potential of supramolecular synthon in designing coamorphous systems with tailored physical stability, and further provided a deeper insight into the mechanisms of physical stability of coamorphous systems.

Ordering of ionic liquids at a charged sapphire interface: Evolution with cationic chain length

HYPOTHESIS

Room Temperature Ionic Liquids (RTILs) bulk's molecular layering dominates their structure also at the RTIL/sapphire interface, increasing the layer spacing with the cationic alkyl chain length n. However, the negatively-charged sapphire surface compresses the layers, increases the layering range, and affects the intra-layer structure in yet unknown ways.

EXPERIMENTS

X-ray reflectivity (XR) off the RTIL/sapphire interface, for a broad homologous RTIL series 1-alkyl-3-methylimidazolium bis(trifluoromethansulfonyl)imide, hitherto unavailable for any RTIL.

FINDINGS

RTIL layers against the sapphire, exhibit two spacings: da and db. da is n-varying, follows the behavior of the bulk spacing but exhibits a downshift, thus showing significant layer compression, and over twofold polar slab thinning. The latter suggests exclusion of anions from the interfacial region due to the negative sapphire charging by x-ray-released electrons. The layering range is larger than the bulk's. db is short and near n-independent, suggesting polar moieties' layering, the coexistence mode of which with the da-spaced layering is unclear. Comparing the present layering with the bulk's and the RTIL/air interface's provides insight into the Coulomb and dispersion interaction balance dominating the RTIL's structure and the impact thereon of the presence of a charged solid interface.

Mechanistic insights into the crystallization of coamorphous drug systems

In our previous study, the coamorphous formulation of lurasidone hydrochloride (LH) with saccharin (SAC) showed significantly enhanced dissolution and physical stability compared to crystalline/amorphous LH. However, the coamorphous system is still in amorphous state, and has the tendency to recrystallization, which will in turn result in the loss of above advantages. In this study, the crystallization kinetics under isothermal and non-isothermal conditions was investigated. Compared to amorphous LH, coamorphous LH-SAC showed 68.3-361.2 and 2.6-6.1 times lower crystallization rates in glassy state and supercooled liquid state, respectively. After co-amorphization, the addition of SAC changed the crystallization mechanism of amorphous LH from nucleation-controlled to diffusion-controlled manner. Amorphous LH followed the site-saturated nucleation, whereas the coamorphous system exhibited a fixed number of nuclei. The non-isothermal crystallization indicated amorphous LH and coamorphous LH-SAC showed two-dimensional (JMAEK 2) and three-dimensional (JMAEK 3) growth of nuclei, respectively. Furthermore, coamorphous LH-SAC exhibited higher molecular mobility and dynamic fragility (m_D) than amorphous LH, which is kinetically unfavorable for its physical stability. However, from thermodynamic perspective, coamorphous LH-SAC had a higher configurational entropy, i.e., a higher entropy barrier for crystallization, which is beneficial to hinder its crystallization. Therefore, it was concluded that the higher configurational entropy rather than the molecular mobility was proposed to be responsible for its improved stability. In addition, molecular dynamics simulations with miscibility, radial distribution function and binding energy calculations suggested coamorphous components exhibited good miscibility and strong intermolecular interactions, which was also conductive to the enhancement in its stability. This study offers an in-depth understanding about the effect of the coformer on the crystallization kinetics of coamorphous systems, and points out the important contribution of the configurational entropy in stabilizing the coamorphous systems.

Interfacial Forces in Free-Standing Layers of Melted Polyethylene, from Critical to Nanoscopic Thicknesses

Molecular dynamics simulations of ultrathin free-standing layers made of melted (373.15–673.15 K) polyethylene chains, which exhibit a lower melting temperature (compared to the bulk value), were carried out to investigate the dominant pressure forces that shape the conformation of chains at the interfacial and bulk liquid regions. We investigated layer thicknesses, tL, from the critical limit of mechanical stability up to lengths of tens of nm and found a normal distribution of bonds dominated by slightly stretched chains across the entire layer, even at large temperatures. In the bulk region, the contribution of bond vibrations to pressure was one order of magnitude larger than the contributions from interchain interactions, which changed from cohesive to noncohesive at larger temperatures just at a transition temperature that was found to be close to the experimentally derived onset temperature for thermal stability. The interchain interactions produced noncohesive interfacial regions at all temperatures in both directions (normal and lateral to the surface layer). Predictions for the value of the surface tension, γ, were consistent with experimental results and were independent of tL. However, the real interfacial thickness—measured from the outermost part of the interface up to the point where γ reached its maximum value—was found to be dependent on tL, located at a distance of 62 Å from the Gibbs dividing surface in the largest layer studied (1568 chains or 313,600 bins); this was ~4 times the length of the interfacial thickness measured in the density profiles.

Gelation switch of polyamorphic indomethacin depending on the thermal procedure

Amorphous indomethacin (IMC) prepared under different thermal procedures via melt quenching method showed significantly different dissolution behaviors. This study aims to investigate the influence of thermal procedures on the formation of IMC polyamorphism and to explore the mechanism for their different dissolution behaviors. Amorphous IMC samples were prepared by melting crystalline IMC under a series of temperatures (160-195 °C), respectively, followed by quenching in liquid nitrogen. Samples obtained under 170 °C exhibited bi-halo shapes at ∼15° and ∼26° (2θ), while the ones above 175 °C showed a single halo at ∼21° (2θ), suggesting amorphous IMC prepared under different thermal procedures probably have different local molecular arrangements. In comparison to crystalline IMC, amorphous IMC obtained under 170 °C showed significantly higher dissolution profiles with good dispersibility in aqueous medium, however, all amorphous IMC samples prepared above 175 °C demonstrated much lower dissolution with significant gelation, which seemed like a gelation switch existed for polyamorphic IMC when the preparation temperature was between 170 and 175 °C. Based on physicochemical characterizations, amorphous IMC prepared under 170 °C had higher surface free energy, more surficial hydrophilic groups and better wettability than the ones made above 175 °C. Molecular dynamics simulations revealed that the amorphous samples prepared below 170 °C had similar binding energy values in the range of 310.045-325.479 kcal/mol, while those prepared above 175 °C were significantly lower within 212.193-235.073 kcal/mol. Such binding energy difference might be responsible for their different local molecular arrangements after different thermal procedures. The current study deeply reminds us that the thermal procedure of preparation methods may significantly affect the physicochemical properties of amorphous materials, which should be paid special attention to the polymorphic selection during pharmaceutical development.

Critical Thickness of Free-Standing Nanothin Films Made of Melted Polyethylene Chains via Molecular Dynamics

The mechanical stability of nanothin free-standing films made of melted polyethylene chains was predicted via molecular dynamics simulations in the range of 373.15-673.15 K. The predicted critical thickness, tc, increased with the square of the temperature, T, with additional chains needed as T increased. From T = 373.15 K up to the thermal limit of stability for polyethylene, tc values were in the range of nanothin thicknesses (3.42-5.63 nm), which approximately corresponds to 44-55 chains per 100 nm2. The density at the center of the layer and the interfacial properties studied (density profiles, interfacial thickness, and radius of gyration) showed independence from the film thickness at the same T. The polyethylene layer at its tc showed a lower melting T (<373.15 K) than bulk polyethylene.

Cation folding and the thermal stability limit of the ionic liquid [BMIM+][BF4 -] under total vacuum

Molecular dynamics simulations reveal the behavior of the bimodal distribution of cation conformations (folded/unfolded) in ionic liquids based on alkylated imidazoles, such as [BMIM+][BF4 -]. The alkyl chains of the cations can fold and block interactions between the cations and anions, thereby reducing the cohesivity of the liquid. At room temperature, the folded conformations represent less than one-third of the total conformations. In contrast to the behavior observed during the thermal denaturation of proteins, in ionic liquids, the concentration of folded cations grows when the temperature increases. At the equimolar concentration, the system reaches the reported experimental temperature of thermal stability (similar to the thermal denaturation behavior). There is an outermost layer of cations at the interface that can tilt toward the interface and cover a layer of anions adsorbed at the interface. This interfacial conformation makes the system stable in transverse directions and unstable in the normal direction at temperatures in the region of thermal instability, limiting the rate of vaporization of neutral ion pairs, which are observed as rare events at temperatures as low as 773.15 K.

Comment on "Bi-layering at ionic liquid surfaces: a sum - frequency generation vibrational spectroscopy - and molecular dynamics simulation-based study" by T. Iwahashi, T. Ishiyama, Y. Sakai, A. Morita, D. Kim and Y. Ouchi, Phys. Chem. Chem. Phys., 2020, 22, 12565

This Comment raises several questions concerning the surface structure concluded in the paper referenced in the title. Specifically, that paper ignores previous experiments and simulations which demonstrate for the same ionic liquids depth-decaying, multilayered surface-normal density profiles rather than the claimed molecular mono- or bi-layers. We demonstrate that the claimed structure does not reproduce the measured X-ray reflectivity, which probes directly the surface-normal density profile. The measured reflectivities are found, however, to be well-reproduced by a multilayered density model. These results, and previous experimental and simulation results, cast severe doubt on the validity of the surface structure claimed in the paper referenced in the title.

Evidence of Direct Dissolution of CO2 into the Ionic Liquid [C4min] [NTf2] during Their Initial Interaction

Ionic liquids (ILs) are known for their high ability to capture CO₂. However, the mechanism of CO₂ solubility into ILs during their initial interaction remains controversial. In this study, we analyzed the initial dissolution of CO₂ into an IL 1-butyl-3-methylimidazolium bis-(trifluoromethylsulfonyl)imide ([C₄min] [NTf₂]) by measuring its solubility using a combination of a molecular beam and a flowing liquid jet sheet beam (FJSB) and the King and Wells method (KW method). The temperature dependence of the initial dissolution probability indicates that the solubility of CO₂ in the IL [C₄min] [NTf₂] increases with increasing temperature. This result is not consistent with what has been reported in an equilibrium state. The initial dissolution probability was well-fitted by the Vogel-Fulcher-Tammann (VFT) equation, which describes the dynamical cage structure in IL systems. We also find that the initial dissolution probability was correlated to the cage lifetime and correlation length. The simple model of CO₂ dissolution into an IL with the cage model was implemented to explain the experimental results in this study. Our results indicate that the initial dissolution of CO₂ into the IL corresponds to a direct solution and not an uptake process.

Probing the electronic structure of ether functionalised ionic liquids using X-ray photoelectron spectroscopy

The charge distribution associated with individual components in functionalised ionic liquids (ILs) can be tuned by careful manipulation of the substituent groups incorporated into the ions. Here we use X-ray photoelectron spectroscopy to investigate the impact of substituent atoms on the electronic structure of similar imidazolium-based systems each paired with a common anion, [Tf₂N]^-. The experimental measurements revealed an unexpected variation in the charge density distribution within the IL cation when the oxygen atom in a poly-ether containing side chain is moved by just one atomic position. This surprising observation is supported by density functional theory calculations.

Reducing uncertainty in simulation estimates of the surface tension through a two-scale finite-size analysis: thicker is better

Recent simulation studies of the surface tension γ, and other properties of thin free-standing films, have revealed unexpected finite size effects in which the variance of the properties vary monotonically with the in-plane width of the films, complicating the extrapolation of estimates of film properties to the thermodynamic limit. We carried out molecular dynamics simulations to determine the origin of this phenomenon, and to address the practical problem of developing a more reliable methodology for estimating γ in the thermodynamic limit. We find that there are two distinct finite size effects that must be addressed in a finite size analysis of γ in thin films. The first finite size scale is the in-plane width of the films and the second scale is the simulation cell size in the transverse direction. Increasing the first scale enhances fluctuations in γ, measured by the standard deviation of their distribution, while increasing the second reduces γ fluctuations due to a corresponding increased 'freedom' of the film to fluctuate out of plane. We find that using progressively large simulation cells in the transverse direction, while keeping the film width fixed to an extent in which the full bulk liquid zone is developed, allows us to obtain a smooth extrapolation to the thermodynamic limit, enabling a reduction of the γ uncertainty to a magnitude on the order of 1% for systems having a reasonably large size, i.e., O (1 μm).