A structural and theoretical study of the alkylammonium nitrates forefather: Liquid methylammonium nitrate

Ionic Twin Nanostructural Comparison: Propylammonium Butanoate vs. Butylammonium Propanoate and Their Interactions with Water

This study investigates the nanostructure of two protic ionic liquids (PILs), [N0 0 0 3][C3CO2] and [N0 0 0 4][C2CO2], with similar polar head groups but varying alkyl chain lengths. An X-ray scattering technique and molecular dynamics simulations have been utilized to characterize the bulk and interfacial properties of these PILs. The findings suggest that the nanostructure of the PILs is primarily determined by the electrostatic forces between charged functional groups playing a dominant role. Despite differences in the alkyl chain lengths, the PILs possess remarkably similar nanostructures. Extending our investigation, we report the impact of water on the nanostructure. Our findings reveal that the addition of water disrupts interactions between cations and anions, weakening Coulombic forces. The disruptive behavior is attributed to the establishment of hydrogen bonds between water and ions. This comprehensive approach provides valuable insights into the nuanced factors shaping the nanostructure of these PILs, which are crucial for tailoring their applications in synthetic chemistry, catalysis, and beyond.

Highlighting the difference in nanostructure between domain-forming and domainless protic ionic liquids

Nanoheterogeneity in some ionic liquids is a known phenomenon, but quantifying or sometimes even identifying it is not a straightforward task. We compared several known and suggested some novel approaches to identify and characterize domain segregation using the results of atomistic simulations. 10 ammonium-based protic ionic liquids with different propensity to form segregated polar and apolar domains as suggested by experimental studies were considered. They include butyl-, propyl-, 2-methoxyethylammonium nitrate, butyl- and propylammonium hydrogen sulfate, butylammonium thiocyanate (domain-forming liquids), ethylammonium and pyrrolidinium nitrate (weakly pronounced segregation), methylammonium and 2-hydroxyethylammonium nitrate (domainless liquids). Molecular dynamics simulations were performed using models based on the OPLS-AA force field with scaled ion charges. Results show that domains can be recognized and the characteristic domain length scale can be determined from peaks of Ripley's functions, peaks and large-period oscillations of finite-volume radial distribution function integral, or difference of such integrals for polar and apolar atoms, and peaks of local atom density variance. These peaks disappear with increasing temperature due to the disruption of segregated domains. In domain-forming liquids, apolar atoms are more homogeneously distributed in space than polar atoms. In addition, the probability of molecular-sized cavity formation is significantly higher in apolar domains, which determines better solubility of apolar species in domain-forming ILs. The suggested approaches can be applied to various nanostructured liquids including both ionic and molecular solvents and mixtures, as well as other systems with mesoscale ordering.

Automated parameterization of quantum-mechanically derived force-fields including explicit sigma holes: A pathway to energetic and structural features of halogen bonds in gas and condensed phase

In classical molecular dynamics, general purpose atomistic force-fields (FFs) often deliver inaccurate results when dealing with halogen bonds (XBs), notwithstanding their crucial role in many fields of science, ranging from material design to drug development. Given the large dimensions of the systems of interest, it would be therefore desirable to increase the FF accuracy maintaining the simplicity of the standard Lennard-Jones (LJ) plus point charge description to avoid an excessive computational cost. A simple yet effective strategy consists in introducing a number of virtual sites able to mimic the so-called "explicit σ-hole." In this work, we present an automated FF parameterization strategy based on a global optimization of both LJ and charge parameters with respect to accurate quantum mechanical data, purposely computed for the system under investigation. As a test case, we report on two homologue series, characterized either by weak or strong XBs, namely, the di-halogenated methanes and the mono-, di-, and tri-substituted acetonitriles, taking into consideration Cl, Br, and I substituents. The resulting quantum mechanically derived FFs are validated for each compound in the gas and in the condensed phase by comparing them to general purpose and specific FFs without virtual sites and to highly accurate reference quantum mechanical data. The results strongly support the adoption of the specific FFs with virtual sites, which overcome the other investigated models in representing both gas phase energetics and the structural patterns of the liquid phase structure related to the presence of XBs.

Coupled hydroxyl and ether functionalisation in EAN derivatives: the effect of hydrogen bond donor/acceptor groups on the structural heterogeneity studied with X-ray diffractions and fixed charge/polarizable simulations

We present a study by energy-dispersive X-ray diffraction of liquid 2-(2-hydroxyethoxy)ethan-1-ammonium nitrate, NH₃CH₂CH₂(OCH₂CH₂OH)⁺NO₃^- (22HHEAN). This ionic liquid is derived from the parent ethylammonium nitrate (EAN) with an ether link in the chain and a hydroxyl group in the terminal position. The absence of peaks at low-q values in the experimental diffraction curve indicates that the added polar groups and the high conformational isomerism of the cations alter strongly the nanosegregation of the parent EAN liquid. Aggregation between ionic species may involve hydrogen bonding between cations and anions and a variety of intermolecular hydrogen bonds between cations. Diffraction patterns are compared with the results of molecular dynamics simulations with two different force fields: the fixed point charge force field (GAFF) with different charge scaling protocols and the polarizable AMOEBA force field. Most point charge models lead to the appearance of a quite evident low q-peak which decreases gradually, when the percentage and type of the scaling (uniform vs. non-uniform) are increased. In the polarisable model and in the model where only anion charges are scaled to 20%, instead, the pre-peak is absent in agreement with our experiments.

Structure and dynamics of propylammonium nitrate-acetonitrile mixtures: An intricate multi-scale system probed with experimental and theoretical techniques

In this article, we report the study of structural and dynamical properties for a series of acetonitrile/propylammonium nitrate mixtures as a function of their composition. These systems display an unusual increase in intensity in their X-ray diffraction patterns in the low-q regime, and their ¹H-NMR diffusion-ordered NMR spectroscopy (DOSY) spectra display unusual diffusivities. However, the magnitude of both phenomena for mixtures of propylammonium nitrate is smaller than those observed for ethylammonium nitrate mixtures with the same cosolvent, suggesting that the cation alkyl tail plays an important role in these observations. The experimental X-ray scattering data are compared with the results of molecular dynamics simulations, including both ab initio studies used to interpret short-range interactions and classical simulations to describe longer range interactions. The higher level calculations highlight the presence of a strong hydrogen bond network within the ionic liquid, only slightly perturbed even at high acetonitrile concentration. These strong interactions lead to the symmetry breaking of the NO₃^- vibrations, with a splitting of about 88 cm^-1 in the ν₃ antisymmetric stretch. The classical force field simulations use a greater number of ion pairs, but are not capable of fully describing the longest range interactions, although they do successfully account for the observed concentration trend, and the analysis of the models confirms the nano-inhomogeneity of these kinds of samples.