Source author record

Soumya Bera

Soumya Bera appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2025arXiv

Critical States of Fermions with ${\mathbb{Z}}_2$ Flux Disorder

Motivated by many contemporary problems in condensed matter physics where matter particles experience random gauge fields, we investigate the physics of fermions on a square lattice with $π$-flux (that realizes Dirac fermions at low energies), subjected to flux disorder arising from a random ${\mathbb{Z}}_2$ gauge field that results from the presence of flux defects (plaquettes with zero flux). At half-filling, where the system possesses BDI symmetry, we show that a new class of critical phases is realized, with the states at zero energy showing a multifractal character. The multifractal properties depend on the concentration $\mathfrak{c}$ of the $π$-flux defects and spatial correlations between the flux defects. These states are characterized by the singularity spectrum, Lyapunov exponents, and transport properties. For any concentration of flux defects, we find that the multi-fractal spectrum shows termination, but $\textit{not freezing}$. We characterize this class of critical states by uncovering a relation between the conductivity and the Lyapunov exponent, which is satisfied by the states irrespective of the concentration or the local correlations between the flux defects. We demonstrate that renormalization group methods, based on perturbing the Dirac point, fail to capture this new class of critical states. This work not only offers new challenges to theory, but is also likely to be useful in understanding a variety of problems where fermions interact with discrete gauge fields.

preprint2025arXiv

Diffusive metal in a percolating Chern insulator

Two-dimensional non-interacting fermions without any anti-unitary symmetries generically get Anderson localized in the presence of disorder. In contrast, topological superconductors with their inherent particle-hole symmetry can host a thermal metallic phase, which is non-universal and depends on the nature of microscopic disorder. In this work, we demonstrate that in the presence of geometric disorders, such as random bond dilution, a robust metal can emerge in a Chern insulator with particle-hole symmetry. The metallic phase is realized when the broken links are weakly stitched via concomitant insertion of $π$ fluxes in the plaquettes. These nucleate low-energy manifolds, which can provide percolating conduction pathways for fermions to elude localization. This diffusive metal, unlike those in superconductors, can carry charge current and even anomalous Hall current. We investigate the transport properties and show that while the topological insulator to Anderson insulator transition exhibits the expected Dirac universality, the metal insulator transition displays a different critical exponent $ν\approx 2$ compared to a disordered topological superconductor, where $ν\approx 1.4$. Our work emphasizes the unique role of geometric disorder in engineering novel phases and their transitions in topological quantum matter.

preprint2020arXiv

Anomalous diffusion in particle-hole symmetric many-body localized systems

In this work we probe the dynamics of the particle-hole symmetric many-body localized (MBL) phase. We provide numerical evidence that it can be characterized by an algebraic propagation of both entanglement and charge, unlike in the conventional MBL case. We explain the mechanism of this anomalous diffusion through a formation of bound states, which coherently propagate via long-range resonances. By projecting onto the two-particle sector of the particle-hole symmetric model, we show that the formation and observed subdiffusive dynamics is a consequence of an interplay between symmetry and interactions.

preprint2020arXiv

High-temperature spin dynamics in the Heisenberg chain: Magnon propagation and emerging KPZ-scaling in the zero-magnetization limit

The large-scale dynamics of quantum integrable systems is often dominated by ballistic modes due to the existence of stable quasi-particles. We here consider as an archetypical example for such a system the spin-$\frac{1}{2}$ XXX Heisenberg chain that features magnons and their bound states. An interesting question, which we here investigate numerically, arises with respect to the fate of ballistic modes at finite temperatures in the limit of zero magnetization $m{=}0$. At a finite magnetization density $m$, the spin autocorrelation function $Π(x,t)$ (at high temperatures) typically exhibits a trimodal behavior with left- and right-moving quasi-particle modes and a broad center peak with slower dynamics. The broadening of the fastest propagating modes exhibits a sub-diffusive $t^{1/3}$ scaling at large magnetization densities, $m {\rightarrow} \frac{1}{2}$, familiar from non-interacting models; it crosses over into a diffusive scaling $t^{1/2}$ upon decreasing the magnetization to smaller values. The behavior of the center peak appears to exhibit a crossover from transient super-diffusion to ballistic relaxation at long times. In the limit $m{\to}0$, the weight carried by the propagating peaks tends to zero; the residual dynamics is carried only by the central peak; it is sub-ballistic and characterized by a dynamical exponent $z$ close to the value $\frac{3}{2}$ familiar from Kardar-Parisi-Zhang (KPZ) scaling. We confirm, employing elaborate finite-time extrapolations, that the spatial scaling of the correlator $Π$ is in excellent agreement with KPZ-type behavior and analyze the corresponding corrections.

preprint2019arXiv

Apparent slow dynamics in the ergodic phase of a driven many-body localized system without extensive conserved quantities

We numerically study the dynamics on the ergodic side of the many-body localization transition in a periodically driven Floquet model with no global conservation laws. We describe and employ a numerical technique based on the fast Walsh-Hadamard transform that allows us to perform an exact time evolution for large systems and long times. As in models with conserved quantities (e.g., energy and/or particle number) we observe a slowing down of the dynamics as the transition into the many-body localized phase is approached. More specifically, our data is consistent with a subballistic spread of entanglement and a stretched-exponential decay of an autocorrelation function, with their associated exponents reflecting slow dynamics near the transition for a fixed system size. However, with access to larger system sizes, we observe a clear flow of the exponents towards faster dynamics and can not rule out that the slow dynamics is a finite-size effect. Furthermore, we observe examples of non-monotonic dependence of the exponents with time, with dynamics initially slowing down but accelerating again at even larger times, consistent with the slow dynamics being a crossover phenomena with a localized critical point.

preprint2019arXiv

Insights on heterogeneity in blinking mechanisms and non-ergodicity using sub-ensemble statistical analysis of single quantum-dots

Photo-luminescence intermittency (blinking) in semiconductor nanocrystals (NCs), a phenomenon ubiquitous to single-emitters, is generally considered to be temporally random intensity fluctuations between bright (On) and dark (Off) states. However, individual quantum-dots (QDs) rarely exhibit such telegraphic signal, and yet, the vast majority of single-NC blinking data are analyzed using a single fixed threshold, which generates binary trajectories. Further, blinking dynamics can vary dramatically over NCs in the ensemble, and it is unclear whether the exponents (m) of single-particle On-/Off-time distributions (P(t)-On/Off), which are used to validate mechanistic models of blinking, are narrowly distributed or not. Here, we sub-classify an ensemble based on the emissivity of QDs, and subsequently compare the (sub)ensemble behaviors. To achieve this, we analyzed a large number (>1000) of intensity trajectories for a model system, Mn+2 doped ZnCdS QDs, which exhibits diverse blinking dynamics. An intensity histogram dependent thresholding method allowed us to construct distributions of relevant blinking parameters (such as m). Interestingly, we find that single QD P(t)-On/Off s follow either truncated power law or power law, and their relative proportion vary over sub-populations. Our results reveal a remarkable variation in m(On/Off) amongst as well as within sub-ensembles, which implies multiple blinking mechanisms being operational among various QDs. We further show that the m(On/Off) obtained via cumulative single-particle P(t)-On/Off is clearly distinct from the weighted mean value of all single-particle m(On/Off), an evidence for the lack of ergodicity. Thus, investigation and analyses of a large number of QDs, albeit for a limited time-span of few decades, is crucial to characterize possible blinking mechanisms and heterogeneity therein

preprint2019arXiv

Sub-diffusion in the Anderson model on random regular graph

We study the finite-time dynamics of an initially localized wave-packet in the Anderson model on the random regular graph (RRG). Considering the full probability distribution $Π(x,t)$ of a particle to be at some distance $x$ from the initial state at time $t$, we give evidence that $Π(x,t)$ spreads sub-diffusively over a range of disorder strengths, wider than a putative non-ergodic phase. We provide a detailed analysis of the propagation of $Π(x,t)$ in space-time $(x,t)$ domain, identifying four different regimes. These regimes in $(x,t)$ are determined by the position of a wave-front $X_{\text{front}}(t)$, which moves sub-diffusively to the most distant sites $X_{\text{front}}(t) \sim t^β$ with an exponent $β< 1$. We support our numerical results by a self-consistent semiclassical picture of wavepacket propagation relating the exponent $β$ with the relaxation rate of the return probability $Π(0,t) \sim e^{-Γt^β}$. Importantly, the Anderson model on the RRG can be considered as proxy of the many-body localization transition (MBL) on the Fock space of a generic interacting system. In the final discussion, we outline possible implications of our findings for MBL.

preprint2018arXiv

Dynamical quantum phase transitions in Weyl semimetals

The quench dynamics in type-I inversion symmetric Weyl semimetals (WSM) are explored in this work which, due to the form of the Hamiltonian, may be readily extended to two-dimensional Chern insulators. We analyze the role of equilibrium topological properties characterized by the Chern number of the pre-quench ground state in dictating the non-equilibrium dynamics of the system, specifically, the emergence of dynamical quantum phase transitions (DQPT). By investigating the ground state fidelity, it is found that a change in the signed Chern number constitutes a sufficient but not necessary condition for the occurrence of DQPTs. Depending on the ratio of the transverse and longitudinal hopping parameters, DQPTs may also be observed for quenches lying entirely within the initial Chern phase. Additionally, we analyze the zeros of the boundary partition function discovering that while the zeros generally form two-dimensional structures resulting in one-dimensional critical times, infinitesimal quenches may lead to one-dimensional zeros with zero-dimensional critical times provided the quench distance scales appropriately with the system size. This is strikingly manifested in the nature of non-analyticies of the dynamical free energy, revealing a logarithmic singularity. In addition, following recent experimental advances in observing the dynamical Fisher zeros of the Loschmidt overlap amplitude through azimuthal Bloch phase vortices by Bloch-state tomography, we rigorously investigate the same in WSMs. Finally, we establish the relationship between the dimension of the critical times and the presence of dynamical vortices, demonstrating that only one-dimensional critical times arising from two-dimensional manifolds of zeros of the boundary partition function lead to dynamical vortices.

preprint2016arXiv

Dirty Weyl semimetals: Stability, phase transition and quantum criticality

We study the stability of three-dimensional incompressible Weyl semimetals in the presence of random quenched charge impurities. Combining numerical analysis and scaling theory we show that in the presence of sufficiently weak randomness (i) Weyl semimetal remains stable, while (ii) double-Weyl semimetal gives rise to compressible diffusive metal where the mean density of states at zero energy is finite. At stronger disorder, Weyl semimetal undergoes a quantum phase transition and enter into a metallic phase. Mean density of states at zero energy serves as the order parameter and displays single-parameter scaling across such disorder driven quantum phase transition. We numerically determine various exponents at the critical point, which appear to be insensitive to the number of Weyl pairs. We also extract the extent of the quantum critical regime in disordered Weyl semimetal and the phase diagram of dirty double Weyl semimetal at finite energies.

preprint2016arXiv

Dynamics of a Qubit in a High-Impedance Transmission Line from a Bath Perspective

We investigate quantum dynamics of a generic model of light-matter interaction in the context of high impedance waveguides, focusing on the behavior of the emitted photonic states, in the framework of the spin-boson model Quantum quenches as well as scattering of an incident coherent pulse are studied using two complementary methods. First, we develop an approximate ansatz for the electromagnetic waves based on a single multimode coherent state wavefunction; formally, this approach combines ideas from adiabatic renormalization, the Born-Markov approximation, and input-output theory. Second, we present numerically exact results for scattering of a weak intensity pulse by using NRG calculations. NRG provides a benchmark for any linear response property throughout the ultra-strong coupling regime. We find that in a sudden quantum quench, the coherent state approach produces physical artifacts, such as improper relaxation to the steady state. These previously unnoticed problems are related to the simplified form of the ansatz that generates spurious correlations within the bath. In the scattering problem, NRG is used to find the transmission and reflection of a single photon, as well as the inelastic scattering of that single photon. Simple analytical formulas are established and tested against the NRG data that predict quantitatively the transport coefficients for up to moderate environmental impedance. These formulas resolve pending issues regarding the presence of inelastic losses in the spin-boson model near absorption resonances, and could be used for comparison to experiments in Josephson waveguide QED. Finally, the scattering results using the coherent state wavefunction approach are compared favorably to the NRG results for very weak incident intensity. We end our study by presenting results at higher power where the response of the system is nonlinear.

preprint2016arXiv

Generalized Dyson model: nature of zero mode and its implication in dynamics

We study the role of the anomalous $E=0$ state in dynamical properties of non-interacting fermionic chains with chiral symmetry and correlated bond disorder in one dimension. These models posses a diverging density of states at zero energy leading to a divergent localization length at the band center. By analytically calculating the localization length for a finite system, we show that correlations in the disorder modify the spatial decay of the $E = 0$ state from being quasilocalized to extended. We numerically simulate charge and entanglement propagation and provide evidence that states close to $E=0$ dominate the dynamical properties. Remarkably, we find that correlations lead to subdiffusive charge propagation, whereas the growth of entanglement is logarithmically slow. A logarithmic scaling of entanglement saturation with system size is also observed, which indicates a behavior akin to quantum critical glasses.

preprint2016arXiv

Graphene with vacancies: supernumerary zero modes

The density of states, $\varrho(E)$, of graphene is investigated within the tight binding (Hückel) approximation in the presence of vacancies. They induce a non-vanishing density of zero modes, $n_\text{zm}$, that act as midgap states: $\varrho(E)=n_\text{zm}δ(E) + \text{smooth}$. As is well known, the actual number of zero modes per sample can in principle exceed the sublattice imbalance: $N_\text{zm}\geq |N_\text{A}-N_\text{B}|$, where $N_\text{A}$, $N_\text{B}$ denote the number of carbon atoms in each sublattice. In this work, we establish a stronger relation that is valid in the thermodynamic limit and that involves the concentration of zero-modes: $n_\text{zm}>|c_\text{A}-c_\text{B}|$, where $c_\text{A}$ and $c_\text{B}$ denote the concentration of vacancies per sublattice; in particular, $n_\text{zm}$ is non-vanishing even in the case of balanced disorder, $N_\text{A}/N_\text{B}=1$. Adopting terminology from benzoid graph theory, the excess modes associated with the current carrying backbone (percolation cluster) are called supernumerary. In the simplest cases such modes can be associated with structural elements of internal boundaries like dangling bonds. Our result suggest that the continuum limit of bipartite hopping models supports nontrivial "supernumerary" terms that escape the present continuum descriptions.

preprint2015arXiv

Local entanglement structure across a many-body localization transition

Local entanglement between pairs of spins, as measured by concurrence, is investigated in a disordered spin model that displays a transition from an ergodic to a many-body localized phase in excited states. It is shown that the concurrence vanishes in the ergodic phase and becomes nonzero and increases in the many-body localized phase. This happen to be correlated with the transition in the spectral statistics from Wigner to Poissonian distribution. A scaling form is found to exist in the second derivative of the concurrence with the disorder strength. It also displays a critical value for the localization transition that is close to what is known in the literature from other measures. An exponential decay of concurrence with distance between spins is observed in the localized phase. Nearest neighbor spin concurrence in this phase is also found to be strongly correlated with the disorder configuration of onsite fields: nearly similar fields implying larger entanglement.

preprint2015arXiv

Many-body localization characterized from a one-particle perspective

We show that the one-particle density matrix $ρ$ can be used to characterize the interaction-driven many-body localization transition in closed fermionic systems. The natural orbitals (the eigenstates of $ρ$) are localized in the many-body localized phase and spread out when one enters the delocalized phase, while the occupation spectrum (the set of eigenvalues of $ρ$) reveals the distinctive Fock-space structure of the many-body eigenstates, exhibiting a step-like discontinuity in the localized phase. The associated one-particle occupation entropy is small in the localized phase and large in the delocalized phase, with diverging fluctuations at the transition. We analyze the inverse participation ratio of the natural orbitals and find that it is independent of system size in the localized phase.

preprint2014arXiv

A generalized multi-polaron expansion for the spin-boson model: Environmental entanglement and the biased two-state system

We develop a systematic variational coherent state expansion for the many-body ground state of the spin-boson model, in which a quantum two-level system is coupled to a continuum of harmonic oscillators. Energetic constraints at the heart of this technique are rationalized in terms of polarons (displacements of the bath states in agreement with classical expectations) and antipolarons (counter-displacements due to quantum tunneling effects). We present a comprehensive study of the ground state two-level system population and coherence as a function of tunneling amplitude, dissipation strength, and bias (akin to asymmetry of the double well potential defining the two-state system). The entanglement among the different environmental modes is investigated by looking at spectroscopic signatures of the bipartite entanglement entropy between a given environmental mode and all the other modes. We observe a drastic change in behavior of this entropy for increasing dissipation, indicative of the entangled nature of the environmental states. In addition, the entropy spreads over a large energy range at strong dissipation, a testimony to the wide entanglement window characterizing the underlying Kondo state. Finally, comparisons to accurate numerical renormalization group calculations and to the exact Bethe Ansatz solution of the model demonstrate the rapid convergence of our variationally-optimized multi-polaron expansion, suggesting that it should also be a useful tool for dissipative models of greater complexity, as relevant for numerous systems of interest in quantum physics and chemistry.

preprint2014arXiv

Impact of quantized vibrations on the efficiency of interfacial charge separation in photovoltaic devices

We demonstrate that charge separation at donor-acceptor interfaces is a complex process that is controlled by the combined action of Coulomb binding for electron-hole pairs and partial relaxation due to quantized phonons. A joint electron-vibration quantum dynamical study reveals that high energy vibrations sensitively tune the charge transfer probability as a function of time and injection energy, due to polaron formation. These results have bearings for the optimization of energy transfer both in organic and quantum dot photovoltaics, as well as in biological light harvesting complexes.

preprint2014arXiv

Stabilizing Spin Coherence Through Environmental Entanglement in Strongly Dissipative Quantum Systems

The key feature of a quantum spin coupled to a harmonic bath---a model dissipative quantum system---is competition between oscillator potential energy and spin tunneling rate. We show that these opposing tendencies cause environmental entanglement through superpositions of adiabatic and antiadiabatic oscillator states, which then stabilizes the spin coherence against strong dissipation. This insight motivates a fast-converging variational coherent-state expansion for the many-body ground state of the spin-boson model, which we substantiate via numerical quantum tomography.

preprint2013arXiv

Statistics of Conductances and Subleading Corrections to Scaling near the Integer Quantum Hall Plateau Transition

We study the critical behavior near the integer quantum Hall plateau transition by focusing on the multifractal (MF) exponents $X_q$ describing the scaling of the disorder-average moments of the point contact conductance $T$ between two points of the sample, within the Chalker-Coddington network model. Past analytical work has related the exponents $X_q$ to the MF exponents $Δ_q$ of the local density of states (LDOS). To verify this relation, we numerically determine the exponents $X_q$ with high accuracy. We thereby provide, at the same time, independent numerical results for the MF exponents $Δ_q$ for the LDOS. The presence of subleading corrections to scaling makes such determination directly from scaling of the moments of $T$ virtually impossible. We overcome this difficulty by using two recent advances. First, we construct pure scaling operators for the moments of $T$ which have precisely the same leading scaling behavior, but no subleading contributions. Secondly, we take into account corrections to scaling from irrelevant (in the renormalization group sense) scaling fields by employing a numerical technique ("stability map") recently developed by us. We thereby numerically confirm the relation between the two sets of exponents, $X_q$ (point contact conductances) and $Δ_q$ (LDOS), and also determine the leading irrelevant (corrections to scaling) exponent $y$ as well as other subleading exponents. Our results suggest a way to access multifractality in an experimental setting.

preprint2013arXiv

Unveiling environmental entanglement in strongly dissipative qubits

The coupling of a qubit to a macroscopic reservoir plays a fundamental role in understanding the complex transition from the quantum to the classical world. Considering a harmonic environment, we use both intuitive arguments and numerical many-body quantum tomography to study the structure of the complete wavefunction arising in the strong-coupling regime, reached for intense qubit-environment interaction. The resulting strongly-correlated many-body ground state is built from quantum superpositions of adiabatic (polaron-like) and non-adiabatic (antipolaron-like) contributions from the bath of quantum oscillators. The emerging Schrödinger cat environmental wavefunctions can be described quantitatively via simple variational coherent states. In contrast to qubit-environment entanglement, we show that non-classicality and entanglement among the modes in the reservoir are crucial for the stabilization of qubit superpositions in regimes where standard theories predict an effectively classical spin.