Cationic vacancies as defects in honeycomb lattices with modular symmetries.

Layered materials tend to exhibit intriguing crystalline symmetries and topological characteristics based on their two dimensional (2D) geometries and defects. We consider the diffusion dynamics of positively charged ions (cations) localized in honeycomb lattices within layered materials when an external electric field, non-trivial topologies, curvatures and cationic vacancies are present. The unit (primitive) cell of the honeycomb lattice is characterized by two generators, [Formula: see text] of modular symmetries in the special linear group with integer entries, corresponding to discrete re-scaling and rotations respectively. Moreover, applying a 2D conformal metric in an idealized model, we can consistently treat cationic vacancies as topological defects in an emergent manifold. The framework can be utilized to elucidate the molecular dynamics of the cations in exemplar honeycomb layered frameworks and the role of quantum geometry and topological defects not only in the diffusion process such as prediction of conductance peaks during cationic (de-)intercalation process, but also pseudo-spin and pseudo-magnetic field degrees of freedom on the cationic honeycomb lattice responsible for bilayers.

Introduction

Layered materials tend to exhibit a myriad of intriguing crystalline symmetries and topological characteristics based on their two-dimensional (2D) geometries and defects.[1,2] Since the discovery of graphene-based systems[3] and layered materials[1] such as honeycomb layered materials[2], a great deal of experimental and theoretical studies has been dedicated to illuminating the role the honeycomb lattice plays in the dynamics of electron quasi-particles and spin degrees of freedom, enriching our understanding of phenomena in materials ranging from high-temperature superconductors and 2D quantum hall systems to topological insulators and Kitaev materials.[2,4-6] Despite their crystal-structural versatility and compositional tuneability attracting interest in various realms of solid-state (electro)chemistry, materials science, condensed matter physics and pioneering the discovery of next-generation energy storage materials[2,7-12], theoretical and experimental studies centered on the honeycomb lattice have focused mainly on its effect on 2D electron and spin dynamics, thus rendering the behavior of larger particles such as positively-charged ions (cations) on the lattice understudied.

In particular, while the honeycomb lattice of graphene is formed by the carbon atoms, the electrodynamics studied is often centered around electron quasi-particles and spin degrees of freedom even when curvatures and topological defects are considered.[3,13,14] By contrast, cations in layered materials can not only form the honeycomb lattice, but are also responsible for the electrodynamics during cation (de-)intercalation process.[15,16] This poses a unique challenge to identifying the relevant quantum electrodynamics of cations in layered materials, not necessarily faced by other 2D systems. Intuitively (and in the continuum limit of the lattice), the quantum problem of electron dynamics in graphene with curvatures and defects is analogous to the problem of 1 + 2 dimensional quantum electrodynamics in a curved spaces, whereas the quantum problem of cations (charged) or their vacancies (uncharged) is analogous to the problem of 2D quantum gravity.[17,19,20]

A specific class of honeycomb layered materials based on transition metal and semi-metal oxides has recently emerged adopting, inter alia, chemical compositions embodied mainly by , or wherein A represents an alkali-ion (Li, Na, K, etc.) or coinage metal ions such as Ag, whereas M is mainly a transition metal species such as Co, Ni, Cu, Zn, etc. and D depicts a pnictogen or chalcogen metal species such as Sb, Bi, As and Te.[21-40] These materials are often referred to as honeycomb layered oxides.[2] The structure of these honeycomb layered oxides is comprised of localized and de-localized A cations sandwiched between slabs entailing M atoms coordinated with oxygen around D atoms, arranged in a honeycomb fashion (as shown in Fig. 1a). Thus, of specific interest to our work is the diffusion dynamics of the de-localized cations when an external electric field and non-trivial topology, curvatures and specific defects are present.[41,42]

Figure 1

(a) A schematic representation of the structure of exemplar honeycomb layered materials , or wherein A represents an alkali ion (Li, Na, K, etc.) or coinage metal ions such as Ag, whereas M is mainly a transition metal species such as Co, Ni, Cu and Zn, and D depicts a pnictogen or chalcogen metal species such as Sb, Bi and Te. The red rectangle at the base indicates the location of the schematic in Fig. 1b; (b) Schematic depicting honeycomb-shaped diffusion pathways with vacant cationic sites in exemplar honeycomb layered materials. The maroon line shows a possible random diffusion pathway of the cation.

Experimental investigations reveal the diffusion of cations to be largely restricted along honeycomb pathways in honeycomb layered tellurates such as (as shown in Fig. 1b).[15] While computational studies are consistent with this observation, they further suggest diffusion of cations in honeycomb pathways is restricted in honeycomb layered oxides for other materials as well, such as and where and Cs are cations with a large ionic radius, exhibiting a prismatic coordination with oxygen atoms of the Ni and Te octahedra forming the inter-layers.[16,43,44] In particular, Van der Waals forces and Coulomb repulsive forces tend to localize the cations in honeycomb lattices, creating a loosely-bound 2D non-Bravais hexagonal lattice with a two-cation basis known as the honeycomb lattice, which favors de-localization leaving vacancies in the hexagonal vertices only when sufficient activation energy from thermal fluctuations or the electric field, is present.[45,46] Consequently, provided the ground state of the system devoid of activation energies was initially vacancy-free, the number of vacancies, h in the honeycomb lattice is expected to closely correlate with the number of de-localized (mobile) cations, . Whence, the number of vacancies directly impacts the performance of the material as an effective cathode.[47]

Meanwhile, in thin layers of superfluids, superconductors and liquid crystals deposited on curved 2D surfaces, topological defects are known to couple to 2D curvature degrees of freedom, leading to the identification of the number of topological defects as the Euler characteristic of the surface.[48-54] This lends credence to analogous treatments for cationic vacancies in layered materials.[2,41] Moreover, layered materials demonstrating a bilayer arrangement of metal atoms (with each layer arranged in a triangular lattice) have been found, a vast majority being Ag-based layered oxides and halides such as (), , (or equivalently as ), , and more recently (where ).[55-66] Preliminary experimental and computational studies reveal that the bilayers represent a monolayer-bilayer phase transition of the honeycomb lattice, with the bifurcation mechanism not clearly understood.[63,67]

Herein, we consider how the honeycomb lattice of cations and emergent geometry constrains the model of cationic diffusion in such honeycomb layered oxides, leading to a rich topological description.[41] Consistently treating the number of vacancies, h as the number of holes and handles (genus) of an emergent 2D manifold with a conformal metric, we conclude that the primitive basis and the corresponding Euler characteristic must obey the modular transformation[68], invariant under the special linear group with integer entries, () with and , where is the total number of cation sites enclosed within the parallelogram with the primitive basis labeled by dx and dy as shown in Fig. 2.

Figure 2

Unit (primitive) cell () of localized cations arranged in a honeycomb lattice with lattice constant, a showing the primitive vectors, dx and dy. Every primitive cell is a parallelogram engulfing two (un-)occupied cationic sites, spanning partial interiors of four cationic-site hexagons.

In particular, the thermodynamics of the diffusive system of cations can be described by the partition function, of the form, which is invariant under the discrete re-scaling (, for large k) and discrete rotations () of the primitive cell for small k, generated respectively by,when the Euler characteristic can be written in a Fourier series,, is the p-th Betti number[69] associated with the n-th primitive cell and topology labeled by h, is an integer, M is the average potential energy of the cations (Strictly speaking, it can be taken as the inverse of the Green’s function of the system dominated by its potential energy, and hence can be complex-valued when the system has a finite life-time.[70]), is the ‘reduced’ inverse temperature with Boltzmann’s constant and is a constant independent of k but otherwise is dependent on the topology, h and temperature, .

The conformal geometry and resultant Betti numbers can be taken to emerge from the underlying matrix group field theory of the self-interactions and quantum correlations of the cations in each primitive cell. In particular, a matrix large group field theory can be considered, where is the size of the group.[17] Consequently, this construction effectively treats cationic vacancies in honeycomb layered materials as topological defects in the underlying field theory, related to modular symmetries of the honeycomb lattice in the context of emergent 2D quantum geometries.[19] For instance, we show that Eq. (2a) follows from the partition function of k pairs of cations forming the primitive cells of the honeycomb lattice, whereby each pair interacts via the Ising Hamiltonian due to emergent pseudo-spin and pseudo-magnetic degrees of freedom associated with the modular symmetry and broken conformal symmetry respectively. Analogous pseudo-degrees of freedom have also been considered in graphene-based systems.[13,14] Thus, the framework can be utilized to elucidate the molecular dynamics of the cations in exemplar honeycomb layered frameworks and the role of quantum geometry and topological defects not only in the diffusion process such as prediction of conductance peaks during the cation (de-)intercalation process, but also pseudo-spin and pseudo-magnetic field degrees of freedom on the cationic honeycomb lattice whose interactions predict cationic bilayered frameworks.[67]

Hereafter, we shall set Planck’s constant, the speed of electromagnetic waves in the material, , Boltzmann’s constant, and the elementary charge of the cations, to unity, , and employ Einstein summation convention unless explicitly stated otherwise.

The model

To build an intuitive geometric and topological picture of cationic diffusion, we begin by summarizing crucial results and clues from an idealized model, previously considered by the present authors in a separate publication.[41] Whilst the cations are positively charged, charge conservation requires the cationic vacancies created through de-localization by the electric field, present in the 2D plane to be considered electrically neutral. Nonetheless, the vacancies can be treated as possessing a fictitious ‘magnetic moment’ given by, where is the unit normal vector to the 2D layer comprised of the honeycomb lattice of cations. Thus, the cations diffusing through and around the vacancies created by the electric field, within the honeycomb lattice will introduce the ‘Aharonov-Casher’ phase[71], where is the boundary of a 2D patch, of the manifold associated with the 2D layer spanning the primitive cells, where locally the Cartesian coordinates are given by . Aligning the manifold with the plane, the unit normal vector becomes .

Applying Stokes’ theorem and substituting respectively Maxwell’s equations (Gauss’ law) and the normalization, where , is the 2D cation number density, is set as the integration cut-off scale in the z-direction and is the total number of mobile cations, yields the condition, where is also the number of vacancies.

Moreover, assuming the diffusion paths trace arc lengths defined by the 2D conformal metric, where is the 2D metric tensor, must satisfy Liouville’s equation[41], with K(x, y) the Gaussian curvature of the manifold. In fact, one can relate the two scalar functions, and by requiring that,Whence, by Stokes’ theorem, the Euler characteristic of the manifold is given by the Gauss-Bonnet/Poincaré-Hopf theorem,where we have used Eqs. (5), (6) and (7) to arrive at our result.

For instance, for a compact orientable 2D manifold homeomorphic to h number of simply-connected 2-tori, the Euler characteristic is given by, where h is the genus of the surface given by,which satisfies for a large number of diffusing cations, . Thus, this avails the avenue to treat the number of cationic vacancies as the genus, h which uniquely defines the emergent topology of the manifold. Moreover, using and with the 3D Levi-Civita symbol normalized as , Eqs. (3) and (7) follow from the phase equations of motion[72,73], on the Minkowski metric,where is the Minkowski metric tensor,are Maxwell’s equations, and are time-like and space-like unit normal four-vectors respectively, is the electromagnetic field strength, is the electromagnetic (U(1)) gauge field, is the dual field strength with the 4D Levi-Civita symbol normalized as and is the current density of the cations. Thus, by Eq. (10), the 2D charge density is related to the Gaussian curvature by[41],To incorporate diffusion in the formalism, we introduce the diffusion current given by, corresponding to Fick’s first law with D the diffusion coefficient and the center of mass momentum of the cations. Taking the cation number density to satisfy Boltzmann distribution at equilibrium, with a ‘gravitational’ potential energy governing the diffusion dynamics and a peculiarly defined center of mass effective mass using the diffusion coefficient, and applying Eq. (10) yields, Thus, Eq. (14) correspond to the Hamilton-Jacobi equations for the cations with corresponding to Hamilton’s principal function, the second equation to the 2D Langevin equation[74,75] and to the mean-free time/path between collisions (friction term). Thus, this serves as the motivation for appearing as the cut-off time and length scale in Eqs. (4) and (10). Moreover, the peculiar relation, can be better understood by applying the Virial theorem[76], where the averages are evaluated at equilibrium using,In this study, we shall consider the particular Hamiltonian for the cations, with momenta, , displacement vectors, , a mass per cation parameter defined as the inverse of the mean time/path between collisions, and the leading interaction term in the potential energy defined proportional to , the mobility of the cations. Typically, other terms such as the Vashishta-Rahman potential[77], which capture interactions of the cations with the slabs atoms especially oxygen, contribute higher order terms neglected herein. This requires that the diffusion coefficient, including cation-cation correlation terms[78], satisfy the Einstein-Smoluchowski relation, as , where we have used the result in Eq. (15a). Thus, requires we have .

Observe that, when cation-cation correlation () terms vanish, the diffusion coefficient becomes the self-diffusion coefficient, whereas takes the role of frequency of the harmonic oscillator. Moreover, since the mobility is a constant, we can re-define it as , where and a is taken to be the lattice constant with dimensions of length. We thus have, and , where G is a gravitational constant, in obvious comparison with Schwarzschild black hole thermodynamics.[18]

Results

Conductance spikes

We note that, Eqs. (10) and (13) require that the diffusion current takes the Chern-Simons form[79], where is the conductivity of a single primitive cell, is the mobility (Einstein-Smoluchowski equation) and,is the Chern-Simons level. However, the current (density) we are interested in is not necessarily the Hall current, but the spatial part of , which couples to the electromagnetic tensor, in Maxwell’s equations given in Eq. (11). Consequently, this predicts conductance spikes whenever primitive cells containing cation sites are activated in the (de-)intercalation process. However, due to typically low measurement sensitivity in existing experimental data, the integer nature of the conductance cannot be ascertained. Nonetheless, distinct conductance spikes can be observed at low resolution (i.e with ) at specific voltage values as shown Fig. 3.

Figure 3

(a) Current, I–Voltage, V characteristics derived from cyclic voltammetry experiment of the K cation extraction (charging) process with as the cathode in a two-electrode setup. K metal was used as the counter electrode and the scanning rate was set to 0.1 mVs[80]. The current peaks at varied spread-out voltage values reflective of the small activation energy of K cations in the honeycomb layered oxide material. (b) The conductance, of during the charging process displaying conductance spikes at varied spread-out voltage values corresponding to extractions of a large number of K cations from multiple honeycomb lattice layers, , where the values of and cannot be separately determined from the results. Nonetheless, the sharp conductance spikes occur evenly distributed at a rough interval of 0.1 V within the voltage interval, 3.2 V to 4.0 V (gray/shaded region)

The low resolution is an artifact of the multi-layered nature of the materials, with each honeycomb lattice not only contributing active primitive cells with cationic sites during the extraction but also multiple sites getting activated at once at spread out external voltage values. Moreover, experimental data suggests that a large activation energy, meV, where is the activation energy of potassium, K cations in the honeycomb lattice[46] and the presence of cationic vacancies before the extraction process would tend to disfavor the cation extraction process from occurring at evenly spread-out (low to high) voltage values during cycling, often leading to a solitary broad current peak centered at the high voltage regime, for instance, as can be seen in the I–V cycling characteristics for with , since the activation energy for Li and Na is vastly greater than that of K, i.e. .[2] Consequently, we have plotted only the Current, I – Voltage, V characteristics of the extraction process for in Fig. 3, which exhibits several distinct current peaks across varied voltage values.

Moreover, we can employ Eq. (8), which suggests the partition function,is given by the sum over different geometries of the manifold with distinct topology, where is the 2D Ricci scalar, is the 2D Ricci tensor and is the general form of the 2D Riemann tensor in Riemannian geometry with the equivalence, assumed to be valid. Evidently, Eq. (19) is the partition function of 2D quantum geometry in Euclidean signature[19], where the coupling constant corresponds to . It is worth noting that, considering emergent geometries within crystals to describe defects is not entirely a novel idea, since it has been considered in great detail for disclinations and dislocations within the context of classical geometries with torsion.[20,81-85]

Finally, we are left to show that Eq. (19) (and equivalently, Eq. (20)) obeys the necessary modular symmetries defined in Eq. (2b), which are imposed by the primitive cell of the cations in honeycomb layered materials. Nonetheless, Eq. (2) approaches Eqs. (19) and (20) in the limit , as required.

Modular symmetries

The honeycomb lattice is spanned by the primitive basis, dx and dy defining a parallelogram enclosing pair of cation sites. Since the unit cell in Fig. 2 is a rhombus, we have . Moreover, transforming the basis by the matrix, where, we find that and , given in Eq. (2b), correspond to the re-scaling, and the discrete rotation, , as illustrated in Figs. 4 and 5 respectively. Moreover, we shall consider the modular form[68] defined as, to completely characterize the honeycomb lattice in an invariant manner, under with . By definition, is invariant under the modular transformations, i.e. . Consequently, f(k) transforms as a modular form of weight 2,Proceeding, we must take the large limit, , which spans the entire honeycomb lattice. Moreover, assuming , we obtain, .

Figure 4

The honeycomb lattice of cations depicting the action of the generator of on the primitive vectors. (a) The primitive vectors dx and dy of the primitive cell of the honeycomb lattice, where is the number of pairs of cations enclosed within the primitive cell; (b) transformation corresponding to the re-scaling of the dx primitive vector and hence an expansion of the unit cell, .

Figure 5

The honeycomb lattice and discrete rotations () generated by acting on the primitive cell. The primitive cell is rotated as shown in (a), (b), (c) and (d) by application of the transformation, such that corresponds to inversion of (a) to (d) and (b) to (c), where is the identity matrix, requiring that .

Now, consider the diffusion dynamics of the cations given in Eqs. (6) and (14). Defining the velocity, , where is the arc length interval along the proper length, ds, we obtain, where we have used from Eqs. (17b) and (8). Thus, setting defines the modular form as the action of a particle of mass, M/2 in 2D Riemannian geometry, or equivalently the exponent of Eq. (19) (or Eq. (20)).

Liouville conformal field theory

To elucidate further properties of the formalism, we shall consider the Liouville action[86],where is the 2D metric tensor with given in Eq. (6), K is the Gaussian curvature associated with and is the Ricci scalar associated with , is a parameter dependent on k and genus h. Setting and , the metric in Eq. (24) reduces to the 2D identity matrix, requiring that the Ricci scalar vanishes, getting rid of the last term even when .

The Liouville action reduces to,where . Thus, for arbitrary , Eq. (25) can be varied with respect to to yield Eq. (6). Moreover, we define the path integral as, which, after summation over j and functional integration over yields, where, is the divergent vacuum energy of with and approximated as separate non-interacting fields, C is a constant that will be set to vanish by regularization, is the area element, , are the allowed momenta/energies of the bosonic field, , , and . Thus, Riemann zeta function regularization requires[87], where .

To elucidate the nature of Eq. (26), we consider the Virasoro algebra[88], spanned by two copies of commuting generators, and for all integers , where , is a real-valued constant (the central charge), satisfying and is the Kronecker delta. Meanwhile, the representations are characterized by a highest weight primary state, satisfying, , or , , , and , . The rest, , and , can be computed by applying the Virasoro algebra in Eq. (28).

Proceeding, it is known that the field is primary when the conformal dimension is given by[89,90], while the marginal condition for the primary field that guarantees conformal invariance of the theory is with which yields,with and . The central charge is given by, which can be written as, where is the temperature, is a defined energy in order for as required. The conformal dimension and spin are given by, respectively. However, the theory is known to be unitary only for and corresponding to and respectively. Since k is considered real, we must have . Thus, not only breaks discrete rotation symmetry but also breaks unitarity. Indeed, unlike Eq. (30), is not invariant under , which breaks the discrete rotation symmetry generated by , where is the number of primitive cells in the honeycomb lattice.

Now, we consider the partition function,where , and , , arriving at the energy, with from Eq. (17b). Consequently,where we have used Eq. (32) and defined the momentum of the primary field, as, Recall that the marginal condition is guaranteed by , which now translates to,for which yields, or equivalently .

In this case, unitarity is only achieved for , which corresponds to the two-sphere satisfying . Moreover, to reconcile with Eq. (26) the extra factor of must correspond to the vacuum energy, after regularization. This requires the momenta be positive integers, which yields the expression by regularization. Note that, for , all other positive values of h cannot satisfy the marginal condition and hence must break conformal invariance.

Discussion

If the Euler characteristic of the n-th primitive cell can be associated with a manifold defined by the Poincaré polynomial,where is the p-th Betti number, the Euler characteristic of emergent manifold can be calculated as the Euler characteristic of the connected sum of emergent manifolds corresponding to the primitive cells, with topology, h, which can be decomposed into a Fourier series given by Eq. (2c), where with and we have used the fact that only the -st Betti numbers are additive in a connected sum. Consequently, the simplest real-valued partition function that respects the modular symmetries, and in their respective limits of k and is proportional to the partition function in Eq. (19) (and hence Eq. (20)) corresponds to Eq. (2), where the limit breaks the discrete rotation symmetry, of the partition function, whilst promoting scale invariance, . Mathematically, this is required since there exists no holomorphic modular forms of weight 2, invariant under both transformations.[68]

Nonetheless, this can be remedied by considering the Euler characteristic proportional to the almost holomorphic modular form of weight 2 in the large genus limit ()[68],where ,with the second Eisenstein series, and the sum of the divisors of the positive integer, n. Thus, the Betti numbers for the emergent geometry of the honeycomb lattice primitive cells correspond to the two-torus, where for such h primitive cells in a connected sum, we set, Thus, unlike in Liouville CFT, modular invariance and hence conformal invariance is guaranteed for all integer values of h and k. However, it is broken for the primitive cell when , corresponding to a phase transition.[67]

To discuss the effects of such a phase transition, we note that, , where,which implies the unit basis acquires a minus sign under . Since exchanges the basis dx with dy and vice-versa, it corresponds to a discrete rotation when acting on the primitive cell. There are 4 such discrete rotations such that correspond to complete rotations of the primitive cell, where is a real number. In addition, exchanges one cationic site in the primitive cell with the other. Thus, under the exchange of two cations belonging to the same primitive cell, the transformation picks up a minus sign. This can be understood as the origin of the pseudo-degree of freedom we shall refer to as pseudo-spin, which distinguishes the two sub-lattices of the honeycomb lattice.[13,67]

Under specific conditions, the pseudo-spin of the graphene lattice can be linked to the spin of the electrons localized on the carbon atoms in the sub-lattice.[13] However, for cations in honeycomb layered oxides, no such identification can be affirmed. Nonetheless, the pseudo-spin degree of freedom, coupled with the group imply the partition function given in Eq. (2a) that the underlying theory of cations is a conformal field theory whose ground state must avoid pseudo-spin frustration by pseudo-spin anti-ferromagnetic behavior, but nonetheless prevents the cations from forming a stable honeycomb lattice due to a repulsive exchange interaction which can be offset by pairing of opposite pseudo-spin degrees of freedom.[67]

The effective theory for two pseudo-spin cations () in non-interacting honeycomb primitive cells corresponds to the 1D Ising Hamiltonian[91],where, is the pseudo-magnetic field[14,67] in the z-direction interacting with the pseudo-spins, which is taken to be proportional to the Euler characteristic, , while A(h) is the Heisenberg term representing the exchange interaction, assumed to depend on the genus, with the cationic vacancy number.

This Ising model is exactly solvable, where standard calculation for the partition function yields[92],where is the trace over the genus h and spins , and are the eigenvalues of the transfer matrix,given by,with,Thus, the non-interacting system of k primitive cells occupied by pairs of cations interacting via their pseudo-spins and the pseudo-magnetic field yields the partition function in Eq. (2a) where , which takes on varied values for different topology configurations, h. Moreover, the exponents of components of the transfer matrix correspond to the pseudo-spin energy states, where corresponds to a gapped phase where the honeycomb lattice bifurcates into bilayers with energies and due to a finite pseudo-magnetic field, , whereas correspond to the ferromagnetic () and anti-ferromagnetic () alignment of the pseudo-spins.

To avoid pseudo-spin frustration when (two-torus, ), the honeycomb lattice must be anti-ferromagnetic described by the singlet bound state, , (), with , which disfavors the ferromagnetic condition.[93] Nonetheless, for a finite pseudo-magnetic field, () the triplet bound state, , ) is allowed, corresponding to other topology configurations. We are interested in , corresponding to the unitarity condition for marginal fields in Liouville conformal field theory. In this case, is the potential energy of the resultant bond, is the 3D relative position of the pseudo-spin up and down cations, , is the Gaussian curvature of the two-sphere () and due to the point like nature of the Fermi surface, is taken to be the non-oscillatory Ruderman-Kittel-Kasuya-Yosida (indirect exchange) interaction mediated by the conduction electrons of the cations, where is a distance scale to be determined.[94,95] Thus, the energy density (integrand), corresponds to a metallophilic interaction[67] between the pseudo-spins in the primitive cell separated by a distance, r apart in a stable bond forming bilayers. The bilayers are stable when and , corresponding to the separation distance between the honeycomb sub-lattices, which can be determined experimentally.[63,67] Consequently, this finite distance scale breaks scale/conformal invariance of the theory.

In conclusion, we have constructed a consistent framework to treat cationic vacancies in honeycomb layered materials as topological defects, h, by relating the Euler characteristic of the manifold to modular symmetries and 2D quantum geometries.[19] The framework predicts integer conductance spikes during (de-)intercalation process, proportional to the number of active cation sites, participating in the diffusion process at high resolution, , which remain unobserved. Nonetheless, the framework greatly elucidates the geometric nature of the diffusion process which occur in these novel materials, and the crucial role played by cationic vacancies as topological defects, and hence should find great utility in finding avenues for performance optimization of such cathode materials for energy storage.[42,43] Further theoretical, computational and experimental treatments and applications are beyond the scope of the present work.[2,17,67]