We study the one-dimensional S = 1/2 Heisenberg model with an uniform and a staggered magnetic fields, using the dynamical density-matrix renormalization group (DDMRG) technique. The DDMRG enables us to investigate the dynamical properties of chain with lengths up to a few hundreds, and the results are numerically exact in the same sense as 'exact diagonalization' results are. Thus, we can analyze the low-energy spectrum almost in the thermodynamic limit. In this work, we calculate the dynamical spin structure factor and demonstrate the performance of the DDMRG method applying the open-end boundary conditions as well as the periodic boundary conditions.
The magnetic and thermodynamic properties of spin-
1∕2 Heisenberg diamond chains are investigated in three different cases: (a)
J1,
J2,
J3>0 (frustrated), (b)
J1,
J3<0,
J2>0 (frustrated), and (c)
J1,
J2>0,
J3<0 (nonfrustrated), where the bond coupling
Ji (i=1,2,3)>0 stands for an antiferromagnetic (AF) interaction, and
<0 for a ferromagnetic (F) interaction. The density-matrix renormalization-group (DMRG) technique is invoked to study the properties of the system in the ground state, while the transfer-matrix renormalization-group (TMRG) technique is applied to explore the thermodynamic properties. The local magnetic moments, spin-correlation functions, and static structure factors are discussed in the ground state for the three cases. It is shown that the static structure factor
S(q) shows peaks at wave vectors $q=a\ensuremath{\pi}∕3$
(a=0,1,2,3,4,5) for different couplings in a zero magnetic field, which, however, in the magnetic fields where the magnetization plateau with
m=1∕6 pertains, exhibits the peaks only at
q=0, $2\ensuremath{\pi}∕3$, and $4\ensuremath{\pi}∕3$, which are found to be couplings independent. The DMRG results of the zero-field static structure factor can be nicely fitted by a linear superposition of six modes, where two fitting equations are proposed. It is observed that the six modes are closely related to the low-lying excitations of the system. At finite temperatures, the magnetization, susceptibility, and specific heat show various behaviors for different couplings. The double-peak structures of the susceptibility, and specific heat against temperature are obtained, where the peak positions and heights are found to depend on the competition of the couplings. It is also uncovered that the
XXZ anisotropy of F and AF couplings leads the system of case (c) to display quite different behaviors. In addition, the experimental data of the susceptibility, specific heat, and magnetization for the compound
Cu3(CO3)2(OH)2 are fairly compared with our TMRG results.
The spin-1/2 Heisenberg antiferromagnet on the frustrated diamond-decorated square lattice is known to feature various zero-field ground-state phases, consisting of extended monomer-dimer and dimer-tetramer ground states as well as a ferrimagnetic regime. Using a combination of analytical arguments, density matrix renormalization group (DMRG), exact diagonalization as well as sign-problem-free quantum Monte Carlo (QMC) calculations, we investigate the properties of this system and the related Lieb lattice in the presence of a finite magnetic field, addressing both the ground-state phase diagram as well as several thermodynamic properties. In addition to the zero-field ground states, we find at high magnetic field a spin-canted phase with a continuously rising magnetization for increasing magnetic field strength as well as the fully polarized paramagnetic phase. At intermediate field strength, we identify a first-order quantum phase transition line between the ferrimagnetic and the monomer-dimer regime. This first-order line extends to finite temperatures, terminating in a line of critical points that belong to the universality class of the two-dimensional Ising model.
This thesis investigates the ground-state properties and quantum phase diagrams of low-dimensional spin systems, with a focus on spin–1/2 ladder chains and related models. These systems provide a fertile ground for exploring the interplay between magnetic frus- tration, quantum correlations, topology, and external fields. By combining advanced nu- merical methods, primarily the density matrix renormalization group (DMRG), with ana- lytical approaches, we analyze how different coupling geometries and magnetic fields give rise to distinct quantum phases, including magnetization plateaus, symmetry-protected and topological phases, and first-order transitions. The study covers both uniform and alternating ladders, as well as coupled ladder systems that interpolate between one- and two-dimensional behavior. In particular, we examine how frustration influences the stabil- ity of collective spin states, and how interladder couplings can drive transitions between resonating valence bond (RVB) states, Haldane phases, and trivial phases hosting collec- tive spin excitations. Although approximate analytical models, such as hard-core boson mappings and spin-wave theory, capture some qualitative aspects of these systems, our results highlight their limitations in accurately describing critical behavior and phase boundaries, reinforcing the importance of high-precision numerical analysis.
In this thesis, we deal with 1-dimensional anisotropic spin models more specifically with the chain and ladder geometries. These models are fundamental to the study of condensed matter theory and quantum magnetism because of their simplicity. Research into the chain model dates back to 1930s with Hans Bethe. One method to study these systems is through numerical analysis including renormalization group techniques. These allow for the diagonalization of the Hamiltonian without sacrificing computation to a large virtual Hilbert space. In this thesis, we present results for both the spin chain and spin ladder geometries in open boundary conditions. This work uses the density matrix renormalization group technique to calculate system energy. Initially we will present a study into the spin energy gap of the spin-\sfrac{1}{2} anisotropic (XXZ) Heisenberg chain in open boundary conditions (OBCs). The energy gap is shown to be reduced to half in the ground sector when compared to periodic boundary conditions (PBCs) due to edge effects. Secondly, a full phase diagram for the spin-\sfrac{1}{2} anisotropic Heisenberg ladder is presented which shows the emergence of a rich schematic with a variety of phases. Lastly, weak coupling limit maps are compared to quantum field theory phase transition predictions. This research was motivated by the lack of comprehensive phase diagram results for ladders and by real materials being investigated at the University of Kent.
This thesis is concerned with the numerical study of one-dimensional (1D) spin-1/2 quantum magnets and related method development. Its focus is on the calculation of dynamical spin correlation functions both at zero and finite temperature. This is motivated by the accessibility of dynamical quantities in experiments such as inelastic neutron scattering (INS) and electron spin resonance (ESR). The numerical methods used in this thesis are based on extensions of the density-matrix renormalization group (DMRG) and are formulated in the framework of matrix product states (MPS). While zero-temperature dynamical correlation functions are computed with existing MPS frequency-domain methods, an MPS frequency-domain approach for their calculation at finite temperature is developed in this thesis. The new method combines the Liouville-space formulation of the dynamics with a moment expansion of the dynamical correlation function. The majority of the results are obtained via MPS-based Chebyshev expansions. These numerical techniques are applied to two different model systems describing real materials. The first one is the material copper pyrimidine dinitrate (Cu-PM) which is modeled by a 1D spin-1/2 Heisenberg antiferromagnet with Dzyaloshinskii-Moriya interactions. The spin dynamics of this model is studied in an applied magnetic field and compared to ESR experiments. Zero-temperature calculations for momentum- and frequency-resolved dynamical quantities give direct access to the intensity of the elementary excitations and go beyond the low-energy description by the quantum sine-Gordon model. Thus, a deviation from the Lorentz invariant dispersion for the single-soliton resonance is found. The presence of the strongest boundary bound state previously predicted from a boundary sine-Gordon field theory is confirmed, while composite boundary-bulk excitations have too low intensities to be found in the numerical results. At finite temperature, there is a temperature-induced crossover of the soliton. Moreover, additional temperature effects such as interbreather transitions emerge, which is confirmed by accompanying ESR experiments on Cu-PM over a wide range of the applied field strength. The second system studied in this thesis is the compound BaCu
2V
2O
8. It is shown that the magnetic properties of this quasi-1D material can be described by a strongly alternating antiferromagnetic-ferromagnetic spin-1/2 Heisenberg chain. As found for other dimerized systems, the strong correlations in BaCu
2V
2O
8 persist even at elevated temperatures. Moreover, these correlations lead to an asymmetric lineshape broadening of the magnetic excitations at finite temperature. Upon raising the temperature, an increasingly asymmetric lineshape is observed in high-resolution INS experiments for BaCu
2V
2O
8. In this thesis, the lineshape is calculated by the MPS finite-temperature method developed as a part of this work. Comparing these results to the INS data and the lineshape obtained by a diagrammatic approach, excellent agreement is found over a broad temperature range. This demonstrates that coherent quantum behavior persists at elevated temperatures in BaCu
2V
2O
8 and that it can be predicted quantitatively.