Dry-Jet Wet Spinning of Thermally Stable Lignin-Textile Grade Polyacrylonitrile Fibers Regenerated from Chloride-Based Ionic Liquids Compounds

In this paper, we report on the use of amorphous lignin, a waste by-product of the paper industry, for the production of high performance carbon fibers (CF) as precursor with improved thermal stability and thermo-mechanical properties. The precursor was prepared by blending of lignin with polyacrylonitrile (PAN), which was previously dissolved in an ionic liquid. The fibers thus produced offered very high thermal stability as compared with the fiber consisting of pure PAN. The molecular compatibility, miscibility, and thermal stability of the system were studied by means of shear rheological measurements. The achieved mechanical properties were found to be related to the temperature-dependent relaxation time (consistence parameter) of the spinning dope and the diffusion kinetics of the ionic liquids from the fibers into the coagulation bath. Furthermore, thermogravimetric analysis (TGA), differential scanning calorimetry (DSC), and dynamic mechanical tests (DMA) were utilized to understand in-depth the thermal and the stabilization kinetics of the developed fibers and the impact of lignin on the stabilization process of the fibers. Low molecular weight lignin increased the thermally induced physical shrinkage, suggesting disturbing effects on the semi-crystalline domains of the PAN matrix, and suppressed the chemically induced shrinkage of the fibers. The knowledge gained throughout the present paper allows summarizing a novel avenue to develop lignin-based CF designed with adjusted thermal stability.


Introduction
Carbon fiber (CF) reinforced thermoplastics and thermosets combine high stiffness and high strengths with low density, enabling them as light-weight construction materials for space and aviation sectors, wind energy, and automotive industry.
The high tensile strength and high tensile modulus of CF are owing to the unique preferably along with the fiber axis oriented turbostratic or graphite-like carbonaceous crystal structure. However, the characteristic and unavoidable porosity, as well as the pore dimensions and misorientation of the pores with respect to the fiber main axis, govern, to a wide extent, the mechanical performance of the fibers. In order to achieve certain tensile strengths and tensile moduli of CF, defined structural parameters should be maintained. Some of these parameters are as follows:  crystal thickness Lc;  crystal width parallel La// and perpendicular La┴ to the fiber axis in meridional and equatorial modes;  interlayer spacing of the basal planes d002;  relative orientation degree and Herman's orientation factors of the crystals with respect to the fiber main axis;  porosity;  true pore length, overall tilt angle (misorientation) of the pores with respect to the fiber main axis;  specific surface of the pores;  intersection length and radius of gyration;  orientation degree of the pores alongside the fiber main axis.
The structural characteristics can be adjusted utilizing over 200 manufacturing parameters of the precursor system used and the manufacturing process. In one of our last papers, we reported in detail about the structure-property relation of CF and the impact of the manufacturing process on the resulting properties [1].
The global demand of CF increases yearly at a rate of 10-13% (2010-2017) [2][3][4][5]. Co-polymers of polyacrylonitrile constitute as precursors ca. 90% of the global CFCF market. Ten percent of the CF is made from mesophase pitch and isotropic pitch. One of the critical drawbacks of CF was high costs, making the breakthrough of this material in the low-market segments, that is, the automotive industry, more difficult. The precursor costs of polyacrylonitrile account for slightly more than 50% of the costs of PAN-based CF [5]. CF with reduced costs should also maintain a certain level of mechanical properties, that is, tensile strength >1.7 GPa, tensile modulus >210 GPa, as well as elongation at break >1% [6].
From the economic viewpoint, different research methodologies were recently reported [6]. For instance, utilization of heavy tow textile grade acrylic fibers prepared by pseudo-melt spinning technology and stabilization via advanced oxidative stabilization technique, that is, the microwave assisted plasma was reported [7].
Another important aspect of CF production is the low carbon yield (up to 50%) of the polyacrylonitrile. This means that the production capacity of a PAN-based CF plant is one half of that of PAN plant.
Lignin is a by-product of the pulping process and a low-cost material. The high abundance, onethird of the dry weight of the lignocellulosic biomass, and the aromatic character of this material are of a paramount importance for further economic utilization for high value-added products, such as CF precursors. Caused by the pulping process, Lignin showed consistency problems, making it unsuitable for fiber formation, that is, very low molecular weight values and the high polydispersity lead to brittleness of the fiber structure, and the high ash content leads to defects and flaws in the fine fiber structure.
Utilizing lignin as a precursor material for high value-added products like CF is not a new research topic. Different research routes were pursued to formulate lignin-based precursor fibers and, subsequently, lignin-based CF with high mechanical properties. To date, lignin-based CF show low mechanical performance, that is, low tensile strength < 1 GPa, low tensile modulus <150 GPa, and low elongation at break <1% [6]. Additionally, the X-ray porosity of lignin-based CF is high and exceeds 30-38% depending on the thermal conversion [8]. The low molecular weight decreases the probability of conjunction and joining zones between the carbonaceous crystals in the subsequent lignin based CF. The high ash content increases the possibility of the development of flaw defects in the fine structure of CF, and hence deteriorates the mechanical performance of the fibers. The high content of hydroxyl groups hinders the fusibility of lignin and the concentration and constitution of methoxy groups increase the glass transition region of the lignin [9,10]. This hinders in turn the processability of the material utilizing eco-friendly processes like melt spinning. In addition, the high content of hydroxyl groups and the low molecular weight (as a molecular origin) increase the agglomeration tendency of the lignin fractions even in strong solvents, that is, in the water-based binary solvent system with N, N-dimethylformamide (DMF), dimethyl sulfoxide (DMSO), or dimethylacetamide (DMAc). This leads in the wet-spinning process with polyacrylonitrile to leaching phenomena in the coagulation bath [11].
Additionally, the high hydroxyl group content of lignin increases the hydrogen bonding and hinders the fusibility of the polymers. This is the main reason why lignin cannot be melt-processed without any chemical modifications. The melt-spinning and compatibilization of lignin with different polyolefin polymers and liquid crystal polymers as a supporting matrix require a chemical derivatization, that is, alkylation or acetylation [12][13][14].
Most of the research works regarding lignin-based CF mainly considered the utilization of softwood Kraft lignin types for the following shortly summarized reasons.
Softwood Kraft lignin can be technically extracted from the black liquors with a relatively low (chemically bounded) sulfur content (1-3 wt. %), for instance, utilizing the Lignoboost ® -Process [15], compared with the lignin extracted from the sulphite-process, for example, magnesium or sodiumbased sulphite process.
Organosolv hardwood lignin is very pure (in terms of ash content) and almost sulfur free owing to the organosolv extraction process; however, the high methoxy group content was found to reduce the thermal stability of the lignin compared with the softwood lignin [16,17].
Softwood lignin is highly susceptible to any thermal processing, which initiates radical polymerization and increases the molecular weight owing to the cross-linking. This mechanism shortens the stabilization time required during the oxidative stabilization of softwood lignin-based precursor fibers [18]. Different studies report the compounding of softwood and hard wood lignin to increase the thermal stability of the mixture and enhance the thermal oxidative stabilization [19,20].
To date, lignin-based CF, regardless of the PAN/lignin-based CF, show low mechanical properties, and the achieved properties are inferior to the set targets [6].
The absence of the semi-crystalline structure of lignin prevents the orientation of the molecules alongside the fiber axis during the high shear stress extrusion through fine spinning nozzles.
Recently, Clauss et al. reported the possibility of increasing the molecular weight of lignin utilizing the reactive compounding. It is assumed that, the higher the molecular weight, the better the orientation of lignin molecules alongside the fiber axis would be. Residual chars up to 49% @ 800 °C under nitrogen purge were achieved [13].
A crucial drawback is the extremely low thermal stability of lignin [21,22] and the undefined thermal decomposition owing to the lack of semi-crystalline sequences or domains. This becomes visible in terms of higher porosity of the residual chars at elevated temperatures (above 2000 °C) and can be evaluated utilizing the µ CT-measurements. Comparatively, Steudle et al. reported in their work [8] the increasing X-ray determined porosity of the lignin-based CF at high treatment temperatures. The measured porosity values reached 30-38% at a heat treatment temperature of 2000 °C. This is owing to the fact that lignin in general consists of non-repetitive or random chemical bonds, which was found to enhance the release of volatile products in the early pyrolysis region between 200 and 600 °C and to produce weakly cross-linked clusters, which are furthermore thermally unstable [21]. On the basis of our recent work, it is well known that the origin of the pores in CF is mainly caused by the amorphous regions in the precursor fibers [1,6]. The randomly constituted amorphous lignin oligomers inevitably lead to an undefined pore structure in the subsequent CF.
Recently, different research works aimed to develop continuous precursor fibers based on different trinary systems, for example, lignin/pure cellulose/ionic liquids or lignin/Kraft pulp/ionic liquids, utilizing the dry-jet wet spinning technology [23][24][25][26]. The utilization of cellulose aimed to over-compensate the low molecular weight of lignin. Strong and homogenous precursor fibers with increased thermal stability and higher carbon yield for the systems (lignin/cellulose) could be achieved. Nevertheless, not enough structural data of the fiber fine structure were reported. Moreover, the impact of lignin on the fine structure and the thermo-mechanical behavior of the fibers were not discussed in depth.
Byrne et al. [27] studied the stabilization mechanisms of cellulose/lignin fibers regenerated from the ionic liquid: 1,5-diazabicyclo(4.3.0)non-5-ene-1-ium acetate (DBNH + ) (OAc -) and the impact of lignin on the stabilization using FT-IR techniques. Unfortunately, no structural data regarding the impact of lignin on the thermo-mechanical behavior of the fibers, that is, thermally induced entropic shrinkage, oxidative shrinkage, and the influence of lignin on the fine structure of the fibers, were reported.
In order to gain CF based on lignin/PAN in combination, main attention should be paid to the adjustment of the thermal stability of the starting material, as well as the ordering and orientation of the lignin oligomers alongside the fiber axis. This fundamental research methodology will be the scope of the present research paper. Additionally, the present paper focuses on the impact of ionic liquids, softwood lignin oligomers, and the processing parameters on the thermal, thermomechanical, and mechanical properties of the precursor fibers. Furthermore, the impact of the processing parameters and rheological data will be correlated to the fine structure of the fibers, that is, crystal sizes, interlayer spacing, relative orientation degree, and Herman's orientation factor. The impact of lignin on the stabilization kinetics and the kinetics of the thermo-mechanical behavior of the fibers will also be discussed in depth throughout the present paper.
Textile grade homo-polyacrylonitrile (HPAN) with acrylonitrile content greater than 99.5 wt. % and molecular weight of approximately 200,000 g·mol −1 was supplied by DOLAN GmbH, D-93309 Kelheim, Germany. Softwood BioChoice TM Kraft Lignin (SBKL) was obtained from UPM Biochemicals, FI-00101 Helsinki, Finland. Softwood Kraft lignin was selected, as Softwood Kraft lignin is industrially available as lignin-type and, upon the thermal conversion process, shows lower void formation tendency, as the content of the methoxy groups is lower than the Hardwood organosolv lignin. SBKL was dried at 80 °C under vacuum for 24 h prior to further characterization, purification, and further processing. The molecular characteristics of the utilized lignin are summarized in the supplementary data.
Ash removal from SBKL: For CF precursors [6], lignin should be extensively purified to remove any contaminating alkali metals or free sulfur. This is of importance to avoid any possible flaws generated in the fibers by the ash, as reported by Baker et al. [28], Klett et al. [29], and Ishikawa et al. [30].
The impurities in lignin including acid soluble carbohydrate substitutes and ash have been extracted via the dialysis process [31]. Around 200 g of dried SBKL powder was dispersed under vigorous stirring at 1200 rpm in a two-liter solution of hydrochloric acid with the molarity of 2.5 M at 60 °C for one hour. The pH was controlled between 1 and 3. After washing, the SBKL was precipitated and separated from the acidic solution via a common filtration set-up using a ceramic Brüchner funnel and flask connected to a vacuum pump. Ash-free filter papers (Grade 390, Pore retention 3-5 µ m) were supplied by Sartorius Stedim Biotech GmbH, D-37,079 Göttingen, Germany.
The as-separated SBKL cake was washed under rigorous stirring with four liters of hot (T = 60 °C) deionized water. After successive washing, the cake was washed off with deionized water and dried at 70 °C under vacuum for 24 h.
After drying, the remained ash content of lignin was investigated according to the norm DIN EN ISO 3451-1:2008. The ash content of the neat SBKL was 1.29 wt. % and, after the acidic washing, the ash content was reduced to 0.15 wt. %.

Micro-Compounder
In order to evaluate the miscibility of the trinary system textile grade polyacrylonitrile, Biochoice TM Softwood Lignin and ionic liquids the DSM Xplore ® micro compounder, Netherland, was utilized. The DSM micro compounder was equipped with conveying fully intermeshing twin screws and a chamber volume of 15 cm 3 , and was used to develop different compounds of lignin, textile grade PAN, and ionic liquids.

Shear Rheological Measurements
The shear rheological measurements of the developed trinary system dopes (PAN/Softwood BioChoice TM /Ionic Liquids) and the binary system dopes (PAN/Ionic Liquids) were investigated utilizing the ARES G2 strain-controlled oscillatory rheometer (TA Instruments, 159 Lukens Drive, New Castle, UK) to describe their viscoelastic pseudo-flow behavior and their thermal stability. The strain resolution was 0.04 µ rad and the minimum angular displacement oscillation was 1 µ rad. The samples were heated in a forced convection furnace using nitrogen as a heating medium. Cone-plate geometry made from stainless steel with the diameter of 25 mm and cone angle of 0.0997 rad was used. The gap offset was 0.0559 mm, giving a sample volume of 0.4091 cm 3 .

Oscillatory Strain Sweep Test
To identify the spinning dopes for being linear or non-linear viscoelastic, the strain dependency of both the storage moduli (G') and loss moduli (G") was evaluated at a constant frequency and two different temperatures (85 °C and 105 °C). Additionally, the yield stresses of the spinning dopes were estimated for the studied dopes under the same conditions.

Oscillatory Frequency Sweep Test
To ensure the comparability of the data of different spinning dopes, the deformation amplitude of 6% in the linear viscoelastic range was selected for all dopes from the strain sweep tests. At this strain amplitude, frequency sweep measurements at five different temperatures with temperature increments of 10 K in the temperature range 65 °C to 105 °C were performed within a modest frequency range from 100 Hz to 0.1 Hz and a resolution of nine points per decade.
The master curves of the complex viscosities as a function to the angular frequency were displayed for five different reference temperatures.
Assuming the applicability [32] of the Cox-Merz empirical relationship [33] to the binary or trinary polymer system of the spinning dopes at low frequencies, the plot of the steady state shear viscosity versus the shear rate ̇ was assumed to correspond to the magnitude of the complex viscosity | * | versus the angular frequency : From the plot of the steady state shear viscosity versus shear rate , the Carreau-Yasuda fit model was used to describe the pseudo-flow behavior and to calculate the zero shear rate viscosity as well as the consistency (relaxation time) of the spinning dope at the related five reference temperatures.
where 0 is the zero shear rate viscosity, ∞ is the infinite shear rate viscosity, k is the consistency parameter (characteristic relaxation time), ̇ is the shear rate, n is the power low index, and a is a parameter describing the flow transition between the Newtonian plateau and the power law region.
The TTS was evaluated from the logarithmic value of the horizontal shift factors αT of the storage modulus at 85 °C versus the term 1/T. If there was no break in the TTS plot, the system was considered miscible and the activation energy of the flow was calculated from the slope of the linear regression of the TTS plot using the Arrhenius model.
where is the flow activation energy, R is the ideal gas constant (8.31448 J·K −1 ·mol −1 ), and T0 is the reference temperature.

Dry-Jet Wet Spinning
Selected spinning dopes were compounded on a self-constructed vertical kneader under nitrogen purge. Prior to the spinning trails, the dopes were filtered and degassed overnight. Dry-Jet Wet spinning trails were performed with a modular laboratory scale Dry-Jet Wet spinning system at Tampere University, Finland. A titanium spinneret nozzle made from with a 300 µ m cylindrical nozzle with the aspect ratio L/D of 2 and an internal opening cone angle of 60° was used. The degassed and filtered spinning dope was permanently forwarded from the sealed tank to the Zeith ® gear pump by nitrogen gas with the pressure of 4 bar at 85 °C. The gear pump achieved a minimum rotational speed of 0.5 rpm. This in turn allows (depending on the spinning dope viscosity) a minimum throughput of 0.08-0.09 g·min −1 . The Dry-Jet Wet Spinning machine allowed to increase speed of the first take-up godet between 2.5 and 50 m·min −1 . The maximum possible overall in-line stretching of the attenuated fibers was 10 times.
The spinning set-up (spinning dope tank, spinning gear, spinning tubes, and spinning nozzle) was permanently heated by a silicon heating tape system, type KM-HT-203-SO, SAF Wärmetechnik GmbH, D-69509 Mörlenbach). The temperature of the spinning set-up and spinning nozzle was permanently controlled and adjusted to 85 °C for all spinning trails. The coagulation (phaseseparation process) of the attenuated filaments was performed in a deionized water bath with an initial specific conductivity of 0.5 µ S·cm −1 at 25 °C. All spinning trials were performed at a constant air gap of 5 cm between the spinning nozzle and the surface of the water bath and a constant length of the coagulation line in the water bath of 80 cm.
The temperature of the coagulation bath was controlled via an ethylene glycol cooling unit and high surface area stainless heat exchanging unit.
The take-up zone consisted of three glass godets with corrugated surfaces. Between the first and second godets, a heated water bath (60 °C) with deionized water was used to stretch the fibers in-line. All fibers in this work were collected on the second godets under permanent washing with deionized water (23 °C).

Single Fiber Tensile Test
Single fiber tensile test was conducted to evaluate the mechanical properties of the fibers, that is, tenacity, Young's modulus, the elongation at break, and work to rupture, using a FAVIGRAPH from Textechno H. Stein GmbH (Mönchengladbach, Germany). The Young's modulus was determined by the slope of the linear regression of the tenacity-strain curves at terminal values of the strain (0.05-0.5%). The measuring force cell was 100 N, the clamping length was 30 mm, and the constant strain rate was selected to be 1.39%·s −1 (25 mm·min −1 ) at a relative humidity of 65% and ambient temperature (23 °C). For each developed fiber type, 20 sample measurements were performed.
The FAVIGRAPH was also equipped with a vibration cell to determine the linear density (fineness) of each single fiber prior to the tensile test using the vibration method according to ASTM D 1577.

Wide Angle X-Ray Scattering (WAXS)
The supramolecular structure of the developed as-spun fibers and the off-line hot-stretched fibers was studied by the wide angle X-ray scattering (WAXS) using the X-ray diffractometer D8 DISCOVER (Bruker AXS GmbH, Karlsruhe, Germany) with pine-hole collimation system with 0.5 mm diameter and a Cu Kα monochromatized radiation source ( = 0.15418 nm). Large MICROGAPTM Detector (VÅ NTEC-500) with a 10,000 mm 2 gap-free active area was used. The fibers were mounted on a self-constructed Euler Cradle system in order to achieve very wide diffraction angles and to avoid any shadow effects of the clamping system during the measurements.
The crystal size of the develop fibers was calculated, depending on the orientation status of the fibers, from the meridional and equatorial reflection at 2θ~17° using the Scherrer equation: where ℎ is the crystal size related to the (hkl) lattice plane and K is a dimensionless shape factor. For polymeric materials, it has a typical value of (0.87-1.00). In this study, K was 0.9; is the X-ray wavelength and is equal to 0.15418 nm; is the corrected line broadening at half the maximum intensity (FWHM).
can be calculated from the following equation: where is the apparent FWHM of the (hkl) reflection. The instrumental line broadening was 4.538·10 −3 rad for the reflection at 2θ ~ 17° and 6.685·10 −3 rad for the reflection at 2θ~30°. is the Bragg angle and was calculated from the fitting related to (hkl) reflection with the pseudo-Voigt cross function using OriginLab2019.
The interlayer spacing of the reflection planes can be calculated from the Bragg's law as follows: where ℎ is the interlayer spacing of the (hkl) lattice planes, is the wavelength of the incident wave, and is the scattering angle. The average number of planes N in the crystalline domain can be calculated as follows: From the corrected azimuthal intensity distribution of the equatorial (200) reflection, the relative orientation degree (OD) and the Herman's orientation factor fH parallel to the fiber axis were evaluated for all fibers. The azimuthal intensity distribution of the equatorial reflection at 2θ~17° was corrected by the underground from the meridional scan.
The relative orientation degree can be calculated from the corrected plot of the azimuthal intensity distribution as follows: The term (200) is the full width at half intensity maximum of the corrected azimuthal intensity of the reflection at 2θ~17° and can be calculated from a suitable fit function, that is, Pearson VII. The fibers were fully oriented if OD was near to 1 and non-oriented if OD was near to 0. The Herman's orientation factor was calculated from the following equations: The crystals in the fiber structure were perfectly oriented alongside the fiber axis if was near to 1 and random for = 0. If was negative or near to −0.5, this means the crystals were perpendicular or perfectly perpendicular, respectively, to the symmetry axis of the fibers.

Thermogravimetric Analysis (TGA)
The TGA measurements were performed using the TGA Q5000 from TA Instruments (Water Corporation, New Castle, USA) under nitrogen and air. The initial sample weight was ~6 mg and the gas flow was adjusted to be 50 mL·min −1 . The measurements were performed under the heating rates of 5, 10, 20, and 50 K·min −1 .

Shrinkage Measurements
Shrinkage measurements were performed on bundles of precursor fibers consisting of 20 single filaments using DMA Q800 from TA Instruments (Water Corporation, New Castle, USA) at an initial gauge length of 16.20 mm. The force resolution of the device was 0.01 mN, the displacement resolution was 1 nm, and the sensitivity of the loss factor tan δ was 0.0001. The shrinkage tests at different pre loads and different heating rates (1, 2, and 5 K·min −1 ) were performed on a bundle of 20 filaments in the temperature range between 50 and 325 °C to understand in-depth the impact of lignin oligomers on the kinetics of the heat-induced physical entropic shrinkage and on the kinetics of the chemically induced shrinkage during the oxidative stabilization at elevated temperatures under air (oxidative environment) under different loads.

Differential Scanning Calorimetry (DSC)
DSC measurements under nitrogen and under air were conducted for the regenerated neat PAN and PAN/lignin fibers using Discovery DSC 2500 from TA Instruments (TA Instruments (Water Corporation, New Castle, USA) with Fusion Cell TM technology for flat baseline. The stabilization procedures were performed under air to mimic the stabilization process and under nitrogen to reveal the stabilizing role of lignin on the stabilization process of the developed fibers. The weight of the samples was between ~3 and 3.5 mg. All samples were pre-dried under vacuum for 24 at 80 °C and stored in a sealed desiccator over silica gel pellets prior to any measurement. Heating rates of 1, 2, 5, 7, and 10 K·min −1 were selected for the DSC measurements.

Scanning Electron Microscopy (SEM) and Energy-Dispersive X-Ray Spectroscopy (EDX)
For morphology and structure evaluation of the fibers, SEM analyses were carried out on the highly stretched fibers using a high resolution CrossBeam ® NEON ® 40 EsB system with FEG (Field Emission Gun) from Carl Zeiss Microscopy GmbH, Germany and an SEM ULTRA plus with FEG from Carl Zeiss Microscopy GmbH, Germany, and equipped with energy-dispersive X-ray spectroscopy (EDX) XFlash ® FlatQUAD Detector from Bruker Nano GmbH, Germany.
SEM of the fiber surface and cross sections as well as EDX scans of the cross section of the highly stretched fibers embedded in epoxy were performed to reveal topographic information and the existence of chemical elements in the fiber structure. It was important to maintain that the detection threshold of the elements was given by 0.1 wt. %. Any traces under this limit were not able to be detected using this method. SEM images were taken at 3 keV and EDX measurements at 6 keV primary beam energy. All samples were coated with carbon prior to analysis with SEM and EDX to prevent charging artefacts.

Microscopy
For the evaluation of fiber diameter and cross section of the fibers, the fibers were embedded in an epoxy resin and polished to 0.1 µ m. The polished surfaces were then investigated using (VHX 100, Keyence Corporation, Osaka, Japan) polarized reflected-light microscope with magnification up to 2000×. Table 1 shows the composition of the three investigated dopes, which were prepared utilizing the micro-compounder DSM MC-15.  Figure 1 shows the procedure of the preparation of the dopes ND and L1 and the impact of lignin on the reduction of the axial force during the compounding.

Micro-Compounding
All dopes show relaxation behavior after the increase of the compounding speed from 100 to 250 rpm and holding the speed at 250 rpm. For lignin samples, the relaxation rate was higher than that without lignin, implying the softening effect of the lignin oligomers during compounding. This could be induced by the low molecular weight of lignin. Figure 2a shows the linear viscoelastic (LVE) range of the three spinning dopes at 85 °C. The system ND (PAN/IL) showed the lowest critical deformation strain of 16.6%. Spinning dopes with lignin show higher critical deformation strains, which are enhanced with the increasing lignin content, and the values of the critical deformation strain are 21.4% and 27.7% for the dope L1 and L2, respectively. These results suggest that the presence of lignin promotes the resistance of these dopes to deformation, regardless of the low molecular weight of lignin.

Rheological Behavior
In the LVE region, the loss modulus was greater than the storage modulus for all dopes (PAN/IL and PAN/lignin/IL) in the entire deformation range, indicating a dominant viscoelastic liquid-like behavior of all spinning dopes. It was also clear that the lignin-free spinning dope (ND) showed a higher storage modulus in the LVE region than the system L1. With any further increase in the lignin content, a decrement of the storage modulus value can be found in the LVE region. The same behavior was observed for the loss modulus. This implies that lignin acts as a softening agent in the dopes owing to the low molecular weight of lignin.
For further oscillatory investigations, a deformation strain of 6% was adjusted, thus this value was high enough concerning the rotational moment resolution of the utilized rheometer and fulfils the criteria of the LVE region for all dopes. The storage modulus G' and complex viscosity * versus the angular frequency as a log-log plot obtained from the frequency sweep tests are illustrated in Supplementary Data Part II ( Figure  S1). All spinning dopes show strong dependence of the complex viscosity on the temperature, and a general shear thinning behavior at increasing angular frequency and increasing temperatures was detected.
In the terminal region (at low angular frequencies), linear polymer melts typically show the response of a viscoelastic liquid with the scaling property of approximately G′ ∝ ω2 and G″ ∝ ω, indicating a liquid-like behavior.
In the case of the pseudo-molten spinning dopes, at, for example, 85 °C, however, the PAN/IL (ND) dope showed a slope of G' versus the angular frequency of 1.04 and a slope of G" versus of 0.76 at small terminal values of . The dope L1 (PAN/lignin is 2:1) showed corresponding slopes of 1.36 (G' versus ) and 0.95 (G" versus ), while the dope L2 (PAN/lignin was 1:1) showed slopes of 1.03 and 0.92, respectively.
The storage and loss moduli of the dopes show, in general, frequency dependency, and are thus considered to be viscoelastic, although the scaling properties are different from polymer melt, where G′ ∝ ω2 and G″ ∝ ω.
Adding lignin in small concentrations, that is, dope L1, leads to an increase of both the slope of G' versus ω and the slope of G'' versus ω in comparison with the case of neat PAN dopes (ND). For dope L1, the increase of the slope of G' versus ω was higher than the increase of G'' versus ω, implying that lignin oligomers promote more elasticity to the dopes. This suggests, furthermore, some molecular interaction between the PAN and lignin oligomers. Further increase of the share of lignin in the blend, that is, the case of dope L2, however, reduces the slope to values comparable of that of PAN. This suggests the softening effect of lignin at high lignin shares.
From the frequency sweep tests, the master curves of the complex viscosity versus the angular frequency with a reference temperature of 85 °C, for all dope systems, can be seen from Figure 2b. By assuming the applicability of the Cox-Merz rule (|η* (ω)|→η()) for the developed dopes, the steady shear viscosity as a function of the shear rate ̇ can be calculated. All dopes showed a temperaturedependent shear-thinning behavior.
Low molecular lignin reduced the zero-rate viscosity, that is, in the case of dope L1. However, when increasing the lignin share, that is, dope L2, a slight increase of the zero-shear viscosity over that of dope L1 was observed, obviously owing to the increased overall polymer content in the dope. From the transformed curves of the steady shear viscosity versus the shear rate and applying the Carreau-Yasuda fit function, the specific temperature-dependent and composition-dependent characteristics of all the dopes, that is, the zero-rate viscosity and relaxation time, were calculated. Figure 2c shows the dependence of the zero-shear viscosity on temperature for all dopes. In general, the dopes with lignin shares of 33.5 wt. % in the PAN matrix showed lowered zero-shear viscosity over the studied temperature range than that of the neat PAN dopes. Higher lignin shares, that is, dope L2, slightly increased the zero-shear viscosity compared with that of dope L1, owing to the higher overall polymer content in the dopes (please refer to Table 1) Figure 2d illustrates the influence of the temperature and the composition of the dopes on the calculated relaxation time (the consistency parameter). A higher temperature reduced the relaxation time of all dopes and the neat PAN dopes showed the lowest relaxation times. Introducing lignin to the dopes, that is, dope L1, increases the longest relaxation time of the dopes. This effect became strongly visible at lower softening temperatures <85 °C. This implies that the low molecular weight lignin oligomers considerably decelerate the relaxation process of the PAN chains. This effect is very important to be taken into consideration, as the coagulation process will be performed at much lower temperatures.
The compatibility of the dopes was determined using the Han plot, which is temperature invariant and showed a linear correlation in the plot log G' versus log G". The system can be considered as compatible if nearly the same slope of the pure matrix can be observed at different lignin shares in the blend; otherwise, the system was considered to be phase-separated [41]. Figure 3a shows the Han plot of the three dopes at different temperatures (over the frequency range of 0.1 to 100 Hz). The neat PAN dope (ND) showed a slope of 1.61. Introducing lignin to the dope systems at different lignin shares, that is, dope L1 and L2, slightly increased the slope to 1.72 (+ 6.2%) in both cases. Although lignin consists of low molecular oligomers and showed no repeated monomer sequences, the spinning dopes with different lignin content showed temperatureindependent and almost composition-independent slope values of 1.61-1.72, indicating the compatible composition and single-phase consistency between the PAN and low molecular weight purified softwood Kraft lignin BioChoice TM .
Additionally, the miscibility/immiscibility of the PAN, lignin, and ionic liquids was also investigated using the time-temperature-superposition plot. The G' ( ) curves for thermorheologically simple pseudo-melts were superimposed by horizontal shifts (αT) along the angular frequency axis. The vertical shifts (bT) were set to unity. Figure 3b showed the TTS-Arrhenius plot of the three investigated dopes at a reference temperature of 84.9 °C. Regardless of how strong or weak the intra-molecular interactions between lignin oligomers and PAN chains were, the TTS principle was valid for the three spinning dopes. Certainly, lignin oligomers caused no break to the TTS plot at shares up to 50 wt. %, hence the compositions can be miscible.    (Figure 2b). In general, if no molecular changes in the dopes under isothermal conditions take place, the storage modulus should remain constant, and thus time-independent. The system can be, hence, termed as thermally stable. Any decrease in the storage modulus over time will be caused by the molecular degradation affecting the relaxation process, whereas the increase of the storage modulus over the time signifies many molecular changes, for example, cross-linking reaction or gelation. In the case of the filled system, this would mean development in the structural rigidity or elasticity. Figure 3c illustrates the time sweep test of the dopes at 85 °C and a deformation of 6%. After 30,000 s, the neat PAN dopes show an increase of the storage modulus of 32.6% relative to the initial value. Lignin-containing dopes L1 and L2 showed comparatively similar behavior and an increase of the storage modulus by 27.9% and 33.3%, respectively. Figure 3d shows the development of the loss factor in the time sweep tests at 85 °C over the predefined period. Obviously, the neat PAN and lignin-containing dopes show a similar tendency. The loss factors decreased by 11.5% for the neat PAN dope system and by 9% for lignin-based dopes over time, where the increasing rate of the storage modulus was higher than the increasing rate of the loss modulus. This means that the material became, under the applied isothermal conditions over the time, more elastic than viscous. The complex viscosity of the neat PAN system showed an increase of 15%. The complex viscosity of the lignin-containing dopes L1 and L2 increased by 14.6% and 18.7%, respectively, as illustrated in Figure 3d.

Spinning by Dry-Jet Wet Method
The spinning trails were conducted for the dope of neat PAN (ND) and the lignin-containing spinning dope (L1). The spinnability was investigated and the spinnability charts were summarized in Supplementary Data Part III. Depending on the spinnability and the stability in the stretching zone, further investigation matrix was constructed in order to pursue the study of (a) the take up speed at the first godets, (b) the temperature of the coagulation bath temperature and (c) the in-line and post stretching in solid-state on the morphology and fine fiber structure as well as on the mechanical properties of the fibers. Table 2 summarizes the investigated parameters for both dopes PAN/IL (ND) and PAN/lignin/IL (L1). Table 2. Conducted spinning trails to evaluate the impact of different processing parameters, that is, the temperature of the coagulation bath (TC_bath), the take-up speed of the first and the second godets, as well as the off-line drawing on the structure and mechanical properties of the regenerated PAN and PAN/lignin fibers. All spinning trails were conducted at a constant rotational speed of the spinning pump @ 0.5 rpm. *) Post-drawn (2.5 times) in solid-state at 170 °C under air. DR, draw ratio; HPAN, homo-polyacrylonitrile.    Figure 5a showed the influence of the take-up speed of the first godet (V1) on the diameter of the as-spun fibers and on the draw down ratio. Obviously, thicker filaments were obtained from lignin-containing dope in comparison with the neat PAN dope under similar processing conditions (spinning temperature 85 °C; TC_bath = 6 °C).

Evaluation of the Drawing Dependent Fineness and Diameter
In the Dry-Jet Wet spinning process, the flow of the extruded mass in the air gap can be considered as a non-isothermal viscoelastic extensional flow assuming a simplified cylindrical symmetry along the spinning line. Ionic liquids have in general low-volatility [42]. Taking the masstransfer or continuity equation, the momentum equation, the energy equation, and the constitutive equation into consideration, the diameter of the fibers at the first godet is related to the following parameters [43]: where D is the diameter of the fiber in µ m collected on the first godet, Q is the mass throughput in g·min −1 , ρ is the density of the dope in g·cm −3 , c is the polymer fraction in the spinning dope in wt. %, At a constant take-up speed set-up and the spinning nozzle (D0 = 300 µ m), the ratio of the final fiber diameter of lignin-containing fibers D2 to the diameter of the neat PAN fibers D1 was related to the following parameters: The mass throughput of the lignin-containing dopes and the neat PAN was 0.09 g·min −1 and 0.11 g·min −1 , respectively; the polymer content ratio c2/c1 was equal to 18.37/13.04 = 1.41. Taking these parameters into account, the density ratio of the spinning dopes 1 / 2 must be greater than 1 in order to analyze the experimental data of the fiber diameter. The dynamic simulation of diameter formation in the Dry-Jet wet spinning of lignin-containing dopes will be conducted in another paper and was not within the scope of the present paper. A discrepancy in the draw down ratio between the neat PAN fibers and PAN/lignin fibers was detected and the discrepancy was increased by increasing the take-up speed. This effect became intense by in-line stretching and solid state post drawing the filaments, as shown in the Figure 5a.
At a low in-line draw ratio (DR = 2), the discrepancy was ca. 30%. At the highest solid-state drawing (DR = 10), the difference was roughly 51%. This means that the draw ratio-dependent decrement of the diameter (lignin-containing fiber) is much lower than that of neat PAN fibers. This effect can also be seen from the values of the linear density.
The fineness of the developed fibers can be seen in Figure 5b. At a total draw ratio of 10, neat PAN fibers showed linear density of 2.46 ± 0.14 dtex and diameter values of 16.5 ± 0.16 µ m. Neat PAN fibers were thinner than the lignin-containing fibers 5.46 ± 0.37 dtex with diameter values of 24.9 ± 0.27 µ m. From the linear density and the fiber diameter, the theoretically determined density ρ of the fibers can be calculated from the following equation: where T tex is the fineness of the fibers in tex, ρ is the density of the fibers in g·cm −3 , and D is the diameter of the fiber in µ m. The theoretically determined density of highly drawn (DR = 10) neat PAN fibers and highly drawn (DR = 10) lignin-containing fibers was 1.15 g·cm −3 and 1.12 g·cm −3 , respectively.

Mechanical Properties
The PAN and PAN/lignin fibers show typical tenacity-elongation curves, as shown in Figure 6a. Lignin-containing fibers show lower tenacity values under the same manufacturing conditions. This general tenacity effect is owing to (a) the higher fineness values (lower diameters) of the neat PAN fibers and (b) the disturbed fine structure of the lignin-containing fibers by the low molecular weight lignin oligomers. Nonetheless, lignin-containing fibers in general show higher elongation at break values. The tenacity of the neat PAN fibers and lignin-containing fibers increases with higher draw ratio, as seen in Figure 6b,c. As it was possible to produce roughly two times finer neat PAN fibers than the lignin-containing fibers under the same conditions (owing to the higher polymer content and lower density of the lignin-containing dopes), better mechanical properties of the neat PAN fibers were reached. Table 3 illustrates in detail the impact of two important processing parameters, that is, the draw ratio and the temperature of the coagulation bath, on the fineness and mechanical performance of the neat PAN and lignin-containing fibers. Obviously, in both cases of neat PAN fibers and lignin-containing fibers coagulated at 6 °C, the tenacity increases with increased drawing. Generally, lignin-containing fibers show lower tenacities and tensile moduli. The highest tenacity of the neat PAN fibers reached an average value of 99.18 ± 2.99 cN·tex −1 at a total draw ratio of 10. On the basis of the diameter calculations (D = 16.5 ± 0.16 µ m), the tensile strength and the tensile modulus   10. This implies that low molecular lignin oligomers do not generally affect the elongation at break of the fibers. Additional work regarding the critical fiber radius and its dependency on the fineness will be conducted in the future to understand in-depth the impact of lignin on the brittleness of the fiber structure. Furthermore, Table 3 showed the impact of the coagulation bath temperature on the mechanical performance of the developed neat PAN and lignin-containing fibers. Clearly, the higher the coagulation temperature, the thicker the developed fibers offered lower the mechanical properties of the fibers. This tendency was observed for neat PAN fibers as well as for lignin-containing fibers. For instance, in the case of neat PAN fibers drawn four times in-line, the tenacity was reduced from 41.38 ± 1.88 cN·tex −1 to 24.57 ± 0.96 cN·tex −1 (reduction by ca. 68.4%). For the lignin-containing fibers, the reduction was roughly 38.1% from 26.13 ± 0.68 cN·tex −1 to 18.92 ± 0.55 cN·tex −1 . This important effect can be explained by the influence of the diffusion rate of the ionic liquids from the high concentration region (the attenuated fiber being drawn in the coagulation zone) to the low-concentration region (the deionized water in the coagulation bath) being higher at higher temperatures. This shortens in turn the available period to draw and thus orient the polymer chains. This was also obvious in terms of higher fineness values and lower orientation parameters (please refer to the crystallographic fine structure section). The orientation is furthermore only possible if the ionic liquids can act as long as possible as a "lubricant-agent" in the fiber structure while orienting the fibers via the extensional forces being applied in the coagulation zone. Another important parameter, which influences the orientation, is the temperature-dependent relaxation time of the spinning dope. At higher temperatures, the dopes, in general, show shorter relaxation times than at lower temperatures. This was of rather importance in the entry zone (after the air-gap) of the attenuated fibers in the coagulation bath. This effect together with the influence of the thermodynamic diffusion rate can be pursued in detail by the orientation values of the fibers, that is, the relative orientation degree and Herman's orientation factor, as described in the following section. A similar effect of the relaxation time and temperature-dependent diffusion rate of the ionic liquids on the mechanical properties was found in the study of the Dry-Jet Wet spinning of cellulose fibers from ionic liquid dopes [48].

Crystallographic Fine Structure
The structure-process-property relations of the developed fibers were influenced to a wide extent by (a) the material composition, that is, low molecular lignin oligomers, on the deformation behavior of the fibers; and (b) the processing parameters, that is, the coagulation bath temperature, and in-line and solid state stretching. It is well known that the mechanical performance and deformation behavior of the fibers are influenced to a wide extent by the fine structure of the fibers, that is, the crystal structure and the orientation parameters of the crystalline domains. For each type of fiber, equatorial and meridional scans were performed at different drawing ratios and at different coagulation bath temperatures. The azimuthal scans of the intensity of the PAN reflection at 2θ~17° were recorded.  Figure 7a shows the equatorial (e) and meridional (m) wide-angle X-ray diffraction patterns of the as-spun fibers, in-line drawn fibers, and solid-state post drawn fibers. As-spun neat PAN and lignin-containing fibers (DR = 0) and in-line drawn fibers (DR to 4) show diffuse reflection at 2θ~17° as well as a broad diffuse scattering maximum 2θ values between 22° and 30°. The peak intensity of the reflection at 2θ~17° in the equatorial mode increased and the full width at half maximum decreased with the increased stretching, as shown in Figure 7c,d. Only at high stretching values (DR = 10), neat PAN and lignin-containing fibers showed the highest intensity of the equatorial reflection at 2θ~17° and an equatorial peak at 2θ values ~30° became clear in both cases (Figure 7a,c,d).  Table 4 summarizes the crystallographic fine structural data of the neat PAN fibers and lignincontaining fibers. Additional information about the impact of some processing parameters was also listed. As lignin consists of low molecular oligomers with no defined semi crystalline structure, the semi-crystalline structure of the lignin-containing fibers can only be originally given by the PAN chains. At low temperatures (6 °C) of the coagulation bath, neat PAN fibers show increased equatorial cryptal size L2θ~17°(e) with increasing in-line drawing. In contrast, the meridional crystal size L2θ~17°(m) decreased with increasing drawing. The crystal size ratio L2θ~17°(e)/L2θ~17°(m) increased around 39% at inline draw ratio of 4. The interlayer spacing d2θ~17° = 5.3 Å of the PAN chains was found to be similar to that of the original powder 5.26 Å at all drawing ratios. At a total drawing ratio of 10, the equatorial crystal size L2θ~17°(e) was 12.83 nm. The equatorial crystal size of the reflection at 2θ~30° L2θ~30°(e) was 7.70 nm and the associated interlayer spacing d2θ~30° was 3.0 Å . The aspect ratio of the unit cell b/a (Figure 7b) was equal to 1.732 and the crystal unit is thus hexagonal [49].
Lignin deteriorated to a wide extent the size of the crystal domains, while the characteristics of the unit cell by PAN were still preserved (Figure 7b). In other words, the interlayer spacing d2θ~17° was 5.3 Å and the d2θ~30° was 3.0 Å , hence the aspect ratio of the unit cell b/a was 1.732. The crystal size of the highly oriented lignin-containing fibers was less than that of PAN at the same draw ratio. At 10 times drawing, the crystal size L2θ~17°(e) was 7.16 nm and L2θ~30°(e) was 4.92 nm. This gave a lower number of reflection plains N2θ~17° = 15 and N2θ~30° = 17 than that of neat PAN fibers (N2θ~17° = 26 and N2θ~30° = 26). This implies that lignin did not affect the unit cell characteristics, but disturbed the number of PAN chains contributed in the crystal domain, as seen illustratively in Figure 7b. By increasing the temperature of the coagulation bath, the crystal size was reduced in both cases of neat PAN and lignin-containing fibers as a result of the fast diffusion of ionic liquids from the extruded fibers in the coagulation bath, as seen in Table 4 in both PAN and PAN/lignin-series. The aspect ratio of the L2θ~17°(e)/L2θ~17°(m) at a defined drawing was reduced by increasing the coagulation temperature. For instance, this ratio was reduced from 1.53 for neat PAN fibers drawn at DR = 4 to 1.28 for neat PAN fibers coagulated at 20 °C. This suggests that, at higher diffusion rates of the ionic liquids at the equivalent extensional deformation in the coagulation bath (V1 was 10 m·min −1 for both cases), the ionic liquids left the fiber at a very early stage, and the chains became unsoftened and thus less able to be oriented. Lowered crystal size was a result of the less orientation, as seen in the Table 4. The similar effect was found in the case of lignin-containing materials. In summary, it was obvious that lignin does not affect the increasing rate of the ratio L2θ~17°(e)/L2θ~17°(m), but deteriorates to a wide extent the crystal size L2θ~17°. From the azimuthal distribution of the reflection at 2θ~17°, the relative orientation degree (OD) and the Herman's orientation factor of all developed fibers were calculated.
As shown in Table 4, OD increased by increasing the drawing ratio. For the fiber series coagulated at 6 °C, the highest relative OD was 94.7% for neat PAN fibers and 88.9 for lignincontaining fibers. This implies that lignin reduced and blocked the ability of the PAN chains to be oriented during stretching. The same effect was also visible in terms of Herman's orientation effect, which was 0.73 and 0.69 for neat PAN and lignin-containing fibers at DR of 10, respectively.
Increasing the coagulation temperature to 20 °C reduced the orientation degree for both series of the developed fibers compared with those fibers coagulated at 6 °C and at the same in-line draw ratio of 2 and 4, as seen in Table 4. These findings support the explanations of the generally lowered mechanical properties at a higher coagulation temperature discussed in the section of the mechanical properties.
Additional research work was also required to investigate the influence of lignin on the periodical structure length (ordered domain length) of the crystal units parallel to the fibers axis as well as the true void length and tilt angle of the voids with respect to the fiber axis. This can be done using the small angle X-ray scattering and will be discussed in another extended research work.  Obviously, neat PAN fibers offered lower char residue values than PAN/lignin fibers at a similar heating rate. This implies that PAN/lignin compounds regenerated from ionic liquids show lower thermal degradation rates, and thus possessed higher thermal stability. This could be initiated by cross-linking reactions caused by the ionic liquids. The related DTG curves showed significantly lowered degradation peaks of lignin-containing fibers. The maximum degradation rate was less than 0.3%·K −1 within the temperature range between 250 and 500 °C in the case of lignin-containing fibers and around 1.5%·K −1 for neat PAN fibers in the same temperature range.
The incorporation of lignin into the PAN system and regeneration from ionic liquids dopes decreased to a wide extent the degradation rate, suggesting that the used softwood lignin participates in and influences to a large extent the degradation and thermally induced cross-linking mechanisms of PAN. Considering that the degradation takes place at DTG values above 0.05%·K −1 , the degradation of lignin-containing fibers starts significantly at 220 °C and ended at 500 °C under nitrogen. Comparatively, the degradation of neat PAN fibers starts at 250 °C and ceased at around 450 °C.
For a heating rate of 5 K·min −1 and at 800 °C under nitrogen, the char residue value of nonstabilized lignin-containing fibers was 56%. This high residue value is important from two points of view: (a) Neat lignin cannot maintain such values under comparable conditions. Usually, the char residue of lignin lies between 20 and 30% at 800 °C under nitrogen. (b) For CF, it is important to maintain high carbon yields to increase the productivity, as described in the introduction of this paper.
To further investigate the degradation activation energies from TGA measurements under nitrogen, peak temperatures of different degradation processes were collected at different heating rates. On the basis of these data, the activation energies of different processes were calculated by the Kissinger method from the following plot: ( ) The activation energy of the related process was then calculated from the slope by the following equation: where is the heating rate in K·min −1 , is the ideal gas constant R = 8.3144598 J·K −1 ·mol −1 , is the temperature at the peak maximum in K, From the intercept, the pre-exponential factor A can be calculated. The rate constant kTm at a defined temperature Tm can be calculated from the following equation: For comparison purposes, the activation energy was also calculated using the Flynn-Wall-Ozawa (FWO) method from the following plot: The activation energy of the related process was then calculated from the slope by the following equation: Only degradation peaks from DTG curves greater than 0.05%·K −1 of neat PAN and PAN/lignin fibers under nitrogen were considered for this evaluation (Figure 8c,d). The activation energies were listed in Table 5.
Neat PAN fibers showed three main degradation processes, as seen in the Figure 8c,d. This thermally induced degradation behavior is comparable to PAN fibers formed via the wet spinning technique, as described elsewhere [50,51]. Table 5. Activation energy of different degradation peaks of neat PAN fibers and PAN/lignin fibers under nitrogen atmosphere calculated from the Costs-Redfern method and the Flynn-Wall-Ozawa method; * error of the linear regression. For more details about the peak ID assignments, please refer to Supplementary Data Figure S3. For neat PAN fibers, the first degradation under nitrogen is very fast and had a maximum of around 1.5%·K −1 with the activation energy of 132 ± 16 kJ·mol −1 . In conjunction with the following DSC measurements, in this degradation zone, the conversion from nitrile groups to -C≡N to imine groups -C = N-could cyclize in the same chain with the adjacent nitrile group or crosslinks the polymer chains [52]. This reaction is followed by a consequent degradation shoulder with a lowered degradation rate between 0.75 and 1.00%·K −1 in the temperature range between 250 °C and 375 °C. The activation energy of this reaction was found to be 112 ± 15 kJ·mol −1 according to the Kissinger method. Intramolecular cross-linking reactions take place at this stage according to Rahaman et al. [52]. In the third degradation process, at higher temperatures, polymer chain scission reactions and gaseous volatile by-products, for example, methane, hydrogen, hydrogen cyanide, ammonia, and tarry molecules, could be ejected by the fibers according to the literature [52][53][54]. The highest activation energy was found to be at this stage and was 364 kJ·mol −1 according to the Flynn-Wall-Ozawa method. Introducing lignin in the PAN fibers increased the activation energy of the first (179 kJ·mol −1 ) and the second (211 kJ·mol −1 ) reaction processes and reduced the activation energy at elevated temperatures (268 kJ·mol −1 ). This is owing to the fact that softwood lignin is highly susceptible to any thermal processing, which initiates radical polymerization and increases the molecular weight owing to the cross-linking according to the available literature [18,19,55].

Shrinkage Measurements
During the stabilization process of the polyacrylonitrile fibers, two different structural changes take place in a clear sequence [54]. The first one is the entropic physical shrinkage. The origin of this shrinkage is mainly caused by the amorphous region of the polymer, which showed higher chain mobility than the crystalline region. Polyacrylonitrile showed single-hybrid-phase morphology [56,57]. For moderately drawn fibers, the morphology consists of ordered domains with their own glass transition region and amorphous disordered domains with a lower glass transition region [49,56,58,59]. The second part of the heat treatment induced shrinkage of the polyacrylonitrile fibers is a chemical shrinkage, which is mainly caused by the stabilization reactions and the conversion reaction of the nitrile groups -C≡N to imine groups -C = N-, which crosslinks the intra-molecular structure and forms a shorter ladder structure. The crosslinking reactions presume an oxidative environment [54]. Under a constant tensile load of 2.03 mN·tex −1 , the shrinkage of the fibers was recorded at different heating rates to understand the kinetics of the shrinkage process. The tensile load was generated in the case of neat PAN fibers by applying 0.01 N for 20 filaments of the mean fineness 4.92 tex and by 0.022 N for 20 filaments of the mean fineness 10.94 tex in the case of PAN/lignin fibers, respectively. Figure 9a,b illustrates the longitudinal shrinkage behavior of the fibers and the first derivative of the strain versus the temperature at different heating rates and at an adequate low specific tension force for both fibers. Neat PAN fibers show in general lower overall shrinkage compared with the lignin-containing fibers. At very low heating rate of 1 K·min −1 under air, the overall shrinkage was 20.8% and 23.5 for the neat PAN fibers and lignin-containing fibers, respectively. For the highly drawn PAN fibers, the chemical shrinkage showed a higher contribution to the overall shrinkage than the physical shrinkage. At a heating rate of 2 K·min −1 , the chemical shrinkage contributed 53.6% to the overall shrinkage. In contrast, the physical shrinkage of lignin-containing fibers contributed 58.1% (@ 1 K·min −1 ), 63.3% (@ 2 K·min −1 ), and 72.3% (@ 5 K·min −1 ) to the overall shrinkage. This phenomenon is owed to the fact that, the higher the heating rate, the lower the period given to a material to allow it to cross-link and oxidize during the heat treatment.
PAN fibers (DR = 10) showed an activation energy of the physical shrinkage of 92 ± 14 kJ·mol −1 (Figure 9c, Table 6). Introducing lignin to the fibers increased the activation energy of the physical shrinkage by around 3.14 times to 289 ± 28 kJ·mol −1 (Figure 9d, Table 6). From the crystallographic analysis, lignin decreased the crystalline domain size and increased the amorphous volume part in the fiber structure. This strengthened the assumption of a new amorphous phase of PAN with a higher volume percentage. Knowing that the activation energy is a specific value and represents per definition the minimum amount of energy required to initiate a process or reaction, some new amorphous phase physically bounded or crosslinked to lignin oligomers can be assumed to be responsible for such behavior. Lignin-containing fibers showed furthermore two associated chemical shrinkage phases (Figure 9b). The first chemical shrinkage phase (phase I) of PAN/lignin fibers showed a higher activation energy compared with that of PAN fibers, which were 147 ± 11 kJ·mol −1 and 135 ± 3 kJ·mol −1 , respectively. Phase I is mostly caused by conversion reactions of the nitrile groups -C≡N to -C = N-groups and the intra-chain crosslinking reactions in the ordered domains of the PAN phase. The activation energy of the second shrinkage phase (phase II) was 100 ± 7 kJ·mol −1 . This shrinkage is assumed to be caused by the cross-linking reactions of lignin. With the increase of the specific tension by 3.7 times, as seen from Figure 9e,f, a reduction of the overall shrinkage of the neat PAN fibers at a heating rate of 5 K·min −1 at 325 °C under air from 11.6% to 6.8% can be observed. Obviously, the peak maximum (from the first derivative of the strains versus the temperature) of the chemical shrinkage was found to be load-independent for the studied cases and was at 290 °C. The peak position of the physical shrinkage was shifted to a lower temperature. This change was very small for the studied load cases. Under a specific tension of 7.5 mN·tex −1 , PAN fibers started to be extended at 175 °C. This assumes that, at this temperature, the glass transition temperature of some amorphous part of the fibers was reached and the fibers became deformable at the applied tension load. This phenomenon was not detectable at lower tension loads (2.03 mN·tex −1 ). This extension took place over a range of 175 °C to 252 °C, at which the chemical shrinkage started to occur, as illustrated in Figure 9e.
In the case of lignin-containing fibers, however, the onset of glass transition was found at 163 °C and ended at 240 °C with a shift to the left side by roughly ~10 K compared with the neat PAN fibers under the same higher specific tension load. Additionally, the higher load remarkably reduced the peak maximum of the physical shrinkage (obtained from the first derivative of the strain) from 146 °C to 139 °C, as seen from Figure 9f. The chemical shrinkage of phase I occurred at a constant temperature at 266 °C for different loads and the first derivative of the strain versus the temperature was negative, meaning that shrinkage took place. The chemical shrinkage of phase II (lignin-induced) was also at 297 °C in both cases. Remarkably, at this temperature positive, values of the first derivative of the strain versus the temperature were recorded, which means that the fibers at this load level started to be stretched again. This could be initiated by the low molecular volatiles being gassed off by lignin molecules during the early pyrolysis stage (starting from 290 °C, as recorded by the TGA curves). These findings suggest for the future work a detailed evaluation of the activation energies of the physical and chemical shrinkage of the fibers at different specific tension loads and higher drawing ratios.

Thermal Stabilization Kinetics
In order to convert the precursor fibers into CF, a critical thermo-oxidative stabilization step is indispensable. In this step, the fibers must be converted into a thermally more stable state. The mechanisms of the thermo-oxidative stabilization of homo-polyacrylonitrile fibers and copolymers of the acrylonitrile were extensively studied in the literature [52,54,60,61].
In summary, the stabilization reactions of the PAN fibers can be described by conversion of the nitrile groups -C≡N to imine groups -C = N-, dehydration reactions of the PAN backbone chain, oxidation reactions (under air atmosphere), and crosslinking reactions, which take place in an overlapped fashion. Different volatiles, for example, H2O, CO2, CO, HCN, and NH3, were reported as a result of the stabilization reaction [52]. Figure 10 illustrates the DSC measurements of different heating rates for neat PAN fibers and PAN/lignin fibers conducted under air or nitrogen atmosphere at different heating rates. The measurements were conducted in order to reveal the influence of both oxygen (oxidative air atmosphere) and softwood lignin on the stabilization behavior of the developed fibers.
Under nitrogen (Figure 10a), neat PAN fibers showed single and sharp exothermal peaks of the conversion reaction of nitrile groups to imine groups. This peak shifted to higher temperatures by increasing the heating rate. Introducing lignin to the fibers (Figure 10b) resulted in the overall reduction of the released heat during the exothermal conversion reaction of nitrile groups to imine groups. An additional associated exothermal shoulder was also detected, suggesting exothermal crosslinking reactions induced by lignin. Furthermore, the onset and peak temperatures in the case of lignin-containing fibers shifted (compared with the neat PAN fibers) to lower temperatures by ~10 K for heating rates between 5 and 10 K·min −1 and by ~7 K for heating rates between 1 and 2 K·min −1 .
Neat PAN fibers showed, under oxidative atmosphere (air, Figure 10c), an additional exothermal oxidation shoulder followed by a subsequent exothermal crosslinking reaction peak. Increasing the heating rate reduced the FWHM of the oxidation and crosslinking shoulders and shifted the peaks to higher temperatures. In comparison with the neat fibers stabilized under nitrogen, the peak maxima of the conversion reactions of nitrile groups to imine groups were very slightly shifted to lower temperatures (Figure 10a,c). The addition of lignin reduced the overall released exothermal heat compared with that of neat PAN under the same conditions, suggesting that lignin oligomers actively acted as a stabilizer for the stabilization process.
Utilizing the Kissinger method, the activation energies and the kinetic parameters of the thermal stabilization processes of neat PAN fibers and lignin-containing fibers were calculated under air and nitrogen atmospheres and are summarized in Table 7. The activation energy of the conversion reactions of the nitrile groups -C≡N to imine groups -C = N-in the case of neat PAN fibers was 146 ± 2 kJ·mol −1 . This activation energy is comparable-within the linear regression error-to the activation energy of the neat PAN fibers under air, which was 151 ± 5 kJ·mol −1 . Under oxidative atmosphere, neat PAN fibers showed oxidation and crosslinking reactions with activation energies of 212 ± 29 kJ·mol −1 and 117 ± 7 kJ·mol −1 , respectively. Introducing lignin to the PAN fibers under oxidative atmosphere, the activation energy of the conversion reactions from nitrile to imine groups was increased by 12.30%. Furthermore, the rate constant calculated at 270 °C from the pre-exponential factor A ( Table 7) increased from 0.133 to 0.250 min −1 by ca. 88%. This implies that lignin oligomers accelerate the conversion reactions of nitrile to imine groups. Comparatively, higher activation energy (163 ± 1 kJ·mol −1 ) of the conversion reaction of nitrile to imine groups was detected for lignin-containing fibers under nitrogen. The rate constant was also increased from 0.124 min −1 (under nitrogen) to 0.232 min −1 . Similar behavior for homopolyacrylonitrile was reported by Heine et al. [60] at low heating rates up to 10 K·min −1 . The incorporation of lignin decreased the rate constant at 270 °C of the oxidation reaction by 43.1% from 0.072 to 0.041 min −1 . The activation energy of the oxidation process was also reduced by 34.1% to 139 ± 60 kJ·mol −1 . Under air, the constant rate of the crosslinking reaction in the case of the lignincontaining fibers increased by 30%, suggesting the supportive role of lignin to crosslink the structure of PAN. Under nitrogen atmosphere, the activation energy of the crosslinking reaction of PAN/lignin fibers was 236 ± 19 kJ·mol −1 . This activation energy was reduced by 90% to 124 ± 4 kJ·mol −1 for the same fibers under oxidative atmosphere. The constant rate of this reaction was increased 2.43 times    Figure 11a,b showed the EDX curves of the cross section of the fibers embedded in an epoxy matrix. High resolution images of the SEM of the cross section and the EDX images of the elements nitrogen, chlorine, oxygen, and sulphur are provided in Supplementary Data Figure S3.

Scanning Electron Microscopy (SEM) and Energy-Dispersive X-Ray Spectroscopy (EDX)
The evidence threshold of the atoms lies at 0.1 at. %. Chlorine could not be detected by this method, suggesting the high regeneration rate of the ionic liquids from the fiber structure.
The quantitatively detected portions of C/N/O/S/Cl were 66.17:28.07:5.72:0.03:0.01 in at. %, respectively. The amount of detected chlorine traces was around 0.01 at. %, suggesting a high regeneration rate of the ionic liquids out of the fiber structure. In the case of PAN/lignin fibers, more sulfur (0.16 at. %) was detected (Figure 11b) owing to the chemically bounded sulfur to the lignin caused by the early pulping process. The quantitatively detected portions of C/N/O/S/Cl were 64.99:22.31:12.53:0.16:0.00 in wt. %, respectively. SEM images (Figure 11c,d) of the fiber surfaces showed, in both cases (PAN and PAN/lignin fibers), smooth and homogenous surface with longitudinal grooves caused by the flow lines in the nozzle during the extrusion process. Neither pores over the cross section of the fibers could be detected nor did a lobulated surface occur, suggesting a controllable and homogeneous coagulation process.

Conclusions
In the present study, we explored the shear rheological flow behavior, as well as the miscibility, compatibility, and rheological determined thermal stability, of the mixture of PAN, purified softwood lignin, and 1-ethyl-3-methylimidazolium chloride. All investigated dope systems were found to be miscible in term of the Arrhenius plot and compatible for different lignin content in terms of the temperature-invariant Han plot.
The developed fiber based on PAN and lignin offered homogeneous and smooth-surface fibers with fineness values ~6 dtex, tensile strengths ~616 MPa, tensile moduli ~30 GPa, as well as elongation at break ~13.27%.
The coagulation process and coagulation bath temperature extensively influenced the fine crystallographic structure of the fibers. Lowered crystal size, orientation degree, and Herman's orientation factors were a result of increased coagulation temperatures. No spinning was possible at high coagulation temperatures, that is, at 40 °C. This was found to be correlated to the temperaturedependent diffusion rate of the ionic liquids from the fiber structure in the coagulation bath and to the temperature-dependent relaxation times of the spinning dope. At low coagulation temperatures of 6 °C, lignin was found to deteriorate the crystal size of the PAN phase in the fiber structure while maintaining the dimensions of the unit cell of the PAN phase. Furthermore, introducing the amorphous low molecular weight lignin oligomers to the fibers reduced the orientation degree and Herman's orientation factors at a given coagulation bath temperature and draw ratio. The essential thermal properties of the developed fibers were investigated. The TGA analysis of the regenerated fibers under inert atmosphere showed unusual and better thermal stability of the lignin-containing fibers than the neat PAN fibers and powders as well as neat lignin powder. This could also be detected at the lowered DTG peaks under the same conditions (heat rate, atmosphere, and temperature range) in comparison with the neat PAN fibers. Lignin increased the activation energy of the characteristic first and second decomposition peaks and lowered the activation energy of the thermal degradation at higher temperatures. The thermo-mechanical shrinkage tests under air supported the findings in the crystallographic fine structure of the fibers. Lignin increased the amorphous part in the fiber structure and introduced a physical interaction with the amorphous part. The existence of a single glass transition of the fibers with a slight positive shift was found in the case of the PAN fibers. Additional research work will be conducted in the future to understand the impact of lignin on the domain length of the PAN phase in the neat PAN fibers and lignin-containing fibers using the small angle X-ray scattering technique. In order to investigate the degree of crystallinity, the ultrafast DSC technique is required, as shown by Furushima et al. in [62]. The crystal conversion mechanisms, orientation recovery mechanisms, and structure formation of the subsequent CF on the base of the developed fibers will also be pursued. It can be also important to mention that no corrosion reactions (caused by the ionic liquids) were found to take place, neither for the measuring geometry of the ARES G2 Rheometer during the long-term rheological characterization nor during the spinning trails.

Supplementary Materials:
The following are available online at www.mdpi.com/1996-1944/13/17/3687/s1, Figure S1: Frequency sweep tests of (a) dope ND; (b) dope L1; (c) dope L2 at different temperatures., Figure S2: (a) Plot of the throughput of the up-scaled spinning dopes ND and L1 versus the rotational speed of the spinning pump; (b) the spinnability chart of the spinning dope ND and (c) the spinnability chart of the spinning dope L1 (round green circles mean the system was found to be spinnable; and round tacked red circles for un-spinnable set-up. Figure S3: DTG curves of (a) PAN fibers and (b) PAN/lignin fibers (DR =2) for detailed assignment of the Peak ID 1, 2 and 3 for both fiber types. Figure