Advances in Condensed Matter Physics

Volume 2015 (2015), Article ID 607021, 10 pages

http://dx.doi.org/10.1155/2015/607021

## Noise from Glauber Dynamics: Self-Consistent Interaction and Temperature Dependent Correlations

Department of Physics and Astronomy, University of California, Riverside, CA 92521, USA

Received 10 December 2014; Accepted 18 December 2014

Academic Editor: Sergei Sergeenkov

Copyright © 2015 Amrit De. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Low frequency flicker noise has been argued to occur in spatially extended metastable systems near a critical point (Bak et al., 1987). An Ising-Glauber model based method is suggested here to systematically obtain temperature dependent *n*th-order correlation functions for arbitrary interacting two-level systems (TLSs). This model is fully consistent with existing methods to calculate noise spectra from TLSs and complements them. However, with as such no *a priori* assumptions on the typical log normal distribution of fluctuation rates, it is shown that noise manifests in two different cases: first in the thermodynamic limit on a 2D lattice with long range antiferromagnetic interactions at low temperatures and second in the case of a statistical ensemble of finite-sized spin clusters representing disorder, but where each cluster is ordered due to ferromagnetic interactions.

#### 1. Introduction

Low frequency flicker noise or noise is seen in all sorts of systems such as in solid-state devices, road traffic, and biological, geological, and financial systems. Despite its widespread occurrence, there is no common physical underlying mechanism that gives rise to all the different manifestations of this noise [1]. However, it has been argued that this flicker noise will occur in dynamical systems with extended spatial degrees of freedom that are barely stable, which evolve into self-organized critical structures [2].

Ever since its first measurement in the 1920s [3], noise has been seen in a wide variety of solid state systems [1, 4, 5]. In the case of spin glasses, noise is related to the magnetic fluctuations [6, 7]. This is directly related to the low temperature kinetics of glassy systems and is often associated with the onset of a spin-glass phase [8].

In many solid state devices the presence of parasitic two-level systems (TLSs), possibly due to the presence of defects, generates random-telegraphic noise (RTN) [9, 10]. Typically, if there are a large number of RTN sources (same as flip-flopping Ising spins or fluctuating TLSs), then a log normal distribution of their switching rates gives rise to a type power spectrum within some range of frequencies.

Understanding and reducing noise are extremely important for improving the performance of various semiconductor devices [11] such as field-effect transistors (FETs) [12, 13], core-shell nanowire FETs [14], GaN/AlGaN heterostructures [15], and graphene based electronics [16, 17]. noise is also seen in metal-insulator composites and tunnel junctions [18–20].

For solid state quantum computing, noise is a major problem as it is a significant source of decoherence [5, 21–24]. In semiconductor quantum dots (QDs), RTN is observed when the potential in the dot lines up with the chemical potential in the reservoir, causing electrons to randomly tunnel back and forth [25–27]. In the case of superconducting flux and phase qubits, experiments have revealed that the magnetic flux noise in these devices have a -type power spectrum [28, 29]. Though this type of magnetic flux noise was first observed in SQUIDs in the 80s [30, 31], its origins were never fully explained. Recent activity in quantum computing has however revived interest in this subject [5, 32–43]. Overall, a better understanding of this noise can tremendously aid in making qubits more robust [10, 44–46].

Various theoretical models have been proposed to explain noise in SQUIDs. Examples include a model proposed by Koch et al. [47], where unpaired and noninteracting electrons randomly hop between traps with different spin orientations and a distribution of trap energies. noise can also be explained using a dangling bond model [48]. Choi et al. [49] recently explained this in terms of metal induced gap states that arise due to the potential disorder at the metal-insulator interface. In general, hopping conductivity models are typically used to explain noise in electron glasses [50–53]. Sometime a model was proposed, where the amplitude of noise scaled as [54]; however, this model did not consider interactions between TLSs.

Interactions are important for understanding noise spectra in different types of systems such as in spin glasses [8], where sometimes the temperature dependence of the noise can even be conflicting [6, 55]. Coulomb interactions and correlated hopping are particularly important in the case of Coulomb glasses [53, 56]. Spatially correlated fluctuations and the temperature dependence are also very important for understanding noise spectrum in metal films [57]. In the case of noise in SQUIDs, experimental evidence suggests that the TLSs are strongly interacting and that there is a net spin polarization [34, 42, 58], which cannot be explained by a spin glass phase [59]. Another example where both interactions and disorder also play an important role is in the case of metal-insulator transitions in 2D systems [60], where noise has been observed [61, 62].

The purpose of this paper is to introduce a general fully self-consistent model that takes the temperature dependence, TLS-TLS interactions, and spatial correlations into account while calculating the noise spectra. This is done within the framework of Ising-Glauber dynamics. A straightforward method is suggested to systematically extract any th correlation function for interacting Ising spins. This method is well suited for numerics as well as for analytics in smaller systems. It also allows open quantum systems methods such as the quasi-Hamiltonian formalism [24, 40, 63–67] to connect to more complicated underlying temperature dependent noise microscopics. A temperature dependent spin-cluster model with ferromagnetic interactions, based on this method, was proposed recently to explain various puzzling features of magnetization noise in SQUIDs [68].

In the case of the interaction dependent model presented here, as such no* a priori* assumptions are needed for a log normal distribution of fluctuation rates in order to obtain noise. For long range antiferromagnetic interactions, there are signatures of type noise on an infinite 2D lattice, without making any heuristic assumption for switching rates, whereas, for long range ferromagnetic interactions, noise is shown to manifest from multiple clusters with a normal distribution of lattice constants. This spin-cluster model represents a disordered system. The overall suggested method is fully self-consistent: in the infinite temperature limit, the model can reduce to usual case with some assumed switching rate distribution.

This paper is organized as follows. In Section 2, the model and the correlation function calculation technique are discussed. An explicit example is given for a pair of interacting TLSs, where analytical expressions are derived for the correlation functions and the noise power spectra. More general results are discussed in Section 3, where this model can be used to show the manifestation of noise in two different cases. This is followed by Section 4.

#### 2. Model for Noise Dynamics

##### 2.1. Uncorrelated Two-Level System Fluctuations

A fluctuating two-level system is an Ising spin whose spin flips have a random time dependence, which is a stochastic process. Its dynamics is therefore governed by the master equation, [69], where is a matrix of transition rates (such that the sum of each of its columns is zero) and is the flipping probability matrix for the TLS. If the average occupation of the two states is the same (for unbiased fluctuators), then for spin- and the corresponding flipping probability matrix is
where is the th spin’s relaxation rate. For a single Ising spin, then the th order autocorrelation function can be obtained from as follows:
where is the same as in (1) but with and and are the initial and final state vectors for the Ising spin that satisfy . For an unbiased spin (*i.e.,* with equal flipping probabilities), and is the Pauli matrix.

For uncorrelated fluctuating Ising spins, the probability matrix is simply and the matrix , which for the case of two fluctuating Ising spins is where is a identity matrix. Alternatively, if one started with the matrix of transition rates, then the explicit flipping probability matrix is explicitly .

For the two uncorrelated Ising spins (or TLS), the two-point autocorrelation function and cross-correlation function can be obtained as follows: where , , and the projection vectors and . As expected, the cross-correlation term is zero as we only considered uncorrelated fluctuating Ising spins. The noise power spectrum can be obtained by taking the Fourier transform of the above correlation functions.

In general for spins (interacting or noninteracting), any th order correlation function between arbitrary spins can be calculated as follows: and it is implied that And and are in the same lexicographically ordered Ising spin basis. The matrices for interacting Ising spins are discussed in the next section.

##### 2.2. Correlated Two-Level System Fluctuations with Ising-Glauber Dynamics

The method described in the previous section can also be used to obtain arbitrary -point correlation functions for correlated spin fluctuations. The system’s overall time development is governed by the master equation as usual. The nonequilibrium spin dynamics for a system of fluctuating spins can be handled using the Ising-Glauber model [70, 71]. Single-site Glauber dynamics (also known as the kinetic Ising model) is a Markov process where the new distribution of spins depends only on the current spin configuration. A single spin is flipped at a given site, and the new configuration agrees with the old one everywhere except at the site where the spin has been flipped. In Glauber dynamics, the conditional probability for a single spin to flip is determined by the Boltzmann factor. For correlated spin fluctuations, the matrix elements of the matrix of transition rates are Here, is a vector denoting the present spin configuration of the lattice, denotes the spin configuration of the lattice at an earlier instance of time, and is the relaxation rate of the spin that is flipped. The nonnegative off-diagonal matrix elements in (10) satisfy the detailed balance condition and the diagonal terms are just the negative sum of the off-diagonal column elements so that all column elements sum up to zero, which ensures the conservation of probability. The systems temporal dynamics is then governed by the flipping probability matrix, which is . The eigenvalues of are either zero, which corresponds to the equilibrium distribution, or real-negative which also eventually tend to the equilibrium distribution as [69].

The Hamiltonian in the Boltzmann factor in (10) is where is the magnetic field and is the spin-spin interaction between the th and th spins.

##### 2.3. Two Correlated Fluctuators

In order to obtain the power spectrum of a given pair of interacting spins, one has to first obtain the two-point correlation functions. Consider two fluctuating Ising spins that are correlated. If , then the matrix in the basis has the following form:where . From this form, it is apparent that if , then (4) is recovered. The flipping probability matrix is then which can be obtained explicitly. The autocorrelation functions and the cross-correlation functions can be calculated using (6) and (7): where and and The projection vectors corresponding to the zero eigenvalue solution of are

In the infinite temperature limit () the results of (6) and (7) are recovered, whereas, in the limit of ,
where . Note that here , which implies that two interacting TLSs can be expressed as a single TLS with effective flipping rate . This could lead to a simple way to* renormalize* and reduce large dimensional matrices for more conducive lattice structures. The effective flipping rates can then be introduced in the form of an effective .

The temperature dependent autocorrelation functions are shown in Figure 1(a) and the cross-correlation functions are shown for ferromagnetic interactions in Figure 1(b) and antiferromagnetic interactions in Figure 1(c). As expected for antiferromagnetic interactions, the cross-correlation function is anticorrelated for a pair of fluctuators at initial times, whereas the autocorrelation function does not depend on the sign of . Also note that the maximum amplitude of the cross-correlation function is inversely proportional to the temperature as expected, thereby ensuring that the cross-correlation terms vanish at high temperatures.