Structure Optimization of Some Single-Ion Conducting Polymer Electrolytes with Increased Conductivity Used in “Beyond Lithium-Ion” Batteries

Simulation techniques implemented with the HFSS program were used for structure optimization from the point of view of increasing the conductivity of the batteries’ electrolytes. Our analysis was focused on reliable “beyond lithium-ion” batteries, using single-ion conducting polymer electrolytes, in a gel variant. Their conductivity can be increased by tuning and correlating the internal parameters of the structure. Materials in the battery system were modeled at the nanoscale with HFSS: electrodes–electrolyte–moving ions. Some new materials reported in the literature were studied, like poly(ethylene glycol) dimethacrylate-x-styrene sulfonate (PEGDMA-SS) or PU-TFMSI for the electrolyte; p-dopable polytriphenyl amine for cathodes in Na-ion batteries or sulfur cathodes in Mg-ion or Al-ion batteries. The coarse-grained molecular dynamics model combined with the atomistic model were both considered for structural simulation at the molecular level. Issues like interaction forces at the nanoscopic scale, charge carrier mobility, conductivity in the cell, and energy density of the electrodes were implied in the analysis. The results were compared to the reported experimental data, to confirm the method and for error analysis. For the real structures of gel polymer electrolytes, this method can indicate that their conductivity increases up to 15%, and even up to 26% in the resonant cases, via parameter correlation. The tuning and control of material properties becomes a problem of structure optimization, solved with non-invasive simulation methods, in agreement with the experiment.


Introduction
The "beyond lithium-ion" batteries implemented with single-ion conducting polymer electrolytes (SIPEs) present a higher energy density when the materials for electrolyte, electrodes, and moving ions are properly chosen.This represents a complex task, due to the fact that these materials are strongly correlated, and internal interactions occur at the molecular level [1][2][3].
The solid SIPEs have no more disadvantages like electrolyte leakage or flammability, which are characteristics of batteries with liquid electrolytes.For these conducting polymers, anions are fixed to the polymeric chain and only cations are capable of migrating in the structure, from one chain to another, but bind and unbind with the chains by weak molecular forces.Thus, only one type of ion, the cation, participates in the electrical conduction inside the battery cell.This behavior presents some advantages, namely that concentration polarization does not appear in the SIPEs with immobile anions, which implies the absence of cation depletion.For this type of SIPEs, the Sand and diffusion coherences are invalid, and cation conduction occurs by migration and not by diffusion, like in the dual-ion batteries [4][5][6].
Usually, electrodes chosen for different categories of electrolytes can be an insertion or a conversion material type.In the case of SIPE electrodes, the proper choice is represented by the conversion electrode type, usually carbonates or ethereal solvents.Other categories of materials can be an option, which involve bond breakage and significant structural changes in the electrode materials during redox reactions, like light chalcogens, sulfur, and oxygen, considering that the electrode material has to be compatible with the movable cations' nature.
In the electrode material, redox reactions occur, and the conversion anode stores and releases cations, storing charge via cation replacement in an ionic compound by the active ion.The negative charge state can be stabilized by conjugated structures via charge delocalization [7,8].The energy density generated by an electrode is directly correlated with the charge capacity and the redox potential of the active species involved.The compatibility between cathode material and movable cation species should be chosen using the periodical system, analyzing valences and redox potential [7].Redox-active polymers and also carbonyl-containing polymers, particularly polyquinones and polyimides, are also good candidates for electrodes, also in correlation with the polymer electrolyte.These structures present the advantage of higher energy density due to the complex phenomena occurring inside, which facilitates the conductivity increase.In these cases, the possibility of tuning the structural parameters that control the conductivity in the cell is also considerable.
The gel polymer electrolytes are a good option for performance battery cells because a solvent inserted in a polymer helps in ion transport across the polymer matrix.The solvent has to be chosen with a suitable dielectric constant, which is noted to impact ion conductivity.The gel polymers are typically prepared via thermal or UV curing methods by initially impregnating liquid precursors inside the electrode [1,[9][10][11].It is compulsory to apply an in situ polymerization method, which presents the advantage of decreasing the interfacial resistance between the solid electrolyte and electrodes [12,13].
In the new "beyond lithium-ion" batteries, the Li + ions are replaced with cations like Na + -which has similar physico-chemical properties to lithium-or with K + , Mg 2+ and Ca 2+ , or Al 3+ .The interest in these materials is for their accessibility, cost, correlation with electrode materials, easy conduction in the electrolyte, etc., all at an even lower cost.The valences, ionic radii, and mobilities of the alternative cations are of interest when the proper electrodes-electrolyte materials are chosen.Different variants are tested in the literature and new phenomenological valences occur when studies are evolving.
Consequently, the purpose of the new SIPE variants for electrolytes is to have an increased conductivity (the goal is a conductivity higher than 10 −3 S/m), which can be achieved by increasing their chain mobility and tuning the strength of the interaction between the polymer polar groups and the active ion [6,14], which is no more the classical lithium.This represents the main interest of our study, where the idea is to analyze and optimize different SIPE variants.
A non-invasive analysis method, which is very flexible and facilitates the variation in internal and external parameters, is of much interest in practice.The structure optimization for deserving a specific task (e.g., improved conductivity in the case of SIPEs) can be efficiently performed when more than one parameter is varied and this variation has to be correlated.Following this idea, the simulation method of the intimate structure at the molecular level, with fair consideration of the internal interaction between the constituents, can provide us with information about the composite material behavior when different stimuli are applied.
We performed our study for the case of "Beyond lithium-ion" batteries with improved cations conductivity, comprising a gel single-ion polymer electrolyte.The study target was "something else than Lithium cations", with possibilities of indicating results for new Polymers 2024, 16, 368 3 of 21 cation species, from an open list, and also "a larger variety of materials for the medium where conduction occurs", such as, for now, polymers, but the material list is also open.
The internal structure was described by simulation based on the coarse-grained molecular dynamics (CGMD) model, modified for the gel polymer electrolyte.The mesoscopic particle-based model or dissipative particle dynamics (DPD) model has been used for conductivity determination of the phases inside the polymeric gel.The electrolyte conductivity was calculated based on our method and also the electrodes' energy density.Structure optimization techniques via more-parameter correlation were revealed.The possible resonances corresponding to the conductivity maxima were identified when different structural parameters were correlated.Our new approach consists of physical parameters maximization and discovery of the resonances, by internal parameters correlation, and more than one internal parameter has to be modified to obtain a maximum (e.g., electrical conductivity maxima).
Our study was dedicated to improving the model for the new materials' analysis for batteries.

Materials and Methods
For analyzing the "beyond lithium-ion" batteries with SIPE, the structure of interest is represented by the electrolyte-electrodes-moving ions system.These elements are interdependent and the phenomena implied are correlated.
Different variants of structures have been analyzed to find possibilities to increase the conductivity inside.The gel polymer variant is recommended to be better, for which an amount of solvent is incorporated into the SIPE, compatible with the electrode material (carbonates, ethereal solvents).Conversion electrode materials (e.g., halides, light chalcogens) have also been considered and the results for the conductivity in the cell have been compared with those reported in the literature [15,16].
The single-ion versus dual-ion conducting electrolytes [5,17] were adopted for analysis.In the dual-ion batteries case, cations have to be inserted into the anode, and at the same time, anions are inserted into the cathode.These cations can be metal ions (Li + , Mg + , Ca 2+ , Al 3+ , etc.), while the anion is usually the counterion in the salt containing the cation, such as PF 6 − or TFSI − , etc.In dual-ion batteries, the internal processes are more complicated due to the two transport ion species inside, but they present the advantage of high energy densities, much easier to obtain in comparison with single-ion electrolytes.This represents no more an obstacle for implementing new variants, as single-ion batteries using polymers follow the tendency to be cheaper, safer, and easier to be synthesized.The cation migration in the single-ion structure is illustrated in Figure 1.
Table 1.The chemical structure of the molecular clusters (monomers and polymers) for the single-ion conducting polymer electrolytes considered for the analysis.
Table 1.The chemical structure of the molecular clusters (monomers and polymers) for the singleion conducting polymer electrolytes considered for the analysis.
The cathode represents the source of cations in the cell.The chosen conversion cathodes have been reported with high energy density for battery applications [26].The calculation of the density of charge by applying the conservation of charge law for the electrode system is of interest.The conversion anode stores and releases cations; in fact, it stores charge via cation replacement in an ionic compound via the active ion [27]: The cathode represents the source of cations in the cell.The chosen conversion cathodes have been reported with high energy density for battery applications [26].The calculation of the density of charge by applying the conservation of charge law for the electrode system is of interest.The conversion anode stores and releases cations; in fact, it stores charge via cation replacement in an ionic compound via the active ion [27]: where M is the cation, X is an anion, and A is the active ion implied in the reaction.Phenomenologically, a high ionic conductivity in SIPE can be achieved above the glass transition temperature in different conditions [28,29].For increasing the ionic conductivity in the cell, the polymeric chain mobility can be increased and the strength of the interaction between the polymer polar groups and the active ion can be tuned when the polymeric chain suffers modifications [30].The ion solubility has to be increased, but the ions must not become trapped [30].
Considering the internal phenomena, a mixture of ionic and electronic conductivity occurs in the cell, which can be characterized using a simulation model, developed for non-invasive testing.We focused on modeling the current transport and the potential field inside this complex structure of electrodes-electrolyte-conducting ions.We have to apply the charge conservation law for the three-dimensional electrical conduction inside the structure, illustrated by Equation ( 2): where i represents the electrical current, σ is the electrical conductivity, and ϕ is the electrical potential.The quantities i = i(x,y,z), that is, V = V(x,y,z), respectively, were considered in 3D representation, describing a volumetric conduction in the polymeric gel.The structure was described considering the intermolecular interactions of the micro-components, and at the macroscopic level, we had to take into account the ohmic losses in all the conducting materials (electrolyte, electrodes, and current collectors) and the contact resistance at the pairs of interfaces.Thus, the Laplace equation was solved.The effective conductivity of the polymeric gel was calculated based on the coarse-grained molecular dynamics (CGMD) model, modified for the gel polymer electrolyte, with details presented in the following.If for electrical potential, an evolution in an interval is considered (e.g., 0-5 V), then Equation (1) interconnects the current with the obtained potential, and the internal parameters of the structure and parametrical evolution of these quantities can be studied.Structural parameters like conducting ions radius or solvent concentration in SIPE were considered for analysis.Equation (1) was graphically solved using the Mathcad program, with the help of which the 3D gradient could be represented like a surface in space and intersected with the (0, 0) coordinates plane.The obtained values for i, when different internal parameters were changed, were represented on 3D graphs, which are available for structure optimization.
Each material was described at the molecular level in HFSS and then the proper mesh was set, in agreement with the molecular dimensions.The mesh was set for each particular analysis, due to the fact that the coarse-grained molecular dynamics (CGMD) model combined with the atomistic model was the basis of our simulation method.The model was developed by calculating the implied parameters' values in the case of a solid-state polymer electrolyte, and then the model was modified for the gel polymer electrolyte case.
The coarse-grained molecular scale was combined in analysis with the atomic scale, meaning an interval of 10 −10 to 10 −6 m and for each case a proper mesh.The polymer macromolecules are represented as beads.Each CG bead represents a group of monomers in a whole polymer structure.Internal interaction was considered between atoms in the chain, where the parameters of interest were the bond length, bond angle, and bond dihedral potentials, but a few degrees of freedom were neglected in the chain to simplify the global chain structure, less important for the interactions with other species.Thus, the beads for each cell configuration were selected.We had to consider as well the interactions with other molecular species at interfaces (interactions between beads and beads or between beads and other molecules (firstly, the solvent-polymer interactions).
In our case, the interacting species are electrolytes (polymer), electrodes, and moving ions.The parameters for calculation of these interacting forces are mainly the bead pair distance (r ij ), bond length (l ij ), bond angles (θ ijk ), and dihedral angles (ψ ijk ) (a 3D spatial configuration of the beads was considered).
Considering the polymer beads, the mesoscopic particle-based model or dissipative particle dynamics (DPD) model [31] was used for the conductivity determination of the phases inside the polymeric gel.The CG spheres interact with each other through purely repulsive soft potentials.These interactions between beads can be fine-tuned to capture the macroscopic phenomena on larger time scales.Polymeric coarse-grained beds interact with forces given by [31][32][33], illustrated in Equation (3): where i and j represent the indices of the two interacting beads, respectively, F C ij are conservative forces (non-bonded forces, electric, elastic, describing the repulsive properties between coarse-grained beads), F D ij are dissipative forces (friction, describing the friction dissipation between the structural system in the simulated bead), and F R ij are random forces (describes Brownian motion at ambient temperature) [29,30].
The most important term in expression ( 1) is determined by the conservative forces, which can be estimated with an expression as follows (Equation ( 4)), describing a soft pure repulsive (excluded volume) force through a distance r ij between beads i and j [6,33]: where a ij is the maximum repulsion force between beads I and j; ω C r ij is the weight function of the conservative force, ranging from 0 to 1, which represents a simple decaying function of the distance related with the cutoff distance r C , which is the interaction radius on the dimensions of the simulation system, correlated with the simulation mesh; e ij = r i −r j r ij is the unit vector pointing from particle j to particle i. a ij represents the dynamic interaction parameter that contains the physical-chemical information relevant to the atomic group and depends on the type of conservative force (repulsion force in our case) [6,33].Conservative forces were considered to vanish at distances greater than r C .Dissipative forces in expression (1) take into account the effects of viscosity, which slows down the particles' motion with respect to each other and is given in Equation ( 5) [6]: where γ ij is the friction coefficient; ω D r ij is the weight function of dissipative force; and v ij = v i − v j represents the particles' relative velocity.The random forces in expression (1) can be characterized considering the thermal or vibrational energy of the system [6], as it is indicated in Equation ( 6): where σ ij represents the magnitude of the noise; ω R r ij is the weight function of the stochastic force; ∆t is the integral time step and ζ ij is a random number between 0 and 1 which ensure the conservation of the total momentum.For each considered polymeric gel system, the bead structure was defined and the parameters necessary for interacting forces calculation were estimated (the bead pair distance (r ij ), bonds length (l ij ), bond angles (θ ijk ), and dihedral angles (ψ ijk )).
An important issue is the crosslink between the polymeric chain structural groups.Each type of crosslink can be described by the model, but the number and density of the crosslinks vary in the same structure depending on different parameters, like synthesis method, synthesis ambient conditions (e.g., temperature, pH, etc.), presence of a secondary adjuvant reactant, interaction forces with the gel, nature of the contact electrodes.During the operating process of the battery, the cations' dynamic at discharge and the voltage applied at charging can influence also the crosslinks panel.In this case, only a mediated crosslink density, decided after a great number of experiments, can offer us information about the matter.The presence of the crosslinks influences the cations' mobility inside the SIPEs, but is not the decisional factor.For deciding a result for the quantities describing the cations' possibilities of moving inside the crosslinked SIPEs, we considered a crosslink density justified by the values of the physical parameters of the polymer, meaning the molecular mass, the volumetric density, the polarity, and the elasticity coefficient.
A calculation of the ionic conductivity was performed using the interaction force inside the polymeric gel.The ions' mobility, µ(E, t), is decided by the interaction forces calculated above.The conductivity σ = e•p(x, y, z, t)•µ(E, t) where is e the elementary charge, p(x,y,z,t) is the charge carrier concentration, and E is the electric field.
The determined quantities were used as parameters inserted in the HFSS set-up, where the materials have been described.The results for conductivity of the medium components (polymeric gel) were obtained via simulation.
Generally, the electrolyte conductivity can be increased by increasing the chain mobility and tuning the strength of the interaction between the polymer polar groups and the active ion.These structural details can be implemented in the simulation set-up.
The effective conductivity formula that was used here was the general effective medium (GEM) equation [28,34], valuable for a medium composed of conducting and insulating phases inside the polymeric gel and is illustrated by Equation ( 7 where σ 1 , σ 2 represent the electric conductivities of the conducting and insulating phase, respectively, and σ ef is the effective conductivity of the polymeric gel; A is a constant that depends on the phases' nature; the exponent 1 ases inside the polymeri ases inside the polymeri depends on the filler volume fraction f and to the grain (beads) shape.The A and 1 ases inside the polymeri ases inside the polymeri parameters were calculated based on the results indicated in the literature for particular values of the effective conductivity determined in specific conditions for different polymeric gels and values of similar magnitude order have been used in the formula for new polymer variants until a confirmed result for the effective conductivity was obtained.An iterative process of approximation was applied until the results have been validated.The formula has to be applied for each conduction mechanism in the polymeric gel.
Then the current density j = σ ef •E was computed.The purpose was also to obtain the energy density of the gel polymer battery [35,36].This quantity can be estimated with the Formulas ( 8): where C cell is the capacity of the electrode; V cell is the operating cell voltage; m components represents the mass of individual components of the electrode; C theor is the theoretical capacity of anode, N/P is the excess anode content in the system, and N/P area is a constant that represents the fraction of excess area of the anode [35,36].

Results for the Material Parameters
The applied working strategy was: -Establishing the bead groups in the polymer structure; -Calculating all the interaction forces based on the theoretical model and using data from reference values like those in the literature; -Establishing the conduction mechanism; -Determining the charge carriers' mobility using the simulation set-up when the electrical stimuli were set; -Calculating the conductivity in the polymeric gel; -Calculating the current density; -Calculating the energy density corresponding to different electrode-electrolyte combinations; -Three-dimensional graphical representation in function of different parameters: conducting ions radius, solvent concentration in SIPE, energy density generated by the conversion electrode (correlated with the charge capacity and redox potential of the active species involved), strain and stress in the polymeric electrolyte; parameters corresponding to the beads: the bead pair distance, bonds length, bond angles, and dihedral angles; choosing of the simulation mesh; and external parameters: electric field for polymeric gel testing, temperature.
The current density throughout the domain was determined when different internal and external parameters of the exploited battery were varied, like the conducting ions radius, solvent concentration in SIPE, respectively the obtained electric potential.The considered electrode thickness was of 400 micrometers (µm), with particles inside on the order of µm, with deferential particle sizes (~7×), higher dimensions for active material particles, and lower for the binder.The ionic radii for the considered cations were: Na + , Mg 2+ , and Al 3+ : 116 pm (1.16 Å), 86 pm (0.86 Å), respectively 67.5 pm (0.675 Å).
The method was initially tested on the usually used conductive polymers (CPs) based on different polyethylene oxide, polycarbonates and polysiloxanes.The obtained results were in good agreement with those presented in the literature [11,14,15,18].
The current density versus battery potential were represented on double-parametrical graphs, in order to observe and optimize their dependence on the internal parameters of the structure, characterizing the mobile cations (Na + , K + , Mg 2+ and Ca 2+ , Al 3+ ) with different valences, the SIPE chain, the solvent incorporated into the SIPE or the conversion electrode.A set of results for the three SIPE mentioned above are illustrated in Figures 2-5.Simulations were performed with HFSS program (by Ansys: Ansys HFSS 2022 R1).
One observes that for higher current densities, ionic radius is recommended to be lower, the current density decreasing for bigger cations and also the battery voltage.
One observes that the solvent concentration has an optimum value for which the conduction in the polymer is maximized.Simulations have indicated that this maximum depends in the same time on solvent and movable cation nature.
One observes that a couple of maxima occur on the current density graphs, illustrating internal resonances between a couple of geometrical parameters (conducting ions radius with bead pair distances in our case), which can be identified with the help of the simulation method.
For the solid SIPEs that can be synthesized with increased conductivity, different methods of analysis were reported in the literature.A comparative analysis is imposed in order to reveal the advantages of our method of study.We have considered for illustration the three polymers species analyzed in the 3D parametrical graphs: PEGDMA-SS, PU-TFMSI, and PDMS-poly(MPA-Na + -r-PEGMA) 20 -PDMS 10 .The comparative analysis is presented in Table 2.We have presented in the table the current density evolution in the gelled SIPEs when different control parameters were varied in specific intervals, imposed by the SIPE nature.In these intervals, the current density of the SIPEs changes, but this property is tuned by the control parameters.These control parameters were enumerated in the table, specific for each analyzing method.One observes that our method offers the possibility of changing the greater number of parameters and their influence to be monitored simultaneously.The current density versus battery potential were represented on double-parametrical graphs, in order to observe and optimize their dependence on the internal parameters of the structure, characterizing the mobile cations (Na + , K + , Mg 2+ and Ca 2+ , Al 3+ ) with different valences, the SIPE chain, the solvent incorporated into the SIPE or the conversion electrode.A set of results for the three SIPE mentioned above are illustrated in Figures 2-5.Simulations were performed with HFSS program (by Ansys: Ansys HFSS 2022 R1).
One observes that for higher current densities, ionic radius is recommended to be lower, the current density decreasing for bigger cations and also the battery voltage.
One observes that the solvent concentration has an optimum value for which the conduction in the polymer is maximized.Simulations have indicated that this maximum depends in the same time on solvent and movable cation nature.
One observes that a couple of maxima occur on the current density graphs, illustrating internal resonances between a couple of geometrical parameters (conducting ions

Conclusions
The new SIPE variants present increased conductivity that can be tuned and improved by different methods, like increasing their chain mobility and tuning the strength of the interaction between the polymer polar groups and the active ion, which is a cation that migrates in the structure.These were the basic ideas in our effort to analyze and synthesize different SIPE variants with improved conductivity.
New results were obtained for different single-ion CPs like poly(ethylene glycol) dimethacrylate-x-styrene sulfonate (PEGDMA-SS), PU-TFMSI, or PDMS-poly(MPA-Na +r-PEGMA) 20 -PDMS 10 , but not only these.The current density versus battery potential was represented on double-parametrical graphs in order to observe and optimize their dependence on the internal parameters of the structure, characterizing the mobile cations (Na + , K + , Mg 2+ , and Ca 2+ , Al 3+ ), the SIPE chain, the solvent incorporated into the SIPE, and the conversion electrode.A set of results for the three SIPEs mentioned above was illustrated in this paper.The simulations were performed with the HFSS program (by Ansys: Ansys HFSS 2022 R1) based on the modified (CGMD) model for the structural description of the cell constituents.
The conductivity of the gel polymeric electrolyte depends on different parameters: conducting ions radius, solvent concentration in SIPE, energy density generated by the conversion electrode (correlated with the charge capacity and redox potential of the active species involved), strain and stress in the polymeric electrolyte; parameters corresponding to the beads: the bead pair distance, bonds length, bond angles, and dihedral angles; choosing of the simulation mesh; and external parameters: electric field for polymeric gel testing, temperature.
If we consider the matter of an optimal choice of the materials for batteries (SIPEs, gel, and cations) among the analyzed structures, the most important conclusion is that better results for the SIPE conductivity can be obtained only by parameter correlation, e.g., due to the fact that the current density decreases for bigger cations and the battery voltage is important, the Al 3+ cation with a smaller radius can be a choice.Also, the increased valence number (+3) of a cation offers more possibilities to interact with the bounded anions for increased mobility.Another example: when the solvent is chosen, the solvent concentration has an optimum value for which the conduction in the polymer is maximized, depending on solvent and movable cation nature.The gelled PU-TFMSI (+(BzO) 2 or THF) can be a choice in this case.
The influence of the cathode material has to be also considered, considering the energy density generated by the conversion cathode.The results illustrated for the Na-ion batteries recommend the poly(anthraquinonyl imide) cathode, with better results, rather than an p-dopable polytriphenyl amine cathode associated with the same cation type.
Another important aspect revealed by our study is the occurrence of the SIPE's current density maxima, illustrated on some double-parametrical graphs.These maxima are of great interest in practice and their occurrence is an event that has to be caused inside the battery physical system.These maxima can be obtained in results of the internal resonances determined by correlation of specific geometrical parameters (conducting ions radius with bead pair distances in our case), which can be identified with the help of the simulation method.Mg-cation batteries and also Al-cation batteries present these maxima when we adjust the energy density generated by the conversion cathode.Our study performed for more battery constituents revealed the presence of these maxima for every type of cation, but their occurrence is difficult to obtain due to the fact that it implies a very fine correlation of the structural parameters (the phenomenon is a resonant one).From this point of view, the role of the simulation is essential.
The idea of current density maxima opens perspectives and a field for future work in order to identify the possible resonances occurring in the material sample.If these resonances are caused in practice by tuning properly and correlating internal parameters, a consistent increase in conduction can be caused at the structure level.

Figure 1 .
Figure 1.Single-ion conductive polymer working principle in the battery cell.(a) Structure of electrodes-electrolyte; (b) cation migration in the polymer matrix, from one chain to another; fixed anions can be observed.

Figure 1 .
Figure 1.Single-ion conductive polymer working principle in the battery cell.(a) Structure of electrodes-electrolyte; (b) cation migration in the polymer matrix, from one chain to another; fixed anions can be observed.

Figure 2 . 2 2Figure 2 .Figure 3 .Figure 3 .Figure 4 .
Figure 2.(a,b) Current density in function of the electric potential, for three species of SIPEs (1), (2), and (3) used for "beyond lithium-ion" batteries electrolyte, considering the ionic radius of the mobile cation like parameter; plots are represented in two positions for clarity; (c) Cross-section of the surface plots for constant ionic radius of 0.675 Å (continuous curves), respectively cross-section of the surface plots for constant ionic radius of 2.5 Å (dotted curves), represented on the same 2D graph; (d) Cross-section of the surface plots for constant electric potential of 1 V (continuous curves), respectively cross-section of the surface plots for constant electric potential of 2.3 V (dotted curves), represented on the same 2D graph.

Figure 4 . 21 Figure 5 .
Figure 4. (a,b) Current density in function of the electric potential for Na-ion batteries, considering the energy density generated by the conversion cathode like parameter.Surface graph (cat1) corresponds to the p-dopable polytriphenyl amine cathode; graph (cat2) is for poly(anthraquinonyl imide) cathode; plots are represented in two positions for clarity; (c) Cross-section of the surface plots for constant energy density of 40 Wh kg −1 (continuous curves), respectively cross-section of the surface plots for constant energy density of 180 Wh kg −1 (dotted curves), represented on the same 2D graph; (d) Crosssection of the surface plots for constant electric potential of 1 V (continuous curves), respectively cross-section of the surface plots for constant electric potential of 2.3 V (dotted curves), represented on the same 2D graph.

Figure 5 .
Figure 5. (a,b) Current density in function of the electric potential for Mg-ion batteries (Mg batt) versus Al-ion batteries (Al batt), considering the energy density generated by the conversion cathode like parameter.The sulfur cathode Mo 6 S 8 was considered for Mg-ion batteries (orange graph), respectively the Co 3 S 4 cathode for Al-ion batteries (green graph); plots are represented in two positions for clarity; (c) Cross-section of the surface plots for constant energy density of 80 W•h•kg −1 (continuous curves), respectively cross-section of the surface plots for constant energy density of 320 W•h•kg −1 (dotted curves), represented on the same 2D graph.Curves for the first maxima, occurring at ca. 220 W•h•kg −1 for Al batteries, respectively at 182 W•h•kg −1 for Mg batteries, were also represented on the graph.(d) Cross-section of the surface plots for constant electric potential of 1 V (continuous curves), respectively cross-section of the surface plots for constant electric potential of 2.3 V (dotted curves), represented on the same 2D graph.

Table 1 .
The chemical structure of the molecular clusters (monomers and polymers) for the singleion conducting polymer electrolytes considered for the analysis.

Table 2 .
Method of analysis applied for the current density determination in the "beyond lithium-ion" batteries' structures with the considered SIPEs and comparison of the results and methods indicated in the literature with our method.