1. Introduction
In recent years, increasing concern has been raised about the growing environmental burden caused by the use of fossil fuels and the overall increase in energy consumption, especially in EU countries [
1,
2,
3]. Searching for clean energy, particularly decarbonized sources, is relevant, and according to photovoltaics, PV [
4,
5], Lithium-ion batteries, LIB [
6], wind mills [
7], and fuel cells [
8,
9], these appear as innovating and useful means for overcoming such challenges. The fuel cell is a device that converts the electrochemical energy of a fuel’s oxidation-reduction into electricity, heat, and water. One of the most important components is the membrane, and among them, perfluorosulfonic acid ones (PFSA) as Nafion®, Flemion®, Aquivion®, and 3M ®membranes are commercially available from Chemours, Fumatech, Solvay, and 3M Innovative, respectively [
10]. Many articles have been published on the synthesis, characterization, processing, and electrochemical application into portable, small, and bigger embedded devices (vehicles that are already produced from Toyota [
11], Honda, or Hyundai). It has been found that proton transport properties of PFSA membranes strongly depend on their morphology [
12]. The self-assembly of fluorinated side-groups determines the volume fraction, topology, and size of water channels. Despite numerous publications on the correlation between chemical structure, preparation conditions, and final supramolecular morphology of polymeric membranes, this question remains open due to polymer specificity, such as slow structure formation and polydispersity of model ionomers [
13,
14,
15]. To better understand the self-assembly process in ion-transporting membranes, it would be of interest to propose models (i.e. small molecules) capable of self-assembling under conditions similar to those for fluorinated comb-like copolymers. Such semi-rigid amphiphilic compounds, known as calamitics, are known to form layered, cylindrical, or cubic superstructures as a function of molecular shape, temperature, humidity, etc. [
16,
17,
18]. For example, Berrod et al. [
19] suggested original perfluorosulfonic acid molecules that can reproduce phase-separated morphology of PFSA ionomer membranes used as electrolytes in fuel cells.
Investigation of thermotropic mesogens are also of great interest from the fundamental point of view as the peculiarities of structure and assembly mechanisms of liquid crystals is yet to be fully understood [
20,
21,
22]. Specific details on molecular level arrangements, conditions of transitions between different phases and so on can be resolved through various simulations techniques [
23,
24], and comparison between experimental and theoretical data can be used for the refinements of force fields and models, leading to better understating of physical reasons of mentioned assemblies. Up today there are numerous modeling studies related to thermotropic mesogens [
23,
24,
25,
26,
27] which show the high ability of modern methods to predict and reveal some features of the structuring of mesogens qualitatively as well as quantitatively.
2. Materials and Methods
Synthesis
1,1,2,2-tetrafluoro-2-(1,1,2,2-tetrafluoro-2-iodoetoxy)ethanesulfonyl fluoride was supplied from Appollo while ethylene provided by Air Liquide, France; Perkadox 16S was purchased from Akzo Chemicals and tertbutanol and Li2CO3 from Sigma Aldrich.
Ethylenation of the iodofluorocompounds.
Ethylenation of 1,1,2,2-tetrafluoro-2-(1,1,2,2-tetrafluoro-2-iodoetoxy)ethanesulfonyl fluoride occurred in 50 mL Hastelloyl (HC276) autoclave equipped with inlet and outlet valves, a manometer, a rupture disk, a mechanical stirrer and a controller to check the speeding rate and the temperature. A solution composed of di(tert-butylcycloheylperoxy dicarbonate (Perkadox 16S) (0.247 g, 0.63 mmol), 28 mL of tert-butanol, and 16.82 g (39.4 mmol) of 1,1,2,2-tetrafluoro-2-(1,1,2,2-tetrafluoro-2-iodoetoxy)ethanesulfonyl fluoride was degassed by N2 bubbling for 30 min. During this time, the autoclave was put under vacuum (40 × 10−3 bar) for 30 min to remove the residual traces of oxygen. The above solution was transferred into the autoclave under vacuum via a funnel tightly connected to the introduction valve of the autoclave. Then, the vessel was cooled in a liquid nitrogen bath, and ethylene (2 g, 45.5 mmol) was introduced through the autoclave inlet valve while monitoring its amount by double weighing (i.e., the difference of masses before and after filling the autoclave with ethylene). The autoclave was then allowed to warm to ambient temperature and heated to the target reaction temperature (60 °C) under mechanical stirring and the evolutions of pressure and temperature were recorded. Suddenly, ca. after one hour-reaction at 60 °C, a fast high increase of pressure up to 70 bar occurred immediately followed by a drop of pressure as low as 5 bar, considering a complete reaction. The autoclave was then cooled in an ice bath for 30 minutes and no released of unreacted gas was noted, followed by the opening of the vessel. The total product mixture (yellow solid) was recovered from the autoclave and characterized by 19F NMR. The evidence by the absence of signal at -59 ppm assigned to CF2CF2I end-group indicated the quantitative conversion of 5-iodo-3-oxaperfluoropentane-1-sulfonyl fluoride. The crude product was then dissolved in acetone and precipitated from water and finally dried under vacuum (20 × 10−3 bar, 50 °C) for 8 h yielding 17.55 g of a yellowish powder (yield=94%).
Oxidation of the sulfonyl fluoride in I-CH2CH2CF2CF2OCF2CF2SO2F into SO3H was carried out in slightly basic conditions in presence of K2CO3 to lead to ICH2CH2CF2CF2OCF2CF2SO3Li (yield not determined). This latter was then acidified to chemically change the -SO3Li end-group into sulfonic acid, thus producing ICH2CH2CF2CF2OCF2CF2SO3H as a white powder. The overall yield from 1,1,2,2-tetrafluoro-2-(1,1,2,2-tetrafluoro-2-iodoetoxy)ethanesulfonyl fluoride was ca. 78%.
Methods of Characterization
The morphology and physical properties of the sample were studied through thermogravimetric analysis, differential scanning calorimetry, polarized optical microscopy, temperature-resolved small- (SAXS) and wide-angle X-ray scattering (WAXS), and computer simulation.
Thermogravimetric experiments were conducted using TGA Q500 from TA Instruments from 25 °C to 900 °C at a heating rate of 10 °C under a nitrogen flow 100 ml/min.
Differential scanning calorimetry (DSC) measurements were carried out using a Polyma 214, Netzsch. Samples were placed in pierced aluminum pans and subjected to two heating − cooling cycles from 20 to 260 °C at a rate of 10 K.min−1.
Microscopic images of thin films were captured using a Carl Zeiss AxioScope A1 POL optical microscope in polarized light with 100x objective. The films, prepared from the melt between two cover glasses, were positioned on a Linkam LTS heating stage. A small amount of the sample was heated to 260 °C, above melting temperature, and then crystallized at two different cooling rates: 1 and 10 K.min−1. Images were captured with a 5 MP CMOS camera.
Temperature-resolved SAXS and WAXS measurements were carried out at the ID02 beamline of the European Synchrotron Radiation Facility (ESRF) in Grenoble (France). The measurements were performed in transmission geometry with X-ray photons at an energy of 11.9 keV. The calibration of the norm of the reciprocal vector s ( where θ is the Bragg angle and λ-wavelength) was done using several diffraction orders of silver behenate (for SAXS) and corundum (for WAXS).
All-atom molecular dynamics (AAMD) simulations of I(CH
2)
2(CF
2)
2O(CF
2)
2SO
2F melt under different temperature conditions were carried out as follows. The initial unit cell configuration was established through Monte Carlo annealing among a set of most common symmetry groups, with geometry optimization using the COMPASS force field [
28]. The elementary configuration with the lowest energy was then replicated to form a periodic box measuring 100 Å x 30 Å x 30 Å, containing 288 I(CH2)
2(CF
2)
2O(CF
2)
2SO
2F molecules. The COMPASS force field was utilized for the AAMD simulations, and electrostatic interactions were computed using the PPPM method [
29].
The system was equilibrated in the
isobaric-isothermal ensemble (NPT) for 10 ns with a Berendsen barostat [
30] and a Nosé -Hoover thermostat [
31], at a temperature of T = 25 °C and a pressure of P = 1 atm. Subsequently, the simulation box was heated sequentially by 5 °C every 10 ns up to 280 °C. Statistics were gathered over the last 1 ns of every heating step across three independent heating procedures.
Radial distribution function (RDF) and static structure factor S(k) were utilized for charactering the translational order of molecules within the box. Orientational ordering was monitored by calculating the first eigenvalue of nematic Q-tensor,
[
32], where γ is the angle between molecules' first principal axis and the director. The size of the molecules was estimated based on the calculation of the gyration radius R
g.
3. Results and discussion
1,1,2,2-tetrafluoro-2-(1,1,2,2-tetrafluoro-4-iodobutoxy)ethanesulfonic acid was synthesized in three steps from 1,1,2,2-tetrafluoro-2-(1,1,2,2-tetrafluoro-2-iodoetoxy)ethanesulfonyl fluoride (
Scheme 1): i) by a quantitative ethylenation initiated by
tert-butyl cyclohexyl peroxydicarbonate at 60 °C (that generates only the monoadduct as in previous works [
33]. The vanishing of the signal at -59 ppm (highlighted by
19F NMR spectrum) and the presence of two complex multiplets at 2.3 and 3.0 ppm in the
1H NMR spectrum evidenced the formation of the ethylenated product (with the insertion of one ethylene unit only); ii) the second step enabled to modify the sulfonyl fluoride end group into SO
3Li under mild basic conditions to avoid the dehydroiodination. iii) the final steps dealt with the acidification of that latter into sulfonic acid. These reactions were monitored by
19F,
1H and IR spectroscopies. The O-F frequencies (bending modes) centered at 1470 and 810 cm
-1 (noted in the spectrum of the product bearing SO
3Li end-group) were absent after the hydrolysis reaction while the -S(O)(O)OH stretching vibrations of the sulfonic acid groups at 3377 cm
-1 and 1026 cm
-1 were present, indicating the successful hydrolysis reaction. Additionally, the hydrolysis of the -SO
2F groups was confirmed by the
19F NMR spectrum that does not display any signal at +45 ppm characteristic of the sulfonyl fluoride group [
34]. The absence of this peak indicates that the hydrolysis-reaction yield was quasi quantitative.
The overall yield was 78 % from 5-iodo-3-oxaperfluoropentane-1-sulfonyl fluoride.
Structural Characteristics of I(CH2)2(CF2)2O(CF2)2SO2OH
The first DSC heating thermogram shows a double endothermic peak within the temperature range of 80 to 110 °C, which is attributed to the sample’s dehydration or dehydroiodination. The peak at 150 °C is likely due to a crystal-to-crystal transition, typical for amphiphilic calamitic molecules [
35]. An intense peak at 218 °C corresponds to the sample’s isotropization [
36]. The second heating curve shows only isotropization, indicating irreversibility of the solid-to-solid transition.
Figure 1.
DSC thermograms of I(CH2)2(CF2)2O(CF2)2SO2OH: 1st heating (yellow), cooling (blue) and 2nd heating (red). Exo processes are down.
Figure 1.
DSC thermograms of I(CH2)2(CF2)2O(CF2)2SO2OH: 1st heating (yellow), cooling (blue) and 2nd heating (red). Exo processes are down.
Polarized optical microscopy of thin films of I(CH
2)
2(CF
2)
2O(CF
2)
2SO
2OH at 216 °C reveals crystalline spherulites and rod-like domains of the smectic liquid crystalline (LC) phase (
Figure 2a). It is plausible that the LC phase is thermotropic, which signifies that it is in equilibrium in a certain temperature range. During cooling, the nucleation rate increases, leading to the emergence of a grainy polycrystalline texture at room temperature (
Figure 2b).
Structural analysis was performed using temperature-resolved small- (SAXS) and wide-angle (WAXS) X-ray scattering. At room temperature, powder diffraction indicates a highly ordered crystalline structure with numerous peaks in the WAXS region (
Figure 3). The SAXS region shows peaks typical for a layered structure. The complete table of the peaks is presented in
Table S1. According to SAXS, the longest inter-layer distance is 2.72 nm. The high-intensity peaks at s=4.3 and 4.9 nm
-1 are associated with aluminum foil used as a sample container.
The crystalline structure observed can be assigned to a monoclinic unit cell with lattice parameters: a = 27.5 Å, b = 5.2 Å, c = 5.0 Å, and angles α = 81°, β = 90°, and γ = 90°. Temperature-dependent SAXS and WAXS measurements are presented in
Figure 4(a) and
Figure 4(b), respectively. During heating, distinct changes in both peak positions and intensities are noticeable in the SAXS and WAXS regions. Specifically, shifts observed at 90 °C and 150 °C are attributed to dehydration and a crystal-to-crystal transition, respectively. The isotropization process at 220 °C becomes evident with the fading of sharp crystalline peaks.
Indexation of the high-temperature phase indicates the formation of a triclinic unit cell with the following parameters: a=25.6 Å, b=5.12 Å, c=4.96 Å, α=95°, β=97°, γ=88°.
To gain insights into the arrangement of I(CH2)2(CF2)2O(CF2)2SO2OH at different temperatures, all-atom molecular dynamics simulations of a box sized 100 Å x 30 Å x 30 Å, filled with 288 molecules, were conducted at room temperature, followed by a heating procedure with a temperature increase up to 280 °C. A detailed description of the modeling methodology is provided in the Materials and Methods section.
Initially, after system equilibration at T=25 °C, the molecules are elongated and retain crystal ordering. Approximation of the lattice parameters, obtained through identification of the minimal symmetric parallelepiped [
37] and averaged over the system, yielded the following values:
a = 27.2 Å,
b = 5.0 Å,
c = 5.0 Å, 𝛼 =81.1°, 𝛽 = 89.6°, 𝛾 = 90.4°. These parameters correspond to a cell with monoclinic symmetry and are in good agreement with experimental data (
Figure S1).
After equilibration at room temperature, the box of I(CH
2)
2(CF
2)
2O(CF
2)
2SO
2OH was slowly heated in the NPT (isobaric-isothermal ensemble) ensemble from T=25 °C to T = 280 °C. Within this temperature range, the system undergoes a series of transitions, analyzed using nematic order parameter S [
32], which indicates whether the molecules’ main axis is oriented in one direction (S=1) or their orientations are distributed isotropically (S=0). The presence of a layered structure was monitored through the calculation of the static structure factor S(k/k*).
Figure 5 shows the dependency of the nematic order parameter S on temperature, the insets provide typical snapshots of the system under different states and the normalized static structure factor for the chosen T. The latter was plotted against k/k*, where k* is the location of the first maximum. For the layered structure, positions of consequent peaks k(1):k(2):k(3): k(4) should satisfy the relations 1:2:3:4 [
38], indicating the periodicity of layers. This relation breaks when layered ordering is disrupted. Positions of these peaks are marked with vertical lines.
Overall, the dependency of the nematic order parameter on temperature is qualitatively similar to what is expected for ellipsoidal mesogens [
39,
40], with S close to 1 at low temperatures (crystal) and decreasing monotonically as T increases. Sharp changes in the order parameter S indicate transitions between different states (crystal, smectic, nematic, isotropic), with the most pronounced drop at the nematic-isotropic transition.
In the interval between 25 °C and ~ 170 °C, the molecules maintain a layered arrangement (peaks at k/k* = 1, 2, 3, 4 at T = 150 °C,
Figure 5 upper right,
Figure S2(left)) and retain a high orientational ordering parameter, which decreases from S = 0.9 to S = 0.64 upon heating. At T ~ 170 °C, layers begin to overlap, leading to the destruction of the smectic phase - at T = 180 °C, only the first two peaks of S(k/k*) remain (
Figure 5 bottom right,
Figure S2(middle)). However, the system maintains high orientational order for I(CH
2)
2(CF
2)
2O(CF
2)
2SO
2OH (S ~ 0.64). Thus, at T ~ 170 °C the smectic to nematic transition can be observed. The model temperature for this transition to the nematic phase was found to lie within the interval between 170 and 245 °C, with a monotonic decrease of S upon heating from 0.64 to 0.42. A subsequent increase in T leads to a sharp drop in the nematic order parameter to 0.19 at T = 250 °C, with the nematic-to-isotropic transition temperature occurring approximately at T ~ 245 °C,
Figure S2 (right).
Therefore, within the all-atom molecular dynamics framework, two pronounced transition temperatures were identified: smectic-to-nematic at T ~170 °C and nematic-to-isotropic at T ~ 245 °C, both of which are in good agreement with experimental data (
Figure 3). The fact that the values derived from simulations are slightly higher than those observed experimentally (i.e., 150 and 230 °C, respectively) can be attributed to finite-size effect or to inaccuracies within the chosen force field.