Isotropic liquid state of cocoa butter.

The complex crystal structure of coca butter (CB) is responsible for the unique melting behavior, surface gloss, and mechanical properties of chocolate. While most studies concentrated on the crystalline state of CB, few studied the isotropic liquid state, which has a major impact on the crystallization process and the characteristics of the resulting crystals. In this study, the molecular organizations of the main CB triacylglycerols (TAGs; 1,3-dipalmitoyl-2-oleoylglycerol, palmitoyl-oleoyl-stearoylglycerol, POS, and 1,3-distearoyl-2-oleoylglycerol) were studied. The findings revealed the tunning-fork (Tf) conformation, commonly found in the crystalline state, is the least abundant in the isotropic liquid state of CB and pure TAGs. Notably, POS was found to interact with itself in CB, while its molecules with Tf conformation, although in small amounts in the mixture, tend to pair with each other at lower temperatures. These results highlight the significance of POS in CB crystallization and provide insights for developing CB alternatives.

Dehydroroyleanone as a Building Block for a Drug Delivery Platform Based on Self-Assembled Nanoparticles: Structural Studies and Chemical Modification.

6,7-Dehydroroyleanone (DHR) is a caspase-induced cytotoxic abietane diterpene, frequently found on Plectranthus spp. A pharmaceutical formulation consisting of a DHR-squalene conjugate was synthesized and analyzed by different techniques such as scanning electron microscopy (SEM). The facile production of the dispersion of DHR-squalene conjugate nanoparticles in phosphate buffer (pH 7.4) suggests that this nanodelivery platform may be an effective system to improve the solubility and bioavailability of DHR, so that therapeutical systemic levels may be achieved.

High-Performance Ionanofluids from Subzipped Carbon Nanotube Networks.

Investments in the transfer and storage of thermal energy along with renewable energy sources strengthen health and economic infrastructure. These factors intensify energy diversification and the more rapid post-COVID recovery of economies. Ionanofluids (INFs) composed of long multiwalled carbon nanotubes (MWCNTs) rich in sp2-hybridized atoms and ionic liquids (ILs) display excellent thermal conductivity enhancement with respect to the pure IL, high thermal stability, and attractive rheology. However, the influence of the morphology, physicochemistry of nanoparticles and the IL-nanostructure interactions on the mechanism of heat transfer and rheological properties of INFs remain unidentified. Here, we show that intertube nanolayer coalescence, supported by 1D geometry assembly, leads to the subzipping of MWCNT bundles and formation of thermal bridges toward 3D networks in the whole INF volume. We identified stable networks of straight and bent MWCNTs separated by a layer of ions at the junctions. We found that the interactions between the ultrasonication-induced breaking nanotubes and the cations were covalent in nature. Furthermore, we found that the ionic layer imposed by close MWCNT surfaces favored enrichment of the cis conformer of the bis(trifluoromethylsulfonyl)imide anion. Our results demonstrate how the molecular perfection of the MWCNT structure with its supramolecular arrangement affects the extraordinary thermal conductivity enhancement of INFs. Thus, we gave the realistic description of the interactions at the IL-CNT interface with its (super)structure and chemistry as well as the molecular structure of the continuous phase. We anticipate our results to be a starting point for more complex studies on the supramolecular zipping mechanism. For example, ionically functionalized MWCNTs toward polyionic systems─of projected and controlled nanolayers─could enable the design of even more efficient heat-transfer fluids and miniaturization of flexible electronics.

DLPGEN: Preparing Molecular Dynamics Simulations with Support for Polarizable Force Fields.

Preparing input files for molecular dynamics (MD) simulations can be a tedious task, particularly if different MD programs need to be used. Most simulation packages are accompanied by applications that handle this task, and, in some cases, software to perform interconversion between MD programs can be found. However, the conversion between different types of files is not always foolproof or the force field may not be fully supported, as quite often observed with polarization models. This work describes the program DLPGEN, which produces input files for the MD programs GROMACS, CHARMM, DL_POLY, and LAMMPS. The program can prepare polarizable force fields using a self-consistent field approach or with a dual thermostat-extended Lagrangian model. For GROMACS, a new polarization scheme suitable for the simulation of molecules containing more than one virtual particle is described. Results obtained for ethanol in the liquid state revealed that the system configurational energy, liquid density, and self-diffusion coefficient obtained with GROMACS are in good agreement with the data found with LAMMPS and CHARMM. In the case of DL_POLY, a problem with the Shells temperature control was found, suggesting that this program may not be suitable for the simulation of molecules containing multiple Drude particles.

Water Solubility Trends in Ionic Liquids: The Quantitative Structure-Property Relationship Model versus Molecular Dynamics.

The knowledge of water solubility in ionic liquids (ILs) is an important property with an impact on the design of many physical and chemical processes, like the purification of organic compounds or the establishment of decontamination procedures. The development of methods to predict or establish solubility trends in ILs is, therefore, extremely relevant, as it may avoid expensive and time-consuming experimental procedures. In this work, we compare results of water solubility in ILs predicted by a quantitative structure-property relationship (QSPR) model with trends found using aggregation studies in molecular dynamics (MD) simulation results. This study was performed for ILs combining the cations 1-butyl-1-methylpyrrolidinium and 1-butyl-1-methylmorpholinium, with the anions bis(pentafluoroethylsulfonyl)imide (BETI-), trifluoromethanesulfonate (TF-), and tetrafluoroborate (BF4-). Both methods indicated that, at 298.15 K, the water solubility in ILs was almost independent of the investigated cations. However, if the IL is composed of a hydrophobic anion, a slight increase in the mixability of the IL with water may be observed if the cation can form H-bonds. The QSPR model indicated that the hydrophobic BETI- anion leads to solubilities (xH2O ∼ 0.33), approximately half of those predicted when the cations are combined with TF- and BF4- anions (xH2O ∼ 0.60). The MD results suggested that this difference is essentially related to the ability of the water molecules to interact with the anion. This interaction involves the formation of networks of molecules, where H2O is completely solvated by anions. These structures make the formation of interactions between water molecules difficult, which are responsible for their segregation from solution and, therefore, to liquid-liquid phase separation. For the investigated ILs, the MD data also suggest that the solubility trends are inversely proportional to the number of "isolated" anions relative to ···AN-H2O-AN-H2O··· networks.

First and Second Dissociation Enthalpies in Bi-Component Crystals Consisting of Maleic Acid and L-Phenylalanine.

The energetics of the stepwise dissociation of a A:B2 bi-component crystal, according to A:B2(cr) → A:B(cr) + B(cr) and A:B(cr) → A(cr) + B(cr), was investigated using MA:Phe2 and MA:Phe (MA = maleic acid; Phe = L-phenylalanine) as model systems. The enthalpy changes associated with these sequential processes and with the overall dissociation reaction A:B2(cr) → A(cr) + 2B(cr) were determined by solution calorimetry. It was found that they are all positive, indicating that there is a lattice enthalpy gain when MA:Phe2 is formed, either from the individual precursors or by adding Phe to MA:Phe. Single-crystal X-ray diffraction (SCXRD) analysis showed that MA:Phe2 is best described as a protic salt containing a maleate anion (MA-) and two non-equivalent L-phenylalanine units, both linked to MA- by NH···O hydrogen bonds (H-bond): one of these units is protonated (HPhe+) and the other zwitterionic (Phe±). Only MA- and HPhe+ molecules are present in the MA:Phe lattice. In this case, however, NH···O and OH···O H-bonds are formed between each MA- unit and two HPhe+ molecules. Despite these structural differences, the enthalpy cost for the removal of the zwitterionic Phe± unit from the MA:Phe2 lattice to yield MA:Phe is only 0.9 ± 0.4 kJ mol-1 higher than that for the dissociation of MA:Phe, which requires a proton transfer from HPhe+ to MA- and the rearrangement of L-phenylalanine to the zwitterionic, Phe±, form. Finally, a comparison of the dissociation energetics and structures of MA:Phe and of the previously reported glycine maleate (MA:Gly) analogue indicated that parameters, such as the packing coefficient, density, hydrogen bonds formed, or fusion temperature, are not necessarily good descriptors of dissociation enthalpy or lattice enthalpy trends when bi-component crystals with different molecular composition are being compared, even if the stoichiometry is the same.

The Solubility of Gases in Ionic Liquids: A Chemoinformatic Predictive and Interpretable Approach.

This work comprises the study of solubilities of gases in ionic liquids (ILs) using a chemoinformatic approach. It is based on the codification, of the atomic inter-component interactions, cation/gas and anion/gas, which are used to obtain a pattern of activation in a Kohonen Neural Network (MOLMAP descriptors). A robust predictive model has been obtained with the Random Forest algorithm and used the maximum proximity as a confidence measure of a given chemical system compared to the training set. The encoding method has been validated with molecular dynamics. This encoding approach is a valuable estimator of attractive/repulsive interactions of a generical chemical system IL+gas. This method has been used as a fast/visual form of identification of the reasons behind the differences observed between the solubility of CO2 and O2 in 1-butyl-3-methylimidazolium hexafluorophosphate (BMIM PF6 ) at identical temperature and pressure (TP) conditions, The effect of variable cation and anion effect has been evaluated.

Tuning the miscibility of water in imide-based ionic liquids.

Liquid-liquid phase behavior measurements were performed for binary mixtures of water and ionic liquids (ILs) containing the same 1-ethyl-3-methylimidazolium ([C2mim]+) cation and different imide-based anions, having symmetric (bis(fluorosulfonyl)imide ([FSI]-)) or asymmetric structures (2,2,2-trifluoromethylsulfonyl-N-cyanoamide ([TFSAM]-) and 2,2,2-trifluoro-N-(trifluoromethylsulfonyl)acetamide ([TSAC]-)). An inversion of phase behavior was observed: while below ∼298 K, the miscibility of water in the studied ILs increases according to the order [C2mim][TSAC] < [C2mim][FSI] < [C2mim][NTf2], for temperatures above ∼303 K, the reverse trend is observed [C2mim][NTf2] < [C2mim][FSI] < [C2mim][TSAC]. In turn, above ∼306 K the [C2mim][TFSAM] is completely miscible with H2O in all ranges of concentrations. The obtained results also revealed an unusual water solubility variation of 11% in [C2mim][FSI], and 20% in [C2mim][TSAC], when the system temperature was changed by less than 1 K, around 298 K and 301 K, respectively. Molecular Dynamics (MD) simulations were used to understand the IL-water interactions and rationalize the experimental observations. These results suggested that the miscibility trends are mainly related to the ability of the water molecules to form water-anion and water-water aggregates in solution.

C13 - a new empirical force field to characterize the mechanical behavior of carbyne chains.

An accurate prediction of the mechanical behavior of long carbyne chains depends on the suitable modeling of bond alternation in these chains. While first-principles methods are a good approach, less computationally demanding empirical potentials are preferable for large carbyne-containing systems. AIREBO and Reax empirical potentials have extensively and successfully been used for simulating the mechanical behavior of graphene and carbon nanotubes. However, it remains unclear if these potentials can be directly applied in the accurate mechanical modeling of carbon nanostructures with sp hybridization, without re-parameterization. Here, a new force-field for carbyne, designated as C13 potential, that takes bond alternation into account, is presented. This new empirical potential was parameterized from ab initio calculations. Molecular dynamics (MD) simulations using the developed force-field are then conducted to determine the mechanical properties of carbyne chains under tensile loading, namely to assess their dependence on chain length and temperature. The bending stiffness of carbyne and its persistence length are also calculated. The results obtained are validated through comparison with results available in the literature. Lastly, the C13 potential is employed to model, for the first time, the tensile and the compressive behaviors of the hybrid system composed of carbon nanotubes infilled with carbyne chains.

Neat ionic liquids versus ionic liquid mixtures: a combination of experimental data and molecular simulation.

Simple mixtures of ionic liquids (IL-IL mixtures) can become a promising approach for the substitution of task-specific ILs. Such a concept was explored in this article by comparison of the thermophysical properties of neat 1-ethyl-3-methylimidazolium 2,2,2-trifluoromethylsulfonyl-N-cyanoamide, [C2mim][TFSAM], and equimolar mixtures of two structurally similar ILs having more common ions. Molecular dynamics (MD) simulations were additionally used to further highlight structural aspects of these systems at a molecular level.

Polymorphism in Simvastatin: Twinning, Disorder, and Enantiotropic Phase Transitions.

Simvastatin is one of the most widely used active pharmaceutical ingredients for the treatment of hyperlipidemias. Because the compound is employed as a solid in drug formulations, particular attention should be given to the characterization of different polymorphs, their stability domains, and the nature of the phase transitions that relate them. In this work, the phase transitions delimiting the stability domains of three previously reported simvastatin forms were investigated from structural, energetics, and dynamical points of view based on single crystal X-ray diffraction (SCXRD), hot stage microscopy (HSM), and differential scanning calorimetry (DSC) experiments (conventional scans and heat capacity measurements), complemented with molecular dynamics (MD) simulations. Previous assignments of the crystal forms were confirmed by SCXRD: forms I and II were found to be orthorhombic ( P212121, Z'/ Z = 1/4) and form III was monoclinic ( P21, Z'/ Z = 2/4). The obtained results further indicated that (i) the transitions between different forms are observed at 235.9 ± 0.1 K (form III → form II) and at 275.2 ± 0.2 K (form II → form I) in DSC runs carried out at 10 K min-1 and close to these values when other types of techniques are used (e.g., HSM). (ii) They are enantiotropic (i.e., there is a transition temperature relating the two phases before fusion at which the stability order is reversed), fast, reversible, with very little hysteresis between heating and cooling modes, and occur under single crystal to single crystal conditions. (iii) A nucleation and growth mechanism seems to be followed since HSM experiments on single crystals evidenced the propagation of an interface, accompanied by a change of birefringence and crystal contraction or expansion (more subtle in the case of form III → form II), when the phase transitions are triggered. (iv) Consistent with the reversible and small hysteresis nature of the phase transitions, the SCXRD results indicated that the molecular packing is very similar in all forms and the main structural differences are associated with conformational changes of the "ester tail". (v) The MD simulations further suggested that the tail is essentially "frozen" in two conformations below the III → II transition temperature, becomes progressively less hindered throughout the stability domain of form II, and acquires a large conformational freedom above the II → I transition. Finally, the fact that these transitions were found to be fast and reversible suggests that polymorphism is unlikely to be a problem for pharmaceutical formulations employing crystalline simvastatin because, if present, the III and II forms will readily convert to form I at ambient temperature.

Extraction Optimization and Structural and Thermal Characterization of the Antimicrobial Abietane 7α-Acetoxy-6ß-hydroxyroyleanone.

The abietane 7α-acetoxy-6ß-hydroxyroyleanone (AHR), obtained from plant extracts, is an attractive lead for drug development, given its known antimicrobial properties. Two basic requirements to establish any compound as a new drug are the development of a convenient extraction process and the characterization of its structural and thermal properties. In this work seven different methods were tested to optimize the extraction of AHR from Plectranthus grandidentatus. Supercritical fluid extraction (SFE) proved to be the method of choice, delivering an amount of AHR (57.351 µg·mg-1) approximately six times higher than the second best method (maceration in acetone; 9.77 µg·mg-1). Single crystal X-ray diffraction analysis of the ARH molecular and crystal structure carried out at 167 ± 2 K and 296 ± 2 K showed only a single phase, here dubbed form III (orthorhombic space group P21212), at those temperatures. The presence of two other polymorphs above room temperature was, however, evidenced by differential scanning calorimetry (DSC). The three forms are enantiotropically related, with the form III → form II and form II → form I transitions occurring at 333.5 ± 1.6 K and 352.0 ± 1.6 K, respectively. The fact that the transitions are reversible suggests that polymorphism is not likely to be an issue in the development pharmaceutical formulations based on ARH. DSC experiments also showed that the compound decomposes on melting at 500.8 ± 0.8 K. Melting should therefore be avoided if, for example, strategies to improve solubility based on the production of glassy materials or solid dispersions are considered.

Solvation of alcohols in ionic liquids - understanding the effect of the anion and cation.

In this work, we studied the effect of anion and cation properties on the interaction of alcohols with ionic liquids (ILs), using propan-1-ol as a molecular probe. The enthalpies of solution at infinite dilution of propan-1-ol in several ILs were measured by isothermal titration calorimetry (ITC). The calorimetric results were analysed together with molecular dynamics simulation and quantum chemical calculations of the interaction of the hydroxyl group of propan-1-ol with the anions. The results evidenced the role of the anion's basicity in the intermolecular interactions of alcohols and ionic liquids and further revealed a secondary effect of the cation nature on the solvation process. The effect of the anion basicity on the strength of the interaction of alcohols with ionic liquids was evaluated by comparing the results obtained for ILs with the same cation and different anions, [C4C1im][anion] (anions NTf2, PF6, FAP, DCA and TFA). The effect of the cation (size, aromaticity, charge distribution, and acidity) was explored using five different cations of the NTf2 series, [cation][NTf2] (cations C4C1im, C4C1pirr, C4py, C4C1pip, and C3C1C1im).

Modeling Halogen Bonds in Ionic Liquids: A Force Field for Imidazolium and Halo-Imidazolium Derivatives.

In this work, a force field for molecular dynamics and Monte Carlo simulations of ionic liquids containing imidazolium and halo-imidazolium derivatives is presented. This force field is an extension of the well-known CL&P and OPLS-AA models and was validated by comparing predicted crystalline structures for 22 ionic liquid compounds with the corresponding data deposited at the Cambridge Structural Database. The obtained results indicate that the proposed force field extension allows the reproduction of the crystal data with an absolute average deviation lower than 2.4%. Finally, it was also established that the halogen atoms covalently bound to the studied imidazolium cations are positively charged and do not exhibit a so-called σ-hole feature. For this reason, the formation of halogen bonds in the proposed force field appears naturally from the parametrized atomic point-charge distribution, without the necessity of any extra interaction sites.

Structure-property relationships in protic ionic liquids: a study of solvent-solvent and solvent-solute interactions.

The ionic nature of a functionalized protic ionic liquid cannot be rationalized simply through the differences in aqueous proton dissociation constants between the acid precursor and the conjugate acid of the base precursor. The extent of proton transfer, i.e. the equilibrium ionicity, of a tertiary ammonium acetate protic ionic liquid can be significantly increased by introducing an hydroxyl functional group on the cation, compared to the alkyl or amino-functionalized analogues. This increase in apparent ionic nature correlates well with variations in solvent-solute and solvent-solvent interaction parameters, as well as with physicochemical properties such as viscosity.

Ionic liquids with anions based on fluorosulfonyl derivatives: from asymmetrical substitutions to a consistent force field model.

Herein, seven anions including four imide-based, namely bis[(trifluoromethyl)sulfonyl]imide (TFSI), bis(fluorosulfonyl)imide (FSI), bis[(pentafluoroethyl)sulfonyl]imide (BETI), 2,2,2-trifluoromethylsulfonyl-N-cyanoamide (TFSAM) and 2,2,2-trifluoro-N-(trifluoromethylsulfonyl) acetamide (TSAC), and two sulfonate anions, trifluoromethanesulfonate (triflate, TF) and nonafluorobutanesulfonate (NF), are considered and compared. The volumetric mass density and dynamic viscosity of five ionic liquids containing these anions combined with the commonly used 1-ethyl-3-methylimidazolium cation (C2C1im), [C2C1im][FSI], [C2C1im][BETI], [C2C1im][TFSAM], [C2C1im][TSAC] and [C2C1im][NF] are measured in the temperature range of 293.15 ≤ T/K ≤ 353.15 and at atmospheric pressure. The results show that [C2mim][FSI] and [C2mim][TFSAM] exhibit the lowest densities and viscosities among all the studied ionic liquids. The experimental volumetric data is used to validate a more consistent re-parameterization of the CL&P force field for use in MD simulations of ionic liquids containing the ubiquitous bis[(trifluoromethyl)sulfonyl]imide and trifluoromethanesulfonate anions and to extend the application of the model to other molten salts with similar ions.

Structure-property relationships in protic ionic liquids: a thermochemical study.

How does cation functionality influence the strength of intermolecular interactions in protic ionic liquids (PILs)? Quantifying the energetics of PILs can be an invaluable tool to answer this fundamental question. With this in view, we have determined the standard molar enthalpy of vaporization, , and the standard molar enthalpy of formation, , of three tertiary ammonium acetate PILs with varying cation functionality, and of their corresponding precursor amines, through a combination of Calvet-drop microcalorimetry, solution calorimetry, and ab initio calculations. The obtained results suggest that these PILs vaporize as their neutral acid and base precursors. We also found a strong correlation between of the PILs and of their corresponding amines. This suggests that, within this series of PILs, the influence of cation modification on their cohesive energies follows a group additivity rule. Finally, no correlation between the of PILs and the extent of proton transfer, as estimated from the difference in aqueous pKa between the precursor acid and the conjugate acid of the precursor base, was observed.

ForConX: A forcefield conversion tool based on XML.

The force field conversion from one MD program to another one is exhausting and error-prone. Although single conversion tools from one MD program to another exist not every combination and both directions of conversion are available for the favorite MD programs Amber, Charmm, Dl-Poly, Gromacs, and Lammps. We present here a general tool for the force field conversion on the basis of an XML document. The force field is converted to and from this XML structure facilitating the implementation of new MD programs for the conversion. Furthermore, the XML structure is human readable and can be manipulated before continuing the conversion. We report, as testcases, the conversions of topologies for acetonitrile, dimethylformamide, and 1-ethyl-3-methylimidazolium trifluoromethanesulfonate comprising also Urey-Bradley and Ryckaert-Bellemans potentials. © 2017 Wiley Periodicals, Inc.

AGGREGATES: Finding structures in simulation results of solutions.

Molecular Dynamic and Monte-Carlo simulations are widely used to investigate the structure and physical properties of solids and liquids at a molecular level. Tools to extract the most relevant information from the obtained results are, however, in considerable demand. One such tool, the program AGGREGATES, is described in this work. Based on distance criteria, the program searches trajectory files for the presence of molecular clusters and computes several statistical and shape properties for these structures. Tools designed to investigate the local organization and the molecular conformations in the clusters are also available. Among these, it is introduced a new approach to perform a First Shell Analysis, by looking for the presence of atomic contacts between molecules. These elements are particularly useful to obtain information on molecular assembly processes (such as the nucleation of crystals or colloidal particles) or to investigate polymorphism in organic compounds. The program features are illustrated here through the investigation of the 4'-hydroxyacetophenone + ethanol system. © 2017 Wiley Periodicals, Inc.

Mixtures of the 1-ethyl-3-methylimidazolium acetate ionic liquid with different inorganic salts: insights into their interactions.

In this work, we explore the interactions between the ionic liquid 1-ethyl-3-methylimidazolim acetate and different inorganic salts belonging to two different cation families, those based on ammonium and others based on sodium. NMR and Raman spectroscopy are used to screen for changes in the molecular environment of the ions in the ionic liquid + inorganic salt mixtures as compared to pure ionic liquid. The ion self-diffusion coefficients are determined from NMR data, allowing the discussion of the ionicity values of the ionic liquid + inorganic salt mixtures calculated using different methods. Our data reveal that preferential interactions are established between the ionic liquid and ammonium-based salts, as opposed to sodium-based salts. Computational calculations show the formation of aggregates between the ionic liquid and the inorganic salt, which is consistent with the spectroscopic data, and indicate that the acetate anion of the ionic liquid establishes preferential interactions with the ammonium cation of the inorganic salts, leaving the imidazolium cation less engaged in the media.