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 [CIm][NO] rather than [CMIm][NO].
2. Text Corrections
The functional form of the OPLS-AA force field is corrected below. Replace
with
In addition, with the corrected force field, the nitrate anions are no longer rigid [
3], so remove the following sentence from line 110:
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 , before increasing the temperature by discrete steps of 25 until either the internal temperature is 600 or the material melts into an isotropic liquid. The latter condition was utilized for the 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 , before increasing the temperature by discrete steps of 25 up to 600 .
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,
with
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 and . 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 and (d) 600 . 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 and (d) 600 . 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).
with
Figure 3.
Observed phases in the ionic liquid crystal material, [C
16MIm][NO
3]. Typical configurations are presented from temperatures ranging between (
a) 300
and (
d) 600
. 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, [C
16MIm][NO
3]. Typical configurations are presented from temperatures ranging between (
a) 300
and (
d) 600
. 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).
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.
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.
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.
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.
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.
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.
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 [C
12MIm][NO
3], but observed the isotropic phase above 318
. 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 [C
12MIm][NO
3], but observed the isotropic phase above 318
. As a function of alkyl tail length, the lower transition temperature exhibits a slight non-monotonicity for these species.
Species | Transition | Temperature, K |
---|
This Work | Experimental [5] | Coarse-Grained [6] |
---|
C12 | Crystal–Smectic-A | 425 | - | - |
Smectic-A–Isotropic | 475 | - | - |
C14 | Crystal–Smectic-A | 400 | 336 | - |
Smectic-A–Isotropic | 550 | 356 | - |
C16 | Crystal–Smectic-A | 425 | 339 | 500 |
Smectic-A–Isotropic | 575 | 404 | 560 |
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 [C
12MIm][NO
3], but observed the isotropic phase above 302
. 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 [C
12MIm][NO
3], but observed the isotropic phase above 302
. As a function of alkyl tail length, the trends of transition temperatures are consistent with experimental results [
7].
Species | Transition | Temperature, K |
---|
This Work | Experimental [7] | Coarse-Grained [6] |
---|
C12 | Crystal–Smectic-A | 450 | - | - |
Smectic-A–Isotropic | 500 | - | - |
C14 | Crystal–Smectic-A | 475 | 315 | - |
Smectic-A–Isotropic | 550 | 401 | - |
C16 | Crystal–Smectic-A | 500 | 320 | 505 |
Smectic-A–Isotropic | 575 | 457 | 560 |
Replace
Figure 7.
Mean-squared displacement of the species in [C16MIm][NO3] as a function of time at temperatures of (a) 300 in the crystalline state, (b) 400 just before the transition to smectic, (c) 425 in the smectic-A phase, and (d) 600 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 in the crystalline state, (b) 400 just before the transition to smectic, (c) 425 in the smectic-A phase, and (d) 600 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.)
with
Figure 7.
Mean-squared displacement of the species in [C14MIm][NO3] as a function of time at temperatures of (a) 400 in the crystalline state, (b) 450 just before the transition to smectic, (c) 500 in the smectic-A phase, and (d) 600 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 in the crystalline state, (b) 450 just before the transition to smectic, (c) 500 in the smectic-A phase, and (d) 600 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.)
Replace
Figure 8.
In-plane radial distribution function, of the species in [C16MIm][NO3], trajectory-averaged over the last 200 ns, at temperatures of (a) 300 in the crystalline state, (b) 400 just before the transition to the smectic-A state, (c) 425 in the smectic-A phase, and (d) 600 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 -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, of the species in [C16MIm][NO3], trajectory-averaged over the last 200 ns, at temperatures of (a) 300 in the crystalline state, (b) 400 just before the transition to the smectic-A state, (c) 425 in the smectic-A phase, and (d) 600 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 -plane. The phase transition to the isotropic liquid loses almost all structure, except for pairing of the opposite charges.
with
Figure 8.
In-plane radial distribution function, of the species in [C14MIm][NO3], trajectory-averaged over the last 200 ns, at temperatures of (a) 400 in the crystalline state, (b) 450 just before the transition to the smectic-A state, (c) 500 in the smectic-A phase, and (d) 600 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 -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, of the species in [C14MIm][NO3], trajectory-averaged over the last 200 ns, at temperatures of (a) 400 in the crystalline state, (b) 450 just before the transition to the smectic-A state, (c) 500 in the smectic-A phase, and (d) 600 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 -plane. The phase transition to the isotropic liquid loses almost all structure, except for pairing of the opposite charges.
The authors regret any inconvenience these changes cause to readers, but seek to correct the record in the name of accurate and reproducible science.