Search

This tutorial is a pedagogical introduction to recent methods of computing quantum spin entanglement witnesses from spectroscopy, with a special focus on neutron scattering on quantum spin systems. We offer a brief introduction to the concepts and equations, define a data analysis protocol, and discuss the interpretation of three entanglement witnesses: one-tangle, two-tangle, and Quantum Fisher Information. We also discuss practical experimental considerations, and give three examples of extracting entanglement witnesses from experimental data: Copper Nitrate, KCuF3, and NiPS3., Comment: 11 pages, 5 figures

We study the impact of hole-doping on the Kitaev-Heisenberg model on the honeycomb lattice. We investigate the pairing tendencies and correlation functions in the framework of a $t-J-K$ model using density matrix renormalization group calculations on three-leg cylinders. In the case of the pure Kitaev model, which realizes a quantum spin-liquid phase at half-filling, we find that binding of two holes only occurs at low values of the hopping, where the holes are slow. We have theoretically verified that pair formation occurs in the limit of immobile holes, where the pure Kitaev model remains exactly solvable. When we instead fix the hopping at an intermediate, more realistic, value, and vary the Heisenberg and Kitaev interaction strengths, we find pairing tendencies only in the N\'eel phase. This is in contrast to prior mean-field calculations, highlighting the importance of accounting for the kinetic energy of dopants in generalized Kitaev models. Interestingly, we also find signatures of pair-density wave formation over the studied range of model parameters, namely a periodic modulation of the charge density as well as the spin-spin and pair-pair correlations in real space. Moreover, we present a comparative study of the different correlations as a function of doping. We finally discuss the potential for experimentally observing the studied physics in quantum materials and heterostructures., Comment: Main text: 12 pages, 7 figures. Supplemental material: 5 pages, 1 figure

The detection and certification of entanglement and quantum correlations in materials is of fundamental and far-reaching importance, and has seen significant recent progress. It impacts both our understanding of the basic science of quantum many-body phenomena as well as the identification of systems suitable for novel technologies. Frameworks suitable to condensed matter that connect measurements to entanglement and coherence have been developed in the context of quantum information theory. These take the form of entanglement witnesses and quantum correlation measures. The underlying theory of these quantities, their relation to condensed matter experimental techniques, and their application to real materials are comprehensively reviewed. In addition, their usage in e.g. protocols, the relative advantages and disadvantages of witnesses and measures, and future prospects in, e.g., correlated electrons, entanglement dynamics, and entangled spectroscopic probes, are presented. Consideration is given to the interdisciplinary nature of this emerging research and substantial ongoing progress by providing an accessible and practical treatment from fundamentals to application. Particular emphasis is placed on quantities accessible to collective measurements, including by susceptibility and spectroscopic techniques. This includes the magnetic susceptibility witness, one-tangle, concurrence and two-tangle, two-site quantum discord, and quantum coherence measures such as the quantum Fisher information., Comment: 41 pages, 13 figures, 2 tables

We investigate the pairing tendencies in the hole-doped Haldane spin-1 chain. To allow for doping, we extend the original spin chain Hamiltonian into a fermionic model involving a two-orbital Hubbard chain at intermediate or strong repulsive interaction strengths $U$, and for degenerate orbitals. At half-filling and large $U$, the ferromagnetic Hund's coupling, $J_\mathrm{H}$, generates effective spin-$1$ moments, with antiferromagnetic correlations between sites. Using large-scale density matrix renormalization group calculations, we study accurately the system's behavior under light hole-doping. For $U=1.6$ in units of the non-interacting bandwidth and for $J_\mathrm{H}/U\gtrsim 0.275$ we find that singlet pairing dominates the long-distance physics, establishing this system as a promising platform for repulsively mediated superconductivity. We provide concrete examples of materials that could realize the physics described here. We also provide evidence that the system approaches a Luther-Emery liquid state at large system sizes, reminiscent of the behavior of doped one-orbital two-leg ladders at weak coupling, which also have superconducting tendencies. The numerically calculated central charge approaches one in the thermodynamic limit, indicating a single gapless mode as is expected for the Luther-Emery state. Exponents characterizing the power-law decays of singlet pair-pair and charge density-density correlations are determined, and found to approximately satisfy the Luther-Emery identity., Comment: Main text: 10 pages, 9 figures. Supplementary: 4 pages, 2 figures

Magnetic excitations in the spin chain candidate Sr$_2$V$_3$O$_9$ have been investigated by inelastic neutron scattering on a single crystal sample. A spinon continuum with a bandwidth of $\sim22$ meV is observed along the chain formed by alternating magnetic V$^{4+}$ and nonmagnetic V$^{5+}$ ions. Incipient magnetic Bragg peaks due to weak ferromagnetic interchain couplings emerge when approaching the magnetic transition at $T_N\sim 5.3$ K while the excitations remain gapless within the instrumental resolution. Comparisons to the Bethe ansatz, density matrix renormalization group (DMRG) calculations, and effective field theories confirm Sr$_2$V$_3$O$_9$ as a host of weakly coupled $S = 1/2$ chains dominated by antiferromagnetic intrachain interactions of $\sim7.1$(1) meV., Comment: 19 pages, 9 figures

Quantum correlations are a fundamental property of quantum many-body states. Yet they remain experimentally elusive, hindering certification of genuine quantum behavior, especially in quantum materials. Here we show that the momentum-dependent dynamical susceptibility measured via inelastic neutron scattering enables the systematic reconstruction of a general family of quantum correlation functions, which express the degree of quantum coherence in the fluctuations of two spins at arbitrary mutual distance. Using neutron scattering data on the compound KCuF$_3$ $\unicode{x2014}$ a system of weakly coupled $S=1/2$ Heisenberg chains $\unicode{x2014}$ and of numerically exact quantum Monte Carlo data, we show that quantum correlations possess a radically different spatial structure with respect to conventional correlations. Indeed, they exhibit a new emergent length scale $\unicode{x2014}$ the quantum coherence length $\unicode{x2014}$ which is finite at any finite temperature (including when long-range magnetic order develops). Moreover, we show theoretically that coupled Heisenberg spin chains exhibit a form of quantum monogamy, with a trade-off between quantum correlations along and transverse to the spin chains. These results highlight real-space quantum correlators as an informative, model-independent means of probing the underlying quantum state of real quantum materials., Comment: 12 pages, 10 figures

Fully gapped, spin singlet superconductors with antisymmetric spin-orbit coupling in a Zeeman magnetic field provide a promising route to realize superconducting states with non-Abelian topological order and therefore fault-tolerant quantum computation. Here we use a quantum Monte Carlo dynamical cluster approximation to study the superconducting properties of a doped two-dimensional attractive Hubbard model with Rashba spin-orbit coupling in a Zeeman magnetic field. We generally find that the Rashba coupling has a beneficial effect towards s-wave superconductivity. In the presence of a finite Zeeman field, when superconductivity is suppressed by Pauli pair-breaking, the Rashba coupling counteracts the spin imbalance created by the Zeeman field by mixing the spins, and thus restores superconductivity at finite temperatures. We show that this favorable effect of the spin-orbit coupling is traced to a spin-flip driven enhancement of the amplitude for the propagation of a pair of electrons in time-reversed states. Moreover, by inspecting the Fermi surface of the interacting model, we show that for sufficiently large Rashba coupling and Zeeman field, the superconducting state is expected to be topologically non-trivial.

We calculate the dynamical spin structure factor of the generalized spin-$1/2$ compass spin chain using the density matrix renormalization group. The model, also known as the twisted Kitaev spin chain, was recently proposed to be relevant for the description of the spin chain compound CoNb$_2$O$_6$. It features bond-dependent interactions and interpolates between an Ising chain and a one-dimensional variant of Kitaev's honeycomb spin model. The structure factor, in turn, is found to interpolate from gapped and non-dispersive in the Ising limit to gapless with non-trivial continua in the Kitaev limit. In particular, the component of the structure factor perpendicular to the Ising directions changes abruptly at the Kitaev point into a dispersionless continuum due to the emergence of an extensive groundstate degeneracy. We show this continuum is consistent with analytical Jordan-Wigner results. We also discuss implications for future inelastic scattering experiments and applications to materials, particularly CoNb$_2$O$_6$., Comment: Main text: 10 pages, 6 figures. Supplemental material: 5 pages, 2 figures

Recent work has demonstrated that quantum Fisher information (QFI), a witness of multipartite entanglement, and magnetic Van Hove correlations $G(r,t)$, a probe of local real-space real-time spin dynamics, can be successfully extracted from inelastic neutron scattering on spin systems through accurate measurements of the dynamical spin structure factor $S(k,\omega)$. Here we apply theoretically these ideas to the half-filled Hubbard chain with nearest-neighbor hopping, away from the strong-coupling limit. This model has nontrivial redistribution of spectral weight in $S(k,\omega)$ going from the non-interacting limit ($U=0$) to strong coupling ($U\rightarrow \infty$), where it reduces to the Heisenberg quantum spin chain. We use the density matrix renormalization group (DMRG) to find $S(k,\omega)$, from which QFI is then calculated. We find that QFI grows with $U$. With realistic energy resolution it becomes capable of witnessing bipartite entanglement above $U=2.5$ (in units of the hopping), where it also changes slope. This point is also proximate to slope changes of the bandwidth $W(U)$ and the half-chain von Neumann entanglement entropy. We compute $G(r,t)$ by Fourier-transforming $S(k,\omega)$. The results indicate a crossover in the short-time short-distance dynamics at low $U$ characterized by ferromagnetic lightcone wavefronts, to a Heisenberg-like behavior at large $U$ featuring antiferromagnetic lightcones and spatially period-doubled antiferromagnetism. We find this crossover has largely been completed by $U=3$. Our results thus provide evidence that, in several aspects, the strong-coupling limit of the Hubbard chain is reached qualitatively already at a relatively modest interaction strength. We discuss experimental candidates for observing the $G(r,t)$ dynamics found at low $U$., Comment: 15 pages, 10 figures

Single crystal inelastic neutron scattering data contain rich information about the structure and dynamics of a material. Yet the challenge of matching sophisticated theoretical models with large data volumes is compounded by computational complexity and the ill-posed nature of the inverse scattering problem. Here we utilize a novel machine-learning-assisted framework featuring multiple neural network architectures to address this via high-dimensional modeling and numerical methods. A comprehensive data set of diffraction and inelastic neutron scattering measured on the Kitaev material $\alpha-$RuCl$_3$ is processed to extract its Hamiltonian. Semiclassical Landau-Lifshitz dynamics and Monte-Carlo simulations were employed to explore the parameter space of an extended Kitaev-Heisenberg Hamiltonian. A machine-learning-assisted iterative algorithm was developed to map the uncertainty manifold to match experimental data; a non-linear autoencoder used to undertake information compression; and Radial Basis networks utilized as fast surrogates for diffraction and dynamics simulations to predict potential spin Hamiltonians with uncertainty. Exact diagonalization calculations were employed to assess the impact of quantum fluctuations on the selected parameters around the best prediction., Comment: 22 pages, 18 figures

Traditional spectroscopy, by its very nature, characterizes properties of physical systems in the momentum and frequency domains. The most interesting and potentially practically useful quantum many-body effects however emerge from the deep composition of local, short-time correlations. Here, using inelastic neutron scattering and methods of integrability, we experimentally observe and theoretically describe a local, coherent, long-lived, quasiperiodically oscillating magnetic state emerging out of the distillation of propagating excitations following a local quantum quench in a Heisenberg antiferromagnetic chain. This "quantum wake" displays similarities to Floquet states, discrete time crystals and nonlinear Luttinger liquids., Comment: 7 pages 4 figures, with 4 pages supplemental information

We investigate the magnetic excitations of elemental gadolinium (Gd) using inelastic neutron scattering, showing that Gd is a Dirac magnon material with nodal lines at $K$ and nodal planes at half integer $\ell$. We find an anisotropic intensity winding around the $K$-point Dirac magnon cone, which is interpreted to indicate Berry phase physics. Using linear spin wave theory calculations, we show the nodal lines have non-trivial Berry phases, and topological surface modes. We also discuss the origin of the nodal plane in terms of a screw-axis symmetry, and introduce a topological invariant characterizing its presence and effect on the scattering intensity. Together, these results indicate a highly nontrivial topology, which is generic to hexagonal close packed ferromagnets. We discuss potential implications for other such systems., Comment: 4 pages, 4 figures, plus four pages supplemental information. Accompanying paper to arXiv:2107.11370

We present a comprehensive study of the magnetic exchange Hamiltonian of elemental Gadolinium. We use neutron scattering to measure the magnon spectrum over the entire Brillouin zone, and fit the excitations to a spin wave model to extract the first 26 nearest neighbor magnetic exchange interactions with rigorously defined uncertainty. We find these exchange interactions to follow RKKY behavior, oscillating from ferromagnetic to antiferromagnetic as a function of distance. Finally, we discuss the topological features and degeneracies in Gd, and HCP ferromagnets in general. We show theoretically how, with asymmetric exchange, topological properties could be tuned with a magnetic field., Comment: 6 pages, 10 figures, 5 pages appendices. Accompanying paper: arXiv:2107.11372

We demonstrate how quantum entanglement can be directly witnessed in the quasi-1D Heisenberg antiferromagnet KCuF$_3$. We apply three entanglement witnesses --- one-tangle, two-tangle, and quantum Fisher information --- to its inelastic neutron spectrum, and compare with spectra simulated by finite-temperature density matrix renormalization group (DMRG) and classical Monte Carlo methods. We find that each witness provides direct access to entanglement. Of these, quantum Fisher information is the most robust experimentally, and indicates the presence of at least bipartite entanglement up to at least 50 K, corresponding to around 10% of the spinon zone-boundary energy. We apply quantum Fisher information to higher spin-S Heisenberg chains, and show theoretically that the witnessable entanglement gets suppressed to lower temperatures as the quantum number increases. Finally, we outline how these results can be applied to higher dimensional quantum materials to witness and quantify entanglement., Comment: 9 pages and 5 figures, four pages and six figures of appendices, and five pages supplemental information

The lack of methods to experimentally detect and quantify entanglement in quantum matter impedes our ability to identify materials hosting highly entangled phases, such as quantum spin liquids. We thus investigate the feasibility of using inelastic neutron scattering (INS) to implement a model-independent measurement protocol for entanglement based on three entanglement witnesses: one-tangle, two-tangle, and quantum Fisher information (QFI). We perform high-resolution INS measurements on Cs$_2$CoCl$_4$, a close realization of the $S=1/2$ transverse-field XXZ spin chain, where we can control entanglement using the magnetic field, and compare with density-matrix renormalization group calculations for validation. The three witnesses allow us to infer entanglement properties and make deductions about the quantum state in the material. We find QFI to be a particularly robust experimental probe of entanglement, whereas the one- and two-tangles require more careful analysis. Our results lay the foundation for a general entanglement detection protocol for quantum spin systems., Comment: Main text: 7 pages, 4 figures. Supplementary Information: 15 pages, 15 figures

Operator square-roots are ubiquitous in theoretical physics. They appear, for example, in the Holstein-Primakoff representation of spin operators and in the Klein-Gordon equation. Often the use of a perturbative expansion is the only recourse when dealing with them. In this work we show that under certain conditions differential equations can be derived which can be used to find perturbatively inaccessible approximations to operator square-roots. Specifically, for the number operator $\hat n=\hat a^\dag a$ we show that the square-root $\sqrt{\hat n}$ near $\hat n=0$ can be approximated by a polynomial in $\hat n$. This result is unexpected because a Taylor expansion fails. A polynomial expression in $\hat n$ is possible because $\hat n$ is an operator, and its constituents $a$ and $a^\dag$ have a non-trivial commutator $[a,a^\dag]=1$ and do not behave as scalars. We apply our approach to the zero mass Klein-Gordon Hamiltonian in a constant magnetic field, and as a main application, the Holstein-Primakoff representation of spin operators, where we are able to find new expressions that are polynomial in bosonic operators. We prove that these new expressions exactly reproduce spin operators. Our expressions are manifestly Hermitian, which offer an advantage over other methods, such as the Dyson-Maleev representation.

The magnetic excitation spectrum of the quantum magnet YbCl$_3$ is studied with inelastic neutron scattering. The spectrum exhibits an unusually sharp feature within a broad continuum, as well as conventional spin waves. By including both transverse and longitudinal channels of the neutron response, linear spin wave theory with a single Heisenberg interaction on the honeycomb lattice reproduces all of the key features in the spectrum. In particular, the broad continuum corresponds to a two-magnon contribution from the longitudinal channel, while the sharp feature within this continuum is identified as a Van Hove singularity in the joint density of states, which indicates the two-dimensional nature of the two-magnon continuum. We term these singularities magneto-caustic features in analogy with caustic features in ray optics where focused envelopes of light are generated when light passes through or reflects from curved or distorted surfaces. The experimental demonstration of a sharp Van Hove singularity in a two-magnon continuum is important because analogous features in potential two-spinon continua could distinguish quantum spin liquids from merely disordered systems. These results establish YbCl$_3$ as a nearly ideal two-dimensional honeycomb lattice material hosting strong quantum effects in the unfrustrated limit., Comment: 17 pages, 5 figures

We numerically investigate the momentum-space entanglement entropy and entanglement spectrum of the random-dimer model and its generalizations, which circumvent Anderson localization, after a quench in the Hamiltonian parameters. The type of dynamics that occurs depends on whether or not the Fermi level of the initial state is near the energy of the delocalized states present in these models. If the Fermi level of the initial state is near the energy of the delocalized states, we observe an interesting slow logarithmic-like growth of the momentum-space entanglement entropy followed by an eventual saturation. Otherwise, the momentum-space entanglement entropy is found to rapidly saturate. We also find that the momentum-space entanglement spectrum reveals the presence of delocalized states in these models for long times after the quench and the many-body entanglement gap decays logarithmically in time when the Fermi level is near the energy of the delocalized states., Comment: 4+e pages, 3 figures

What is the correct low-energy spin Hamiltonian description of $\alpha$-RuCl$_3$? The material is a promising Kitaev spin liquid candidate, but is also known to order magnetically, the description of which necessitates additional interaction terms. The nature of these interactions, their magnitudes and even signs, remain an open question. In this work we systematically investigate dynamical and thermodynamic magnetic properties of proposed effective Hamiltonians. We calculate zero-temperature inelastic neutron scattering (INS) intensities using exact diagonalization, and magnetic specific heat using a thermal pure quantum states method. We find that no single current model satisfactorily explains all observed phenomena of $\alpha$-RuCl$_3$. In particular, we find that Hamiltonians derived from first principles can capture the experimentally observed high-temperature peak in the magnetic specific heat, while overestimating the magnon energy at the zone center. In contrast, other models reproduce important features of the INS data, but do not adequately describe the magnetic specific heat. To address this discrepancy we propose a modified ab initio model that is consistent with both magnetic specific heat and low-energy features of INS data., Comment: Main text: 12 pages, 8 figures. Supplemental Information: 16 pages, 18 figures

We use the constrained random phase approximation (cRPA) to derive from first principles the Ru-$t_{2g}$ Wannier function based model for the Kitaev spin-liquid candidate material $\alpha$-RuCl$_3$. We find the non-local Coulomb repulsion to be sizable compared to the local one. In addition we obtain the contribution to the Hamiltonian from the spin-orbit coupling and find it to also contain non-negligible non-local terms. We invoke strong coupling perturbation theory to investigate the influence of these non-local elements of the Coulomb repulsion and the spin-orbit coupling on the magnetic interactions. We find that the non-local Coulomb repulsions cause a strong enhancement of the magnetic interactions, which deviate from experimental fits reported in the literature. Our results contribute to the understanding and design of quantum spin liquid materials via first principles calculations.

In this paper we develop an analogue of Hamilton-Jacobi theory for the time-evolution operator of a quantum many-particle system. The theory offers a useful approach to develop approximations to the time-evolution operator, and also provides a unified framework and starting point for many well-known approximations to the time-evolution operator. In the important special case of periodically driven systems at stroboscopic times, we find relatively simple equations for the coupling constants of the Floquet Hamiltonian, where a straightforward truncation of the couplings leads to a powerful class of approximations. Using our theory, we construct a flow chart that illustrates the connection between various common approximations, which also highlights some missing connections and associated approximation schemes. These missing connections turn out to imply an analytically accessible approximation that is the "inverse" of a rotating frame approximation and thus has a range of validity complementary to it. We numerically test the various methods on the one-dimensional Ising model to confirm the ranges of validity that one would expect from the approximations used. The theory provides a map of the relations between the growing number of approximations for the time-evolution operator. We describe these relations in a table showing the limitations and advantages of many common approximations, as well as the new approximations introduced in this paper., Comment: 17 pages, 5 figures, 1 table

We present a theoretical method to generate a highly accurate {\em time-independent} Hamiltonian governing the finite-time behavior of a time-periodic system. The method exploits infinitesimal unitary transformation steps, from which renormalization group-like flow equations are derived to produce the effective Hamiltonian. Our tractable method has a range of validity reaching into frequency regimes that are usually inaccessible via high frequency $\omega$ expansions in the parameter $h/\omega$, where $h$ is the upper limit for the strength of local interactions. We demonstrate our approach on both interacting and non-interacting many-body Hamiltonians where it offers an improvement over the more well-known Magnus expansion and other high frequency expansions. For the interacting models, we compare our approximate results to those found via exact diagonalization. While the approximation generally performs better globally than other high frequency approximations, the improvement is especially pronounced in the regime of lower frequencies and strong external driving. This regime is of special interest because of its proximity to the resonant regime where the effect of a periodic drive is the most dramatic. Our results open a new route towards identifying novel non-equilibrium regimes and behaviors in driven quantum many-particle systems., Comment: 25 pages, 14 figures

We theoretically study magnetic and topological properties of antiferromagnetic kagome spin systems in the presence of both in- and out-of-plane Dzyaloshinskii-Moriya interactions. In materials such as the iron jarosites, the in-plane interactions stabilize a canted noncollinear "umbrella" magnetic configuration with finite scalar spin chirality. We derive expressions for the canting angle, and use the resulting order as a starting point for a spin-wave analysis. We find topological magnon bands, characterized by non-zero Chern numbers. We calculate the magnon thermal Hall conductivity, and propose the iron jarosites as a promising candidate system for observing the magnon thermal Hall effect in a noncollinear spin configuration. We also show that the thermal conductivity can be tuned by varying an applied magnetic field, or the in-plane Dzyaloshinskii-Moriya strength. In contrast with previous studies of topological magnon bands, the effect is found to be large even in the limit of small canting., Comment: 11 pages, 11 figures

We theoretically study the magnetic properties of pyrochlore iridate bilayer and trilayer thin films grown along the $[111]$ direction using a strong coupling approach. We find the ground state magnetic configurations on a mean field level and carry out a spin-wave analysis about them. In the trilayer case the ground state is found to be the all-in-all-out (AIAO) state, whereas the bilayer has a deformed AIAO state. For all parameters of the spin-orbit coupled Hamiltonian we study, the lowest magnon band in the trilayer case has a nonzero Chern number. In the bilayer case we also find a parameter range with nonzero Chern numbers. We calculate the magnon Hall response for both geometries, finding a striking sign change as function of temperature. Using a slave-boson mean-field theory we study the doping of the trilayer system and discover an unconventional time-reversal symmetry broken $d+id$ superconducting state. Our study complements prior work in the weak coupling limit and suggests that the $[111]$ grown thin film pyrochlore iridates are a promising candidate for topological properties and unconventional orders., Comment: Main text: 6 pages, 4 figures. Supplementary material: 7 pages, 6 figures. Revised version implementing corrections from Erratum Phys. Rev. Lett. 121, 019901 (2018)

We advocate that in critical spin chains, and possibly in a larger class of 1D critical models, a gap in the momentum-space entanglement spectrum separates the universal part of the spectrum, which is determined by the associated conformal field theory, from the non-universal part, which is specific to the model. To this end, we provide affirmative evidence from multicritical spin chains with low energy sectors described by the SU(2)$_2$ or the SU(3)$_1$ Wess-Zumino-Witten model., Comment: 5 pages, 5 figures

We study the electronic contribution to the thermal conductivity and the thermopower of Weyl and Dirac semimetals using a semiclassical Boltzmann approach. We investigate the effect of various relaxation processes including disorder and interactions on the thermoelectric properties, and also consider doping away from the Weyl or Dirac point. We find that the thermal conductivity and thermopower have an interesting dependence on the chemical potential that is characteristic of the linear electronic dispersion, and that the electron-electron interactions modify the Lorenz number. For the interacting system, we also use the Kubo formalism to obtain the transport coefficients. We find exact agreement between the Kubo and Boltzmann approaches at high temperatures. We also consider the effect of electric and magnetic fields on the thermal conductivity in various orientations with respect to the temperature gradient. Notably, when the temperature gradient and magnetic field are parallel, we find a large contribution to the longitudinal thermal conductivity that is quadratic in the magnetic field strength, similar to the magnetic field dependence of the longitudinal electrical conductivity due to the presence of the chiral anomaly when no thermal gradient is present., Comment: 17 pages, 2 Figures