\section{Introduction}
\label{sec:intro}

\begin{figure}[th]
\centering
\includegraphics[width=1.0\linewidth]{images/teaser_compressed.pdf}
\caption{Our TediGAN is the first method that unifies text-guided image generation and manipulation into one same framework, leading to naturally continuous operations from generation to manipulation (a), and inherently supports image synthesis with multi-modal inputs (b), such as sketches or semantic labels with or without instance (texts or real images) guidance.
}
\label{fig:teaser}
\end{figure}
How to create or edit an image of the desired content without tedious manual operations is a difficult but meaningful task.
To make image generation and manipulation more readily and user-friendly, recent studies have been focusing on the image synthesis conditioned on a variety of guidance, such as sketch~\cite{ghosh2019isketchnfill,xia2019sketch}, semantic label~\cite{isola2017image,wang2018high}, or textual description~\cite{nam2018tagan,xu2018attngan}.
Despite the success of its label and sketch counterparts, most state-of-the-art text-guided image generation and manipulation methods are only able to produce low-quality images~\cite{reed2016generative,dong2017semantic}.
Those aiming at generating high-quality images from texts typically design a multi-stage architecture and train their models in a progressive manner.
To be more specific, there are usually three stages in the main module, and each stage contains a generator and a discriminator. Three stages are trained at the same time, and progressively generate images of three different scales, \emph{i.e.}, \(\text{64}^2 \to \text{128}^2 \to \text{256}^2\).
The initial image with rough shape and color would be refined to a high-resolution one.
However, the multi-stage training process is time-consuming and cumbersome, making the aforementioned methods unfeasible for higher resolution.
Furthermore, the pretrained text-image matching model they used fails to exploit attribute-level cross-modal information and leads to mismatched attributes when generating images from texts~\cite{xu2018attngan,li2019control,zhu2019dmgan,cheng2020rifegan}, or undesired changes of irrelevant attributes when manipulating images~\cite{dong2017semantic,nam2018tagan,li2020manigan,li2020lightweight}.

\begin{figure}[t]
\centering
\includegraphics[width=0.95\linewidth]{images/mapping.pdf}
\caption{{Projecting Multi-Modal Embedding into the $\mathcal{W}$ Space of StyleGAN.}
Taking visual and linguistic embedding for example, the left illustrates visual-linguistic similarity learning, where the visual embedding $\w^v$ and linguistic embedding $\w^l$ are expected to be close enough.
The right demonstrates text-guided image manipulation.
Given a source image and a text guidance, we first get their embedding $\w^v$ and $\w^l$ in $\mathcal{W}$ space through corresponding encoders. We then perform style mixing for target layers and get the target latent code $\w^t$. The final $\w^{t*}$ is obtained through instance-level optimization. The edited image can be generated from the StyleGAN generator.
}
\label{fig:mapping}
\end{figure}

Recent progress on generative adversarial networks (GANs) has established an entirely different image generation paradigm that achieves phenomenal quality, fidelity, and realism.
StyleGAN~\cite{karras2019style}, one of the most notable GAN frameworks, introduces a novel style-based generator architecture and can produce high-resolution images with unmatched photorealism.
Some recent work~\cite{karras2019style} has demonstrated that the intermediate latent space $\W$ of StyleGAN, inducted from a learned piece-wise continuous mapping, yields less entangled representations and offers more feasible manipulation.
The superior characteristics of $\W$ space appeal to numerous researchers to develop advanced GAN inversion techniques~\cite{xia2021survey,bau2019inverting,abdal2019image2stylegan} to invert real images back into the StyleGAN's latent space and perform meaningful manipulation.
The most popular way~\cite{zhu2020indomain,richardson2020encoding} is to train an additional encoder to map real images into the $\W$ space, which leads to not only faithful reconstruction but also semantically meaningful editing.
Furthermore, it is easy to introduce the hierarchically semantic property of the $\W$ space to any GAN model by simply learning an extra mapping network before a fixed, pretrained StyleGAN generator.
We thoroughly investigated the existing GAN inversion methods, and found all is about how to map images into the latent space of a well-trained GAN model.
The other modalities like texts, however, have not received any attention.

In this paper, for the first time, we propose a GAN inversion technique that can map multi-modal information, \emph{e.g.}, texts, sketches, or labels, into a common latent space of a pretrained StyleGAN.
Based on that, we propose a very simple yet effective method for \textcolor{blue}{Te}xt-guided \textcolor{blue}{di}verse image generation and manipulation via \textcolor{blue}{GAN} (abbreviated \textcolor{blue}{TediGAN}).
Our proposed method introduces three novel modules. The first StyleGAN inversion module learns the inversion where an image encoder can map a real image to the $\W$ space,
while the second visual-linguistic similarity module learns linguistic representations that are consistent with the visual representations by projecting both into a common $\W$ space, as shown in Figure~\ref{fig:mapping}.
The third instance-level optimization module is to preserve the identity after editing, which can precisely manipulate the desired attributes consistent with the texts while faithfully reconstructing the unconcerned ones.
Our proposed method can generate diverse and high-quality results
with a resolution up to $\text{1024}^2$,
and inherently support image synthesis with multi-modal inputs, such as sketches or semantic labels with or without instance (texts or real images) guidance.
Due to the utilization of a pretrained StyleGAN model, our method can provide the lowest effect guarantee, \emph{i.e.}, our method can always produce pleasing results no matter how uncommon the given text or image is.
Furthermore, to fill the gaps in the text-to-image synthesis dataset for faces, we create the Multi-Modal CelebA-HQ dataset to facilitate the research community.
Following the format of the two popular text-to-image synthesis datasets, \emph{i.e.}, CUB~\cite{wah2011caltech} for birds and COCO~\cite{lin2014microsoft} for natural scenes, we create ten unique descriptions for each image in the CelebA-HQ~\cite{karras2017progressive}.
Besides real faces and textual descriptions, the introduced dataset also contains the label map and sketch for the text-guided generation with multi-modal inputs.

In summary, this work has the following contributions:
\begin{itemize}
\item
We propose a unified framework that can generate diverse images given the same input text, or text with image for manipulation, allowing the user to edit the appearance of different attributes interactively.

\item We propose a GAN inversion technique that can map multi-modal information into a common latent space of a pretrained StyleGAN where the instance-level image-text alignment can be learned.
\item We introduce the Multi-Modal CelebA-HQ dataset, consisting of multi-modal face images and corresponding textual descriptions, to facilitate the community.
\end{itemize}

\begin{figure*}[th]
\centering
\begin{minipage}[t]{0.33\textwidth}
\centering
\includegraphics[width=0.95\linewidth]{images/high-res-gene.jpg}\\
\vspace{-5.0pt}\fontsize{7.0pt}{\baselineskip}\selectfont{\textit{a smiling young woman with short blonde hair}}
\end{minipage}
\begin{minipage}[t]{0.33\textwidth}
\centering
\includegraphics[width=0.95\linewidth]{images/high-res-lab.jpg}\\
\vspace{-5.0pt}\fontsize{7.0pt}{\baselineskip}\selectfont{\textit{he is young and wears beard}}
\end{minipage}
\begin{minipage}[t]{0.33\textwidth}
\centering
\includegraphics[width=0.95\linewidth]{images/high-res-skt.jpg}\\
\vspace{-5.0pt}\fontsize{7.0pt}{\baselineskip}\selectfont{\textit{a young woman with long black hair}}
\end{minipage}
\caption{Diverse High-Resolution Results from Text.
Our method can achieve text-guided diverse image generation and manipulation up to an unprecedented resolution at $\text{1024}^2$.}
\label{fig:high-res-gene}
\end{figure*}

\section{Related Work}
\label{sec:related_work}

\paragraph{Text-to-Image Generation.}
There are basically two categories of GAN-based text-to-image generation methods.
The first category produces images from texts directly by one generator and one discriminator.
For example, Reed~\emph{et al.}~\cite{reed2016generative} propose to use conditional GANs to generate plausible images from given text descriptions.
Tao~\emph{et al.}~\cite{tao2020dfgan} propose a simplified backbone that generates high-resolution images directly by Wasserstein distance and fuses the text information into visual feature maps to improve the image quality and text-image consistency.
Despite the plainness and conciseness, the one-stage models produce dissatisfied results in terms of both photo-realism and text-relevance in some cases.
Thus, another thread of research focuses on multi-stage processing.
Zhang~\emph{et al.}~\cite{zhang2017stackgan} stack two GANs to generate high-resolution images from text descriptions through a sketch-refinement process.
They further propose a three-stage architecture~\cite{zhang2018stackgan++} that stacks multiple generators and discriminators, where multi-scale images are generated progressively in a course-to-fine manner.
Xu~\emph{et al.}~\cite{xu2018attngan} improve the work of~\cite{zhang2018stackgan++} from two aspects.
First, they introduce attention mechanisms to explore fine-grained text and image representations.
Second, they propose a Deep Attentional Multimodal Similarity Model (DAMSM) to compute the similarity between the generated image and the sentence.
The subsequent studies basically follow the framework of~\cite{xu2018attngan} and have proposed several variants by introducing different mechanisms like attention~\cite{li2019control} or memory writing gate~\cite{zhu2019dmgan}.
However, the multi-stage frameworks produce results that look like a simple combination of visual attributes from different image scales.

\vspace{-5pt}
\paragraph{Text-Guided Image Manipulation.}
Similar to text-to-image generation, manipulating given images using texts also produces results that contain desired visual attributes.
Differently, the modified results should only change certain parts and preserve text-irrelevant contents of the original images.
For example, Dong~\emph{et al.}~\cite{dong2017semantic} propose an encoder-decoder architecture to modify an image according to a given text.
Nam~\emph{et al.}~\cite{nam2018text} disentangle different visual attributes by introducing a text-adaptive discriminator, which can provide finer training feedback to the generator.
Li \emph{et al.}~\cite{li2020manigan} introduce a multi-stage network with a novel text-image combination module to produce high-quality results.
Similar to text-to-image generation, the text-based image manipulation methods with the best performance are basically based on the multi-stage framework.
Different from all existing methods, we propose a novel framework that unifies text-guided image generation and manipulation methods and can generate high-resolution and diverse images \texttt{directly} without multi-stage processing.

\vspace{-5pt}
\paragraph{Image-Text Matching.}
One key of text-guided image generation or manipulation is to match visual attributes with corresponding words. To do this, current methods usually provide explicit word-level training feedback from the elaborately-designed discriminator~\cite{li2020manigan,li2020lightweight}.
There is also a rich line of work proposed to address a related direction named image-text matching, or visual-semantic alignment, aiming at exploiting the matching relationships and making the corresponding alignments between text and image.
Most of them can be categorized into two-branch deep architecture according to the granularity of representations for both modalities, \emph{i.e.}, global~\cite{kiros2014unifying,mao2014deep,ma2015multimodal} or local~\cite{karpathy2014deep,karpathy2015deep,lee2018stacked} representations.
The first category employs deep neural networks to extract the global features of both modalities, based on which their similarities are measured~\cite{mao2014deep}.
Another thread of work performs instance-level image-text matching~\cite{nam2017dual,lee2018stacked,song2019polysemous}, learning the correspondences between words and image regions~\cite{karpathy2015deep}.

\begin{figure*}[ht]
\centering
\includegraphics[width=1.0\linewidth]{images/comparison_generation_compressed.pdf}
\caption{Comparison of Text-to-Image Generation on Our Multi-modal CelebA-HQ dataset.}
\label{fig:comp_gen}
\end{figure*}

\section{The TediGAN Framework}
\label{sec:method}
We first learn the inversion, \emph{i.e.}, training an image encoder to map the real images to the latent space such that all codes produced by the encoder can be recovered at both the pixel-level and the semantic-level.
We then use the hierarchical characteristic of $\W$ space to learn the text-image matching by mapping the image and text into the same joint embedding space.
To preserve identity in manipulation, we propose an instance-level optimization, involving the trained encoder as a regularization to better reconstruct the pixel values without affecting the semantic property of the inverted code.

\subsection{StyleGAN Inversion Module}
\label{subsec:gan-inversion}
The inversion module aims at training an image encoder that can map a real face image to the latent space of a fixed StyleGAN model pretrained on the FFHQ dataset~\cite{karras2019style}.
The reason we invert a trained GAN model instead of training one from scratch is that, in this way, we can go beyond the limitations of a paired text-image dataset.
The StyleGAN is trained in an unsupervised setting and covers much higher quality and wider diversity, which makes our method able to produce satisfactory edited results with images in the wild.
In order to facilitate subsequent alignment with text attributes, our goal for inversion is not only to reconstruct the input image by pixel values but also to acquire the inverted code that is semantically meaningful and interpretable~\cite{shen2020interpreting,yang2019semantic}.

Before introducing our method, we first briefly establish problem settings and notations.
A GAN model typically consists of a generator $G(\cdot): \Z\rightarrow\X$ to synthesize fake images and a discriminator $D(\cdot)$ to distinguish real data from the synthesized.
In contrast, GAN inversion studies the reverse mapping, which is to find the best latent code $\z^{*}$ by inverting a given image $\x$ to the latent space of a well-trained GAN.
A popular solution is to train an additional encoder $E_v(\cdot): \X\rightarrow\Z$~\cite{zhu2016generative,bau2019seeing} (subscript $v$ means visual).
To be specific, a collection of latent codes $\z^{s}$ are first randomly sampled from a prior distribution, \emph{e.g.}, normal distribution, and fed into $G(\cdot)$ to get the synthesis $\x^{s}$ as the training pairs.
The introduced encoder $E_v(\cdot)$ takes $\x^{s}$ and $\z^{s}$ as inputs and supervisions respectively and is trained with
\begin{align}
  \min_{\Theta_{E_v}}\Loss_{E_v} = ||\z^{s} - E_v(G(\z^{s}))||_2^2,
  \label{eq:conventional-encoder}
\end{align}
where $||\cdot||_2$ denotes the $l_2$ distance and $\Theta_{E_v}$ represents the parameters of the encoder $E_v(\cdot)$.

Despite of its fast inference, the aforementioned procedure simply learns a deterministic model with no regard to whether the codes produced by the encoder align with the semantic knowledge learned by $G(\cdot)$.
The supervision by only reconstructing $\z^{s}$ is not powerful enough to train $E_v(\cdot)$, and $G(\cdot)$ is actually not fully used to guide the training of $E_v(\cdot)$, leading to the incapability of inverting real images.
To solve these problems, we use a totally different strategy to train an encoder for GAN inversion as in~\cite{zhu2020indomain}.
There are two main differences compared with the conventional framework:
(a) the encoder is trained with real images rather than with synthesized images, making it more applicable to real applications;
(b) the reconstruction is at the image space instead of latent space, which provides semantic knowledge and accurate supervision and allows integration of powerful image generation losses such as perceptual loss~\cite{johnson2016perceptual} and LPIPS~\cite{zhang2018unreasonable}.
Hence, the training process can be formulated as
\begin{small}
\begin{align}
  &\begin{aligned}
    \min_{\Theta_{E_v}}\Loss_{E_v} \!=\! ||\x - G(E_v(\x))||_2^2\ &\!+\!\lambda_{1} ||F(\x) \!-\! F(G(E_v(\x)))||_2^2\ \\
    &\!-\!\lambda_{2}\E[D_v(G(E_v(\x)))], \label{eq:encoder}
  \end{aligned} \\
  &\begin{aligned}
    \min_{\Theta_{D_v}}\Loss_{D_v} = \!\E[D_v(G(E_v(\x)))] \!- \!\E[D_v(\x)] \! + \!\frac{\lambda_3}{2} {\E}[||\nabla_{{\x}}D_v(\x)||_2^2],
  \end{aligned} \label{eq:discriminator}
\end{align}
\end{small}
where $\Theta_{E_v}$ and $\Theta_{D_v}$ are learnable parameters,
$\lambda_{1}$, $\lambda_{2}$, and $\lambda_{3}$ are the hyper-parameters,
and $F(\cdot)$ denotes the VGG feature extraction model.

Through the learned image encoder, we can map a real image into the $\mathcal{W}$ space.
The obtained code is guaranteed to align with the semantic domain of the StyleGAN generator and can be further utilized to mine cross-modal similarity between the image-text instance pairs.

\begin{figure*}[ht]
\centering
\includegraphics[width=0.90\linewidth]{images/comparison_manipulation_compressed.pdf}
\caption{Qualitative Comparison of Image Manipulation using Natural Language Descriptions.}
\label{fig:comp_man}
\end{figure*}

\subsection{Visual-Linguistic Similarity Learning}
\label{subsec:vls-model}

Once the inversion module is trained, given a real image, we can map it into the $\mathcal{W}$ space of StyleGAN.
The next problem is how to train a text encoder that learns the associations and alignments between image and text.
Instead of training a text encoder in the same way as the image encoder or the aforementioned DAMSM, we propose a visual-linguistic similarity module to project the image and text into a common embedding space, \emph{i.e.}, the $\W$ space, as shown in Figure~\ref{fig:mapping}.
Given a real image and its descriptions, we encode them into the $\W$ space by using the previously trained image encoder and a text encoder.
The obtained latent code is the concatenation of $L$ different $C$-dimensional $\w$ vectors, one for each input layer of StyleGAN.
The multi-modal alignment can be trained with
\begin{align}
  \min_{\Theta_{E_l}}\Loss_{E_l} = ||\sum_{i=1}^{L} p_i (\w^v_i - \w^l_i)||_2^2,
  \label{eq:reg}
\end{align}
where $\Theta_{E_l}$ represents the parameters of the text encoder $E_l(\cdot)$ and subscript $l$ means linguistic;
$\w^v, \w^l \in \W^{L \times C}$ are the obtained image embedding and text embedding;
$\w^v = f(E_v(\x))$ is the projected code of the image embedding $\z$ in the input latent space $\Z$ using a non-linear mapping network $f:\Z \to \W$; $\w^l$ shares a similar definition;
$\w^v$ and $\w^l$ are with the same shape $L \times C$, meaning to have $L$ \textit{layers} and each with a $C$-\textit{dimensional} latent code;
and $p_i$ is the weight of $i$-th layer in the latent code.

Compared with DAMSM, our proposed module is lightweight and easy to train. More importantly, this module achieves instance-level alignment~\cite{Wang2020CVSE}, \emph{i.e.}, learning correspondences between visual and linguistic attributes, by leveraging the disentanglability of StyleGAN.
The text encoder is trained with the proposed visual-linguistic similarity loss together with the pairwise ranking loss~\cite{kiros2014unifying,dong2017semantic}, which is omitted from Equation~\ref{eq:reg}.

\subsection{Instance-Level Optimization}
\label{subsec:instance-level-optimization}
One of the main challenges of face manipulation is the identity preservation. Due to the limited representation capability, learning a perfect reverse mapping with an encoder alone is not easy.
To preserve identity, some recent methods~\cite{ververas2020slidergan,richardson2020encoding} incorporate a dedicated face recognition loss~\cite{deng2019arcface} to measure the cosine similarity between the output image and its source.
Different from their methods, for text-guided image manipulation, we implement an instance-level optimization module to precisely manipulate the desired attributes consistent with the descriptions while faithfully reconstructing the unconcerned ones.
We use the inverted latent code $\z$ as the initialization, and the image encoder is included as a regularization to preserve the latent code within the semantic domain of the generator.
To summarize, the objective function for optimization is
\begin{align}
  \begin{aligned}
    \z^{*} = \arg\min_{\z}\ ||\x - G(\z)||_2^2\ &+ \lambda_{1}^{\prime}||F(\x) - F(G(\z))||_2^2\ \\
    &+ \lambda_{2}^{\prime}||\z - E_v(G(\z))||_2^2,
  \end{aligned} \label{eq:optimization}
\end{align}
where $\x$ is the original image to manipulate, $\lambda_{1}^{\prime}$ and $\lambda_{2}^{\prime}$ are the loss weights corresponding to the perceptual loss and the encoder regularization term, respectively.

\begin{figure}[t]
\centering
\includegraphics[width=1.0\linewidth]{images/control_mechanism.pdf}
\caption{Control Mechanism of Our TediGAN Framework.
Different layer in the StyleGAN generator represents different attributes.
Changing the value of a certain layer would change the corresponding attributes of the image.
Since the texts and images are mapped into the common latent space, we can synthesize images with certain attributes by selecting attribute-specific layers.
The control mechanism mixes layers of the style code $\w^s$ by partially replacing corresponding layers of the content $\w^c$.
When $\w^s$ is a randomly sampled latent code, it is the text-to-image generation and when $\w^s$ is the image embedding, it performs text-guided image manipulation.
}
\label{fig:control_mechanism}
\end{figure}

\subsection{Control Mechanism}
\label{subsec:control_mechanism}
\paragraph{Attribute-Specific Selection.}
The two different tasks, \emph{i.e.}, text-to-image generation and text-guide image manipulation, are unified into one framework by our proposed control mechanism.
Our mechanism is based on the style mixing of StyleGAN.
The layer-wise representation of StyleGAN learns disentanglement of semantic fragments (attributes or objects).
In general, different layer $\w_i$ represents different attributes and is fed into the $i$-th layer of the generator.
Changing the value of a certain layer would change the corresponding attributes of the image.
As shown in Figure~\ref{fig:mapping}, given two codes with the same size $\w^c, \w^s \in \W^{\;L \times C}$ denoting content code and style code, this control mechanism selects attribute-specific layers and mixes those layers of $\w^s$ by partially replacing corresponding layers of $\w^c$.
For text-to-image generation, the produced images should be consistent with the textual description, thus $\w^c$ should be the linguistic code, and randomly sampled latent code with the same size acts as $\w^s$ to provide diversity (results are shown in Figure~\ref{fig:diverse_image}).
For text-guided image manipulation, $\w^c$ is the visual embedding while $\w^s$ is the linguistic embedding, the layers for mixing should be relevant to the text, for the purpose of modifying the relevant attributes only and keeping the unrelated ones unchanged.

\vspace{-5pt}
\paragraph{Supported Modality.}
The style code $\w^s$ and content code $\w^c$ could be sketch, label, image, and noise, as shown in Figure~\ref{fig:control_mechanism}, which makes our TediGAN feasible for multi-modal image synthesis.
The control mechanism provides high accessibility, diversity, controllability, and accurateness for image generation and manipulation.
Due to the control mechanism, as shown in Figure~\ref{fig:teaser}, our method inherently supports continuous operations and multi-modal synthesis for sketches and semantic labels with descriptions.
To produce the diverse results, all we need to do is to keep the layers related to the text unchanged and replace the others with the randomly sampled latent code.
If we want to generate images from other modality with text guidance, take the sketch as an example, we can train an additional sketch image encoder $E_{vs}$ in the same way as training the real image encoder and leave the other parts unchanged.

\vspace{-5pt}
\paragraph{Layerwise Analysis.}
The pre-trained StyleGAN we used in most experiments is to generate images of 256 $\times$ 256 (\emph{i.e.}, size 256), whose has 14 layers of the intermediate vector.
For a synthesis network trained to generate images of 512 $\times$ 512, the intermediate vector would be of shape (16, 512) (and (18, 512) for 1024 $\times$ 1024), where the number of the layers $L$ is determined by $2\log_2\mathrm{R}-2$ and $\mathrm{R}$ is the image size.
In general, layers in the generator at lower resolutions (\emph{e.g.}, 4 $\times$ 4 and 8 $\times$ 8) control high-level styles such as eyeglasses and head pose,
layers in the middle (\emph{e.g.}, as 16 $\times$ 16 and 32 $\times$ 32) control hairstyle and facial expression, while the final layers (\emph{e.g.}, 64 $\times$ 64 to 1024 $\times$ 1024) control color schemes and fine-grained details.
Based on empirical observations, we list the attributes represented by different layers of a 14-layer StyleGAN in Table~\ref{tab:layerwise_analysis}.
The layers from 11-14 represent micro features or fine structures, such as stubble, freckles, or skin pores, which can be regarded as the stochastic variation.
High-resolution images contain lots of facial details and cannot be obtained by simply upsampling from the lower-resolutions, making the stochastic variations especially important as they improve the visual perception without affecting the main structures and attributes of the synthesized image.

\begin{table}[th]
\centering
\caption{The Empirical Layerwise Analysis of a 14-layer StyleGAN Generator. The 13-th and 14-th layers are omitted since there is basically no visible difference.}
\scalebox{0.85}{
\begin{tabular}{cc|cc}
\toprule
$n$-th & attribute & $n$-th & attribute \\
\midrule
1 &eye glasses & 7 & hair color  \\
2 &head pose & 8 & face color  \\
3 &face shape & 9 & age \\
4 &hair length, nose, lip & 10 & gender \\
5 &cheekbones  & 11 & micro features \\
6 &chin & 12 & micro features \\
\bottomrule
\end{tabular}
}
\label{tab:layerwise_analysis}
\end{table}

\section{Experiments}
\label{sec:experiments}

\subsection{Experiments Setup}
\label{subsec:setup}

\begin{figure}[th]
\includegraphics[width=0.9\linewidth]{images/diverse_image_compressed.pdf}
\centering
\caption{Diverse Text-to-image Generation.}
\label{fig:diverse_image}
\end{figure}

\begin{table}[th]
\centering
\caption{{Quantitative Comparison of Text-to-Image Generation.} We use FID, LPIPS, accuracy (Acc.), and realism (Real.) to compare the state of the art and our method on the proposed Multi-modal CelebA-HQ dataset. $\downarrow$ means the lower the better while $\uparrow$ means the opposite.}
\scalebox{0.8}{
\begin{tabular}{c|ccccc}
\toprule
Method &FID $\downarrow$ &LPIPS $\downarrow$ & Acc. (\%) $\uparrow$ &Real. (\%) $\uparrow$\\
\midrule
AttnGAN~\cite{xu2018attngan} &125.98 &0.512 &14.2 &20.3 \\
ControlGAN~\cite{li2019control} &116.32 &0.522 &18.2 &22.5 \\
DFGAN~\cite{tao2020dfgan} &137.60 &0.581 &22.8 &25.5 \\
DM-GAN~\cite{zhu2019dmgan} &131.05 &0.544 &19.5 &12.8 \\
TediGAN &\textbf{106.37} &\textbf{0.456} &\textbf{25.3} &\textbf{31.7} \\
\bottomrule
\end{tabular}
}
\label{tab:quan_gen}
\end{table}

\paragraph{Datasets and Baseline Models.}
To achieve text-guided image generation and manipulation, the first step is to build a dataset that contains photo-realistic facial images and corresponding descriptions.
We introduce the Multi-Modal CelebA-HQ dataset, a large-scale face image dataset that has 30,000 high-resolution face images, each having a high-quality segmentation mask, sketch, and descriptive text.
We evaluate our proposed method on text and image partitions, comparing with state-of-the-art approaches AttnGAN~\cite{xu2018attngan}, ControlGAN~\cite{li2019control}, DM-GAN~\cite{zhu2019dmgan}, and DFGAN~\cite{tao2020dfgan} for image generation, and comparing with ManiGAN~\cite{li2020manigan} for image manipulation using natural language descriptions.
All methods are retrained with the default settings on the proposed Multi-Modal CelebA-HQ dataset.

\vspace{-5pt}
\paragraph{Evaluation Metric.}
For evaluation, there are four important aspects: image quality, image diversity, accuracy, and realism~\cite{li2019control,li2020lightweight}.
The quality of generated or manipulated images is evaluated through Fr\'echet Inception Distance (FID)~\cite{heusel2017gans}.
The diversity is measured by the Learned Perceptual Image Patch Similarity (LPIPS)~\cite{zhang2018unreasonable}.
For image generation, the accuracy is evaluated by the similarity between the text and the corresponding generated image.
For manipulation, the accuracy is evaluated by whether the modified visual attributes of the synthetic image are aligned with the given description and text-irrelevant contents are preserved.
The accuracy and realism are evaluated through a user study, where the users are asked to judge which one is more photo-realistic, and more coherent with the given texts.
We test accuracy and realism by randomly sampling 50 images with the same conditions and collect more than 20 surveys from different people with various backgrounds.

\begin{table}[th]
\centering
\caption{{Quantitative Comparison of Text-Guided Image Manipulation.} We use FID, accuracy (Acc.), and realism (Real.) to compare with the state of the art ManiGAN~\cite{li2020manigan}.}
\scalebox{0.8}{
\begin{tabular}{c|cc|cc}
\toprule
&\multicolumn{2}{c}{CelebA} &\multicolumn{2}{c}{Non-CelebA} \\
Method
&ManiGAN~\cite{li2020manigan} &Ours
&ManiGAN~\cite{li2020manigan} &Ours \\
\midrule
FID $\downarrow$ &117.89 &\textbf{107.25} &143.39 &\textbf{135.47}  \\
Acc. (\%) $\uparrow$ &40.9 &\textbf{59.1} &12.8 &\textbf{87.2}  \\
Real. (\%) $\uparrow$ &36.2 &\textbf{63.8} &21.7 &\textbf{78.3}  \\
\bottomrule
\end{tabular}
}
\label{tab:quan_man}
\end{table}

\subsection{Comparison with State-of-the-Art Methods}
\label{subsec:comparison}

\subsubsection{Text-to-Image Generation}
\label{subsec:exp_gen}

\paragraph{Quantitative Comparison.}
In our experiments, we evaluate the FID and LPIPS on a large number of samples generated from randomly selected text descriptions.
To evaluate accuracy and realism, we generate images from 50 randomly sampled texts using different methods. In a user study, users are asked to judge which one is the most photo-realistic and most coherent with the given texts.
The results are demonstrated in Table~\ref{tab:quan_gen}.
Compared with the state-of-the-arts, our method achieves better FID, LPIPS, accuracy, and realism values, which proves that our methods can generate images with the highest quality, diversity, photorealism, and text-relevance.

\vspace{-5pt}
\paragraph{Qualitative Comparison.}
Most existing text-to-image generation methods, as shown in Figure~\ref{fig:comp_gen}, can generate photo-realistic and text-relevant results.
However, some attributes contained in the text do not appear in the generated image, and the generated image looks like featureless paint and lacks details.
This ``featureless painterly'' look~\cite{karras2019style} would be significantly obvious and irredeemable when generating higher resolution images using the multi-stage training methods~\cite{xu2018attngan,li2019control,zhu2019dmgan}.
Furthermore, most existing solutions have limited diversity of the outputs, even if the provided conditions contain different meanings. For example, ``\textit{has a beard}" might mean a goatee, short or long beard, and could have different colors.
Our method can not only generate results with diversity but also realise the expectation to change where you want by using the control mechanism.
To produce diverse results, with the layers related to the text unchanged, the other layers could be replaced by any values sampled from the prior distribution.
For example, as shown in the first row of Figure~\ref{fig:diverse_image}, the key visual attributes (\textit{women, black long hair, earrings, and smiling}) are preserved, while the other attributes, like haircuts, makeups, face shapes, and head poses, show a great degree of diversity.
The images in the second row illustrate more precise control ability. We keep the layers representing face shape and head pose the same and change the others.
Figure~\ref{fig:high-res-gene} shows high-quality and diverse results with resolution at 1024 $\times$ 1024.

\vspace{-5pt}
\subsubsection{Text-Guided Image Manipulation}
\label{subsec:exp_man}

\vspace{-5pt}
\paragraph{Quantitative Comparison.}
In our experiments, we evaluate the FID and conduct a user study on randomly selected images from both CelebA and Non-CelebA datasets with randomly chosen descriptions.
The results are shown in Table~\ref{tab:quan_man}. Compared with ManiGAN~\cite{li2020manigan}, our method achieves better FID, accuracy, and realism.
This indicates that our method can produce high-quality synthetic images, and the modifications are highly aligned with the given descriptions, while preserving other text-irrelevant contents.

\vspace{-5pt}
\paragraph{Qualitative Comparison.}
Figure~\ref{fig:comp_man} shows the visual comparisons between the recent method ManiGAN~\cite{li2020manigan} and ours.
As shown, the second row is to add earrings and change the face shape and hair style of the woman, our method completes this difficult case while ManiGAN fails to produce required attributes.
ManiGAN produces less satisfactory modified results: in some cases, the text-relevant regions are not modified and the text-irrelevant ones are changed.
Furthermore, since the StyleGAN we used is pretrained on a very large face dataset~\cite{karras2019style}, the latent space almost covers the full space of facial attributes, which makes our method robust for real images in the wild.
The images in last two columns are results of out-of-distribution (Non-CelebA), \emph{i.e.}, images from other face dataset such as~\cite{chelnokova2014rewards, courset2018caucasian, yi2019apdrawinggan}, which illustrate that our method is prepared to produce pleasing results with images in the wild.

\section{Ablation Study and Discussion}
\label{sec:ablation}

\paragraph{Instance-Level Optimization.}
The comparison of with or without instance-level optimization is shown in Figure~\ref{fig:inversion_result}.
As shown, the inversion results of the image encoder preserve all attributes of the original images, which is sufficient for text-to-image generation since there is no identity to preserve (Figure~\ref{fig:inversion_result} (c)).
Manipulating a given image according to a text, however, should not change the unrelated attributes especially one's identity, which is preserved after the instance-level optimization (Figure~\ref{fig:inversion_result} (d)).
We also compare with
a recent inversion-based image synthesis method pSp~\cite{richardson2020encoding} that incorporates a dedicated recognition loss~\cite{deng2019arcface} during training.
Despite both preserving the identity, the optional instance-level optimization provides a non-deterministic way to refine the final results accordingly.

\vspace{-5pt}
\paragraph{Visual-Linguistic Similarity.}
The text encoder is trained using our visual-linguistic similarity and a very simple pairwise ranking loss~\cite{kiros2014unifying,dong2017semantic} to align text and image embedding.
Although the learned text embedding can handle near-miss cases, as shown in Figure~\ref{fig:near-miss}, we found this plain strategy sometimes may lead to insufficient disentanglement of attributes and mismatching of image-text alignment, leaving some room for improvement.

\begin{figure}[t]
\includegraphics[width=0.95\linewidth]{images/inversion_result_compressed.pdf}
\centering
\caption{{Inversion Results.}
(a) original image;
(b) inversion result of pSp~\cite{richardson2020encoding};
(c) inversion result of our image encoder (Section~\ref{subsec:gan-inversion});
(d) inversion results after optimization (Section~\ref{subsec:instance-level-optimization}).
}
\label{fig:inversion_result}
\end{figure}

\begin{figure}[t]
\begin{center}
\includegraphics[width=0.95\linewidth]{images/near-miss_compressed.pdf}
\end{center}
\caption{Illustration of Near-miss Cases.}
\label{fig:near-miss}
\end{figure}

\vspace{-5pt}
\paragraph{Potential Issue with StyleGAN.}
In our experiments, we found that some unrelated attributes are unwantedly changed when we manipulate a given image according to a text description.
We thought it might be the problem of visual-linguistic similarity learning in the first place.
However, when performing layer-wise style mixing on the inverted codes of two real images, the interference still occurs.
This means some facial attributes remain entangled in the $\mathcal{W}$ space, where different attributes should be orthogonal (meaning without affecting other attributes).
Another inherent defect of StyleGAN is that some attributes, such as hats, necklaces and earrings, are not well represented in its latent space. This makes our method perform less satisfactorily sometimes when adding or removing jewelry or accessories through natural language descriptions.

\section{Conclusion}
\label{sec:conclusion}
We have proposed a novel method for image synthesis using textual descriptions, which unifies two different tasks (text-guided image generation and manipulation) into the same framework and achieves high accessibility, diversity, controllability, and accurateness for facial image generation and manipulation.
Through the proposed multi-modal GAN inversion and large-scale multi-modal dataset, our method can effectively synthesize images with unprecedented quality.
Extensive experimental results demonstrate the superiority of our method, in terms of the effectiveness of image synthesis, the capability of generating high-quality results, and the extendability for multi-modal inputs.

{\small
\bibliographystyle{ieee}


\section{Introduction}

The enigmatic quantum spin liquid~(QSL) has drawn a lot of attention
ever since the seminal work of Anderson in 1973 \cite{Anderson1973}.
In 2006, Kitaev proposed an exactly solvable spin-$1/2$ model on the honeycomb lattice
and demonstrated that its ground state is an exotic QSL
with emergent Majorana fermion excitations \cite{Kitaev2006}.
The past decade has witnessed a surge of interest in realization of the Kitaev honeycomb model
on real materials with $4d$ or $5d$ magnetic ions, which includes iridates and $\alpha$-RuCl$_3$
(see Refs.~\cite{Jackeli2009,RauLeeKee2016,TakagiTJ2019} and referees therein).
However, because of the inevitable non-Kitaev interactions,
e.g., the Heisenberg interaction and a symmetric off-diagonal exchange $\Gamma$-interaction \cite{RanLeeKeePRL2014},
these materials are shown to display magnetic orders at lowest temperatures
\cite{LiuBYetal2011,ChaloupkaJH2013,PlumbCSetal2014,JohnsonWHetal2015}.
Nevertheless, it is believed that the effective $K$-$\Gamma$ model is the dominant ingredient
to describe $\alpha$-RuCl$_3$ \cite{WangDYLi2017}.

From a theoretical point of view,
although the quantum phase diagram of the $K$-$\Gamma$ model on a honeycomb lattice is elusive,
several magnetically ordered phases and distinct QSLs are demonstrated to exist
\cite{RanLeeKeePRL2014,CatunYWetal2018,GohlkeWYetal2018,WangBL2019,LuoZhaoKeeWang2019,YamadaSS2020},
indicating the strong quantum fluctuation enhanced by competing interactions.
Given the notorious difficulty in two dimension,
it is beneficial and constructive to reduce the dimensionality where many full-fledged analytical and numerical methods
capable of addressing problems in one-dimensional~(1D) quantum spin chains are available.
Recently, the phase diagram of the isotropic $K$-$\Gamma$ chain has been studied
by the density-matrix renormalization group~(DMRG) method
and the non-Abelian bosonization technique \cite{YangKG2020,YangSN2021}.
It is shown that about 2/3 of the phase diagram is occupied by a gapless Luttinger liquid (LL).
The ferromagnetic~(FM) Kitaev limit is merely a transition point,
while a critical segment near the antiferromagnetic~(AFM) Kitaev limit is identified.
Two symmetry-breaking phases termed the FM$_{U_6}$ phase and the $M_2$ phase
(see Fig.~\ref{FIG-GSPD} for the nomenclature of the magnetically ordered phases) are also reported.
Later on, it is found that FM and AFM Heisenberg interactions could open up a wide region
of the FM$_{U_6}$ phase and the LL, respectively \cite{YangJKG2020}.
However, how to enlarge the territory of the puzzling $M_2$ phase is still unclear.
Aligning with this effort, a two-leg $K$-$\Gamma$ ladder under a [111] magnetic field is also studied,
revealing a rich phase diagram with several emergent phases \cite{SorenseCGK2021}.

Aside from the exotic phases and quantum criticality,
quantum spin chains also provide an excellent platform for theoretical studies of various quantum phase transitions~(QPTs) \cite{Sachdev2011},
of particular interest is the topological QPT that is beyond Landau's paradigm.
The topological QPT occurs between two different phases
without any explicit symmetry breaking \cite{KT1973,TsuiSG1982,HaldanePRL1983,WenRMP2017}.
The Haldane phase is such an example of symmetry-protected topological~(SPT) phase \cite{PollmannSPT2012},
which possesses a nonlocal string order parameter~(SOP) due to a hidden $\mathbb{Z}_2\times\mathbb{Z}_2$ symmetry breaking \cite{denNijsRom1989,KennedyTasaki1992},
Dating back to 1992, Hida originally pointed out that the bond-alternating spin-$1/2$ Heisenberg chain could host the Haldane phase
due to the imbalance of the neighboring coupling intensities,
leading to the formation of either total spin 0 or 1 out of the two spin-$1/2$ degrees of freedom \cite{Hida1992}.
Therefore, bond alternation is a practical route to legalize the validity of the SPT phase in spin-$1/2$ chains \cite{ChenGuWen2011}.
We also note that the anisotropic Kitaev spin chain hosts two disordered phases which undergo a direct transition at the isotropic point \cite{BrzezickiDO2007,YouTian2008}.
These observations motivate us to investigate the ground-state properties of the $K$-$\Gamma$ chain
by altering the bond strength of adjacent sites.

In this paper, we study the phases and QPTs of a bond-alternating $S$ = $1/2$ $K$-$\Gamma$ chain.
When $\Gamma = 0$, it is the Kitaev spin chain,
otherwise known as the exactly solvable 1D quantum compass model~(QCM) \cite{BrzezickiDO2007,YouTian2008}.
Beyond that it is nonintegrable except for some special points and lines when $|K|$ = $|\Gamma|$.
Therefore, we resort to the DMRG method \cite{White1992,Peschel1999,Schollwock2005}
to map out the quantum phase diagram.
The phase boundaries are determined by various quantities including the energy gap and entanglement entropy.
The central charge is calculated to distinguish the universality class of a continuous QPT.

The structure of the paper is as follows.
In Sec.~\ref{SEC:Model} we introduce the theoretical model under investigation,
analyse the symmetry properties, and present the phase diagram of interest.
Following this, we study two topological QPTs in Secs.~\ref{SEC:EHOH} and \ref{SEC:QCM}.
Section~\ref{SEC:SSB} presents the magnetic order parameters of symmetry-breaking phases.
In Sec.~\ref{SEC:FMKtv} we study the transitions between the left and right panels of the phase diagram.
We conclude with a summary in Sec.~\ref{SEC:Conclusion}.
Finally, a brief review of the diagonalization of QCM
and some other useful contents are presented in the Supplemental Material \cite{SuppMat}.

\section{Model and Method}\label{SEC:Model}
We consider a bond-alternating spin-$1/2$ $K$-$\Gamma$ chain with
\begin{align}\label{J1J2KG-Ham}
\mathcal{H} = \sum_{l=1}^{L/2} g_x\mathcal{H}_{2l-1,2l}^{(x)}(\theta) + g_y\mathcal{H}_{2l,2l+1}^{(y)}(\theta)
\end{align}
where $L$ is the chain length, $g_x$~($g_y$) is the odd~(even) bond strength, and
\begin{align}\label{KtvGam-Ham}
\mathcal{H}_{i,j}^{(\gamma)}(\theta) = K S_i^{\gamma}S_j^{\gamma} + \Gamma (S_i^{\alpha}S_j^{\beta}+S_i^{\beta}S_j^{\alpha}).
\end{align}
Here, $K$ and $\Gamma$ are the Kitaev interaction and the off-diagonal exchange interaction, respectively.
$\gamma$ could be either $x$ or $y$ and it specifies the spin direction associated with the referred bond, see Fig.~\ref{FIG-Bond}(a).
For each $\gamma$-bond, $\alpha$ and $\beta$ are the two remaining mutually exclusive spin directions.
In what follows we parametrize $K=\sin\theta$ and $\Gamma=\cos\theta$ with $\theta \in (-\pi, \pi]$.

\begin{figure}[!ht]
\centering
\includegraphics[width=0.99\columnwidth, clip]{Fig01.pdf}\\
\caption{(a) Sketch of the bond structure in the original form.
    Here, $x$~(red) and $y$~(green) stand for the $\gamma$-index
    and the width is proportional to the bond strength.
    (b) Pictorial bond structure of the Hamiltonian in the rotated basis.
    The overhanging bond at each site is determined by the remaining one along the chain.
    The shaded region represents the six-site unit cell.}
    \label{FIG-Bond}
\end{figure}

Before carrying out numerical calculation,
let us analyze the symmetries in the parameter space
which will reduce the computational cost.
Akin to the transverse field Ising model \cite{Pfeuty1970},
the model \eqref{J1J2KG-Ham} owns a duality relation which could be seen by applying the spin rotation transformation
$(S_i^x, S_i^y, S_i^z) \to (-S_i^y, -S_i^x, -S_i^z)$.
This implies that each eigenvalue $E$ of $H$ satisfies the relation
\begin{equation}\label{SelfDual}
E(g) = gE(1/g)
\end{equation}
where $g \equiv g_y/g_x$ is the relative bond strength.
On the other hand, by virtue of a global spin rotation around the $z$-axis by $\pi$,
$(S_i^x, S_i^y, S_i^z) \to (S_i^y, -S_i^x, S_i^z)$,
the Kitaev interaction remains uninfluenced whereas the sign of $\Gamma$-interaction is altered \cite{YangKG2020}.
We thus instantly find that
\begin{equation}\label{GammaSym}
E(K, \Gamma) = E(K, -\Gamma),
\end{equation}
or equivalently, $\theta \mapsto \pi-\theta$.
These relations in Eq.~\eqref{SelfDual} and Eq.~\eqref{GammaSym}
allow us to focus on the phase diagram primarily in the reduced parameter range
$\theta \in [-\pi/2, \pi/2]$ and $g = g_y/g_x \in [0, 1]$
and then map out the whole phase diagram shown in Fig.~\ref{FIG-GSPD}.

Using a site-ordering cross decimation rotation with a periodicity of six sites,
all the cross terms of $S_i^{\alpha}S_j^{\beta}$ with $\alpha \neq \beta$ in Eq.~\eqref{KtvGam-Ham} will vanish \cite{YangKG2020}.
This $U_6$ transformation is given by
\begin{eqnarray}\label{EQ:SOCD}
\text{sublattice $1$}: & (x,y,z) & \rightarrow (\tilde{x},\tilde{y},\tilde{z}),      \nonumber\\
\text{sublattice $2$}: & (x,y,z) & \rightarrow (-\tilde{x},-\tilde{z},-\tilde{y}),   \nonumber\\
\text{sublattice $3$}: & (x,y,z) & \rightarrow (\tilde{y},\tilde{z},\tilde{x}),      \nonumber\\
\text{sublattice $4$}: & (x,y,z) & \rightarrow (-\tilde{y},-\tilde{x},-\tilde{z}),   \nonumber\\
\text{sublattice $5$}: & (x,y,z) & \rightarrow (\tilde{z},\tilde{x},\tilde{y}),      \nonumber\\
\text{sublattice $6$}: & (x,y,z) & \rightarrow (-\tilde{z},-\tilde{y},-\tilde{x}),
\end{eqnarray}
where $\gamma \big(= x (\tilde{x}), y (\tilde{y}), z (\tilde{z})\big)$ denotes the spin component of $S^{\gamma}$ ($\tilde{S}^{\gamma}$).
Under this transformation the original Hamitonian acquires the following form \cite{YangKG2020}
\begin{align}\label{U1XYZ-Ham}
\tilde{\mathcal{H}}_{i,j}^{(\gamma)}(\theta) = -K \tilde{S}_i^{\gamma} \tilde{S}_j^{\gamma}
- \Gamma (\tilde{S}_i^{\alpha} \tilde{S}_j^{\alpha}+ \tilde{S}_i^{\beta} \tilde{S}_j^{\beta})
\end{align}
in which the bonds $\gamma$ = $\tilde{x}$ (red), $\tilde{z}$ (blue), and $\tilde{y}$ (green) circularly, as depicted in Fig.~\ref{FIG-Bond}(b).
$\tilde{\mathbf{S}} = (\tilde{S}_i^{x}, \tilde{S}_i^{y}, \tilde{S}_i^{z})$ is the spin operator in the rotated basis.
Such a $U_6$ transformation does not alter the energy spectra (i.e., energy and its degeneracy)
but simplifies the spin-spin correlation functions.
Therefore, we will preferentially focus on the rotated Hamiltonian in Eq.~\eqref{U1XYZ-Ham} unless stated explicitly otherwise.
The exceptions are Secs.~\ref{SEC:QCM} and \ref{SEC:FMKtv}
where it is convenient to calculate the correlation functions in the original basis.
In addition, combining Eq.~\eqref{GammaSym},
it is apparently that Eq.~\eqref{U1XYZ-Ham} has a $SU(2)$ symmetric structure
when $|K| = |\Gamma|$.
Specifically, In the range $\theta \in [-\pi/2, \pi/2]$,
the point $\theta = -\pi/4$ and $\pi/4$ corresponds to an
AFM and FM Heisenberg chain, respectively.

\begin{figure}[!ht]
\centering
\includegraphics[width=0.90\columnwidth, clip]{Fig02.pdf}\\
\caption{Quantum phase diagram of the bond-alternating spin-$1/2$ $K$-$\Gamma$ chain with $K=\sin\theta$ and $\Gamma=\cos\theta$.
    The green thick line marked by $\textcircled{1}$ is the isotropic $K$-$\Gamma$ chain,
    and the black solid circles at $\theta = \pi/4$ and $-\pi/4$ represent the hidden $SU(2)$ FM and AFM Heisenberg chains, respectively.
    The asterisk ($\ast$) in the vertical line represents a multicritical point.
    There are seven distinct phases in the right panel which are the main focus of the paper.
    The EH--OH transition and $A_x$--$A_y$ transition are continuous with a central charge of $c = 1$ and $1/2$, respectively.
    The nature of the magnetically ordered states, the FM$_{U_6}$ phase, the $M_1$ phase, and the $M_2$ phase,
    are clarified in a six-site $U_6$ rotated basis.
    See the main text for details.}
    \label{FIG-GSPD}
\end{figure}

The numerical calculations are performed by the DMRG method
\cite{White1992,Peschel1999,Schollwock2005},
which is a powerful technique for 1D many-body problems.
Periodic boundary condition~(PBC) is preferred to weaken the finite-size effect
and open boundary condition~(OBC) is also adopted occasionally for comparison.
We keep up to 2000 states so as to ensure a typical truncated error of $\sim\!10^{-7}$ or less.
The chain length $L$ is strictly considered to be the multiple of 6,
consistent with the structure of unit cell and the $U_6$ transformation.

The resultant phase diagram is shown in Fig.~\ref{FIG-GSPD},
which has a salient feature of mirror (left-right) symmetry.
Focusing on the right half circle, it harbours \textit{seven} distinguishing phases.
Four of them, i.e., the even-Haldane~(EH) and odd-Haldane~(OH) phases
and the $A_x$ and $A_y$ phases \cite{Kitaev2006},
are disordered and could be characterized by nonlocal SOPs of different kinds.
The rest are three magnetically ordered phases named FM$_{U_6}$ phase and $M_1$ and $M_2$ phases.
The FM$_{U_6}$ phase is collinear in the rotated basis and exhibits six-fold degeneracy.
The $M_1$ and $M_2$ phases show more complex spin patterns where all their three spin components are finite.
For the $M_1$ phases, it is stabilized at the region where $g \lesssim 1/\sqrt3$ or $g \gtrsim \sqrt3$
and \textit{one} of the spin components is dominantly the biggest.
For the $M_2$ phases, it locates around the very isotropic line of $g \simeq 1$ where $K/\Gamma > 1$
and \textit{two} of its spin components are almost the same and is larger than the third.

\section{EH-OH topological QPT}\label{SEC:EHOH}

Straightforwardly, when $\theta = -\pi/4$ Eq.~\eqref{J1J2KG-Ham} and Eq.~\eqref{U1XYZ-Ham}
turn out to be a bond-alternating AFM Heisenberg chain \cite{Hida1992,Barnes1999,Johnston2000,WangLiCho2013,Haghshenas2014}.
It is well-established that there is a topological EH-OH transition at $g = 1$ with a central charge $c = 1$ \cite{WangLiCho2013}.
For either $g < 1$ or $g > 1$, the ground state could be characterized by a SOP
which is nonzero inside the phase but vanishes otherwise (see Sec. II in the Supplemental Material \cite{SuppMat}).
Specifically, the two phases could be distinguished by the even- and odd-SOPs which are defined as \cite{Hida1992}
\begin{equation}\label{EQ:SOPev}
\mathcal{O}_{e}^{\alpha} = \lim_{|j-i|\to\infty}\mathcal{O}^{\alpha}(2i,2j+1)
\end{equation}
and
\begin{equation}\label{EQ:SOPod}
\mathcal{O}_{o}^{\alpha} = \lim_{|j-i|\to\infty}\mathcal{O}^{\alpha}(2i-1,2j)
\end{equation}
where
\begin{equation}\label{EQ:SOPij}
\mathcal{O}^{\alpha}(p,q) = -4\left\langle \tilde{S}_p^{\alpha} \Big(\prod\limits_{p<r<q} e^{i\pi \tilde{S}_r^{\alpha}}\Big) \tilde{S}_q^{\alpha} \right\rangle.
\end{equation}
Here, $\alpha = x, y, z$.
The even-SOP $\mathcal{O}_{e}^{\alpha}$ is calculated from an even site ($2i$) to an odd site ($2j+1$)
while the odd-SOP $\mathcal{O}_{o}^{\alpha}$ is measured from an odd site ($2i-1$) to an even site ($2j$).
At the critical point $g = 1$, both SOPs are identical and decay as $\mathcal{O}_{e/o} \sim L^{-1/4}$ \cite{BortzSS2007}.
From a topological perspective, the EH phase is trivial while the OH is a SPT phase
which is isomorphic to the ground state of the spin-1 Heisenberg chain \cite{HaldanePRL1983}.
For the OH phase, its ground state is unique under PBC
but has a fourfold degeneracy under OBC because of two edge spin-$1/2$s.
In addition, the degeneracy of the lowest-lying entanglement spectrum
is twofold~(fourfold) under OBC~(PBC) \cite{Pollmann2010}.

\begin{figure}[!ht]
\centering
\includegraphics[width=0.95\columnwidth, clip]{Fig03.pdf}\\
\caption{(a) Entanglement entropy $\mathcal{S}_1$ (open symbols) and $\mathcal{S}_2$ (filled symbols) of different cuts
    for the $g_x$-$g_y$ $\Gamma$-chain with $\theta= 0.00\pi$.
    The chain length $L$ is 48~(red circle), 96~(green triangle), and 144~(blue square).
    (b) Bulk entanglement gap $\Delta\mathcal{S}$ in the same region as (a).
    $W = -1$ and 1 are the sign of $\Delta\mathcal{S}$ when $g < 1$ and $g > 1$, respectively.}
    \label{FIG-BulkDltVNE}
\end{figure}

As shown in Fig.~\ref{FIG-GSPD},
the EH and OH phases extend to a large region of the parameter space.
To demonstrate it, we focus on the line of $\theta = 0$, which is the $\Gamma$-chain limit.
We begin by studying a so-called \textit{bulk entanglement gap} $\Delta\mathcal{S}$ \cite{Tan2020},
which comes from the even-odd oscillation of the entanglement entropy
$\mathcal{S}_L(l) = -\mathrm{Tr}(\rho_l\ln\rho_l)$
where $\rho_l$ is the reduced density matrix of the subsystem
with a contiguous spatial segment $l$ \cite{VidalLRK2003}.
Depending on whether $l \gg 1$ is odd or even,
$\mathcal{S}_L(l)$ saturates to a constant value of $\mathcal{S}_1$ and $\mathcal{S}_2$, respectively.
The bulk entanglement gap is thus defined as $\Delta\mathcal{S} = \mathcal{S}_2 - \mathcal{S}_1$ \cite{Tan2020}.
Figure~\ref{FIG-BulkDltVNE}(a) shows $\mathcal{S}_1$~(open symbols) and $\mathcal{S}_2$~(fill symbols)
for the chain length $L$ = 48, 96, and 144.
For these lengths chosen, $\mathcal{S}_1$ corresponds to cut one strong valence bond consistently,
while $\mathcal{S}_2$ stands for cutting a strong valence bond
zero time or twice when $g < 1$ or $g > 1$, respectively.
This implies that the  bulk entanglement gap $\Delta\mathcal{S}/\ln2$ tends to be $-1$ or $1$ at the limit case
where $g\to0$ or $g\to\infty$, respectively~(see Fig.~\ref{FIG-BulkDltVNE}(b)).
Near the critical point, $\mathcal{S}_1 \simeq \mathcal{S}_2$ and thus $\Delta\mathcal{S} \simeq 0$.
Moreover, defining $\delta = \frac{g-1}{g+1}$,
the quantity scales as $\Delta\mathcal{S} \sim \ln2-(1-|\delta|)^2$ when away from criticality,
whereas $\Delta\mathcal{S} \sim -\delta\ln|\delta|$ when close to criticality \cite{Tan2020}.
Therefore, the value of $\Delta\mathcal{S}$ is bounded to $\pm\ln2$
and its sign $W = \textrm{sgn}\big(\Delta\mathcal{S}\big)$
could be used to characterize the corresponding topological sector.
The sign change~(e.g., from $W = -1$ to 1) is a signal for the topological QPT.

\begin{figure}[!ht]
\centering
\includegraphics[width=0.95\columnwidth, clip]{Fig04.pdf}\\
\caption{(a) SOPs of the even type $\mathcal{O}_{e}^z$ (open symbols) and odd type $\mathcal{O}_{o}^z$ (filled symbols)
    for the $g_x$-$g_y$ $\Gamma$-chain with $\theta= 0.00\pi$.
    The inset shows the asymptotic decay of SOP $\mathcal{O}$ at $g_x = g_y$.
    (b) Energy gap $\Delta_L$ in the same region as (a).
    Inset shows the linear extrapolation of the energy gap at $g_x = g_y$.}\label{FIG-GamSOPGap}
\end{figure}

To further reveal the nature of phases at different topological sectors,
we measure the even-SOP $\mathcal{O}_{e}^{z}$ (see Eq.~\eqref{EQ:SOPev})
and odd-SOP $\mathcal{O}_{o}^{z}$ (see Eq.~\eqref{EQ:SOPod}).
It is clearly shown in Fig.~\ref{FIG-GamSOPGap}(a) that
$\mathcal{O}_{e}^{z}$~($\mathcal{O}_{o}^{z}$) is finite when $g < 1$~($g > 1$) and is vanishingly small otherwise.
The finite-size effects of $\mathcal{O}_{e/o}^{z}$ are very weak,
except for a narrow window that is close to the critical region.
As shown in the inset, both types of SOPs $\mathcal{O}_{e/o}^{z}$ decay algebraically as $L^{-\alpha}$
where the critical exponent $\alpha \approx 0.26$, which is fairly close to the value of $1/4$ at $\theta = -\pi/4$ \cite{BortzSS2007}.
For an infinite size system, SOPs $\mathcal{O}_{e/o}^{z}$ scale as $\delta^{1/6}$ \cite{Hida1992}.
As a result, the critical exponent is given as $\beta = 1/12$ because $\mathcal{O} \propto \delta^{2\beta}$.
Hence, this topological QPT belongs to the Gaussian universality class.
We also calculate the excitation gap $\Delta_T$,
which is defined as the energy difference between the first excited state and the ground state.
Figure~\ref{FIG-GamSOPGap}(b) shows that $\Delta_T$ is very robust when $g \neq 1$.
Near $g = 1$, it has a pronounced drop with size increased.
As shown in the inset, $\Delta_T$ is zero when $L\to\infty$,
showing that the ground state of the isotropic $\Gamma$-chain is critical.
To extract the central charge $c$, we calculate the von Neumann entanglement entropy $\mathcal{S}_L(L/2)$
for a series of chain length $L$ and the central charge is fitted by
$\mathcal{S}_L = \frac{c}{3}\ln(L/\pi) + c'$.
Our best fitting suggests that $c \simeq 0.997(5)$ (not shown),
which is very close to 1 of the LL.

\begin{figure}[!ht]
\centering
\includegraphics[width=0.95\columnwidth, clip]{Fig05.pdf}\\
\caption{(a) Behaviors of the energy density $w = E_g(L)/L$ for the isotropic $\Gamma$-chain under OBC~(black) and PBC~(red).
    The inset shows the six-site periodicity of $w$ under PBC.
    (b) and (c) show the estimate of $e_g$ for the isotropic $\Gamma$-chain under OBC and PBC, respectively.
    The roman numerals \upperRomannumeral{1} (green and red) and \upperRomannumeral{2} (blue)
    mark two different methods illustrated in the text.}\label{FIG-GamEg}
\end{figure}

Figure~\ref{FIG-GamEg}(a) displays the energy density $E_g/L$ of the isotropic $\Gamma$-chain
under OBC (black triangle) and PBC (red circle).
For the OBC case, $E_g/L$ decreases smoothly and saturates around $-0.30$ as $L$ increases.
In contrast, it is not monotonically increasing but exhibits an oscillation with six-site periodicity for the PBC (see inset).
As a comparison, we note that such an abnormal energy density behavior is absent in the isotropic Kitaev spin chain
(see Fig. 1 in the Supplemental Material \cite{SuppMat}).
This phenomenon in the $\Gamma$-chain is striking and may be related to the unusual energy behavior
of the $\Gamma$ model on the honeycomb lattice \cite{LuoZhaoKeeWang2019}.
In Ref.~[\onlinecite{LuoZhaoKeeWang2019}] the total energy $E_g$ is calculated on a series of honeycomb clusters
where OBC~(PBC) is utilized on the $L_x$~($L_y$) direction of cylinders.
For any cylinder with fixing $L_y$, the energy density $E_g/N$ ($N=L_xL_y$) varies linearly with $1/L_x$.
However, by increasing the circumference of the cylinders with $L_x/L_y = 2$,
the energy density $E_g/N$ is no longer monotonous and exhibits a skew sawtooth behavior.

To round off the calculation, we give an estimate of the ground-state energy per-site $e_g$ of the isotropic $\Gamma$-chain.
We note that our $\Gamma$-chain contains both $x$ and $y$ bonds~(see Eq.~\eqref{J1J2KG-Ham}),
and there is no analytical solution so far.
It is fundamentally different from a $z$-bond $\Gamma$-chain
which could be solved exactly via the Jordan-Wigner transformation \cite{YouGam2020}.
At the quantum critical point, the finite-size scaling of the ground-state energy $E_g(L)$ is known to be \cite{BloteCN1986,Affleck1986}
\begin{equation}\label{EQ:EgFSS}
E_g(L) = Le_g + \varepsilon_b - \frac{\Delta_b}{L} + \mathcal{O}(L^{-2}),
\end{equation}
where $e_g$ is the average bulk energy per-site,
$\varepsilon_b$ is the size-independent surface energy which vanishes in the case of PBC,
and $\Delta_b$ is the subleading correlation term.
It is found that $\Delta_b = \pi c/6$~($\pi c/24$) for PBC~(OBC) where $c$ is the central charge \cite{BloteCN1986,Affleck1986}.
By definition we have $e_g = \lim_{L\to\infty} e_L$
where $e_L$ is the energy per-site of the chain with length $L$.
For the energy obtained under the OBC,
there are two ways to extrapolate it to the thermodynamic limit;
one is $e_L^{\textrm{I}} = E_g(L)/L$ and the other is $e_L^{\textrm{II}} = \big(E_g(L)-E_g(L\!-\!2)\big)/2$.
It is easy to check that convergence speed of the latter is faster than the former.
As shown in Fig.~\ref{FIG-GamEg}(b), the quadratic fittings of the two give that
$e_g^{\textrm{I}} = -0.29959362$ and $e_g^{\textrm{II}} = -0.29959375$,
yielding an estimate for the ground-state energy per-site in the thermodynamic limit
of $e_g = -0.2995937(1)$ with seven significant digits.
Meanwhile, we also extrapolate the energy under PBC
by using the solid points in Fig.~\ref{FIG-GamEg}(a) where $L$ is a multiple of six~(see Fig.~\ref{FIG-GamEg}(c)).
Our result suggests that $e_g \approx -0.299594$,
which is fairly consistent with the high-precision value revealed by the calculation under OBC.

\section{Extended quantum compass model}\label{SEC:QCM}

In the absence of $\Gamma$-interaction,
Eq.~\eqref{J1J2KG-Ham} and Eq.~\eqref{KtvGam-Ham} are reduced to the 1D Kitaev spin chain,
which is also known as the 1D QCM in some other context \cite{BrzezickiDO2007,YouTian2008}.
The QCM could be solved exactly by Jordan-Wigner transformation
and its dispersion relation is almost the same as that of the transverse field Ising model.
There is a topological QPT between the gapped $A_x$ and $A_y$ phases at $g = 1$ \cite{Kitaev2006}.
However, because of intermediate symmetries,
the ground state of the QCM possesses a huge number of degeneracy $2^{N/2-1}$~($2^{N/2}$) under PBC~(OBC)
where $N$ is the total number of sites \cite{BrzezickiDO2007,NussinovBrink2015}.
Equivalently, the QCM could be rewritten as a Majorana fermion chain complemented by $N$ decoupled Majorana fermions.
Since each Majorana fermion has $\sqrt2$ degrees of freedom,
the redundant Majorana fermions thus contribute a ground-state degeneracy of $\mathcal{O}(2^{N/2})$ \cite{FengZX2007}.
These degenerate ground states are vulnerable and can be totally lifted by an infinitesimal transverse field \cite{SunChen2009}.
The entanglement \cite{Subrahmanyam2013,YouWangYi2018,VimalSubra2018,WuYou2019},
energy dynamics \cite{SteinBrenig2016}, and the dissipative behavior \cite{ShibataKat2019}
of the QCM have been studied over the years.

Using the spin duality transformation, the topological $A_x$ and $A_y$ phases, respectively,
are shown to possess nonlocal string correlators \cite{FengZX2007},
\begin{equation}\label{EQ:SOPx}
\mathcal{O}_K^x(2r) = \lim_{r\to\infty}\left\langle\prod_{k=1}^{2r} \sigma_k^x\right\rangle
\end{equation}
and
\begin{equation}\label{EQ:SOPy}
\mathcal{O}_K^y(2r) = \lim_{r\to\infty}\left\langle\prod_{k=2}^{2r+1} \sigma_k^y\right\rangle,
\end{equation}
where $\sigma_k^x$ and $\sigma_k^y$ are Pauli matrices, i.e., twice the spin-1/2 operators in the \textit{original} basis.
Here, the nonlocal SOPs are defined in the original basis embedded in Eq.~\eqref{KtvGam-Ham}.
Generalization of nonlocal SOPs to a two-leg Kitaev ladder has been discussed in a recent work \cite{Catuneanu2019}.
At the critical point $g = 1$, these Kitaev-type SOPs vanish in an algebraic behavior at long-distance limit $n\gg1$ \cite{Pfeuty1970},
\begin{equation}\label{EQ:KSOPAsy}
\mathcal{O}_K^{x/y}(n) =e^{1/4}2^{1/12}A^{-3} n^{-1/4}\Big(1-\frac{1}{64}n^{-2} + \cdots\Big),
\end{equation}
where $A \simeq 1.2824$.
For infinite-size case, they obey a scaling law and $\mathcal{O}_K^{x} \sim (1-g^2)^{1/4}$ when $g\to1^{-}$.

\begin{figure}[!ht]
\centering
\includegraphics[width=0.95\columnwidth, clip]{Fig06.pdf}\\
\caption{Kitaev-type SOP $\mathcal{O}_K^{x/y}$ of the $g_x$-$g_y$ $K$-$\Gamma$ chain for $\theta=0.48\pi$
    with chain length $L = 48$~(red circle) and 72~(green triangle).
    The inset shows the behavior of $\mathcal{O}_K^{x}$ at the isotropic point $g = 1$.
    Values at $\theta = 0.50\pi$~(blue square) and $\theta = 0.48\pi$~(red cross) are shown
    for several chain length $L$ = 48, 72, 96, 144, and 216.
    The sold line is the correlation function defined in Eq.~\eqref{EQ:KSOPAsy}.}
    \label{FIG-Tht048}
\end{figure}

Hereafter we show numerically that,
the topological $A_x$ and $A_y$ phases are extended
when $\theta$ is slightly deviated from $\pi/2$ (AFM Kitaev point).
We demonstrate this by calculating the SOPs shown in Eq.~\eqref{EQ:SOPx} and Eq.~\eqref{EQ:SOPy} for $\theta = 0.48\pi$ with $2r = L/2$.
As presented in Fig.~\ref{FIG-Tht048},
the SOPs change smoothly and are very robust in each corresponding phase,
showing the validity of them in this region.
When $g = 1$, the two have the same value due to the self-dual relation and they decrease visibly as $L$ grows.
To measure how the SOPs vary at this point,
we calculate $\mathcal{O}_K^{x}$ for $L$ up to 216 sites and the results are shown in the inset.
The values at $\theta = 0.50\pi$ are also shown for comparison.
As revealed by Eq.~\eqref{EQ:KSOPAsy},
the leading term of $\mathcal{O}_K^{x}$ is $\sim1/L^{1/4}$,
so its decay ratio is not very rapid for modest chain length $L$ as shown by the solid line in the inset.
However, it is constructive to note that $\mathcal{O}_K^{x}$ at $\theta = 0.48\pi$ and $0.50\pi$ are very close
but the curve of $\theta = 0.48\pi$ is shifted downward slightly when compared with the latter.
For $\theta = 0.50\pi$ it is known that $\mathcal{O}_K^{x}$ vanishes at $g = 1$ when $L\to\infty$ \cite{FengZX2007}.
So it is reasonable for us to believe that it will also go to zero ultimately for $\theta = 0.48\pi$.
As a result, $g = 1$ is still inferred as the critical point for the $A_x$--$A_y$ transition.

\begin{figure}[!ht]
\centering
\includegraphics[width=0.95\columnwidth, clip]{Fig07.pdf}\\
\caption{Estimate of the central charge $c$ for the critical $K$-$\Gamma$ chain at ($g = 1.00$, $\theta = 0.48\pi$).}\label{FIG-Tht048CC}
\end{figure}

To confirm the criticality at $\theta = 0.48\pi$, we now turn to calculate the central charge.
The central charge is usually extracted from the coefficient of the logarithmic correlation
in the entanglement entropy \cite{VidalLRK2003}.
However, this method is not optimal for the critical Kitaev phase
because of the macroscopic ground-state degeneracy \cite{BrzezickiDO2007}.
As a result, it is challenging to get a minimally entangled state
which is essential for a reliable estimate of the central charge.
The practical way to handle this problem is by the energy scaling as shown in Eq.~\eqref{EQ:EgFSS}.
For the PBC, the central charge is given by the following formula
\begin{equation}\label{EQ:CCEgFit}
c_L \simeq \frac{6}{\pi}\big(Le_g-E_g(L)\big)L
\end{equation}
where $e_g$ is the only relevant parameter.
Following a similar procedure illustrated in Fig.~\ref{FIG-GamEg}(b),
we find $e_g\approx-0.1591092$ for $\theta = 0.48\pi$,
which is only slightly larger than that of $-1/(2\pi) = -0.1591549\cdots$ for $\theta = \pi/2$.
We have also calculated the ground-state energy $E_g(L)$ for a series of chain length $L$ ranging from 24 to 144.
The fitting central charge via Eq.~\eqref{EQ:CCEgFit} is shown in Fig.~\ref{FIG-Tht048CC}.
It can be found that the central charge is very close to $1/2$ and suffers from a tiny finite-size effect.
Therefore, we draw the conclusion that the central charge $c = 1/2$,
and the transition belongs to the same universality class as that at $\theta = \pi/2$ \cite{WangCho2015,YangPH2020},
confirming the existence of an extended region of $A_x$ and $A_y$ phases and the critical transition line between them.

\section{The symmetry breaking phases}\label{SEC:SSB}

\subsection{Degeneracy and spin patterns}
Like the AFM case shown in Sec.~\ref{SEC:EHOH},
Eq.~\eqref{J1J2KG-Ham} and Eq.~\eqref{U1XYZ-Ham} could be reduced to
the bond-alternating FM Heisenberg chain when $\theta = \pi/4$ (i.e., $K = \Gamma$).
For this model its ground-state energy $E_g = -(1+g)KL/8$ with a $(L+1)$-fold ground-state degeneracy \cite{KomaNach1997}.
Although it is gapless inherently, the system is not conformally invariant.
Specially, when $g = 1$ it is shown that there is an effective central charge $c_{\rm eff} = 3/2$ \cite{ChenXMetal2013}.
Around the isotropic $SU(2)$ FM point by tuning $\theta$ along the line of $g = 1$, there is an $O_h\to D_4$ symmetry-breaking phase
which has six-fold degenerate ground states along the $\pm\hat{x}$, $\pm\hat{y}$, and $\pm\hat{z}$ spin directions \cite{YangKG2020}.
In addition, the local magnetization, say $\langle \tilde{S}_i^z\rangle$,
shows a three-site periodicity where two of them are equal,
\begin{equation}\label{EQ:FMG1.00}
\langle \tilde{\mathbf{S}}_1\rangle = {c}\hat{z},\; \langle \tilde{\mathbf{S}}_2\rangle = {a}\hat{z},\; \langle \tilde{\mathbf{S}}_3\rangle = {a}\hat{z},
\end{equation}
in which ${a}$ and ${c}$ are the strengths of the local magnetization.
With the $U_6$ transformation shown in Eq.~\eqref{EQ:SOCD} in mind,
it is easy to check that spins in the original basis are
($|S_1^z|$, $|S_2^y|$, $|S_3^y|$; $|S_4^z|$, $|S_5^x|$, $|S_6^x|$) = $(c, a, a; c, a, a)$.
The inherent frustration in model \eqref{U1XYZ-Ham} is accidentally eliminated when $\theta = \pi/4$.
Away from this line, the interplay of bond anisotropy and competing interactions would enhance quantum fluctuations,
giving rise to new type of magnetic orderings.
It is shown in Fig.~\ref{FIG-GSPD} that there are three distinct magnetically ordered states in the middle area
where one is a collinear FM$_{U_6}$ phase while the other two are dubbed $M_1$ and $M_2$ phases.

\begin{figure}[!ht]
\centering
\includegraphics[width=0.95\columnwidth, clip]{Fig08.pdf}\\
\caption{Low-lying energy levels $E_{\upsilon}$ for $\theta = 0.30\pi$ with chain length $L = 24$.
    The ground-state degeneracy is sixlet and octuplet for $g$ = 0.80 (red circle, FM$_{U_6}$ phase) and $g$ = 0.50 (green triangle, $M_1$ phase), respectively.}
    \label{FIG-FMDegeneracy}
\end{figure}

To begin with, by reducing the strength of $g$ along the line of $\theta = 0.30\pi$,
we find that FM$_{U_6}$ phase could survive until $g\simeq\sqrt3/3$
where the ground-state degeneracy changes from sixlet to octuplet.
Figure~\ref{FIG-FMDegeneracy} show the first ($L$+2) energy levels $E_{\upsilon}$ ($\upsilon$ = $0\to25$) of a 24-site chain
at $g$ = 0.80 (red circle) and 0.50 (green triangle).
One can readily recognize that there is a energy step at the sixth (eighth) energy level for $g$ = 0.80~(0.50).
The energy barrier is $\sim\!10^{-3}$, which is several orders larger than the energy splitting within the degenerate ground states.
The energy step at the 24-th energy level is extremely steep,
which is a reminiscence of the $(L+1)$-fold degeneracy at the $SU(2)$ FM line of $\theta = \pi/4$.
We have checked the ground-state degeneracies under spin chains of $L = 48$ and 72 as well,
and the results remain unchanged in the DMRG calculation with up to $m = 4000$ states kept.

\begin{figure}[!ht]
\centering
\includegraphics[width=0.95\columnwidth, clip]{Fig09.pdf}\\
\caption{Local magnetization $\tilde{S}_i^{\alpha}$ ($\alpha$ = $x, y, z$) as a function of site index $i$ for a 24-site chain.
    (a) $\tilde{S}_i^z$ in the FM$_{U_6}$ phase at ($g = 0.80$, $\theta = 0.30\pi$).
    (b) $\tilde{S}_i^x$ (red circle), $\tilde{S}_i^y$ (green triangle), and $\tilde{S}_i^z$~(blue square) in the $M_1$ phase at ($g = 0.50$, $\theta = 0.30\pi$).}
    \label{FIG-FMPattern}
\end{figure}

We then study magnetization distributions of the symmetry breaking phases.
For the FM$_{U_6}$ phase, the spin ordering is very similar to what is shown in Eq.~\eqref{EQ:FMG1.00}
but with a six-site periodicity, see Fig.~\ref{FIG-FMPattern}(a).
It is observed that
\begin{equation}\label{EQ:FMG0.80}
\big(\langle \tilde{\mathbf{S}}_1\rangle, \langle \tilde{\mathbf{S}}_2\rangle, \langle \tilde{\mathbf{S}}_3\rangle; \langle \tilde{\mathbf{S}}_4\rangle, \langle \tilde{\mathbf{S}}_5\rangle, \langle \tilde{\mathbf{S}}_6\rangle\big) \!=\! \big({c},{b},{b};{c},{a},{a}\big)\hat{z}
\end{equation}
where ${a}$, ${b}$ and ${c}$ are magnitudes of the spin orderings along the $\hat{z}$ direction and $a, b, c \leq S$.
There is a slight difference between ${a}$ and ${b}$ when $g \neq 1$.
That is, ${b} < {a}$~(${b} > {a}$) when $g < 1$~($g > 1$).
They are equal at the isotropic case, consistent with the group-theoretical argument \cite{YangKG2020}.
As can be seen from Fig.~\ref{FIG-FMPattern}(a), ${c}$ is the smallest value of the three albeit
its difference to the penultimate value (it is ${b}$ when $g < 1$) becomes negligible as $g$ is decreased.
For the $M_1$ phase shown in Fig.~\ref{FIG-FMPattern}(b),
the $z$ component of the magnetization still shows the pattern in Eq.~\eqref{EQ:FMG0.80},
except that $c$ and $b$ are very close in value but are visibly smaller than $a$.
Most importantly, the $x$ and $y$ components of the spins in the $M_1$ phase also become nonzero
and shows the permutation relation within each even and odd sublattice.
Following the $\eta$-notation of Rousochatzakis and Perkins \cite{RousochatzakisPerkins2017}, we find that
\begin{eqnarray}\label{EQ:N1Spin135}
\langle \tilde{\mathbf{S}}_1\rangle \!=\!
    \left(
    \begin{array}{c}
        \eta_x a \\
        \eta_y b \\
        \eta_z c
    \end{array}
    \right),
\langle \tilde{\mathbf{S}}_3\rangle \!=\!
    \left(
    \begin{array}{c}
        \eta_x c \\
        \eta_y a \\
        \eta_z b
    \end{array}
    \right),
\langle \tilde{\mathbf{S}}_5\rangle \!=\!
    \left(
    \begin{array}{c}
        \eta_x b \\
        \eta_y c \\
        \eta_z a
    \end{array}
    \right)
\end{eqnarray}
and
\begin{eqnarray}\label{EQ:N1Spin246}
\langle \tilde{\mathbf{S}}_2\rangle \!=\!
    \left(
    \begin{array}{c}
        \eta_x a \\
        \eta_y c \\
        \eta_z b
    \end{array}
    \right),
\langle \tilde{\mathbf{S}}_4\rangle \!=\!
    \left(
    \begin{array}{c}
        \eta_x b \\
        \eta_y a \\
        \eta_z c
    \end{array}
    \right),
\langle \tilde{\mathbf{S}}_6\rangle \!=\!
    \left(
    \begin{array}{c}
        \eta_x c \\
        \eta_y b \\
        \eta_z a
    \end{array}
    \right).
\end{eqnarray}
Here, $a, b, c$ $( \geq 0)$ are the intensities of the magnetization,
while $\eta_{x}, \eta_y, \eta_z$ $( = \pm 1)$ are the Ising variables.
It is worth noting that $a$, $b$, and $c$ in the $M_1$ phase satisfy the restriction $\sqrt{a^2+b^2+c^2} \leq S$,
and it is quite different from these in the FM$_{U_6}$ phase (see Eq.~\eqref{EQ:FMG0.80}).
All the three $\eta$'s are free to choose any of 1 or $-1$,
accounting for the eight-fold degeneracy of the $M_1$ phase shown in Fig.~\ref{FIG-FMDegeneracy}.
In addition, by applying the inversion $U_6$ transformation,
the spins in the original basis have the following relation
$(|S_{\upsilon}^x|, |S_{\upsilon}^y|, |S_{\upsilon}^z|) = (a, b, c)$ for $1 \leq \upsilon \leq L$.
The spin structure of the $M_1$ phase is noncoplanar in the rotated basis
and it could be verified by the scalar spin chirality defined as
\begin{equation}\label{EQ:ChiDef}
\hat{\chi}_{ijk} = \tilde{\mathbf{S}}_i \cdot (\tilde{\mathbf{S}}_j \times \tilde{\mathbf{S}}_k).
\end{equation}
It is easy to check that $\hat{\chi}_{135} = \hat{\chi}_{246} \equiv \hat{\chi}_{0}$ and
\begin{align}\label{EQ:ChiNonzero}
\hat{\chi}_{0} &= \eta_x\eta_y\eta_z \big(a^3 + b^3 + c^3 - 3abc\big) \nonumber\\
&= \frac{\eta}{2} (a + b + c)\big[(a - b)^2 + (b - c)^2 + (c - a)^2\big]
\end{align}
with $\eta \equiv \eta_x\eta_y\eta_z$.
Eq.~\eqref{EQ:ChiNonzero} suggests that as long as $a$, $b$, and $c$ are not all the same,
which always holds as observed from Fig.~\ref{FIG-FMPattern}(b),
$\hat{\chi}_{0}$ will be nonzero, in line with the noncoplanar pattern of the $M_1$ phase.

\begin{figure}[!ht]
\centering
\includegraphics[width=0.95\columnwidth, clip]{Fig10.pdf}\\
\caption{The $a$ (red circle), $b$ (green triangle), and $c$ (blue square) magnetization components of the FM$_{U_6}$ and $M_1$ phases
    for $\theta = 0.30\pi$ with chain length $L = 48$.
    The middle region is the FM$_{U_6}$ phase ($0.5778 \lesssim g \lesssim 1.731$)
    while the side ones are the self-dual $M_1$ phases.}\label{FIG-FMTht030}
\end{figure}

Figure~\ref{FIG-FMTht030} shows the $(a, b, c)$ components of the magnetization along the line of $\theta = 0.30\pi$.
In the wide region of $0.5778 \lesssim g \lesssim 1.731$,
the ground state is the FM$_{U_6}$ phase where all the spins point along the $z$ direction with a almost saturated moment.
The three species $a$, $b$, and $c$ are totally different as long as $g \neq 1$.
The $M_1$ phase takes over when $0.43 < g < 0.5778$,
and magnitudes of the magnetization are suppressed approximately to $3/4$ (for $a$) or 1/2 (for $b, c$) of the saturated value.

For the $M_2$ phase, the local magnetization is fragile and we thus extract their values by calculating the spin-spin correlation functions defined as
\begin{equation}\label{EQ:CorrFunc}
C_{\upsilon}^{\alpha}(l) = \langle \tilde{S}_{\upsilon}^{\alpha} \tilde{S}_{\upsilon+l}^{\alpha} \rangle
\end{equation}
where $\alpha = x, y, z$ and $\upsilon$ is the reference site.
For simplicity we firstly consider the isotropic case ($g$ = 1) which shows a three-site periodicity in the rotated basis.
The correlators $C^{x/y/z}$ at $(g = 1.00, \theta = 0.42\pi)$ are calculated based on a 48-site chain, see Fig.~\ref{FIG-CxCyCz}.
These values are very stable when the site distance $l$ is larger than 10,
and we estimate the local magnetization as $\langle \tilde{S}_{\upsilon}^{\alpha} \rangle = \sqrt{C_{\upsilon}^{\alpha}(L/2)}$ with $L = 48$.
The local magnetization $\langle \tilde{\mathbf{S}}_{\upsilon} \rangle = (\langle \tilde{S}_{\upsilon}^x\rangle, \langle \tilde{S}_{\upsilon}^y\rangle, \langle \tilde{S}_{\upsilon}^z\rangle)^T$
within the 3-site unit cell is
\begin{eqnarray*}
\left(\langle \tilde{\mathbf{S}}_1 \rangle,\langle \tilde{\mathbf{S}}_2 \rangle,\langle \tilde{\mathbf{S}}_3 \rangle\right) \!=\!
    \left(
    \left[\!
    \begin{array}{c}
        0.214 \\
        0.214 \\
        0.103
    \end{array}
    \!\right]\!,
    \left[\!
    \begin{array}{c}
        0.214 \\
        0.103 \\
        0.214
    \end{array}
    \!\right]\!,
    \left[\!
    \begin{array}{c}
        0.103 \\
        0.214 \\
        0.214
    \end{array}
    \!\right]
    \right)\!.
\end{eqnarray*}
When away from the isotropic line where $g = 1$, the unit cell is doubled
and there is a same magnetization distribution to that one shown in Eq.~\eqref{EQ:N1Spin135} and Eq.~\eqref{EQ:N1Spin246}.
It should be noted that both $M_1$ and $M_2$ phases are eight-fold degenerate and their difference lies in the relative values among $a$, $b$, and $c$.
For $M_1$ phase we have $c \simeq b < a$ (when $g < 1$) or $c \simeq a < b$ (when $g > 1$),
while for $M_2$ phase $c$ is much smaller than $a, b$.

\begin{figure}[!ht]
\centering
\includegraphics[width=0.95\columnwidth, clip]{Fig11.pdf}\\
\caption{Spin correlators $C_{\upsilon}^{\alpha}$ ($\alpha$ = $x, y, z$) of a segment $l$ for a 48-site chain.
    The reference site $\upsilon$ could be I (site 1), II (site 2), III (site 3).
    (a), (b), and (c) are for $C^{x}$, $C^{y}$, and $C^{z}$ in the $M_2$ phase at ($g = 1.00$, $\theta = 0.42\pi$), respectively.}
    \label{FIG-CxCyCz}
\end{figure}

\subsection{Magnetic orderings of $M_1$ and $M_2$ phases}
This section is devoted to study the transitions to $M_1$ and $M_2$ phases.
We begin by considering the transitions to $M_1$ phase along the path of $g = 0.5$.
The SOPs of the Haldane-type $\mathcal{O}_H = \mathcal{O}_e^z$ (via Eq.~\eqref{U1XYZ-Ham}) and
the Kitaev-type $\mathcal{O}_K = \mathcal{O}_K^{x}$ (via Eq.~\eqref{KtvGam-Ham}) are plotted in Fig.~\ref{FIG-OPWpJ2050}(a).
At $\theta_H \approx 0.215\pi$ and $\theta_K \approx 0.383\pi$, the two SOPs are discontinuous,
indicating of first-order transitions between the EH~($A_x$) phase and the intermediate $M_1$ phase.
The order parameter $O_{M_1}$ is shown in Fig.~\ref{FIG-OPWpJ2050}(b).
Here, only $c$ flavor of $\langle S_i^z\rangle$ is chosen for the sake of brevity.
We find that it is very robust with a negligible finite-size effect.
In addition, there is also a nonvanishing correlation of a rank-2 spin-nematic (SN) ordering defined in Eq.~\eqref{EQ:SNOrder} (not shown).
Notably, when crossing the line of $\theta = \pi/4$ where the ground state is the $SU(2)$ FM phase,
$O_{M_1}$ has a discontinuity because of the inherent difference of the spin orientations.

\begin{figure}[!ht]
\centering
\includegraphics[width=0.95\columnwidth, clip]{Fig12.pdf}\\
\caption{(a) SOPs of the Haldane-type $\mathcal{O}_H$~(filled symbols) and Kitaev-type $\mathcal{O}_K$~(open symbols)
    for $g = 0.5$ with chain length $L = 48$~(green triangle) and 72~(blue square).
    (b) Order parameter of the $M_1$ phase $O_{M_1}$~($c$ component only) and (c) Flux density $\langle W_p\rangle$
    in the same region as (a).}
    \label{FIG-OPWpJ2050}
\end{figure}

Meanwhile, it is appealing to know how the flux density $\langle\hat{W}_p\rangle$,
\begin{equation}\label{EQ:Wp}
\hat{W}_p = 2^6 \tilde{S}_1^z \tilde{S}_2^y \tilde{S}_3^x \tilde{S}_4^z \tilde{S}_5^y \tilde{S}_6^x,
\end{equation}
evolves in each different phase. Similar to the two-dimensional counterpart \cite{Kitaev2006},
the quantity in Eq.~\eqref{EQ:Wp} is the product of spin operators on consecutive overhanging bonds
within the six-site unit cell~(see Fig.~\ref{FIG-Bond}(b)).
In Fig.~\ref{FIG-OPWpJ2050}(c), we plot the flux density $\langle\hat{W}_p\rangle$ versus $\theta$
in the whole region of $\theta \in [0, \pi/2]$.
It is clearly shown that $\langle\hat{W}_p\rangle$ is zero in the EH phase.
In the $M_1$ phase, $\langle\hat{W}_p\rangle$ starts from a nonzero value and goes up with the increasing of $\theta$
except for $\theta = \pi/4$ where $\langle\hat{W}_p\rangle \simeq 0$.
In the $A_x$ phase, however, $\langle\hat{W}_p\rangle$ decreases from 0.17 or so
and does not stop dropping until $\theta \simeq 0.50$ where $\langle\hat{W}_p\rangle < 0$.
Without doubt, the flux density $\langle\hat{W}_p\rangle$ shows a crucial difference among the three distinct phases.
The jump and kink are excellent probes for phase transitions involved.

We now turn to the transition around the $M_2$ phase.
It is shown in the isotropic $K$-$\Gamma$ chain that there is an intermediate phase when $0.40\pi \lesssim \theta \lesssim 0.466\pi$ \cite{YangSN2021}.
This phase is now recognized as the $M_2$ phase with a nonzero magnetization,
and it could survive against small anisotropy where $|g-1| \ll 1$, and then gives way to the conventional collinear FM$_{U_6}$ phase.
Although both phases are magnetically ordered,
we appreciate the rank-2 SN ordering as a sensitive probe to capture the phase transition.
It is finite in the $M_2$ phase while vanishes in the FM$_{U_6}$ phase since the latter is collinear.
The SN correlation function in the long-distance limit is known as \cite{YangSN2021}
\begin{equation}\label{EQ:SNOrder}
O_{SN}^2 = \lim_{|j-i|\gg1}\langle \tilde{S}_{3i+1}^{y}\tilde{S}_{3i+2}^{z} \cdot \tilde{S}_{3j+1}^{y}\tilde{S}_{3j+2}^{z} \rangle.
\end{equation}
Since the species of spins in the correlator of Eq.~\eqref{EQ:SNOrder} always come in pair,
the sign of the correlation remains uninfluenced.
In addition, on can infer from Eq.~\eqref{EQ:N1Spin135} and Eq.~\eqref{EQ:N1Spin246}
that there is no difference of the magnetization among the degenerate ground states,
indicating that all of them will produce a same value of the order parameter.
As a result, we do not need to distinguish these states and
the value of the order parameter could be safely obtained via Eq.~\eqref{EQ:SNOrder}.

\begin{figure}[!ht]
\centering
\includegraphics[width=0.95\columnwidth, clip]{Fig13.pdf}\\
\caption{(a) Correlations of the SN order (see Eq.~\eqref{EQ:SNOrder}) for chain length $L = 48, 96, 144$, and 192.
    The selected point at ($g = 1.00$, $\theta = 0.42\pi$) is deep in the $M_2$ phase.
    Inset: Extrapolation of the SN order to infinite-size system.
    (b) Order parameter $O_{SN}$ of the $M_2$ phase for $\theta = 0.42\pi$
    with chain length $L = 48$~(red circle) and 72~(blue square).}
    \label{FIG-SN}
\end{figure}

Figure~\ref{FIG-SN} shows the order parameter $O_{SN}$ for $\theta = 0.42\pi$.
To check for the finite-size effect, we consider the isotropic case $g = 1$ and
calculate the correlation function in Eq.~\eqref{EQ:SNOrder}
with $i = 0$ and $j$ =0, 1, 2, $\cdots$, $L/6$, see Fig.~~\ref{FIG-SN}(a).
It is found that the correlators saturate to a finite value after several times of oscillation.
In what follows we shall define $O_{SN} \equiv O_{SN}(i=0, j=L/6)$.
The inset shows the extrapolation of the SN order parameter for chain length $L$ = 48, 72, 96, 144 and 192,
from which we can clearly find that $O_{SN}$ is very robust against $L$.
We then extend the calculation of $O_{SN}$ for $0.8 \leq g \leq 1.2$ and the results are summarized in Fig.~~\ref{FIG-SN}(b).
Deep into the $M_2$ phase, $O_{SN}$ is very stable although there is a modest suppression near the boundaries.
The transitions between the $M_2$ phase and the collinear FM$_{U_6}$ phase are accompanied by the jumps of $O_{SN}$,
from which the transition points are determined as $g_t \approx$ 0.945 and 1.065, respectively.
It is worth mentioning that the transition points satisfy the self-dual relation shown in Eq.~\eqref{SelfDual}
since they are related as $ 0.945 \simeq 1/1.065$.

Empirically, the ground-state energy per-site of the $FM_{U_6}$ phase is given by
\begin{equation}\label{EQ:EgFMU6}
e_g^{\textrm{FM}_{U_6}} = -\frac{1}{6}\big[K(a^2+gb^2) + 2c\Gamma(b+ga)\big]
\end{equation}
where $a$, $b$, $c$ are almost saturated (see Fig.~\ref{FIG-FMPattern}(a) and Fig.~\ref{FIG-FMTht030}).
For example, at the hidden $SU(2)$ FM point where $K = \Gamma$ and $g = 1$,
we have $a = b = c = 1/2$ and the energy inferred from Eq.~\eqref{EQ:EgFMU6} is $-K/4$,
consistent with the analytical result \cite{KomaNach1997}.
For the $M_2$ (and also $M_1$) phase, the energy displays a very similar form except that
the prefactor ($1/6$) in Eq.~\eqref{EQ:EgFMU6} should be $1/2$.
In addition, $(a, b, c)$ subject to the constraint $\overline{M} \equiv \sqrt{a^2 + b^2 + c^2} \leq S$.
However, the total magnetization $\overline{M}$ of the $M_2$ phase is far from saturated,
and it is only 0.320 with $(a, b, c) = (0.214, 0.214, 0.103)$ at $(g = 1.00, \theta = 0.42\pi)$.
By adding the bond alternation with $g \neq 1$,
there is a slight enhancement of $\overline{M}$, lowering the ground-state energy and thus opening a finite region of $M_2$ phase.

\section{Transitions around the FM Kitaev limit}\label{SEC:FMKtv}

So far, we have mainly concentrated on the right panel of the phase diagram shown in Fig.~\ref{FIG-GSPD}.
Phases in the left panel could be obtained from the right part after a mirror operation.
However, little is known about the transition types of the adjacent phases near the axis of symmetry.
After an inspection of the first-order energy derivative $\partial e_g/\partial\theta$ along the line of $g = 1$,
Yang \textit{et. al.} claimed that the transition at the FM Kitaev point is of first order \cite{YangKG2020}.
A variational Monte Carlo study, amazingly, suggests that the $Z_2$ QSL at that point
could survive up to a small $\Gamma$-interaction \cite{WangLiu2020}.
Herein, we find that the FM Kitaev point is a confluence point of two transition lines,
i.e., the $A_x$-$A_y$ and the LL-LL$'$ transition lines.
It is thus a multicritical point which accounts for the difficulty in determining the nature of transition
(for extended discussion, see Sec. III in the Supplemental Material \cite{SuppMat}).
By virtue of an efficient bond-reversal method \cite{LuoQPT2019},
we argue in the following that the aforesaid topological QPT between the two LLs is continuous.
Nevertheless, away from the symmetric line of $g = 1$, transitions at $\theta = -\pi/2$ are of first order
without closing the gap at the transition points.

Here we illustrate how to use the bond-reversal method to determine the transition type around the FM Kitaev point.
By tuning $\Gamma$ from negative to positive,
the ground-state energy $e_g$ must be symmetric with respect to the $\Gamma = 0$ line
(i.e., $\theta = -\pi/2$) due to the symmetry relation of Eq.~\eqref{GammaSym}.
When $\theta$ is slightly away from $-\pi/2$,
the sign of $\Gamma$-interaction is changed,
and local expectations of $\Gamma_x = \langle S_1^yS_2^z + S_1^zS_2^y\rangle$
and $\Gamma_y = \langle S_2^zS_3^x + S_2^xS_3^z\rangle$ in the \textit{original} basis must be reversed.
In this regard, we thus define the difference of bond strength~(DBS) $\mathcal{D}$ as
\begin{equation}\label{EQ:DBS}
\mathcal{D} = \frac{1}{2}(\Gamma_x + \Gamma_y).
\end{equation}
The DBS $\mathcal{D}$ is a sensitive probe for a first-order QPT because it have a jump at the transition point.
Oppositely, there is a continuous QPT if $\mathcal{D}$ is smoothly changed \cite{LuoQPT2019}.
We note in passing that, physically, $\mathcal{D}$ is equivalent to
the first-order derivative of ground-state energy $\partial e_g/\partial\theta$ when $g = 1$.
However, the energy derivative depends on the increment $\delta\theta$ which may cause artificial oscillation.
In this sense, the DBS $\mathcal{D}$ is obviously superior and is more reliable.

\begin{figure}[!ht]
\centering
\includegraphics[width=0.95\columnwidth, clip]{Fig14.pdf}\\
\caption{(a) DBS $\mathcal{D}$ of the $g_x$-$g_y$ $K$-$\Gamma$ chain for $g$ = 1.0 with
    chain length $L = 24$~(red circle), 48~(green triangle) and 72~(blue square).
    Inset: Zoom in of the DBS near the FM Kitaev limit with $\Gamma > 0$.
    (b) DBS $\mathcal{D}$ for $g$ = 0.5~(red circle), 1.5~(green triangle) and 2.0~(blue square) with chain length $L = 24$.}
    \label{FIG-KtvDBS}
\end{figure}

Figure~\ref{FIG-KtvDBS}(a) shows the DBS $\mathcal{D}$ of the LL--LL$'$ transition when $g$ = 1.0.
It is rather smooth without any jump in a wide region of $-0.55\pi \leq \theta \leq -0.45\pi$.
The size-dependent behavior is insignificant except for a narrow slit near the FM Kitaev point.
As can be seen from the inset of Fig.~\ref{FIG-KtvDBS}(a),
the DBS $\mathcal{D}$ shows a well-controlled scaling behavior for different chain length $L$,
and it does not have a jump although its slope becomes sharp as $L$ increased,
indicating of a multicritical behavior.
We recall that such a multicritical point is analogous to the one existing in a 1D transverse $XY$ spin chain
which owns an intersection point of two transition lines of different universality classes
\cite{LiebSM1961,DamleSach1996,YangZS2019}.
When shifting away from $g$ = 1.0, the ground state is occupied by gapped EH or OH phase.
The DBS $\mathcal{D}$ for three selected EH--EH$'$ transition~($g$ = 0.5) and OH--OH$'$ transitions~($g$ = 1.5 and 2.0)
are shown in Fig.~\ref{FIG-KtvDBS}(b).
For all the cases there are appreciable jumps of $\mathcal{D}$ at $\theta = -\pi/2$,
representing the hallmark character of the first-order QPT.

\section{Summary and Discussion}\label{SEC:Conclusion}
We have numerically studied the phases and phase transitions in a bond-alternating spin-$1/2$ $K$-$\Gamma$ chain,
which is an excellent platform to reveal many aspects of one-dimensional quantum magnetism.
By calculating various conventional symmetry-breaking order parameters and nonlocal SOPs,
we unveil a rich quantum phase diagram which contains seven different phases.
Near the AFM Kitaev spin chain limit, there is a critical segment with a macroscopic ground-state degeneracy.
It is unstable against bond alternation, resulting in two gapped disordered $A_x$ and $A_y$ phases characterized by nonlocal SOPs.
The $A_x$--$A_y$ topological QPT falls in the Ising universality class with a central charge $c = 1/2$.
On the other hand, starting from the FM Kitaev spin chain limit by increasing the $\Gamma$ interaction,
there are EH and OH phases in the inner circle ($g < 1$) and outer circle ($g < 1$), respectively.
The EH--OH transition is determined by the SOPs which vanish algebraically at the transition boundary.
It could also be captured by the entanglement gap which undergoes a sign change when crossing the critical point.
This transition belongs to the Gaussian universality class with a central charge $c = 1$,
identical to that of the bond-alternating spin-$1/2$ AFM Heisenberg chain.
The FM Kitaev point is recognized as a multicritical point converging several different phases.
In addition, there are also three distinct magnetically ordered states, named FM$_{U_6}$, $M_1$, and $M_2$ phases,
in the presence of AFM Kitaev interaction.
The FM$_{U_6}$ phase has a six-fold degeneracy and is situated in a wide region around the isotropic line of $g = 1$.
The $M_1$ and $M_2$ phases are highly spatially modulated and could have a rank-2 spin-nematic ordering.

The isotropic $\Gamma$-chain is conformally invariant with a central charge $c = 1$.
While its ground-state energy smoothly varies with the chain length under OBC,
it surprisingly shows an unconventional six-site periodicity under PBC.
We remark that this phenomenon has a profound relation to the abnormal energy scaling in two dimensional honeycomb lattice \cite{LuoZhaoKeeWang2019}.
Given that there is an emergent $SU(2)$ symmetry at this $\Gamma$ limit \cite{YangKG2020},
we conjecture that the versatile Bethe ansatz may be capable to give an exact solution of the isotropic $\Gamma$-chain.

In closing, our work demonstrates the essential role played by the bond alternation in enriching the underlying phase diagram.
The bond alternation is a relevant perturbation to either open up the energy gap or rearrange the distribution of magnetization,
leaving the possibility for the emergence of novel phases.
Our study also highlights the richness of Kitaev systems with AFM exchange interaction.
Although the $K$-$\Gamma$ model is widely recognized as a cornerstone
to describe candidate Kitaev materials like $\alpha$-RuCl$_3$,
much less attention has been paid to $K > 0$ as Kitaev interaction is likely negative in these materials.
A theoretical proposal for the AFM Kitaev interaction in $f$-electron based magnets has been proposed \cite{JangSKM2019}.
Our study thus corroborates a new direction to hunt for exotic phases in a less explored area.

\begin{acknowledgments}
Q.L. would like to thank W.-L. You for fruitful discussions on the quantum compass model.
X.W. was supported by the National Program on Key Research Project (Grant No. 2016YFA0300501)
and the National Natural Science Foundation of China (Grant No. 11974244),
He also acknowledged the support from a Shanghai talent program.
J.Z. was supported by the National Natural Science Foundation of China (Grant No. 11874188).
H.-Y.K. was supported by the NSERC Discovery Grant No. 06089-2016,
the Centre for Quantum Materials at the University of Toronto,
the Canadian Institute for Advanced Research,
and also a funding from the Canada Research Chairs Program.
Computations at early stages were performed on the Tianhe-2JK at the Beijing Computational Science Research Center (CSRC).
Computations were mostly performed on the Niagara supercomputer at the SciNet HPC Consortium.
SciNet is funded by: the Canada Foundation for Innovation under the auspices of Compute Canada;
the Government of Ontario; Ontario Research Fund - Research Excellence; and the University of Toronto.
\end{acknowledgments}