Source author record

Zoltán G. Soos

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

Catalog footprint

Research overview

9works
3topics
4close collaborators

Published work

9 publications

preprint2022arXiv

Density matrix renormalization group approach to the low temperature thermodynamics of correlated 1D fermionic models

The low temperature thermodynamics of correlated 1D fermionic models with spin and charge degrees of freedom is obtained by exact diagonalization (ED) of small systems and followed by density matrix renormalization group (DMRG) calculations that target the lowest hundreds of states {E(N)}\{E(N)\} at system size NN instead of the ground state. Progressively larger NN reaches T<0.05tT < 0.05t in correlated models with electron transfer tt between first neighbors and bandwidth 4t4t. The size dependence of the many-fermion basis is explicitly included for arbitrary interactions by scaling the partition function. The remaining size dependence is then entirely due to the energy spectrum {E(N)}\{E(N)\} of the model. The ED/DMRG method is applied to Hubbard and extended Hubbard models, both gapped and gapless, with Ne=NN_e = N or N/2N/2 electrons and is validated against exact results for the magnetic susceptibility χ(T)χ(T) and entropy S(T)S(T) per site. Some limitations of the method are noted. Special attention is given to the bond-order-wave phase of the extended Hubbard model with competing interactions and low TT thermodynamics sensitive to small gaps.

preprint2022arXiv

Spin Peierls transition of J1−J2J_{1}-J_{2} and extended models with ferromagnetic J1J_{1}.Sublattice dimerization and thermodynamics of zigzag chains in ββ-TeVO4_{4}

The spin−1/2-1/2 chain with ferromagnetic exchange J1<0J_1 < 0 between first neighbors and antiferromagnetic J2>0J_2 > 0 between second neighbors supports two spin-Peierls (SP) instabilities depending on the frustration α=J2/∣J1∣α= J_2/\vert J_1\vert. Instead of chain dimerization with two spins per unit cell, J1−J2J_1-J_2 models with α>0.65α> 0.65 and linear spin-phonon coupling are unconditionally unstable to sublattice dimerization with four spins per unit cell. Unequal J1J_1 to neighbors to the right and left extends the model to gapped (γ>0γ> 0) chains with conditional SP transitions at TSPT_{SP} to dimerized sublattices and a weaker specific heat C(T)C(T) anomaly. The spin susceptibility χ(T)χ(T) and C(T)C(T) are obtained in the thermodynamic limit by a combination of exact diagonalization of small systems with α>0.65α> 0.65 and density matrix renormalization group (DMRG) calculations of systems up to N∼100N \sim 100 spins. Both J1−J2J_1-J_2 and γ>0γ> 0 models account quantitatively for χ(T)χ(T) and C(T)C(T) in the paramagnetic phase of ββ-TeVO4_{4} for T>8T > 8 K, but lower TT indicates a gapped chain instead of a J1−J2J_1-J_2 model as previously thought. The same parameters and TSP=4.6T_{SP} = 4.6 K generate a C(T)/TC(T)/T anomaly that reproduces the anomaly at the 4.64.6 K transition of ββ-TeVO4_{4}, but not the weak χ(T)χ(T) signature.

preprint2020arXiv

Modeling the spin-Peierls transition of spin-1/21/2 chains with correlated states: J1−J2J_1-J_2 model, CuGeO3_3 and TTF-CuS4_4C4_4(CF3_3)4_4

The spin-Peierls transition at TSPT_{SP} of spin-1/21/2 chains with isotropic exchange interactions has previously been modeled as correlated for T>TSPT > T_{SP} and mean field for T<TSPT < T_{SP}. We use correlated states throughout in the J1−J2J_1-J_2 model with antiferromagnetic exchange J1J_1 and J2=αJ1J_2 = αJ_1 between first and second neighbors, respectively, and variable frustration 0≤α≤0.500 \leq α\leq 0.50. The thermodynamic limit is reached at high TT by exact diagonalization of short chains and at low TT by density matrix renormalization group calculations of progressively longer chains. In contrast to mean field results, correlated states of 1D models with linear spin-phonon coupling and a harmonic adiabatic lattice provide an internally consistent description in which the parameter TSPT_{SP} yields both the stiffness and the lattice dimerization δ(T)δ(T). The relation between TSPT_{SP} and Δ(δ,α)Δ(δ,α), the T=0T = 0 gap induced by dimerization, depends strongly on αα and deviates from the BCS gap relation that holds in uncorrelated spin chains. Correlated states account quantitatively for the magnetic susceptibility of TTF-CuS4_4C4_4(CF3_3)4_4 crystals (J1=79J_1 = 79 K, α=0α= 0, TSP=12T_{SP} = 12 K) and CuGeO3_3 crystals (J1=160J_1 = 160 K, α=0.35α= 0.35, TSP=14T_{SP} = 14 K). The same parameters describe the specific heat anomaly of CuGeO3_3 and inelastic neutron scattering. Modeling the spin-Peierls transition with correlated states exploits the fact that δ(0)δ(0) limits the range of spin correlations at T=0T = 0 while T>0T > 0 limits the range at δ=0δ= 0.

preprint2016arXiv

Boundary-induced spin density waves in linear Heisenberg antiferromagnetic spin chains with S≥1\mathbf{S \ge 1}

Linear Heisenberg antiferromagnets (HAFs) are chains of spin-SS sites with isotropic exchange JJ between neighbors. Open and periodic boundary conditions return the same ground state energy in the thermodynamic limit, but not the same spin SGS_G when S≥1S \ge 1. The ground state of open chains of N spins has SG=0S_G = 0 or SS, respectively, for even or odd N. Density matrix renormalization group (DMRG) calculations with different algorithms for even and odd N are presented up to N = 500 for the energy and spin densities ρ(r,N)ρ(r,N) of edge states in HAFs with S=1S = 1, 3/2 and 2. The edge states are boundary-induced spin density waves (BI-SDWs) with ρ(r,N)∝(−1)r−1ρ(r,N)\propto(-1)^{r-1} for r=1,2,…Nr=1,2,\ldots N. The SDWs are in phase when N is odd, out of phase when N is even, and have finite excitation energy Γ(N)Γ(N) that decreases exponentially with N for integer SS and faster than 1/N for half integer SS. The spin densities and excitation energy are quantitatively modeled for integer SS chains longer than 5ξ5 ξ spins by two parameters, the correlation length ξξ and the SDW amplitude, with ξ=6.048ξ= 6.048 for S=1S = 1 and 49.0 for S=2S = 2. The BI-SDWs of S=3/2S = 3/2 chains are not localized and are qualitatively different for even and odd N. Exchange between the ends for odd N is mediated by a delocalized effective spin in the middle that increases ∣Γ(N)∣|Γ(N)| and weakens the size dependence. The nonlinear sigma model (NLσσM) has been applied the HAFs, primarily to S=1S = 1 with even N, to discuss spin densities and exchange between localized states at the ends as Γ(N)∝(−1)Nexp⁡(−N/ξ)Γ(N) \propto (-1)^N \exp(-N/ξ)...

preprint2016arXiv

Numerical study of incommensurate and decoupled phases of spin-1/2 chains with isotropic exchange J1, J2 between first and second neighbors

The spin-1/2 chain with isotropic exchange J1, J2 > 0 between first and second neighbors is frustrated for either sign of J1 and has a singlet ground state (GS) for J1/J2 ≥−4\ge - 4. Its rich quantum phase diagram supports gapless, gapped, commensurate (C), incommensurate (IC) and other phases. Critical points J1/J2 are evaluated using exact diagonalization (ED) and density matrix renormalization group (DMRG) calculations. The wave vector qGq_G of spin correlations is related to GS degeneracy and obtained as the peak of the spin structure factor S(q)S(q). Variable qGq_G indicates IC phases in two J1/J2J1/J2 intervals, [ -4, -1.24] and [0.44, 2], and a C-IC point at J1/J2 = 2. The decoupled C phase in [-1.24, 0.44] has constant qG=π/2q_G = π/2, nondegenerate GS, and a lowest triplet state with broken spin density on sublattices of odd and even numbered sites. The lowest triplet and singlet excitations, EmE_m and EσE_σ, are degenerate in finite systems at specific frustration J1/J2J1/J2. Level crossing extrapolates in the thermodynamic limit to the same critical points as qGq_G. The S(q)S(q) peak diverges at qG=πq_G = π in the gapless phase with J1/J2>4.148J1/J2 > 4.148 and quasi-long-range order (QLRO({ππ})). S(q)S(q) diverges at ±π/2\pm π/2 in the decoupled phase with QLRO({ππ}/2), but is finite in gapped phases with finite range correlations. Numerical results and field theory agree at small J2/J1J2/J1 but disagree for the decoupled phase with weak exchange J1J1 between sublattices. Two related models are summarized: one has an exact gapless decoupled phase with QLRO({ππ}/2) and no IC phases; the other has a single IC phase without a decoupled phase in between.

preprint2015arXiv

Level crossing, spin structure factor and quantum phases of the frustrated spin-1/2 chain with first and second neighbor exchange

The spin-1/2 chain with isotropic Heisenberg exchange J1J_1, J2>0J_2 > 0 between first and second neighbors is frustrated for either sign of J1. Its quantum phase diagram has critical points at fixed J1/J2J_1/J_2 between gapless phases with nondegenerate ground state (GS) and quasi-long-range order (QLRO) and gapped phases with doubly degenerate GS and spin correlation functions of finite range. In finite chains, exact diagonalization (ED) estimates critical points as level crossing of excited states. GS spin correlations enter in the spin structure factor S(q)S(q) that diverges at wave vector qmq_m in QLRO(qmq_m) phases with periodicity 2π/qm2π/q_m but remains finite in gapped phases. S(qm)S(q_m) is evaluated using ED and density matrix renormalization group (DMRG) calculations. Level crossing and the magnitude of S(qm)S(q_m) are independent and complementary probes of quantum phases, based respectively on excited and ground states. Both indicate a gapless QLRO(π/2π/2) phase between −1.2<J1/∣J2∣<0.45-1.2 < J_1/|J_2| < 0.45. Numerical results and field theory agree well for quantum critical points at small frustration J2J_2 but disagree in the sector of weak exchange J1J_1 between Heisenberg antiferromagnetic chains on sublattices of odd and even-numbered sites.

preprint2012arXiv

Spin parity and broken symmetry in finite spin-1/2 chains with frustrated exchange: quantum transition from high to low spin

Exact diagonalization of finite spin-1/2 chains with periodic boundary conditions is applied to the ground state (gs) of chains with ferromagnetic (F) exchange J1<0J_1 < 0between first neighbors, antiferromagnetic (AF) exchange J2=αJ1>0J_2 = αJ_1 > 0between second neighbors, and axial anisotropy 0≤Δ≤10 \le Δ\le 1. In zero field, the gs is in the Sz=0S_z = 0 sector for the relevant parameters and is doubly degenerate at multiple points γm=(αm,Δm)γ_m = (α_m, Δ_m) in the αα, ΔΔ plane. Degeneracy under inversion at sites or spin parity or both leads, respectively, to a bond order wave (BOW), to staggered magnetization or to vector chiral (VC) order. Exact results up to N=28N = 28 spins directly yield order parameters and spin correlation functions whose weak N dependencies allow inferences about infinite chains. The high-spin gs at J2=0J_2 = 0 changes discontinuously at γ1=(−1/4,1)γ_1 = (-1/4, 1) to a singlet in the isotropic (Δ=1Δ= 1) chain. The transition from high to low spin S(α,Δ)S(α, Δ) is continuous for Δ<ΔB=0.95±0.01 Δ< Δ_B = 0.95 \pm 0.01 on the degeneracy line α1(Δ)α_1(Δ). The gs has staggered magnetization between ΔA=0.72Δ_A = 0.72 and ΔBΔ_B, and a BOW for Δ<ΔAΔ< Δ_A. When both inversion and spin parity are reversed at γmγ_m, the correlation functions C(p)C(p) for spins separated by pp sites are identical. C(p)C(p) minima are shifted by π/2π/2 from the minima of VC order parameters at separation pp, consistent with right and left-handed helices along the z axis and spins in the xy plane. Degenerate gs of finite chains are related to quantum phase diagrams of extended αα, ΔΔ chains, with good agreement for order parameters along the line α1(Δ)α_1(Δ).

preprint2011arXiv

Surface dipole of F4TCNQ films: Collective charge transfer and dipole-dipole repulsion in submonolayers

A charge transfer (CT) model is introduced for strong acceptors like A = F4TCNQ that are ionized on surfaces at low coverage θ. Each A forms a CT dimer with a surface state S. Dipole-dipole repulsion grows as θ^(3/2) up to θ = 1 in a full monolayer. Electron transfer ρ(θ) within dimers is found self-consistently and decreases with increasing coverage. The surface dipole and work function shift ΔΦ are proportional to θρ(θ). The CT model has implications for photoemission and accounts for ΔΦ(d) profiles of F4TCNQ films of thickness d on Cu(1,1,1) or on hydrogenated diamond (1,0,0). The CT model also describes ΔΦ(θ) of organic donors on metals and is contrasted to previous treatments of dipoles on a surface.

preprint2010arXiv

Ionization potentials of crystalline organic thin films: Position dependence due to molecular shape and charge redistribution

In addition to electronic polarization or charge redistribution, the shape of neutral conjugated molecules yields position-dependent ionization potentials and electron affinities in organic thin films. Self-consistent I(n) and A(n) are computed in each layer n of 10-layer films of prototypical organics on a metal. The depth dependence of I(n) is discussed at surfaces of anthracene, C60 and PTCDA. The shape contribution can be substantial, up to 0.5 eV, and comes primarily from charge-quadrupole interactions.

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.