
We study the quantum phases driven by interaction in a semimetal with a
quadratic band touching at the Fermi level. By combining the density matrix
renormalization group (DMRG), analytical power expanded Gibbs potential method,
and the weak coupling renormalization group, we study a spinless fermion system
on a checkerboard lattice at halffilling, which has a quadratic band touching
in the absence of interaction. In the presence of strong nearestneighbor
($V_1$) and nextnearestneighbor ($V_2$) interactions, we identify a site
nematic insulator phase, a stripe insulator phase, and a phase separation
region, in agreement with the phase diagram obtained analytically in the strong
coupling limit (i.e. in the absence of fermion hopping). In the intermediate
interaction regime, we establish a quantum anomalous Hall phase in the DMRG as
evidenced by the spontaneous timereversal symmetry breaking and the appearance
of a quantized Chern number $C = 1$. For weak interaction, we utilize the power
expanded Gibbs potential method that treats $V_1$ and $V_2$ on equal footing,
as well as the weak coupling renormalization group. Our analytical results
reveal that not only the repulsive $V_1$ interaction, but also the $V_2$
interaction (both repulsive and attractive), can drive the quantum anomalous
Hall phase. We also determine the phase boundary in the $V_1$$V_2$ plane that
separates the semimetal from the quantum anomalous Hall state. Finally, we show
that the nematic semimetal, which was proposed for $V_2 \gg V_1$ at weak
coupling in a previous study, is absent, and the quantum anomalous Hall state
is the only weak coupling instability of the spinless quadratic band touching
semimetal.

We study dynamical spin structure factor (DSSF) of S = 1/2 Heisenberg model
on the kagome lattice (KAFM) by means of densitymatrix renormalization group.
By comparison with the welldefined magnetic ordered state and chiral spin
liquid sitting nearby in the phase diagram, the KAFM with the nearestneighbor
interaction shows distinct dynamical response behaviors. First of all, the DSSF
displays important spectral intensity predominantly at low frequency region
around Q = M point in momentum space, and shows a broad spectral distribution
at high frequency region for momenta along the boundary of the extended
Brillouin zone. Secondly, spinon excitation spectrum is identified from
momentum and energy resolved DSSF, which shows critical behavior with much
reduced spectrum intensity comparing to the neighboring chiral spin liquid. By
adding a weak DzyaloshinkiiMoriya interaction, the DSSF demonstrates a strong
sensitivity to the boundary condition smore consistent with a gapless spin
liquid. These results capture the main observations in the inelastic neutron
scattering measurements of herbertsmithite, and indicate the spin liquid nature
of the ground state with fractionalized spinon excitations. By following the
DSSF crossing the quantum phase transition between the CSL and the magnetical
ordered phase, we identify the spinon condensation driving the quantum phase
transition.

We explore the frustrated spin$1/2$ Heisenberg model on the star lattice
with antiferromagnetic (AF) couplings inside each triangle and ferromagnetic
(FM) intertriangle couplings ($J_e<0$), and calculate its magnetic and
thermodynamic properties. We show that the FM couplings do not sabotage the
magnetic disordering of the ground state due to the frustration from the AF
interactions inside each triangle, but trigger a fully gapped
inversionsymmetrybreaking trimerized valence bond crystal (TVBC) with
emergent spin1 degrees of freedom. We discover that with strengthening $J_e$,
the system scales exponentially, either with or without a magnetic field $h$:
the order parameter, the five critical fields that separate the $J_e$$h$
groundstate phase diagram into six phases, and the excitation gap obtained by
lowtemperature specific heat, all depend exponentially on $J_e$. We calculate
the temperature dependence of the specific heat, which can be directly compared
with future experiments.

We study the spin$1$ Heisenberg model on the square lattice with the
antiferromagnetic nearestneighbor $J_1$ and the nextnearestneighbor $J_2$
couplings by using the infinite projected entangledpair state (iPEPS) ansatz
and density matrix renormalization group (DMRG) calculation. The iPEPS
simulation, which studies the model directly in the thermodynamic limit, finds
a crossing of the ground state from the N\'eel magnetic state to the stripe
magnetic state at $J_2/J_1 \simeq 0.549$, showing a direct phase transition. In
the finitesize DMRG calculation on the cylinder geometry up to the cylinder
width $L_y = 10$, we find a very small intermediate regime $\sim 0.005 J_1$
between the two magnetic order phases, which may imply the absent intermediate
phase. Both calculations identify that the stripe order comes with a
firstorder transition at $J_2/J_1 \simeq 0.549$. Our results indicate that
unlike the spin$1/2$ $J_1J_2$ square model, quantum fluctuations in the
spin$1$ model may be not strong enough to stabilize an intermediate
nonmagnetic phase.

We study the effects of bond randomness in the spin$1/2$ $J_1J_2$
triangular Heisenberg model using exact diagonalization and density matrix
renormalization group. With increasing bond randomness, we identify a
randomness induced disordered (RID) phase without any magnetic order, dimer
order, spin glass order, or valencebond glass order. The finitesize scaling
of gaps suggests the gapless nature of the spin triplet and singlet
excitations, which is supported by the broad continuum of dynamical spin
structure factor. By studying the bipartite entanglement spectrum on cylinder
system, we find the features of the lowlying spectrum in the RID phase, which
may distinguish the RID phase and the $J_1J_2$ spin liquid phase. We discuss
the nature of this spinliquidlike phase and the indication of our results for
understanding triangular spinliquidlike materials.

The exotic normal state of iron chalcogenide superconductor FeSe, which
exhibits vanishing magnetic order and possesses an electronic nematic order,
triggered extensive explorations of its magnetic ground state. To understand
its novel properties, we study the ground state of a highly frustrated spin$1$
system with bilinearbiquadratic interactions using unbiased largescale
density matrix renormalization group. Remarkably, with increasing biquadratic
interactions, we find a paramagnetic phase between N\'eel and stripe magnetic
ordered phases. We identify this phase as a candidate of nematic quantum spin
liquid by the compelling evidences, including vanished spin and quadrupolar
orders, absence of lattice translational symmetry breaking, and a persistent
nonzero lattice nematic order in the thermodynamic limit. The established
quantum phase diagram natually explains the observations of enhanced spin
fluctuations of FeSe in neutron scattering measurement and the phase transition
with increasing pressure. This identified paramagnetic phase provides a new
possibility to understand the novel properties of FeSe.

We study the spin$1/2$ Heisenberg model on the triangular lattice with the
nearestneighbor $J_1 > 0$, the nextnearestneighobr $J_2 > 0$ Heisenberg
interactions, and the additional scalar chiral interaction $J_{\chi}(\vec{S}_i
\times \vec{S}_j) \cdot \vec{S}_k$ for the three spins in all the triangles
using largescale density matrix renormalization group calculation on cylinder
geometry. With increasing $J_2$ ($J_2/J_1 \leq 0.3$) and $J_{\chi}$
($J_{\chi}/J_1 \leq 1.0$) interactions, we establish a quantum phase diagram
with the magnetically ordered $120^{\circ}$ phase, stripe phase, and
noncoplanar tetrahedral phase. In between these magnetic order phases, we find
a chiral spin liquid (CSL) phase, which is identified as a $\nu = 1/2$ bosonic
fractional quantum Hall state with possible spontaneous rotational symmetry
breaking. By switching on the chiral interaction, we find that the previously
identified spin liquid in the $J_1  J_2$ triangular model ($0.08 \lesssim
J_2/J_1 \lesssim 0.15$) shows a phase transition to the CSL phase at very small
$J_{\chi}$. We also compute spin triplet gap in both spin liquid phases, and
our finitesize results suggest large gap in the odd topological sector but
small or vanishing gap in the even sector. We discuss the implications of our
results to the nature of the spin liquid phases.

Iron chalcogenides display a rich variety of electronic orders in their phase
diagram. A particularly enigmatic case is FeTe, a metal which possesses
coexisting hole and electron Fermi surfaces as in the iron pnictides but has a
distinct ($\pi$/2,$\pi/2$) bicollinear antiferromagnetic order in the Fe square
lattice. While localmoment physics has been recognized as essential for
understanding the electronic order, it has been a longstanding challenge to
understand how the bicollinear antiferromagnetic ground state emerges in a
proper quantum spin model. We show here that a bilinearbiquadratic spin$1$
model on a square lattice with nonzero ringexchange interactions exhibits the
bicollinear antiferromagnetic order over an extended parameter space in its
phase diagram. Our work shows that frustrated magnetism in the quantum spin
model provides a unified description of the electronic orders in the iron
chalcogenides and iron pnictides.

By using Gutzwiller projected fermionic wave functions and variational Monte
Carlo technique, we study the spin$1/2$ Heisenberg model with the
firstneighbor ($J_1$), secondneighbor ($J_2$), and additional scalar chiral
interaction $J_{\chi}{\bf S}_i \cdot ({\bf S}_j \times {\bf S}_k)$ on the
triangular lattice. In the nonmagnetic phase of the $J_1J_2$ triangular model
with $0.08 \lesssim J_2/J_1 \lesssim 0.16$, recent densitymatrix
renormalization group (DMRG) studies [Zhu and White, Phys. Rev. B {\bf 92},
041105 (2015); Hu, Gong, Zhu, and Sheng, Phys. Rev. B {\bf 92}, 140403 (2015)]
find a possible gapped spin liquid with the signal of a competition between a
chiral and a $Z_2$ spin liquid. Motivated by the DMRG results, we consider the
chiral interaction $J_{\chi}{\bf S}_i \cdot ({\bf S}_j \times {\bf S}_k)$ as a
pertubation for this nonmagnetic phase. We find that with growing $J_{\chi}$,
the gapless U(1) Dirac spin liquid, which has the best variational energy for
$J_{\chi}=0$, exhibits the energy instability towards a gapped spin liquid with
nontrivial magnetic fluxes and nonzero chiral order. We calculate topological
Chern number and groundstate degeneracy, both of which identify this flux
state as the chiral spin liquid with fractionalized Chern number $C=1/2$ and
twofold topological degeneracy. Our results indicate a positive direction to
stabilize a chiral spin liquid near the nonmagnetic phase of the $J_1J_2$
triangular model.

We study the ground state phase diagram of the quantum spin$1/2$ Heisenberg
model on the kagom\'{e} lattice with first ($J_1 < 0$), second ($J_2 < 0$),
and thirdneighbor interactions ($J_d > 0$) by means of analytical lowenergy
field theory and numerical densitymatrix renormalization group (DMRG) studies.
The results offer a consistent picture of the $J_d$dominant regime in terms of
three sets of spin chains weakly coupled by the ferromagnetic interchain
interactions $J_{1,2}$. When either $J_1$ or $J_2$ is dominant, the model is
found to support one of two cuboctohedral phases, cuboc1 and cuboc2. These
cuboc states host noncoplanar longranged magnetic order and possess finite
scalar spin chirality. However, in the compensated regime $J_1 \simeq J_2$, a
valence bond crystal phase emerges between the two cuboc phases. We find
excellent agreement between an analytical theory based on coupled spin chains
and unbiased DMRG calculations, including at a very detailed level of
comparison of the structure of the valence bond crystal state. To our
knowledge, this is the first such comprehensive understanding of a highly
frustrated twodimensional (2d) quantum antiferromagnet. We find no evidence of
either the onedimensional (1d) gapless spin liquid or the chiral spin liquids,
which were previously suggested by parton mean field theories.

Unconventional superconductivity often emerges near antiferromagnetic order.
This paradigm applies to many ironbased superconductors, but the iron
chalcogenide FeSe represents a major puzzle. Its normal state has a nematic
order (spontaneously broken lattice rotational symmetry) which, in contrast to
the iron pnictides, is not accompanied by an antiferromagnetic order. In one
proposal, such properties of FeSe originate from the orbital degrees of freedom
associated with the Fe's 3d electrons. In a competing proposition, spin physics
is central to FeSe but magnetic frustration has replaced antiferromagnetism
with a more exotic spin state. Here, we show that a frustrated spin1
bilinearbiquadratic model on the square lattice displays quantum transitions
between nonmagnetic and antiferromagnetic phases, each of which is also
nematic. Our largescale density matrix renormalization group calculations show
that tuning the ratio of competing interactions in the model induces a
transition from a $(\pi,0)$ antiferroquadrupolar order to a $(\pi,0)$
antiferromagnetic state, either directly or through a $(\pi/2,\pi)$
antiferromagnetic order. Our findings explain the recent dramatic experimental
observations of an orthorhombic antiferromagnetic order in the pressurized
FeSe, and suggest that superconductivity in a wide range of ironbased
materials has a common origin.

By using the variational Monte Carlo technique, we study the spin$1/2$ XXZ
antiferromagnetic model (with easyplane anisotropy) on the kagome lattice. A
class of Gutzwiller projected fermionic states with a spin Jastrow factor is
considered to describe either spin liquids (with $U(1)$ or $Z_2$ symmetry) or
magnetically ordered phases (with ${\bf q}=(0,0)$ or ${\bf q}=(4\pi/3,0)$). We
find that the magnetic states are not stable in the thermodynamic limit.
Moreover, there is no energy gain to break the gauge symmetry from $U(1)$ to
$Z_2$ within the spinliquid states, as previously found in the Heisenberg
model. The best variational wave function is therefore the $U(1)$ Dirac state,
supplemented by the spin Jastrow factor. Furthermore, a vanishing $S=2$ spin
gap is obtained at the variational level, in the whole regime from the $XY$ to
the Heisenberg model.

Strongly correlated systems with geometric frustrations can host the emergent
phases of matter with unconventional properties. Here, we study the spin $S =
1$ Heisenberg model on the honeycomb lattice with the antiferromagnetic first
($J_1$) and secondneighbor ($J_2$) interactions ($0.0 \leq J_2/J_1 \leq 0.5$)
by means of density matrix renormalization group (DMRG). In the parameter
regime $J_2/J_1 \lesssim 0.27$, the system sustains a N\'{e}el
antiferromagnetic phase. At the large $J_2$ side $J_2/J_1 \gtrsim 0.32$, a
stripe antiferromagnetic phase is found. Between the two magnetic ordered
phases $0.27 \lesssim J_2/J_1 \lesssim 0.32$, we find a \textit{nonmagnetic}
intermediate region with a plaquette valencebond order. Although our
calculations are limited within $6$ unitcell width on cylinder, we present
evidence that this plaquette state could be a strong candidate for this
nonmagnetic region in the thermodynamic limit. We also briefly discuss the
nature of the quantum phase transitions in the system. We gain further insight
of the nonmagnetic phases in the spin$1$ system by comparing its phase
diagram with the spin$1/2$ system.

We study the spin$1/2$ Heisenberg model on the triangular lattice with the
antiferromagnetic first ($J_1$) and second ($J_2$) nearestneighbor
interactions using density matrix renormalization group. By studying the spin
correlation function, we find a $120^{\circ}$ magnetic order phase for $J_2
\lesssim 0.07 J_1$ and a stripe antiferromagnetic phase for $J_2 \gtrsim 0.15
J_1$. Between these two phases, we identify a spin liquid region characterized
by the exponential decaying spin and dimer correlations, as well as the large
spin singlet and triplet excitation gaps on finitesize systems. We find two
near degenerating ground states with distinct properties in two sectors, which
indicates more than one spin liquid candidates in this region. While the sector
with spinon is found to respect the time reversal symmetry, the even sector
without a spinon breaks such a symmetry for finitesize systems. Furthermore,
we detect the signature of the fractionalization by following the evolution of
different ground states with inserting spin flux into the cylinder system.
Moreover, by tuning the anisotropic bond coupling, we explore the nature of the
spin liquid phase and find the optimal parameter region for the gapped $Z_2$
spin liquid.

We study the quantum phase diagram of the spin$1/2$ Heisenberg model on the
kagom\'e lattice with first, second, and thirdneighbor interactions $J_1$,
$J_2$, and $J_3$ by means of density matrix renormalization group. For small
$J_2$ and $J_3$, this model sustains a timereversal invariant quantum spin
liquid phase. With increasing $J_2$ and $J_3$, we find in addition a $q=(0,0)$
N\'{e}el phase, a chiral spin liquid phase, a valencebond crystal phase, and a
complex noncoplanar magnetically ordered state with spins forming the vertices
of a cuboctahedron known as a cuboc1 phase. Both the chiral spin liquid and
cuboc1 phase break time reversal symmetry in the sense of spontaneous scalar
spin chirality. We show that the chiralities in the chiral spin liquid and
cuboc1 are distinct, and that these two states are separated by a strong first
order phase transition. The transitions from the chiral spin liquid to both the
$q=(0,0)$ phase and to timereversal symmetric spin liquid, however, are
consistent with continuous quantum phase transitions.

We study the spin1/2 Heisenberg model on the square lattice with first and
secondneighbor antiferromagnetic interactions J1 and J2, which possesses a
nonmagnetic region that has been debated for many years and might realize the
interesting Z2 spin liquid. We use the density matrix renormalization group
approach with explicit implementation of SU(2) spin rotation symmetry and study
the model accurately on open cylinders with different boundary conditions. With
increasing J2, we find a Neel phase, a plaquette valencebond (PVB) phase with
a finite spin gap, and a possible spin liquid in a small region of J2 between
these two phases. From the finitesize scaling of the magnetic order parameter,
we estimate that the Neel order vanishes at J2/J1~0.44. For 0.5<J2/J1<0.61, we
find dimer correlations and PVB textures whose decay lengths grow strongly with
increasing system width, consistent with a longrange PVB order in the
twodimensional limit. The dimerdimer correlations reveal the swave character
of the PVB order. For 0.44<J2/J1<0.5, spin order, dimer order, and spin gap are
small on finitesize systems and appear to scale to zero with increasing system
width, which is consistent with a possible gapless SL or a nearcritical
behavior. We compare and contrast our results with earlier numerical studies.

The fractional quantum Hall effect (FQHE) realized in twodimensional
electron systems under a magnetic field is one of the most remarkable
discoveries in condensed matter physics. Interestingly, it has been proposed
that FQHE can also emerge in timereversal invariant spin systems, known as the
chiral spin liquid (CSL) characterized by the topological order and the
emerging of the fractionalized quasiparticles. A CSL can naturally lead to the
exotic superconductivity originating from the condense of anyonic
quasiparticles. Although CSL was highly sought after for more than twenty
years, it had never been found in a spin isotropic Heisenberg model or related
materials. By developing a densitymatrix renormalization group based method
for adiabatically inserting flux, we discover a FQHE in a spin$\frac{1}{2}$
isotropic kagome Heisenberg model. We identify this FQHE state as the
longsought CSL with a uniform chiral order spontaneously breaking time
reversal symmetry, which is uniquely characterized by the halfinteger
quantized topological Chern number protected by a robust excitation gap. The
CSL is found to be at the neighbor of the previously identified $Z_2$ spin
liquid, which may lead to an exotic quantum phase transition between two gapped
topological spin liquids.

We use density matrix renormalization group (DMRG) algorithm to study the
phase diagram of the spin1/2 Heisenberg model on honeycomb lattice with first
(J1) and second (J2) neighbor antiferromagnetic interactions, where a Z2 spin
liquid region has been proposed. By implementing SU(2) symmetry in the DMRG
code, we are able to obtain accurate results for long cylinders with width
slightly over 15 lattice spacings and torus up to the size N=2*6*6. With
increasing J2, we find a Neel phase with vanishing spin gap and a plaquette
valencebond (PVB) phase with nonzero spin gap. By extrapolating the square of
the staggered magnetic moment ms^2 on finitesize cylinders to thermodynamic
limit, we find the Neel order vanishing at J2/J1~0.22. For 0.25<J2/J1<0.35, we
find a possible PVB order, which shows a fast growing PVB decay length with
increasing system width. For 0.22<J2/J1<0.25, both spin and dimer orders are
weak for all systems we have studied, which is consistent with a possible spin
liquid phase or a deconfined quantum critical point. We present calculations of
the topological entanglement entropy, compare the DMRG results with the
variational Monte Carlo, and discuss possible scenarios in the thermodynamic
limit for this region.

Motivated by the crystal structures of [(CuCl2tachH)3Cl]Cl2 and Ca3Co2O6, we
develop a lowenergy effective theory using the bosonization technique for a
spin1/2 frustrated threeleg spin tube with trigonal prism units in two limit
cases. The features obtained with the effective theory are numerically
elucidated by the density matrix renormalization group method. Three different
quantum phases in the ground state of the system, say, one gapped dimerized
phase and two distinct gapless phases, are identified, where the two gapless
phases are found to have the conformal central charge c=1 and 3/2,
respectively. Spin gaps, spin and dimer correlation functions, and the
entanglement entropy are obtained. In particular, it is disclosed that the
critical phase with c=3/2 is the consequence of spin frustrations, which might
belong to the SU(2)k=2 WessZuminoWittenNovikov universality class, and is
induced by the twist term in the bosonized Hamiltonian density.

A linearized tensor renormalization group (LTRG) algorithm is proposed to
calculate the thermodynamic properties of onedimensional quantum lattice
models, that is incorporated with the infinite timeevolving block decimation
technique, and allows for treating directly the twodimensional transfermatrix
tensor network. To illustrate its feasibility, the thermodynamic quantities of
the quantum XY spin chain are calculated accurately by the LTRG, and the
precision is shown to be comparable with (even better than) the transfer matrix
renormalization group (TMRG) method. Unlike the TMRG scheme that can only deal
with the infinite chains, the present LTRG algorithm could treat both finite
and infinite systems, and may be readily extended to boson and fermion quantum
lattice models.

The twodimensional Ising model on a distorted Kagom\'{e} lattice is studied
by means of exact solutions and the tensor renormalisation group (TRG) method.
The zerofield phase diagrams are obtained, where three phases such as
ferromagnetic, ferrimagnetic and paramagnetic phases, along with the
secondorder phase transitions, have been identified. The TRG results are quite
accurate and reliable in comparison to the exact solutions. In a magnetic
field, the magnetization ($m$), susceptibility and specific heat are studied by
the TRG algorithm, where the $m=1/3$ plateaux are observed in the magnetization
curves for some couplings. The experimental data of susceptibility for the
complex Co(N$_3$)$_2$(bpg)$\cdot$ DMF$_{4/3}$ are fitted with the TRG results,
giving the couplings of the complex $J=22K$ and $J'=33K$.

The magnetic and thermodynamic properties of a new ferrimagnetic decorated
spin(1/2,1) Heisenberg chain with spin1 pendant spins are investigated for
three cases: (A) J1,J2>0; (B) J1>0, J2<0; and (C) J1<0, J2>0, where J1 and J2
are the exchange couplings between spins in the chain and along the rung,
respectively. The lowlying and magnetic properties are explored jointly by the
realspace renormalization group, spin wave, and densitymatrix renormalization
group methods, while the transfermatrix renormalization group method is
invoked to study the thermodynamics. It is found that the magnon spectra
consist of a gapless and two gapped branches. Two branches in case (C) have
intersections. The coupling dependence of lowenergy gaps are analyzed. In a
magnetic field, an m=3/2 (m is the magnetization per unit cell) plateau is
observed for case (A), while two plateaux at m=1/2 and 3/2 are observed for
cases (B) and (C). Between the two plateaux in cases (B) and (C), the
sublattice magnetizations for the spins coupled by ferromagnetic interactions
have novel decreasing regions with increasing the magnetic field. At finite
temperature, the zerofield susceptibility temperature product chi*T and
specific heat exhibit distinct exotic features with varying the couplings and
temperature for different cases. chi*T is found to converge as T approaches
zero, which is different from the divergent behavior in the spin(1/2,1)
mixedspin chain without pendants. The observed thermodynamic behaviors are
also discussed with the help of their lowlying excitations.

The spin1/2 Heisenberg antiferromagnet on the distorted honeycomb (DHC)
lattice is studied by means of the tensor renormalization group method. It is
unveiled that the system has a quantum phase transition of secondorder between
the gapped quantum dimer phase and a collinear Neel phase at the critical point
of coupling ratio \alpha_{c} = 0.54, where the quantum critical exponents \nu =
0.69(2) and \gamma = 1.363(8) are obtained. The quantum criticality is found to
fall into the O(3) universality class. A groundstate phase diagram in the
fieldcoupling ratio plane is proposed, where the phases such as the dimer,
semiclassical Neel, and polarized phases are identified. A link between the
present spin system to the boson Hubbard model on the DHC lattice is also
discussed.

The magnetic properties of a trigonal prism unit of the spin2 frustrated
compound Ca3Co2O6 are studied by means of the densitymatrix renormalization
group method. A magnetization plateau at $ms/3$ ($ms$ is the saturation
magnetization) with ferrimagnetic structure is observed. By fitting the
experimental data of magnetic curve, an estimation of the couplings gives
J1=26.84K, J_{2}=0.39K, and J_{3}=0.52K. The local magnetic moments are
unveiled to exhibit an incommensurate sinusoidally modulation along the three
chains of the trigonal prism, which gives a strong theoretical support to the
experimentally observed incommensurate partially disordered antiferromagnetic
state for Ca3Co2O6. The present result suggests that the modulation indeed
originates from the competition of antiferromagnetic and ferromagnetic
couplings.

The thermal pairwise entanglement (TE) of the S=1/2 XY chain in a transverse
magnetic field is exactly resolved by means of the JordanWigner
transformation. It is found that the TE vanishes at a common temperature
Tc~0.4843J, which is irrelevant to the field. A thermal quantity is proposed to
witness the entangled state. Furthermore, the TE of the S=1/2
antiferromagneticferromagnetic (AFF) Heisenberg chain is studied by the
transfermatrix renormalization group method.The TEs of the spins coupled by AF
and F interactions are found to behave distinctively. The vanishing temperature
of the fieldinduced TE of the spins coupled by F interactions is observed to
change with the magnetic field. The results are further confirmed and analyzed
by the meanfield theory.