Modeling and Simulations of the Sulfur Inﬁltration in Activated Carbon Fabrics during Composite Cathode Fabrication for Lithium-Sulfur Batteries

: During the manufacture of a composite cathode for lithium-sulfur (Li-S) batteries it is important to realize homogeneous inﬁltration of a speciﬁed amount of sulfur, targeted to be at least 5 mg cm − 2 to achieve good battery performance in terms of high energy density. A model of the sulfur inﬁltration is presented in this study, taking into account the pore size distribution of the porous cathode host, phase transitions in sulfur, and formation of different sulfur allotropes, depending on pore size, formation energy and available thermal energy. Simulations of sulfur inﬁltration into an activated carbon fabric at a hot-plate temperature of 175 ◦ C for two hours predicted a composite cathode with 41 wt% sulfur (8.3 mg cm − 2 ), in excellent agreement with the experiment. The pore size distribution of the porous carbon host proved critical for both the extent and form of retained sulfur, where pores below 0.4 nm could not accommodate any sulfur, pores between 0.4 and 0.7 nm retained S 4 and S 6 allotropes, and pores between 0.7 and 1.5 nm contained S 8 .


Introduction
Lithium-sulfur (Li-S) batteries are one of the latest generation lithium batteries, pursued due to their high theoretical energy density, 2510 Wh kg −1 [1], compared to an expected maximum of 250 Wh kg −1 for conventional Li-ion batteries [2][3][4]. Li-S batteries typically comprise a lithium anode and a composite sulfur cathode, consisting of sulfur as the active cathode material and a porous host, which provides electronic conductivity and structural integrity. Such hosts include activated carbon fabrics (ACF) [5], hollow nanofibers [6], hierarchically porous carbon coatings [7][8][9][10], graphene, and graphene oxide [11,12], where oxygen groups and functional, transition metal oxide groups [9,10] trap sulfur and polysulfides and thus reduce the "shuttling" effect during battery operation, and highly conductive materials such as carbon black and multiwall carbon nanotubes [13,14]. Furthermore, the porous carbon cathode host contributes a supercapacitor element [12,[14][15][16][17][18][19][20][21][22][23], hence, increasing the power density of the Li-S battery. In this manner, a Li-S battery with such composite cathode is a hybridized battery/Li-ion capacitor [24], combining the high energy density of the Li-S battery and the high power density of the Li-ion capacitor (consisting of the lithium anode and porous carbon cathode), similarly to composite supercapacitors [25,26] but with much higher energy density, as the theoretical capacity of a Li-S battery reaches 1670 Ah g s −1 (where g s −1 denotes per gram of sulfur). A Li-S battery with a composite cathode can be easily recycled by disassembly [27] and liquid processing, dissolution of binder and sulfur [28,29], dielectrophoresis [30], and filtration techniques.
In order for the Li-S battery to achieve the expected high capacity and energy density, the amount of sulfur needs to be maximized in the composite cathode, recently specified above 5 mg cm −2 of cathode [31], with at least 70 wt% sulfur in the composite cathode

Materials and Experimental Methods
The cathode host was ACF Kynol ® ACC-507-15 (Kynol Europe GmbH (Hamburg, Germany): areal density of 12 mg cm −2 and specific surface area BET = 1461 m 2 g −1 ), presented in Figure 1a,b, with a discretized pore size distribution (PSD) (Figure 1c) derived from data of nitrogen adsorption experiments [24]. The PSD started at 0.3 nm and exhibited peaks at 0.69 nm and 1.31 nm, which were in the micropore range. Micropores were expected to be of slit shape in phenolic-derived ACFs, as was the Kynol ® fabric [24,38]. Figure 1b displays that macro-and mesopores of approximately circular cross section were also present and, during sulfur infiltration or battery operation, provided the transport channels for the transported species, sulfur, or ions [24], respectively, to reach the large proportion of micropores that dominated the PSD of the ACF. Sulfur powder (Sigma Aldrich, Gillingham, UK) was used as the active cathode material, via an infiltration stage into the porous ACF host. channels for the transported species, sulfur, or ions [24], respectively, to reach the large proportion of micropores that dominated the PSD of the ACF. Sulfur powder (Sigma Aldrich, Gillingham, UK) was used as the active cathode material, via an infiltration stage into the porous ACF host. ACF samples of 1 cm 2 (12 mg) were placed on 1-cm 2 aluminum foil (of 30-µm thickness, utilized as current collector in the battery cathode). Twelve mg of sulfur powder was sprinkled on the free top surface of each ACF sample. The assembly was wrapped by a single layer of aluminum foil (to entrap any sulfur vapor as much as possible), taking care for the wrap not to touch the sulfur powder. Finally, each wrapped assembly was placed in a sealed pressure vessel from stainless steel of a circular base of 160-mm diameter and height of 110 mm, as illustrated in Figure 2a. The pressure vessel was placed on a hot plate for a specified period of time. Different experiments were carried out at three different hot-plate temperatures: Tb = 150, 175, and 200 °C. External temperature monitoring was carried out during each experiment using an infrared (IR) camera ( Figure 2b). On the basis of the ideal gas law, a gauge pressure range of 0-0.6 bar was estimated, where the safety valve in the pressure vessel was set to 1 bar and was never reached, i.e., the safety valve did not open in any of the experiments. It must be mentioned that the dimensions of the pressure vessel were designed so if its base allocated a maximum number of eight samples, the estimated maximum gauge pressure was below the safety valve setting of 1 bar. After 40 min and after 2 h at Tb, followed by 12 h cooling to room temperature, the samples were removed from the vessel; the 40-and 120-min experiments were independent experiments, where the 2-h experiments were not interrupted. Each sample on the current collector foil was weighed before and after sulfur powder was added, as well as at the end of the specified period of sulfur infiltration processing. Samples were examined under scanning electron microscope HR-SEM JEOL-7100 F (JEOL Ltd., Tokyo, Japan) in terms of SEM and EDX (energy-dispersive X-ray spectroscopy). ACF samples of 1 cm 2 (12 mg) were placed on 1-cm 2 aluminum foil (of 30-µm thickness, utilized as current collector in the battery cathode). Twelve mg of sulfur powder was sprinkled on the free top surface of each ACF sample. The assembly was wrapped by a single layer of aluminum foil (to entrap any sulfur vapor as much as possible), taking care for the wrap not to touch the sulfur powder. Finally, each wrapped assembly was placed in a sealed pressure vessel from stainless steel of a circular base of 160-mm diameter and height of 110 mm, as illustrated in Figure 2a. The pressure vessel was placed on a hot plate for a specified period of time. Different experiments were carried out at three different hot-plate temperatures: T b = 150, 175, and 200 • C. External temperature monitoring was carried out during each experiment using an infrared (IR) camera ( Figure 2b). On the basis of the ideal gas law, a gauge pressure range of 0-0.6 bar was estimated, where the safety valve in the pressure vessel was set to 1 bar and was never reached, i.e., the safety valve did not open in any of the experiments. It must be mentioned that the dimensions of the pressure vessel were designed so if its base allocated a maximum number of eight samples, the estimated maximum gauge pressure was below the safety valve setting of 1 bar. After 40 min and after 2 h at T b , followed by 12 h cooling to room temperature, the samples were removed from the vessel; the 40-and 120-min experiments were independent experiments, where the 2-h experiments were not interrupted. Each sample on the current collector foil was weighed before and after sulfur powder was added, as well as at the end of the specified period of sulfur infiltration processing. Samples were examined under scanning electron microscope HR-SEM JEOL-7100 F (JEOL Ltd., Tokyo, Japan) in terms of SEM and EDX (energy-dispersive X-ray spectroscopy).

Modeling of the Sulfur Infiltration in the Porous Cathode Host
A numerical model was developed for the heat transfer and sulfur impregnation in the cathode host. On the basis of the diagram presented in Figure 2a, one-dimensional equations were assembled in the y-axis, given the small thickness of each component in the assembly.
All possible sulfur allotropes were considered for melt or low-temperature vapor [36] and the lowest energy structures were selected, as predicted by density functional theory (DFT) simulations by Jackson et al. [39]. Further molecular simulations were carried out in this study for geometrical optimization, using Avogadro, and the resulted molecular models are displayed in Figure 3a, together with the measured minimum dimension for each allotrope model and their formation energy, as calculated by Jackson et al. [39] employing DFT simulations. On the basis of the minimum dimension of each sulfur allotrope, it was deducted which allotrope may be accommodated in each pore size, as is presented for the peaks of the PSD of the porous cathode host ACF in Figure 3b. Pores smaller than 0.4 nm could not accommodate any sulfur in any form and pores smaller than 0.7 nm could not accommodate S8 and S7.

Modeling of the Sulfur Infiltration in the Porous Cathode Host
A numerical model was developed for the heat transfer and sulfur impregnation in the cathode host. On the basis of the diagram presented in Figure 2a, one-dimensional equations were assembled in the y-axis, given the small thickness of each component in the assembly.
All possible sulfur allotropes were considered for melt or low-temperature vapor [36] and the lowest energy structures were selected, as predicted by density functional theory (DFT) simulations by Jackson et al. [39]. Further molecular simulations were carried out in this study for geometrical optimization, using Avogadro, and the resulted molecular models are displayed in Figure 3a, together with the measured minimum dimension for each allotrope model and their formation energy, as calculated by Jackson et al. [39] employing DFT simulations. On the basis of the minimum dimension of each sulfur allotrope, it was deducted which allotrope may be accommodated in each pore size, as is presented for the peaks of the PSD of the porous cathode host ACF in Figure 3b. Pores smaller than 0.4 nm could not accommodate any sulfur in any form and pores smaller than 0.7 nm could not accommodate S 8 and S 7 .

Modeling of the Sulfur Infiltration in the Porous Cathode Host
A numerical model was developed for the heat transfer and sulfur impregnation in the cathode host. On the basis of the diagram presented in Figure 2a, one-dimensional equations were assembled in the y-axis, given the small thickness of each component in the assembly.
All possible sulfur allotropes were considered for melt or low-temperature vapor [36] and the lowest energy structures were selected, as predicted by density functional theory (DFT) simulations by Jackson et al. [39]. Further molecular simulations were carried out in this study for geometrical optimization, using Avogadro, and the resulted molecular models are displayed in Figure 3a, together with the measured minimum dimension for each allotrope model and their formation energy, as calculated by Jackson et al. [39] employing DFT simulations. On the basis of the minimum dimension of each sulfur allotrope, it was deducted which allotrope may be accommodated in each pore size, as is presented for the peaks of the PSD of the porous cathode host ACF in Figure 3b. Pores smaller than 0.4 nm could not accommodate any sulfur in any form and pores smaller than 0.7 nm could not accommodate S8 and S7.
of thickness of 30 µm; (j) air (at temperature T (K) and pressure P (Pa): The heat transfer equation for each single material layer, i, is given by: where T is the temperature, t is the time, y is the vertical axis, ρ is density, c p is specific heat capacity, k is thermal conductivity, subscript i denotes layer i, and ∑ j R j DH f,Sj is the power lost due to the formation of new sulfur allotropes j. The air pressure in the container was updated according to the updated air temperature, employing the ideal gas law. The specific heat capacity of sulfur, c p,S , is given as a function of the absolute temperature, T (in K), by the following relation with the values of the constants given in Table 1 [40]. For the composite ACF-S, the thermal diffusivity, k i ρ i c p,i , is given by the rule of mixtures, where the thermal diffusivity of the sulfur is weighted by the volume fraction of the impregnated sulfur in the ACF. At the boundaries between adjacent layers, there is heat flux continuity boundary condition, with a contact resistance, 1/h c , to take into account any non-perfect contact at the ACF/current collector interface and previous interfaces (layers (b)-(e)). A fixed boundary condition, T b , is used at the bottom of layer (b), in contact with the hot plate. The time-changing temperature at the ceiling of the pressure container is given by the temperature at that height from the IR images, monitoring the experiment, and used as input data in the heat transfer simulation. A heat transfer coefficient of 390 W m −2 K −1 was used for the heat transfer via natural convection between the S-ACF and air [41,42] and of 2280 W m −2 K −1 between air and aluminum foil or stainless steel top plate [43].
At the sulfur melt-solid interface, with sulfur melt point at 119.7 • C, there is the boundary condition: where ∂S m ∂t is the rate of sulfur melt generation, ∆H m is the heat of fusion of sulfur, ∆H m = 54 kJ kg −1 [40], and, for the rest of properties, subscripts s and m denote solid and melt, respectively.
Sublimation (solid-to-vapor for temperatures below 119.7 • C) and vaporization of sulfur (melt-to-vapor for temperatures above 119.7 • C) occur from the top surface of sulfur layer (g) ((g) top subscript), where a heat flux loss occurs at that boundary with the air layer boundary (h), according to the relation: where the heat of vaporization is ∆H S,v = 306.25 kJ kg −1 and the heat of sublimation is ∆H S,v = 360.31 kJ kg −1 [40]. The rate of generation of S vapor sulfur moles, dN S,v /dt, is given by the equation [44,45]: where A is the area of the top surface of layer (g), U transfer is the mass transfer coefficient, P S,v is the sulfur vapor pressure and R is the gas constant R = 8.3145 J mol −1 K −1 . The mass transfer coefficient is given by the relation [45]: (6) where g is gravity (g = 9.81 m s −2 ), µ S,v is the sulfur vapor viscosity given by the linear relation against temperature fitting the experimental data from [46], as presented in Figure 4a, and D S,v is the diffusion coefficient of the sulfur vapor, given as a function of temperature, T (in K), by the following relation [44], assuming binary coefficients with air being the second vapor-phase material [44]:  The sulfur vapor pressure is considered a function of temperature from the polynomial data fit of experimental data from [47], as is presented in Figure 4b.
The second system of equations describes the mass transport of each sulfur allotrope, Sj for j = 2-8 (Figure 3a), and infiltration through the porous ACF in the average y-axis and in each pore size, p, from the discretized PSD in Figure 1c: where Cj,p is the concentration of sulfur allotrope j in pore size p, up is the velocity in pore size p, Dj,p is the diffusion coefficient of sulfur allotrope j in pore size p, Rj is the rate of formation of new allotrope Sj in pore size dp, rj is the rate of loss of allotrope Sj due to the formation of lower order allotropes for smaller pores, and Fdecay,p is a decay factor where Fdecay,j,p = 1 for the S8 allotrope (j = 8) but depends on local thermal energy and allotrope formation energy for the lower order sulfur allotropes. Parallel infiltration paths are considered for the different pore sizes, with the mean velocity, up, for each pore size calculated according to Darcy's law in the same manner as matrix infiltration in a fibrous medium in composites' manufacturing [37,48]: where εp is the pore fraction for pores sized p, Lf is the filled length across the thickness of the porous ACF, Kp is the permeability of the porous medium for pore size dp given by the Carman-Kozeny equation for flow transverse to fibers with Kozeny constant of 3.4 [49]: Pc,p is the capillary pressure given by the Young-Laplace equation [37,50]: where σ is the surface tension of the infiltrating sulfur fluid, θ is its contact angle with the carbon pore walls (σ = 61 mN m −1 and θ = 4.3° [51] in this study), and µS,m is the viscosity of the sulfur melt, its value taken from experimental data depending on temperature [33] as:  The sulfur vapor pressure is considered a function of temperature from the polynomial data fit of experimental data from [47], as is presented in Figure 4b.
The second system of equations describes the mass transport of each sulfur allotrope, S j for j = 2-8 (Figure 3a), and infiltration through the porous ACF in the average y-axis and in each pore size, p, from the discretized PSD in Figure 1c: where C j,p is the concentration of sulfur allotrope j in pore size p, u p is the velocity in pore size p, D j,p is the diffusion coefficient of sulfur allotrope j in pore size p, R j is the rate of formation of new allotrope S j in pore size d p , r j is the rate of loss of allotrope S j due to the formation of lower order allotropes for smaller pores, and F decay,p is a decay factor where F decay,j,p = 1 for the S 8 allotrope (j = 8) but depends on local thermal energy and allotrope formation energy for the lower order sulfur allotropes. Parallel infiltration paths are considered for the different pore sizes, with the mean velocity, u p , for each pore size calculated according to Darcy's law in the same manner as matrix infiltration in a fibrous medium in composites' manufacturing [37,48]: where ε p is the pore fraction for pores sized p, L f is the filled length across the thickness of the porous ACF, K p is the permeability of the porous medium for pore size d p given by the Carman-Kozeny equation for flow transverse to fibers with Kozeny constant of 3.4 [49]: P c,p is the capillary pressure given by the Young-Laplace equation [37,50]: where σ is the surface tension of the infiltrating sulfur fluid, θ is its contact angle with the carbon pore walls (σ = 61 mN m −1 and θ = 4.3 • [51] in this study), and µ S,m is the viscosity of the sulfur melt, its value taken from experimental data depending on temperature [33] as: Sulfur allotropes are generated (R j ) and infiltrate the pores in pore size hierarchy, as follows. All pores with d p ≥ L min,S8 are infiltrated by the S 8 allotrope. Pores with L min,S8 > d p ≥ L min,S7 are first given the chance to be infiltrated with the S 7 allotrope, if there is sufficient local thermal energy, E Th , available to overcome the formation energy of S 7 , DH f,S7 . Otherwise, formation and infiltration of S 6 is considered and consecutive smaller allotropes, pores with L min,S7 > d p ≥ L min,S6, are first given the chance to be infiltrated with the S 6 allotrope, if there is sufficient local thermal energy available to overcome the formation energy of S 6 , DH f,S6 . Otherwise, formation and infiltration of S 5 is considered and consecutive smaller allotropes. The procedure continues in the same manner until pores of dp < L min,S2 are encountered, which cannot be infiltrated by any form of sulfur. The rate of formation of a sulfur allotrope, R j , encompasses the full concentration of the parent allotrope (Figure 3a) for that pore size, multiplied by the decay factor, F decay,j,p , given by the equation [24]: The diffusion coefficient of sulfur allotrope j in pore size d p is given by the following relation [17] (where k B is Boltzmann's constant), which is a modified Einstein's equation to take into account tortuosity and constrictivity effects [24]: The system of the heat transfer and sulfur transport equations was solved using the implicit time-dependent finite volume/finite difference technique [52,53]. Computer simulations were performed for the experimental process described in Section 2, at T b = 150 • C or 175 • C or 200 • C for two alternative, specified infiltration durations: 40 min or 2 h, employing a timestep ∆t = 0.1 s and a y-increment ∆y = 5 µm in layers (a)-(i).

Results
The monitoring IR camera showed that after the hot plate reached the set temperature and the steel container with the samples was placed on the hot plate, the hot-plate temperature was disturbed and it took about 30-40 min to reach the set temperature again. From then on, the samples were left for sulfur infiltration for either 40 min or 2 h. To ensure that the correct experimental temperature profile was used, both the bottom and top temperatures of the steel container were inputted in the computer simulations from the monitoring IR camera. For hot plate-temperature T b = 150 • C, the sulfur powder did not melt and was still on top of the ACF samples when the steel container was opened, after either 40 min or 2 h. The computational simulations predicted an end, steady-state temperature of 109.88 • C at the top of the ACF sample, which was below the sulfur melting temperature, which justified the experimental result.
The computational simulations predicted an end, steady-state temperature of 125.55 • C and 146.8 • C at the top of the sulfur layer for initially set, experimental hot-plate temperatures of 175 and 200 • C, respectively. These steady state temperatures at the top of ACF were then above the sulfur melting point. Figure 5 displays the predicted total mass of sulfur (in the form of any allotrope) that infiltrated the porous ACF and also and in top sulfur layer as a function of the infiltration time, for the two, set, hot-plate temperatures, 175 and 200 • C. It also compares the predictions with the experimental data of the weight gain after sulfur infiltration of 40 min and 2 h. Initially, the sulfur was approximately 12 mg cm −2 , same as the amount of sulfur added initially. Sulfur evaporation or sublimation seemed to start slowly and then accelerated as the sample was heated and continued until the sulfur started impregnating meso-and micropores, where the numerical model considered that it was trapped without vapor escape from the ACF. In the simulations, sulfur seemed to impregnate and stay in the macropores for longer in the samples at set, hot-plate temperature, T b = 175 • C, from which it may have evaporated faster due to the larger surface area of macropores, compared with the assumed flat surface of the sulfur layer at the top of ACF. Hence, there was faster initial sulfur evaporation for T b = 175 • C, compared to the samples at T b = 200 • C. However, overall, more sulfur evaporated until the end for the samples at T b = 200 • C, which was attributed to the higher temperature resulting in higher vapor pressure, as seen in Figure 4b. As a result, the set hot-plate temperature T b = 175 • C is recommended, provided that the sulfur is distributed homogeneously through the thickness of the ACF sample and there is no accumulated sulfur layer at the top or bottom surface on the sample. This was investigated, both in the predictions and experimentally, in the rest of this study. added initially. Sulfur evaporation or sublimation seemed to start slowly and then accelerated as the sample was heated and continued until the sulfur started impregnating meso-and micropores, where the numerical model considered that it was trapped without vapor escape from the ACF. In the simulations, sulfur seemed to impregnate and stay in the macropores for longer in the samples at set, hot-plate temperature, Tb = 175 °C, from which it may have evaporated faster due to the larger surface area of macropores, compared with the assumed flat surface of the sulfur layer at the top of ACF. Hence, there was faster initial sulfur evaporation for Tb = 175 °C, compared to the samples at Tb = 200 °C. However, overall, more sulfur evaporated until the end for the samples at Tb = 200 °C, which was attributed to the higher temperature resulting in higher vapor pressure, as seen in Figure 4b. As a result, the set hot-plate temperature Tb = 175 °C is recommended, provided that the sulfur is distributed homogeneously through the thickness of the ACF sample and there is no accumulated sulfur layer at the top or bottom surface on the sample. This was investigated, both in the predictions and experimentally, in the rest of this study.   Figure 8, where S, C, and mixed S/C element maps depict that there was no sulfur layer surrounding the carbon fiber and, in fact, the sulfur had higher concentration in the central region of the fiber drawn by capillary flow into the micropores there. Any sulfur and carbon in the background of the fiber in Figure 8 was from the neighboring fibers, as shown in Figure 7b.  Figure 7 are from samples cut after their sulfur impregnation. Figure 7 illustrates that sulfur impregnated homogeneously across the whole thickness of the sample, without any accumulation at the top or the bottom of the sample. No leftover, unimpregnated sulfur was detected in any of the case studies: T b = 200 or 175 • C, impregnation time 40 or 120 min. The cross section of one of the cut fibers of Figure 7b is presented in Figure 8, where S, C, and mixed S/C element maps depict that there was no sulfur layer surrounding the carbon fiber and, in fact, the sulfur had higher concentration in the central region of the fiber drawn by capillary flow into the micropores there. Any sulfur and carbon in the background of the fiber in Figure 8 was from the neighboring fibers, as shown in Figure 7b.      Figure 9 presents the predicted evolution of different sulfur allotropes near the bottom of the ACF. Sulfur infiltration took place in a hierarchical manner, where the large pores were infiltrated first, followed by mesopores and, lastly, by micropores. S8 entered rapidly and hierarchically the large pores that could accommodate it, i.e., dp > 0.7 nm. For pores smaller than 0.7 nm, S8 converted to S6, and then to S4 for dp < 0.65 nm.  Figure 10 presents the predicted contour plots of the sulfur concentration at the end of the 40 min and the 2 h infiltration periods, for set hot-plate temperature Tb = 175 °C . Micropores greater than 0.7 nm were impregnated with S8. Pores smaller than 0.8 nm and larger than 0.6 nm contained S6, and pores smaller than 0.64 nm and larger than 0.4 nm contained S4. The sulfur distribution was homogeneous through the ACF thickness at the end of the 2 h infiltration period. After 40 min sulfur infiltration, there was some low S8 concentration in the meso-and macropores near the top ACF surface, which was depleted in the subsequent infiltration time for the sulfur infiltration into the lower parts, while the initial top reservoir of sulfur over the ACF was exhausted. At the end of the 40 min infiltration period, S4 had not fully infiltrated yet the bottom 20% of the ACF but the S4 infiltration was completed through the whole ACF thickness at the end of the 2-h period.   Figure 9 presents the predicted evolution of different sulfur allotropes near the bottom of the ACF. Sulfur infiltration took place in a hierarchical manner, where the large pores were infiltrated first, followed by mesopores and, lastly, by micropores. S8 entered rapidly and hierarchically the large pores that could accommodate it, i.e., dp > 0.7 nm. For pores smaller than 0.7 nm, S8 converted to S6, and then to S4 for dp < 0.65 nm.  Figure 10 presents the predicted contour plots of the sulfur concentration at the end of the 40 min and the 2 h infiltration periods, for set hot-plate temperature Tb = 175 °C. Micropores greater than 0.7 nm were impregnated with S8. Pores smaller than 0.8 nm and larger than 0.6 nm contained S6, and pores smaller than 0.64 nm and larger than 0.4 nm contained S4. The sulfur distribution was homogeneous through the ACF thickness at the end of the 2 h infiltration period. After 40 min sulfur infiltration, there was some low S8 concentration in the meso-and macropores near the top ACF surface, which was depleted in the subsequent infiltration time for the sulfur infiltration into the lower parts, while the initial top reservoir of sulfur over the ACF was exhausted. At the end of the 40 min infiltration period, S4 had not fully infiltrated yet the bottom 20% of the ACF but the S4 infiltration was completed through the whole ACF thickness at the end of the 2-h period.  Figure 10 presents the predicted contour plots of the sulfur concentration at the end of the 40 min and the 2 h infiltration periods, for set hot-plate temperature T b = 175 • C. Micropores greater than 0.7 nm were impregnated with S 8 . Pores smaller than 0.8 nm and larger than 0.6 nm contained S 6, and pores smaller than 0.64 nm and larger than 0.4 nm contained S 4 . The sulfur distribution was homogeneous through the ACF thickness at the end of the 2 h infiltration period. After 40 min sulfur infiltration, there was some low S 8 concentration in the meso-and macropores near the top ACF surface, which was depleted in the subsequent infiltration time for the sulfur infiltration into the lower parts, while the initial top reservoir of sulfur over the ACF was exhausted. At the end of the 40 min infiltration period, S 4 had not fully infiltrated yet the bottom 20% of the ACF but the S 4 infiltration was completed through the whole ACF thickness at the end of the 2-h period. mpos. Sci. 2021, 5, x FOR PEER REVIEW 13 of 17

Discussion
It is most interesting that in both the experiment and the computational simulation a sulfur loss of 59-62 wt% occurred, depending on hot-plate temperature. Sulfur vapor may have escaped outside the aluminum wrap, as after several such experiments the inner walls of the stainless steel pressure vessel obtained a shiny, yellowish hue. The heat transfer characteristics of the experimental rig greatly affected the temperature distribution and, hence, sulfur vaporization and also infiltration. Alternative experimental setups in the literature include sulfur impregnation with the sample in a convection oven [12] where, in fact, even higher degree of sulfur vaporization was observed [12]. The methodology of a sealed container on a hot plate with set hot-plate temperature, employed in this study, may represent continuous, large-scale sulfur deposition on a substrate tape moving

Discussion
It is most interesting that in both the experiment and the computational simulation a sulfur loss of 59-62 wt% occurred, depending on hot-plate temperature. Sulfur vapor may have escaped outside the aluminum wrap, as after several such experiments the inner walls of the stainless steel pressure vessel obtained a shiny, yellowish hue. The heat transfer characteristics of the experimental rig greatly affected the temperature distribution and, hence, sulfur vaporization and also infiltration. Alternative experimental setups in the literature include sulfur impregnation with the sample in a convection oven [12] where, in fact, even higher degree of sulfur vaporization was observed [12]. The methodology of a sealed container on a hot plate with set hot-plate temperature, employed in this study, may represent continuous, large-scale sulfur deposition on a substrate tape moving continually on a conveyor belt. The substrate may be ACF on aluminum foil current collector, their assembly moving on the conveyor belt with sulfur powder being sprayed on top and the bed of the moving belt being heated to the set temperature. Hence, the heat transfer and sulfur infiltration model presented in this study is applicable to a scaled fabrication of composite S-ACF cathode. As sulfur loss through vaporization and sublimation is generally undesirable, process design simulations may try to optimize the temperature and external pressure to maintain sulfur infiltration in melt state rather than vapor state. This would mean that only S 8 allotrope would be present, as lower-order sulfur allotropes are formed in vapor state only. As a result, cathode host micropores below 0.7 nm would not contain any sulfur and such unfilled pores would, instead, absorb and waste the electrolyte solvent of the Li-S battery without any prospect of active participation in the redox battery reactions in the absence of sulfur.
It was predicted that in pores smaller than 0.7 nm (and greater than 0.4 nm) the sulfur resides in the form of mainly S 6 and S 4 allotropes. The so confined sulfur in small carbon micropores cannot move to the anode so easily in a Li-S battery, which offers the advantage of limiting the "shuttling" effect in the Li-S battery [54].
However, S 4 in micropores smaller than 0.64 nm, which occupy 16% of the total volume of pores sized 0.46-119 nm in the ACF of this study, reduces the capacity of the first potential plateau at 2.3-2.4 V in the discharge of a Li-S battery for some composite cathodes, as this plateau corresponds to the conversion of Li 2 S 8 to Li 2 S 6 [5,54]. Therefore, it can be seen that the pore size distribution of the porous carbon host is critical for both the amount and form of infiltrated sulfur, which are going to affect the capacity of the Li-S battery, depending on the amount and allotrope type of sulfur and its energy stored and delivered, which depends on the integrated product of potential and capacity.
Hence, the design of future cathode hosts for Li-S batteries is a complex exercise. Micropores smaller than 0.7 nm are detrimental for the Li-S battery energy, as they can only accommodate lower-order sulfides that reduce the first voltage plateau in the battery while these micropores consume electrolyte, which is sometimes a waste as the absorbed electrolyte remains passive if there is no sufficient sulfur to utilize it in redox reactions. However, these pores filled with the electrolyte offer a high-power supercapacitor component and the device is then a hybrid Li-S battery/Li-ion capacitor. Furthermore, larger micropores might be beneficial to the Li-S battery as they can accommodate S 8 and restrict escape of polysulfides, reducing the "shuttling" effect during battery cycling.
The numerical model of heat transfer and sulfur infiltration outlined in this paper may also be applied to other porous cathode hosts, such as activated carbon powder [7,9,10] graphene [11,12] or other porous carbon material [6] with measured pore size distribution that can be used as input data in the simulations. The material properties and porosity would be used to determine the thermal properties (heat capacity and thermal conductivity) for the heat transfer, Equation (1). The pore size distribution of the porous host controls the sulfur impregnation in Equations (8)- (11) and the type of sulfur allotrope that may be formed in each pore size. In general, simulations of the sulfur infiltration may comprise the first stage, in a comprehensive numerical model of cathode fabrication and Li-S battery operation. Combining simulations of the composite cathode fabrication and battery cycling will help in the design of optimized porous hosts and composite cathodes.

Conclusions
The pore size distribution of a porous carbon host for a composite cathode is very important in a Li-S battery, both in cathode manufacture and during battery operation. A numerical model has been developed in this study for the description of the sulfur infiltration into the porous host while heated above the melting point of sulfur in cathode manufacture, taking into account the pore size distribution of the porous host and including heat transfer, sulfur sublimation, and evaporation, as well as formation of lower-order sulfur allotropes, depending on pore size and available thermal energy compared to their energy of formation. Computational simulations predicted a composite cathode with 38 and 41 wt% sulfur for set hot-plate temperatures of 200 and 175 • C, respectively, with homogeneous sulfur distribution through the ACF thickness after 2 h of infiltration. In terms of the porous carbon host, pores smaller than 0.4 nm could not accommodate any sulfur allotrope. For pores smaller than 0.7 nm, sulfur resided mainly in the form of S 4 and of S 6 . S 8 resided in micropores between 0.7 and 1.5 nm.

Conflicts of Interest:
The authors declare no conflict of interest.