Time-reversal symmetry breaking in a square lattice

The bulk conductivity of a two-dimensional system is studied assuming that quantum interference effects break time-reversal symmetry in the presence of strong spin-orbit interaction and strong lattice potential. The study is carried out by direct diagonalization in order to explore the nonlinear-response regime. The system displays a quantized conductivity that depends on the intensity of the electric field and under specific conditions the conductivity limit at zero electric field shows a nonvanishing value.


I. INTRODUCTION
An essential result from quantum mechanics prescribes that when two operators commute there exists an eigenbasis that diagonalizes them simultaneously, so that the elements of such an eigenbasis conform at the same time to both operators. It is however important to highlight that this principle does not dictate that any eigenbasis of the first operator is also an eigenbasis of the other, which would automatically imply that any eigenstate of the first operator should conform to the second one. Nevertheless, such an implication takes place in one instance: When the spectrum of the first operator is nondegenerate. These facts are at the center of the theory of (standard) phase transitions: In a scenario where the Hamiltonian commutes with an unitary operator (generated by the symmetry) a critical point separates a trivial or symmetric phase, where the state conforms to both the Hamiltonian and the symmetry, from a non-trivial or broken phase, where the physical state is no longer a symmetry eigenstate. A transition of this kind is only possible when the Hamiltonian spectrum goes from nondegenerate to degenerate, being the latter case the only one where the physical state can break the symmetry. From this perspective a phase transition is essentially the arising (or suppression) of the Hamiltonian's degeneracy. The mechanisms by which the symmetry is broken must be on the one hand irreversible [1], since otherwise equilibrium states would retain the Hamiltonian's symmetries, and on the other hand global, so that they affect the state as a whole and the symmetry be broken everywhere. Without these mechanisms the symmetry would not break and the phase transition would not take place. A typical example of this kind of transition is the change from paramagnetic (symmetric) to ferromagnetic (broken) in spin systems. In the case of electron systems, a paramount result known as the Kramers degeneracy [2] has significant implications in connection to phase transitions: Time-reversal spin-systems have de- * Electronic address: kfjimenezfals@gmail.com † Electronic address: reslenjo@yahoo.com generate spectra. As such, the notion of phase transition in this kind of systems cannot be accommodated in the standard symmetry-breaking paradigm associated with a transition from degenerate to non-degenerate or vice versa. However, it has been observed that a phase transition can take place whereupon the symmetry is broken locally in the non-trivial phase, i.e., some parts of the state, usually those associated with the system's bulk, display symmetry, but others, like those associated with the system's boundary, do not. This phase is known as the topological insulator [3][4][5]. In contrast, the trivial phase lacks symmetry-breaking mechanisms entirely, thus being known as the standard insulator. The symmetry associated with the transition is time-reversal (TR) [2], while the symmetry-breaking mechanism is quantum interference, triggered by spin-orbit interaction. In order to appreciate the mentioned mechanism with some degree of formality, let us imagine a situation in which a spin-up electron moving forward bumps elastically into an obstacle, as portrayed in figure 1. The incident particle can be scattered in different ways but the case is such that the only state available for backscattering corresponds to spin down. The obstacle effect can be modeled by the following term Notice that the effect is such that rotations in opposite directions by and angle θ are equally considered. This means the potential is TR invariant, i.e., it does not mutate under the change t → −t [17]. This case corresponds to a non-magnetic obstacle. The problem can be approached by considering solutions on each side of the obstacle following the decomposition shown in figure 1.
Demanding wave function continuity, ψ I (0) = ψ II (0), it is found that A = A ′ and V = V ′ . Integration of the Schrodinger equation around the origin yields For finite values of U the only nontrivial solution is for V = 0 and θ = π, which represents a spin-up state with perfect conduction in spite of the particle hitting an obstacle. Contrary to potential (1), this perfect conduction state is not TR invariant. This effect is produced by the destructive interference of backscattering paths and as such is a quantum-mechanical effect [4]. This result is sustained by the fact fact that backscattering is only possible via spin inversion. The existence of spinflipped channels in opposite directions is guaranteed by the Kramers degeneracy as long as the Hamiltonian be TR invariant. The other determining factor is that the number of equal-energy moving-channels in a given direction be odd, i.e., energy degeneracy must be an odd multiple of two. Contrariwise, the suppression of backscattering paths due to quantum interference is unfeasible when the number of moving channels in one direction is even. This can be seen considering a Hamiltonian that displays the common form (the constant term σ 2 is included only to keep a reference to spin states) It then followsĤ| − p =Ĥ| p . As a result, given an eigenstate with average momentum p it is always possible to construct another eigenstate with the same energy, spin and spatial distribution, but opposite momentum − p , providing in this way direct backscattering channels. In this example the number of moving channels on each direction is always even, because every solution admits spin up and down. The correspondence between dissipation when there is an even number of moving channels and conduction when there is an odd number of moving channels can be represented using the members of the group Z 2 = {0, 1}. This equivalence has prompted the use of the adjective "topological" when referring to the case of nonvanishing conductivity. Also, it has been shown that it is possible to formally establish the classification of a given system from its Block energy-structure using topology methods [6]. One way of inhibiting even numbers of moving channels is to deliberately break the model's TR invariance. The simplest strategy consists in applying a magnetic field to a spinless electron In a two-dimensional space with perpendicular magnetic field the vector potential is given by Â = (0, eBx). The associated term breaks the parity symmetry since it is not invariant under the change x → − x. It can therefore be said that neither parity nor TR are independently preserved. However, assuming that U (− x) = U ( x), it is noticeable that the Hamiltonian is invariant under parity and TR applied simultaneously. In absence of symmetry breaking mechanisms, it is at least possible to build eigenstates that be invariant under the joined effect of both symmetries. These states take place both in the bulk and on the edges, as shown in figure 2. The interesting thing about the states that develop on the surface is that backward-moving states shift from forwardmoving states, leaving only one moving channel on each edge and allowing the arising of an energy gap associated with backscattering as other elements of the problem, as for example an electric field or electron-electron collisions, are taken into account. As backscattering channels become sufficiently suppressed, motion on the edges becomes dissipationless. As can be seen, the conduction mechanism is rooted in the local breaking of the TR symmetry rather than in quantum interference, since in this example the Hamiltonian is not TR invariant. It is known that at low temperatures the transverse conductivity of (6) comes in integer multiples n of e 2 /h [7,8]. This phenomenon is known as the integer quantum Hall effect [9]. Number n has been shown to correspond to a topological invariant [10], which explains the notable robustness of the effect observed in experimental measurements. In a topological insulator dissipation channels are displaced just as in the Hall effect, but the role of the magnetic field is taken over by the spin, such that A. In the system's boundary the number of moving channels in one direction is odd. And B. The Hamiltonian is TR invariant. These facts combined lead to the quantum interference that provokes the state to break TR symmetry and so become conducting [11,12]. It is often said in this regard that conduction is "protected" by the TR symmetry. Topological insulators are peculiar in that they can conduct even though their bulk spectra are gaped.
The goal of the present study is to provide a numerical analysis of a single-body model where the symmetrybreaking takes place in the bulk and is generated by a strong lattice potential that acts in a way analogous to the edge of a topological insulator. It is of interest to con-sider the electric field as an integral part of the problem and to observe how the conductivity depends on this field beyond the linear approximation. This approach intends to shed insight by helping visualize the system's response as a complement to the more abstract analytical formulation often found in related studies. Interestingly, this procedure yields a quantized conductivity that shows a dependence with the number of bands below the Fermi energy and in some cases this conductivity remains finite as the electric field goes to zero, suggesting in this way a superconducting state.

II. THE MODEL AND ITS EIGENSYSTEM
The model corresponds to an electron that moves on a two-dimensional potential U (x, y) under the action of an electric field E in the y direction. Spin-orbit interaction arises as a coupling between the z-components of spin and angular momentum. The Hamiltonian is written aŝ Constants m and e represent mass and charge respectively. The intensity of the spin-orbit interaction is mediated by constant λ. The potential is written as Current technology allows for a high degree of control over the model's parameters in optical lattices or superlattices [13,14], being these scenarios where the effects reported further ahead are more likely to be observed. For a numerical analysis it is necessary to bound the system in order to provide a compact Hilbert space, hence periodic boundary conditions are imposed on the x-axis over a square lattice of side L. Lattice constant a is such that L = N a, being N the square root of the total number of real unit-cells in the lattice. It can be noticed that due to the terms of spin-orbit and electric field, Hamiltonian (7) does not display translational invariance in neither axis and therefore it does not admit a treatment in terms of Bloch functions. However, it is possible to consider an alternative symmetry arising from simultaneous translations of space and momentum, but for this it is necessary to add a term, as followŝ The extra-term can be considered either as a physical confining potential, in which case it becomes an integral part of the Hamiltonian, or as a perturbation. Both (7) and (9) are TR invariant, but only (9) commutes with the following symmetry operator ( = 1) This can be confirmed usingTxT −1 =x + a and Tp yT −1 =p y − aλσ z . Fundamental results dictate that there exists a common basis forĤ 0 andT . The most general way of writing an eigenfunction ofT with eigenvalue e ika is ψ(x, p y , σ z ) = e ikx u k (x, p y , σ z ), (11) subject to the condition u k (x + a, p y − aλσ z , σ z ) = u k (x, p y , σ z ). These requirements are met for functions defined as insofar as being n an arbitrary integer. Solving for w yields being q = nN − j. Integer j is not bounded, since the corresponding momentum eigenvalue p x = 2πj L in (12) remains always consistent with boundary conditions. This is not the case for the position eigenvalue y = 2πw because the system is bounded on the y axis, therefore Using (14) it then follows Since q in (14) can take negative values, the total number of position states is given by Q = 2q max + 1. Inserting Q and solving for L gives As a consequence, the system length depends on the number of states and the interaction constant. This conditioning certainty arises from the symmetry and seems to be related to the fact that Heisenberg's uncertainty principle establishes a phase space grating. Another result is that the system's eigenfunctions are periodic in the p y -space with period λL, as can be seen from (14) and (12). Besides, periodic boundary conditions on the x axis ψ(x+L, p y ) = ψ(x, p y ) determine as valid values of k = 2πl L , for l integer. The size of a unit cell in k-space is 2π a . The eigenvalue problem can be formulated in terms of the symmetry functions as wherein being E the system's energy. As the problem is separable with respect to spin, it is valid to setσ z = 1 in order to focus on the spin-up case. Technically, the resulting expression breaks TR symmetry, but the operation is justified by the fact that quantum interference breaks such a symmetry through the same quantum interference effect that takes place in the boundary of a topological insulator, with the difference that in this case the role of the boundary is taken over by the periodic potential, which must in addition be sufficiently strong to induce trajectories surrounding potential peaks in the bulk, as sketched in figure 3. This mechanism must derive in an energy spectrum whose degeneracy comes in odd multiples of two, which is to be checked further ahead. Sincê x andp y are to be used as a complete set of commuting observables, the following results must be considered, p x = −i∂ x , y = i∂ py . Basis functions are taken as eigenfunctions ofp x andŷ normalized over the xp y -cell Calculation of the matrix elements, j ′ , q ′ |Ĥ ↑ k |j, q , yields if n = n ′ , and g(n, n ′ ) = 0 if n = n ′ , Interestingly, matricial elements are all real even though the basis functions are complex. This helps reduce computation costs. When this Hamiltonian is numerically diagonalized the respective eigenfunctions take the next form The momentum integer is j = N n − q and the total number of momentum states is J = 2n max + 1. Figure 4 presents the first bands of the spin-up-part of th Hamiltonian for a set of reasonable parameters and strong lattice potential. In addition to being flat, the band pattern is in general nondegenerate, except for specific values of E. This can be better seen on the left side of figure 5, which shows energy as a function of E. Given the TR symmetry of the whole Hamiltonian together with its degeneracy profile, it is viable to assume that the quantum interference mechanisms that take place in the boundary of a topological insulator also take place in this system.

III. CONDUCTIVITY
The mean values of momentum and position over state (11) are Taking the effect of electron-electron collisions as a perturbation, it can be said that at zero temperature the contribution of a given band to the x-conductivity is at first order proportional to the sum of momentum mean values over that band The transverse conductivity is proportional to the sum of contributions from all the bands below the Fermi level Being α a proportionality constant. Figure 5 shows energy as well as conductivity against electric field. It can be seen that σ xy comes in integer multiples of a constant that shows no dependence with the electric field. This conductivity displays a stair pattern, being constant over intervals of different extension and increasing by integer steps of different size in the same points where a given band crosses the corresponding Fermi energy as the electric field grows. The stair pattern became less appreciable in simulations with smaller λ. In the quantum Hall effect, the integer that determines the conductivity is given by the number of times the wave-function phase winds around the boundary of a two-dimensional Brillouin zone [7,10]. Such an integer is known in topology as the Chern invariant [3]. This parallel does not apply here since the inclusion of the electric field in the Hamiltonian breaks translational invariance in the y-axis and the reciprocal lattice becomes one-dimensional. Whether there are additional topology constructs that apply in this context or the system's discreteness can be ascribed to deeper precursors remains to be seen. Whatever the case, this result shows that the conductivity's quantization does not depend on the linear-response assumptions necessary to obtain the Kubo formula [7,15]. Another curious trait of figure 5 is that for the lowest two values of Fermi energy conductivity features finite limits at zero electric field. For this to happen the Fermi energy must equal one of the system's energy values for vanishing electric field. This limit is consistent with a superconductor state because it features dissipationless conduction on account of the quantum interference effect already discussed in this document. It would then be interesting to add a magnetic field to Hamiltonian (9) and see whether this magnetic field is offset by the induced magnetization, giving in this way a solid evidence of the superconducting state. Contrary to transversal conductivity, the first order longitudinal conductivity vanishes for any electric field. This happens because the system displays backscattering channels in the y-axis, These channels are nonetheless displaced, so that TR symmetry is broken locally just as in the Hall effect. By a similar mechanism, longitudinal conductivity could be induced adding an electric field in the x-axis. Likewise, the inclusion of interaction terms directly in the Hamiltonian would open a gap between the states involved in (28), since dissipation channels are spatially separated and a particle would normally experience collisions in going from x to −x in proportion to the magnitude of x. The same phenomenon can also explain the longitudinal conduction in topological insulators [16]. Although electron-electron interaction was not considered in this work, the single body functions obtained in section II are the starting point to build a second-quantization Hamiltonian giving a more accurate representation of the system, in which case it is reasonable to expect nonvanishing longitudinal-conduction.

IV. CONCLUSIONS
A single-body spin-orbit interaction model has been used to study the conductivity pattern produced by a symmetry-breaking mechanism that takes place in the system's bulk. The electric field responsible for charge transport has been included in the Hamiltonian and the study has been carried out by direct diagonalization in order to explore the system's response beyond the linear approximation. Such a response displays a discrete pattern that is consistent with the quantization of conductivity over the range of fields considered in the study, showing in this way that a quantized conductivity does no depend on linearity approximations. The quantization value is found to depend on the number of energy bands located below the system's Fermi energy and also on the intensity of the electric field. In the particular case of the Fermi energy being equal to an eigenenergy of the zero-field Hamiltonian, the limit of conductivity for vanishing electric fields proves to be finite. This feature is consistent with a superconducting state but additional tests are necessary to confirm such a hypothesis. Were it verified, it would mean the mechanism by which the TR symmetry is broken in a superconductor can be related to the one taking place on the edges of topological insulators. Overall, both the system's physics as well as the perspective granted by the numerical analysis display interesting features.

V. ACKNOWLEDGMENTS
Financial support from Vicerrectoría de Investigaciones, Extensión y Proyección Social from Universidad del Atlántico through the program Investigación Formativa is gratefully thanked by both authors.