Random physics

Alberto Verga, research notebook

\(\newcommand{\I}{\mathrm{i}} \newcommand{\E}{\mathrm{e}} \newcommand{\D}{\mathop{}\!\mathrm{d}} \newcommand{\Di}[1]{\mathop{}\!\mathrm{d}#1\,} \newcommand{\Dd}[1]{\frac{\mathop{}\!\mathrm{d}}{\mathop{}\!\mathrm{d}#1}}\)

»quantum chaos»kicked rotator»dynamical localization»random matrices

Quantum walk

In a classical random walk the moving particle explores the configuration space at a rate growing as the square root of the number of steps \(\sim \sqrt{t}\). Instead, a quantum walk, as first defined by Aharonov et al (1993), can explore its configuration space at much faster rate \(\sim t\). This is a consequence of the superposition principle, at each time step, after the result of a coin toss, the moving particle initially at site \(n\), simultaneously occupies the two neighboring sites with probabilities given by the square of the respective site amplitudes.

A step of a quantum walk is the result of the composition of two unitary operators, which gives an unitary evolution operator and a final measurement step:

  • Coin toss: coin tossing corresponds, in the simplest case, to a two state quantum system, \(|\uparrow\rangle\) (up) and \(|\downarrow>\) (down); commonly used coin \(C\) operators are Hadamard \(A\) and spin rotation \(R\):

    $$A = \frac{1}{\sqrt{2}} \begin{pmatrix} 1 & 1 \\ 1 & -1 \end{pmatrix} \,, \quad R(\theta) = \E^{-\I \sigma_y \theta/2} = \begin{pmatrix} \cos (\theta/2) & -\sin (\theta/2) \\ \sin (\theta/2) & \cos (\theta/2) \end{pmatrix}\,. $$
    The action of both coin operators results in a superposition of the two up and down basis states, from an initial up or down state.

  • Shift: the motion step corresponds to the application of a shift operator (one step translation in the appropriated space), conditioned to the result of the coin tossing:

    $$S=\sum_x \left( |x+1\rangle \langle x|\otimes |\uparrow\rangle\langle\uparrow| + |x-1\rangle \langle x|\otimes |\downarrow\rangle\langle\downarrow| \right) $$
    where \(\otimes\) denotes the Kronecker product; \(S\) shifts the position \(|x\rangle\) to the right \(|x+1\rangle\), if the coin state is \(|\uparrow\rangle\) or to the left \(|x-1\rangle\), if it is \(|\downarrow\rangle\).

  • Evolution: the walk dynamics is governed by the unitary operator,

    $$U = S(I \otimes C)$$
    that is the composition of \(S\), a shift operator acting on a discrete space (position, energy, angular momentum states, etc.), and \(C\) the coin operator acting on an internal degree of freedom; \(I\) is the identity. Repeated application of \(U\) leads to the history
    $$|\psi_t\rangle = U^t |\psi_0\rangle$$
    of the walker.

  • Measurement: the last step of the quantum walk algorithm consists in some projection to a specific state \(|s\rangle\). The measure of \(|s\rangle\) provides the probability of this state \(|\langle s| \psi_t \rangle|^2\).


Write a code to simulate the Hadamard random walk. Show that the evolution of a state initially at \(x=0\) depends on the spin state: compare the sequences with initial coin states \(|\uparrow\rangle\) and \((|\uparrow\rangle + \I |\downarrow\rangle)/\sqrt{2}\). Verify that \(\langle x_t| x^2 |x_t\rangle\sim t^2\), where \(|x_t\rangle\) is the state \(|\psi_t\rangle\) projected to any coin state (summed over coin states), at time \(t\) for large \(t\).

Topology of the quantum walk

One of the goals of condensed matter physics is to characterize the different phases or states of matter: insulator, conductor, ferromagnetic, superfluid, etc. Landau establishes a theory that can account of these states and its changes, using the concept of order parameter, associated with the symmetries of the systems and the mechanisms that can brake these symmetries. A symmetry is broken, for instance, under the disorder-order transition, between a paramagnetic state and a ferromagnetic state when the temperature decreases. The isotropy of the magnetization field (the order parameter) in the paramagnetic state, is replaced by a magnetization directed along some specific direction, characteristic of the ferromagnetic state.

However, after the discovery of the quantum Hall effect, it became clear that other order states exists, not related to symmetries, but with the topology of, for instance, the electronic bands. Edge states in the integer and fractional Hall effect, solitons in one dimensional molecular chains (polyacetylene), or spin-polarized currents in topological insulators, are manifestations of some topological property of matter.

In order to identify the topology of a material, we use as an example the Bloch theory of electronic bands applied to a two dimensional material. We consider \(|u_{\boldsymbol k}\rangle\), the Bloch function of quasi-momentum \(\boldsymbol k\), and the system’s Hamiltonian \(H(\boldsymbol k)\) which is a function of the quasi-momentum vector inside the Brillouin zone, \(\boldsymbol k \in \mathrm{BZ}\). We assume that the material is an insulator whose energy bands, the set of eigenvalues of \(H\), \(E(\boldsymbol k)\), are separated by a gap. The topology of the band structure is related to the mapping from the crystal momentum space (the Brillouin zone), a torus in two dimensions, to the Hilbert space of the gapped Hamiltonian. If the gap remains open under changes of the quasi-momentum, the Hamiltonians belongs to the same topological class. These topological classes possess a Chern invariant \(C \in \mathbb{Z}\), given by the (gauge independent) Berry flux integral:

\begin{equation} C = \frac{1}{2\pi} \oint \D^2 \boldsymbol k \cdot \boldsymbol F\,, \quad \boldsymbol F = \frac{\partial }{\partial \boldsymbol k} \wedge \boldsymbol A\,,\quad \boldsymbol A = -\I \langle u_{\boldsymbol k}| \frac{\partial }{\partial \boldsymbol k} |u_{\boldsymbol k}\rangle \end{equation}

where the Berry connexion \(\boldsymbol A\) and the curvature \(\boldsymbol F\) are the analogs of the vector potential and the magnetic field; this analogy comes from the indetermination of the eigenvector’s phase \(|u_{\boldsymbol k}\rangle \rightarrow \E^{\I \phi(\boldsymbol k)} |u_{\boldsymbol k}\rangle\), which translates into the gauge transformation \(\boldsymbol A \rightarrow \boldsymbol A + \partial_{\boldsymbol k} \phi(\boldsymbol k)\). The Chern number itself, appears as the charge of a magnetic monopole, in this analogy.

Su Schrieffer Heeger model

The simplest non trivial topological system is a two-band model with hopping amplitudes between sites of type A and B:

\begin{equation} H = \sum_{x=1}^{L} (\epsilon + \Delta) c^\dagger_{\mathrm{A},x} c_{\mathrm{B},x} + \sum_{x=1}^{L-1} (\epsilon - \Delta) c^\dagger_{\mathrm{A},x+1} c_{\mathrm{B},x} + \text{hc}\,, \end{equation}

we take the length of the unit cell equal to one. This model was first introduced by Su, Schrieffer and Haager (1979) to investigate the charge transfer in polyacetylene, a (CH) chain with alternating double bonds. \(H\) takes into account, in the tight binding approximation, the two kind of bonds \(\epsilon + \Delta\) (double bond) in the unit cell (AB) and \(\epsilon - \Delta\) (single bond) between cells. Two equivalent phases are possible by exchange of single and double bonds: for example a chain with an initial double bound (phase A) and other with an initial single bond (phase B); putting both together introduces a defect, two successive single bonds, for instance. Exchange of the bond pattern corresponds to a sign change of the parameter \(\Delta\), which is also proportional to the energy gap between the two electronic bands \(\Delta\).

Using a momentum representation \(k\), and introducing a spinor \(\psi_k\),

$$ \psi_k = \begin{pmatrix} a_k \\ b_k \end{pmatrix} \,, \quad a_k = \frac{1}{\sqrt{N}}\sum_n \E^{-\I x k} c_{\mathrm{A},x}\,, \quad b_k = \frac{1}{\sqrt{N}}\sum_n \E^{-\I x k} c_{\mathrm{B},x} $$

the Hamiltonian can be written in the diagonal form:

\begin{equation} H = \sum_k \psi^\dagger_k \boldsymbol d(k) \cdot \boldsymbol \sigma \psi_k \,,\quad d_x = (1 + \Delta) + (1 - \Delta) \cos k\,, \quad d_y = (1 - \Delta) \sin k \end{equation}

(taking \(\epsilon=1\), the energy scale) from which we readily find the two eigenvalues

\begin{equation} E_\pm = \pm \sqrt{ d_x^2 + d_y^2} = \pm d \end{equation}

and the eigenvectors,

\begin{equation} u_+ = \frac{1}{\sqrt{2}} \begin{pmatrix} 1 \\ \E^{\I \varphi_k} \end{pmatrix}\,, \quad u_- = \frac{1}{\sqrt{2}} \begin{pmatrix} 1 \\ -\E^{\I \varphi_k} \end{pmatrix} \end{equation}

sshd sshe

Varying the value of \(k\) in the circle \((-\pi,\pi)\), the vector \(\boldsymbol d(k)\) describes a circle; if \(\Delta>0\) the circle do not contain the origin, but if \(\Delta<0\) it will contain the origin. The winding number can be computed from the logarithm of the complex function \(h(k) = d_x(k) - \I d_y(k)= |E_k|\exp[\I \mathrm{Arg}\,h(k)]\):

\begin{equation} C = \frac{\I}{2\pi}\int_{-\pi}^{\pi} \Di{k} \frac{\partial }{\partial k}\log q(k)\,, \quad q = \frac{h}{|h|}\,. \end{equation}

Using the above definitions, one gets

\begin{equation} C = \begin{cases} 0, &\text{ if } \Delta > 0 \\ 1, &\text{ if } \Delta < 0 \end{cases} %> \end{equation}

The split step walk

We consider a walker on a one dimensional lattice of points \(x\), whose motion is determined by the orientation of its spin. A single step decomposes in:

  • a spin rotation \(R(\theta) = \E^{-\I \sigma_y \theta/2}\),
  • and a shift \(S=S_\downarrow + S_\uparrow\) towards the left \(S_\downarrow\) for the spin down, and to the right \(S_\uparrow\) for the spin up,
  • the evolution operator is
    $$U = S [I \otimes R(\theta)] \equiv \E^{-\I H_{\text{eff}}}$$
    where we defined an effective Hamiltonian.

Using the momentum representation

$$S = \int_{-\pi}^{\pi} \Di{k} |k\rangle\langle k| \otimes \E^{\I k \sigma_z}$$

the effective Hamiltonian becomes diagonal

\begin{equation} H_{\text{eff}} = \int_{-\pi}^{\pi} \Di{k} |k\rangle\langle k| \otimes E_k \boldsymbol n(k) \cdot \boldsymbol \sigma \end{equation}

with the energy spectrum

\begin{equation} \cos E_k = \cos (\theta/2) \cos k \end{equation}

and the spin quantization axis \(\boldsymbol n(k)\) of the spinor eigenstates is,

\begin{equation} \boldsymbol n(k) = \frac{1}{\sin E_k} \begin{pmatrix} \sin (\theta/2) \sin k\\ \sin (\theta/2) \cos k\\ -\cos (\theta/2) \sin k \end{pmatrix} \end{equation}
def rotation_1(N, theta):
    q = 0.5*theta*ones(2*N + 1)
    return array([[cos(q), -sin(q)], 
                  [sin(q), cos(q)]])
    return kron(eye(2*N+1), r1)

def rotation_2(N, theta_m, theta_p):
    x = arange(-N, N+1)
    delta = 0.01
    q = 0.5*((theta_p + theta_m)/2 + 
        (theta_p - theta_m)/2*tanh((x + 0.5)/delta))
    return array([[cos(q), -sin(q)], 
                [sin(q), cos(q)]])

def init_psi(N, a, b):
    psi = zeros((2, 2*N + 1))
    psi[0,N] = a
    psi[1,N] = b
    return psi

def mult(m, v):
    return transpose(einsum('ij...,j...', m, v),(1,0))

def qw_split(N, theta, theta_m, theta_p, a, b):
    """Reference: T. Kitagawa (2012) 
        "Topological phenomena in quantum walks",
        Quantum Information Processing 11, p. 1107
    r1 = rotation_1(N, theta)
    r2 = rotation_2(N, theta_m, theta_p)
    psi = init_psi(N, a, b)
    psi_t = zeros((2, 2*N + 1, N+1))
    psi_t[:,:,0] = init_psi(N, a, b)
    for n in arange(1, N+1):
        psi1 = mult(r1, psi) # rotation theta1
        psi2 = psi1
        psi2[0,:] = roll(psi1[0,:], 1) # shift up
        psi3 = mult(r2, psi2) # rotation theta2
        psi = psi3
        psi[1,:] = roll(psi3[1,:], -1) # shift down
        psi_t[:,:,n] = psi
    return psi_t

def measure(psi):
    return abs(psi[0,:])**2 + abs(psi[1,:])**2

free loc

»quantum chaos»kicked rotator»dynamical localization»random matrices