Localization of weakly disordered flat band states

Certain tight binding lattices host macroscopically degenerate flat spectral bands. Their origin is rooted in local symmetries of the lattice, with destructive interference leading to the existence of compact localized eigenstates. We study the robustness of this localization to disorder in different classes of flat band lattices in one and two dimensions. Depending on the flat band class, the flat band states can either be robust, preserving their strong localization for weak disorder W, or they are destroyed and acquire large localization lengths ξ that diverge with a variety of unconventional exponents ν, ξ ~ 1 /W ^ν .     

A new criterion for localization in Anderson's model for disordered lattices

A new criterion for localization is obtained based on the Economou—Cohen approach to the problem of localization in Anderson's model for disordered lattices. The new criterion is remarkably successful in all cases where independent checks are available.     

Non-Abelian gauge potentials driven localization transition in quasiperiodic optical lattices

Derived from quantum waves immersed in an Abelian gauge potential, the quasiperiodic Aubry-André-Harper (AAH) model is a simple yet powerful Hamiltonian to study the Anderson localization of ultracold atoms. Here, we investigate the localization properties of ultracold atoms in quasiperiodic optical lattices subject to a non-Abelian gauge potential, which are depicted by non-Abelian AAH models. We identify that the non-Abelian AAH models can bear the self-duality. We analyze the localization of such non-Abelian self-dual optical lattices, revealing a rich phase diagram driven by the non-Abelian gauge potential involved: a transition from a pure delocalization phase, then to coexistence phases, and finally to a pure localization phase. This is in stark contrast to the Abelian counterpart that does not support the coexistence phases. Our results establish the connection between localization and gauge symmetry, and thus comprise a new insight on the fundamental aspects of localization in quasiperiodic systems, from the perspective of non-Abelian gauge potential.     

Quantum subdiffusion with two- and three-body interactions

We study the dynamics of a few-quantum-particle cloud in the presence of two- and three-body interactions in weakly disordered one-dimensional lattices. The interaction is dramatically enhancing the Anderson localization length ξ ₁ of noninteracting particles. We launch compact wave packets and show that few-body interactions lead to transient subdiffusion of wave packets, m ₂ ~ t ^α , α< 1, on length scales beyond ξ ₁. The subdiffusion exponent is independent of the number of particles. Two-body interactions yield α ≈ 0.5 for two and three particles, while three-body interactions decrease it to α ≈ 0.2. The tails of expanding wave packets exhibit exponential localization with a slowly decreasing exponent. We relate our results to subdiffusion in nonlinear random lattices, and to results on restricted diffusion in high-dimensional spaces like e.g. on comb lattices.     

<Emphasis Type="Italic">q</Emphasis>-Breathers in discrete nonlinear Schrödinger arrays with weak disorder

Nonlinearity and disorder are key players in vibrational lattice dynamics, responsible for localization and derealization phenomena. q-Breathers—periodic orbits in nonlinear lattices, exponentially localized in the reciprocal linear mode space—is a fundamental class of nonlinear oscillatory modes, currently found in disorder-free systems. In this paper we generalize the concept of q-breathers to the case of weak disorder, taking the Discrete Nonlinear Schrödinger chain as an example. We show that g-breathers retain exponential localization near the central mode, provided that disorder is sufficiently small. We analyze statistical properties of the instability threshold and uncover its sensitive dependence on a particular realization. Remarkably, the threshold can be intentionally increased or decreased by specifically arranged inhomogeneities. This effect allows us to formulate an approach to controlling the energy flow between the modes. The relevance to other model arrays and experiments with miniature mechanical lattices, light and matter waves propagation in optical potentials is discussed.     

On the criterion of the crystal-liquid phase transition

A localization criterion is proposed for the crystal-liquid phase transition (PT). According to this criterion, the PT begins when the E _d/k _b T ratio reaches a boundary value E _d(s)/k _b T _m such that a solid phase is present above it and a liquid phase is present below it in a phase diagram. Here, E _d is the energy of atom delocalization, k _b is the Boltzmann constant, T is the temperature, and E _d(s) is the delocalization energy for a solid phase at melting point T _m. This criterion is shown to generalize the Lindemann criterion of melting to the case of crystallization and the Löven criterion of crystallization to the case of melting. This localization criterion is found to be applicable for both normally melting substances and substances that melt with a decrease in the specific volume upon the transition into a liquid phase. The relation of the localization criterion to the vacancy and diffusional criteria of the crystal-liquid PT has been studied. The inequality T _N < T _m, where T _N is the temperature of the onset of crystallization, is explained using the localization criterion. The calculated values of the T _N /T _m ratio coincide well with the experimental estimates. The maximum value of T _N /T _m is likely to be most probable in crystals with a bcc structure and a small value of the Grüneisen parameter. The T _N /T _m ratio is analyzed at the points in the PT where no change in the specific volume occurs and an entropy jump is nonzero.     

Percolation on Infinite Graphs and Isoperimetric Inequalities

We consider the Bernoulli bond percolation process (with parameter p) on infinite graphs and we give a general criterion for bounded degree graphs to exhibit a non-trivial percolation threshold based either on a single isoperimetric inequality if the graph has a bi-infinite geodesic, or two isoperimetric inequalities if the graph has not a bi-infinite geodesic. This new criterion extends previous criteria and brings together a large class of amenable graphs (such as regular lattices) and non-amenable graphs (such trees). We also study the finite connectivity in graphs satisfying the new general criterion and show that graphs in this class with a bi-infinite geodesic always have finite connectivity functions with exponential decay when p is sufficiently close to one. On the other hand, we show that there are graphs in the same class with no bi-infinite geodesic for which the finite connectivity decays sub-exponentially (down to polynomially) in the highly supercritical phase even for p arbitrarily close to one.     

Inhomogeneity-induced second-order phase transitions in the Potts model on hierarchical lattices

Thermodynamics of the Potts model with an arbitrary number of states is analyzed for a class of hierarchical lattices of fractal dimension d > 1. In contrast to the case of crystal lattice, it is shown that all phase transitions on lattices of this type are of the second order. Critical exponents are determined, their dependence on structural parameters is examined, and scaling relations between them are established. A structural criterion for change in transition order is discussed for inhomogeneous systems. Application of the results to critical phenomena in phase transitions in dilute crystals and porous media is discussed.     

Maximum coefficient of light extinction in a nonabsorbing dispersive medium

The maximum value of the light extinction coefficient μ, which can be observed in a dispersive medium with a relative refractive index n of the scattering particles, is studied within the framework of a quasi-crystalline approximation for nonabsorbing dispersive media consisting of monodisperse spherical scatterers. A change in the diffraction parameter x of the scattering particles and their volume concentration c _v is accompanied by nonmonotonic variations of the extinction coefficient, and the function μ(x, c _v ) exhibits several maxima. The dimensions and concentrations of particles are determined, for which the extinction coefficient reaches the absolute maximum μ^max. The μ^max value exhibits a monotonic growth with increasing relative refractive index n of the scattering particles. The conditions of validity of the Ioffe-Regel criterion of radiation localization have been studied. It is established that the localization in nonabsorbing dispersive media can be observed only for n ? 2.7. The intervals of x and c _v in which the criterion of radiation localization is satisfied in dispersive media consisting of particles with n = 3.0 and 3.5 are determined.     

Bootstrap Percolation and Kinetically Constrained Models on Hyperbolic Lattices

We study bootstrap percolation (BP) on hyperbolic lattices obtained by regular tilings of the hyperbolic plane. Our work is motivated by the connection between the BP transition and the dynamical transition of kinetically constrained models, which are in turn relevant for the study of glass and jamming transitions. We show that for generic tilings there exists a BP transition at a nontrivial critical density, 0<ρ _c <1. Thus, despite the presence of loops on all length scales in hyperbolic lattices, the behavior is very different from that on Euclidean lattices where the critical density is either zero or one. Furthermore, we show that the transition has a mixed character since it is discontinuous but characterized by a diverging correlation length, similarly to what happens on Bethe lattices and random graphs of constant connectivity.     

DNA分子链电子结构特性研究

在单电子紧束缚近似下，建立了一维无序二元DNA分子链模型，计算了链长为2×10⁴个碱基对的DNA分子链的电子态密度、局域化特性，并探讨了碱基对的不同组分、格点能量无序度对电子局域态的影响.结果表明：由于DNA分子链中格点能量无序及碱基对的不同组分的存在，其电子波函数呈现出局域化的特性，而局域长度作为衡量电子局域化程度的一个尺度，受碱基对的组分及格点能量无序度的影响. 关键词： DNA分子链 电子结构 电子局域态 局域长度     

Numerical study of localization in Anderson model for disordered systems

The localization properties in the Anderson model of two dimensional square lattices are investigated numerically. Pretty large lattices composed of 10⁴ (= 100 × 100) sites are dealt with, and the overall behaviors of the eigenvectors near the band center are directly examined. Fairly sharp transition and the exponentially decaying localized states are visualized.     

Non-interacting dimer kinetics in hypercubic lattices

The exact formulation of the kinetic of dimer in hypercubic lattices is developed in the framework of the kinetic lattice gas model. The so-called local evolution rules are used to obtain the hierarchy of equation of motion for the correlation functions where processes like adsorption and desorption are included. The hierarchy of equations are truncated using a mean field (m, n) closures which allows the analytical treatment of the system. A general expression for non-interacting dimer isotherm and two particle correlation functions are obtained in hypercubic lattices.     

Pairing enhancement in Betts lattices with next nearest neighbor couplings: Exact results

Electron instabilities in the Hubbard model with the next nearest neighbor coupling are calculated by exact diagonalization in finite, two-dimensional Betts cells (lattices). A viable spin and charge coherent pairing, signaled by quantum critical points and a negative charge gap region, is found in 8- and 10-site Betts lattices at small and moderate U regions consistent with our exact results in elementary bipartite geometries [Phys. Rev. B 78 (2008) 075431]. The contour isolines for continuous temperature driven-crossover between the Mott-Hubbard insulating and coherent pairing phases are demonstrated. The criteria for smooth and abrupt phase transitions are found for systematic enhancement of coherent pairing by optimization of the next nearest neighbor coupling parameter.     

Quantization and confinement effects in superconducting nanostructures

The confinement of the flux lines by a lattice of submicron holes (‘antidots’) has been studied in nanostructured superconducting Pb/Ge multilayers. By introducing regular arrays of sufficiently large antidots, multi-quanta vortex lattices have been stabilized. Sharp cusp-like magnetization (M) anomalies, appearing at matching fieldsH_m=mφ₀/Sin superconducting films with the antidot lattices having a unit cell areaS, are successfully explained. These anomalies are, analogues of the well-knownM(H) cusp atH_c1, but for the onset of multi-quanta (m+1)φ₀-vortices penetration at each subsequent matching fieldH_m. It is shown that theM(H) curve between the matching fieldsH_m<H<H_m+1follows a simpleM∝ln(H−H_m) dependence. These experimental observations have revealed an unusual expansion of validity of the London limit in superconductors with lattices of relatively large antidots. The successful high quality fit of theM(H,T) curves convincingly demonstrates that a new type of the critical stateB=const (‘single-terrace critical state’) can be realized in superconductors with the antidot lattices.     

Ising ferromagnet with biquadratic exchange

A spin one Ising system with biquadratic exchange, is investigated, using Green's function technique in random phase approximation (RPA). Transition temperature T_c and <(S^z)²> at T_c, are found to increase with biquadratic exchange parameter α for sc, bcc and fcc lattices. The variation of <(S^z)²> at T_c with α is found to be the same for the above lattices.     

Kauffman Cellular Automata on Quasicrystal Topology

In this paper, we perform numerical simulations to study Kauffman cellular automata (KCA) on quasiperiod lattices. In particular, we investigate phase transition, magnetic entropy, and propagation speed of the damage on these lattices. Both the critical threshold parameter \(p_{c}\) and the critical exponents are estimated with good precision. In order to investigate the increase of statistical fluctuations and the onset of chaos in the critical region of the model, we have also defined a magnetic entropy to these systems. It is seen that the magnetic entropy behaves in a different way when one passes from the frozen regime (p < p_c) to the chaotic regime (p > p_c). For a further analysis, the robustness of the propagation of failures is checked by introducing a quenched site dilution probability q on the lattices. It is seen that the damage spreading is quite sensitive when a small fraction of the lattice sites are disconnected. A finite-size scaling analysis is employed to estimate the critical exponents. From these numerical estimates, we claim that on both pure (q =?0) and diluted (q =?0.05) quasiperiodic lattices, the KCA model belongs to the same universality class than on square lattices. Furthermore, with the aim of comparing the dynamical behavior between periodic and quasiperiodic systems, the propagation speed of the damage is also calculated for the square lattice assuming the same conditions. It is found that on square lattices the propagation speed of the damage obeys a power law as \(v\sim (p-p_{c})^{\alpha }\), whereas on quasiperiod lattices, it follows a logarithmic law as \(v \sim \ln (p-p_{c})^{\alpha }\).     

Coupled matter-wave solitons in optical lattices

We make use of a potential model to study the dynamics of two coupled matter-wave or Bose-Einstein condensate (BEC) solitons loaded in optical lattices. With separate attention to linear and nonlinear lattices we find some remarkable differences for response of the system to effects of these lattices. As opposed to the case of linear optical lattice (LOL), the nonlinear lattice (NOL) can be used to control the mutual interaction between the two solitons. For a given lattice wave number k, the effective potentials in which the two solitons move are such that the well (V_eff(NOL)), resulting from the juxtaposition of soliton interaction and nonlinear lattice potential, is deeper than the corresponding well V_eff(LOL). But these effective potentials have opposite k dependence in the sense that the depth of V_eff(LOL) increases as k increases and that of V_eff(NOL) decreases for higher k values. We verify that the effectiveness of optical lattices to regulate the motion of the coupled solitons depends sensitively on the initial locations of the motionless solitons as well as values of the lattice wave number. For both LOL and NOL the two solitons meet each other due to mutual interaction if their initial locations are taken within the potential wells with the difference that the solitons in the NOL approach each other rather rapidly and take roughly half the time to meet as compared with the time needed for such coalescence in the LOL. In the NOL, the soliton profiles can move freely and respond to the lattice periodicity when the separation between their initial locations are as twice as that needed for a similar free movement in the LOL. We observe that, in both cases, slow tuning of the optical lattices by varying k with respect to a time parameter τ drags the oscillatory solitons apart to take them to different locations. In our potential model the oscillatory solitons appear to propagate undistorted. But a fully numerical calculation indicates that during evolution they exhibit decay and revival.     

On stochastic spatial patterns and neuronal polarity

To investigate the statistical behavior in the sizes of finite clusters for percolation, cluster size distribution n _s (p) for site and bond percolations at different lattices and dimensions was simulated using a modified algorithm. An equation to approximate the finite cluster size distribution n _s (p) was obtained and expressed as: log?(n _s (p)) = as ? b log?s + c. Based on the analysis of simulation data, we found that the equation is valid for p from 0 to 1 on site and for the bond percolation of two-dimensional (2D) and three-dimensional (3D) lattices. Furthermore, the relationship between the coefficients of the equation and the occupied ratio p was studied using the finite-size scaling method. When \(x = D(p - p_c )L^{y_t }\) , p < p _c , and D was a nonuniversal metric factor. a was found to be related only to p, and the a-x curves of different lattices were nearly overlapped; b was related to the dimensions and p, and the scaled data of the b of all lattices with the same dimension tended to fall on the same curves. Unlike a and b, c apparently had a quadratic relation with x in 2D lattices and linear relation with x in 3D lattices. The results of this paper could significantly reduce the amount of tasks required to obtain numerical data of on the cluster size distribution for p from 0 to p _c .