Tuning Solvent Miscibility: A Fundamental Assessment on the Example of Induced Methanol/n-Dodecane Phase Separation
- Roberto MacchieraldoRoberto MacchieraldoMulliken Center for Theoretical Chemistry, University of Bonn, Beringstrasse 4+6, D-53115 Bonn, GermanyMore by Roberto Macchieraldo
- ,
- Sascha GehrkeSascha GehrkeMulliken Center for Theoretical Chemistry, University of Bonn, Beringstrasse 4+6, D-53115 Bonn, GermanyMax Planck Institute for Chemical Energy Conversion, Stiftstrasse 34-36, D-45413 Mülheim an der Ruhr, GermanyMore by Sascha Gehrke
- ,
- Nagaphani K. BatchuNagaphani K. BatchuDepartment of Chemistry, KU Leuven, Celestijnenlaan 200F, bus 2404, B-3001 Heverlee, BelgiumMore by Nagaphani K. Batchu
- ,
- Barbara Kirchner*Barbara Kirchner*E-mail: [email protected] (B.K.).Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstrasse 4+6, D-53115 Bonn, GermanyMore by Barbara Kirchner
- , and
- Koen Binnemans*Koen Binnemans*E-mail: [email protected] (K.B.).Department of Chemistry, KU Leuven, Celestijnenlaan 200F, bus 2404, B-3001 Heverlee, BelgiumMore by Koen Binnemans
Abstract

In this work, we assess the fundamental aspects of mutual miscibility of solvents by studying the mixing of two potential candidates, methanol and n-dodecane, for nonaqueous solvent extraction. To do so, 1H NMR spectroscopy and molecular dynamics simulations are used jointly. The NMR spectra show that good phase separation can be obtained by adding LiCl and that the addition of a popular extractant (tri-n-butyl phosphate) yields the opposite effect. It is also demonstrated that in a specific case the poor phase separation is not due to the migration of n-dodecane into the more polar phase, but due to the transfer of the extractant into it, which is especially relevant when considering industrial applications of solvent extraction. With the aid of molecular dynamics simulations, explanations of this behavior are given. Specifically, an increase of all hydrogen-bond lifetimes is found to be consequent to the addition of LiCl which implies an indirect influence on the methanol liquid structure, by favoring a stronger hydrogen-bond network. Therefore, we found that better phase separation is not directly due to the presence of LiCl, but due to the “hardening” of the hydrogen-bond network.
Introduction
Figure 1

Figure 1. Schematic representation of the procedure described in this work.
Methods
Chemicals and Instrumentation
Experimental Details
Computational Details
Results
Volume Analysis
MeOH + LiCl | DD | DD + TBP |
---|---|---|
neat MeOH | 1:1.1 | 1:2.3 |
0.5 | 1:1 | 1:2 |
1 | 1:1 | 1:1.9 |
2 | 1:1 | 1:1.9 |
3 | 1:1 | 1:1.9 |
Data are referred to the volume ratio between the top phase (apolar) and the bottom phase (polar) obtained after mixing and phase separation.
DD + TBP | MeOH + LiCl |
---|---|
neat DD | 1:1 |
0.25 | 1:1.4 |
0.5 | 1:1.8 |
1 | 1:1.9 |
2 | 1:4.5 |
Data reported is referred to the volume ratio between the top phase (apolar) and the bottom phase (polar) after mixing and phase separation.
1H NMR Spectra
Systems Simulated via Molecular Dynamics
system | MD | MDT | MLD | MLDT |
---|---|---|---|---|
MeOH | 3000 | 3000 | 3000 | 3000 |
DD | 32 | 100 | 2 | 10 |
TBP | 120 | 120 | ||
LiCl | 361 | 361 |
Data refer to the number of molecules in each system.
Solvent Partitioning
system | MeOH | MeOH + LiCl | LiCl | DD | DD + TBP | TBP |
---|---|---|---|---|---|---|
MD | 1.0 | 12.9 (32) | ||||
MLD | 1.0 | 1.0 | 291.8 (361) | 1.5 (2) | ||
MDT | 2.7 | 9.3 (100) | 3.9 | 6.4 (120) | ||
MLDT | 1.2 | 1.2 | 274.0 (361) | 7.0 (10) | 3.4 | 3.6 (120) |
The values refer to the average number of domains of the subsets into the systems (within parentheses is the number of molecules for that component. For MeOH, the number of molecules is always 3000).
Figure 2

Figure 2. Snapshot taken from the MLDT system. All atoms are represented by their van der Waals radii. Color scheme: yellow for LiCl, red for MeOH, green for TBP, and blue for DD.
Mobility
Figure 3

Figure 3. Mean square displacements of the bottom phase components. Top left: MeOH; top right: DD; bottom left: Li+ and Cl–; bottom right: DD and TBP.
component | MeOH | Li+ | Cl– | DD | TBP |
---|---|---|---|---|---|
MD | 2.53 | 1.74 | |||
MLD | 0.91 | 0.38 | 0.48 | 0.47 | |
MDT | 2.27 | 0.71 | 0.71 | ||
MLDT | 0.82 | 0.36 | 0.39 | 0.69 | 0.39 |
Structural Properties
Figure 4

Figure 4. RDFs of the interactions of MeOH molecules in the bottom phases.
Figure 5

Figure 5. Example of the liquid structure of MeOH in the presence of LiCl taken from the MLDT system. LiCl is represented with van der Waals radii while MeOH is displayed with the ball-and-stick models. Color scheme: brown for Li, green for Cl, orange for C, red for O, and white for H.
Figure 6

Figure 6. Intramolecular RDFs of the distance between the terminal carbons of DD molecules.
Hydrogen-Bond (HB) Network
MeOH | Cl– | TBP | |
---|---|---|---|
MD | 13.8 | ||
MLD | 28.8 | 15.5 | |
MDT | 16.5 | 10.2 | |
MLDT | 32.6 | 18.7 | 20.1 |
Conclusions
Supporting Information
The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpcb.9b00839.
1H NMR spectra; radial distribution function; RDFs of the interactions of MeOH(O) molecules (PDF)
Terms & Conditions
Most electronic Supporting Information files are available without a subscription to ACS Web Editions. Such files may be downloaded by article for research use (if there is a public use license linked to the relevant article, that license may permit other uses). Permission may be obtained from ACS for other uses through requests via the RightsLink permission system: http://pubs.acs.org/page/copyright/permissions.html.
Acknowledgments
B.K., K.B., and R.M. acknowledge the MSCA ETN project SOCRATES (http://etn-socrates.eu/), which has received funding from the European Union’s EU Framework Program for Research and Innovation Horizon 2020 under Grant Agreement No. 721385. S.G. was supported by the International Max Planck Research School for Reactive Structure Analyses for Chemical Reactions. K.B. and N.K.B. acknowledge funding from the European Research Council (ERC) under the European Union’s Horizon 2020 Research and Innovation Program: Grant Agreement 694078—Solvometallurgy for critical metals (SOLCRIMET).
References
This article references 50 other publications.
- 1Hiskey, B. Kirk-Othmer Encyclopedia of Chemical Technology; Wiley Blackwell: Hoboken, New Jersey, 2000.Google ScholarThere is no corresponding record for this reference.
- 2Anderson, C. G.; Dunne, R. C.; Uhrie, J. L. Mineral Processing and Extractive Metallurgy: 100 years of Innovation; Society for Mining, Metallurgy, and Exploration: Englewood, 2014.Google ScholarThere is no corresponding record for this reference.
- 3Hayes, P. C. Process Principles in Minerals and Materials Production; Hayes Publishing Co: Sherwood, 2003.Google ScholarThere is no corresponding record for this reference.
- 4Geissdoerfer, M.; Savaget, P.; Bocken, N. M.; Hultink, E. J. The circular economy - a new sustainability paradigm?. J. Clean. Prod. 2017, 143, 757– 768, DOI: 10.1016/j.jclepro.2016.12.048Google ScholarThere is no corresponding record for this reference.
- 5Binnemans, K.; Jones, P. T. Solvometallurgy: an emerging branch of extractive netallurgy. J. Sustainable Metall. 2017, 3, 570– 600, DOI: 10.1007/s40831-017-0128-2Google ScholarThere is no corresponding record for this reference.
- 6Abbott, A. P.; Boothby, D.; Capper, G.; Davies, D. L.; Rasheed, R. K. Deep eutectic solvents formed between choline chloride and carboxylic acids: versatile alternatives to ionic liquids. J. Am. Chem. Soc. 2004, 126, 9142– 9147, DOI: 10.1021/ja048266jGoogle Scholar6https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2cXlt1Gns70%253D&md5=5eba50624615fcd43af5afb2b11c0808Deep Eutectic Solvents Formed between Choline Chloride and Carboxylic Acids: Versatile Alternatives to Ionic LiquidsAbbott, Andrew P.; Boothby, David; Capper, Glen; Davies, David L.; Rasheed, Raymond K.Journal of the American Chemical Society (2004), 126 (29), 9142-9147CODEN: JACSAT; ISSN:0002-7863. (American Chemical Society)Deep Eutectic Solvents (DES) can be formed between a variety of quaternary ammonium salts and carboxylic acids. The phys. properties are significantly affected by the structure of the carboxylic acid but the phase behavior of the mixts. can be simply modeled by taking account of the mole fraction of carboxylic acid in the mixt. The phys. properties such as viscosity, cond., and surface tension of these DES are similar to ambient temp. ionic liqs. and insight into the cause of these properties is gained using hole-theory. It is shown that the cond. and viscosity of these liqs. is controlled by ion mobility and the availability of voids of suitable dimensions, and this is consistent with the fluidity of other ionic liqs. and molten salts. The DES are also shown to be good solvents for metal oxides, which could have potential application for metal extn.
- 7Rogers, R. D.; Seddon, K. R. Ionic liquids-solvents of the future?. Science 2003, 302, 792– 793, DOI: 10.1126/science.1090313Google Scholar7https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A280%3ADC%252BD3srhsV2jug%253D%253D&md5=0d242bb94d4f3a63249d106c6dab9e79Chemistry. Ionic liquids--solvents of the future?Rogers Robin D; Seddon Kenneth RScience (New York, N.Y.) (2003), 302 (5646), 792-3 ISSN:.There is no expanded citation for this reference.
- 8Gutowski, K. E.; Broker, G. A.; Willauer, H. D.; Huddleston, J. G.; Swatloski, R. P.; Holbrey, J. D.; Rogers, R. D. Controlling the aqueous miscibility of ionic liquids: aqueous biphasic systems of Weter-miscible ionic liquids and water-structuring salts for recycle, metathesis, and separations. J. Am. Chem. Soc. 2003, 125, 6632– 6633, DOI: 10.1021/ja0351802Google Scholar8https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3sXjs12ks7k%253D&md5=c46b031e62907c751321aeab5c7cee1fControlling the Aqueous Miscibility of Ionic Liquids: Aqueous Biphasic Systems of Water-Miscible Ionic Liquids and Water-Structuring Salts for Recycle, Metathesis, and SeparationsGutowski, Keith E.; Broker, Grant A.; Willauer, Heather D.; Huddleston, Jonathan G.; Swatloski, Richard P.; Holbrey, John D.; Rogers, Robin D.Journal of the American Chemical Society (2003), 125 (22), 6632-6633CODEN: JACSAT; ISSN:0002-7863. (American Chemical Society)Hydrophilic ionic liqs. can be salted-out and concd. from aq. soln. upon addn. of salts forming aq. biphasic systems as illustrated by the phase behavior of mixts. of 1-butyl-3-methylimidazolium chloride and K3PO4.
- 9Earle, M. J.; Esperança, J. M.; Gilea, M. A.; Lopes, J. N. C.; Rebelo, L. P.; Magee, J. W.; Seddon, K. R.; Widegren, J. A. The distillation and volatility of ionic liquids. Nature 2006, 439, 831– 834, DOI: 10.1038/nature04451Google Scholar9https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD28XhsVSkurk%253D&md5=9e8d177c2100d4ecc45c62f3ca92c357The distillation and volatility of ionic liquidsEarle, Martyn J.; Esperanca, Jose M. S. S.; Gilea, Manuela A.; Canongia Lopes, Jose N.; Rebelo, Luis P. N.; Magee, Joseph W.; Seddon, Kenneth R.; Widegren, Jason A.Nature (London, United Kingdom) (2006), 439 (7078), 831-834CODEN: NATUAS; ISSN:0028-0836. (Nature Publishing Group)It is widely believed that a defining characteristic of ionic liqs. (or low-temp. molten salts) is that they exert no measurable vapor pressure, and hence cannot be distd. Here we demonstrate that this is unfounded, and that many ionic liqs. can be distd. at low pressure without decompn. Ionic liqs. represent matter solely composed of ions, and so are perceived as non-volatile substances. During the last decade, interest in the field of ionic liqs. has burgeoned, producing a wealth of intellectual and technol. challenges and opportunities for the prodn. of new chem. and extractive processes, fuel cells and batteries, and new composite materials. Much of this potential is underpinned by their presumed involatility. This characteristic, however, can severely restrict the attainability of high purity levels for ionic liqs. (when they contain poorly volatile components) in recycling schemes, as well as excluding their use in gas-phase processes. We anticipate that our demonstration that some selected families of commonly used aprotic ionic liqs. can be distd. at 200-300 °C and low pressure, with concomitant recovery of significant amts. of pure substance, will permit these currently excluded applications to be realized.
- 10Fukumoto, K.; Yoshizawa, M.; Ohno, H. Room temperature ionic liquids from 20 Natural amino acids. J. Am. Chem. Soc. 2005, 127, 2398– 2399, DOI: 10.1021/ja043451iGoogle Scholar10https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2MXptlSgsA%253D%253D&md5=b8ffaba7d45996888c6c26708d46cfa9Room Temperature Ionic Liquids from 20 Natural Amino AcidsFukumoto, Kenta; Yoshizawa, Masahiro; Ohno, HiroyukiJournal of the American Chemical Society (2005), 127 (8), 2398-2399CODEN: JACSAT; ISSN:0002-7863. (American Chemical Society)The authors have succeeded in synthesizing room temp. ionic liqs. from 20 natural amino acids (AA). These amino acid ionic liqs. are H-AA-O-·(1-ethyl-3-methylimidazolium)+ salts, which dissolved native amino acids despite water-free conditions. Furthermore, these ionic liqs. are sol. in various org. solvents, such as chloroform. Effects of acidity, hydrogen bonding ability, and steric factors on the properties of these ionic liqs. were analyzed as the function of side groups.
- 11Wang, Y.; Voth, G. A. Unique spatial heterogeneity in ionic liquids. J. Am. Chem. Soc. 2005, 127, 12192– 12193, DOI: 10.1021/ja053796gGoogle Scholar11https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2MXnslOrtbc%253D&md5=139435d4120f25a320b1cae5af1e0fabUnique Spatial Heterogeneity in Ionic LiquidsWang, Yanting; Voth, Gregory A.Journal of the American Chemical Society (2005), 127 (35), 12192-12193CODEN: JACSAT; ISSN:0002-7863. (American Chemical Society)A multiscale coarse-graining model for ionic liqs. has been extended to investigate the unique aggregation of cations in ionic liqs. through computer simulation. It has been found that, with sufficiently long side chains, the tail groups of cations aggregate to form spatially heterogeneous domains, while headgroups of the cations and the anions distribute as uniformly as possible. This is understood as the result of competition between the charged electrostatic interactions between headgroups and anions and the collective short-range interactions between the neutral tail groups. This aggregation can help to explain a no. of exptl. obsd. phys. phenomena in ionic liqs.
- 12Scurto, A. M.; Aki, S. N. V. K.; Brennecke, J. F. CO2 as a separation switch for ionic liquid/organic mixtures. J. Am. Chem. Soc. 2002, 124, 10276– 10277, DOI: 10.1021/ja0268682Google Scholar12https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD38XlvFGlt7w%253D&md5=a7afd1d2406e4de0c70dbbfc7c3dba23CO2 as a Separation Switch for Ionic Liquid/Organic MixturesScurto, Aaron M.; Aki, Sudhir N. V. K.; Brennecke, Joan F.Journal of the American Chemical Society (2002), 124 (35), 10276-10277CODEN: JACSAT; ISSN:0002-7863. (American Chemical Society)A novel technique to sep. ionic liqs. from org. compds. is introduced which uses carbon dioxide to induce the formation of an ionic liq.-rich phase and an org.-rich liq. phase in mixts. of methanol and 3-butyl-1-methyl-imidazolium hexafluorophosphate ([C4mim][PF6]). If the temp. is above the crit. temp. of CO2 then the methanol-rich phase can become completely miscible with the CO2-rich phase, and this new phase is completely ionic liq.-free. Since CO2 is nonpolar, it is not equipped to solvate ions. As the CO2 dissolves in the methanol/[C4mim][PF6] mixt., the solvent power of the CO2-expanded liq. is significantly reduced, inducing the formation of the second liq. phase that is rich in ionic liq. This presents a new way to recover products from ionic liq. mixts. and purify org. phases that have been contaminated with ionic liq. Moreover, these results have important implications for reactions done in CO2/ionic liq. biphasic mixts.
- 13Sastri, V. Selective leaching of metals from ores with organic solvents and chelating agents. J. Sci. Ind. Res. 1975, 34, 663– 665Google Scholar13https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaE28Xks1Kiur8%253D&md5=b045526d899e5f5fb1ee928d8dd423a8Selective leaching of metals from ores with organic solvents and chelating agentsSastri, V. S.Journal of Scientific & Industrial Research (1975), 34 (12), 663-5CODEN: JSIRAC; ISSN:0022-4456.A review with 16 refs. is given on the use of org. solvents and chelating agents in the selective leaching of metals from ores.
- 14Raghavan, S.; Fuerstenau, D. A Lyometallurgical Process for Leaching Copper from Chrysocolla. In Hydrometallurgy Fundamentals, Technology and Innovations , Proceedings of the Milton E., Wordsworth (IV) Int. Symposium on Hydrometallurgy, 1993; pp 283– 297.Google ScholarThere is no corresponding record for this reference.
- 15Chizhevskaya, S. V.; Chekmarev, A. M.; Klimenko, O. M.; Povetkina, M. V.; Sinegribova, O. A.; Cox, M. In Nontraditional Methods of Treating High-Silicon Ores Containing Rare Elements , Hydrometallurgy ’94: Papers presented at the international symposium ‘Hydrometallurgy ’94’ organized by the Institution of Mining and Metallurgy and the Society of Chemical Industry, and held in Cambridge, England, from 11 to 15 July, 1994; pp 219– 228.Google ScholarThere is no corresponding record for this reference.
- 16Hiroshi, T.; Kazushige, T.; Yoshimi, O. Direct Dissolution of Water-Insoluble Uranium Compounds by Contact with Neutral Organic Solvents Pretreated with Nitric Acid. U.S. Patent US3288568, 1966.Google ScholarThere is no corresponding record for this reference.
- 17Batchu, N. K.; Vander Hoogerstraete, T.; Banerjee, D.; Binnemans, K. Non-aqueous solvent extraction of rare-earth nitrates from ethylene glycol to n-dodecane by Cyanex 923. Sep. Purif. Technol. 2017, 174, 544– 553, DOI: 10.1016/j.seppur.2016.10.039Google Scholar17https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28XhslyjsL%252FP&md5=6f6ebbff4bbd5464d240c22512f0cf07Non-aqueous solvent extraction of rare-earth nitrates from ethylene glycol to n-dodecane by Cyanex 923Batchu, Nagaphani Kumar; Vander Hoogerstraete, Tom; Banerjee, Dipanjan; Binnemans, KoenSeparation and Purification Technology (2017), 174 (), 544-553CODEN: SPUTFP; ISSN:1383-5866. (Elsevier B.V.)A solvent extn. process comprising two immiscible org. phases was developed for the extn. of rare earths. The more polar org. phase was ethylene glycol with dissolved rare-earth nitrate salts and lithium nitrate, while the less polar phase was a soln. of the neutral extractant Cyanex 923 dissolved in n-dodecane. The solvent extn. mechanism was detd. by slope anal. and the main species in the org. phase was identified by Extended x-ray Absorption Fine Structure (EXAFS) studies. The extn. from the ethylene glycol soln. was compared with extn. from an aq. feed soln. When compared to aq. feed solns., the light rare-earth elements (LREEs) are less efficiently extd. and the heavy rare-earth elements (HREEs) more efficiently extd. from an ethylene glycol feed soln., resulting into the easy sepn. of HREEs from LREEs. The sepn. factors between neighboring elements are higher for this non-aq. solvent extn. process than for extn. from an aq. feed soln.
- 18Li, Z.; Li, X.; Raiguel, S.; Binnemans, K. Separation of transition metals from rare earths by non-aqueous solvent extraction from ethylene glycol solutions using Aliquat 336. Sep. Purif. Technol. 2018, 201, 318– 326, DOI: 10.1016/j.seppur.2018.03.022Google Scholar18https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXlslKgs7k%253D&md5=b1a5cd1a9bc8ea4fa342fa21a0c64f00Separation of transition metals from rare earths by nonaqueous solvent extraction from ethylene glycol solutions using Aliquat 336Li, Zheng; Li, Xiaohua; Raiguel, Stijn; Binnemans, KoenSeparation and Purification Technology (2018), 201 (), 318-326CODEN: SPUTFP; ISSN:1383-5866. (Elsevier B.V.)Solvent extn. is a widely used sepn. technique in extractive metallurgy. A conventional solvent extn. system consists of an aq. phase and an immiscible org. phase. Replacement of H2O by a polar org. solvent can lead to superior metal sepns. Co(II) and Sm(III) chlorides dissolved in H2O and ethylene glycol (EG), resp., with LiCl as chloride source, were extd. by Aliquat 336 dild. in toluene. Both Co and Sm were extd. from the aq. solns., but Co was extd. more efficiently from the ethylene glycol solns. than from the aq. solns., whereas Sm was not extd. at all from the ethylene glycol solns. As a result, Co and Sm could be sepd. completely in a single extn. step from ethylene glycol solns. The mechanisms of Co extn. by Aliquat 336 from the ethylene glycol and aq. solns. are similar, as validated by slope anal. and UV-visible absorption spectroscopy. Sm was extd. from the aq. soln. through the salting-out effect of LiCl. LiCl has a much lower salting-out effect for Sm in ethylene glycol than in H2O due to the lower dielec. const. of ethylene glycol and the lower soly. of LiCl in ethylene glycol. Consequently, Sm is not salted out from ethylene glycol, leading to a very efficient sepn. of Co and Sm. This sepn. effect can also be applied to the sepn. of other transition metal and rare-earth metal pairs, including Fe/Nd and Zn/Eu.
- 19Misawa, M.; Yoshida, K.; Maruyama, K.; Munemura, H.; Hosokawa, Y. Salt-induced phase separation in aqueous solution. J. Phys. Chem. Solids 1999, 60, 1301– 1306, DOI: 10.1016/S0022-3697(99)00108-0Google Scholar19https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK1MXls1equrY%253D&md5=16a5c9ae3269dc0a0c8cf3ca669cca89Salt-induced phase separation in aqueous solutionMisawa, M.; Yoshida, K.; Maruyama, K.; Munemura, H.; Hosokawa, Y.Journal of Physics and Chemistry of Solids (1999), 60 (8-9), 1301-1306CODEN: JPCSAW; ISSN:0022-3697. (Elsevier Science Ltd.)KCl-induced phase sepn. phenomenon in 1-propanol aq. soln. has been studied by means of small angle (SANS) and quasi-elastic neutron scattering measurements. The significant increase of SANS intensity by the addn. of KCl is attributed to the aggregate of clusters, of which the fractal nature of dimension 1.8-1.9 is found. The motion of water mols. is much suppressed in the 1-propanol-water soln., while it is slightly recovered by the addn. of KCl. The aggregate of clusters and liq.-liq. phase sepn. are discussed in relation to the motion of water mols.
- 20Batchu, N. K.; Vander Hoogerstraete, T.; Banerjee, D.; Binnemans, K. Separation of rare-earth ions from ethylene glycol (+LiCl) solutions by non-aqueous solvent extraction with Cyanex 923. RSC Adv. 2017, 7, 45351– 45362, DOI: 10.1039/C7RA09144CGoogle Scholar20https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXhsFCktL7I&md5=4e3bd3a8ac008e12abeb98934d7b2d9fSeparation of rare-earth ions from ethylene glycol (+LiCl) solutions by non-aqueous solvent extraction with Cyanex 923Batchu, Nagaphani Kumar; Vander Hoogerstraete, Tom; Banerjee, Dipanjan; Binnemans, KoenRSC Advances (2017), 7 (72), 45351-45362CODEN: RSCACL; ISSN:2046-2069. (Royal Society of Chemistry)The sepn. of a mixt. of rare earths by non-aq. solvent extn. with two immiscible org. phases has been studied. The more polar org. phase was ethylene glycol with dissolved lithium chloride and the less polar org. phase was the extractant dild. in n-dodecane. Cyanex 923 was found to be the most performant extractant amongst the investigated acidic, basic and solvating extractants: Cyanex 272, Cyphos IL 101, Aliquat 336, bis(2-ethylhexyl)amine, trioctylphosphine oxide (TOPO) and Cyanex 923. The replacement of the aq. chloride feed solns. by non-aq. ethylene glycol feed solns. had a profound effect on the distribution ratios and sepn. factors. The sepn. factors for extn. of pairs of rare earths from aq. chloride solns. by Cyanex 923 are too low to be of practical use. On the contrary, a mixt. of rare earths can be sepd. conveniently in four different groups by extn. with Cyanex 923 from ethylene glycol (+LiCl) solns. The influence of several parameters such as the chloride concn., the type of chloride salt, the addn. of other polar solvents to the ethylene glycol phase, the addn. of second extractant to the less polar org. phase, and the addn. of complexing agents to the ethylene glycol phase has been studied. The extn. mechanism for extn. of ytterbium(III) was studied by slope anal. expts. The ytterbium(III) species in the ethylene glycol phase and the extd. species in the n-dodecane phase were detd. by EXAFS. Furthermore, a conceptual flow sheet for the fractionation of rare earths from an ethylene glycol (+LiCl) feed soln. into different groups by extn. with Cyanex 923 has been proposed. The new extn. system is useful for extn. of scandium and for sepn. of scandium from the other REEs.
- 21Casás, L. M.; Touriño, A.; Orge, B.; Marino, G.; Iglesias, M.; Tojo, J. Thermophysical properties of acetone or methanol + n–alkane (C9 to C12) mixtures. J. Chem. Eng. Data 2002, 47, 887– 893, DOI: 10.1021/je0103059Google Scholar21https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD38XktV2gsLw%253D&md5=dc84bff396c2f3d1c9ef3216dbafb093Thermophysical Properties of Acetone or Methanol + n-Alkane (C9 to C12) MixturesCasas, Lidia M.; Tourino, Antonio; Orge, Beatriz; Marino, Gonzalo; Iglesias, Miguel; Tojo, JoseJournal of Chemical and Engineering Data (2002), 47 (4), 887-893CODEN: JCEAAX; ISSN:0021-9568. (American Chemical Society)Liq.-liq. equil. for methanol + n-alkanes (C9 to C12) were detd. in the temp. range from 278.15 to 308.15 K, using a visual static method. Densities, refractive indexes, and speeds of sound were also measured at 288.15, 298.15, and 308.15 K for acetone + n-alkane mixts. and at 298.15 and 308.15 K for methanol + n-alkane mixts. to study the temp. influence on mixing phenomena. The results were used to model liq.-liq. coexistence and derived properties trend. The calcd. values based on the UNIQUAC equation were found to be similar to those based on the NRTL model. The Prigogine-Flory-Paterson theory was applied to est. the excess molar volumes and the UNIFAC group contribution model for phase equil.
- 22Casás, L. M.; Orge, B.; Díaz, C.; Tojo, J. Liquid-liquid equilibria for mixtures of methyl acetate+methanol+n-alkane (C10–C12) at several temperatures and 1 atm. J. Chem. Thermodyn. 2004, 36, 237– 243, DOI: 10.1016/j.jct.2003.11.009Google Scholar22https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2cXhtVensL0%253D&md5=98c0ac16ca68e3c3d525015980b6cdebLiquid-liquid equilibria for mixtures of {methyl acetate + methanol + n-alkane (C10-C12)} at several temperatures and 1 atmCasas, Lidia M.; Orge, Beatriz; Diaz, Concepcion; Tojo, JoseJournal of Chemical Thermodynamics (2004), 36 (3), 237-243CODEN: JCTDAF; ISSN:0021-9614. (Elsevier Science Ltd.)The liq.-liq. equil. of the mixts. {methyl acetate + methanol + n-alkane (C10-C12)} at atm. pressure in the temp. range (278.15 to 308.15) K are presented. The equil. compns. were measured by gas chromatog. from samples splitted isothermally into a glass-stirred device and phase diagrams are reported at each temp. The exptl. detd. liq.-liq. equil. were satisfactorily correlated by the NRTL and UNIQUAC equations. A poor description has been obtained using the UNIFAC method and their modifications.
- 23Touriño, A.; Casás, L.; Marino, G.; Iglesias, M.; Orge, B.; Tojo, J. Liquid phase behaviour and thermodynamics of acetone+methanol+n-alkane (C9–C12) mixtures. Fluid Phase Equilib. 2003, 206, 61– 85, DOI: 10.1016/S0378-3812(02)00306-0Google Scholar23https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3sXjtl2qu70%253D&md5=0ff46624a178726290d28c9d57a0666fLiquid phase behaviour and thermodynamics of acetone+methanol+n-alkane (C9-C12) mixturesTourino, A.; Casas, L. M.; Marino, G.; Iglesias, M.; Orge, B.; Tojo, J.Fluid Phase Equilibria (2003), 206 (1-2), 61-85CODEN: FPEQDT; ISSN:0378-3812. (Elsevier Science B.V.)This paper reports the results of a new exptl. study of two-liq. phases equil. as a function of temp. and one-liq. phase thermodn. magnitudes (densities, refractive indexes and speeds of sound), covering the compn. diagrams for the mixts. acetone + methanol + n-alkane (C9-C12). The obtained coexistence curves were very asym. with respect to equimolar alc. + n-alkane compn., such effect increasing with the length of the aliph. chain and temp. A comparative anal. was performed by application of different methods to predict exptl. liq.-liq. equil. (LLE) behavior and thermodn. of these ternary mixts. The obtained exptl. results let us know the potential role of the linear aliph. alkanes as sepn. agents for heterogeneous modified distn. of the azeotrope acetone + methanol.
- 24Martínez, L.; Andrade, R.; Birgin, E. G.; Martínez, J. M. A package for building initial configurations for molecular dynamics simulations. J. Comput. Chem. 2009, 30, 2157– 2164, DOI: 10.1002/jcc.21224Google Scholar24https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1MXptleqsb8%253D&md5=2a76255c873b866a26540f7e84496272PACKMOL: A package for building initial configurations for molecular dynamics simulationsMartinez, L.; Andrade, R.; Birgin, E. G.; Martinez, J. M.Journal of Computational Chemistry (2009), 30 (13), 2157-2164CODEN: JCCHDD; ISSN:0192-8651. (John Wiley & Sons, Inc.)Adequate initial configurations for mol. dynamics simulations consist of arrangements of mols. distributed in space in such a way to approx. represent the system's overall structure. In order that the simulations are not disrupted by large van der Waals repulsive interactions, atoms from different mols. must keep safe pairwise distances. Obtaining such a mol. arrangement can be considered a packing problem: Each type mol. must satisfy spatial constraints related to the geometry of the system, and the distance between atoms of different mols. must be greater than some specified tolerance. We have developed a code able to pack millions of atoms, grouped in arbitrarily complex mols., inside a variety of three-dimensional regions. The regions may be intersections of spheres, ellipses, cylinders, planes, or boxes. The user must provide only the structure of one mol. of each type and the geometrical constraints that each type of mol. must satisfy. Building complex mixts., interfaces, solvating biomols. in water, other solvents, or mixts. of solvents, is straightforward. In addn., different atoms belonging to the same mol. may also be restricted to different spatial regions, in such a way that more ordered mol. arrangements can be built, as micelles, lipid double-layers, etc. The packing time for state-of-the-art mol. dynamics systems varies from a few seconds to a few minutes in a personal computer. The input files are simple and currently compatible with PDB, Tinker, Molden, or Moldy coordinate files. The package is distributed as free software and can be downloaded from . © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009.
- 25Plimpton, S. Fast parallel algorithms for short-range molecular dynamics. J. Comput. Phys. 1995, 117, 1– 19, DOI: 10.1006/jcph.1995.1039Google Scholar25https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK2MXlt1ejs7Y%253D&md5=715052332237e4cf8501f8fb81234017Fast parallel algorithms for short-range molecular dynamicsPlimpton, SteveJournal of Computational Physics (1995), 117 (1), 1-19CODEN: JCTPAH; ISSN:0021-9991.Three parallel algorithms for classical mol. dynamics are presented. The first assigns each processor a fixed subset of atoms; the second assigns each a fixed subset of inter-at. forces to compute; the third assigns each a fixed spatial region. The algorithms are suitable for mol. dynamics models which can be difficult to parallelize efficiently - those with short-range forces where the neighbors of each atom change rapidly. They can be implemented on any distributed-memory parallel machine which allows for message-passing of data between independently executing processors. The algorithms are tested on a std. Lennard-Jones benchmark problem for system sizes ranging from 500 to 100,000,000 atoms on several parallel supercomputers - the nCUBE 2, Intel iPSC/860 and Paragon, and Cray T3D. Comparing the results to the fastest reported vectorized Cray Y-MP and C90 algorithm shows that the current generation of parallel machines is competitive with conventional vector supercomputers even for small problems. For large problems, the spatial algorithm achieves parallel efficiencies of 90% and a 1840-node Intel Paragon performs up to 165 faster than a single Cray C90 processor. Trade-offs between the three algorithms and guidelines for adapting them to more complex mol. dynamics simulations are also discussed.
- 26Jorgensen, W. L.; Maxwell, D. S.; Tirado-Rives, J. Development and testing of the OPLS all-atom force field on conformational energetics and properties of organic liquids. J. Am. Chem. Soc. 1996, 118, 11225– 11236, DOI: 10.1021/ja9621760Google Scholar26https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK28XmtlOitrs%253D&md5=fef2924a69421881390282aa309ae91bDevelopment and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic LiquidsJorgensen, William L.; Maxwell, David S.; Tirado-Rives, JulianJournal of the American Chemical Society (1996), 118 (45), 11225-11236CODEN: JACSAT; ISSN:0002-7863. (American Chemical Society)The parametrization and testing of the OPLS all-atom force field for org. mols. and peptides are described. Parameters for both torsional and nonbonded energetics have been derived, while the bond stretching and angle bending parameters have been adopted mostly from the AMBER all-atom force field. The torsional parameters were detd. by fitting to rotational energy profiles obtained from ab initio MO calcns. at the RHF/6-31G*//RHF/6-31G* level for more than 50 org. mols. and ions. The quality of the fits was high with av. errors for conformational energies of less than 0.2 kcal/mol. The force-field results for mol. structures are also demonstrated to closely match the ab initio predictions. The nonbonded parameters were developed in conjunction with Monte Carlo statistical mechanics simulations by computing thermodn. and structural properties for 34 pure org. liqs. including alkanes, alkenes, alcs., ethers, acetals, thiols, sulfides, disulfides, aldehydes, ketones, and amides. Av. errors in comparison with exptl. data are 2% for heats of vaporization and densities. The Monte Carlo simulations included sampling all internal and intermol. degrees of freedom. It is found that such non-polar and monofunctional systems do not show significant condensed-phase effects on internal energies in going from the gas phase to the pure liqs.
- 27Das, A.; Sahu, P.; Ali, S. M. Molecular dynamics simulation for the calibration of the OPLS force field using DFT derived partial charges for the screening of alkyl phosphate ligands by studying structure, dynamics, and thermodynamics. J. Chem. Eng. Data 2017, 62, 2280– 2295, DOI: 10.1021/acs.jced.7b00096Google Scholar27https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXht1Slt73E&md5=2ec0cfaf7877976bc614ed86dea8dd79Molecular Dynamics Simulation for the Calibration of the OPLS Force Field Using DFT Derived Partial Charges for the Screening of Alkyl Phosphate Ligands by Studying Structure, Dynamics, and ThermodynamicsDas, Arya; Sahu, Pooja; Ali, Sk. MusharafJournal of Chemical & Engineering Data (2017), 62 (8), 2280-2295CODEN: JCEAAX; ISSN:0021-9568. (American Chemical Society)Mol. dynamics (MD) simulations were performed to calibrate the all-atom optimized potential for liq. simulations (OPLS-AA) force field using partial quantum charges calcd. from four different population anal. methods: Mulliken, Lowdin, NPA, and ChelpG for predicting the thermophys. properties of pure liqs. like tri-n-butylphosphate (TBP), tri-isoamylphosphate (TiAP), triethylphosphate (TEP), and dodecane to det. a potential solvent for the nuclear fuel cycle. The structural, dynamic, and thermodn. properties were calcd. in NVT ensembles by introducing the partial charges on each atom calcd. from d. functional theory (DFT). The calcd. structural and dynamic properties were affected by the different partial charges on TBP, TiAP, and TEP. The estd. liq. d. employing partial charges obtained from Mulliken population anal. with OPLS force field leads to an excellent agreement with the exptl. data (within 0.36-1.41%). The diffusivity and the pair correlation function (PCF) for all of the ligands have been calcd. and validated wherever literature data is available. The free energies of hydration and solvation for all of the ligands were evaluated using thermodn. integration technique and the hydration free energy for TEP is within 8.3% of the exptl. value, and for other properties they are not available in the literature for comparison. Furthermore, the partition coeff. of the ligands calcd. using MD derived free energy difference between the water-dodecane system resembles the trend predicted by DFT/COSMO-RS calcns. which is in qual. agreement with the exptl. results. Among the four-charge model, the computed dipole moment of TBP and TEP using the Mulliken charge is found to be in good agreement with the exptl. results. Finally, the superiority of TiAP over TBP as an extg. agent for the UO22+ ion has been demonstrated by a higher calcd. free energy of extn., ΔGext, over TBP using DFT. Overall the Mulliken charge embedded calibrated OPLS-AA force field is perhaps the most reliable one as it does not incorporate any arbitrary scaling in the force field or Lennard-Jones parameters and thus can be used indubitably to evaluate the liq. state properties of alkyl phosphates and n-alkanes which eventually assist in the invent of future generation extractants.
- 28Lennard-Jones, J. E. Cohesion. Proc. Phys. Soc. 1931, 43, 461– 482, DOI: 10.1088/0959-5309/43/5/301Google Scholar28https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaA38XjtVGl&md5=de693d597acc93abfe649daa7b96fd88CohesionLennard-Jones, J. E.Proceedings of the Physical Society, London (1931), 43 (), 461-82CODEN: PPSOAU; ISSN:0370-1328.A review. The various types of cohesive forces, such as homeopolar, ionic and metallic, are considered from the viewpoint of wave mechanics; each type is discussed. The nature of the van der Waals fields is treated and previous work, including the author's, is discussed.
- 29Eastwood, J.; Hockney, R.; Lawrence, D. P3M3DP-The three-dimensional periodic particle-particle/particle-mesh program. Comput. Phys. Commun. 1984, 35, 618– 619, DOI: 10.1016/S0010-4655(84)82783-6Google ScholarThere is no corresponding record for this reference.
- 30Nosé, S. Constant-temperature molecular dynamics. J. Phys.: Condens. Matter 1990, 2, SA115, DOI: 10.1088/0953-8984/2/S/013Google ScholarThere is no corresponding record for this reference.
- 31Hoover, W. G.; Posch, H. A.; Holian, B. L.; Gillan, M. J.; Mareschal, M.; Massobrio, C. Dissipative irreversibility from Nosé’s reversible mechanics. Mol. Simul. 1987, 1, 79– 86, DOI: 10.1080/08927028708080932Google ScholarThere is no corresponding record for this reference.
- 32Brehm, M.; Kirchner, B. TRAVIS - A free analyzer and visualizer for monte carlo and molecular dynamics trajectories. J. Chem. Inf. Model. 2011, 51, 2007– 2023, DOI: 10.1021/ci200217wGoogle Scholar32https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3MXptleqs7s%253D&md5=68eac025b6aaefd7961cba05b68e7ca3TRAVIS - A Free Analyzer and Visualizer for Monte Carlo and Molecular Dynamics TrajectoriesBrehm, Martin; Kirchner, BarbaraJournal of Chemical Information and Modeling (2011), 51 (8), 2007-2023CODEN: JCISD8; ISSN:1549-9596. (American Chemical Society)We present TRAVIS ("TRajectory Analyzer and VISualizer"), a free program package for analyzing and visualizing Monte Carlo and mol. dynamics trajectories. The aim of TRAVIS is to collect as many analyses as possible in one program, creating a powerful tool and making it unnecessary to use many different programs for evaluating simulations. This should greatly rationalize and simplify the work-flow of analyzing trajectories. TRAVIS is written in C++, open-source free-ware and licensed under the terms of the GNU General Public License (GPLv3). It is easy to install (platform independent, no external libraries) and easy to use. In this article, we present some of the algorithms that are implemented in TRAVIS - many of them widely known for a long time, but some of them also to appear in literature for the first time. All shown analyses only require a std. MD trajectory as input data.
- 33Brehm, M.; Weber, H.; Thomas, M.; Hollòczki, O.; Kirchner, B. Domain analysis in nanostructured liquids: A post-molecular dynamics study at the example of ionic liquids. ChemPhysChem 2015, 16, 3271– 3277, DOI: 10.1002/cphc.201500471Google Scholar33https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2MXhsFajtLjJ&md5=4723ce7f422ad2886440a6146c7f55afDomain Analysis in Nanostructured Liquids: A Post-Molecular Dynamics Study at the Example of Ionic LiquidsBrehm, Martin; Weber, Henry; Thomas, Martin; Holloczki, Oldamur; Kirchner, BarbaraChemPhysChem (2015), 16 (15), 3271-3277CODEN: CPCHFT; ISSN:1439-4235. (Wiley-VCH Verlag GmbH & Co. KGaA)In the present computational work, we develop a new tool for our trajectory anal. program TRAVIS to analyze the well-known behavior of liqs. to sep. into microphases. The dissection of the liq. into domains of different subsets, for example, in the case of fluorinated ionic liqs. with anionic and cationic head groups (forming together the polar domain), fluorous, and alkyl subsets is followed by radical Voronoi tessellation. This leads to useful av. quantities of the subset neighbor count, i.e., the domain count that gives the amt. of particular domains in the liq., the domain vol. and surface, as well as the isoperimetric quotient, which provides a measure of the deviation of the domains from a spherical shape. Thus, the newly implemented method allows anal. of the domains in terms of their nos. and shapes on a qual. and also quant. basis. It is a simple, direct, and automated anal. that does not need evaluation of the structure beforehand in terms of, for example, first solvent shell min. In the microheterogeneous ionic liqs. that we used as examples, the polar subsets always form a single domain in all investigated liqs. Although the fluorous side chains are also more or less connected in one or, maximally, two domains, the alkyl phases are dispersed.
- 34Luzar, A.; Chandler, D. Effect of environment on hydrogen bond dynamics in liquid water. Phys. Rev. Lett. 1996, 76, 928– 931, DOI: 10.1103/PhysRevLett.76.928Google Scholar34https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK28Xotl2iuw%253D%253D&md5=5f3b5cab4da60239ec27fae8d34b9a21Effect of environment on hydrogen bond dynamics in liquid waterLuzar, Alenka; Chandler, DavidPhysical Review Letters (1996), 76 (6), 928-31CODEN: PRLTAO; ISSN:0031-9007. (American Physical Society)In a series of mol.-dynamics calcns., we simulate the dynamics of forming and breaking a hydrogen bond in liq. water at room temp. We show that this dynamics is clearly nonexponential, yet virtually uncorrelated with the fluctuations of neighboring bonds.
- 35Luzar, A.; Chandler, D. Hydrogen-bond kinetics in liquid water. Nature 1996, 379, 55– 57, DOI: 10.1038/379055a0Google Scholar35https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK28XivVKhtQ%253D%253D&md5=282f1fbda5394e9ebd896db3892a005cHydrogen-bond kinetics in liquid waterLuzar, Alenka; Chandler, DavidNature (London) (1996), 379 (6560), 55-7CODEN: NATUAS; ISSN:0028-0836. (Macmillan Magazines)Hydrogen bonds play a crucial role in the behavior of water; their spatial patterns and fluctuations characterize the structure and dynamics of the liq. The processes of breaking and making hydrogen bonds in the condensed phase can be probed indirectly by a variety of exptl. techniques, and more quant. information can be obtained from computer simulations. In particular, simulations have revealed that on long timescales the relaxation behavior of hydrogen bonds in liq. water exhibit non-exponential kinetics, suggesting that bond making and breaking are not simple processes characterized by well defined rate consts. Here we show that these kinetics can be understood in terms of an interplay between diffusion and hydrogen-bond dynamics. In our model, which can be extended to other hydrogen-bonded liqs., diffusion governs whether a specific pair of water mols. are near neighbors, and hydrogen bonds between such pairs form and persist at random with av. lifetimes detd. by rate consts. for bond making and breaking.
- 36Luzar, A. Resolving the hydrogen bond dynamics conundrum. J. Chem. Phys. 2000, 113, 10663– 10675, DOI: 10.1063/1.1320826Google Scholar36https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3cXosFantb4%253D&md5=bc078176c816b6143f48bac25241748cResolving the hydrogen bond dynamics conundrumLuzar, AlenkaJournal of Chemical Physics (2000), 113 (23), 10663-10675CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)This paper analyzes dynamic properties of hydrogen bonds in liq. water. We use mol. dynamics simulation to calc. different probability densities that govern the time evolution of the formation and rupture of hydrogen bonds. We provide anal. connections between these functions. Excellent agreement with our simulation results is obsd. We prove transition state theory rate const. to be identical to the inverse of the assocd. mean first passage time (hydrogen bond lifetime). Hence, the anal. establishes its Arrhenius temp. dependence. We give the explicit relation between reactive flux correlation function for the relaxation dynamics of hydrogen bonds, and their first passage time probability densities. All the different observations in the existing literature, assocd. with various ests. of hydrogen bonding times in liq. water that are affected (or not affected) by particular bond criteria, as well as by different definitions of hydrogen bond lifetimes applied in simulation, can be easily reconciled within the framework of reactive flux correlation function approach.
- 37Gehrke, S.; von Domaros, M.; Clark, R.; Hollòczki, O.; Brehm, M.; Welton, T.; Luzar, A.; Kirchner, B. Structure and lifetimes in ionic liquids and their mixtures. Faraday Discuss. 2018, 206, 219– 245, DOI: 10.1039/C7FD00166EGoogle Scholar37https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXhtFCltbjM&md5=156128b1e8a1f4c6214e7e3b193e5719Structure and lifetimes in ionic liquids and their mixturesGehrke, Sascha; von Domaros, Michael; Clark, Ryan; Holloczki, Oldamur; Brehm, Martin; Welton, Tom; Luzar, Alenka; Kirchner, BarbaraFaraday Discussions (2018), 206 (Ionic Liquids: From Fundamental Properties to Practical Applications), 219-245CODEN: FDISE6; ISSN:1359-6640. (Royal Society of Chemistry)With the aid of mol. dynamics simulations, we study the structure and dynamics of different ionic liq. systems, with focus on hydrogen bond, ion pair and ion cage formation. To do so, we report radial distribution functions, their no. integrals, and various time-correlation functions, from which we ext. well-defined lifetimes by means of the reactive flux formalism. We explore the influence of polarizable force fields vs. non-polarizable ones with downscaled charges (±0.8) for the example of 1-butyl-3-methylimidazolium bromide. Furthermore, we use 1-butyl-3-methylimidazolium trifluoromethanesulfonate to investigate the impact of temp. and mixing with water as well as with the chloride ionic liq. Smaller coordination nos., larger distances, and tremendously accelerated dynamics are obsd. when the polarizable force field is applied. The same trends are found with increasing temp. Adding water decreases the ion-ion coordination nos. whereas the water-ion and water-water coordination is enhanced. A domain anal. reveals that the nonpolar parts of the ions are dispersed and when more water is added the water clusters increase in size. The dynamics accelerate in general upon addn. of water. In the ionic liq. mixt., the coordination no. around the cation changes between the two anions, but the no. integrals of the cation around the anions remain const. and the dynamics slow down with increasing content of the chloride ionic liq.
- 38Brehm, M.; Weber, H.; Thomas, M.; Hollòczki, O.; Kirchner, B. Domain analysis in nanostructured liquids: a post-molecular dynamics study at the example of ionic liquids. ChemPhysChem 2015, 16, 3271– 3277, DOI: 10.1002/cphc.201500471Google Scholar38https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2MXhsFajtLjJ&md5=4723ce7f422ad2886440a6146c7f55afDomain Analysis in Nanostructured Liquids: A Post-Molecular Dynamics Study at the Example of Ionic LiquidsBrehm, Martin; Weber, Henry; Thomas, Martin; Holloczki, Oldamur; Kirchner, BarbaraChemPhysChem (2015), 16 (15), 3271-3277CODEN: CPCHFT; ISSN:1439-4235. (Wiley-VCH Verlag GmbH & Co. KGaA)In the present computational work, we develop a new tool for our trajectory anal. program TRAVIS to analyze the well-known behavior of liqs. to sep. into microphases. The dissection of the liq. into domains of different subsets, for example, in the case of fluorinated ionic liqs. with anionic and cationic head groups (forming together the polar domain), fluorous, and alkyl subsets is followed by radical Voronoi tessellation. This leads to useful av. quantities of the subset neighbor count, i.e., the domain count that gives the amt. of particular domains in the liq., the domain vol. and surface, as well as the isoperimetric quotient, which provides a measure of the deviation of the domains from a spherical shape. Thus, the newly implemented method allows anal. of the domains in terms of their nos. and shapes on a qual. and also quant. basis. It is a simple, direct, and automated anal. that does not need evaluation of the structure beforehand in terms of, for example, first solvent shell min. In the microheterogeneous ionic liqs. that we used as examples, the polar subsets always form a single domain in all investigated liqs. Although the fluorous side chains are also more or less connected in one or, maximally, two domains, the alkyl phases are dispersed.
- 39Elfgen, R.; Hollòczki, O.; Kirchner, B. A molecular level understanding of template effects in ionic liquids. Acc. Chem. Res. 2017, 50, 2949– 2957, DOI: 10.1021/acs.accounts.7b00436Google Scholar39https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXhsl2ms7%252FL&md5=11aa98fbdf91a41743fabb670bb036a7A Molecular Level Understanding of Template Effects in Ionic LiquidsElfgen, Roman; Holloczki, Oldamur; Kirchner, BarbaraAccounts of Chemical Research (2017), 50 (12), 2949-2957CODEN: ACHRE4; ISSN:0001-4842. (American Chemical Society)The structure-directing or template effect has been invoked several times for ionic liqs. to explain the different outcome in material synthesis, namely, different scaffolds or geometrical arrangements with varying ionic liqs. It is obvious to assume that such an effect can originate from the most likely complex microstructure, being present within the ionic liq. itself. In that regard, ionic liqs. have already been shown to undergo a nanosegregation into polar and nonpolar phases, which is commonly known and denoted as microheterogeneity. In order to provide detailed insight on the mol. level and to understand the effects rising from this structuring, the authors performed mol. dynamics simulations on selected very simple model systems composed of 1-alkyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide, considering Et, Bu, hexyl, and octyl side chains attached to the cations, mixed with either n-dodecanol or n-butanol. By analyzing snapshots of the simulation boxes and calcg. spatial distribution functions, one can visualize that with increasing side chains, the systems show considerable nanosegregation into polar and nonpolar domains. Combined angular and distance distribution functions show that in case of the nanosegregating systems the side chains of the cations are preferentially arranged in a parallel fashion, which indicates a micelle-like structure for the ionic liqs. The alc. mols. participate in and are, therefore, influenced by this microheterogeneity. It can be shown that in the case of the short IL alkyl side chains, the self-aggregation of the nonpolar units of the alcs. is much stronger, while for the long chain cations, the nonpolar entities of the alcs. are most often connected to the nonpolar units of the ionic liqs. Using our domain anal. tool, we can quantify these observations by tracking the no., size, and shape of the polar and nonpolar entities present in the different investigated systems. The aforementioned combined angular-distance distribution functions reveal a structure-directing effect of the ionic liqs. on the alc. mols. within our simple model systems. The ionic liqs. act as template and order the alc. mols. according to their own structure, resulting in a parallel alignment of the alkyl side chains of the alcs. and ionic liq. cations, with both polar groups being at the same side. These observations show that the microheterogeneous structure of ionic liqs. can indeed be applied to order substrates with respect to each other or, for example, to catalysts in a predetd. fashion, opening new possibilities for explaining or enhancing selectivities of chem. reactions in ionic liqs.
- 40Macchieraldo, R.; Esser, L.; Elfgen, R.; Voepel, P.; Zahn, S.; Smarsly, B. M.; Kirchner, B. Hydrophilic ionic liquid mixtures of weakly and strongly coordinating anions with and without water. ACS Omega 2018, 3, 8567– 8582, DOI: 10.1021/acsomega.8b00995Google Scholar40https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXhsVajsrrF&md5=232946b7312d7bca99099aea1a5477b6Hydrophilic Ionic Liquid Mixtures of Weakly and Strongly Coordinating Anions with and without WaterMacchieraldo, Roberto; Esser, Lars; Elfgen, Roman; Voepel, Pascal; Zahn, Stefan; Smarsly, Bernd M.; Kirchner, BarbaraACS Omega (2018), 3 (8), 8567-8582CODEN: ACSODF; ISSN:2470-1343. (American Chemical Society)With the aid of ab initio mol. dynamics simulations, we investigate an ionic liq. mixt. composed of three components 1-butyl-3-methylimidazolium [C4C1Im]+, tetrafluoroborate [BF4]- and chloride [Cl]- without and with water. In the pure IL mixt., we observe an already complex network of interactions between cations and anions, adding water to the system even extends the complexity. Obsd. no. integrals show that the coordination no. between cations and anions is reduced in the system with water compared to the pure system. Further studies show, that the Coulombic network of the strongly coordinating anion [Cl]- is disturbed by water, while the weakly coordinating anion [BF4]- is not. These observations can also be confirmed by the Voronoi polyhedra anal., which shows that the polar network of microheterogeneous IL collapses by the introduction of water. Hydrogen-acceptor interactions revealed, that the [Cl]- anions are transferred from being situated in the IL to the water continuum while [BF4]-is almost unperturbed, these effects mainly influence the interplay of the ionic liq. network.
- 41Ray, P.; Vogl, T.; Balducci, A.; Kirchner, B. Structural investigations on lithium-doped protic and aprotic ionic liquids. J. Phys. Chem. B 2017, 121, 5279– 5292, DOI: 10.1021/acs.jpcb.7b02636Google Scholar41https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXntFeiuro%253D&md5=246661ee345211fbff9b320750518a8dStructural Investigations on Lithium-Doped Protic and Aprotic Ionic LiquidsRay, Promit; Vogl, Thomas; Balducci, Andrea; Kirchner, BarbaraJournal of Physical Chemistry B (2017), 121 (20), 5279-5292CODEN: JPCBFK; ISSN:1520-5207. (American Chemical Society)Solns. of lithium bis(trifluoromethanesulfonyl)imide (LiNTf2), in four different [NTf2]--based ionic liqs., are extensively investigated as potential electrolytes for lithium-ion batteries. Solvation of the [Li]+ ions in the ionic liqs. and its impact on their physicochem. properties are studied herein with the aid of mol. dynamics simulations. The cationic components of the investigated liqs. were systematically varied so as to individually evaluate effects of specific structural changes; increase in ring size, the addn. of an alkyl chain and absence of an acidic proton, on the solvation and mobility of the [Li]+ cations. The studied cations also allow for a direct comparison between solns. of [Li]+ salt in protic and aprotic ionic liqs. Emphasis is laid on elucidating the interactions between the [Li]+ and [NTf2]- ions revealing slightly higher coordination nos. for the aprotic solvent, benchmarked against exptl. measurements. The study suggests that the ionic liqs. largely retain their structure upon salt addn., with interactions within the liqs. only slightly perturbed. The rattling motion of the [Li]+ cations within cages formed by the surrounding [NTf2]- anions is examd. by the anal. of [Li]+ autocorrelation functions. Overall, the solvation mechanism of [Li]+ salt, within the hydrogen-bonded network of the ionic liqs., is detailed from classical and ab initio mol. dynamics simulations.
- 42Kirchner, B.; Hutter, J.; Kuo, I.-F. W.; Mundy, C. J. Hydrophobic hydratation from Car-Parrinello simulations. Int. J. Mod. Phys. B 2004, 18, 1951– 1962, DOI: 10.1142/S0217979204007241Google Scholar42https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2cXmvVaru7Y%253D&md5=0bd7f67008e20970e686112f086857baHydrophobic hydration from Car-Parrinello simulationsKirchner, Barbara; Hutter, Juerg; Kuo, I.-Feng W.; Mundy, Christopher J.International Journal of Modern Physics B: Condensed Matter Physics, Statistical Physics, Applied Physics (2004), 18 (14), 1951-1962CODEN: IJPBEV; ISSN:0217-9792. (World Scientific Publishing Co. Pte. Ltd.)In this work we investigate the fast anomalous diffusion of hydrogen mols. in water using Car-Parrinello mol. dynamics simulations. We employ Voronoi polyhedra anal. to distinguish between void diffusion and void hopping. Our results indicate that a combination of both mechanism is sufficient to explain anomalous diffusion. Furthermore, we investigate the geometry and the electronic structure of the first solvation shell.
- 43Ludwig, R. Water: From Clusters to the Bulk. Angew. Chem., Int. Ed. 2001, 40, 1808– 1827, DOI: 10.1002/1521-3773(20010518)40:10<1808::AID-ANIE1808>3.0.CO;2-1Google Scholar43https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3MXjvFyqtrw%253D&md5=6194b6257b7ef8884216627ec9d53132Water: from clusters to the bulkLudwig, RalfAngewandte Chemie, International Edition (2001), 40 (10), 1808-1827CODEN: ACIEF5; ISSN:1433-7851. (Wiley-VCH Verlag GmbH)A review with 186 refs. Water is of fundamental importance for human life and plays an important role in many biol. and chem. systems. Although water is the most abundant compd. on earth, it is definitely not a simple liq. It possesses strongly polar hydrogen bonds which are responsible for a striking set of anomalous phys. and chem. properties. For more than a century the combined importance and peculiarity of water inspired scientists to construct conceptual models, which in themselves reproduce the obsd. behavior of the liq. The exploration of structural and binding properties of small water complexes provides a key for understanding bulk water in its liq. and solid phase and for understanding solvation phenomena. Modern ab initio quantum chem. methods and high-resoln. spectroscopy methods have been extremely successful in describing such structures. Cluster models for liq. water try to mimic the transition from these clusters to bulk water. The important question is: What cluster properties are required to describe liq.-phase behavior.
- 44Huber, H.; Kirchner, B.; Searles, D. J. Is there an iceberg effect in the water/DMSO mixture? Some information from computational chemistry. J. Mol. Liq. 2002, 98-99, 71– 77, DOI: 10.1016/S0167-7322(01)00310-5Google Scholar44https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD38Xis1ehtbg%253D&md5=479106f6b2917627b4060648be17ac40Is There an Iceberg Effect in the Water/DMSO Mixture? Some Information from Computational Chemistry.Huber, Hanspeter; Kirchner, Barbara; Searles, Debra J.Journal of Molecular Liquids (2002), 98-99 (), 71-77CODEN: JMLIDT; ISSN:0167-7322. (Elsevier Science S.A.)The deuteron quadrupole coupling const. of heavy water in a mixt. with DMSO measured by B. C. Gordalla and M. D. Zeidler (1986)shows an unexpected behavior. A completely different evaluation from computational methods with independent approxns. does not reproduce this finding. A crit. discussion of all approxns. and uncertainties of the computational method is given using the limited results available for this system and results from studies on pure water.
- 45Galamba, N. Water’s structure around hydrophobic solutes and the iceberg model. J. Phys. Chem. B 2013, 117, 2153– 2159, DOI: 10.1021/jp310649nGoogle Scholar45https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3sXhsFGrtr0%253D&md5=4413978b1e29db158e633af3d0a40e03Water's Structure around Hydrophobic Solutes and the Iceberg ModelGalamba, N.Journal of Physical Chemistry B (2013), 117 (7), 2153-2159CODEN: JPCBFK; ISSN:1520-5207. (American Chemical Society)The structure of water in the hydration shells of small hydrophobic solutes was investigated through mol. dynamics. The results show that a subset of water mols. in the first hydration shell of a nonpolar solute have a significantly enhanced tetrahedrality and a slightly larger no. of hydrogen bonds, relative to the mols. in water at room temp., consistent with the exptl. obsd. neg. excess entropy and increased heat capacity of hydrophobic solns. at room temp. This ordering results from the rearrangement of a small no. of water mols. near the nonpolar solutes that occupy one to two vertices of the enhanced water tetrahedra. Although this structuring is not nearly like that often assocd. with a literal interpretation of the term "iceberg" in the Frank and Evans iceberg model, it does support a moderate interpretation of this model. Thus, the tetrahedral orientational order of this ensemble of water mols. is comparable to that of liq. water at ∼10°C, although not accompanied by the small contraction of the O-O distance obsd. in cold water. The structural changes of water in the vicinity of small nonpolar solutes cannot be inferred from the water radial distribution functions, explaining why this increased ordering is not obsd. through neutron diffraction expts. The present results restore a mol. view where the slower translational and reorientational dynamics of water near hydrophobic groups has a structural equiv. resembling water at low temps.
- 46Chowdhary, J.; Ladanyi, B. M. Hydrogen bond dynamics at the water/hydrocarbon interface. J. Phys. Chem. B 2009, 113, 4045– 4053, DOI: 10.1021/jp8061509Google Scholar46https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1cXhtl2ltLfN&md5=8f09abe1c2d114eded38bcc622dd70e7Hydrogen Bond Dynamics at the Water/Hydrocarbon InterfaceChowdhary, Janamejaya; Ladanyi, Branka M.Journal of Physical Chemistry B (2009), 113 (13), 4045-4053CODEN: JPCBFK; ISSN:1520-6106. (American Chemical Society)The dynamics of hydrogen bond formation and breakage for water in the vicinity of water/hydrocarbon liq. interfaces is studied using mol. dynamics simulations. Several liq. alkanes are considered as the hydrocarbon phase to det. the effects of their chain length and extent of branching on the properties of the adjacent water phase. In addn. to defining the interface location in terms of the lab.-frame d. profiles, the effects of interfacial fluctuations are considered by locating the interface in terms of the proximity of the mols. of the other phase. The hydrogen bond dynamics of interfacial water is weakly influenced by the identity of the hydrocarbon phase and by capillary waves. In addn. to calcg. hydrogen bond time correlations, the authors examine how the hydrogen bond dynamics depend on local coordination and det. the extent of cooperativity in the population relaxation of the hydrogen bonds that a given mol. participates in. The contributions of translational diffusion and reorientation of mol. O-H bonds to the mechanism of hydrogen bond breakage and reformation are studied. In previous work, rotation of the principal axes of water is anisotropic at the interface and depends on the initial orientation of the mol. relative to the interface. Here, the authors extend this anal. to the reorientation of the O-H vector and to hydrogen bond time correlation. Hydrogen bond dynamics are also sensitive to the initial orientation of the mols. participating in the hydrogen bond.
- 47Li, Q.; Wu, G.; Yu, Z. The role of methyl groups in the formation of hydrogen bond in DMSO-methanol mixtures. J. Am. Chem. Soc. 2006, 128, 1438– 1439, DOI: 10.1021/ja0569149Google Scholar47https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD28XjsFeltA%253D%253D&md5=5ece10fd302fa06d97bed75146bf6fb0The Role of Methyl Groups in the Formation of Hydrogen Bond in DMSO-Methanol MixturesLi, Qingzhong; Wu, Guoshi; Yu, ZhiwuJournal of the American Chemical Society (2006), 128 (5), 1438-1439CODEN: JACSAT; ISSN:0002-7863. (American Chemical Society)When examg. the formation energetics of a hydrogen-bonded complex R-X-H···Y-R', focus has been almost always on the atoms directly involved, namely the atoms X, Y, and H. Little attention has been paid to the effects of the secondary alkyl groups R and R'. Taking DMSO-methanol binary system as an example, we have studied the roles of the alkyl groups in stabilizing the hydrogen bonds by employing FTIR and NMR techniques and quantum chem. calcns. We found that Me groups play different roles in response to the hydrogen-bonding interactions. The Me groups of DMSO are electron-donating, whereas that of methanol is electron-withdrawing, both making pos. contributions. The findings reveal non-negligible effects of secondary alkyl groups in hydrogen bonding interaction and may shed light on the understanding of other more complicated hydrogen-bonded systems in chem. and biol. systems.
- 48Yao, H.; Ke, H.; Zhang, X.; Pan, S.-J.; Li, M.-S.; Yang, L.-P.; Schreckenbach, G.; Jiang, W. Molecular recognition of hydrophilic molecules in water by combining the hydrophobic effect with hydrogen bonding. J. Am. Chem. Soc. 2018, 140, 13466– 13477, DOI: 10.1021/jacs.8b09157Google Scholar48https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXhslKqt7rE&md5=bb3a11ceab70b41baf08ccdc55b71749Molecular Recognition of Hydrophilic Molecules in Water by Combining the Hydrophobic Effect with Hydrogen BondingYao, Huan; Ke, Hua; Zhang, Xiaobin; Pan, San-Jiang; Li, Ming-Shuang; Yang, Liu-Pan; Schreckenbach, Georg; Jiang, WeiJournal of the American Chemical Society (2018), 140 (41), 13466-13477CODEN: JACSAT; ISSN:0002-7863. (American Chemical Society)During the last half a century, great achievements have been made in mol. recognition in parallel with the invention of numerous synthetic receptors. However, the selective recognition of hydrophilic mols. in water remains a generally accepted challenge in supramol. chem. but is commonplace in nature. In an earlier Communication [Huang et al. J. Am. Chem. Soc. 2016, 138, 14550], we reported a pair of endo-functionalized mol. tubes that surprisingly prefer highly hydrophilic mols. over hydrophobic mols. of a similar size and shape. The hydrophobic effect and hydrogen bonding were proposed to be responsible, but their exact roles were not fully elucidated. In this Article, we present a thorough study on the binding behavior of these mol. tubes toward 44 hydrophilic mols. in water. Principal component anal. reveals that the binding strength is weakly correlated to the hydrophobicity, vol., surface area, and dipole moment of guests. Furthermore, mol. dynamics simulations show the hydrophobic effect through releasing the poorly hydrogen-bonded cavity water contributes to the binding of all the hydrophilic mols., while hydrogen bonding differentiates these mols. and is thus the key to achieve a high selectivity toward certain hydrophilic mols. over other mols. with a similar size and shape. Therefore, a good guest for these mol. tubes should meet the following criteria: the hydrogen-bonding sites should be complementary, and the mol. vol. should be large enough to expel all the cavity water but not too large to cause steric hindrance. This rule of thumb may also be used to design a selective receptor for certain hydrophilic mols. Following these guidelines, a "best-fit" guest was found for the syn-configured mol. tube with a binding const. as high as 106 M-1.
- 49Yamazoe, K.; Higaki, Y.; Inutsuka, Y.; Miyawaki, J.; Cui, Y.-T.; Takahara, A.; Harada, Y. Enhancement of the hydrogen-bonding network of water confined in a polyelectrolyte brush. Langmuir 2017, 33, 3954– 3959, DOI: 10.1021/acs.langmuir.7b00243Google Scholar49https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXltlWhs7o%253D&md5=9b6e57791c57cdd8e371998949ea8b52Enhancement of the Hydrogen-Bonding Network of Water Confined in a Polyelectrolyte BrushYamazoe, Kosuke; Higaki, Yuji; Inutsuka, Yoshihiro; Miyawaki, Jun; Cui, Yi-Tao; Takahara, Atsushi; Harada, YoshihisaLangmuir (2017), 33 (16), 3954-3959CODEN: LANGD5; ISSN:0743-7463. (American Chemical Society)Water existing in the vicinity of polyelectrolytes exhibits unique structural properties, which demonstrate key roles in chem., biol., and geoscience. In this study, X-ray absorption and emission spectroscopy was employed to observe the local hydrogen-bonding structure of water confined in a charged polyelectrolyte brush. Even at room temp., a majority of the water mols. confined in the polyelectrolyte brush exhibited one type of hydrogen-bonding configuration: a slightly distorted, albeit ordered configuration. The findings from this study provide new insights in terms of correlation between the function and local structure of water at the interface of biol. materials under physiol. conditions.
- 50Wallen, S. L.; Palmer, B. J.; Garrett, B. C.; Yonker, C. R. Density and temperature effects on the hydrogen bond structure of liquid methanol. J. Phys. Chem. 1996, 100, 3959– 3964, DOI: 10.1021/jp9524082Google Scholar50https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK28XhtVKrsrs%253D&md5=2414d722bd2b4ab2e4c78644b6a62b8cDensity and Temperature Effects on the Hydrogen Bond Structure of Liquid MethanolWallen, S. L.; Palmer, B. J.; Garrett, B. C.; Yonker, C. R.Journal of Physical Chemistry (1996), 100 (10), 3959-64CODEN: JPCHAX; ISSN:0022-3654. (American Chemical Society)The H bond structure of liq. MeOH was studied as a function of pressure and temp. up to 2.8 kbar and from 297 to 413 K. Chem. shifts of the CH3 and OH groups were monitored throughout this pressure and temp. regime, and the chem. shift difference between these two groups was used to describe changes of the H bond network in MeOH. The H bond equil. was studied using mol. dynamics simulations and a phenomenol. model describing clustering in liq. MeOH. Results are presented concerning the size and distribution of H-bonded clusters in MeOH as a function of pressure and temp. The extent of H bonding decreases upon an increase in temp. The results for pressure are equivocal, the phenomenol. model suggests that H bonding decreases with increasing pressure, which supports earlier interpretations regarding the measured self-diffusion coeffs. in deuterated MeOH as a function of pressure. The mol. dynamics simulations, however, show an increase in H bonding with increasing pressure.
Cited By
This article is cited by 7 publications.
- Mohsen Farshad, Mark J. DelloStritto, Antonio Suma, Vincenzo Carnevale. Detecting Liquid–Liquid Phase Separations Using Molecular Dynamics Simulations and Spectral Clustering. The Journal of Physical Chemistry B 2023, 127
(16)
, 3682-3689. https://doi.org/10.1021/acs.jpcb.3c00805
- Zheng Li, Brecht Dewulf, Koen Binnemans. Nonaqueous Solvent Extraction for Enhanced Metal Separations: Concept, Systems, and Mechanisms. Industrial & Engineering Chemistry Research 2021, 60
(48)
, 17285-17302. https://doi.org/10.1021/acs.iecr.1c02287
- Patryk Sikorski, Wojciech Tomaszewski, Katarzyna Kosiek, Tadeusz Hofman. The phase split forced by salts or carbohydrates in nonaqueous systems. Fluid Phase Equilibria 2023, 570 , 113790. https://doi.org/10.1016/j.fluid.2023.113790
- Maria Atanassova, Vanya Kurteva. Mutual Solubilities between Ethylene Glycol and Organic Diluents: Gas Chromatography and NMR. Molecules 2023, 28
(13)
, 5121. https://doi.org/10.3390/molecules28135121
- Aaron D. Wilson, Zi Hao Foo, Ashini S. Jayasinghe, Caleb Stetson, Hyeonseok Lee, Harry W. Rollins, Akshay Deshmukh, John H. Lienhard. Modeling Henry's law and phase separations of water–NaCl–organic mixtures with solvation and ion-pairing. Physical Chemistry Chemical Physics 2023, 13 https://doi.org/10.1039/D3CP02003G
- Roberto Macchieraldo, Johannes Ingenmey, Barbara Kirchner. Understanding the Complex Surface Interplay for Extraction: A Molecular Dynamics Study. Chemistry – A European Journal 2020, 26
(65)
, 14969-14977. https://doi.org/10.1002/chem.202002744
- M. Brehm, M. Thomas, S. Gehrke, B. Kirchner. TRAVIS—A free analyzer for trajectories from molecular simulation. The Journal of Chemical Physics 2020, 152
(16)
https://doi.org/10.1063/5.0005078
Abstract
Figure 1
Figure 1. Schematic representation of the procedure described in this work.
Figure 2
Figure 2. Snapshot taken from the MLDT system. All atoms are represented by their van der Waals radii. Color scheme: yellow for LiCl, red for MeOH, green for TBP, and blue for DD.
Figure 3
Figure 3. Mean square displacements of the bottom phase components. Top left: MeOH; top right: DD; bottom left: Li+ and Cl–; bottom right: DD and TBP.
Figure 4
Figure 4. RDFs of the interactions of MeOH molecules in the bottom phases.
Figure 5
Figure 5. Example of the liquid structure of MeOH in the presence of LiCl taken from the MLDT system. LiCl is represented with van der Waals radii while MeOH is displayed with the ball-and-stick models. Color scheme: brown for Li, green for Cl, orange for C, red for O, and white for H.
Figure 6
Figure 6. Intramolecular RDFs of the distance between the terminal carbons of DD molecules.
References
ARTICLE SECTIONSThis article references 50 other publications.
- 1Hiskey, B. Kirk-Othmer Encyclopedia of Chemical Technology; Wiley Blackwell: Hoboken, New Jersey, 2000.Google ScholarThere is no corresponding record for this reference.
- 2Anderson, C. G.; Dunne, R. C.; Uhrie, J. L. Mineral Processing and Extractive Metallurgy: 100 years of Innovation; Society for Mining, Metallurgy, and Exploration: Englewood, 2014.Google ScholarThere is no corresponding record for this reference.
- 3Hayes, P. C. Process Principles in Minerals and Materials Production; Hayes Publishing Co: Sherwood, 2003.Google ScholarThere is no corresponding record for this reference.
- 4Geissdoerfer, M.; Savaget, P.; Bocken, N. M.; Hultink, E. J. The circular economy - a new sustainability paradigm?. J. Clean. Prod. 2017, 143, 757– 768, DOI: 10.1016/j.jclepro.2016.12.048Google ScholarThere is no corresponding record for this reference.
- 5Binnemans, K.; Jones, P. T. Solvometallurgy: an emerging branch of extractive netallurgy. J. Sustainable Metall. 2017, 3, 570– 600, DOI: 10.1007/s40831-017-0128-2Google ScholarThere is no corresponding record for this reference.
- 6Abbott, A. P.; Boothby, D.; Capper, G.; Davies, D. L.; Rasheed, R. K. Deep eutectic solvents formed between choline chloride and carboxylic acids: versatile alternatives to ionic liquids. J. Am. Chem. Soc. 2004, 126, 9142– 9147, DOI: 10.1021/ja048266jGoogle Scholar6https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2cXlt1Gns70%253D&md5=5eba50624615fcd43af5afb2b11c0808Deep Eutectic Solvents Formed between Choline Chloride and Carboxylic Acids: Versatile Alternatives to Ionic LiquidsAbbott, Andrew P.; Boothby, David; Capper, Glen; Davies, David L.; Rasheed, Raymond K.Journal of the American Chemical Society (2004), 126 (29), 9142-9147CODEN: JACSAT; ISSN:0002-7863. (American Chemical Society)Deep Eutectic Solvents (DES) can be formed between a variety of quaternary ammonium salts and carboxylic acids. The phys. properties are significantly affected by the structure of the carboxylic acid but the phase behavior of the mixts. can be simply modeled by taking account of the mole fraction of carboxylic acid in the mixt. The phys. properties such as viscosity, cond., and surface tension of these DES are similar to ambient temp. ionic liqs. and insight into the cause of these properties is gained using hole-theory. It is shown that the cond. and viscosity of these liqs. is controlled by ion mobility and the availability of voids of suitable dimensions, and this is consistent with the fluidity of other ionic liqs. and molten salts. The DES are also shown to be good solvents for metal oxides, which could have potential application for metal extn.
- 7Rogers, R. D.; Seddon, K. R. Ionic liquids-solvents of the future?. Science 2003, 302, 792– 793, DOI: 10.1126/science.1090313Google Scholar7https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A280%3ADC%252BD3srhsV2jug%253D%253D&md5=0d242bb94d4f3a63249d106c6dab9e79Chemistry. Ionic liquids--solvents of the future?Rogers Robin D; Seddon Kenneth RScience (New York, N.Y.) (2003), 302 (5646), 792-3 ISSN:.There is no expanded citation for this reference.
- 8Gutowski, K. E.; Broker, G. A.; Willauer, H. D.; Huddleston, J. G.; Swatloski, R. P.; Holbrey, J. D.; Rogers, R. D. Controlling the aqueous miscibility of ionic liquids: aqueous biphasic systems of Weter-miscible ionic liquids and water-structuring salts for recycle, metathesis, and separations. J. Am. Chem. Soc. 2003, 125, 6632– 6633, DOI: 10.1021/ja0351802Google Scholar8https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3sXjs12ks7k%253D&md5=c46b031e62907c751321aeab5c7cee1fControlling the Aqueous Miscibility of Ionic Liquids: Aqueous Biphasic Systems of Water-Miscible Ionic Liquids and Water-Structuring Salts for Recycle, Metathesis, and SeparationsGutowski, Keith E.; Broker, Grant A.; Willauer, Heather D.; Huddleston, Jonathan G.; Swatloski, Richard P.; Holbrey, John D.; Rogers, Robin D.Journal of the American Chemical Society (2003), 125 (22), 6632-6633CODEN: JACSAT; ISSN:0002-7863. (American Chemical Society)Hydrophilic ionic liqs. can be salted-out and concd. from aq. soln. upon addn. of salts forming aq. biphasic systems as illustrated by the phase behavior of mixts. of 1-butyl-3-methylimidazolium chloride and K3PO4.
- 9Earle, M. J.; Esperança, J. M.; Gilea, M. A.; Lopes, J. N. C.; Rebelo, L. P.; Magee, J. W.; Seddon, K. R.; Widegren, J. A. The distillation and volatility of ionic liquids. Nature 2006, 439, 831– 834, DOI: 10.1038/nature04451Google Scholar9https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD28XhsVSkurk%253D&md5=9e8d177c2100d4ecc45c62f3ca92c357The distillation and volatility of ionic liquidsEarle, Martyn J.; Esperanca, Jose M. S. S.; Gilea, Manuela A.; Canongia Lopes, Jose N.; Rebelo, Luis P. N.; Magee, Joseph W.; Seddon, Kenneth R.; Widegren, Jason A.Nature (London, United Kingdom) (2006), 439 (7078), 831-834CODEN: NATUAS; ISSN:0028-0836. (Nature Publishing Group)It is widely believed that a defining characteristic of ionic liqs. (or low-temp. molten salts) is that they exert no measurable vapor pressure, and hence cannot be distd. Here we demonstrate that this is unfounded, and that many ionic liqs. can be distd. at low pressure without decompn. Ionic liqs. represent matter solely composed of ions, and so are perceived as non-volatile substances. During the last decade, interest in the field of ionic liqs. has burgeoned, producing a wealth of intellectual and technol. challenges and opportunities for the prodn. of new chem. and extractive processes, fuel cells and batteries, and new composite materials. Much of this potential is underpinned by their presumed involatility. This characteristic, however, can severely restrict the attainability of high purity levels for ionic liqs. (when they contain poorly volatile components) in recycling schemes, as well as excluding their use in gas-phase processes. We anticipate that our demonstration that some selected families of commonly used aprotic ionic liqs. can be distd. at 200-300 °C and low pressure, with concomitant recovery of significant amts. of pure substance, will permit these currently excluded applications to be realized.
- 10Fukumoto, K.; Yoshizawa, M.; Ohno, H. Room temperature ionic liquids from 20 Natural amino acids. J. Am. Chem. Soc. 2005, 127, 2398– 2399, DOI: 10.1021/ja043451iGoogle Scholar10https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2MXptlSgsA%253D%253D&md5=b8ffaba7d45996888c6c26708d46cfa9Room Temperature Ionic Liquids from 20 Natural Amino AcidsFukumoto, Kenta; Yoshizawa, Masahiro; Ohno, HiroyukiJournal of the American Chemical Society (2005), 127 (8), 2398-2399CODEN: JACSAT; ISSN:0002-7863. (American Chemical Society)The authors have succeeded in synthesizing room temp. ionic liqs. from 20 natural amino acids (AA). These amino acid ionic liqs. are H-AA-O-·(1-ethyl-3-methylimidazolium)+ salts, which dissolved native amino acids despite water-free conditions. Furthermore, these ionic liqs. are sol. in various org. solvents, such as chloroform. Effects of acidity, hydrogen bonding ability, and steric factors on the properties of these ionic liqs. were analyzed as the function of side groups.
- 11Wang, Y.; Voth, G. A. Unique spatial heterogeneity in ionic liquids. J. Am. Chem. Soc. 2005, 127, 12192– 12193, DOI: 10.1021/ja053796gGoogle Scholar11https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2MXnslOrtbc%253D&md5=139435d4120f25a320b1cae5af1e0fabUnique Spatial Heterogeneity in Ionic LiquidsWang, Yanting; Voth, Gregory A.Journal of the American Chemical Society (2005), 127 (35), 12192-12193CODEN: JACSAT; ISSN:0002-7863. (American Chemical Society)A multiscale coarse-graining model for ionic liqs. has been extended to investigate the unique aggregation of cations in ionic liqs. through computer simulation. It has been found that, with sufficiently long side chains, the tail groups of cations aggregate to form spatially heterogeneous domains, while headgroups of the cations and the anions distribute as uniformly as possible. This is understood as the result of competition between the charged electrostatic interactions between headgroups and anions and the collective short-range interactions between the neutral tail groups. This aggregation can help to explain a no. of exptl. obsd. phys. phenomena in ionic liqs.
- 12Scurto, A. M.; Aki, S. N. V. K.; Brennecke, J. F. CO2 as a separation switch for ionic liquid/organic mixtures. J. Am. Chem. Soc. 2002, 124, 10276– 10277, DOI: 10.1021/ja0268682Google Scholar12https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD38XlvFGlt7w%253D&md5=a7afd1d2406e4de0c70dbbfc7c3dba23CO2 as a Separation Switch for Ionic Liquid/Organic MixturesScurto, Aaron M.; Aki, Sudhir N. V. K.; Brennecke, Joan F.Journal of the American Chemical Society (2002), 124 (35), 10276-10277CODEN: JACSAT; ISSN:0002-7863. (American Chemical Society)A novel technique to sep. ionic liqs. from org. compds. is introduced which uses carbon dioxide to induce the formation of an ionic liq.-rich phase and an org.-rich liq. phase in mixts. of methanol and 3-butyl-1-methyl-imidazolium hexafluorophosphate ([C4mim][PF6]). If the temp. is above the crit. temp. of CO2 then the methanol-rich phase can become completely miscible with the CO2-rich phase, and this new phase is completely ionic liq.-free. Since CO2 is nonpolar, it is not equipped to solvate ions. As the CO2 dissolves in the methanol/[C4mim][PF6] mixt., the solvent power of the CO2-expanded liq. is significantly reduced, inducing the formation of the second liq. phase that is rich in ionic liq. This presents a new way to recover products from ionic liq. mixts. and purify org. phases that have been contaminated with ionic liq. Moreover, these results have important implications for reactions done in CO2/ionic liq. biphasic mixts.
- 13Sastri, V. Selective leaching of metals from ores with organic solvents and chelating agents. J. Sci. Ind. Res. 1975, 34, 663– 665Google Scholar13https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaE28Xks1Kiur8%253D&md5=b045526d899e5f5fb1ee928d8dd423a8Selective leaching of metals from ores with organic solvents and chelating agentsSastri, V. S.Journal of Scientific & Industrial Research (1975), 34 (12), 663-5CODEN: JSIRAC; ISSN:0022-4456.A review with 16 refs. is given on the use of org. solvents and chelating agents in the selective leaching of metals from ores.
- 14Raghavan, S.; Fuerstenau, D. A Lyometallurgical Process for Leaching Copper from Chrysocolla. In Hydrometallurgy Fundamentals, Technology and Innovations , Proceedings of the Milton E., Wordsworth (IV) Int. Symposium on Hydrometallurgy, 1993; pp 283– 297.Google ScholarThere is no corresponding record for this reference.
- 15Chizhevskaya, S. V.; Chekmarev, A. M.; Klimenko, O. M.; Povetkina, M. V.; Sinegribova, O. A.; Cox, M. In Nontraditional Methods of Treating High-Silicon Ores Containing Rare Elements , Hydrometallurgy ’94: Papers presented at the international symposium ‘Hydrometallurgy ’94’ organized by the Institution of Mining and Metallurgy and the Society of Chemical Industry, and held in Cambridge, England, from 11 to 15 July, 1994; pp 219– 228.Google ScholarThere is no corresponding record for this reference.
- 16Hiroshi, T.; Kazushige, T.; Yoshimi, O. Direct Dissolution of Water-Insoluble Uranium Compounds by Contact with Neutral Organic Solvents Pretreated with Nitric Acid. U.S. Patent US3288568, 1966.Google ScholarThere is no corresponding record for this reference.
- 17Batchu, N. K.; Vander Hoogerstraete, T.; Banerjee, D.; Binnemans, K. Non-aqueous solvent extraction of rare-earth nitrates from ethylene glycol to n-dodecane by Cyanex 923. Sep. Purif. Technol. 2017, 174, 544– 553, DOI: 10.1016/j.seppur.2016.10.039Google Scholar17https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28XhslyjsL%252FP&md5=6f6ebbff4bbd5464d240c22512f0cf07Non-aqueous solvent extraction of rare-earth nitrates from ethylene glycol to n-dodecane by Cyanex 923Batchu, Nagaphani Kumar; Vander Hoogerstraete, Tom; Banerjee, Dipanjan; Binnemans, KoenSeparation and Purification Technology (2017), 174 (), 544-553CODEN: SPUTFP; ISSN:1383-5866. (Elsevier B.V.)A solvent extn. process comprising two immiscible org. phases was developed for the extn. of rare earths. The more polar org. phase was ethylene glycol with dissolved rare-earth nitrate salts and lithium nitrate, while the less polar phase was a soln. of the neutral extractant Cyanex 923 dissolved in n-dodecane. The solvent extn. mechanism was detd. by slope anal. and the main species in the org. phase was identified by Extended x-ray Absorption Fine Structure (EXAFS) studies. The extn. from the ethylene glycol soln. was compared with extn. from an aq. feed soln. When compared to aq. feed solns., the light rare-earth elements (LREEs) are less efficiently extd. and the heavy rare-earth elements (HREEs) more efficiently extd. from an ethylene glycol feed soln., resulting into the easy sepn. of HREEs from LREEs. The sepn. factors between neighboring elements are higher for this non-aq. solvent extn. process than for extn. from an aq. feed soln.
- 18Li, Z.; Li, X.; Raiguel, S.; Binnemans, K. Separation of transition metals from rare earths by non-aqueous solvent extraction from ethylene glycol solutions using Aliquat 336. Sep. Purif. Technol. 2018, 201, 318– 326, DOI: 10.1016/j.seppur.2018.03.022Google Scholar18https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXlslKgs7k%253D&md5=b1a5cd1a9bc8ea4fa342fa21a0c64f00Separation of transition metals from rare earths by nonaqueous solvent extraction from ethylene glycol solutions using Aliquat 336Li, Zheng; Li, Xiaohua; Raiguel, Stijn; Binnemans, KoenSeparation and Purification Technology (2018), 201 (), 318-326CODEN: SPUTFP; ISSN:1383-5866. (Elsevier B.V.)Solvent extn. is a widely used sepn. technique in extractive metallurgy. A conventional solvent extn. system consists of an aq. phase and an immiscible org. phase. Replacement of H2O by a polar org. solvent can lead to superior metal sepns. Co(II) and Sm(III) chlorides dissolved in H2O and ethylene glycol (EG), resp., with LiCl as chloride source, were extd. by Aliquat 336 dild. in toluene. Both Co and Sm were extd. from the aq. solns., but Co was extd. more efficiently from the ethylene glycol solns. than from the aq. solns., whereas Sm was not extd. at all from the ethylene glycol solns. As a result, Co and Sm could be sepd. completely in a single extn. step from ethylene glycol solns. The mechanisms of Co extn. by Aliquat 336 from the ethylene glycol and aq. solns. are similar, as validated by slope anal. and UV-visible absorption spectroscopy. Sm was extd. from the aq. soln. through the salting-out effect of LiCl. LiCl has a much lower salting-out effect for Sm in ethylene glycol than in H2O due to the lower dielec. const. of ethylene glycol and the lower soly. of LiCl in ethylene glycol. Consequently, Sm is not salted out from ethylene glycol, leading to a very efficient sepn. of Co and Sm. This sepn. effect can also be applied to the sepn. of other transition metal and rare-earth metal pairs, including Fe/Nd and Zn/Eu.
- 19Misawa, M.; Yoshida, K.; Maruyama, K.; Munemura, H.; Hosokawa, Y. Salt-induced phase separation in aqueous solution. J. Phys. Chem. Solids 1999, 60, 1301– 1306, DOI: 10.1016/S0022-3697(99)00108-0Google Scholar19https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK1MXls1equrY%253D&md5=16a5c9ae3269dc0a0c8cf3ca669cca89Salt-induced phase separation in aqueous solutionMisawa, M.; Yoshida, K.; Maruyama, K.; Munemura, H.; Hosokawa, Y.Journal of Physics and Chemistry of Solids (1999), 60 (8-9), 1301-1306CODEN: JPCSAW; ISSN:0022-3697. (Elsevier Science Ltd.)KCl-induced phase sepn. phenomenon in 1-propanol aq. soln. has been studied by means of small angle (SANS) and quasi-elastic neutron scattering measurements. The significant increase of SANS intensity by the addn. of KCl is attributed to the aggregate of clusters, of which the fractal nature of dimension 1.8-1.9 is found. The motion of water mols. is much suppressed in the 1-propanol-water soln., while it is slightly recovered by the addn. of KCl. The aggregate of clusters and liq.-liq. phase sepn. are discussed in relation to the motion of water mols.
- 20Batchu, N. K.; Vander Hoogerstraete, T.; Banerjee, D.; Binnemans, K. Separation of rare-earth ions from ethylene glycol (+LiCl) solutions by non-aqueous solvent extraction with Cyanex 923. RSC Adv. 2017, 7, 45351– 45362, DOI: 10.1039/C7RA09144CGoogle Scholar20https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXhsFCktL7I&md5=4e3bd3a8ac008e12abeb98934d7b2d9fSeparation of rare-earth ions from ethylene glycol (+LiCl) solutions by non-aqueous solvent extraction with Cyanex 923Batchu, Nagaphani Kumar; Vander Hoogerstraete, Tom; Banerjee, Dipanjan; Binnemans, KoenRSC Advances (2017), 7 (72), 45351-45362CODEN: RSCACL; ISSN:2046-2069. (Royal Society of Chemistry)The sepn. of a mixt. of rare earths by non-aq. solvent extn. with two immiscible org. phases has been studied. The more polar org. phase was ethylene glycol with dissolved lithium chloride and the less polar org. phase was the extractant dild. in n-dodecane. Cyanex 923 was found to be the most performant extractant amongst the investigated acidic, basic and solvating extractants: Cyanex 272, Cyphos IL 101, Aliquat 336, bis(2-ethylhexyl)amine, trioctylphosphine oxide (TOPO) and Cyanex 923. The replacement of the aq. chloride feed solns. by non-aq. ethylene glycol feed solns. had a profound effect on the distribution ratios and sepn. factors. The sepn. factors for extn. of pairs of rare earths from aq. chloride solns. by Cyanex 923 are too low to be of practical use. On the contrary, a mixt. of rare earths can be sepd. conveniently in four different groups by extn. with Cyanex 923 from ethylene glycol (+LiCl) solns. The influence of several parameters such as the chloride concn., the type of chloride salt, the addn. of other polar solvents to the ethylene glycol phase, the addn. of second extractant to the less polar org. phase, and the addn. of complexing agents to the ethylene glycol phase has been studied. The extn. mechanism for extn. of ytterbium(III) was studied by slope anal. expts. The ytterbium(III) species in the ethylene glycol phase and the extd. species in the n-dodecane phase were detd. by EXAFS. Furthermore, a conceptual flow sheet for the fractionation of rare earths from an ethylene glycol (+LiCl) feed soln. into different groups by extn. with Cyanex 923 has been proposed. The new extn. system is useful for extn. of scandium and for sepn. of scandium from the other REEs.
- 21Casás, L. M.; Touriño, A.; Orge, B.; Marino, G.; Iglesias, M.; Tojo, J. Thermophysical properties of acetone or methanol + n–alkane (C9 to C12) mixtures. J. Chem. Eng. Data 2002, 47, 887– 893, DOI: 10.1021/je0103059Google Scholar21https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD38XktV2gsLw%253D&md5=dc84bff396c2f3d1c9ef3216dbafb093Thermophysical Properties of Acetone or Methanol + n-Alkane (C9 to C12) MixturesCasas, Lidia M.; Tourino, Antonio; Orge, Beatriz; Marino, Gonzalo; Iglesias, Miguel; Tojo, JoseJournal of Chemical and Engineering Data (2002), 47 (4), 887-893CODEN: JCEAAX; ISSN:0021-9568. (American Chemical Society)Liq.-liq. equil. for methanol + n-alkanes (C9 to C12) were detd. in the temp. range from 278.15 to 308.15 K, using a visual static method. Densities, refractive indexes, and speeds of sound were also measured at 288.15, 298.15, and 308.15 K for acetone + n-alkane mixts. and at 298.15 and 308.15 K for methanol + n-alkane mixts. to study the temp. influence on mixing phenomena. The results were used to model liq.-liq. coexistence and derived properties trend. The calcd. values based on the UNIQUAC equation were found to be similar to those based on the NRTL model. The Prigogine-Flory-Paterson theory was applied to est. the excess molar volumes and the UNIFAC group contribution model for phase equil.
- 22Casás, L. M.; Orge, B.; Díaz, C.; Tojo, J. Liquid-liquid equilibria for mixtures of methyl acetate+methanol+n-alkane (C10–C12) at several temperatures and 1 atm. J. Chem. Thermodyn. 2004, 36, 237– 243, DOI: 10.1016/j.jct.2003.11.009Google Scholar22https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2cXhtVensL0%253D&md5=98c0ac16ca68e3c3d525015980b6cdebLiquid-liquid equilibria for mixtures of {methyl acetate + methanol + n-alkane (C10-C12)} at several temperatures and 1 atmCasas, Lidia M.; Orge, Beatriz; Diaz, Concepcion; Tojo, JoseJournal of Chemical Thermodynamics (2004), 36 (3), 237-243CODEN: JCTDAF; ISSN:0021-9614. (Elsevier Science Ltd.)The liq.-liq. equil. of the mixts. {methyl acetate + methanol + n-alkane (C10-C12)} at atm. pressure in the temp. range (278.15 to 308.15) K are presented. The equil. compns. were measured by gas chromatog. from samples splitted isothermally into a glass-stirred device and phase diagrams are reported at each temp. The exptl. detd. liq.-liq. equil. were satisfactorily correlated by the NRTL and UNIQUAC equations. A poor description has been obtained using the UNIFAC method and their modifications.
- 23Touriño, A.; Casás, L.; Marino, G.; Iglesias, M.; Orge, B.; Tojo, J. Liquid phase behaviour and thermodynamics of acetone+methanol+n-alkane (C9–C12) mixtures. Fluid Phase Equilib. 2003, 206, 61– 85, DOI: 10.1016/S0378-3812(02)00306-0Google Scholar23https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3sXjtl2qu70%253D&md5=0ff46624a178726290d28c9d57a0666fLiquid phase behaviour and thermodynamics of acetone+methanol+n-alkane (C9-C12) mixturesTourino, A.; Casas, L. M.; Marino, G.; Iglesias, M.; Orge, B.; Tojo, J.Fluid Phase Equilibria (2003), 206 (1-2), 61-85CODEN: FPEQDT; ISSN:0378-3812. (Elsevier Science B.V.)This paper reports the results of a new exptl. study of two-liq. phases equil. as a function of temp. and one-liq. phase thermodn. magnitudes (densities, refractive indexes and speeds of sound), covering the compn. diagrams for the mixts. acetone + methanol + n-alkane (C9-C12). The obtained coexistence curves were very asym. with respect to equimolar alc. + n-alkane compn., such effect increasing with the length of the aliph. chain and temp. A comparative anal. was performed by application of different methods to predict exptl. liq.-liq. equil. (LLE) behavior and thermodn. of these ternary mixts. The obtained exptl. results let us know the potential role of the linear aliph. alkanes as sepn. agents for heterogeneous modified distn. of the azeotrope acetone + methanol.
- 24Martínez, L.; Andrade, R.; Birgin, E. G.; Martínez, J. M. A package for building initial configurations for molecular dynamics simulations. J. Comput. Chem. 2009, 30, 2157– 2164, DOI: 10.1002/jcc.21224Google Scholar24https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1MXptleqsb8%253D&md5=2a76255c873b866a26540f7e84496272PACKMOL: A package for building initial configurations for molecular dynamics simulationsMartinez, L.; Andrade, R.; Birgin, E. G.; Martinez, J. M.Journal of Computational Chemistry (2009), 30 (13), 2157-2164CODEN: JCCHDD; ISSN:0192-8651. (John Wiley & Sons, Inc.)Adequate initial configurations for mol. dynamics simulations consist of arrangements of mols. distributed in space in such a way to approx. represent the system's overall structure. In order that the simulations are not disrupted by large van der Waals repulsive interactions, atoms from different mols. must keep safe pairwise distances. Obtaining such a mol. arrangement can be considered a packing problem: Each type mol. must satisfy spatial constraints related to the geometry of the system, and the distance between atoms of different mols. must be greater than some specified tolerance. We have developed a code able to pack millions of atoms, grouped in arbitrarily complex mols., inside a variety of three-dimensional regions. The regions may be intersections of spheres, ellipses, cylinders, planes, or boxes. The user must provide only the structure of one mol. of each type and the geometrical constraints that each type of mol. must satisfy. Building complex mixts., interfaces, solvating biomols. in water, other solvents, or mixts. of solvents, is straightforward. In addn., different atoms belonging to the same mol. may also be restricted to different spatial regions, in such a way that more ordered mol. arrangements can be built, as micelles, lipid double-layers, etc. The packing time for state-of-the-art mol. dynamics systems varies from a few seconds to a few minutes in a personal computer. The input files are simple and currently compatible with PDB, Tinker, Molden, or Moldy coordinate files. The package is distributed as free software and can be downloaded from . © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009.
- 25Plimpton, S. Fast parallel algorithms for short-range molecular dynamics. J. Comput. Phys. 1995, 117, 1– 19, DOI: 10.1006/jcph.1995.1039Google Scholar25https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK2MXlt1ejs7Y%253D&md5=715052332237e4cf8501f8fb81234017Fast parallel algorithms for short-range molecular dynamicsPlimpton, SteveJournal of Computational Physics (1995), 117 (1), 1-19CODEN: JCTPAH; ISSN:0021-9991.Three parallel algorithms for classical mol. dynamics are presented. The first assigns each processor a fixed subset of atoms; the second assigns each a fixed subset of inter-at. forces to compute; the third assigns each a fixed spatial region. The algorithms are suitable for mol. dynamics models which can be difficult to parallelize efficiently - those with short-range forces where the neighbors of each atom change rapidly. They can be implemented on any distributed-memory parallel machine which allows for message-passing of data between independently executing processors. The algorithms are tested on a std. Lennard-Jones benchmark problem for system sizes ranging from 500 to 100,000,000 atoms on several parallel supercomputers - the nCUBE 2, Intel iPSC/860 and Paragon, and Cray T3D. Comparing the results to the fastest reported vectorized Cray Y-MP and C90 algorithm shows that the current generation of parallel machines is competitive with conventional vector supercomputers even for small problems. For large problems, the spatial algorithm achieves parallel efficiencies of 90% and a 1840-node Intel Paragon performs up to 165 faster than a single Cray C90 processor. Trade-offs between the three algorithms and guidelines for adapting them to more complex mol. dynamics simulations are also discussed.
- 26Jorgensen, W. L.; Maxwell, D. S.; Tirado-Rives, J. Development and testing of the OPLS all-atom force field on conformational energetics and properties of organic liquids. J. Am. Chem. Soc. 1996, 118, 11225– 11236, DOI: 10.1021/ja9621760Google Scholar26https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK28XmtlOitrs%253D&md5=fef2924a69421881390282aa309ae91bDevelopment and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic LiquidsJorgensen, William L.; Maxwell, David S.; Tirado-Rives, JulianJournal of the American Chemical Society (1996), 118 (45), 11225-11236CODEN: JACSAT; ISSN:0002-7863. (American Chemical Society)The parametrization and testing of the OPLS all-atom force field for org. mols. and peptides are described. Parameters for both torsional and nonbonded energetics have been derived, while the bond stretching and angle bending parameters have been adopted mostly from the AMBER all-atom force field. The torsional parameters were detd. by fitting to rotational energy profiles obtained from ab initio MO calcns. at the RHF/6-31G*//RHF/6-31G* level for more than 50 org. mols. and ions. The quality of the fits was high with av. errors for conformational energies of less than 0.2 kcal/mol. The force-field results for mol. structures are also demonstrated to closely match the ab initio predictions. The nonbonded parameters were developed in conjunction with Monte Carlo statistical mechanics simulations by computing thermodn. and structural properties for 34 pure org. liqs. including alkanes, alkenes, alcs., ethers, acetals, thiols, sulfides, disulfides, aldehydes, ketones, and amides. Av. errors in comparison with exptl. data are 2% for heats of vaporization and densities. The Monte Carlo simulations included sampling all internal and intermol. degrees of freedom. It is found that such non-polar and monofunctional systems do not show significant condensed-phase effects on internal energies in going from the gas phase to the pure liqs.
- 27Das, A.; Sahu, P.; Ali, S. M. Molecular dynamics simulation for the calibration of the OPLS force field using DFT derived partial charges for the screening of alkyl phosphate ligands by studying structure, dynamics, and thermodynamics. J. Chem. Eng. Data 2017, 62, 2280– 2295, DOI: 10.1021/acs.jced.7b00096Google Scholar27https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXht1Slt73E&md5=2ec0cfaf7877976bc614ed86dea8dd79Molecular Dynamics Simulation for the Calibration of the OPLS Force Field Using DFT Derived Partial Charges for the Screening of Alkyl Phosphate Ligands by Studying Structure, Dynamics, and ThermodynamicsDas, Arya; Sahu, Pooja; Ali, Sk. MusharafJournal of Chemical & Engineering Data (2017), 62 (8), 2280-2295CODEN: JCEAAX; ISSN:0021-9568. (American Chemical Society)Mol. dynamics (MD) simulations were performed to calibrate the all-atom optimized potential for liq. simulations (OPLS-AA) force field using partial quantum charges calcd. from four different population anal. methods: Mulliken, Lowdin, NPA, and ChelpG for predicting the thermophys. properties of pure liqs. like tri-n-butylphosphate (TBP), tri-isoamylphosphate (TiAP), triethylphosphate (TEP), and dodecane to det. a potential solvent for the nuclear fuel cycle. The structural, dynamic, and thermodn. properties were calcd. in NVT ensembles by introducing the partial charges on each atom calcd. from d. functional theory (DFT). The calcd. structural and dynamic properties were affected by the different partial charges on TBP, TiAP, and TEP. The estd. liq. d. employing partial charges obtained from Mulliken population anal. with OPLS force field leads to an excellent agreement with the exptl. data (within 0.36-1.41%). The diffusivity and the pair correlation function (PCF) for all of the ligands have been calcd. and validated wherever literature data is available. The free energies of hydration and solvation for all of the ligands were evaluated using thermodn. integration technique and the hydration free energy for TEP is within 8.3% of the exptl. value, and for other properties they are not available in the literature for comparison. Furthermore, the partition coeff. of the ligands calcd. using MD derived free energy difference between the water-dodecane system resembles the trend predicted by DFT/COSMO-RS calcns. which is in qual. agreement with the exptl. results. Among the four-charge model, the computed dipole moment of TBP and TEP using the Mulliken charge is found to be in good agreement with the exptl. results. Finally, the superiority of TiAP over TBP as an extg. agent for the UO22+ ion has been demonstrated by a higher calcd. free energy of extn., ΔGext, over TBP using DFT. Overall the Mulliken charge embedded calibrated OPLS-AA force field is perhaps the most reliable one as it does not incorporate any arbitrary scaling in the force field or Lennard-Jones parameters and thus can be used indubitably to evaluate the liq. state properties of alkyl phosphates and n-alkanes which eventually assist in the invent of future generation extractants.
- 28Lennard-Jones, J. E. Cohesion. Proc. Phys. Soc. 1931, 43, 461– 482, DOI: 10.1088/0959-5309/43/5/301Google Scholar28https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaA38XjtVGl&md5=de693d597acc93abfe649daa7b96fd88CohesionLennard-Jones, J. E.Proceedings of the Physical Society, London (1931), 43 (), 461-82CODEN: PPSOAU; ISSN:0370-1328.A review. The various types of cohesive forces, such as homeopolar, ionic and metallic, are considered from the viewpoint of wave mechanics; each type is discussed. The nature of the van der Waals fields is treated and previous work, including the author's, is discussed.
- 29Eastwood, J.; Hockney, R.; Lawrence, D. P3M3DP-The three-dimensional periodic particle-particle/particle-mesh program. Comput. Phys. Commun. 1984, 35, 618– 619, DOI: 10.1016/S0010-4655(84)82783-6Google ScholarThere is no corresponding record for this reference.
- 30Nosé, S. Constant-temperature molecular dynamics. J. Phys.: Condens. Matter 1990, 2, SA115, DOI: 10.1088/0953-8984/2/S/013Google ScholarThere is no corresponding record for this reference.
- 31Hoover, W. G.; Posch, H. A.; Holian, B. L.; Gillan, M. J.; Mareschal, M.; Massobrio, C. Dissipative irreversibility from Nosé’s reversible mechanics. Mol. Simul. 1987, 1, 79– 86, DOI: 10.1080/08927028708080932Google ScholarThere is no corresponding record for this reference.
- 32Brehm, M.; Kirchner, B. TRAVIS - A free analyzer and visualizer for monte carlo and molecular dynamics trajectories. J. Chem. Inf. Model. 2011, 51, 2007– 2023, DOI: 10.1021/ci200217wGoogle Scholar32https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3MXptleqs7s%253D&md5=68eac025b6aaefd7961cba05b68e7ca3TRAVIS - A Free Analyzer and Visualizer for Monte Carlo and Molecular Dynamics TrajectoriesBrehm, Martin; Kirchner, BarbaraJournal of Chemical Information and Modeling (2011), 51 (8), 2007-2023CODEN: JCISD8; ISSN:1549-9596. (American Chemical Society)We present TRAVIS ("TRajectory Analyzer and VISualizer"), a free program package for analyzing and visualizing Monte Carlo and mol. dynamics trajectories. The aim of TRAVIS is to collect as many analyses as possible in one program, creating a powerful tool and making it unnecessary to use many different programs for evaluating simulations. This should greatly rationalize and simplify the work-flow of analyzing trajectories. TRAVIS is written in C++, open-source free-ware and licensed under the terms of the GNU General Public License (GPLv3). It is easy to install (platform independent, no external libraries) and easy to use. In this article, we present some of the algorithms that are implemented in TRAVIS - many of them widely known for a long time, but some of them also to appear in literature for the first time. All shown analyses only require a std. MD trajectory as input data.
- 33Brehm, M.; Weber, H.; Thomas, M.; Hollòczki, O.; Kirchner, B. Domain analysis in nanostructured liquids: A post-molecular dynamics study at the example of ionic liquids. ChemPhysChem 2015, 16, 3271– 3277, DOI: 10.1002/cphc.201500471Google Scholar33https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2MXhsFajtLjJ&md5=4723ce7f422ad2886440a6146c7f55afDomain Analysis in Nanostructured Liquids: A Post-Molecular Dynamics Study at the Example of Ionic LiquidsBrehm, Martin; Weber, Henry; Thomas, Martin; Holloczki, Oldamur; Kirchner, BarbaraChemPhysChem (2015), 16 (15), 3271-3277CODEN: CPCHFT; ISSN:1439-4235. (Wiley-VCH Verlag GmbH & Co. KGaA)In the present computational work, we develop a new tool for our trajectory anal. program TRAVIS to analyze the well-known behavior of liqs. to sep. into microphases. The dissection of the liq. into domains of different subsets, for example, in the case of fluorinated ionic liqs. with anionic and cationic head groups (forming together the polar domain), fluorous, and alkyl subsets is followed by radical Voronoi tessellation. This leads to useful av. quantities of the subset neighbor count, i.e., the domain count that gives the amt. of particular domains in the liq., the domain vol. and surface, as well as the isoperimetric quotient, which provides a measure of the deviation of the domains from a spherical shape. Thus, the newly implemented method allows anal. of the domains in terms of their nos. and shapes on a qual. and also quant. basis. It is a simple, direct, and automated anal. that does not need evaluation of the structure beforehand in terms of, for example, first solvent shell min. In the microheterogeneous ionic liqs. that we used as examples, the polar subsets always form a single domain in all investigated liqs. Although the fluorous side chains are also more or less connected in one or, maximally, two domains, the alkyl phases are dispersed.
- 34Luzar, A.; Chandler, D. Effect of environment on hydrogen bond dynamics in liquid water. Phys. Rev. Lett. 1996, 76, 928– 931, DOI: 10.1103/PhysRevLett.76.928Google Scholar34https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK28Xotl2iuw%253D%253D&md5=5f3b5cab4da60239ec27fae8d34b9a21Effect of environment on hydrogen bond dynamics in liquid waterLuzar, Alenka; Chandler, DavidPhysical Review Letters (1996), 76 (6), 928-31CODEN: PRLTAO; ISSN:0031-9007. (American Physical Society)In a series of mol.-dynamics calcns., we simulate the dynamics of forming and breaking a hydrogen bond in liq. water at room temp. We show that this dynamics is clearly nonexponential, yet virtually uncorrelated with the fluctuations of neighboring bonds.
- 35Luzar, A.; Chandler, D. Hydrogen-bond kinetics in liquid water. Nature 1996, 379, 55– 57, DOI: 10.1038/379055a0Google Scholar35https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK28XivVKhtQ%253D%253D&md5=282f1fbda5394e9ebd896db3892a005cHydrogen-bond kinetics in liquid waterLuzar, Alenka; Chandler, DavidNature (London) (1996), 379 (6560), 55-7CODEN: NATUAS; ISSN:0028-0836. (Macmillan Magazines)Hydrogen bonds play a crucial role in the behavior of water; their spatial patterns and fluctuations characterize the structure and dynamics of the liq. The processes of breaking and making hydrogen bonds in the condensed phase can be probed indirectly by a variety of exptl. techniques, and more quant. information can be obtained from computer simulations. In particular, simulations have revealed that on long timescales the relaxation behavior of hydrogen bonds in liq. water exhibit non-exponential kinetics, suggesting that bond making and breaking are not simple processes characterized by well defined rate consts. Here we show that these kinetics can be understood in terms of an interplay between diffusion and hydrogen-bond dynamics. In our model, which can be extended to other hydrogen-bonded liqs., diffusion governs whether a specific pair of water mols. are near neighbors, and hydrogen bonds between such pairs form and persist at random with av. lifetimes detd. by rate consts. for bond making and breaking.
- 36Luzar, A. Resolving the hydrogen bond dynamics conundrum. J. Chem. Phys. 2000, 113, 10663– 10675, DOI: 10.1063/1.1320826Google Scholar36https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3cXosFantb4%253D&md5=bc078176c816b6143f48bac25241748cResolving the hydrogen bond dynamics conundrumLuzar, AlenkaJournal of Chemical Physics (2000), 113 (23), 10663-10675CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)This paper analyzes dynamic properties of hydrogen bonds in liq. water. We use mol. dynamics simulation to calc. different probability densities that govern the time evolution of the formation and rupture of hydrogen bonds. We provide anal. connections between these functions. Excellent agreement with our simulation results is obsd. We prove transition state theory rate const. to be identical to the inverse of the assocd. mean first passage time (hydrogen bond lifetime). Hence, the anal. establishes its Arrhenius temp. dependence. We give the explicit relation between reactive flux correlation function for the relaxation dynamics of hydrogen bonds, and their first passage time probability densities. All the different observations in the existing literature, assocd. with various ests. of hydrogen bonding times in liq. water that are affected (or not affected) by particular bond criteria, as well as by different definitions of hydrogen bond lifetimes applied in simulation, can be easily reconciled within the framework of reactive flux correlation function approach.
- 37Gehrke, S.; von Domaros, M.; Clark, R.; Hollòczki, O.; Brehm, M.; Welton, T.; Luzar, A.; Kirchner, B. Structure and lifetimes in ionic liquids and their mixtures. Faraday Discuss. 2018, 206, 219– 245, DOI: 10.1039/C7FD00166EGoogle Scholar37https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXhtFCltbjM&md5=156128b1e8a1f4c6214e7e3b193e5719Structure and lifetimes in ionic liquids and their mixturesGehrke, Sascha; von Domaros, Michael; Clark, Ryan; Holloczki, Oldamur; Brehm, Martin; Welton, Tom; Luzar, Alenka; Kirchner, BarbaraFaraday Discussions (2018), 206 (Ionic Liquids: From Fundamental Properties to Practical Applications), 219-245CODEN: FDISE6; ISSN:1359-6640. (Royal Society of Chemistry)With the aid of mol. dynamics simulations, we study the structure and dynamics of different ionic liq. systems, with focus on hydrogen bond, ion pair and ion cage formation. To do so, we report radial distribution functions, their no. integrals, and various time-correlation functions, from which we ext. well-defined lifetimes by means of the reactive flux formalism. We explore the influence of polarizable force fields vs. non-polarizable ones with downscaled charges (±0.8) for the example of 1-butyl-3-methylimidazolium bromide. Furthermore, we use 1-butyl-3-methylimidazolium trifluoromethanesulfonate to investigate the impact of temp. and mixing with water as well as with the chloride ionic liq. Smaller coordination nos., larger distances, and tremendously accelerated dynamics are obsd. when the polarizable force field is applied. The same trends are found with increasing temp. Adding water decreases the ion-ion coordination nos. whereas the water-ion and water-water coordination is enhanced. A domain anal. reveals that the nonpolar parts of the ions are dispersed and when more water is added the water clusters increase in size. The dynamics accelerate in general upon addn. of water. In the ionic liq. mixt., the coordination no. around the cation changes between the two anions, but the no. integrals of the cation around the anions remain const. and the dynamics slow down with increasing content of the chloride ionic liq.
- 38Brehm, M.; Weber, H.; Thomas, M.; Hollòczki, O.; Kirchner, B. Domain analysis in nanostructured liquids: a post-molecular dynamics study at the example of ionic liquids. ChemPhysChem 2015, 16, 3271– 3277, DOI: 10.1002/cphc.201500471Google Scholar38https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2MXhsFajtLjJ&md5=4723ce7f422ad2886440a6146c7f55afDomain Analysis in Nanostructured Liquids: A Post-Molecular Dynamics Study at the Example of Ionic LiquidsBrehm, Martin; Weber, Henry; Thomas, Martin; Holloczki, Oldamur; Kirchner, BarbaraChemPhysChem (2015), 16 (15), 3271-3277CODEN: CPCHFT; ISSN:1439-4235. (Wiley-VCH Verlag GmbH & Co. KGaA)In the present computational work, we develop a new tool for our trajectory anal. program TRAVIS to analyze the well-known behavior of liqs. to sep. into microphases. The dissection of the liq. into domains of different subsets, for example, in the case of fluorinated ionic liqs. with anionic and cationic head groups (forming together the polar domain), fluorous, and alkyl subsets is followed by radical Voronoi tessellation. This leads to useful av. quantities of the subset neighbor count, i.e., the domain count that gives the amt. of particular domains in the liq., the domain vol. and surface, as well as the isoperimetric quotient, which provides a measure of the deviation of the domains from a spherical shape. Thus, the newly implemented method allows anal. of the domains in terms of their nos. and shapes on a qual. and also quant. basis. It is a simple, direct, and automated anal. that does not need evaluation of the structure beforehand in terms of, for example, first solvent shell min. In the microheterogeneous ionic liqs. that we used as examples, the polar subsets always form a single domain in all investigated liqs. Although the fluorous side chains are also more or less connected in one or, maximally, two domains, the alkyl phases are dispersed.
- 39Elfgen, R.; Hollòczki, O.; Kirchner, B. A molecular level understanding of template effects in ionic liquids. Acc. Chem. Res. 2017, 50, 2949– 2957, DOI: 10.1021/acs.accounts.7b00436Google Scholar39https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXhsl2ms7%252FL&md5=11aa98fbdf91a41743fabb670bb036a7A Molecular Level Understanding of Template Effects in Ionic LiquidsElfgen, Roman; Holloczki, Oldamur; Kirchner, BarbaraAccounts of Chemical Research (2017), 50 (12), 2949-2957CODEN: ACHRE4; ISSN:0001-4842. (American Chemical Society)The structure-directing or template effect has been invoked several times for ionic liqs. to explain the different outcome in material synthesis, namely, different scaffolds or geometrical arrangements with varying ionic liqs. It is obvious to assume that such an effect can originate from the most likely complex microstructure, being present within the ionic liq. itself. In that regard, ionic liqs. have already been shown to undergo a nanosegregation into polar and nonpolar phases, which is commonly known and denoted as microheterogeneity. In order to provide detailed insight on the mol. level and to understand the effects rising from this structuring, the authors performed mol. dynamics simulations on selected very simple model systems composed of 1-alkyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide, considering Et, Bu, hexyl, and octyl side chains attached to the cations, mixed with either n-dodecanol or n-butanol. By analyzing snapshots of the simulation boxes and calcg. spatial distribution functions, one can visualize that with increasing side chains, the systems show considerable nanosegregation into polar and nonpolar domains. Combined angular and distance distribution functions show that in case of the nanosegregating systems the side chains of the cations are preferentially arranged in a parallel fashion, which indicates a micelle-like structure for the ionic liqs. The alc. mols. participate in and are, therefore, influenced by this microheterogeneity. It can be shown that in the case of the short IL alkyl side chains, the self-aggregation of the nonpolar units of the alcs. is much stronger, while for the long chain cations, the nonpolar entities of the alcs. are most often connected to the nonpolar units of the ionic liqs. Using our domain anal. tool, we can quantify these observations by tracking the no., size, and shape of the polar and nonpolar entities present in the different investigated systems. The aforementioned combined angular-distance distribution functions reveal a structure-directing effect of the ionic liqs. on the alc. mols. within our simple model systems. The ionic liqs. act as template and order the alc. mols. according to their own structure, resulting in a parallel alignment of the alkyl side chains of the alcs. and ionic liq. cations, with both polar groups being at the same side. These observations show that the microheterogeneous structure of ionic liqs. can indeed be applied to order substrates with respect to each other or, for example, to catalysts in a predetd. fashion, opening new possibilities for explaining or enhancing selectivities of chem. reactions in ionic liqs.
- 40Macchieraldo, R.; Esser, L.; Elfgen, R.; Voepel, P.; Zahn, S.; Smarsly, B. M.; Kirchner, B. Hydrophilic ionic liquid mixtures of weakly and strongly coordinating anions with and without water. ACS Omega 2018, 3, 8567– 8582, DOI: 10.1021/acsomega.8b00995Google Scholar40https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXhsVajsrrF&md5=232946b7312d7bca99099aea1a5477b6Hydrophilic Ionic Liquid Mixtures of Weakly and Strongly Coordinating Anions with and without WaterMacchieraldo, Roberto; Esser, Lars; Elfgen, Roman; Voepel, Pascal; Zahn, Stefan; Smarsly, Bernd M.; Kirchner, BarbaraACS Omega (2018), 3 (8), 8567-8582CODEN: ACSODF; ISSN:2470-1343. (American Chemical Society)With the aid of ab initio mol. dynamics simulations, we investigate an ionic liq. mixt. composed of three components 1-butyl-3-methylimidazolium [C4C1Im]+, tetrafluoroborate [BF4]- and chloride [Cl]- without and with water. In the pure IL mixt., we observe an already complex network of interactions between cations and anions, adding water to the system even extends the complexity. Obsd. no. integrals show that the coordination no. between cations and anions is reduced in the system with water compared to the pure system. Further studies show, that the Coulombic network of the strongly coordinating anion [Cl]- is disturbed by water, while the weakly coordinating anion [BF4]- is not. These observations can also be confirmed by the Voronoi polyhedra anal., which shows that the polar network of microheterogeneous IL collapses by the introduction of water. Hydrogen-acceptor interactions revealed, that the [Cl]- anions are transferred from being situated in the IL to the water continuum while [BF4]-is almost unperturbed, these effects mainly influence the interplay of the ionic liq. network.
- 41Ray, P.; Vogl, T.; Balducci, A.; Kirchner, B. Structural investigations on lithium-doped protic and aprotic ionic liquids. J. Phys. Chem. B 2017, 121, 5279– 5292, DOI: 10.1021/acs.jpcb.7b02636Google Scholar41https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXntFeiuro%253D&md5=246661ee345211fbff9b320750518a8dStructural Investigations on Lithium-Doped Protic and Aprotic Ionic LiquidsRay, Promit; Vogl, Thomas; Balducci, Andrea; Kirchner, BarbaraJournal of Physical Chemistry B (2017), 121 (20), 5279-5292CODEN: JPCBFK; ISSN:1520-5207. (American Chemical Society)Solns. of lithium bis(trifluoromethanesulfonyl)imide (LiNTf2), in four different [NTf2]--based ionic liqs., are extensively investigated as potential electrolytes for lithium-ion batteries. Solvation of the [Li]+ ions in the ionic liqs. and its impact on their physicochem. properties are studied herein with the aid of mol. dynamics simulations. The cationic components of the investigated liqs. were systematically varied so as to individually evaluate effects of specific structural changes; increase in ring size, the addn. of an alkyl chain and absence of an acidic proton, on the solvation and mobility of the [Li]+ cations. The studied cations also allow for a direct comparison between solns. of [Li]+ salt in protic and aprotic ionic liqs. Emphasis is laid on elucidating the interactions between the [Li]+ and [NTf2]- ions revealing slightly higher coordination nos. for the aprotic solvent, benchmarked against exptl. measurements. The study suggests that the ionic liqs. largely retain their structure upon salt addn., with interactions within the liqs. only slightly perturbed. The rattling motion of the [Li]+ cations within cages formed by the surrounding [NTf2]- anions is examd. by the anal. of [Li]+ autocorrelation functions. Overall, the solvation mechanism of [Li]+ salt, within the hydrogen-bonded network of the ionic liqs., is detailed from classical and ab initio mol. dynamics simulations.
- 42Kirchner, B.; Hutter, J.; Kuo, I.-F. W.; Mundy, C. J. Hydrophobic hydratation from Car-Parrinello simulations. Int. J. Mod. Phys. B 2004, 18, 1951– 1962, DOI: 10.1142/S0217979204007241Google Scholar42https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2cXmvVaru7Y%253D&md5=0bd7f67008e20970e686112f086857baHydrophobic hydration from Car-Parrinello simulationsKirchner, Barbara; Hutter, Juerg; Kuo, I.-Feng W.; Mundy, Christopher J.International Journal of Modern Physics B: Condensed Matter Physics, Statistical Physics, Applied Physics (2004), 18 (14), 1951-1962CODEN: IJPBEV; ISSN:0217-9792. (World Scientific Publishing Co. Pte. Ltd.)In this work we investigate the fast anomalous diffusion of hydrogen mols. in water using Car-Parrinello mol. dynamics simulations. We employ Voronoi polyhedra anal. to distinguish between void diffusion and void hopping. Our results indicate that a combination of both mechanism is sufficient to explain anomalous diffusion. Furthermore, we investigate the geometry and the electronic structure of the first solvation shell.
- 43Ludwig, R. Water: From Clusters to the Bulk. Angew. Chem., Int. Ed. 2001, 40, 1808– 1827, DOI: 10.1002/1521-3773(20010518)40:10<1808::AID-ANIE1808>3.0.CO;2-1Google Scholar43https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3MXjvFyqtrw%253D&md5=6194b6257b7ef8884216627ec9d53132Water: from clusters to the bulkLudwig, RalfAngewandte Chemie, International Edition (2001), 40 (10), 1808-1827CODEN: ACIEF5; ISSN:1433-7851. (Wiley-VCH Verlag GmbH)A review with 186 refs. Water is of fundamental importance for human life and plays an important role in many biol. and chem. systems. Although water is the most abundant compd. on earth, it is definitely not a simple liq. It possesses strongly polar hydrogen bonds which are responsible for a striking set of anomalous phys. and chem. properties. For more than a century the combined importance and peculiarity of water inspired scientists to construct conceptual models, which in themselves reproduce the obsd. behavior of the liq. The exploration of structural and binding properties of small water complexes provides a key for understanding bulk water in its liq. and solid phase and for understanding solvation phenomena. Modern ab initio quantum chem. methods and high-resoln. spectroscopy methods have been extremely successful in describing such structures. Cluster models for liq. water try to mimic the transition from these clusters to bulk water. The important question is: What cluster properties are required to describe liq.-phase behavior.
- 44Huber, H.; Kirchner, B.; Searles, D. J. Is there an iceberg effect in the water/DMSO mixture? Some information from computational chemistry. J. Mol. Liq. 2002, 98-99, 71– 77, DOI: 10.1016/S0167-7322(01)00310-5Google Scholar44https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD38Xis1ehtbg%253D&md5=479106f6b2917627b4060648be17ac40Is There an Iceberg Effect in the Water/DMSO Mixture? Some Information from Computational Chemistry.Huber, Hanspeter; Kirchner, Barbara; Searles, Debra J.Journal of Molecular Liquids (2002), 98-99 (), 71-77CODEN: JMLIDT; ISSN:0167-7322. (Elsevier Science S.A.)The deuteron quadrupole coupling const. of heavy water in a mixt. with DMSO measured by B. C. Gordalla and M. D. Zeidler (1986)shows an unexpected behavior. A completely different evaluation from computational methods with independent approxns. does not reproduce this finding. A crit. discussion of all approxns. and uncertainties of the computational method is given using the limited results available for this system and results from studies on pure water.
- 45Galamba, N. Water’s structure around hydrophobic solutes and the iceberg model. J. Phys. Chem. B 2013, 117, 2153– 2159, DOI: 10.1021/jp310649nGoogle Scholar45https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3sXhsFGrtr0%253D&md5=4413978b1e29db158e633af3d0a40e03Water's Structure around Hydrophobic Solutes and the Iceberg ModelGalamba, N.Journal of Physical Chemistry B (2013), 117 (7), 2153-2159CODEN: JPCBFK; ISSN:1520-5207. (American Chemical Society)The structure of water in the hydration shells of small hydrophobic solutes was investigated through mol. dynamics. The results show that a subset of water mols. in the first hydration shell of a nonpolar solute have a significantly enhanced tetrahedrality and a slightly larger no. of hydrogen bonds, relative to the mols. in water at room temp., consistent with the exptl. obsd. neg. excess entropy and increased heat capacity of hydrophobic solns. at room temp. This ordering results from the rearrangement of a small no. of water mols. near the nonpolar solutes that occupy one to two vertices of the enhanced water tetrahedra. Although this structuring is not nearly like that often assocd. with a literal interpretation of the term "iceberg" in the Frank and Evans iceberg model, it does support a moderate interpretation of this model. Thus, the tetrahedral orientational order of this ensemble of water mols. is comparable to that of liq. water at ∼10°C, although not accompanied by the small contraction of the O-O distance obsd. in cold water. The structural changes of water in the vicinity of small nonpolar solutes cannot be inferred from the water radial distribution functions, explaining why this increased ordering is not obsd. through neutron diffraction expts. The present results restore a mol. view where the slower translational and reorientational dynamics of water near hydrophobic groups has a structural equiv. resembling water at low temps.
- 46Chowdhary, J.; Ladanyi, B. M. Hydrogen bond dynamics at the water/hydrocarbon interface. J. Phys. Chem. B 2009, 113, 4045– 4053, DOI: 10.1021/jp8061509Google Scholar46https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1cXhtl2ltLfN&md5=8f09abe1c2d114eded38bcc622dd70e7Hydrogen Bond Dynamics at the Water/Hydrocarbon InterfaceChowdhary, Janamejaya; Ladanyi, Branka M.Journal of Physical Chemistry B (2009), 113 (13), 4045-4053CODEN: JPCBFK; ISSN:1520-6106. (American Chemical Society)The dynamics of hydrogen bond formation and breakage for water in the vicinity of water/hydrocarbon liq. interfaces is studied using mol. dynamics simulations. Several liq. alkanes are considered as the hydrocarbon phase to det. the effects of their chain length and extent of branching on the properties of the adjacent water phase. In addn. to defining the interface location in terms of the lab.-frame d. profiles, the effects of interfacial fluctuations are considered by locating the interface in terms of the proximity of the mols. of the other phase. The hydrogen bond dynamics of interfacial water is weakly influenced by the identity of the hydrocarbon phase and by capillary waves. In addn. to calcg. hydrogen bond time correlations, the authors examine how the hydrogen bond dynamics depend on local coordination and det. the extent of cooperativity in the population relaxation of the hydrogen bonds that a given mol. participates in. The contributions of translational diffusion and reorientation of mol. O-H bonds to the mechanism of hydrogen bond breakage and reformation are studied. In previous work, rotation of the principal axes of water is anisotropic at the interface and depends on the initial orientation of the mol. relative to the interface. Here, the authors extend this anal. to the reorientation of the O-H vector and to hydrogen bond time correlation. Hydrogen bond dynamics are also sensitive to the initial orientation of the mols. participating in the hydrogen bond.
- 47Li, Q.; Wu, G.; Yu, Z. The role of methyl groups in the formation of hydrogen bond in DMSO-methanol mixtures. J. Am. Chem. Soc. 2006, 128, 1438– 1439, DOI: 10.1021/ja0569149Google Scholar47https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD28XjsFeltA%253D%253D&md5=5ece10fd302fa06d97bed75146bf6fb0The Role of Methyl Groups in the Formation of Hydrogen Bond in DMSO-Methanol MixturesLi, Qingzhong; Wu, Guoshi; Yu, ZhiwuJournal of the American Chemical Society (2006), 128 (5), 1438-1439CODEN: JACSAT; ISSN:0002-7863. (American Chemical Society)When examg. the formation energetics of a hydrogen-bonded complex R-X-H···Y-R', focus has been almost always on the atoms directly involved, namely the atoms X, Y, and H. Little attention has been paid to the effects of the secondary alkyl groups R and R'. Taking DMSO-methanol binary system as an example, we have studied the roles of the alkyl groups in stabilizing the hydrogen bonds by employing FTIR and NMR techniques and quantum chem. calcns. We found that Me groups play different roles in response to the hydrogen-bonding interactions. The Me groups of DMSO are electron-donating, whereas that of methanol is electron-withdrawing, both making pos. contributions. The findings reveal non-negligible effects of secondary alkyl groups in hydrogen bonding interaction and may shed light on the understanding of other more complicated hydrogen-bonded systems in chem. and biol. systems.
- 48Yao, H.; Ke, H.; Zhang, X.; Pan, S.-J.; Li, M.-S.; Yang, L.-P.; Schreckenbach, G.; Jiang, W. Molecular recognition of hydrophilic molecules in water by combining the hydrophobic effect with hydrogen bonding. J. Am. Chem. Soc. 2018, 140, 13466– 13477, DOI: 10.1021/jacs.8b09157Google Scholar48https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXhslKqt7rE&md5=bb3a11ceab70b41baf08ccdc55b71749Molecular Recognition of Hydrophilic Molecules in Water by Combining the Hydrophobic Effect with Hydrogen BondingYao, Huan; Ke, Hua; Zhang, Xiaobin; Pan, San-Jiang; Li, Ming-Shuang; Yang, Liu-Pan; Schreckenbach, Georg; Jiang, WeiJournal of the American Chemical Society (2018), 140 (41), 13466-13477CODEN: JACSAT; ISSN:0002-7863. (American Chemical Society)During the last half a century, great achievements have been made in mol. recognition in parallel with the invention of numerous synthetic receptors. However, the selective recognition of hydrophilic mols. in water remains a generally accepted challenge in supramol. chem. but is commonplace in nature. In an earlier Communication [Huang et al. J. Am. Chem. Soc. 2016, 138, 14550], we reported a pair of endo-functionalized mol. tubes that surprisingly prefer highly hydrophilic mols. over hydrophobic mols. of a similar size and shape. The hydrophobic effect and hydrogen bonding were proposed to be responsible, but their exact roles were not fully elucidated. In this Article, we present a thorough study on the binding behavior of these mol. tubes toward 44 hydrophilic mols. in water. Principal component anal. reveals that the binding strength is weakly correlated to the hydrophobicity, vol., surface area, and dipole moment of guests. Furthermore, mol. dynamics simulations show the hydrophobic effect through releasing the poorly hydrogen-bonded cavity water contributes to the binding of all the hydrophilic mols., while hydrogen bonding differentiates these mols. and is thus the key to achieve a high selectivity toward certain hydrophilic mols. over other mols. with a similar size and shape. Therefore, a good guest for these mol. tubes should meet the following criteria: the hydrogen-bonding sites should be complementary, and the mol. vol. should be large enough to expel all the cavity water but not too large to cause steric hindrance. This rule of thumb may also be used to design a selective receptor for certain hydrophilic mols. Following these guidelines, a "best-fit" guest was found for the syn-configured mol. tube with a binding const. as high as 106 M-1.
- 49Yamazoe, K.; Higaki, Y.; Inutsuka, Y.; Miyawaki, J.; Cui, Y.-T.; Takahara, A.; Harada, Y. Enhancement of the hydrogen-bonding network of water confined in a polyelectrolyte brush. Langmuir 2017, 33, 3954– 3959, DOI: 10.1021/acs.langmuir.7b00243Google Scholar49https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXltlWhs7o%253D&md5=9b6e57791c57cdd8e371998949ea8b52Enhancement of the Hydrogen-Bonding Network of Water Confined in a Polyelectrolyte BrushYamazoe, Kosuke; Higaki, Yuji; Inutsuka, Yoshihiro; Miyawaki, Jun; Cui, Yi-Tao; Takahara, Atsushi; Harada, YoshihisaLangmuir (2017), 33 (16), 3954-3959CODEN: LANGD5; ISSN:0743-7463. (American Chemical Society)Water existing in the vicinity of polyelectrolytes exhibits unique structural properties, which demonstrate key roles in chem., biol., and geoscience. In this study, X-ray absorption and emission spectroscopy was employed to observe the local hydrogen-bonding structure of water confined in a charged polyelectrolyte brush. Even at room temp., a majority of the water mols. confined in the polyelectrolyte brush exhibited one type of hydrogen-bonding configuration: a slightly distorted, albeit ordered configuration. The findings from this study provide new insights in terms of correlation between the function and local structure of water at the interface of biol. materials under physiol. conditions.
- 50Wallen, S. L.; Palmer, B. J.; Garrett, B. C.; Yonker, C. R. Density and temperature effects on the hydrogen bond structure of liquid methanol. J. Phys. Chem. 1996, 100, 3959– 3964, DOI: 10.1021/jp9524082Google Scholar50https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK28XhtVKrsrs%253D&md5=2414d722bd2b4ab2e4c78644b6a62b8cDensity and Temperature Effects on the Hydrogen Bond Structure of Liquid MethanolWallen, S. L.; Palmer, B. J.; Garrett, B. C.; Yonker, C. R.Journal of Physical Chemistry (1996), 100 (10), 3959-64CODEN: JPCHAX; ISSN:0022-3654. (American Chemical Society)The H bond structure of liq. MeOH was studied as a function of pressure and temp. up to 2.8 kbar and from 297 to 413 K. Chem. shifts of the CH3 and OH groups were monitored throughout this pressure and temp. regime, and the chem. shift difference between these two groups was used to describe changes of the H bond network in MeOH. The H bond equil. was studied using mol. dynamics simulations and a phenomenol. model describing clustering in liq. MeOH. Results are presented concerning the size and distribution of H-bonded clusters in MeOH as a function of pressure and temp. The extent of H bonding decreases upon an increase in temp. The results for pressure are equivocal, the phenomenol. model suggests that H bonding decreases with increasing pressure, which supports earlier interpretations regarding the measured self-diffusion coeffs. in deuterated MeOH as a function of pressure. The mol. dynamics simulations, however, show an increase in H bonding with increasing pressure.
Supporting Information
Supporting Information
ARTICLE SECTIONSThe Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpcb.9b00839.
1H NMR spectra; radial distribution function; RDFs of the interactions of MeOH(O) molecules (PDF)
Terms & Conditions
Most electronic Supporting Information files are available without a subscription to ACS Web Editions. Such files may be downloaded by article for research use (if there is a public use license linked to the relevant article, that license may permit other uses). Permission may be obtained from ACS for other uses through requests via the RightsLink permission system: http://pubs.acs.org/page/copyright/permissions.html.