Defect production due to quenching through a multicritical point

We study the generation of defects when a quantum spin system is quenched through a multicritical point by changing a parameter of the Hamiltonian as $t/\tau$, where $\tau$ is the characteristic time scale of quenching. We argue that when a quantum system is quenched across a multicritical point, the density of defects ($n$) in the final state is not necessarily given by the Kibble-Zurek scaling form $n \sim 1/\tau^{d \nu/(z \nu +1)}$, where $d$ is the spatial dimension, and $\nu$ and $z$ are respectively the correlation length and dynamical exponent associated with the quantum critical point. We propose a generalized scaling form of the defect density given by $n \sim 1/\tau^{d/(2z_2)}$, where the exponent $z_2$ determines the behavior of the off-diagonal term of the $2 \times 2$ Landau-Zener matrix at the multicritical point. This scaling is valid not only at a multicritical point but also at an ordinary critical point.Comment: 4 pages, 2 figures, updated references and added one figur

Quenching through Dirac and semi-Dirac points in optical Lattices: Kibble-Zurek scaling for anisotropic Quantum-Critical systems

We propose that Kibble-Zurek scaling can be studied in optical lattices by creating geometries that support, Dirac, Semi-Dirac and Quadratic Band Crossings. On a Honeycomb lattice with fermions, as a staggered on-site potential is varied through zero, the system crosses the gapless Dirac points, and we show that the density of defects created scales as $1/\tau$, where $\tau$ is the inverse rate of change of the potential, in agreement with the Kibble-Zurek relation. We generalize the result for a passage through a semi-Dirac point in $d$ dimensions, in which spectrum is linear in $m$ parallel directions and quadratic in rest of the perpendicular $(d-m)$ directions. We find that the defect density is given by $1 /{\tau^{m\nu_{||}z_{||}+(d-m)\nu_{\perp}z_{\perp}}}$ where $\nu_{||}, z_{||}$ and $\nu_{\perp},z_{\perp}$ are the dynamical exponents and the correlation length exponents along the parallel and perpendicular directions, respectively. The scaling relations are also generalized to the case of non-linear quenching

Quenching across quantum critical points: role of topological patterns

We introduce a one-dimensional version of the Kitaev model consisting of spins on a two-legged ladder and characterized by Z_2 invariants on the plaquettes of the ladder. We map the model to a fermionic system and identify the topological sectors associated with different Z_2 patterns in terms of fermion occupation numbers. Within these different sectors, we investigate the effect of a linear quench across a quantum critical point. We study the dominant behavior of the system by employing a Landau-Zener-type analysis of the effective Hamiltonian in the low-energy subspace for which the effective quenching can sometimes be non-linear. We show that the quenching leads to a residual energy which scales as a power of the quenching rate, and that the power depends on the topological sectors and their symmetry properties in a non-trivial way. This behavior is consistent with the general theory of quantum quenching, but with the correlation length exponent \nu being different in different sectors.Comment: 5 pages including 2 figures; this is the published versio

Adiabatic multicritical quantum quenches: Continuously varying exponents depending on the direction of quenching

We study adiabatic quantum quenches across a quantum multicritical point (MCP) using a quenching scheme that enables the system to hit the MCP along different paths. We show that the power-law scaling of the defect density with the rate of driving depends non-trivially on the path, i.e., the exponent varies continuously with the parameter $\alpha$ that defines the path, up to a critical value $\alpha= \alpha_c$; on the other hand for $\alpha \geq \alpha_c$, the scaling exponent saturates to a constant value. We show that dynamically generated and {\it path($\alpha$)-dependent} effective critical exponents associated with the quasicritical points lying close to the MCP (on the ferromagnetic side), where the energy-gap is minimum, lead to this continuously varying exponent. The scaling relations are established using the integrable transverse XY spin chain and generalized to a MCP associated with a $d$-dimensional quantum many-body systems (not reducible to two-level systems) using adiabatic perturbation theory. We also calculate the effective {\it path-dependent} dimensional shift $d_0(\alpha)$ (or the shift in center of the impulse region) that appears in the scaling relation for special paths lying entirely in the paramagnetic phase. Numerically obtained results are in good agreement with analytical predictions.Comment: 5 pages, 4 figure

Quenching Dynamics of a quantum XY spin-1/2 chain in presence of a transverse field

We study the quantum dynamics of a one-dimensional spin-1/2 anisotropic XY model in a transverse field when the transverse field or the anisotropic interaction is quenched at a slow but uniform rate. The two quenching schemes are called transverse and anisotropic quenching respectively. Our emphasis in this paper is on the anisotropic quenching scheme and we compare the results with those of the other scheme. In the process of anisotropic quenching, the system crosses all the quantum critical lines of the phase diagram where the relaxation time diverges. The evolution is non-adiabatic in the time interval when the parameters are close to their critical values, and is adiabatic otherwise. The density of defects produced due to non-adiabatic transitions is calculated by mapping the many-particle system to an equivalent Landau-Zener problem and is generally found to vary as $1/\sqrt{\tau}$, where $\tau$ is the characteristic time scale of quenching, a scenario that supports the Kibble-Zurek mechanism. Interestingly, in the case of anisotropic quenching, there exists an additional non-adiabatic transition, in comparison to the transverse quenching case, with the corresponding probability peaking at an incommensurate value of the wave vector. In the special case in which the system passes through a multi-critical point, the defect density is found to vary as $1/\tau^{1/6}$. The von Neumann entropy of the final state is shown to maximize at a quenching rate around which the ordering of the final state changes from antiferromagnetic to ferromagnetic.Comment: 8 pages, 6 figure

Defect generation in a spin-1/2 transverse XY chain under repeated quenching of the transverse field

We study the quenching dynamics of a one-dimensional spin-1/2 $XY$ model in a transverse field when the transverse field $h(=t/\tau)$ is quenched repeatedly between $-\infty$ and $+\infty$. A single passage from $h \to - \infty$ to $h \to +\infty$ or the other way around is referred to as a half-period of quenching. For an even number of half-periods, the transverse field is brought back to the initial value of $-\infty$; in the case of an odd number of half-periods, the dynamics is stopped at $h \to +\infty$. The density of defects produced due to the non-adiabatic transitions is calculated by mapping the many-particle system to an equivalent Landau-Zener problem and is generally found to vary as $1/\sqrt{\tau}$ for large $\tau$; however, the magnitude is found to depend on the number of half-periods of quenching. For two successive half-periods, the defect density is found to decrease in comparison to a single half-period, suggesting the existence of a corrective mechanism in the reverse path. A similar behavior of the density of defects and the local entropy is observed for repeated quenching. The defect density decays as $1/{\sqrt\tau}$ for large $\tau$ for any number of half-periods, and shows a increase in kink density for small $\tau$ for an even number; the entropy shows qualitatively the same behavior for any number of half-periods. The probability of non-adiabatic transitions and the local entropy saturate to 1/2 and $\ln 2$, respectively, for a large number of repeated quenching.Comment: 5 pages, 3 figure

A random fiber bundle with many discontinuities in the threshold distribution

We study the breakdown of a random fiber bundle model (RFBM) with $n$-discontinuities in the threshold distribution using the global load sharing scheme. In other words, $n+1$ different classes of fibers identified on the basis of their threshold strengths are mixed such that the strengths of the fibers in the $i-th$ class are uniformly distributed between the values $\sigma_{2i-2}$ and $\sigma_{2i-1}$ where $1 \leq i \leq n+1$. Moreover, there is a gap in the threshold distribution between $i-th$ and $i+1-th$ class. We show that although the critical stress depends on the parameter values of the system, the critical exponents are identical to that obtained in the recursive dynamics of a RFBM with a uniform distribution and global load sharing. The avalanche size distribution (ASD), on the other hand, shows a non-universal, non-power law behavior for smaller values of avalanche sizes which becomes prominent only when a critical distribution is approached. We establish that the behavior of the avalanche size distribution for an arbitrary $n$ is qualitatively similar to a RFBM with a single discontinuity in the threshold distribution ($n=1$), especially when the density and the range of threshold values of fibers belonging to strongest ($n+1$)-th class is kept identical in all the cases.Comment: 6 pages, 4 figures, Accepted in Phys. Rev.

Quantum Discord in a spin-1/2 transverse XY Chain Following a Quench

We report a study on the zero-temperature quantum discord as a measure of two-spin correlation of a transverse XY spin chain following a quench across a quantum critical point and investigate the behavior of mutual information, classical correlations and hence of discord in the final state as a function of the rate of quenching. We show that though discord vanishes in the limit of very slow as well as very fast quenching, it exhibits a peak for an intermediate value of the quenching rate. We show that though discord and also the mutual information exhibit a similar behavior with respect to the quenching rate to that of concurrence or negativity following an identical quenching, there are quantitative differences. Our studies indicate that like concurrence, discord also exhibits a power law scaling with the rate of quenching in the limit of slow quenching though it may not be expressible in a closed power law form. We also explore the behavior of discord on quenching linearly across a quantum multicritical point (MCP) and observe a scaling similar to that of the defect density.Comment: 6 pages, 5 figure

Effect of discontinuity in threshold distribution on the critical behaviour of a random fiber bundle

The critical behaviour of a Random Fiber Bundle Model with mixed uniform distribution of threshold strengths and global load sharing rule is studied with a special emphasis on the nature of distribution of avalanches for different parameters of the distribution. The discontinuity in the threshold strength distribution of fibers non-trivially modifies the critical stress as well as puts a restriction on the allowed values of parameters for which the recursive dynamics approach holds good. The discontinuity leads to a non-universal behaviour in the avalanche size distribution for smaller values of avalanche size. We observe that apart from the mean field behaviour for larger avalanches, a new behaviour for smaller avalanche size is observed as a critical threshold distribution is approached. The phenomenological understanding of the above result is provided using the exact analytical result for the avalanche size distribution. Most interestingly,the prominence of non-universal behaviour in avalanche size distribution depends on the system parameters.Comment: 6 pages, 4 figures, text and figures modifie

Adiabatic dynamics of an inhomogeneous quantum phase transition: the case of z > 1 dynamical exponent

We consider an inhomogeneous quantum phase transition across a multicritical point of the XY quantum spin chain. This is an example of a Lifshitz transition with a dynamical exponent z = 2. Just like in the case z = 1 considered in New J. Phys. 12, 055007 (2010) when a critical front propagates much faster than the maximal group velocity of quasiparticles vq, then the transition is effectively homogeneous: density of excitations obeys a generalized Kibble-Zurek mechanism and scales with the sixth root of the transition rate. However, unlike for z = 1, the inhomogeneous transition becomes adiabatic not below vq but a lower threshold velocity v', proportional to inhomogeneity of the transition, where the excitations are suppressed exponentially. Interestingly, the adiabatic threshold v' is nonzero despite vanishing minimal group velocity of low energy quasiparticles. In the adiabatic regime below v' the inhomogeneous transition can be used for efficient adiabatic quantum state preparation in a quantum simulator: the time required for the critical front to sweep across a chain of N spins adiabatically is merely linear in N, while the corresponding time for a homogeneous transition across the multicritical point scales with the sixth power of N. What is more, excitations after the adiabatic inhomogeneous transition, if any, are brushed away by the critical front to the end of the spin chain.Comment: 10 pages, 6 figures, improved version accepted in NJ