Sigma-Holes in Battery Materials Using Iso-Electrostatic Potential Surfaces

The presence of highly electronegative atoms in Li-ion batteries anticipates the formation of σ-hole regions that may strongly affect the ionic conductivity. The σ-hole consists of a region of positive electrostatic potential extending in the direction of the covalent bond between atoms of groups IV–VII due to anisotropic charge distribution. Graphite electrodes in Li-ion batteries that become halogenated due to the electrolyte, as well as some solid electrolyte materials, can exhibit these σ-holes. Since Li-ions should be able to drift in any part of the battery, the fact that they can be attracted and eventually absorbed by regions of strong negative potentials produced by high-electronegativity counterions becomes detrimental to ionic conductivity. Therefore, the presence of positive well-defined regions, repulsive to the Li-ions, might act as lubricant for Li-ions drifting through electrolytes, thus improving the Li-ion conductivity. In addition, the σ-holes might also have a strong effect on the formation of the passivating layer, known as the solid electrolyte interphase (SEI) at electrode surfaces, which is of paramount importance for the performance of rechargeable batteries. Here we investigate the existence of σ-holes on surfaces of graphite anodes and of a few solid electrolytes by examining the electrostatic potentials calculated using density functional theory.


Introduction
σ-holes are particular electronic arrangements where a region of positive electrostatic potential extends in the direction of the covalent bond between atoms of groups IV-VII due to anisotropic charge distribution.Graphite electrodes in Li-ion batteries that become halogenated due to electrolyte decomposition, as well as of some solid electrolyte materials, can exhibit these σ-holes.The importance of these electronic distributions is that they can favor ionic or electronic motion in directions that would not be possible in their absence.Here, we study them in relation to battery materials, where ionic and electronic transport is of utmost importance.
In particular, the formation of a solid electrolyte interphase (SEI) between anode and electrolyte in Li-ion batteries takes place due to the leakage of electrons from the negative anode (generally graphite) to the electrolyte (solutions based on ethylene-carbonate) in a complex set of reactions involving the counterions and any other additives inserted in the electrolyte to improve their chemical behavior.The leaked electron reacts with the components of the electrolyte, forming a variety of moieties that settle down on the surface of the anode.The SEI thus grows as more electrons are leaked; however, the rate of leakage decreases as the thickness of the SEI increases.Usually, a commercial battery reaches such passivation in the very first cycle of its use or, at most, in the initial cycles (a few hours).This explains why drift velocities of Li-ions in rechargeable batteries are extremely slow, ~1 mm/h.Electrolytes of present batteries are usually in liquid phase.This is usually the state of matter that allows the best conductivity of positive ions while maintaining good contact to the electrodes.Thus, a typical Li-ion rechargeable battery uses LiCoO 2 as the positive electrode (cathode) [1], graphite as the negative electrode (anode), and a liquid electrolyte solution of LiPF 6 dissolved in ethylene carbonate (EC) [2,3].However, present requirements of higher energy and power for transportation and storage backup for electrical grids demand a radical change of these typical components of a battery.To satisfy energy density requirements, proposed materials include sulfur for cathodes [4] along with lithium metal for anodes [5].Both materials bring new challenges; one of them is the formation of dendrites [6], which may cause internal short-circuits in the battery, which is a severe safety issue.One possible solution to avoid the formation of dendrites is the use of solid electrolytes simply because they offer better resistance to the growth of dendrites than liquids [7].Interfacial behavior of the solid electrolytes with anode materials, such as silicon, have been studied [8][9][10][11][12].A downside of solid electrolytes is their ionic conductivity, which is on the order of 10 −3 S cm −1 [13], lower than the ones obtained for organic liquid electrolytes, which are on the order of 10 −2 S cm −1 [14].In spite of the vast research in this field, practical solid electrolytes still do not meet the desired level of ionic conductivities and have the problem of poor interfacial contact (solid-to-solid) to the electrodes.However, a few solid electrolyte materials have shown promising high ionic conductivity, becoming possible options for higher-capacity batteries, including, for example: Li 3 PS 4 [15], Li 6 PS 5 Cl [16], and Li 7 P 2 S 8 I [17].Li 7 P 2 S 8 I (LPSI), was intensively studied as a solid electrolyte because it exhibits a high Li ionic conductivity of 0.64 mS cm −1 , and it has high electrochemical stability, being able to stand voltages up to 10 V (vs.Li/Li + ) at room temperature [17,18].Therefore, understanding and being able to predict their ionic conductivity mechanisms as a function of their electronic structure is necessary.
The σ-hole designates the electron-deficient outer lobe of a half-filled p orbital involved in forming a covalent bond between atoms of groups IV-VII.It involves a region of positive electrostatic potential which can interact attractively (non-covalently) with negative sites on other molecules.In this work, we report the presence of σ-holes or, in general, exposed regions of positive potential in places that were previously thought to exhibit negative charge, including, for example, the outer portions of halogen atoms opposite a bond with another nonmetal element [19][20][21][22][23].The susceptibility of this positive charge to nonbonding interactions with its neighbors [24][25][26], as well as the possibility of its repulsion of other positively charged species, make these examinations necessary for further innovations in materials.

Methodology
The term ∑ A Z A |R A −r| is known in DFT as the external potential, v(r), which uniquely determines the density ρ(r), or vice versa, in Coulomb's Law, Here, Z A represents charge on nucleus A at position R A .Each term contributes independently for the nuclei and electrons, respectively [27].Thus, as concluded earlier [28], the electrostatic potential function is a transformation of the electron density and therefore of the wavefunction, containing all information available in both, and with the additional feature of presenting such information in a more suitable way to study interactions between molecules (intermolecular), as well as interactions in molecules (intramolecular).Poisson's equation, beyond previously-explored classical methods, relates electrostatic potential with surrounding electron density through the equation , where δ represents the three-dimensional Dirac delta function.Important interactions include the phenomena at the boundary between positive and negative regions of a molecule, cluster, or crystal, which are of utmost interest for all kinds of physical and chemical interactions with other systems.Interactions that are almost otherwise impossible to track with electron density and wavefunction properties can be revealed with these DFT models [29].Molecular electrostatic potentials (MEPs), in recent applications, Crystals 2018, 8, 33 3 of 13 have become powerful means of analyzing sensitive binding mechanisms, such as the pairing of DNA bases [28,30] and the handling of volatile materials that require precise detection of minute traces [31].MEP's ability to serve as a logic gate has also found use in information processing as a basis for data operations [32][33][34][35][36].
All calculations of extended systems (crystals and surfaces) were performed using the program Vienna ab initio simulation package (VASP) [37][38][39] and those of finite systems (molecules) using the Gaussian-09 program [40].For the extended systems, DFT with periodic boundary conditions and the projector augmented wave (PAW) [41] method were used.In this method, the generalized pseudopotential [42] and the linear augmented-plane-wave (LAPW) [43] are combined, allowing high-quality first-principle calculations.The PBE [44] functional (for extended systems) and the Becke three-parameter hybrid exchange functional [45] with the first-principles generalized gradient approximation correlation functional Perdew-Wang-91 [46], B3PW91 (for molecular systems).An augmented plane-wave basis set with maximum energy of 300-400 eV was used for extended systems and the double-ξ plus polarization 6-31G(d) basis set was used for molecules [47,48].
The material used to study the crystallographic surfaces is modeled as a slab of the crystal structure (graphite, Li 7 P 2 S 8 I [17], and Li 3 PS 4 [15]) with a vacuum above periodicity parallel to the slab on the order of 18 Å, large enough to isolate the surface from its periodic images.The materials chosen for this study are: graphite, the most used material in commercial Li-ion batteries, and the solid electrolytes Li 7 P 2 S 8 I and Li 3 PS 4 , both of which are of interest given their systems containing σ-holes; the electrostatic potentials of the solid electrolytes are calculated at (010) and (001) surfaces, respectively (Figure 1).

Results and Discussion
The iso-MEP surfaces of the benzene molecule are illustrated in Figure 2. The positive region (blue) expands on the neighborhood of all atoms due to the strong effect of the nuclei.By contrast, the negative regions above and below the plane of the molecule are due to the pi-electrons, leading to a net negative surface.The two cases shown in Figures 2 and 3 illustrate very clearly the effect of the isopotential surfaces.The typical way to display a MEP is to use isodensity surfaces mapped with MEP values.These colorful pictures provide a good qualitative idea of the MEP distribution around a molecule.From the figures, we can see the new approach can be used to show the domains at which the σ-hole is active in halogenated compounds, which the MEP spectrum tends to miss by illustrating the potential on isodensities.A second case is given by chlorobenzene C6H5Cl, shown in Figure 3.For the halogen substituent chlorine, a negative surface covers the positive at low voltages close to 0, but at ±0.3 V, the presence of a σ-hole exposes the positive region underneath.We can attribute this effect to anisotropic distribution of charge leading to accumulation along the equatorial regions of chlorine when it participates in a covalent bond.Therefore, the outermost portion of the halogen atom along the C-Cl bond axis is left vulnerable to possible nonbonding interactions if a negative charge approaches that site (Figure 3).At potentials close to 0, such as 30 mV, the hole is not yet present, while on existing isodensity surfaces of PhCl and halogenated benzenes, it would appear exposed [21].Notable differences for the larger halogens include larger σ-hole circumference at identical isopotential values [24]; the size of the exposed region for phenyl halogens follows the trend PhI > PhBr > PhCl at the ±0.3 V isopotentials.This approach confirms existing literature about increasing halogen capabilities of participating in these types of nonbonding interactions due to differences in size [49,50].Next, we look to take this method to PF  , the anion of lithium hexafluorophosphate, LiPF6, one of the most common salts for Li-ion batteries.Its neutral form has not been reported, but may exist for a very short period of time if removal of an electron from the PF  anion is achieved, upon which the molecule decomposes into PF5 with the release of one of its fluorine substituents [51].Its negative ion exists in solution (Figure 4).The neutral structure of PF6 has lower symmetry than its anion, with each molecule exhibiting D4h and Oh point group symmetries, respectively.Two negative caps of A second case is given by chlorobenzene C 6 H 5 Cl, shown in Figure 3.For the halogen substituent chlorine, a negative surface covers the positive at low voltages close to 0, but at ±0.3 V, the presence of a σ-hole exposes the positive region underneath.We can attribute this effect to anisotropic distribution of charge leading to accumulation along the equatorial regions of chlorine when it participates in a covalent bond.Therefore, the outermost portion of the halogen atom along the C-Cl bond axis is left vulnerable to possible nonbonding interactions if a negative charge approaches that site (Figure 3).At potentials close to 0, such as 30 mV, the hole is not yet present, while on existing isodensity surfaces of PhCl and halogenated benzenes, it would appear exposed [21].Notable differences for the larger halogens include larger σ-hole circumference at identical isopotential values [24]; the size of the exposed region for phenyl halogens follows the trend PhI > PhBr > PhCl at the ±0.3 V isopotentials.This approach confirms existing literature about increasing halogen capabilities of participating in these types of nonbonding interactions due to differences in size [49,50].A second case is given by chlorobenzene C6H5Cl, shown in Figure 3.For the halogen substituent chlorine, a negative surface covers the positive at low voltages close to 0, but at ±0.3 V, the presence of a σ-hole exposes the positive region underneath.We can attribute this effect to anisotropic distribution of charge leading to accumulation along the equatorial regions of chlorine when it participates in a covalent bond.Therefore, the outermost portion of the halogen atom along the C-Cl bond axis is left vulnerable to possible nonbonding interactions if a negative charge approaches that site (Figure 3).At potentials close to 0, such as 30 mV, the hole is not yet present, while on existing isodensity surfaces of PhCl and halogenated benzenes, it would appear exposed [21].Notable differences for the larger halogens include larger σ-hole circumference at identical isopotential values [24]; the size of the exposed region for phenyl halogens follows the trend PhI > PhBr > PhCl at the ±0.3 V isopotentials.This approach confirms existing literature about increasing halogen capabilities of participating in these types of nonbonding interactions due to differences in size [49,50].Next, we look to take this method to PF  , the anion of lithium hexafluorophosphate, LiPF6, one of the most common salts for Li-ion batteries.Its neutral form has not been reported, but may exist for a very short period of time if removal of an electron from the PF  anion is achieved, upon which the molecule decomposes into PF5 with the release of one of its fluorine substituents [51].Its negative ion exists in solution (Figure 4).The neutral structure of PF6 has lower symmetry than its anion, with each molecule exhibiting D4h and Oh point group symmetries, respectively.Two negative caps of Next, we look to take this method to PF 6 − , the anion of lithium hexafluorophosphate, LiPF 6 , one of the most common salts for Li-ion batteries.Its neutral form has not been reported, but may exist for a very short period of time if removal of an electron from the PF 6 − anion is achieved, upon which the molecule decomposes into PF 5 with the release of one of its fluorine substituents [51].Its negative ion exists in solution (Figure 4).The neutral structure of PF 6 has lower symmetry than its anion, We can infer, at least in part, the existence of σ-holes in these molecules, even if only exhibited over a small range.For neutral PF6, moving over a short range from ±0.30 to ±0.35 V gives us a fleeting instance of a σ-hole, which disappears quickly as we reach nuclear domain potentials close to the atoms (Figure 5).The negative ion, PF  , gives us a better look at the σ-hole when moving in isopotential ranges beyond ±5 V (Figure 6).Holes can be seen on the fluorine caps as predicted by the isodensity surface We can infer, at least in part, the existence of σ-holes in these molecules, even if only exhibited over a small range.For neutral PF 6 , moving over a short range from ±0.30 to ±0.35 V gives us a fleeting instance of a σ-hole, which disappears quickly as we reach nuclear domain potentials close to the atoms (Figure 5).We can infer, at least in part, the existence of σ-holes in these molecules, even if only exhibited over a small range.For neutral PF6, moving over a short range from ±0.30 to ±0.35 V gives us a fleeting instance of a σ-hole, which disappears quickly as we reach nuclear domain potentials close to the atoms (Figure 5).The negative ion, PF  , gives us a better look at the σ-hole when moving in isopotential ranges beyond ±5 V (Figure 6).Holes can be seen on the fluorine caps as predicted by the isodensity surface The negative ion, PF 6 − , gives us a better look at the σ-hole when moving in isopotential ranges beyond ±5 V (Figure 6).Holes can be seen on the fluorine caps as predicted by the isodensity surface in Figure 6, and later the phosphorus center becomes visible underneath the negative region, with the equatorial regions of fluorine being the last portions to erode around the ±5.7 V isopotential value.
in Figure 6, and later the phosphorus center becomes visible underneath the negative region, with the equatorial regions of fluorine being the last portions to erode around the ±5.7 V isopotential value.With this in mind, we expand the concept of MEP to extended systems such as crystal electrostatic potentials (CEP) and their surfaces.The first case that calls our attention is the (001) surface of graphite, one of the multiple possible orientations of graphite in batteries.Although in disordered carbons, other graphite facets may be exposed, the one chosen is the most favorable to Li ion intercalation.Even in passivated surfaces such as those H-terminated, the attachment of halogen atoms from the electrolyte is highly possible.To some extent, the σ-hole was observed when fluorine was tested (Figure 7).Fluorine comes directly from the decomposition of PF6 − on the surface or about the electrode, producing a σ-hole from −2.80 to −2.75 V; for other potentials there is no sigma-hole formation.With this in mind, we expand the concept of MEP to extended systems such as crystal electrostatic potentials (CEP) and their surfaces.The first case that calls our attention is the (001) surface of graphite, one of the multiple possible orientations of graphite in batteries.Although in disordered carbons, other graphite facets may be exposed, the one chosen is the most favorable to Li ion intercalation.Even in passivated surfaces such as those H-terminated, the attachment of halogen atoms from the electrolyte is highly possible.To some extent, the σ-hole was observed when fluorine was tested (Figure 7).Fluorine comes directly from the decomposition of PF 6 − on the surface or about the electrode, producing a σ-hole from −2.80 to −2.75 V; for other potentials there is no sigma-hole formation.With this in mind, we expand the concept of MEP to extended systems such as crystal electrostatic potentials (CEP) and their surfaces.The first case that calls our attention is the (001) surface of graphite, one of the multiple possible orientations of graphite in batteries.Although in disordered carbons, other graphite facets may be exposed, the one chosen is the most favorable to Li ion intercalation.Even in passivated surfaces such as those H-terminated, the attachment of halogen atoms from the electrolyte is highly possible.To some extent, the σ-hole was observed when fluorine was tested (Figure 7).Fluorine comes directly from the decomposition of PF6 − on the surface or about the electrode, producing a σ-hole from −2.80 to −2.75 V; for other potentials there is no sigma-hole formation.The electrolyte will react with the anode surface during the first charge [52,53], producing a multicomponent SEI layer.The SEI components may have halogen atoms, such as Cl and I, depending on the electrolyte.A σ-hole was observed when chlorine was on the graphite surface (Figure 8).Chlorine comes from the decomposition of Li 6 PS 5 Cl during the reaction of the first charge.The chlorine produces a small σ-hole from −2.55 to −2.3 V, and for other potentials there is no sigma-hole formation.The electrolyte will react with the anode surface during the first charge [52,53], producing a multicomponent SEI layer.The SEI components may have halogen atoms, such as Cl and I, depending on the electrolyte.A σ-hole was observed when chlorine was on the graphite surface (Figure 8).Chlorine comes from the decomposition of Li6PS5Cl during the reaction of the first charge.The chlorine produces a small σ-hole from −2.55 to −2.3 V, and for other potentials there is no sigmahole formation.In the case of the iodine atom, the formation of the σ-hole was observed (Figure 9).We found that the σ-hole is larger than that of Cl and F, corroborating earlier visualizations of isolated molecules.The iodine atoms might come from the decomposition of small regions of the Li7P2S8I solid electrolyte; however, at this point this is just an assumption.The I atom produces a large σ-hole, while Cl produces a small one, and F an even smaller one.This is directly related to their position in the halogen column on the periodic table.The smallest and most electronegative atom makes a smaller σ-hole.The I atom produces a σ-hole that forms from −2.55 to −2.25 V, and for other potentials there is no σ-hole formation.In the case of the iodine atom, the formation of the σ-hole was observed (Figure 9).We found that the σ-hole is larger than that of Cl and F, corroborating earlier visualizations of isolated molecules.The iodine atoms might come from the decomposition of small regions of the Li 7 P 2 S 8 I solid electrolyte; however, at this point this is just an assumption.The I atom produces a large σ-hole, while Cl produces a small one, and F an even smaller one.This is directly related to their position in the halogen column on the periodic table.The smallest and most electronegative atom makes a smaller σ-hole.The I atom produces a σ-hole that forms from −2.55 to −2.25 V, and for other potentials there is no σ-hole formation.The electrolyte will react with the anode surface during the first charge [52,53], producing a multicomponent SEI layer.The SEI components may have halogen atoms, such as Cl and I, depending on the electrolyte.A σ-hole was observed when chlorine was on the graphite surface (Figure 8).Chlorine comes from the decomposition of Li6PS5Cl during the reaction of the first charge.The chlorine produces a small σ-hole from −2.55 to −2.3 V, and for other potentials there is no sigmahole formation.In the case of the iodine atom, the formation of the σ-hole was observed (Figure 9).We found that the σ-hole is larger than that of Cl and F, corroborating earlier visualizations of isolated molecules.The iodine atoms might come from the decomposition of small regions of the Li7P2S8I solid electrolyte; however, at this point this is just an assumption.The I atom produces a large σ-hole, while Cl produces a small one, and F an even smaller one.This is directly related to their position in the halogen column on the periodic table.The smallest and most electronegative atom makes a smaller σ-hole.The I atom produces a σ-hole that forms from −2.55 to −2.25 V, and for other potentials there is no σ-hole formation.The presence of σ-holes and their nonbonded interactions in sulfur-containing compounds [54,55] leads us to extend this model to possible solid electrolyte materials that can see use in Li-ion batteries.Li3PS4, whose PS4 tetrahedron can feature σ-holes on sulfur, was studied as a crystal on the (001) surface.The formation of σ-holes has a preference to PS4 (Figure 10), the isopotential can be found from −3.3 to −2.8 V, and for other potentials there is no σ-hole formation.The other solid electrolyte material studied is Li7P2S8I on the (010) surface, whose PS4 tetrahedron, like Li3PS4 before it, suggests σ-hole formation.We found σ-holes in the surface along the b-axis, as shown in Figure 11a-c obtained at +0.5 and −3.8 V.The crystal structure of Li7P2S8I, as The presence of σ-holes and their nonbonded interactions in sulfur-containing compounds [54,55] leads us to extend this model to possible solid electrolyte materials that can see use in Li-ion batteries.Li 3 PS 4 , whose PS 4 tetrahedron can feature σ-holes on sulfur, was studied as a crystal on the (001) surface.The formation of σ-holes has a preference to PS 4 (Figure 10), the isopotential can be found from −3.3 to −2.8 V, and for other potentials there is no σ-hole formation.The presence of σ-holes and their nonbonded interactions in sulfur-containing compounds [54,55] leads us to extend this model to possible solid electrolyte materials that can see use in Li-ion batteries.Li3PS4, whose PS4 tetrahedron can feature σ-holes on sulfur, was studied as a crystal on the (001) surface.The formation of σ-holes has a preference to PS4 (Figure 10), the isopotential can be found from −3.3 to −2.8 V, and for other potentials there is no σ-hole formation.The other solid electrolyte material studied is Li7P2S8I on the (010) surface, whose PS4 tetrahedron, like Li3PS4 before it, suggests σ-hole formation.We found σ-holes in the surface along the b-axis, as shown in Figure 11a-c obtained at +0.5 and −3.8 V.The crystal structure of Li7P2S8I, as The other solid electrolyte material studied is Li 7 P 2 S 8 I on the (010) surface, whose PS 4 tetrahedron, like Li 3 PS 4 before it, suggests σ-hole formation.We found σ-holes in the surface along the b-axis, as shown in Figure 11a-c obtained at +0.5 and −3.8 V.The crystal structure of Li 7 P 2 S 8 I, as shown in Figure 11d, reveals the insertion of LiI in the interstitial sites of β-LPS.Additional Li were placed in the Li sites, and the iodine atoms were placed in lattice points of the Pnma space group.Thus, the crystal structure of Pnma is transformed into a similar P4 2 /nmc with a unique b-axis [18].We also determined different isopotentials by keeping the positive potential constant at +0.5 V and increasing the negative from −3.8 to −3.4 V as shown in Figure 12.The σ-hole closes completely by −3.4 V.

Conclusions
Successful demonstration of the σ-hole in molecules allows a larger-scale study by taking advantage of periodic boundary conditions of first-principle programs, such as VASP, when examining isopotentials in crystal surfaces.Previous assumptions of electrostatic behavior in halogen substituents, particularly those of covalently bonded molecules comprised of atoms in groups IV-VII, were previously thought to be uniformly negative, while further examination of both finite and extended systems shows that this might not be the case.By extending crystalline systems and examining isopotential surfaces of periodically repeating molecular structures, we are able to draw more realistic inferences of the molecular environments beyond the vacuum surrounding a single molecule.The systems investigated in this work resemble the molecular environments of graphite anodes and solid electrolytes.Even though the employed models are simplifications of real materials, the electronic structures and their main trends are well-captured, as demonstrated by the vast, well-known, first principles modeling battery materials literature.We found that, for these cases, the σ-holes in the solid electrolytes form preferentially at the PS 4 tetrahedron.For solid electrolytes, we confirm the important role played by the PS 4 in solid electrolytes for the diffusion of lithium ions.In particular, positive regions that would repel the Li-ions might act as driving-forces for Li-ions drifting through electrolytes, thus improving the Li-ion conductivity.It may also have important consequences for the electronic transport, as recently reported [56,57].

Figure 2 .
Figure 2. Isopotential surfaces of benzene, separated between negative (red) and positive (blue) regions corresponding to the same absolute value of potential.Isopotential pairs are ±0.03(second row) and ±0.3 V (bottom).A molecule is depicted for position reference.

Figure 2 .
Figure 2. Isopotential surfaces of benzene, separated between negative (red) and positive (blue) regions corresponding to the same absolute value of potential.Isopotential pairs are ±0.03(second row) and ±0.3 V (bottom).A molecule is depicted for position reference.

Figure 2 .
Figure 2. Isopotential surfaces of benzene, separated between negative (red) and positive (blue) regions corresponding to the same absolute value of potential.Isopotential pairs are ±0.03(second row) and ±0.3 V (bottom).A molecule is depicted for position reference.

Figure 4 .
Figure 4. (a) PF6 around its C4 axis with MEP surface color-coded on the isodensity surface of 0.0004 e/bohr 3 ; and (b) the P F  anion and its MEP surface plotted on the electron density surface of 0.0004 e/bohr 3 .

Figure 4 .
Figure 4. (a) PF 6 around its C 4 axis with MEP surface color-coded on the isodensity surface of 0.0004 e/bohr 3 ; and (b) the PF 6 − anion and its MEP surface plotted on the electron density surface of 0.0004 e/bohr 3 .

Figure 4 .
Figure 4. (a) PF6 around its C4 axis with MEP surface color-coded on the isodensity surface of 0.0004 e/bohr 3 ; and (b) the P F  anion and its MEP surface plotted on the electron density surface of 0.0004 e/bohr 3 .

Figure 10 .
Figure 10.Li3PS4 studied as a crystal on (001) surface (a) Isopotential ±3.1 V surface of Li3PS4 separated between negative (red) and positive (blue) regions; (b) Li3PS4 atomic conformation and isopotential −3.1 V; (c) top view of the crystal of Li3PS4 and isopotential ±3.1 V; and (d) top view of the crystal of Li3PS4.Solid purple spheres represent lithium, yellow for sulfur, and orange represents PS4 tetrahedron.

Figure 10 .
Figure 10.Li3PS4 studied as a crystal on (001) surface (a) Isopotential ±3.1 V surface of Li3PS4 separated between negative (red) and positive (blue) regions; (b) Li3PS4 atomic conformation and isopotential −3.1 V; (c) top view of the crystal of Li3PS4 and isopotential ±3.1 V; and (d) top view of the crystal of Li3PS4.Solid purple spheres represent lithium, yellow for sulfur, and orange represents PS4 tetrahedron.

Figure 10 .
Figure 10.Li 3 PS 4 studied as a crystal on (001) surface (a) Isopotential ±3.1 V surface of Li 3 PS 4 separated between negative (red) and positive (blue) regions; (b) Li 3 PS 4 atomic conformation and isopotential −3.1 V; (c) top view of the crystal of Li 3 PS 4 and isopotential ±3.1 V; and (d) top view of the crystal of Li 3 PS 4 .Solid purple spheres represent lithium, yellow for sulfur, and orange represents PS 4 tetrahedron.

Figure 11 .Figure 11 .
Figure 11.Isopotential surfaces of solid electrolyte Li7P2S8I, separated between negative (red) and positive (blue) regions.(a) Shows the +0.5 and −3.8 V isosurfaces along the b-axis perpendicular to the crystal surface; (b) atomic conformation with the negative isopotential; (c) front view of the isopotential surfaces along b-axis; and (d) atomic conformation front view.Solid purple spheres represent lithium, green for iodine, yellow for sulfur, and orange for the PS4 tetrahedron.

Figure 11 .Figure 12 .
Figure 11.Isopotential surfaces of solid electrolyte Li7P2S8I, separated between negative (red) and positive (blue) regions.(a) Shows the +0.5 and −3.8 V isosurfaces along the b-axis perpendicular to the crystal surface; (b) atomic conformation with the negative isopotential; (c) front view of the isopotential surfaces along b-axis; and (d) atomic conformation front view.Solid purple spheres represent lithium, green for iodine, yellow for sulfur, and orange for the PS4 tetrahedron.