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

In our publication [...].


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 ].

Text Corrections
The functional form of the OPLS-AA force field is corrected below. Replace 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 interlayer 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 [NO 3 ]and [C 14 MIm] + ) begin to exhibit a more smoothed radial distribution function than the aliphatic tails.  Table 1.

Figure Corrections
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 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
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.
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 trajectoryaveraged 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. 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 12 MIm][NO 3 ], 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.  3 ], 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].

Species Transition Temperature, K This Work Experimental [7] Coarse-Grained [6]
C , 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 xy-plane. The phase transition to the isotropic liquid loses almost all structure, except for pairing of the opposite charges.