Next Article in Journal
Phonon Sideband Analysis and Near-Infrared Emission in Heavy Metal Oxide Glasses
Previous Article in Journal
Corrosion Behaviour of 316L Stainless Steel in CNTs–Water Nanofluid: Effect of Temperature
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Correction

Correction: Quevillon, M.J., et al. Charge Transport and Phase Behavior of Imidazolium-Based Ionic Liquid Crystals from Fully Atomistic Simulations. Materials 2018, 11, 64

by
Michael J. Quevillon
and
Jonathan K. Whitmer
*
Department of Chemical and Biomolecular Engineering, University of Notre Dame, Notre Dame, IN 46556, USA
*
Author to whom correspondence should be addressed.
Materials 2021, 14(1), 120; https://doi.org/10.3390/ma14010120
Submission received: 2 May 2020 / Accepted: 8 July 2020 / Published: 30 December 2020
(This article belongs to the Section Materials Simulation and Design)

1. Introduction

In our publication [1], we performed classical molecular dynamics simulations of a system of ionic liquid crystals using a transferable OPLS-AA force field [2] specifically developed for ionic liquids. In performing subsequent research, we became aware of a paper [3] performed by the same authors as reference [2], which updates some of the force field parameters for molecular dynamics simulations. However, the functional form of the torsional angles in the newer paper [3] was listed in a different form than was used in the original paper of the force field [2]. This lead us to review the force field that we had implemented in our study [1] and to determine if this had had significant quantitative or qualitative effects.
With corrected functional form in the subsequent publication [3], the authors wish to make corrections to our publication [1], as well as replace the affected figures. All references to the original force field have been updated to include the subsequent publication.
Additionally, we correct an experimental comparison within Table 1, which used data on [C n Im][NO 3 ] rather than [C n MIm][NO 3 ].

2. Text Corrections

The functional form of the OPLS-AA force field is corrected below. Replace
E = i bonds k r , i r i r o , i 2 + i angles k θ , i θ i θ o , i 2 + i torsions 1 2 V 1 , i ( 1 + cos ϕ i ) + V 2 , i ( 1 + cos 2 ϕ i ) + V 3 , i ( 1 + cos 3 ϕ i ) + V 4 , i ( 1 + cos 4 ϕ i ) + i atoms j > i q i q j e 2 r i j + 4 ε i j σ i j r i j 12 σ i j r i j 6
with
E = i bonds 1 2 k r , i r i r o , i 2 + i angles 1 2 k θ , i θ i θ o , i 2 + i torsions 1 2 V 1 , i ( 1 + cos ϕ i ) + V 2 , i ( 1 cos 2 ϕ i ) + V 3 , i ( 1 + cos 3 ϕ i ) + V 4 , i ( 1 cos 4 ϕ i ) + i atoms j > i q i q j e 2 r i j + 4 ε i j σ i j r i j 12 σ i j r i j 6
In addition, with the corrected force field, the nitrate anions are no longer rigid [3], so remove the following sentence from line 110:
  • Nitrate anions are considered to be rigid bodies, while the imidazolium-based species are able to fluctuate more freely.
With the adjusted force field [3], there were certain aspects of the simulation protocol that were slightly adjusted. Replace the following sentences, starting at line 126,
  • To probe a range of temperatures, the systems were set up in a proposed crystal configuration (see Figure 2) and equilibrated for 20 ns at 300 K , before increasing the temperature by discrete steps of 25 K until either the internal temperature is 600 K or the material melts into an isotropic liquid. The latter condition was utilized for the n = 12 simulations, as the semi-isotropic box exhibited instability within the “isotropic” phase.
with
  • To probe a range of temperatures, the systems were set up in a proposed crystal configuration (see Figure 2) and equilibrated for 20 ns at 300 K , before increasing the temperature by discrete steps of 25 K up to 600 K .
Add the following comment on the low-temperature crystal configurations at line 158,
  • Despite starting from a parallel initial configuration, these systems adopted a tilted, chevron-like ordering upon equilibration. This tilting is consistent with other simulations [4], but the alternating directions of tilt is an effect of our chosen initial configuration. The simplest way for the system to transform into a tilted configuration is by moving the central ionic layer in one direction as a whole, which induces the chevron-like ordering observed here.
Replace the beginning of the sentence starting at line 158,
  • Though we will show later this state exhibits properties suggesting it is a distinct phase …
with
  • Though we will show later this state exhibits structural and dynamic properties suggesting it is a distinct phase …
One important difference between the prior results and those with the corrected force field is addressed below.
Replace the following sentences, starting at line 198,
  • Interestingly, this graph demonstrates a non-monotonicity in transition temperature with alkyl tail length, with a broadened smectic-A phase also appearing for n = 14 and n = 16 . Though a more thorough investigation is required to understand the exact cause for this behavior, we hypothesize that it might involve a trade-off between crystal orderings dominated by electrostatic and van der Waals interactions, as inter-layer Coulombic interactions should be felt more strongly between layers when the aliphatic tails are short, while van der Waals interactions would have increased influence as the tails grow.
with
  • By increasing the length of the alkyl side chain, the relative influence of the inter-layer Coulombic interactions and the van der Waals interactions stabilizes the smectic-A phases, thereby broadening the range of temperatures where the smectic-A is observed.
The long-range ordering behavior somewhat changed with the corrected force field [3], which required altered commentary. In line 266, replace
  • As with the mean-squared displacement, an interesting change occurs before the crystalline–smectic transition, as the charged moieties begin to lose their long-range order, while the longer molecules retain more order.
with
  • As with the mean-squared displacement, an interesting change occurs before the crystalline–smectic transition, where the charged moieties (both [NO3] and [C14MIm]+) begin to exhibit a more smoothed radial distribution function than the aliphatic tails.

3. Figure Corrections

Replace
Figure 3. Observed phases in the ionic liquid crystal material, [C16MIm][NO3]. Typical configurations are presented from temperatures ranging between (a) 300 K and (d) 600 K . Frustrations within the crystal are visible in (a), as some aliphatic tails are bent within the layer and anions exhibit irregular placement. Intermediate temperatures illustrated in (b,c) depict the system on either side of the crystal–smectic phase transition. The isotropic liquid still exhibits segregation, but the layers have lost their long-range order, as evidenced by (d).
Figure 3. Observed phases in the ionic liquid crystal material, [C16MIm][NO3]. Typical configurations are presented from temperatures ranging between (a) 300 K and (d) 600 K . Frustrations within the crystal are visible in (a), as some aliphatic tails are bent within the layer and anions exhibit irregular placement. Intermediate temperatures illustrated in (b,c) depict the system on either side of the crystal–smectic phase transition. The isotropic liquid still exhibits segregation, but the layers have lost their long-range order, as evidenced by (d).
Materials 14 00120 g001
with
Figure 3. Observed phases in the ionic liquid crystal material, [C16MIm][NO3]. Typical configurations are presented from temperatures ranging between (a) 300 K and (d) 600 K . The tilted layers in (a) are consistent with other models [4], though the alternating direction of tilt is most likely due to our simulations’ proposed initial configurations. Frustrations within the crystal are visible in (a), as some aliphatic tails are bent within the layer and anions exhibit irregular placement. Intermediate temperatures illustrated in (b,c) depict the system on either side of the crystal–smectic phase transition. The isotropic liquid still exhibits segregation, but the layers have lost their long-range order, as evidenced by (d).
Figure 3. Observed phases in the ionic liquid crystal material, [C16MIm][NO3]. Typical configurations are presented from temperatures ranging between (a) 300 K and (d) 600 K . The tilted layers in (a) are consistent with other models [4], though the alternating direction of tilt is most likely due to our simulations’ proposed initial configurations. Frustrations within the crystal are visible in (a), as some aliphatic tails are bent within the layer and anions exhibit irregular placement. Intermediate temperatures illustrated in (b,c) depict the system on either side of the crystal–smectic phase transition. The isotropic liquid still exhibits segregation, but the layers have lost their long-range order, as evidenced by (d).
Materials 14 00120 g002
Replace
Figure 4. Intrinsic volume per ion pair as a function of temperature for three different alkyl chain lengths. The most dramatic jump in volume is between the crystalline and smectic phases. The rightmost point on each line is where the material melted into an isotropic liquid. Transition temperatures are summarized in Table 1.
Figure 4. Intrinsic volume per ion pair as a function of temperature for three different alkyl chain lengths. The most dramatic jump in volume is between the crystalline and smectic phases. The rightmost point on each line is where the material melted into an isotropic liquid. Transition temperatures are summarized in Table 1.
Materials 14 00120 g003
with
Figure 4. Intrinsic volume per ion pair as a function of temperature for three different alkyl chain lengths. The most dramatic jump in volume is between the crystalline and smectic phases. The transition between the smectic and isotropic phases is not apparent in this quantity. Transition temperatures are summarized in Table 1.
Figure 4. Intrinsic volume per ion pair as a function of temperature for three different alkyl chain lengths. The most dramatic jump in volume is between the crystalline and smectic phases. The transition between the smectic and isotropic phases is not apparent in this quantity. Transition temperatures are summarized in Table 1.
Materials 14 00120 g004
Replace
Figure 5. Intrinsic enthalpy as a function of temperature for three different alkyl chain lengths. The phase transitions are less pronounced in this system-wide property, but can be identified by the changes in slope of each line.
Figure 5. Intrinsic enthalpy as a function of temperature for three different alkyl chain lengths. The phase transitions are less pronounced in this system-wide property, but can be identified by the changes in slope of each line.
Materials 14 00120 g005
with
Figure 5. Intrinsic enthalpy as a function of temperature for three different alkyl chain lengths. The phase transitions are less pronounced in this system-wide property, but can be identified by the changes in slope of each line.
Figure 5. Intrinsic enthalpy as a function of temperature for three different alkyl chain lengths. The phase transitions are less pronounced in this system-wide property, but can be identified by the changes in slope of each line.
Materials 14 00120 g006
Replace
Figure 6. Nematic order parameter as a function of temperature for three different alkyl chain lengths, averaged over the last 200 ns of the trajectory. Two distinct phase transitions can be observed a discontinuous jumps in this order parameter. The error bars are larger on temperatures near transition points, as the systems fluctuate more there, thus introducing more variance into trajectory-averaged properties.
Figure 6. Nematic order parameter as a function of temperature for three different alkyl chain lengths, averaged over the last 200 ns of the trajectory. Two distinct phase transitions can be observed a discontinuous jumps in this order parameter. The error bars are larger on temperatures near transition points, as the systems fluctuate more there, thus introducing more variance into trajectory-averaged properties.
Materials 14 00120 g007
with
Figure 6. Nematic order parameter as a function of temperature for three different alkyl chain lengths, averaged over the last 200 ns of the trajectory. The layers are calculated separately and averaged, due to the opposite tilting directions in the crystal phase. Two distinct phase transitions can be observed as discontinuous jumps in this order parameter. The error bars are larger on temperatures near transition points, as the systems fluctuate more there, thus introducing more variance into trajectory-averaged properties.
Figure 6. Nematic order parameter as a function of temperature for three different alkyl chain lengths, averaged over the last 200 ns of the trajectory. The layers are calculated separately and averaged, due to the opposite tilting directions in the crystal phase. Two distinct phase transitions can be observed as discontinuous jumps in this order parameter. The error bars are larger on temperatures near transition points, as the systems fluctuate more there, thus introducing more variance into trajectory-averaged properties.
Materials 14 00120 g008
Replace
Table 1. Observed transition temperatures between the various phases for this series of ionic liquid crystals. Lee, et al. [5] did not observe any smectic phase for [C12MIm][NO3], but observed the isotropic phase above 318 K . As a function of alkyl tail length, the lower transition temperature exhibits a slight non-monotonicity for these species.
Table 1. Observed transition temperatures between the various phases for this series of ionic liquid crystals. Lee, et al. [5] did not observe any smectic phase for [C12MIm][NO3], but observed the isotropic phase above 318 K . As a function of alkyl tail length, the lower transition temperature exhibits a slight non-monotonicity for these species.
SpeciesTransitionTemperature, K
This WorkExperimental [5]Coarse-Grained [6]
C12Crystal–Smectic-A425--
Smectic-A–Isotropic475--
C14Crystal–Smectic-A400336-
Smectic-A–Isotropic550356-
C16Crystal–Smectic-A425339500
Smectic-A–Isotropic575404560
with
Table 1. Observed transition temperatures (upper bounds) between the various phases for this series of ionic liquid crystals. Guillet et al. [7] did not observe any smectic phase for [C12MIm][NO3], but observed the isotropic phase above 302 K . As a function of alkyl tail length, the trends of transition temperatures are consistent with experimental results [7].
Table 1. Observed transition temperatures (upper bounds) between the various phases for this series of ionic liquid crystals. Guillet et al. [7] did not observe any smectic phase for [C12MIm][NO3], but observed the isotropic phase above 302 K . As a function of alkyl tail length, the trends of transition temperatures are consistent with experimental results [7].
SpeciesTransitionTemperature, K
This WorkExperimental [7]Coarse-Grained [6]
C12Crystal–Smectic-A450--
Smectic-A–Isotropic500--
C14Crystal–Smectic-A475315-
Smectic-A–Isotropic550401-
C16Crystal–Smectic-A500320505
Smectic-A–Isotropic575457560
Replace
Figure 7. Mean-squared displacement of the species in [C16MIm][NO3] as a function of time at temperatures of (a) 300 K in the crystalline state, (b) 400 K just before the transition to smectic, (c) 425 K in the smectic-A phase, and (d) 600 K isotropic liquid. Reference slopes of 1, indicated as small triangles, help to identify the linear diffusive regime. The colors signify the dimension (red: x, blue: y, black: z); the solid lines are for the long cations and the dotted lines are for the nitrate anions. The most distinctive difference is the divergence of the z-direction from the x- and y-directions in the smectic-A phase, showing that transport within a layer is much simpler than through the layers. (Notes: (c) The x- and y-directions are nearly overlapping. (d) All directions nearly collapse, indicating an isotropic phase.)
Figure 7. Mean-squared displacement of the species in [C16MIm][NO3] as a function of time at temperatures of (a) 300 K in the crystalline state, (b) 400 K just before the transition to smectic, (c) 425 K in the smectic-A phase, and (d) 600 K isotropic liquid. Reference slopes of 1, indicated as small triangles, help to identify the linear diffusive regime. The colors signify the dimension (red: x, blue: y, black: z); the solid lines are for the long cations and the dotted lines are for the nitrate anions. The most distinctive difference is the divergence of the z-direction from the x- and y-directions in the smectic-A phase, showing that transport within a layer is much simpler than through the layers. (Notes: (c) The x- and y-directions are nearly overlapping. (d) All directions nearly collapse, indicating an isotropic phase.)
Materials 14 00120 g009
with
Figure 7. Mean-squared displacement of the species in [C14MIm][NO3] as a function of time at temperatures of (a) 400 K in the crystalline state, (b) 450 K just before the transition to smectic, (c) 500 K in the smectic-A phase, and (d) 600 K isotropic liquid. Reference slopes of 1, indicated as small triangles, help to identify the linear diffusive regime. The colors signify the dimension (red: x, blue: y, black: z); the solid lines are for the long cations and the dotted lines are for the nitrate anions. The most distinctive difference is the divergence of the z-direction from the x- and y-directions in the smectic-A phase, showing that transport within a layer is much simpler than through the layers. (Notes: (c) The x- and y-directions are nearly overlapping. (d) All directions nearly collapse, indicating an isotropic phase.)
Figure 7. Mean-squared displacement of the species in [C14MIm][NO3] as a function of time at temperatures of (a) 400 K in the crystalline state, (b) 450 K just before the transition to smectic, (c) 500 K in the smectic-A phase, and (d) 600 K isotropic liquid. Reference slopes of 1, indicated as small triangles, help to identify the linear diffusive regime. The colors signify the dimension (red: x, blue: y, black: z); the solid lines are for the long cations and the dotted lines are for the nitrate anions. The most distinctive difference is the divergence of the z-direction from the x- and y-directions in the smectic-A phase, showing that transport within a layer is much simpler than through the layers. (Notes: (c) The x- and y-directions are nearly overlapping. (d) All directions nearly collapse, indicating an isotropic phase.)
Materials 14 00120 g010
Replace
Figure 8. In-plane radial distribution function, g ( r ) of the species in [C16MIm][NO3], trajectory-averaged over the last 200 ns, at temperatures of (a) 300 K in the crystalline state, (b) 400 K just before the transition to the smectic-A state, (c) 425 K in the smectic-A phase, and (d) 600 K isotropic liquid. The solid lines signify distribution of a species relative to its own type (black: [C16MIm]+, red: [NO3], blue: imidazolyl ring) while the dotted line denotes placement of [NO3] ions in relation to imidazolyl rings. The first phase transition is characterized by the alkyl tails losing almost all order within the x y -plane. The phase transition to the isotropic liquid loses almost all structure, except for pairing of the opposite charges.
Figure 8. In-plane radial distribution function, g ( r ) of the species in [C16MIm][NO3], trajectory-averaged over the last 200 ns, at temperatures of (a) 300 K in the crystalline state, (b) 400 K just before the transition to the smectic-A state, (c) 425 K in the smectic-A phase, and (d) 600 K isotropic liquid. The solid lines signify distribution of a species relative to its own type (black: [C16MIm]+, red: [NO3], blue: imidazolyl ring) while the dotted line denotes placement of [NO3] ions in relation to imidazolyl rings. The first phase transition is characterized by the alkyl tails losing almost all order within the x y -plane. The phase transition to the isotropic liquid loses almost all structure, except for pairing of the opposite charges.
Materials 14 00120 g011
with
Figure 8. In-plane radial distribution function, g ( r ) of the species in [C14MIm][NO3], trajectory-averaged over the last 200 ns, at temperatures of (a) 400 K in the crystalline state, (b) 450 K just before the transition to the smectic-A state, (c) 500 K in the smectic-A phase, and (d) 600 K isotropic liquid. The solid lines signify the distribution of a species relative to its own type (black: [C14MIm]+, red: [NO3], blue: imidazolyl ring) while the dotted line denotes placement of [NO3] ions in relation to imidazolyl rings. From (a) to (b), a significant decrease in order of tail groups is observed, with the loss of secondary peaks, while the correlations also decay slightly more quickly as a function of distance. The first phase transition is characterized by the alkyl tails losing almost all order within the x y -plane. The phase transition to the isotropic liquid loses almost all structure, except for pairing of the opposite charges.
Figure 8. In-plane radial distribution function, g ( r ) of the species in [C14MIm][NO3], trajectory-averaged over the last 200 ns, at temperatures of (a) 400 K in the crystalline state, (b) 450 K just before the transition to the smectic-A state, (c) 500 K in the smectic-A phase, and (d) 600 K isotropic liquid. The solid lines signify the distribution of a species relative to its own type (black: [C14MIm]+, red: [NO3], blue: imidazolyl ring) while the dotted line denotes placement of [NO3] ions in relation to imidazolyl rings. From (a) to (b), a significant decrease in order of tail groups is observed, with the loss of secondary peaks, while the correlations also decay slightly more quickly as a function of distance. The first phase transition is characterized by the alkyl tails losing almost all order within the x y -plane. The phase transition to the isotropic liquid loses almost all structure, except for pairing of the opposite charges.
Materials 14 00120 g012
The authors regret any inconvenience these changes cause to readers, but seek to correct the record in the name of accurate and reproducible science.

Funding

M.J.Q. and J.K.W. greatly acknowledge computational resources provided by the Notre Dame Center for Research Computing (CRC), startup funds from the University of Notre Dame, and support from the Midwest Integrated Center for Computational Materials (MICCoM), a computational hub supported by the United States Department of Energy, Basic Energy Sciences division.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Quevillon, M.J.; Whitmer, J.K. Charge transport and phase behavior of imidazolium-based ionic liquid crystals from fully atomistic simulations. Materials 2018, 11, 64. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Sambasivarao, S.V.; Acevedo, O. Development of OPLS-AA Force Field Parameters for 68 Unique Ionic Liquids. J. Chem. Theory Comput. 2009, 5, 1038–1050. [Google Scholar] [CrossRef] [PubMed]
  3. Doherty, B.; Zhong, X.; Gathiaka, S.; Li, B.; Acevedo, O. Revisiting OPLS force field parameters for ionic liquid simulations. J. Chem. Theory Comput. 2017, 13, 6131–6145. [Google Scholar] [CrossRef] [PubMed]
  4. Cao, W.; Wang, Y.; Saielli, G. Metastable State during Melting and Solid–Solid Phase Transition of [C nMim][NO3](n = 4 − 12) Ionic Liquids by Molecular Dynamics Simulation. J. Phys. Chem. 2018, 122, 229–239. [Google Scholar] [CrossRef] [PubMed]
  5. Lee, C.K.; Huang, H.W.; Lin, I.J.B. Simple amphiphilic liquid crystalline 1-alkylimidazolium salts. A new solvent system providing a partially ordered environment. Chem. Commun. 2000, 68–74. [Google Scholar] [CrossRef]
  6. Saielli, G. MD simulation of the mesomorphic behaviour of 1-hexadecyl-3-methylimidazolium nitrate: Assessment of the performance of a coarse-grained force field. Soft Matter 2012, 8, 10279–10287. [Google Scholar] [CrossRef]
  7. Guillet, E.; Imbert, D.; Scopelliti, R.; Bünzli, J.C.G. Tuning the emission color of europium-containing ionic liquid-crystalline phases. Chem. Mater. 2004, 16, 4063–4070. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Quevillon, M.J.; Whitmer, J.K. Correction: Quevillon, M.J., et al. Charge Transport and Phase Behavior of Imidazolium-Based Ionic Liquid Crystals from Fully Atomistic Simulations. Materials 2018, 11, 64. Materials 2021, 14, 120. https://doi.org/10.3390/ma14010120

AMA Style

Quevillon MJ, Whitmer JK. Correction: Quevillon, M.J., et al. Charge Transport and Phase Behavior of Imidazolium-Based Ionic Liquid Crystals from Fully Atomistic Simulations. Materials 2018, 11, 64. Materials. 2021; 14(1):120. https://doi.org/10.3390/ma14010120

Chicago/Turabian Style

Quevillon, Michael J., and Jonathan K. Whitmer. 2021. "Correction: Quevillon, M.J., et al. Charge Transport and Phase Behavior of Imidazolium-Based Ionic Liquid Crystals from Fully Atomistic Simulations. Materials 2018, 11, 64" Materials 14, no. 1: 120. https://doi.org/10.3390/ma14010120

APA Style

Quevillon, M. J., & Whitmer, J. K. (2021). Correction: Quevillon, M.J., et al. Charge Transport and Phase Behavior of Imidazolium-Based Ionic Liquid Crystals from Fully Atomistic Simulations. Materials 2018, 11, 64. Materials, 14(1), 120. https://doi.org/10.3390/ma14010120

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop