BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Rajiv R. P. Singh

Rajiv R. P. Singh appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.str-elcond-mat.stat-mechquant-phhep-thBiological PhysicsBiomoleculescond-matcond-mat.mes-hallcond-mat.mtrl-scicond-mat.quant-gascond-mat.supr-conNeurons and CognitionSubcellular Processes

Catalog footprint

What is connected

40works
13topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

40 published item(s)

preprint2022arXiv

Finite Temperature Strong Coupling Expansions for the SU(N) Hubbard Model

We develop finite temperature strong coupling expansions for the SU(N) Hubbard Model in powers of $βt$, $w=\exp{(-βU)}$ and ${1\over βU}$ for arbitrary filling. The expansions are done in the grand canonical ensemble and are most useful at a density of one particle per site, where for $U$ larger than or of order the Bandwidth, the expansions converge over a wide temperature range $t^2/U \ \lesssim \ T \ \lesssim \ 10 U$. By taking the limit $w\to 0$, valid at temperatures much less than $U$, the expansions turn into a high temperature expansion for a dressed SU(N) Heisenberg model that includes nearest-neighbor exchange, further neighbor exchanges and ring exchanges known from the $T=0$ perturbation theory of the SU(2) Hubbard model. Below a filling of one particle per site, the $w\to 0$ limit corresponds to an effective $t-J$ model. The onset of strong correlations can be identified by a plateau-like behavior in the entropy as a function of temperature. At small deviations from one particle per site, the expansions can be arranged in powers of a small parameter $δ=1-n$, the deviation from one particle per site, where the leading $βt$ dependent terms correspond to holes sloshing around in a disordered SU(N) background. We use these expansions to calculate the thermodynamic properties of the model at moderate and high temperatures over a wide parameter range.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Instabilities of spin-1 Kitaev spin liquid phase in presence of single-ion anisotropies

We study the spin-one Kitaev model on the honeycomb lattice in the presence of single-ion anisotropies. We consider two types of single ion anisotropies: A $D_{111}$ anisotropy which preserves the symmetry between $X$, $Y$, and $Z$ bonds but violates flux conservation and a $D_{100}$ anisotropy that breaks the symmetry between $X$, $Y$, and $Z$ bonds but preserves flux conservation. We use series expansion methods, degenerate perturbation theory, and exact diagonalization to study these systems. Large positive $D_{111}$ anisotropy leads to a simple product ground state with conventional magnon-like excitations, while large negative $D_{111}$ leads to a broken symmetry and degenerate ground states. For both signs there is a phase transition at a small $|D_{111}| \approx 0.12$ separating the more conventional phases from the Kitaev spin liquid phase. With large $D_{100}$ anisotropy, the ground state is a simple product state, but the model lacks conventional dispersive excitations due to the large number of conservation laws. Large negative $D_{100}$ leads to decoupled one-dimensional systems and many degenerate ground states. No evidence of a phase transition is seen in our numerical studies at any finite $D_{100}$. Convergence of the series expansion extrapolations all the way to $D_{100}=0$ suggests that the nontrivial Kitaev spin liquid is a singular limit of this type of single-ion anisotropy going to zero, which also restores symmetry between the $X$, $Y$, and $Z$ bonds.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

On the Divergence of the Ferromagnetic Susceptibility in the SU(N) Nagaoka-Thouless Ferromagnet

Using finite temperature strong coupling expansions for the SU(N) Hubbard Model, we calculate the thermodynamic properties of the model in the infinite-$U$ limit for arbitrary density $0\leq ρ\leq 1$ and all $N$. We express the ferromagnetic susceptibility of the model as a Curie term plus a $Δχ$, an excess susceptibility above the Curie-behavior. We show that, on a bipartite lattice, graph by graph the contributions to $Δχ$ are non-negative in the limit that the hole density $δ=1-ρ$ goes to zero. By summing the contributions from all graphs consisting of closed loops we find that the low hole-density ferromagnetic susceptibility diverges exponentially as $\exp{Δ/T}$ as $T \to 0$ in two and higher dimensions. This demonstrates that Nagaoka-Thouless ferromagnetic state exists as a thermodynamic state of matter at low enough density of holes and sufficiently low temperatures. The constant $Δ$ scales with the SU(N) parameter $N$ as $1/N$ implying that ferromagnetism is gradually weakened with increasing $N$ as the characteristic temperature scale for ferromagnetic order goes down.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Onset of charge incompressibility and Mott gaps in the Honeycomb-Lattice SU(4) Hubbard Model: Lessons for Twisted Bilayer Graphene systems

We use finite temperature strong coupling expansions to calculate thermodynamic properties of the Honeycomb-lattice SU(4) Hubbard model. We present numerical results for various properties including chemical potential, compressibility, entropy and specific heat as a function of temperature and density at several $U/t$ values. We study the onset of charge incompressibility and Mott gaps as the temperature is lowered at integer densities. In the incompressible Mott regime, the expansions are recast into a high temperature expansion for a generalized spin model with SU(4) symmetry, which is then used to study the convergence of strong coupling expansions in t/U. We discuss lessons that can be drawn from high temperature properties of a simple Hubbard model regarding Twisted Bilayer Graphene (TBG) and other magic-angle flat-band systems.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Order, Disorder and Monopole Confinement in the Spin-$1/2$ XXZ Model on a Pyrochlore Tube

We study the ground state and thermodynamic properties of the spin-half XXZ model, with an Ising interaction $J_z$ and a transverse exchange interaction $J_{x}$, on a pyrochlore tube obtained by joining together elementary cubes in a one-dimensional array. Periodic boundary conditions in the transverse directions ensure that the bulk of the system consists of corner-sharing tetrahedra, with the same local geometry as the pyrochlore lattice. We use exact diagonalization, the density matrix renormalization group (DMRG), and minimally entangled typical thermal states (METTS) methods to study the system. When $J_z$ is antiferromagnetic ($J_{z}>0$) and $J_x$ is ferromagnetic ($J_{x}<0$), we find a transition from a spin liquid to an XY ferromagnet, which has power-law correlations at $T=0$. For $J_{z}<0$ and $J_{x}>0$, spin-two excitations are found to have lower energy than spin-one at the transition away from the fully polarized state, showing evidence for incipient spin-nematic order. When both interactions are antiferromagnetic, we find a non-degenerate ground state with no broken symmetries and a robust energy gap. The low energy spectra evolve smoothly from predominantly Ising to predominantly XY interactions. In the spin-liquid regime of small $|J_{x}|$, we study the confinement of monopole-anti-monopole pairs and find that the confinement length scale is larger for $J_x<0$ than for $J_{x}>0$, although both length scales are very short. These results are consistent with a local spin-liquid phase for the Heisenberg antiferromagnet with no broken symmetries.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Magnetism of competing high-spin/low-spin states in Ba$_2$NiO$_2$(AgSe)$_2$ and related two-orbital two-electron systems

We discuss the thermodynamic and magnetic properties of a competing high-spin/low-spin two-orbital two-electron model on a square-lattice possibly relevant to the nickelates such as Ba$_2$NiO$_2$(AgSe)$_2$ (BNOAS). We focus on parameter regimes where a high-spin ($S=1$) and a low-spin ($S=0$) state are energetically close to each other and discuss various exchange processes in such a system. The model we study is a variant of, but different from the Kondo-necklace model proposed for the system by Jin et al (H. S. Jin, W. E. Pickett and K. W. Lee, Phys. Rev. Research 2, 033197 (2020)). Although there are similarities between the two models in terms of the ground-state phases and their symmetries, the detailed properties of the different phases and the phase transitions are entirely different and should be easy to distinguish from an experimental point of view.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Thermodynamic behavior of modified integer-spin Kitaev models on the honeycomb lattice

We study the thermodynamic properties of modified spin-$S$ Kitaev models introduced by Baskaran, Sen and Shankar (Phys. Rev. B 78, 115116 (2008)). These models have the property that for half-odd-integer spins their eigenstates map on to those of spin-1/2 Kitaev models, with well-known highly entangled quantum spin-liquid states and Majorana fermions. For integer spins, the Hamiltonian is made out of commuting local operators. Thus, the eigenstates can be chosen to be completely unentangled between different sites, though with a significant degeneracy for each eigenstate. For half-odd-integer spins, the thermodynamic properties can be related to the spin-1/2 Kitaev models apart from an additional degeneracy. Hence we focus here on the case of integer spins. We use transfer matrix methods, high temperature expansions and Monte Carlo simulations to study the thermodynamic properties of ferromagnetic and antiferromagnetic models with spin $S=1$ and $S=2$. Apart from large residual entropies, which all the models have, we find that they can have a variety of different behaviors. Transfer matrix calculations show that for the different models, the correlation lengths can be finite as $T\to 0$, become critical as $T\to 0$ or diverge exponentially as $T\to 0$. There is a conserved $Z_2$ flux variable associated with each hexagonal plaquette which saturates at the value $+1$ as $T\rightarrow0$ in all models except the $S=1$ antiferromagnet where the mean flux remains zero as $T\to 0$. We provide qualitative explanations for these results.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

High-Field Expansion Approach to Kagome Antiferromagnets with Dzyaloshinskii-Moriya Interactions

We apply linked cluster expansion techniques to study the polarized high-field phase of a spin-half antiferromagnet on the Kagome lattice with Heisenberg and Dzyaloshinskii-Moriya interactions (DMI). We find that the Dirac points of the single-magnon spectrum without DMI are robust against arbitrary DMI when the magnetic field lies in the Kagome plane. Unlike the typical case where DMI gaps the spectrum, here we find that varying the DMI merely shifts the location of the Dirac points. In contrast, a magnetic field with a component out of the Kagome plane gaps the spectrum, leading to topological magnon bands. We map out a topological phase diagram as the couplings are varied by computing the band Chern numbers. A pair of phase transitions are observed and we find an enhanced thermal Hall conductivity near the phase boundary.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On two phases inside the Bose condensation dome of Yb$_2$Si$_2$O$_7$

Recent experimental data on Bose-Einstein Condensation (BEC) of magnons in the spin-gap compound Yb$_2$Si$_2$O$_7$ revealed an asymmetric BEC dome arXiv:1810.13096v2. We examine modifications to the Heisenberg model on a breathing honeycomb lattice, showing that this physics can be explained by competing forms of weak anisotropy. We employ a gamut of analytical and numerical techniques to show that the anisotropy yields a field driven phase transition from a state with broken Ising symmetry to a phase which breaks no symmetries and crosses over to the polarized limit.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Confinement, reduced entanglement, and spin-glass order in a random quantum spin-ice model

We study an effective spin model derived perturbatively from random transverse-field Ising model on the pyrochlore lattice. The model consists of spin-configurations on the pyrochlore lattice, restricted to the spin-ice subspace, with spins interacting with random Ising exchange couplings as well as ring exchanges along the hexagons of the lattice. This model is studied by exact diagonalization upto N=64 site systems. We calculate spin-glass correlation functions and local entanglement entropy $S_T$ between spins in a single tetrahedron and the rest of the system. We find that the model undergoes two phase transitions. At weak randomness the model is in a quantum spin-ice phase where $S_T=\ln{6}$. Increasing randomness first leads to a frozen phase, with long-range spin-glass order and $S_T=\ln{2}$ corresponding to the Cat states associated with Ising order. Further increase in randomness leads to a random resonating-hexagon phase with a frozen backbone of spins and a broad distribution of entanglement entropies. The implications of these studies for non-Kramers rare-earth pyrochlores are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Nuclear relaxation rates in the Herbertsmithite Kagome antiferromagnets ZnCu3(OH)6Cl2

Local spectral functions and Nuclear Magnetic Relaxation (NMR) rates, 1/T1, for the spin-half Heisenberg antiferromagnet on the Kagome Lattice are calculated using Moriyas Gaussian approximation, as well as through an extrapolation of multiple frequency moments. The temperature dependence of the calculated rates is compared with the oxygen 1/T1 NMR data in Herbertsmithite. We find that the Gaussian approximation for 1/T1 shows behavior qualitatively similar to experiments with a sharp drop in rates at low temperatures, consistent with a spin-gapped behavior. However, this approximation significantly underestimates the magnitude of 1/T1 even at room temperature. Rates obtained from extrapolation of multiple frequency moments give very good agreement with the room temperature NMR data with J = 195 +/- 20 K and hyperfine couplings determined independently from other measurements. The use of multiple frequency moments also leads to additional low frequency weight in the local structure factors. The convergence of our calculations with higher frequency moments breaks down at low temperatures suggesting the existence of longer range dynamic correlations in the system despite the very short-range static correlations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Spin Generation Via Bulk Spin Current in Three Dimensional Topological Insulators

To date, spin generation in three-dimensional topological insulators is primarily modeled as a single-surface phenomenon, attributed to the momentum-spin locking on each individual surface. In this article we propose a mechanism of spin generation where the role of the insulating yet topologically non-trivial bulk becomes explicit: an external electric field creates a transverse pure spin current through the bulk of a three-dimensional topological insulator, which transports spins between the top and bottom surfaces. Under sufficiently high surface disorder, the spin relaxation time can be extended via the Dyakonov-Perel mechanism. Consequently both the spin generation efficiency and surface conductivity are largely enhanced. Numerical simulation confirms that this spin generation mechanism originates from the unique topological connection of the top and bottom surfaces and is absent in other two dimensional systems such as graphene, even though they possess a similar Dirac cone-type dispersion.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Spiral versus modulated collinear phases in the quantum ANNNH model

Motivated by the discovery of spiral and modulated collinear phases in several magnetic materials, we investigate the magnetic properties of Heisenberg spin S=1/2 antiferromagnets in 2 and 3 dimensions, with frustration arising from 2nd-neighbor couplings in one axial direction (the ANNNH model). Our results clearly demonstrate the presence of an incommensurate spiral phase at T=0 in 2 dimensions, extending to finite temperatures in 3 dimensions. The crossover between Néel and spiral order occurs at a value of the frustration parameter considerably above the classical value 0.25, a sign of substantial quantum fluctuations. We also investigate a possible modulated collinear phase with a wavelength of 4 lattice spacings, and find that it has substantially higher energy and hence is not realized in the model.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Competeing orders in spin-1 and spin-3/2 XXZ Kagome antiferromagnets: A series expansion study

We study the competition between $\sqrt{3} \times \sqrt{3}$ (RT3) and $q=0$ (Q0) magnetic orders in spin-one and spin-$3/2$ Kagome-lattice XXZ antiferromagnets with varying XY anisotropy parameter $Δ$, using series expansion methods. The Hamiltonian is split into two parts: an $H_0$ which favors the classical order in the desired pattern and an $H_1$, which is treated in perturbation theory by a series expansion. We find that the ground state energy series for the RT3 and Q0 phases are identical up to sixth order in the expansion, but ultimately a selection occurs, which depends on spin and the anisotropy $Δ$. Results for ground state energy and the magnetization are presented. These results are compared with recent spin-wave theory and coupled-cluster calculations. The series results for the phase diagram are close to the predictions of spin-wave theory. For the spin-one model at the Heisenberg point ($Δ=1$), our results are consistent with a vanishing order parameter, that is an absence of a magnetically ordered phase. We also develop series expansions for the ground state energy of the spin-one Heisenberg model in the trimerized phase. We find that the ground state energy in this phase is lower than those of magnetically ordered ones, supporting the existence of a spontaneously trimerized phase in this model.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Early breakdown of area-law entanglement at many-body delocalization transition

We introduce the numerical linked cluster (NLC) expansion as a controlled numerical tool for the study of the many-body localization (MBL) transition in a disordered system with continuous non-perturbative disorder. Our approach works directly in the thermodynamic limit, in any spatial dimension, and does not rely on any finite size scaling procedure. We study the onset of many-body delocalization through the breakdown of area-law entanglement in a generic many-body eigenstate. By looking for initial signs of an instability of the localized phase, we obtain a value for the critical disorder, which we believe should be a lower bound for the true value, that is higher than current best estimates from finite size studies. This implies that most current methods tend to overestimate the extent of the localized phase due to finite size effects making the localized phase appear stable at small length scales. We also study the mobility edge in these systems as a function of energy density, and find that our conclusion is the same at all examined energies.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Finite-temperature superconducting correlations of the Hubbard model

We utilize numerical linked-cluster expansions (NLCEs) and the determinantal quantum Monte Carlo algorithm to study pairing correlations in the square lattice Hubbard model. To benchmark the NLCE, we first locate the finite-temperature phase transition of the attractive model to a superconducting state away from half filling. We then explore the superconducting properties of the repulsive model for the d-wave and extended s-wave pairing symmetries. The pairing structure factor shows a strong tendency to d-wave pairing and peaks at an interaction strength comparable to the bandwidth. The extended s-wave structure factor and correlation length are larger at higher temperatures but clearly saturate as temperature is lowered, whereas the d-wave counterparts, which start off lower at high temperatures, continue to rise near half filling. This rise is even more dramatic in the d-wave susceptibility. The convergence of NLCEs breaks down as the susceptibilities and correlation lengths become large, so we are unable to determine the onset of long-range order. However, our results extend the conclusion, previously restricted to only magnetic and charge correlations, that NLCEs offer unique window into pairing in the Hubbard model at strong coupling.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Nonzero temperature Entanglement Negativity of quantum spin models: Area law, Linked Cluster Expansions and Sudden Death

We show that the bipartite logarithmic entanglement negativity (EN) of quantum spin models obeys an area law at all nonzero temperatures. We develop numerical linked cluster (NLC) expansions for the `area-law' logarithmic entanglement negativity as a function of temperature and other parameters. For one-dimensional models the results of NLC are compared with exact diagonalization on finite systems and are found to agree very well. The NLC results are also obtained for two dimensional XXZ and transverse-field Ising models. In all cases, we find a sudden onset (or sudden death) of negativity at a finite temperature above which the negativity is zero. We use perturbation theory to develop a physical picture for this sudden onset (or sudden death). The onset of EN or its magnitude are insensitive to classical finite-temperature phase transitions, supporting the argument for absence of any role of quantum mechanics at such transitions. On approach to a quantum critical point at $T=0$, negativity shows critical scaling in size and temperature.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The many faces of quantum kagome materials: Interplay of further-neighbour exchange and Dzyaloshinskii-Moriya interaction

The field of frustrated magnetism has been enriched significantly by the discovery of various kagome lattice compounds. These materials exhibit a great variety of macroscopic behaviours ranging from magnetic orders to quantum spin liquids. Using large-scale exact diagonalization, we construct the phase diagram of the $S=1/2$ $J_1$-$J_2$ kagome Heisenberg model with $z$-axis Dzyaloshinskii-Moriya interaction $D_z$. We show that this model can systematically account for many of the experimentally observed phases. Small $J_2$ and $D_z$ can stabilize respectively a gapped and a gapless spin liquid. When $J_2$ or $D_z$ is substantial, the ground state develops a $\mathbf{Q}=0$, $120^\textrm{o}$ antiferromagnetic order. The critical strengths for inducing magnetic transition are $D^c_z \sim 0.1\, J_1$ at $J_2=0$, and $J^c_2\sim 0.4\, J_1$ at $D_z=0$. The previously reported values of $D_z$ and $J_2$ for herbertsmithite [ZnCu$_3$(OH)$_6$Cl$_2$] place the compound in close proximity to a quantum critical point.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Unusual Corrections to Scaling and Convergence of Universal Renyi Properties at Quantum Critical Points

At a quantum critical point, bipartite entanglement entropies have universal quantities which are subleading to the ubiquitous area law. For Renyi entropies, these terms are known to be similar to the von Neumann entropy, while being much more amenable to numerical and even experimental measurement. We show here that when calculating universal properties of Renyi entropies, it is important to account for unusual corrections to scaling that arise from relevant local operators present at the conical singularity in the multi-sheeted Riemann surface. These corrections grow in importance with increasing Renyi index. We present studies of Renyi correlation functions in the 1+1 transverse-field Ising model (TFIM) using conformal field theory, mapping to free fermions, and series expansions, and the logarithmic entropy singularity at a corner in 2+1 for both free bosonic field theory and the TFIM, using numerical linked cluster expansions. In all numerical studies, accurate results are only obtained when unusual corrections to scaling are taken into account. In the worst case, an analysis ignoring these corrections can get qualitatively incorrect answers, such as predicting a decrease in critical exponents with the Renyi index, when they are actually increasing. We discuss a two-step extrapolation procedure that can be used to account for the unusual corrections to scaling.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Corner contribution to the entanglement entropy of an O(3) quantum critical point in 2+1 dimensions

The entanglement entropy for a quantum critical system across a boundary with a corner exhibits a sub-leading logarithmic scaling term with a scale-invariant coefficient. Using a Numerical Linked Cluster Expansion, we calculate this universal quantity for a square-lattice bilayer Heisenberg model at its quantum critical point. We find, for this 2+1 dimensional O(3) universality class, that it is thrice the value calculated previously for the Ising universality class. This relation gives substantial evidence that this coefficient provides a measure of the number of degrees of freedom of the theory, analogous to the central charge in a 1+1 dimensional conformal field theory.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Entanglement properties of the antiferromagnetic-singlet transition in the Hubbard model on bilayer square lattices

We calculate the bipartite \Renyi entanglement entropy of an $L\times L\times 2$ bilayer Hubbard model using a determinantal quantum Monte Carlo method recently proposed by Grover [Phys. Rev. Lett. {\bf 111}, 130402 (2013)]. Two types of bipartition are studied: (i) One that divides the lattice into two $L \times L$ planes, and (ii) One that divides the lattice into two equal-size ($L\times L/2\times 2$) bilayers. We compare our calculations with those for the tight-binding model studied by the correlation matrix method. As expected, the entropy for bipartition (i) scales as $L^2$, while the latter scales with $L$ with possible logarithmic corrections. The onset of the antiferromagnet to singlet transition shows up by a saturation of the former to a maximal value and the latter to a small value in the singlet phase. We comment on the large uncertainties in the numerical results with increasing $U$, which would have to be overcome before the critical behavior and logarithmic corrections can be quantified.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Ground state and spectral properties across a Charge Density Wave transition in a triangular-lattice spinless Fermion model

We study ground state properties and particle excitation spectra across a commensurate charge density wave transition in a system of strongly interacting fermions, using series expansion methods and mean-field theory. We consider a 1/3-filled system of spinless fermions on a triangular-lattice, with hopping parameter $t$, nearest-neighbor repulsion $V$, and a sublattice dependent chemical potential $μ_s$. The phase transition is found to be first order for $μ_s=0$, but becomes continuous with increasing $μ_s$. The particle and hole excitation spectra exhibit dramatic changes in the vicinity of the phase transitions and in the charge-density wave ordered phase. We discuss the relevance of this study to the Pinball Fermi liquid phase postulated theoretically in earlier studies as well as to various strongly correlated triangular-lattice materials.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Magnetic correlations and pairing in the 1/5-depleted square lattice Hubbard model

We study the single-orbital Hubbard model on the 1/5-depleted square lattice geometry, which arises in such diverse systems as the spin-gap magnetic insulator CaV4O9 and ordered-vacancy iron selenides, presenting new issues regarding the origin of both magnetic ordering and superconductivity in these materials. We find a rich phase diagram that includes a plaquette singlet phase, a dimer singlet phase, a Neel and a block-spin antiferromagnetic phase, and stripe phases. Quantum Monte Carlo simulations show that the dominant pairing correlations at half filling change character from d-wave in the plaquette phase to extended s-wave upon transition to the Neel phase. These findings have intriguing connections to iron-based superconductors, and suggest that some physics of multiorbital systems can be captured by a single-orbital model at different dopings.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the entanglement across a cubic interface in 3+1 dimensions

We calculate the area, edge and corner Renyi entanglement entropies in the ground state of the transverse-field Ising model, on a simple-cubic lattice, by high-field and low-field series expansions. We find that while the area term is positive and the line term is negative as required by strong subadditivity, the corner contributions are positive in 3-dimensions. Analysis of the series suggests that the expansions converge up to the physical critical point from both sides. The leading area-law Renyi entropies match nicely from the high and low field expansions at the critical point, forming a sharp cusp there. We calculate the coefficients of the logarithmic divergence associated with the corner entropy and compare them with conformal field theory results with smooth interfaces and find a striking correspondence.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Phase diagram of the frustrated quantum-XY model on the honeycomb lattice studied by series expansions: Evidence for proximity to a bicritical point

We study the nearest-neighbor exchange ($J_1$) and second-neighbor exchange ($J_2$) XY antiferromagnet on the honeycomb lattice using ground state series expansions around N{é}el, columnar and dimer phases. The conventional two-sublattice XY N{é}el order at small $J_2$ vanishes at $J_2/J_1=0.22\pm 0.01$ in agreement with results from Density Matrix Renormalization Group (DMRG) studies. Near the transition, we find evidence for an {\it approximate} emergent symmetry between XY and Ising degrees of freedom, namely the nearest-neighbor Ising and XY spin correlations become nearly equal. This suggests that the system is close to a bicritical point separating XY and Ising orders. At still larger $J_2/J_1$ the columnar and dimer energies are found to be nearly degenerate. At even larger $J_2$ the columnar phase is obtained. The ground state energies in all three phases are in good agreement with the values found in the DMRG studies.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Quantum Critical Universality and Singular Corner Entanglement Entropy of Bilayer Heisenberg-Ising model

We consider a bilayer quantum spin model with anisotropic intra-layer exchange couplings. By varying the anisotropy, the quantum critical phenomena changes from XY to Heisenberg to Ising universality class, with two, three and one modes respectively becoming gapless simultaneously. We use series expansion methods to calculate the second and third Renyi entanglement entropies when the system is bipartitioned into two parts. Leading area-law terms and subleading entropies associated with corners are separately calculated. We find clear evidence that the logarithmic singularity associated with the corners is universal in each class. Its coefficient along the Ising critical line is in excellent agreement with those obtained previously for the transverse-field Ising model. Our results provide strong evidence for the idea that the universal terms in the entanglement entropy provide an approximate measure of the low energy degrees of freedom in the system.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Simulated Cytoskeletal Collapse via Tau Degradation

We present a coarse-grained two dimensional mechanical model for the microtubule-tau bundles in neuronal axons in which we remove taus, as can happen in various neurodegenerative conditions such as Alzheimer's disease, tauopathies, and chronic traumatic encephalopathy. Our simplified model includes (i) taus modeled as entropic springs between microtubules, (ii) removal of taus from the bundles due to phosphorylation, and (iii) a possible depletion force between microtubules due to these dissociated phosphorylated taus. We equilibrate upon tau removal using steepest descent relaxation. In the absence of the depletion force, the transverse rigidity to radial compression of the bundle falls to zero at about 60% tau occupancy, in agreement with standard percolation theory results. However, with the attractive depletion force, spring removal leads to a first order collapse of the bundles over a wide range of tau occupancies for physiologically realizable conditions. While our simplest calculations assume a constant concentration of microtubule intercalants to mediate the depletion force, including a dependence that is linear in the detached taus yields the same collapse. Applying percolation theory to removal of taus at microtubule tips, which are likely to be the protective sites against dynamic instability, we argue that the microtubule instability can only obtain at low tau occupancy, from 0.06-0.30 depending upon the tau coordination at the microtubule tips. Hence, the collapse we discover is likely to be more robust over a wide range of tau occupancies than the dynamic instability. We suggest in vitro tests of our predicted collapse.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Accuracy of topological entanglement entropy on finite cylinders

Topological phases are unique states of matter which support non-local excitations which behave as particles with fractional statistics. A universal characterization of gapped topological phases is provided by the topological entanglement entropy (TEE). We study the finite size corrections to the TEE by focusing on systems with $Z_2$ topological ordered state using density-matrix renormalization group (DMRG) and perturbative series expansions. We find that extrapolations of the TEE based on the Renyi entropies with Renyi index $n\geq 2$ suffer from much larger finite size corrections than do extrapolations based on the von Neumann entropy. In particular, when the circumference of the cylinder is about ten times the correlation length, the TEE obtained using von Neumann entropy has an error of order $10^{-3}$, while for Renyi entropies it can even exceed 40%. We discuss the relevance of these findings to previous and future searches for topological ordered phases, including quantum spin liquids.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Anomalies in the Entanglement Properties of the Square Lattice Heisenberg Model

We compute the bipartite entanglement properties of the spin-half square-lattice Heisenberg model by a variety of numerical techniques that include valence bond quantum Monte Carlo (QMC), stochastic series expansion QMC, high temperature series expansions and zero temperature coupling constant expansions around the Ising limit. We find that the area law is always satisfied, but in addition to the entanglement entropy per unit boundary length, there are other terms that depend logarithmically on the subregion size, arising from broken symmetry in the bulk and from the existence of corners at the boundary. We find that the numerical results are anomalous in several ways. First, the bulk term arising from broken symmetry deviates from an exact calculation that can be done for a mean-field Neel state. Second, the corner logs do not agree with the known results for non-interacting Boson modes. And, third, even the finite temperature mutual information shows an anomalous behavior as T goes to zero, suggesting that T->0 and L->infinity limits do not commute. These calculations show that entanglement entropy demonstrates a very rich behavior in d>1, which deserves further attention.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Entanglement at a Two-Dimensional Quantum Critical Point: a Numerical Linked Cluster Expansion Study

We develop a method to calculate the bipartite entanglement entropy of quantum models, in the thermodynamic limit, using a Numerical Linked Cluster Expansion (NLCE) involving only rectangular clusters. It is based on exact diagonalization of all n x m rectangular clusters at the interface between entangled subsystems A and B. We use it to obtain the Renyi entanglement entropy of the two-dimensional transverse field Ising model, for arbitrary real Renyi index alpha. Extrapolating these results as a function of the order of the calculation, we obtain universal pieces of the entanglement entropy associated with lines and corners at the quantum critical point. They show NLCE to be one of the few methods capable of accurately calculating universal properties of arbitrary Renyi entropies at higher dimensional critical points.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Entanglement in ground and excited states of gapped fermion systems and their relationship with fermi surface and thermodynamic equilibrium properties

We study bipartite entanglement entropies in the ground and excited states of model fermion systems, where a staggered potential, $μ_s$, induces a gap in the spectrum. Ground state entanglement entropies satisfy the `area law', and the `area-law' coefficient is found to diverge as a logarithm of the staggered potential, when the system has an extended Fermi surface at $μ_s=0$. On the square-lattice, we show that the coefficient of the logarithmic divergence depends on the fermi surface geometry and its orientation with respect to the real-space interface between subsystems and is related to the Widom conjecture as enunciated by Gioev and Klich (Phys. Rev. Lett. 96, 100503 (2006)). For point Fermi surfaces in two-dimension, the `area-law' coefficient stays finite as $μ_s\to 0$. The von Neumann entanglement entropy associated with the excited states follows a `volume law' and allows us to calculate an entropy density function s_{V}(e), which is substantially different from the thermodynamic entropy density function $s_{T}(e)$, when the lattice is bipartitioned into two equal subsystems but approaches the thermodynamic entropy density as the fraction of sites in the larger subsystem, that is integrated out, approaches unity.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Corrections to Pauling residual entropy and single tetrahedron based approximations for the pyrochlore lattice Ising antiferromagnet

We study corrections to single tetrahedron based approximations for the entropy, specific heat and uniform susceptibility of the pyrochlore lattice Ising antiferromagnet, by a Numerical Linked Cluster (NLC) expansion. In a tetrahedron based NLC, the first order gives the Pauling residual entropy of ${1\over 2}\log{3\over 2}\approx 0.20273$. A 16-th order NLC calculation changes the residual entropy to 0.205507 a correction of 1.37 percent over the Pauling value. At high temperatures, the accuracy of the calculations is verified by a high temperature series expansion. We find the corrections to the single tetrahedron approximations to be at most a few percent for all the thermodynamic properties.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Phase transitions in spin-orbital models with spin-space anisotropies for iron-pnictides: A study through Monte Carlo simulations

The common phase diagrams of superconducting iron pnictides show interesting material specificities in the structural and magnetic phase transitions. In some cases the two transitions are separate and second order, while in others they appear to happen concomitantly as a single first order transition. We explore these differences using Monte Carlo simulations of a two-dimensional Hamiltonian with coupled Heisenberg-spin and Ising-orbital degrees of freedom. In this spin-orbital model, the finite-temperature orbital-ordering transition results in a tetragonal-to-orthorhombic symmetry reduction and is associated with the structural transition in the iron-pnictide materials. With a zero or very small spin space anisotropy, the magnetic transition separates from the orbital one in temperature, and the orbital transition is found to be in the Ising universality class. With increasing anisotropy, the two transitions rapidly merge together and tend to become weakly first order. We also study the case of a single-ion anisotropy and propose that the preferred spin-orientation along the antiferromagnetic direction in these materials is driven by orbital order.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Thermodynamic singularities in the entanglement entropy at a 2D quantum critical point

We study the bipartite entanglement entropy of the two-dimensional (2D) transverse-field Ising model in the thermodynamic limit using series expansion methods. Expansions are developed for the Renyi entropy around both the small-field and large-field limits, allowing the separate calculation of the entanglement associated with lines and corners at the boundary between sub-systems. Series extrapolations are used to extract subleading power laws and logarithmic singularities as the quantum critical point is approached. In 1D, we find excellent agreement with exact results as well as quantum Monte Carlo simulations. In 2D, we find compelling evidence that the entanglement at a corner is significantly different from a free boson field theory. These results demonstrate the power of the series expansion method for calculating entanglement entropy in interacting systems, a fact that will be particularly useful in future searches for exotic quantum criticality in models with and without the sign problem.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Finite Temperature Critical Behavior of Mutual Information

We study mutual information for Renyi entropy of arbitrary index n, in interacting quantum systems at finite-temperature critical points, using high-temperature expansion, quantum Monte Carlo simulations and scaling theory. We find that for n>1, the critical behavior is manifest at two temperatures T_c and n*T_c. For the XXZ model with Ising anisotropy, the coefficient of the area-law has a t*ln(t) singularity, whereas the subleading correction from corners has a logarithmic divergence, with a coefficient related to the exact results of Cardy and Peschel. For T<n*T_c there is a constant term associated with broken symmetries that jumps at both T_c and n*T_c, which can be understood in terms of a scaling function analogous to the boundary entropy of Affleck and Ludwig.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Thermodynamics and phase transitions for the Heisenberg model on the pinwheel distorted kagome lattice

We study the Heisenberg model on the pinwheel distorted kagome lattice as observed in the material Rb_2Cu_3SnF_12. Experimentally relevant thermodynamic properties at finite temperatures are computed utilizing numerical linked-cluster expansions. We also develop a Lanczos-based, zero-temperature, numerical linked cluster expansion to study the approach of the pinwheel distorted lattice to the uniform kagome-lattice Heisenberg model. We find strong evidence for a phase transition before the uniform limit is reached, implying that the ground state of the kagome-lattice Heisenberg model is likely not pinwheel dimerized and is stable to finite pinwheel-dimerizing perturbations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Evaluating force field accuracy with long-time simulations of a tryptophan zipper peptide

We have combined a custom implementation of the fast multiple-time-stepping LN integrator with parallel tempering to explore folding properties of small peptides in implicit solvent on the time scale of microseconds. We applied this algorithm to the synthetic β-hairpin trpzip2 and one of its sequence variants W2W9. Each simulation consisted of over 12 μs of aggregated virtual time. Several measures of folding behavior showed convergence, allowing comparison with experimental equilibrium properties. Our simulations suggest that the electrostatic interaction of tryptophan sidechains is responsible for much of the stability of the native fold. We conclude that the ff99 force field combined with ff96 ϕ and ψ dihedral energies and implicit solvent can reproduce plausible folding behavior in both trpzip2 and W2W9.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Neutron spectra of Herbertsmithite Materials: Observation of a Valence Bond Liquid phase?

We argue that neutron spectra in a short-range ordered Valence Bond state is dominated by two-spinon states localized in small spatial regions such as a pinwheel. These excitations lead to angle averaged dynamic structure factor that is spread over a wide frequency range up to about $2.5 J$, whereas its wavevector dependence at all frequencies remains very close to that of isolated dimers. These results are in excellent agreement with recent Neutron scattering data in the Herbertsmithite materials ZnCu$_3$(OH)$_6$Cl$_2$.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Kagome lattice antiferromagnets and Dzyaloshinsky-Moriya interactions

We study the consequences of in-plane (D_p) and out-of-plane (D_z) Dzyaloshinsky-Moriya (DM) interactions on the thermodynamic properties of spin-1/2} Heisenberg model on the kagome lattice using numerical linked cluster expansions and exact diagonalization, and contrast them with those of other perturbations such as exchange anisotropy and dilution. We find that different combinations of the DM anisotropies lead to a wide variety of thermodynamic behavior, which are quite distinct from those of most other perturbations. We argue that the sudden upturn seen experimentally in the susceptibility of the material ZnCu3(OH)6Cl2 can be understood in terms of Dzyaloshinsky-Moriya anisotropies with D_p>|D_z|. We also show that the measured specific heat of the material puts further constraints on the allowed DM parameters.

0 citations0 reviewsNo code linkedOpen access
preprint1994arXiv

A Plane of Weakly Coupled Heisenberg Chains: Theoretical Arguments and Numerical Calculations

The $S=1/2$, nearest-neighbor, quantum Heisenberg antiferromagnet on the square lattice with spatially anisotropic couplings is reconsidered, with particular attention to the following question: at T=0, does Néel orderdevelop at infinitesimal interchain coupling, or is there a nonzero critical coupling? A heuristic renormalization group argument is presented which suggests that previous theoretical answers to that question are incorrect or at least incomplete, and that the answer is not universal but rather depends on the microscopic details of the model under consideration. Numerical investigations of the nearest-neighbor model are carried out {\it via} zero-temperature series expansions about Ising and dimer Hamiltonians. The results are entirely consistent with a vanishing critical interchain coupling ratio $R_c$; if $R_c$ is finite, it is unlikely to substantially exceed 0.02.

0 citations0 reviewsNo code linkedOpen access