Articles
Physical Chemistry Chemical Physics (14639084)26(40)pp. 26109-26128
In this study, we investigated the effect of DFT density functionals and dispersion correction on an imidazolium-based dicationic ionic liquid (DIL) using ab initio molecular dynamics simulations. To achieve this purpose, the electronic structures, as well as the structural and dynamical properties of [C3(mim)2][NTF2]2 DIL, were obtained using the BLYP and PBE functionals, both with and without D3-correction, and the results were compared with experimental values. Radial distribution functions and structure factors revealed that applying D3-correction increases the interaction between the anion and hydrogen atoms of the rings and side chains. The simulation of the studied DIL with the BLYP-D3 functional depicted lower structural heterogeneity compared to the other functionals. Analysis of Voronoi tessellation and linkage chain conformations showed a reduction in the aggregation of the linkage alkyl chains in the presence of D3-correction, which is more pronounced in the BLYP functional than in PBE. Additionally, it was observed that the probability of forming a hydrogen-bond network depends on both the type of used density functionals and applying dispersion correction. The results of dynamical properties, such as the self-diffusion coefficients, velocity autocorrelation function, and van Hove correlation function, as well as ion pair, ion cage, and hydrogen bond dynamics, indicated that applying D3-correction in both density functionals leads to an increase in the dynamics of the studied DIL. Additionally, the ratio of self-diffusion coefficients of the anion to the cation in the BLYP functional is closer to experimental values compared to the PBE functional. Furthermore, the electronic structure, including dipole moment distribution, and also infrared (IR) and power spectra were studied. Applying D3-correction and the type of density functionals have a significant effect on the dipole moment distribution of ions. Moreover, the results of IR and power spectra demonstrated that only in the BLYP functional, by applying D3-correction, the hydrogen bonding between the anion and the hydrogen atoms of the cation is strengthened at high wavenumbers. Thus, we conclude that applying D3 correction to both the BLYP and PBE density functionals improves the accuracy in describing the various properties of the studied system. Overall, the evaluation of different structural, dynamical, and vibrational properties of [C3(mim)2][NTF2]2 DIL suggests that the BLYP-D3 density functional may be the best choice among the studied density functionals. © 2024 The Royal Society of Chemistry.
Journal of Molecular Liquids (18733166)399
Ionanofluids (INFs) are a novel type of heat transfer fluids consisting of fine nanoparticles suspended in ionic liquids (ILs). They exhibit improved thermal properties, including enhanced thermal conductivity, non-volatility, and non-flammability, making them ideal candidates for various applications, especially solar concentrators. In this work, for the first time, molecular dynamics (MD) simulations are used to investigate the thermophysical, structural, and dynamical properties of a specific INF composed of the [HMIM][BF4] IL and SiC nanoparticles in the temperature range of 298.15–338.15 K and at different nanoparticle volume fractions, i.e. 1.10 %, 2.27 %, and 3.47 %. Radial distribution functions (RDFs) revealed a strong structural correlation between the fluorine atoms of anion and the HR atom of the cation ring and also with the carbon and silicon atoms of the nanoparticle. The results showed that with increasing the volume fraction of nanoparticles, ions tend to have more interactions with nanoparticles, and the interactions between the anions and cations decrease. As a result, the thermal conductivity, viscosity, and density of ions in the INF increases. Finally, we made some comparisons between the simulated thermal conductivities and viscosities of the studied INF with classical models. The results showed that, in some cases, the simulation results were even better than some models compared to the experimental values. The insights obtained from this study not only enhance our understanding of the structure and dynamic behavior of the INF, but also can contribute to the development of innovative materials and more efficient solar thermal systems. © 2024 Elsevier B.V.
Journal of Chemical Physics (10897690)160(15)
This work employs a combination of density functional theory-infrared (IR), conductor-like screening model for real solvents (COSMO-RS), and molecular dynamic (MD) methods to investigate the impact of hydroxyl functional groups on CO2 capture within dicationic ionic liquids (DILs). The COSMO-RS reveals that hydroxyl groups in DILs reduce the macroscopic solubility of CO2 but improve the selectivity of CO2 over CO, H2, and CH4 gases. Quantum methods in the gas phase and MD simulations in the liquid phase were conducted to delve deeper into the underlying mechanisms. The IR spectrum analysis confirms red shifts in CO2’s asymmetric stretching mode and blue shifts in the CR-HR bond of the dication, indicating CO2-DIL interactions and the weakening of the anion-cation interactions caused by the presence of CO2. The results show that the positioning of anions around hydroxyl groups and HR atoms in rings inhibits the proximity of CO2 molecules, causing the hydrogen atoms within methylene groups to accumulate CO2. van der Waals forces were found to dominate the interaction between ions and CO2. The addition of hydroxyl groups strengthens the electrostatic interactions and hydrogen bonds between dications and anions. The stronger interaction energy between ions in [C5(mim)2-(C2)2(OH)2][NTf2]2 limits the displacement of CO2 molecules within this DIL compared to [C5(mim)2-(C4)2][NTf2]2. Compared to [C5(mim)2-(C4)2][NTf2]2, [C5(mim)2-(C2)2(OH)2][NTf2]2 exhibits stronger ion-ion interactions, higher density, and reduced free volume, resulting in a reduction in CO2 capture. These results provide significant insights into the intermolecular interactions and vibrational properties of CO2 in DIL complexes, emphasizing their significance in developing efficient and sustainable strategies for CO2 capture. © 2024 Author(s).
Journal of Molecular Liquids (18733166)393
This study comprehensively investigates the effects of hydroxyl groups in side chains and cation symmetry on intra- and inter-molecular structures, dynamics, and thermodynamic properties of two hydroxyl-functionalized di-cationic ionic liquids (HFDILs). Using a combination of molecular dynamics (MD) simulations and quantum chemistry (QC) calculations, we compare the influence of these effects on various properties with corresponding non-hydroxyl DILs reported in previous works. The structure of HFDILs is explored using angle distribution, radial distribution, spatial distribution, and combined distribution functions. Also, the averaged noncovalent interaction (aNCI) analysis is used to characterize the weak noncovalent interactions in fluctuating environment of HFDILs. Results showed that the oxygen atoms of anion have the strongest interaction with the hydrogen atom of hydroxyl group, followed by the hydrogen atom at the top of the imidazolium ring. Domain analysis and Voronoi calculations are performed to investigate the number of polar and non-polar subunits, revealing that symmetric HFDILs with two hydroxyl groups in their side chains have greater structural heterogeneity than asymmetric ones with only one hydroxyl group due to the larger accumulation of non-polar subunits. To investigate dynamical heterogeneity, we calculate the non-Gaussian parameter, van Hove correlation, hydrogen bond, ion pair, ion cage, and reorientation dynamics. Results indicated that symmetric HFDILs have slower dynamics, greater dynamic heterogeneity, and more stability than asymmetric ones. Finally, the structure and strength of ion triplet interactions are examined using quantum mechanical methods, revealing that symmetric HFDILs are more stable than asymmetric ones. This study provides valuable information for understanding DILs and their potential applications in various fields. © 2023 Elsevier B.V.
In this study, we extended the optimized potentials for liquid simulation-ionic-liquid virtual site (OPLS-VSIL) force field (FF) to imidazolium-based dicationic ionic liquids (DILs) and evaluated the ability of different OPLS-based FFs (i.e., OPLS-2009IL, 0.8*OPLS-2009IL, and OPLS-VSIL) in predicting different properties of the studied DIL by comparing their results with ab initio molecular dynamics (AIMD) simulation and experimental results. To achieve this purpose, MD simulations with three different OPLS-based FFs as well as AIMD simulation were performed for [C3(mim)2][NTF2]2 DIL and its structural, dynamical, vibrational, and volumetric properties were analyzed. Structural properties of the studied DIL, i.e., radial distribution functions (RDFs), structure factor, and hydrogen-bond network, showed that compared to 0.8*OPLS-2009IL FF, there is a much better agreement between the results of both OPLS-2009IL and OPLS-VSIL FFs with the AIMD simulation. On the other hand, the results of dynamical properties, such as mean square displacements, van Hove correlation functions as well as hydrogen bond, ion pair, and ion cage dynamics, depicted that in both 0.8*OPLS-2009IL and OPLS-VSIL FFs, the dynamics of the system is almost similar, and compared to OPLS-2009IL FF, they have better agreements with experimental results where they exist. So, it can be seen that although reducing the total charge of studied DIL by 20% leads to an increase in the dynamics of the system, the type distribution of partial charges on each atom does not significantly affect the system’s dynamics. The calculated infrared (IR) and power spectra showed that the vibrational features of studied DIL in three OPLS-based FFs are mostly the same and reducing total charge and different type distribution of partial charges have no significant effect on the studied system. Furthermore, in volumetric properties, OPLS-VSIL FF shows somehow better agreement with experimental results. Overall, the evaluation of different structural, dynamical, vibrational, and volumetric properties of [C3(mim)2][NTF2]2 DIL shows that the OPLS-VSIL FF may be the best choice among the different studied OPLS FFs. © 2023 Author(s).
