Magnetism, exchange and crystal field parameters in the orbitally unquenched Ising antiferromagnet FePS3

FePS3 is a layered antiferromagnet (T N=123 K) with a marked Ising anisotropy in magnetic properties. The anisotropy arises from the combined effect of the trigonal distortion from octahedral symmetry and spin-orbit coupling on the orbitally degenerate5 T 2g ground state of the Fe2+ ion. The anisotropic paramagnetic susceptibilities are interpreted in terms of the zero field Hamiltonian, ?=?i [?(L iz 2 ?2)+|?|L i .S i ]?? ij J ij S i .S j . The crystal field trigonal distortion parameter ?, the spin-orbit coupling ? and the isotropic Heisenberg exchange,J ij, were evaluated from an analysis of the high temperature paramagnetic susceptibility data using the Correlated Effective Field (CEF) theory for many-body magnetism developed by Lines. Good agreement with experiment were obtained for ?/k=215.5 K; ?/k=166.5 K;J nn k=27.7 K; andJ nnn k=?2.3 K. Using these values of the crystal field and exchange parameters the CEF predicts aT N=122 K for FePS3, which is remarkably close to the observed value of theT N. The accuracy of the CEF approximation was also ascertained by comparing the calculated susceptibilities in the CEF with the experimental susceptibility for the isotropic Heisenberg layered antiferromagnet MnPS3, for which the high temperature series expansion susceptibility is available.

Density-matrix renormalization-group studies of the spin-1/2 Heisenberg system with dimerization and frustration

Using the density-matrix renormalization-group technique, we study the ground-state phase diagram and other low-energy properties of an isotropic antiferromagnetic spin-1/2 chain with both dimerization and frustration, i.e., an alternation delta of the nearest-neighbor exchanges and a next-nearest-neighbor exchange J(2). For delta = 0, the system is gapless for J(2) < J(2c) and has a gap for J(2) > J(2c) where J(2c) is about 0.241. For J(2) = J(2c) the gap above the ground state grows as delta to the power 0.667 +/- 0.001. In the J(2)-delta plane, there is a disorder line 2J(2) + delta = 1. To the left of this line, the peak in the static structure factor S(q) is at q(max) = pi (Neel phase), while to the right of the line, q(max) decreases from pi to pi/2 as J(2) is increased to large values (spiral phase). For delta = 1, the system is equivalent to two coupled chains as on a ladder and it is gapped for all values of the interchain coupling.

Regular versus alternating (FeS4)(n) chains: Magnetism in KFeS2 and CsFeS2

We report crystal magnetic susceptibility results of two S = 1/2 one-dimensional Heisenberg antiferromagnets, KFeS2 and CsFeS2. Both compounds consist of (FeS4)(n) chains with an average Fe-Fe distance of 2.7 Angstrom. In KFeS2, all intrachain Fe-Fe distances are identical. Its magnetic susceptibility is typical of a regular antiferromagnetic chain with spin-spin exchange parameter J = -440.7 K. In CsFeS2, however, the Fe-Fe distances alternate between 2.61 and 2.82 Angstrom. This is reflected in its magnetic susceptibility, which could be fitted with J = -640 K, and the degree of alternation, alpha = 0.3. These compounds form a unique pair, and allow for a convenient experimental comparison of the magnetic properties of regular versus alternating Heisenberg chains.

Ising pyrochlore magnets: Low-temperature properties, "ice rules," and beyond

Pyrochlore magnets are candidates for what Harris et al. [Phys. Rev. Lett. 79, 2554 (1997)] call "spin-ice" behavior. We present theoretical simulations of relevance for the pyrochlore family R2Ti2O7 (R = rare earth) supported by magnetothermal measurements on selected systems. Ey considering long-ranged dipole-dipole as well as short-ranged superexchange interactions, we get three distinct behaviors: (i) an ordered doubly degenerate state, (ii) a highly disordered state with a broad transition to paramagnetism, and (iii) a partially ordered state with a sharp transition to paramagnetism. Closely corresponding behavior is seen in the real compounds.

Two-channel Kondo physics in odd impurity chains

We study odd-membered chains of spin-1/2 impurities, with each end connected to its own metallic lead. For antiferromagnetic exchange coupling, universal two-channel Kondo (2CK) physics is shown to arise at low energies. Two overscreening mechanisms are found to occur depending on coupling strength, with distinct signatures in physical properties. For strong interimpurity coupling, a residual chain spin-1/2 moment experiences a renormalized effective coupling to the leads, while in the weak-coupling regime, Kondo coupling is mediated via incipient single-channel Kondo singlet formation. We also investigate models in which the leads are tunnel-coupled to the impurity chain, permitting variable dot filling under applied gate voltages. Effective low-energy models for each regime of filling are derived, and for even fillings where the chain ground state is a spin singlet, an orbital 2CK effect is found to be operative. Provided mirror symmetry is preserved, 2CK physics is shown to be wholly robust to variable dot filling; in particular, the single-particle spectrum at the Fermi level, and hence the low-temperature zero-bias conductance, is always pinned to half-unitarity. We derive a Friedel-Luttinger sum rule and from it show that, in contrast to a Fermi liquid, the Luttinger integral is nonzero and determined solely by the ``excess'' dot charge as controlled by gate voltage. The relevance of the work to real quantum dot devices, where interlead charge-transfer processes fatal to 2CK physics are present, is also discussed. Physical arguments and numerical renormalization-group techniques are used to obtain a detailed understanding of these problems.

Confinement-deconfinement transition and spin correlations in a generalized Kitaev model

We present a spin model, namely, the Kitaev model augmented by a loop term and perturbed by an Ising Hamiltonian, and show that it exhibits both confinement-deconfinement transitions from spin liquid to antiferromagnetic/spin-chain/ferromagnetic phases and topological quantum phase transitions between gapped and gapless spin-liquid phases. We develop a fermionic resonating-valence-bonds (RVB) mean-field theory to chart out the phase diagram of the model and estimate the stability of its spin-liquid phases, which might be relevant for attempts to realize the model in optical lattices and other spin systems. We present an analytical mean-field theory to study the confinement-deconfinement transition for large coefficient of the loop term and show that this transition is first order within such mean-field analysis in this limit. We also conjecture that in some other regimes, the confinement-deconfinement transitions in the model, predicted to be first order within the mean-field theory, may become second order via a defect condensation mechanism. Finally, we present a general classification of the perturbations to the Kitaev model on the basis of their effect on it's spin correlation functions and derive a necessary and sufficient condition, within the regime of validity of perturbation theory, for the spin correlators to exhibit a long-ranged power-law behavior in the presence of such perturbations. Our results reproduce those of Tikhonov et al. [Phys. Rev. Lett. 106, 067203 (2011)] as a special case.

Monte Carlo simulation of polytypes

Based on the analogy between polytypes and spin-half Ising chains with competing short- and infinite-range interactions, a Monte Carlo simulation of polytypes has been attempted. A general double-layer mechanism connects different states of the polytype chain with about the same probability as the spin-flip mechanism in magnetic Ising chains. It has been possible to simulate various polytypes with periodicities extending up to 12 layers. The Monte Carlo method should be useful in testing different interaction models that may be proposed in the future to describe polytypism.

Effect of dimerization on dynamics of spin-charge separation in Pariser-Parr-Pople model: A time-dependent density matrix renormalization group study

We investigate the effect of static electron-phonon coupling on real-time dynamics of spin and charge transport in pi-conjugated polyene chains. The polyene chain is modeled by the Pariser-Parr-Pople Hamiltonian with dimerized nearest-neighbor parameter t(0)(1 + delta) for short bonds and t(0)(1 - delta) for long bonds, and long-range electron-electron interactions. We follow the time evolution of the spin and charge using time-dependent density matrix renormalization group technique when a hole is injected at one end of the chain in its ground state. We find that spin and charge dynamics followed through spin and charge velocities depend both on chain length and extent of dimerization delta. Analysis of the results requires focusing on physical quantities such as average spin and charge polarizations, particularly in the large dimerization limit. In the dimerization range 0.0 <= delta <= 0.15, spin-charge dynamics is found to have a well-defined behavior, with spin-charge separation (measured as the ratio of charge velocity to spin velocity) as well as the total amount of charge and spin transported in a given time along the chain decreasing as dimerization increases. However, in the range 0.3 <= delta <= 0.5, it is observed that the dynamics of spin and charge transport becomes complicated. It is observed that, for large delta values, spin-charge separation is suppressed and the injected hole fails to travel the entire length of the chain.

Real-time density matrix renormalization group dynamics of spin and charge transport in push-pull polyenes and related systems

In this paper we investigate the effect of terminal substituents on the dynamics of spin and charge transport in donor-acceptor substituted polyenes [D-(CH)(x)-A] chains, also known as push-pull polyenes. We employ a long-range correlated model Hamiltonian for the D-(CH)(x)-A system, and time-dependent density matrix renormalization group technique for time propagating the wave packet obtained by injecting a hole at a terminal site, in the ground state of the system. Our studies reveal that the end groups do not affect spin and charge velocities in any significant way, but change the amount of charge transported. We have compared these push-pull systems with donor-acceptor substituted polymethine imine (PMI), D-(CHN)(x)-A, systems in which besides electron affinities, the nature of p(z) orbitals in conjugation also alternate from site to site. We note that spin and charge dynamics in the PMIs are very different from that observed in the case of push-pull polyenes, and within the time scale of our studies, transport of spin and charge leads to the formation of a ``quasi-static'' state.

Two Novel Heterometallic Chains Featuring Mn-II and Na-I Ions in Trigonal-Prismatic Geometries Alternately Linked to Octahedral Mn-IV Ions: Synthesis, Structures, and Magnetic Behavior

Two new one-dimensional heterometallic complexes, Mn3Na(L)(4)(CH3CO2)(MeOH)(2)]-(ClO4)(2)center dot 3H(2)O (1), Mn3Na(L)(4)(CH3CH2CO2)-(MeOH)(2)](ClO4)(2)center dot 2MeOH center dot H2O (2) LH2 = 2-methyl-2-(2-pyridyl)propane-1,3-diol], have been synthesized and characterized by X-ray crystallography. Both complexes feature Mn-II and Na-I ions in trigonal-prismatic geometries that are linked to octahedral Mn-IV ions by alkoxy bridges. Variable-temperature direct- and alternating-current magnetic susceptibility data indicated a spin ground state of S = 11/2 for both complexes. Density functional theory calculations performed on 1 supported this conclusion.

Quantum Phase Diagram of One-Dimensional Spin and Hubbard Models with Transitions to Bond Order Wave Phases

Similar quantum phase diagrams and transitions are found for three classes of one-dimensional models with equally spaced sites, singlet ground states (GS), inversion symmetry at sites and a bond order wave (BOW) phase in some sectors. The models are frustrated spin-1/2 chains with variable range exchange, half-filled Hubbard models with spin-independent interactions and modified Hubbard models with site energies for describing organic charge transfer salts. In some range of parameters, the models have a first order quantum transition at which the GS expectation value of the sublattice spin < S-A(2)> of odd or even-numbered sites is discontinuous. There is an intermediate BOW phase for other model parameters that lead to two continuous quantum transitions with continuous < S-A(2)>. Exact diagonalization of finite systems and symmetry arguments provide a unified picture of familiar 1D models that have appeared separately in widely different contexts.

Quantum Phase Diagram of One-Dimensional Spin and Hubbard Models with Transitions to Bond Order Wave Phases

Similar quantum phase diagrams and transitions are found for three classes of one-dimensional models with equally spaced sites, singlet ground states (GS), inversion symmetry at sites and a bond order wave (BOW) phase in some sectors. The models are frustrated spin-1/2 chains with variable range exchange, half-filled Hubbard models with spin-independent interactions and modified Hubbard models with site energies for describing organic charge transfer salts. In some range of parameters, the models have a first order quantum transition at which the GS expectation value of the sublattice spin < S-A(2)> of odd or even-numbered sites is discontinuous. There is an intermediate BOW phase for other model parameters that lead to two continuous quantum transitions with continuous < S-A(2)>. Exact diagonalization of finite systems and symmetry arguments provide a unified picture of familiar 1D models that have appeared separately in widely different contexts.

Studies on a frustrated Heisenberg spin chain with alternating ferromagnetic and antiferromagnetic exchanges

We study Heisenberg spin-1/2 and spin-1 chains with alternating ferromagnetic (J(1)(F)) and antiferromagnetic (J(1)(A)) nearest-neighbor interactions and a ferromagnetic next-nearest-neighbor interaction (J(2)(F)). In this model frustration is present due to the non-zero J(2)(F). The model with site spin s behaves like a Haldane spin chain, with site spin 2s in the limit of vanishing J(2)(F) and large J(1)(F)/J(1)(A). We show that the exact ground state of the model can be found along a line in the parameter space. For fixed J(1)(F), the phase diagram in the space of J(1)(A)-J(2)(F) is determined using numerical techniques complemented by analytical calculations. A number of quantities, including the structure factor, energy gap, entanglement entropy and zero temperature magnetization, are studied to understand the complete phase diagram. An interesting and potentially important feature of this model is that it can exhibit a macroscopic magnetization jump in the presence of a magnetic field; we study this using an effective Hamiltonian.

Ising model on a random network with annealed or quenched disorder

We study the equilibrium properties of an Ising model on a disordered random network where the disorder can be quenched or annealed. The network consists of fourfold coordinated sites connected via variable length one-dimensional chains. Our emphasis is on nonuniversal properties and we consider the transition temperature and other equilibrium thermodynamic properties, including those associated with one-dimensional fluctuations arising from the chains. We use analytic methods in the annealed case, and a Monte Carlo simulation for the quenched disorder. Our objective is to study the difference between quenched and annealed results with a broad random distribution of interaction parameters. The former represents a situation where the time scale associated with the randomness is very long and the corresponding degrees of freedom can be viewed as frozen, while the annealed case models the situation where this is not so. We find that the transition temperature and the entropy associated with one-dimensional fluctuations are always higher for quenched disorder than in the annealed case. These differences increase with the strength of the disorder up to a saturating value. We discuss our results in connection to physical systems where a broad distribution of interaction strengths is present.

Efficient density matrix renormalization group algorithm to study Y junctions with integer and half-integer spin

An efficient density matrix renormalization group (DMRG) algorithm is presented and applied to Y junctions, systems with three arms of n sites that meet at a central site. The accuracy is comparable to DMRG of chains. As in chains, new sites are always bonded to the most recently added sites and the superblock Hamiltonian contains only new or once renormalized operators. Junctions of up to N = 3n + 1 approximate to 500 sites are studied with antiferromagnetic (AF) Heisenberg exchange J between nearest-neighbor spins S or electron transfer t between nearest neighbors in half-filled Hubbard models. Exchange or electron transfer is exclusively between sites in two sublattices with N-A not equal N-B. The ground state (GS) and spin densities rho(r) = < S-r(z)> at site r are quite different for junctions with S = 1/2, 1, 3/2, and 2. The GS has finite total spin S-G = 2S(S) for even (odd) N and for M-G = S-G in the S-G spin manifold, rho(r) > 0(< 0) at sites of the larger (smaller) sublattice. S = 1/2 junctions have delocalized states and decreasing spin densities with increasing N. S = 1 junctions have four localized S-z = 1/2 states at the end of each arm and centered on the junction, consistent with localized states in S = 1 chains with finite Haldane gap. The GS of S = 3/2 or 2 junctions of up to 500 spins is a spin density wave with increased amplitude at the ends of arms or near the junction. Quantum fluctuations completely suppress AF order in S = 1/2 or 1 junctions, as well as in half-filled Hubbard junctions, but reduce rather than suppress AF order in S = 3/2 or 2 junctions.