\section{Introduction}

Chirality, representing right- or left-handedness, generates attractive physical properties in material science.
The importance of chirality has been recognized by a relationship between chiral structures and chirality-induced phenomena,
as exemplified by natural optical activity in optical science \cite{Arago1811,Biot1812,Pasteur1,Pasteur2} and chiral helimagnetic ordering,
triggered by an antisymmetric exchange (Dzyaloshinskii-Moriya) interaction \cite{Dzyaloshinskii1958,Moriya1960},
in chiral magnetism \cite{TSi2,Ishida1985,Kishine2015,Togawa2016}.
Indeed, the sign of Dzyaloshinskii-Moriya interaction depends on the handedness of a crystalline structure of chiral helimagnetic materials \cite{Togawa2012}.
In this connection, the importance of synthesizing inorganic enantiopure crystals has been emphasized in material science \cite{Kousaka2017}, which will advance our understanding of a wide range of chirality-induced phenomena in inorganic chiral materials such as magnetochiral dichroism (MChD) \cite{Kibayashi2014,Okumura2015,Nakagawa2017},
nonreciprocal electrical transport \cite{Rikken2001,Yokouchi2017,Aoki2019} and chirality-induced spin selectivity (CISS) effect \cite{Inui2020,Shiota2021,Shishido2021}.

The CISS effect has attracted much attention since a controllability of spin polarization is given in nonmagnetic materials with chiral structures.
The amplitude of the spin polarization becomes as large as that for ferromagnetic materials, while its sign is dependent on the chirality of materials.
It was first observed in chiral molecules \cite{Gohler2011,Xie2011}, and quite recently found in chiral inorganic crystals \cite{Inui2020,Shiota2021,Shishido2021}.
In the latter case, when a charge current flows through the crystal,
the spin-polarized state appears as a linear response of the current-voltage characteristics at room temperature without magnetic field.
The sign of the linear slope is connected to the crystal handedness. In this respect, the CISS signal has a great potential to probe the crystal chirality.
Furthermore, the spin polarization emerges in a length scale over micrometers \cite{Inui2020,Shiota2021} and even over millimeters \cite{Shishido2021} in nonlocal measurements with inorganic chiral crystals.

Note that it is still difficult to evaluate and control the chirality of crystal structures in inorganic materials.
The crystal chirality is generally determined by an absolute structure analysis using single crystal X-ray diffraction measurements.
One of the disadvantages of the X-ray method is that the sample size is limited to a scale of sub-millimeter or smaller because the structure analysis is significantly affected by a large X-ray absorption.
When the size of specimens is in a millimeter scale or larger, such crystals need to be smashed into small portions of sub-millimeter.
In this respect, the single crystal X-ray diffraction measurements correspond to a destructive method for chirality evaluation and probe the chirality in a part of the specimens.

As for controlling the chirality, inorganic chiral compounds normally crystallize racemic twinned crystals, which contain right- and left-handed chirality domains.
One possible method of obtaining a homochiral crystal with the selected handedness is a spontaneous crystallization with stirring solution \cite{Kondepudi1990}.
This technique is applied to a synthesis of water-soluble compounds.
Indeed, the mono-chiral crystals are successfully obtained in NaClO$_3$ \cite{Kondepudi1990} and CsCuCl$_3$ \cite{Kousaka2017,Kousaka2014}, leading to the observation of enhanced MChD signals \cite{Nakagawa2017}.
Obviously, these studies indicate that the enantiopure crystals play a key role in inducing nontrivial physical properties coupled with the handedness of the crystals.

Another example of chirality-controlled synthesis is found in intermetallic monosilicide compounds $T$Si ($T:$ a transition metal) with a B20 type chiral cubic crystal structure.
$T$Si is known as one of the representative chiral magnetic compounds, which exhibit a formation of chiral helimagnetic structure \cite{TSi2,Ishida1985,TSi3} and Skyrmion lattice \cite{Kadowaki1982,Bogdanov1989,Muhlbauer2009,Neubauer2009,Yu2010}.

It is known that single crystals of $T$Si are obtained with the single handedness by use of conventional crystal growth techniques and the crystal chirality depends on the element of $T$ site \cite{TSi2,TSi3,TSi4,Siegfried2015}.
The right-handed crystals are stabilized when $T =$ Fe, while the left-handed ones are grown when $T =$ Mn and Co.
Note that a seed crystal is used in a chirality-selected crystal growth.
Interestingly, the handedness of the grown crystal inherits from the one of the seed crystal.
Namely, when using a seed crystal with the $T$ element appropriate for the desired handedness, the single crystal with the selected handedness is obtained \cite{TSi5,Kousaka2022,Bauer2022}.
For instance, a right-handed CoSi can be grown by using a right-handed FeSi as a seed crystal.
Such ideas of the chirality-controlled crystal growth were indeed demonstrated in a Czochralski method \cite{TSi5} and in a floating zone (FZ) method \cite{Kousaka2022,Bauer2022}.

Intermetallic disillicide compounds $T$Si$_{2}$ with a C40 type chiral hexagonal crystal structure are the next target materials for the chirality-controlled crytal growth.
Recently, $T$Si$_{2}$ is found to show the CISS effect \cite{Shiota2021,Shishido2021}.
The direction of the spin polarization in the right-handed NbSi$_{2}$ is opposite to that in the left-handed TaSi$_{2}$ \cite{Shiota2021}.
These crystals exhibit the spin polarization in a length scale over millimeters \cite{Shishido2021}.

The chiral crystals of centimeters in length are necessary in order to investigate the length scale over which the spin polarization is kept in the chiral crystal.
Moreover, it would be interesting to examine the spin polarization between the same compound $T$Si$_{2}$ of different chirality (the same $T$ disillicide with right-handed and left-handed structures). These studies may reveal a correlation between the crystal chirality and spin polarization via the CISS effect.

In this study, we demonstrate the chirality-controlled crystal growth of NbSi$_{2}$ and TaSi$_{2}$ by using a laser-diode-heated floating zone (LDFZ) method.
With the present technique, we obtained the chirality-selected crystals of NbSi$_{2}$ and TaSi$_{2}$ over centimeters long.
The CISS experiments with NbSi$_{2}$ of the selected chirality unveiled a correlation between the crystal chirality and CISS signals.
Our study exhibits that the CISS measurements work as a non-destructive technique for the chirality evaluation even in centimeter-sized specimens.

\section{Crystal structure and chirality of disillicides}

The disillicides $T$Si$_{2}$ belong to the space group $P6_{2}22$ and $P6_{4}22$, which forms the right- and left-handed crystal structures, respectively, as shown in Fig.~\ref{f-TSi-str}. $T$ and Si atoms occupy Wyckoff $3c$ and $6i$ positions as expressed by $(1/2,0,0)$ and $(x,2x,0)$, respectively.
The parameter $x$ for Si slightly deviates from a value of $1/6$.
Accurately speaking, the $T$ and Si atoms can also locate at Wyckoff $3d$ and $6j$ positions, respectively.
The difference between a pair of the $3c$ and $6i$ positions and that of the $3d$ and $6j$ positions is given by a displacement vector of $(0,0,1/2)$.
When the crystal structures are viewed along the $c$-axis, the right-handed double helices of Si atoms are found in the right-handed crystal and vice versa in the left-handed crystal.

The crystal chirality of the disillicide compounds is examined by an absolute structure analysis using single crystal X-ray diffraction, clarifying that the crystal chirality of $T$Si$_2$ depends on the elements of $T$ site \cite{Sakamoto2005, Onuki2014}.
The right-handed crystals are stabilized when $T =$ Nb, while the left-handed ones are grown when $T =$ V, Cr, and Ta.

\begin{figure}[tb]
\begin{center}
\includegraphics[width=8cm]{fig1_TSi2_str.eps}
\end{center}
\caption{Crystal structures of (a) right-handed and (b) left-handed $T$Si$_{2}$ ($T:$ a transition metal).
The lower schematics show the (001) projection of each structure.
Large green and small blue balls represent $T$ and Si atoms, respectively.
The black arrows are given as an eye guide to recognize the sense of screw alignments of the Si atoms.}
\label{f-TSi-str}
\end{figure}

\begin{figure*}[tb]
\begin{center}
\includegraphics[width=15.0cm]{fig2_FZ_rev7.eps}
\end{center}
\caption{
Schematic drawings of a crystal growth using a laser-diode-heated floating zone (LDFZ) method.
(a) A conventional LDFZ method using a feed rod made of the single compound.
(b) A LDFZ method using a composition-gradient feed rod, which consists of two regions with different compounds.
A combination of the $T$Si$_{2}$ compounds used in the present study is listed.
(c) and (d) The procedure of the LDFZ crystal growth.
The crystal growth starts with the region A, transfers into the region B, and finishes with the edge of the feed rod. As a result, the chirality of the region A is inherited into that of the region B. As for the composition-gradient feed rod (\#3) shown in Fig. \ref{f-FZ}(b), the chirality of TaSi$_{2}$ in the region B of the grown rod is determined by that of NbSi$_{2}$ (left-handed) in the region A, even though TaSi$_{2}$ prefers right-handed chirality in the conventional crystal growth. Various combinations of C40 materials in the feed rod enables the control of chirality of the grown crystals.
}
\label{f-FZ}
\end{figure*}

\begin{figure}[tb]
\begin{center}
\includegraphics[width=8cm]{fig3_chirality-probe_rev5.eps}
\end{center}
\caption{
Schematic drawings of probing the chirality of materials using (a) a single crystal X-ray diffraction method and (b) electrical transport measurements via the CISS effect.
Note that the single crystal X-ray diffraction experiments were carried out by rotating the specimen with $0.5^\circ$ oscillation per frame. Over one thousand frames were collected for a specimen in order to observe thousands of Bragg reflection spots.}
\label{f-chiral-probe}
\end{figure}

\section{Experimental methods}

Single crystals of $T$Si$_{2}$ are obtained by the LDFZ method, as shown in Fig.~\ref{f-FZ}.
The FZ method is a crystal growth technique to obtain rod-shaped single crystals with several millimeters in diameter and several centimeters in length.
Generally speaking, the feed and seed rods for the crystal growth consist of the same compound.
In this study, the feed rods contain only NbSi$_{2}$ or TaSi$_{2}$ for the conventional FZ crystal growth, as shown in Fig.~\ref{f-FZ}(a).
A molten zone is formed by connecting the edges of both the rods that are melted by optical radiation sources.
In the growth process, the rods are pulled down through the optical radiated area.
When the melting part moves out of the radiation, a single crystal is grown at the boundary between the melting and solidifying regions.

However, in some compounds, it is difficult to obtain single crystals via the conventional FZ crystal growth because of a very slow growth rate and/or instability of the molten zone.
Such problems can be solved by preparing a composition-gradient feed rod in the FZ crystal growth \cite{Ikeuchi2003, Enoki2010, Enoki2013}.
This advanced technique was successfully demonstrated for a growth of superconducting copper oxide compounds La$_{2-x}$Sr$_{x}$CuO$_{4}$ (La214) \cite{Ikeuchi2003} and Bi$_{2+x}$Sr$_{2-x}$CuO$_{6}$ \cite{Enoki2010, Enoki2013}.
In the case of La214, the composition-gradient feed rod had a step-like Sr-concentration gradient along the rod. Eventually, the obtained single crystal had an alternation of Sr concentration along the grown direction of the crystal.

Similar technique can be used in a synthesis of chiral inorganic compounds $T$Si.
The LDFZ method using the composition-gradient feed rod of $T$Si sustained the chirality inherited from the seed rod. The grown MnSi and CoSi crystals exhibited the chiral structures, which had the chirality opposite to that obtained by the conventional LDFZ method \cite{Kousaka2022}.

In this study, we performed the crystal growth of $T$Si$_{2}$ by the LDFZ using a composition-gradient feed rod consisting of TaSi$_{2}$ and NbSi$_{2}$, as shown in Fig.~\ref{f-FZ}(b).
To prepare the feed and seed rods, the rod-shaped polycrystalline samples of $T$Si$_{2}$ were synthesized by solid-phase reaction.
A stoichiometric mixture of $T$ and Si in the molar ratio of 1 : 2 was heated at $1100 ^\circ$C in a vacuumed quartz tube for two days.
Single crystals of $T$Si$_{2}$ were grown by the LDFZ furnace (equipped with five 200 W lasers, Crystal Systems Corporation).
The crystal growth was performed in an Ar atmosphere with a flowing rate of 3 $l$$/$min.
To obtain single crystals, traveling speeds of the feed and seed rods were set to be the same, and the maximum speed was 10 mm/hour.
To make a stable molten zone, the feed and seed rods were rotated in the opposite directions at 4 and 5 rpm, respectively.

The crystallographic orientation of the rod-shaped crystals was evaluated by a Laue X-ray backscattering method
using a tungsten X-ray tube and an imaging plate (IPX-LC, IPX Corporation).
The crystal chirality was probed by the absolute structure analysis using a single crystal X-ray diffractometer using a molybdenum X-ray tube (XtaLAB mini II, Rigaku Corporation), as shown in Fig.~\ref{f-chiral-probe}.
In addition, the CISS signals were obtained using electrical transport measurements via an inverse spin Hall effect.

The single crystal X-ray diffraction experiments were performed by rotating the specimen with $0.5^\circ$ oscillation per frame in the $\omega$-axis.
More than one thousand frames were collected to deduct integrated intensities of thousands of the Bragg reflection spots,
which are assigned to the Miller indices $(hkl)$ in a reciprocal hemisphere.
As shown in Fig.~\ref{f-chiral-probe}(a), the specimen size must be smaller than the beam diameter and thus is limited to the product $\mu t$, where $\mu$ and $t$ are X-ray absorption coefficient and path length, respectively.
Resultantly, the absolute structure analysis is applicable only to the specimens of a sub-millimeter scale or smaller one.

Regarding the CISS effect, the direction of spin-polarization depends on the crystal chirality of materials.
In inorganic chiral crystals, the spin-polarization has been probed by an inverse spin Hall effect, where a tungsten (W) or platinum (Pt) electrode, made on the rectangular-shaped specimen, is used for the spin detection, as shown in Fig.~\ref{f-chiral-probe}(b).
There is no limitation on sample size in the CISS measurements. Namely, it ranges from micrometers or millimeters in length \cite{Inui2020,Shiota2021,Shishido2021}.
When a longitudinal electrical current $I_{xx}$ is applied to the $c$-axis of chiral crystal, the spin-polarization occurs in the crystal because of the CISS effect.
The electrode detects the difference of the spin-dependent chemical potential difference between the specimen and the electrode,
and converts the spin polarization into the transverse voltage $V_{xy}$ via the inverse spin Hall effect.
Therefore, the crystal chirality of the specimen can be determined by observing the sign of the slope of the transverse voltage $V_{xy}$.

\section{Experimental results}

\subsection{Chirality-controlled crystal growth}

First, the crystal growth was carried out by the conventional LDFZ method.
In this case, the seed and feed rods were made of the single compound of C40 materials, as shown in Fig.~\ref{f-FZ}(a).
The rod-shaped crystals of NbSi$_{2}$ (\#1) and TaSi$_{2}$ (\#2) were grown by the LDFZ method.

Laue X-ray backscattering experiments showed that the crystallographic $c$-axis of the NbSi$_{2}$ and TaSi$_{2}$ crystals oriented in the growth direction of the obtained rod with a slight deviation of a few degrees.
For chirality determination, a part of the edge of the crystal rod was smashed into some submillimeter pieces and these pieces were evaluated by means of single crystal X-ray diffraction.
The examined pieces of NbSi$_{2}$ and TaSi$_{2}$ crystals exhibited tiny mosaicity and were available for the absolute structure analysis.
With applying an empirical absorption correction to the observed intensity, the absolute structures were solved by the direct method and refined using the SHELXL software package \cite{Shelx}.

In the case of non-centrosymmetric crystal structures, the intensities of a pair of reflections at $(h k l)$ and $(\bar{h} \bar{k} \bar{l})$, termed Bijvoet pairs \cite{Bijvoet1954}, are not equivalent because of anomalous scattering.
The Flack parameter $x$~\cite{Flack1983, Berardinelli1985, Flack1999, Parson2013} is one of the fitting parameters in the structure analysis, and expressed as $I_{obs}(h,k,l)=(1-x)I_{cal}(h, k, l)+xI_{cal}(\bar{h}, \bar{k}, \bar{l})$, where $I_{obs}$ and $I_{cal}$ are the observed and calculated intensities, respectively.
The $I_{cal}(h,k,l)$ corresponds to the intensity calculated from the assumed chiral crystal structure in the analysis, while the $I_{cal}(\bar{h}, \bar{k}, \bar{l})$ is that for the opposite chiral structure.
When the deducted Flack parameter $x$ is refined to be zero, the assumed chiral crystal structure is correct.
When $x$ becomes one, the assumed chiral crystal structure has to be flipped into that with the opposite chirality.

Therefore, the handedness of the crystal structure was deducted by the Flack parameter $x$, which was obtained in the refined process of the analysis.
In the case of the right-handed crystal structure, $x$ becomes zero when the space group of the assumed structure is $P6_{2}22$.
In the case of the left-handed crystal structure, $x$ becomes zero for the assumed structure with the space group of $P6_{4}22$.

The results of the absolute structure analysis are summarized in Table \ref{tbl-NbTaSi2-conventional}.
The handedness of the crystals was left-handed for TaSi$_{2}$, while right-handed for NbSi$_{2}$.
These results are consistent with those reported in the literature \cite{Sakamoto2005}.
The lattice constants $a$ and $c$ of NbSi$_{2}$ are larger (below 0.5 percent) than those of TaSi$_{2}$.
The lattice constants decrease with increasing the atomic number of $T$.
Namely, the values of $a$ and $c$ are proportional to the atomic radius of $T$ and can be used for determining $T$ atoms of the synthesized crystals.

\begin{table}[htbp]
      \caption{Crystal structure analysis of $T$Si$_{2}$, obtained by the LDFZ method using the feed rod made of the single composition (TaSi$_2$ for the rod \#1 and NbSi$_2$ for the rod \#2).}
      \label{tbl-NbTaSi2-conventional}

  \begin{center}

      \begin{tabular}{rrr}
      \hline
      \multicolumn{1}{c}{} & \multicolumn{1}{c}{TaSi$_{2}$} & \multicolumn{1}{c}{NbSi$_{2}$}\\
      \hline
      \multicolumn{1}{c}{Space group} & \multicolumn{1}{c}{$P6_{4}22$} & \multicolumn{1}{c}{$P6_{2}22$}\\
      \multicolumn{1}{c}{$a$ [\AA]} & \multicolumn{1}{c}{$4.7812(3)$} & \multicolumn{1}{c}{$4.8049(7)$}\\
      \multicolumn{1}{c}{$c$ [\AA]} & \multicolumn{1}{c}{$6.5672(4)$} & \multicolumn{1}{c}{$6.6034(5)$}\\
      \multicolumn{3}{c}{}\\
      \multicolumn{1}{c}{$R1$} & \multicolumn{1}{c}{$0.0241$} & \multicolumn{1}{c}{0.0255}\\
      \multicolumn{1}{c}{$wR2$} & \multicolumn{1}{c}{$0.0605$} & \multicolumn{1}{c}{0.0663}\\
      \multicolumn{1}{c}{Flack parameter $x$} & \multicolumn{1}{c}{$-0.076(175)$} & \multicolumn{1}{c}{$-0.109(89)$}\\
      \multicolumn{1}{c}{Absolute structure} & \multicolumn{1}{c}{Left-handed} & \multicolumn{1}{c}{Right-handed}\\
      \hline
      \end{tabular}

  \end {center}
\end{table}

Next, we performed the LDFZ crystal growth using the composition-gradient feed rod in order to obtain the single crystals of $T$Si$_{2}$ with the selected chirality, as shown in Fig.~\ref{f-FZ}(b).
The feed rod (\#3) was prepared for obtaining the left-handed NbSi$_{2}$ crystal.
It consists of TaSi$_{2}$ and NbSi$_{2}$ and thus the crystal growth starts with TaSi$_{2}$ in the region A, as shown in Fig.~\ref{f-FZ}(c).
The seed rod is initially crystallized into the single domain of TaSi$_{2}$.
Then, the feed rod goes into the phase of NbSi$_{2}$ at the beginning of the region B.
In the successive crystal growth, NbSi$_{2}$ is crystallized on the seed crystal of TaSi$_{2}$, as shown in Fig.~\ref{f-FZ}(d).
If the crystallographic structure of NbSi$_{2}$ in the region B is inherited from that of TaSi$_{2}$ in the region A during the crystal growth,
the grown NbSi$_{2}$ crystal should exhibit the chirality opposite to that of the crystal obtained by the conventional LDFZ method (shown in Table~\ref{tbl-NbTaSi2-conventional}).
The same technique was used to obtain the TaSi$_{2}$ crystals with the opposite handedness. In this case (\#4), the feed rod consists of NbSi$_{2}$ and TaSi$_{2}$ for the regions A and B, respectively.

Note that the melting temperature of TaSi$_{2}$ is 200 $^\circ$C higher than that of NbSi$_{2}$.
Thus, a drastic change of the melting temperature inevitably occurs at the boundary between the regions A and B.
In such a situation, it is difficult to stabilize the molten zone using conventional FZ furnaces with ellipsoid mirrors for focusing light emitted from halogen or xenon lamps. The temperature gradient of the conventional FZ furnaces is not as steep as that of the LDFZ furnaces. In this respect,
a precise control of the laser power can keep the molten zone stabilized in the LDFZ method using the well-focused laser.
This is one of the advantages of the crystal growth using the LDFZ method.

Figure~\ref{f-TSi-photo} shows an optical photograph
of the single crystal obtained by using the composition-gradient feed rod (\#3). This rod crystal is made of TaSi$_{2}$ in the initial stage of the crystal growth and followed by the region of NbSi$_{2}$.
The crystalline orientation of the single crystal was evaluated by Laue X-ray backscattering method.
Both the TaSi$_{2}$ and NbSi$_{2}$ regions have the same crystallographic orientation, indicating that the crystallographic structure of TaSi$_{2}$ succeeds to that of NbSi$_{2}$.
The $c$-axis of the crystal is a few degrees tilted from the growth direction of the crystal.

In the case of crystals growth using the composition-gradient feed rod (\#4), the obtained rod kept the composition change from NbSi$_2$ to TaSi$_2$ along the direction of crystal growth as well.
The crystal of TaSi$_{2}$ obtained from the region B includes some grains with different crystallographic orientations.
This misorientation of the crystal growth may be caused by the sudden change of the melting temperature at the boundary between the regions A and B.
Note that the melting temperature in the region A is lower than that in the region B in the \#4 crystal growth.
In this case, when traveling the feed rod at the boundary, it was rather difficult to keep the molten zone stable because it became very slim due to a sudden solidification. More steps of the concentration from NbSi$_2$ to TaSi$_2$ are required in the composition distribution of the feed rod in order to make the melting point change gradually across the boundary.

\begin{figure}[tb]
\begin{center}
\includegraphics[width=8cm]{fig4_TSi2_photo_rev2.eps}
\end{center}
\caption{An optical photograph of the single crystal, obtained by the LDFZ method with the composition-gradient feed rod that consists of TaSi$_{2}$ and NbSi$_{2}$.
}
\label{f-TSi-photo}
\end{figure}

To evaluate the crystal chirality of the grown rods (\#3 and \#4),
small pieces of the single crystals were prepared from both the TaSi$_{2}$ and NbSi$_{2}$ regions.
The deducted Flack parameter $x$ was determined in such pieces by the absolute structure analysis using single crystal X-ray diffraction.
The results are summarized in Tables~\ref{tbl-TaNbSi2-gradient} and \ref{tbl-NbTaSi2-gradient}.
Note that the refined Flack parameters of all the pieces collected from the grown rods are the same within two-sigma error bar.

\begin{table}[htbp]
      \caption{Crystal structure analysis of $T$Si$_{2}$, obtained by the LDFZ method using the composition-gradient feed rod \#3.
                    It consists of TaSi$_{2}$ in the region A and NbSi$_{2}$ in the region B, as shown in Fig.~\ref{f-FZ}(a).}
      \label{tbl-TaNbSi2-gradient}

  \begin{center}

      \begin{tabular}{rrr}
      \hline
      \multicolumn{1}{c}{} & \multicolumn{1}{c}{TaSi$_{2}$} & \multicolumn{1}{c}{NbSi$_{2}$}\\
      \hline
      \multicolumn{1}{c}{Space group} & \multicolumn{1}{c}{$P6_{4}22$} & \multicolumn{1}{c}{$P6_{4}22$}\\
      \multicolumn{1}{c}{$a$ [\AA]} & \multicolumn{1}{c}{$4.7750(3)$} & \multicolumn{1}{c}{4.7904(2)}\\
      \multicolumn{1}{c}{$c$ [\AA]} & \multicolumn{1}{c}{$6.5550(5)$} & \multicolumn{1}{c}{6.5840(4)}\\
      \multicolumn{3}{c}{}\\
      \multicolumn{1}{c}{$R1$} & \multicolumn{1}{c}{$0.0115$} & \multicolumn{1}{c}{0.0074}\\
      \multicolumn{1}{c}{$wR2$} & \multicolumn{1}{c}{$0.0249$} & \multicolumn{1}{c}{0.0164}\\
      \multicolumn{1}{c}{Flack parameter $x$} & \multicolumn{1}{c}{$0.057(136)$} & \multicolumn{1}{c}{$0.021(148)$}\\
      \multicolumn{1}{c}{Absolute structure} & \multicolumn{1}{c}{Left-handed} & \multicolumn{1}{c}{Left-handed}\\
      \hline
      \end{tabular}

  \end {center}
\end{table}

\begin{table}[htbp]
      \caption{Crystal structure analysis of $T$Si$_{2}$, obtained by the LDFZ method with the composition-gradient feed rod \#4.
                    It consists of NbSi$_{2}$ in the region A and TaSi$_{2}$ in the region B as shown in Fig.~\ref{f-FZ}(a).}
      \label{tbl-NbTaSi2-gradient}

  \begin{center}

      \begin{tabular}{rrr}
      \hline
      \multicolumn{1}{c}{} & \multicolumn{1}{c}{NbSi$_{2}$} & \multicolumn{1}{c}{TaSi$_{2}$}\\
      \hline
      \multicolumn{1}{c}{Space group} & \multicolumn{1}{c}{$P6_{2}22$} & \multicolumn{1}{c}{$P6_{2}22$}\\
      \multicolumn{1}{c}{$a$ [\AA]} & \multicolumn{1}{c}{$4.7975(2)$} & \multicolumn{1}{c}{4.7770(4)}\\
      \multicolumn{1}{c}{$c$ [\AA]} & \multicolumn{1}{c}{$6.5837(4)$} & \multicolumn{1}{c}{6.5604(5)}\\
      \multicolumn{3}{c}{}\\
      \multicolumn{1}{c}{$R1$} & \multicolumn{1}{c}{$0.0164$} & \multicolumn{1}{c}{0.0287}\\
      \multicolumn{1}{c}{$wR2$} & \multicolumn{1}{c}{$0.0371$} & \multicolumn{1}{c}{0.0745}\\
      \multicolumn{1}{c}{Flack parameter $x$} & \multicolumn{1}{c}{$-0.077(145)$} & \multicolumn{1}{c}{$-0.021(153)$}\\
      \multicolumn{1}{c}{Absolute structure} & \multicolumn{1}{c}{Right-handed} & \multicolumn{1}{c}{Right-handed}\\
      \hline
      \end{tabular}

  \end {center}
\end{table}

For the crystal growth using the feed rod (\#3), the Flack parameter $x$ of TaSi$_{2}$ in the region A is almost zero with regard to the structure based on the space group of $P6_{4}22$, as shown in Table~\ref{tbl-TaNbSi2-gradient}.
Thus, the obtained crystal of TaSi$_{2}$ has the left-handed crystal structure.
Importantly, the crystal handedness of NbSi$_{2}$ turns to be the same as that of TaSi$_{2}$, while the lattice constants in the NbSi$_{2}$ region reproduces the value obtained by the conventional FZ method.

The handedness control was also performed by using the feed rod (\#4) shown in Table~\ref{tbl-NbTaSi2-gradient}.
For the grown rod, the crystal handedness of NbSi$_{2}$ in the region A and that of TaSi$_{2}$ in the region B were determined to be right-handed.
As for the uniformity of chirality of TaSi$_{2}$ in the region B, it turned out that all the grains examined, including the ones with different crystallographic orientations from that of NbSi$_{2}$ in the region A, exhibit the right-handed crystal structure. In this respect, the crystal handedness of the region B is likely to inherit from that of the region A even though some misorientation occurs during the crystal growth.

These experiments demonstrate that the single crystals of left-handed NbSi$_2$ and right-handed TaSi$_2$ were successfully obtained. These crystals have the handedness opposite to that of the crystals grown by the conventional methods.
Once the single crystals with the opposite chirality are obtained, such crystals can be used as a seed crystal for another crystal growth.
For instance, by feeding polycrystalline NbSi$_{2}$ on the chirality-controlled seed rod, the NbSi$_{2}$ crystals with the opposite chirality are successively obtained without using the composition-gradient feed rod.
This method can be applicable for controlling chirality of $T$Si$_{2}$ with different $T$ element.
We emphasize that a preparation of the chirality-controlled seed crystals is essential in the chirality-controlled crystal growth.

\subsection{CISS measurements}

The CISS measurements were performed with the chiral crystal rods, grown by the LDFZ method, as illustrated in Fig.~\ref{f-CISS-draw}(a).
The spin polarization is monitored by a transverse voltage generated in a tungsten (W) electrode, deposited on a half-cylindrical region of the crystal.

For the CISS measurements, we prepared the chirality-controlled NbSi$_2$ crystals of a cylinder shape with a  dimension of 5 mm in diameter and 70 mm in length, as shown in Fig.~\ref{f-CISS-draw}(b).
More accurately, the LDFZ grown crystals with the right- and left-handedness were cut into a slice with 70 mm long.
The W electrode with 1.5 mm in width and 7 nm in thickness was deposited on the specimen by using a sputtering system (M09-0014, SEED Lab. Corporation).
The electrodes for the current injection are set to be 60 mm apart from each other and contain the W electrode between them.

\begin{figure}[tb]
\begin{center}
\includegraphics[width=8cm]{fig5_CISS-drawing_rev3.eps}
\end{center}
\caption{(a) Schematics of the CISS measurements for a cylinder-shaped specimen. (b) A photograph of the device made for the CISS measurements of the NbSi$_{2}$ crystal with the tungsten (W) electrode. The $I$-$V$ characteristics were collected by applying the charge current into the crystal.
}
\label{f-CISS-draw}
\end{figure}

\begin{figure}[tb]
\begin{center}
\includegraphics[width=8cm]{fig6_CISS_rev3.eps}
\end{center}
\caption{
The $I$-$V$ characteristics of NbSi$_{2}$ collected at room temperature without magnetic field.
The CISS voltage $V_{xy}$ was measured as a function of $I_{xx}$ in the right-handed ($P$6$_{2}$22) and left-handed ($P$6$_{4}$22) NbSi$_{2}$ crystals with the tungsten electrodes.
}
\label{f-CISS}
\end{figure}

The longitudinal electrical current $I_{xx}$ was applied along the long axis of the crystal.
As evaluated by Laue X-ray, the longitudinal direction of the crystals is deviated from the $c$-axis by several degrees and thus the $I_{xx}$ flows in a direction nearly parallel to the $c$-axis.

In a separate experiment, the electrical resistivity $V_{xx}/I_{xx}$ was measured. It was determined to be $39\mu\Omega\cdot$cm at room temperature, which was slightly smaller than those of single crystalline NbSi$_{2}$ obtained by the Czochralski method \cite{Onuki2014}.
This value was an order of magnitude smaller than that of polycrystalline NbSi$_{2}$ \cite{Shishido2021}.

The CISS signals were collected with the NbSi$_{2}$ crystals at room temperature without magnetic field.
The transverse voltage $V_{xy}$ was generated linearly as a function of $I_{xx}$ in the NbSi$_{2}$ crystals, as shown in Fig.~\ref{f-CISS}.

Importantly, the sign of the linear $I$-$V$ characteristics depends on the handedness of the crystals.
The sign reversal occurs because of the spin polarization flipping accompanied by a change of the crystal chirality.
These results indicate that the CISS measurements can probe the handedness of the crystal over centimeters and work as a macroscopic method for chirality determination.

Note that the $V_{xy}$ signals mainly result from the CISS effect rather than longitudinal electrical resistance of NbSi$_{2}$ caused by a misalignment of the W electrodes.
If the signals completely come from the misalignment, they correspond to a resistance of the specimen of 0.3 mm in length.
Taking the width of the W electrode (1.5 mm) into account, it is hard to see such a large misalignment of the electrode deposited on the rod.

Note that the magnitude of the slope $V_{xy} / I_{xx}$ in the present study is three orders magnitude smaller than that obtained in the micrometer-sized single crystals \cite{Shiota2021},
while it is an order of magnitude smaller than that for polycrystals \cite{Shishido2021}.
However, by converting $V_{xy} / I_{xx}$ into the resistivity $\rho_{xy}$ using the specimen dimensions and then normalized by $\rho_{xx}$, the resultant $\rho_{xy} / \rho_{xx}$ exhibits the same order of the magnitude as those for the micrometer-sized single crystals and for the polycrystals.

The spatial-resolved CISS measurements revealed the presence of chirality domains with the opposite chirality in the micrometer-sized specimen \cite{Shiota2021}.
The bulk crystals used in the present study may also contain such chirality domains.
The Flack parameter determined by the single crystal X-ray diffraction measurements indicates the presence of such chirality domains up to a few percent of all the region within its error bar.
However, there is no direct information on a microscopic distribution of the chirality domains in the bulk crystal.
A further quantitative discussion requires for identifying a locational distribution of the chirality domains using the spatial-resolved CISS measurements in a resolution of micrometers or smaller scale.
Such precise CISS measurements for the centimeter-sized crystals will be a future issue.

\section{Summary and Perspective}
In summary, we grew the chirality-selected single crystals of NbSi$_{2}$ and TaSi$_{2}$ using the LDFZ method.
In the case of the conventional LDFZ growth using the seed and feed rods that consist of the same compound, the right-handed NbSi$_{2}$ the left-handed TaSi$_{2}$ crystals were obtained as reported in the literature.
However, in the case of the LDFZ crystal growth using the composition-gradient feed rods,
the chirality of the obtained NbSi$_{2}$ crystals was found to be the left-handed,
while that of the TaSi$_{2}$ crystals was the right-handed.
The chirality is opposite to the reported one in both crystals
and is hardly obtained by the conventional crystal growth technique.
When such crystals with the opposite chirality are used for the seed crystals,
the chiral crystal of the same quality in terms of chirality can be obtained even by using the conventional LDFZ method.
These results indicate that the chirality of the seed crystal is transferred to the grown crystal.

The spin polarization phenomena were observed in the rod-shaped single crystals of NbSi$_{2}$ with millimeters in diameter and centimeters in length.
The sign of the slope in the $I$-$V$ characteristics for the NbSi$_{2}$ bulk rod crystals was qualitatively consistent with the results obtained with the micrometer-sized specimens of NbSi$_{2}$.
The CISS effect appears in such large crystals, and will open up new possibilities for a non-destructive evaluation of chiral crystal structures in inorganic materials.

The next interesting issue will be an investigation of the locational distribution of the crystal chirality over the long crystal.
The distribution of the chirality can also be detected by a circular-polarized resonant X-ray micro diffraction method \cite{Ohsumi2013}.
However, this method proves only on the surface of the specimen because of very large X-ray absorption.
On the other hand, the CISS signal reflects the crystal chirality in the specimen where the current is injected \cite{Inui2020,Shiota2021,Shishido2021}.
When the current is injected in a very narrow region of the specimen, the chirality distribution can be evaluated with a high spatial resolution.

\ack

This work was supported by JSPS KAKENHI Grant Numbers 15H05885, 19KK0070 and 20H02642.

\section*{References}


\section{Introduction}
For a fixed $k\ge 3$, it has long been known that determining the $k$-colourability of a general graph is an NP-complete problem~\cite{Karp1972}. When restricted to families of graphs, however, polynomial-time algorithms can be developed. For example, $k$-colourability can be decided in polynomial-time for all $k$ for perfect graphs~\cite{Grotschel1984} and $P_5$-free graphs~\cite{Hoang2010}. The latter is of particular interest because deciding $k$-colourability of $H$-free graphs remains NP-complete if $H$ contains an induced cycle~\cite{KaminskiLozin2007,Maffray1996} or claw~\cite{Holyer1981,LevenGail1983}. Further, the problem remains NP-complete for $P_6$-free graphs for all $k\ge 5$ and for $P_7$-free graphs for all $k\ge 4$~\cite{Huang2016}. Completing the complexity dichotomy for $4$-colourability of $H$-free graphs when $H$ is connected, is the polynomial-time algorithm to decide $4$-colourability of $P_6$-free graphs~\cite{P6freeconf,P6free1,P6free2}. Thus, $P_5$ is the largest connected graph whose forbidding results in the existence (assuming P$\neq$NP) of polynomial-time $k$-colourability algorithms for all $k$. While many open problems remain on the complexity of $k$-colouring $H$-free graphs, including $3$-colouring $P_8$-free graph~\cite{Bonomo2018} and $6$-colouring $(P_3+P_2)$-free graphs~\cite{HajebiLiSpirkl2021, Golovach2017}, the focus of this paper is to further the study of polynomial-time \textit{certifying} $k$-colourability algorithms. An algorithm is certifying if it returns an easily verifiable witness with each output to certify its correctness. A certifying algorithm offers the considerable advantage of an extra layer of certainty for users that the algorithm was correctly implemented. This has led to an increase in the development in certifying algorithms in the last decade and even to calls from some that only certifying algorithms should be developed for solving complex tasks~\cite{McConnell2011}. For a $k$-colourability algorithm to be certifying, it must return a $k$-colouring if one exists and a $(k+1)$-vertex-critical induced subgraph if the graph is not $k$-colourable. In a family of graphs, if there are only finitely many $(k+1)$-vertex critical graphs, then a polynomial-time $k$-colourability algorithm that certifies negative answers can be developed by determining if any of the vertex-critical graphs are induced subgraphs of the input graph (see~\cite{P5banner2019} for more details).

The classification of all 12 $4$-vertex-critical $P_5$-free graphs was used to develop a linear-time certifying $3$-colourability algorithm for $P_5$-free graphs~\cite{Bruce2009,MaffrayMorel2012}. For all $k\ge 4$, however, there are infinitely $(k+1)$-vertex-critical $P_5$-free graphs~\cite{Hoang2015}. This has led to tremendous interest in determining for which subfamilies of $P_5$-free graphs the polynomial-time $k$-colourability algorithms can be extended to certifying ones. Perhaps the most comprehensive known result is the dichotomy for $H$ of order $4$ that in all cases except $H=2K_2$ or $H=K_3+P_1$, there are only a finite number of $k$-vertex critical $(P_5,H)$-free graphs for all $k$~\cite{KCameron2021}. The open problem of completing such a dichotomy for $H$ of order $5$ was also posed in~\cite{KCameron2021}. Toward resolving this open problem, it is known that there are only finitely many $k$-vertex critical $(P_5,H$)-free for $k\le 5$ when $H=C_5$~\cite{Hoang2015} and for all $k$ when:

\begin{itemize}
\item $H$=banner~\cite[Theorem 3(i)]{Brause2022}
\item $H=K_{2,3}$ and $H=K_{1,4}$~\cite{Kaminski2019}
\item $H=P_2+3P_1$~\cite{CameronHoangSawada2022}
\item $H=P_3+2P_1$~\cite{AbuadasCameronHoangSawada2022}
\item $H=\overline{P_5}$~\cite{Dhaliwal2017}
\item $H=\overline{P_3+P_2}$ and $H=$gem~\cite{CaiGoedgebeurHuang2021}
\end{itemize}

Of particular interest among these results are the last two bullet points, as complete lists of all such $k$-vertex-critical graphs were given for $k=4$ and $k=5$~\cite{Dhaliwal2017,CaiGoedgebeurHuang2021}. These complete lists are essential for implementing the corresponding certifying algorithms. In this paper, we present a new proof that there are only finitely $k$-vertex-critical $(P_5,\text{ gem})$-free graphs. The proof in~\cite{CaiGoedgebeurHuang2021} involves bounding the order of every $k$-vertex-critical graph in one difficult structural case (this case is that the graph belongs to ${\cal H}$, which will be defined later). Our new proof shows the stronger result that there are in fact no vertex-critical graphs in this difficult structural case. Our stronger result allows us to use a simple computer search to completely characterize all such graphs for $k=6$ and $k=7$, allowing for the implementation of polynomial-time certifying $5$- and $6$-colourability algorithms. Beyond their use for implementing certifying colourability algorithms, vertex-critical graphs are of interest in many other areas. In particular, vertex-critical $(P_5,\text{ gem})$-free graphs were very recently used to prove the Borodin-Kostochka Conjecture for $(P_5,\text{ gem})$-free graphs~\cite{CranstonLafayetteRabern2022}. Our proof that there are only finitely many $k$-vertex-critical $(P_5,\text{ gem})$-free graphs for all $k$ would allow for a simpler proof of the Main Theorem in~\cite{CranstonLafayetteRabern2022}, in particular, Case 11, the most involved case, is no longer required.

The rest of the paper is organized as follows: we first introduce some notation and definitions before proving our result on $k$-vertex-critical $(P_5,\text{ gem})$-free graphs in Section~\ref{sec:critP5gem}. In Section~\ref{sec:6and7} we give complete lists of all 19 graphs for $k=6$ and all 46 graphs for $k=7$ and outline our exhaustive computer search (the graphs in graph6 format and our code is available at \cite{graphfiles}). We then conclude by posing some open problems for future research in Section~\ref{sec:conclusion}.

\subsection{Definitions and Notation}

Let $\chi(G)$ denote the \textit{chromatic number} of $G$. A graph $G$ is \textit{$k$-vertex-critical} if $\chi(G)=k$ and $\chi(G-v)<k$ for all $v\in V(G)$. For a given colouring $c$ of a graph $G$ and subset $S\subseteq V(G)$, let $c(S)$ denote the set of all colours used on vertices in $S$. For $S\subseteq V(G)$, let $G[S]$ denote the subgraph of $G$ induced by $S$. For vertices $u,v\in V(G)$, we write $u\sim v$ if $u$ and $v$ are adjacent and $u\nsim v$ is $u$ and $v$ are non-adjacent. A for $S,T\subseteq V(G)$, we say $S$ is \textit{complete} (\textit{anti-complete}) to $T$ if $s\sim t$ ($s\nsim t$) for all $s\in S$ and $t\in T$. A \textit{module} of a graph $G$ is a set $S$ such that for all $v\in V(G)\setminus S$, $v$ is either complete or anti-complete to $S$. For $S\subseteq V(G)$, we say $S$ is a \textit{stable set} if $u\nsim v$ for all $u,v\in S$. A stable set is called \textit{good} if it meets every maximum clique of $G$. Given a graph $G$, a clique ($P_4$-free) expansion of $G$ is a graph obtained from $G$ by replacing some or all of its vertices by cliques ($P_4$-free graphs) and joining every vertex in these cliques ($P_4$-free graphs) according to the adjacencies of $G$.

\section{Vertex-critical ($P_5$, gem)-free graphs}\label{sec:critP5gem}
Throughout this section, we will make extensive use of a structural characterization give in~\cite{Chudetal2020} which requires some further definitions. Let $G_1,G_2,\ldots,G_{10}$ be defined as shown in Figure~\ref{fig:G1G7} and let $\mathcal{G}_i$ and $\mathcal{G}^{\ast}_i$ denote the set of all $P_4$-free and clique expansions of $G_i$, respectively. Throughout, for a graph $G\in \mathcal{G}_i$, we will use $ {Q_j}$ to denote the set of vertices used to replace $v_j$ in the graph $G_i$ according to the labellings in Figure~\ref{fig:G1G7}.
Let ${\cal H}$ be the class of ($P_5$, gem)-free graphs $G$ such that $V(G)$ can be partitioned into seven non-empty sets, $A_1,\ldots,A_7$ such that:
\begin{itemize}
\item $A_1$ is complete to $A_2\cup A_5$ and anti-complete to $A_3\cup A_4\cup A_6\cup A_7$,
\item $A_2$ is complete to $A_1\cup A_3$ and anti-complete to $A_4\cup A_5\cup A_6\cup A_7$,
\item $A_3$ is complete to $A_2\cup A_4\cup A_6$ and anti-complete to $A_1\cup A_5\cup A_7$,
\item $A_4$ is complete to $A_3\cup A_5\cup A_6$ and anti-complete to $A_1\cup A_2\cup A_7$,
\item $A_5$ is complete to $A_1\cup A_4$ and anti-complete to $A_2\cup A_3\cup A_6\cup A_7$,
\item every component of $G[A_7]$ is a module, and
\item every vertex in $A_7$ has a neighbour in $A_6$.
\end{itemize}

\noindent We refer the reader to Figure~\ref{fig:H} for the general form of a graph in ${\cal H}$. Let ${\cal H}^{*}\subseteq {\cal H}$ be defined by all graphs in ${\cal H}$ such that $A_i$ is a clique for all $i\in \{1,2,3,4,5\}$ and each component of $A_7$ is a clique. We are now ready for the structural characterization.
\setcounter{subfigure}{0}
\begin{figure}[htb]
\def0.5{0.5}
\def1{1}
\centering
\subfigure[$G_1$]{
\scalebox{0.5}{
\begin{tikzpicture}
\GraphInit[vstyle=Classic]
\Vertex[Lpos=180,L=\hbox{$v_5$},x=0.0cm,y=3.1235cm]{v0}
\Vertex[L=\hbox{$v_2$},x=5.0cm,y=3.0618cm]{v1}
\Vertex[Lpos=180,L=\hbox{$v_4$},x=0.942cm,y=0.0422cm]{v2}
\Vertex[Lpos=90,L=\hbox{$v_1$},x=2.5253cm,y=5.0cm]{v3}
\Vertex[L=\hbox{$v_3$},x=4.0006cm,y=0.0cm]{v4}

\Edge[](v0)(v2)
\Edge[](v0)(v3)
\Edge[](v1)(v3)
\Edge[](v1)(v4)
\Edge[](v2)(v4)
\end{tikzpicture}}}
\qquad
\subfigure[$G_2$]{
\scalebox{0.5}{
\begin{tikzpicture}
\GraphInit[vstyle=Classic]
\Vertex[Lpos=180,L=\hbox{$v_5$},x=0.0cm,y=3.1235cm]{v4}
\Vertex[L=\hbox{$v_2$},x=5.0cm,y=3.0618cm]{v1}
\Vertex[Lpos=180,L=\hbox{$v_4$},x=0.942cm,y=0.0422cm]{v3}
\Vertex[Lpos=90,L=\hbox{$v_1$},x=2.5253cm,y=5.0cm]{v0}
\Vertex[L=\hbox{$v_3$},x=4.0006cm,y=0.0cm]{v2}
\Vertex[L=\hbox{$v_6$},x=2.5253cm,y=2cm]{v5}
\Edge[](v0)(v1)
\Edge[](v0)(v4)
\Edge[](v0)(v5)
\Edge[](v1)(v2)
\Edge[](v2)(v3)
\Edge[](v2)(v5)
\Edge[](v3)(v4)
\Edge[](v3)(v5)
\end{tikzpicture}}}
\qquad
\subfigure[$G_3$]{
\scalebox{0.5}{
\begin{tikzpicture}
\GraphInit[vstyle=Classic]
\Vertex[Lpos=180,L=\hbox{$v_5$},x=0.0cm,y=3.1235cm]{v4}
\Vertex[L=\hbox{$v_2$},x=5.0cm,y=3.0618cm]{v1}
\Vertex[Lpos=180,L=\hbox{$v_4$},x=0.942cm,y=0.0422cm]{v3}
\Vertex[Lpos=90,L=\hbox{$v_1$},x=2.5253cm,y=5.0cm]{v0}
\Vertex[L=\hbox{$v_3$},x=3cm,y=1.2cm]{v2}
\Vertex[L=\hbox{$v_6$},x=2cm,y=3.1235cm]{v5}
\Vertex[L=\hbox{$v_7$},x=4.0006cm,y=0.0cm]{v6}
\Edge[](v0)(v1)
\Edge[](v0)(v4)
\Edge[](v0)(v5)
\Edge[](v1)(v2)
\Edge[](v1)(v6)
\Edge[](v2)(v3)
\Edge[](v2)(v5)
\Edge[](v3)(v4)
\Edge[](v3)(v5)
\Edge[](v3)(v6)
\end{tikzpicture}}}
\qquad
\subfigure[$G_4$]{
\scalebox{0.5}{
\begin{tikzpicture}
\GraphInit[vstyle=Classic]
\Vertex[Lpos=180,L=\hbox{$v_5$},x=0.0cm,y=3.1235cm]{v4}
\Vertex[L=\hbox{$v_2$},x=5.0cm,y=3.0618cm]{v1}
\Vertex[Lpos=180,L=\hbox{$v_4$},x=0.942cm,y=0.0422cm]{v3}
\Vertex[Lpos=90,L=\hbox{$v_1$},x=2.5253cm,y=5.0cm]{v0}
\Vertex[L=\hbox{$v_3$},x=4.0006cm,y=0.0cm]{v2}
\Vertex[L=\hbox{$v_6$},x=3.5cm,y=1.2cm]{v5}
\Vertex[Lpos=180,L=\hbox{$v_7$},x=2cm,y=3.0618cm]{v6}
\Edge[](v0)(v1)
\Edge[](v0)(v4)
\Edge[](v0)(v5)
\Edge[](v0)(v6)
\Edge[](v1)(v2)
\Edge[](v1)(v6)
\Edge[](v2)(v3)
\Edge[](v2)(v5)
\Edge[](v3)(v4)
\Edge[](v3)(v5)
\Edge[](v3)(v6)
\end{tikzpicture}}}
\qquad
\subfigure[$G_5$]{
\scalebox{0.5}{
\begin{tikzpicture}
\GraphInit[vstyle=Classic]
\Vertex[Lpos=180,L=\hbox{$v_5$},x=0.0cm,y=3.1235cm]{v4}
\Vertex[Lpos=-135,L=\hbox{$v_2$},x=4cm,y=3.0618cm]{v1}
\Vertex[Lpos=180,L=\hbox{$v_4$},x=0.942cm,y=0.0422cm]{v3}
\Vertex[Lpos=90,L=\hbox{$v_1$},x=2.5253cm,y=5.0cm]{v0}
\Vertex[L=\hbox{$v_3$},x=4.0006cm,y=0.0cm]{v2}
\Vertex[L=\hbox{$v_8$},x=5.0cm,y=3.0618cm]{v7}
\Vertex[L=\hbox{$v_6$},x=2.5cm,y=1.6729cm]{v5}
\Vertex[Lpos=-45,L=\hbox{$v_7$},x=1cm,y=3.0618cm]{v6}
\Edge[](v0)(v1)
\Edge[](v0)(v4)
\Edge[](v0)(v5)
\Edge[](v0)(v6)
\Edge[](v0)(v7)
\Edge[](v1)(v2)
\Edge[](v1)(v6)
\Edge[](v2)(v3)
\Edge[](v2)(v5)
\Edge[](v2)(v7)
\Edge[](v3)(v4)
\Edge[](v3)(v5)
\Edge[](v3)(v6)
\end{tikzpicture}}}
\qquad
\subfigure[$G_6$]{
\scalebox{0.5}{
\begin{tikzpicture}
\GraphInit[vstyle=Classic]
\Vertex[Lpos=180,L=\hbox{$v_5$},x=0.0cm,y=3.1235cm]{v4}
\Vertex[L=\hbox{$v_2$},x=5.0cm,y=3.0618cm]{v1}
\Vertex[Lpos=180,L=\hbox{$v_4$},x=0.942cm,y=0.0422cm]{v3}
\Vertex[Lpos=90,L=\hbox{$v_1$},x=3cm,y=3.5cm]{v0}
\Vertex[L=\hbox{$v_3$},x=4.0006cm,y=0.0cm]{v2}
\Vertex[Lpos=90,L=\hbox{$v_8$},x=2.5253cm,y=5.0cm]{v7}
\Vertex[L=\hbox{$v_6$},x=1cm,y=2.5cm]{v5}
\Vertex[Lpos=-90,L=\hbox{$v_7$},x=3.9687cm,y=2.5cm]{v6}
\Edge[](v0)(v1)
\Edge[](v0)(v4)
\Edge[](v0)(v5)
\Edge[](v0)(v6)
\Edge[](v1)(v2)
\Edge[](v1)(v6)
\Edge[](v1)(v7)
\Edge[](v2)(v3)
\Edge[](v2)(v5)
\Edge[](v3)(v4)
\Edge[](v3)(v5)
\Edge[](v3)(v6)
\Edge[](v4)(v7)
\Edge[](v5)(v7)
\end{tikzpicture}}}
\qquad
\subfigure[$G_7$]{
\scalebox{0.5}{
\begin{tikzpicture}
\GraphInit[vstyle=Classic]
\Vertex[Lpos=45,L=\hbox{$v_1$},x=3.5cm,y=1cm]{v0}
\Vertex[Lpos=180,L=\hbox{$v_2$},x=0.942cm,y=0.0422cm]{v1}
\Vertex[L=\hbox{$v_3$},x=4.0006cm,y=0.0cm]{v2}
\Vertex[Lpos=90,L=\hbox{$v_4$},x=2.5253cm,y=5.0cm]{v3}
\Vertex[Lpos=135,L=\hbox{$v_5$},x=1.4765cm,y=2.4803cm]{v4}
\Vertex[L=\hbox{$v_6$},x=5.0cm,y=3.0618cm]{v5}
\Vertex[Lpos=45,L=\hbox{$v_7$},x=1.8cm,y=1cm]{v6}
\Vertex[Lpos=180,L=\hbox{$v_8$},x=0.0cm,y=3.1235cm]{v7}
\Edge[](v0)(v1)
\Edge[](v0)(v2)
\Edge[](v0)(v3)
\Edge[](v0)(v6)
\Edge[](v1)(v2)
\Edge[](v1)(v4)
\Edge[](v1)(v7)
\Edge[](v2)(v5)
\Edge[](v3)(v4)
\Edge[](v3)(v5)
\Edge[](v3)(v7)
\Edge[](v4)(v5)
\Edge[](v4)(v6)
\end{tikzpicture}}}
\subfigure[$G_8$]{
\scalebox{0.5}{
\begin{tikzpicture}
\GraphInit[vstyle=Classic]
\Vertex[Lpos=180,L=\hbox{$v_5$},x=0.0cm,y=3.1235cm]{v4}
\Vertex[L=\hbox{$v_2$},x=5.0cm,y=3.0618cm]{v1}
\Vertex[Lpos=180,L=\hbox{$v_4$},x=0.942cm,y=0.0422cm]{v3}
\Vertex[Lpos=90,L=\hbox{$v_1$},x=2.5253cm,y=5.0cm]{v0}
\Vertex[L=\hbox{$v_3$},x=4.0006cm,y=0.0cm]{v2}

\Vertex[Lpos=135,L=\hbox{$v_6$},x=0.5cm,y=4.3cm]{v5}
\Vertex[Lpos=-90,L=\hbox{$v_7$},x=1.5cm,y=1.4398cm]{v6}
\Vertex[Lpos=90,L=\hbox{$v_8$},x=3.7265cm,y=1.4398cm]{v7}
\Edge[](v0)(v1)
\Edge[](v0)(v4)
\Edge[](v0)(v5)
\Edge[](v0)(v6)
\Edge[](v1)(v2)
\Edge[](v1)(v7)
\Edge[](v2)(v3)
\Edge[](v2)(v6)
\Edge[](v3)(v4)
\Edge[](v3)(v5)
\Edge[](v3)(v7)
\Edge[](v4)(v6)
\Edge[](v5)(v7)
\Edge[](v6)(v7)
\end{tikzpicture}}}
\qquad
\subfigure[$G_9$]{
\scalebox{0.5}{
\begin{tikzpicture}
\GraphInit[vstyle=Classic]
\Vertex[Lpos=180,L=\hbox{$v_5$},x=0.0cm,y=3.1235cm]{v4}
\Vertex[L=\hbox{$v_2$},x=5.0cm,y=3.0618cm]{v1}
\Vertex[Lpos=180,L=\hbox{$v_4$},x=0.942cm,y=0.0422cm]{v3}
\Vertex[Lpos=90,L=\hbox{$v_1$},x=2.5253cm,y=5.0cm]{v0}
\Vertex[L=\hbox{$v_3$},x=4.0006cm,y=0.0cm]{v2}

\Vertex[Lpos=135,L=\hbox{$v_6$},x=0.5cm,y=4.3cm]{v5}
\Vertex[Lpos=45,L=\hbox{$v_7$},x=1.5cm,y=1.4398cm]{v6}
\Vertex[Lpos=90,L=\hbox{$v_8$},x=3.7265cm,y=1.4398cm]{v7}
\Vertex[L=\hbox{$v_9$},x=3.6cm,y=3.0618cm]{v8}
\Edge[](v0)(v1)
\Edge[](v0)(v4)
\Edge[](v0)(v5)
\Edge[](v0)(v6)
\Edge[](v0)(v8)
\Edge[](v1)(v2)
\Edge[](v1)(v7)
\Edge[](v2)(v3)
\Edge[](v2)(v6)
\Edge[](v3)(v4)
\Edge[](v3)(v5)
\Edge[](v3)(v7)
\Edge[](v3)(v8)
\Edge[](v4)(v6)
\Edge[](v5)(v7)
\Edge[](v6)(v7)
\end{tikzpicture}}}
\qquad
\subfigure[$G_{10}$]{
\scalebox{0.5}{
\begin{tikzpicture}
\GraphInit[vstyle=Classic]
\Vertex[Lpos=180,L=\hbox{$v_5$},x=0.0cm,y=3.1235cm]{v4}
\Vertex[L=\hbox{$v_2$},x=5.0cm,y=3.0618cm]{v1}
\Vertex[Lpos=180,L=\hbox{$v_4$},x=0.942cm,y=0.0422cm]{v3}
\Vertex[Lpos=90,L=\hbox{$v_1$},x=2.5253cm,y=5.0cm]{v0}
\Vertex[L=\hbox{$v_3$},x=4.0006cm,y=0.0cm]{v2}

\Vertex[Lpos=135,L=\hbox{$v_6$},x=0.5cm,y=4.3cm]{v5}
\Vertex[Lpos=-90,L=\hbox{$v_7$},x=1.5cm,y=1.4398cm]{v6}
\Vertex[Lpos=90,L=\hbox{$v_8$},x=4cm,y=1.4398cm]{v7}
\Vertex[Lpos=45,L=\hbox{$v_9$},x=2.6cm,y=3.5cm]{v8}
\Edge[](v0)(v1)
\Edge[](v0)(v4)
\Edge[](v0)(v5)
\Edge[](v0)(v6)
\Edge[](v1)(v2)
\Edge[](v1)(v7)
\Edge[](v1)(v8)
\Edge[](v2)(v3)
\Edge[](v2)(v6)
\Edge[](v2)(v8)
\Edge[](v3)(v4)
\Edge[](v3)(v5)
\Edge[](v3)(v7)
\Edge[](v4)(v6)
\Edge[](v4)(v8)
\Edge[](v5)(v7)
\Edge[](v5)(v8)
\Edge[](v6)(v7)
\end{tikzpicture}}}
\caption{Special $(P_5,\text{ gem})$-free graphs used in the structural characterization in \cite{Chudetal2020}.}\label{fig:G1G7}
\end{figure}

\begin{figure}[htb]
\centering
\begin{tikzpicture}
\GraphInit[vstyle=Normal]
\Vertex[L=\hbox{$A_5$},x=0.0cm,y=3.1235cm]{A1}
\Vertex[L=\hbox{$A_3$},x=5.0cm,y=3.0618cm]{A3}
\Vertex[L=\hbox{$A_1$},x=0.942cm,y=0.0422cm]{A5}
\Vertex[L=\hbox{$A_4$},x=2.5253cm,y=5.0cm]{A2}
\Vertex[L=\hbox{$A_2$},x=4.0006cm,y=0.0cm]{A4}
\Vertex[L=\hbox{$A_6$},x=2.4923cm,y=2.2487cm]{A6}
\Vertex[L=\hbox{$A_7$},x=3cm,y=.9cm]{A7}
\SetUpEdge[lw=3pt]

\Edge[](A1)(A2)
\Edge[](A1)(A5)
\Edge[](A2)(A3)
\Edge[](A2)(A6)
\Edge[](A3)(A4)
\Edge[](A3)(A6)
\Edge[](A4)(A5)
\Edge[](A6)(A5)
\SetUpEdge[lw=1pt]
\Edge[style=dashed](A6)(A7)
\end{tikzpicture}
\caption{General form of a graph in ${\cal H}$ where the thick black lines denote sets which are complete to each other and the thin dashed line denotes any adjacency satisfying the requirements in~\cite{Chudetal2020}.}\label{fig:H}
\end{figure}

\begin{theorem}[\cite{Chudetal2020}]\label{thm:P5gemfreestructure}
If $G$ is a $(P_5,\text{ gem})$-free graph, then $G$ is either perfect, in $\mathcal{G}_i$ for some $i\in\{1,\ldots, 10\}$, or $G\in{\cal H}$.
\end{theorem}

Since we are interested in vertex-critical graphs, the following result will be helpful.
\begin{lemma}[\cite{Dhaliwal2017}]\label{lem:module}
If $G$ is $k$-vertex-critical with a non-trivial module $M$, then $M$ is $m$-vertex-critical for some $m<k$.
\end{lemma}

\begin{lemma}\label{lem:kcritinGastorHast}
If $G$ is a $k$-vertex-critical $(P_5,\text{ gem})$-free graph, then $G$ is $K_k$, $G\in {\cal H}^{*}$, or $G\in \mathcal{G}_i^{\ast}$ for some $i=1,\ldots 10$.
\end{lemma}
\begin{proof}
Let $G$ be a $k$-vertex-critical $(P_5,\text{ gem})$-free graph. If $G$ is perfect, then $G=K_k$, otherwise $G\in\mathcal{G}_i$ for some $i\in\{1,\ldots, 10\}$, or $G\in{\cal H}$ from Theorem~\ref{thm:P5gemfreestructure}. Suppose $G\in\mathcal{G}_i$ for some $i\in\{1,\ldots, 10\}$. By definition, every $Q_j$ of $\mathcal{G}_i$ is a non-trivial module and therefore must be $k$-vertex-critical from Lemma~\ref{lem:module}. Further, since each $Q_j$ must induce a $P_4$-free graph, it follows that each must be a clique. Therefore, $G\in\mathcal{G}_i^{\ast}$. Similarly, if $G\in {\cal H}$, then $A_i$ is a non-trivial module that induces a $P_4$-free graph for all $i\in \{1,2,3,4,5\}$. Thus, by Lemma~\ref{lem:module}, $A_i$ is a clique for all $i\in \{1,2,3,4,5\}$ and therefore $G\in{\cal H}^{*}$.
\end{proof}

With Lemma~\ref{lem:kcritinGastorHast} in hand, the most difficult case to deal with is $k$-vertex-critical graphs in ${\cal H}^{*}$. We will show that there are no such graphs in ${\cal H}^{*}$, but it will require several lemmas first to help refine our understanding of the structure.

\begin{lemma}\label{lem:C5sinH}
If $G\in {\cal H}^{*}$ and $S\subset V(G)$ such that $G[S]\cong C_5$, then $S\subseteq \cup_{i=1}^{5}A_i$ or $S\subseteq A_6\cup A_7$.
\end{lemma}
\begin{proof}
Let $G\in {{\cal H}^{*}}$ and $S\subset V(G)$ such that $G[S]\cong C_5$. Clearly, if $S\cap A_7=\emptyset$, then $S\subseteq \cup_{i=1}^{5}A_i$. So suppose $S\cap A_7\neq \emptyset$ and let $a_7\in S\cap A_7$. Since $A_7$ is the disjoint union of cliques
and each component is a module whose neighbours are contained in $A_6$, we must have a vertex $a_6\in A_6 \cap S$. If there is a vertex $a_7'\in S\cap A_7$ such that $a_7'\sim a_7$, then $\{a_6,a_7',a_7\}$ induces a triangle, contradicting $G[S]\cong C_5$. Since $a_7$ must have another neighbour, we must have $a_6'\in A_6$ such that $a_6\neq a_6'$, $a_6'\nsim a_6$, and $a_6'\sim a_7$. Now if $S\not\subseteq A_6\cup A_7$, then there must be a vertex $s\in S\cap (A_1\cup A_3\cup A_4)$. However, now $\{s,a_6,a_6',a_7\}$ induces a $C_4$ in $G$, contradicting $G[S]\cong C_5$. Therefore, $S\subseteq A_6\cup A_7$.
\end{proof}

\begin{lemma}\label{lem:HisG4G10free}
For all $G\in {\cal H}^{*}$, $G$ is $(G_4,G_5,G_6,G_7,G_8,G_9,G_{10})$-free.
\end{lemma}
\begin{proof}

By way of contradiction, suppose $G\in {{\cal H}^{*}}$ and $S\subset V(G)$ such that $G[S]\cong G_i$ for some $i\in \{4,5,\ldots, 10\}$. Thus, there is a $D\subset S$ such that $D$ induces a $C_5$ and it follows from Lemma~\ref{lem:C5sinH} that $D\subseteq \cup_{i=1}^{5}A_i$ or $D\subseteq A_6\cup A_7$.

Suppose that $D\subseteq \cup_{i=1}^{5}A_i$. Then each of the five vertices of $D$ belongs to a distinct $A_i$ for $i = 1,2, \ldots, 5$, since each of these $A_i$'s induces a clique. From the structure of graphs in ${\cal H}^{*}$, every neighbour of a vertex in $D$ must be in $A_6$. Since $A_6$ is a module in $G-A_7$, it must be the case that $N(u)\cap D=N(v)\cap D$ for all $u,v\in S\setminus D$ such that $N(u)\cap D\neq\emptyset$ and $N(v)\cap D\neq\emptyset$.
However, we find that this is not the case for any induced $C_5$ of $G_i$ for any $i\in\{4,5,\ldots,10\}$. For example, for $G_4$, $C=\{v_1,v_2,v_3,v_4,v_5\}$ induces the only $C_5$ in the graph and $N(v_6)\cap C=\{v_1,v_3,v_4\}\neq N(v_7)\cap C=\{v_1,v_2,v_4\}$.  Therefore, $S$ cannot induce $G_i$ for $i = 4, \dots 10$, a contradiction. Thus, $D\not\subseteq \cup_{i=1}^{5}A_i$.

Thus, we may assume that $D\subseteq A_6\cup A_7$. Since each component of $A_7$ is a module and a clique, if $x,y\in D \cap A_7$ such that $x\sim y$, then $N[x]=N[y]$, otherwise a triangle is induced. But  for all $j\in \{4,5,\ldots, 10\}$ and $u,v\in G_j$ such that $u\sim v$, $N[u]\neq N[v]$. Therefore, $D\cap A_7$ induces a stable set in $G$. For any stable set $T$ of $G_j$, it is easily verifiable that $G_j-T$ contains an induced $P_4$. Therefore, $S\setminus A_7$ must contain a subset $P$ such that $G[P]\cong P_4$. Note that $A_6$ is $P_4$-free since $G\in{\cal H}^{*}$, so $P\not\subseteq A_6$ and therefore $S\not\subseteq A_7\cup A_6$. Further, $|P\cap A_7|=1$ since $A_6$ is complete to $A_1\cup A_3\cup A_4$. Thus $P$ must either follow the path $A_5-A_1-A_6-A_3$ or $A_2-A_1-A_6-A_4$. In either case, $P$ starts at a vertex with no neighbours in $D$ (i.e., an induced $C_5$) and ends at a vertex not in $D$ but with a neighbour in $D$. Thus, $S$ must induce $G_7$ as it can readily be checked that this is the only graph with vertices that have no neighbours on an induced $C_5$. In fact, $G_7$ has exactly two such vertices: $v_7$ has no neighbours on the $C_5$ induced by $C=\{v_2,v_3,v_6,v_4,v_8\}$ and $v_8$ has no neighbours on the $C_5$ induced by $C'=\{v_1,v_3,v_6,v_5,v_7\}$. However, all induced $P_4$'s in $G_7$ that start with $v_7$ end on a vertex in the set $\{v_6, v_8, v_3\}\subseteq C$, and all that start with $v_8$ end on a vertex in the set $\{v_6, v_7, v_3\}\subseteq C'$. This contradicts the existence of $P$. Therefore, $D\not\subseteq A_6\cup A_7$ which completes the proof.
\end{proof}

\begin{lemma}[Lemma 3.1 in \cite{Hoang2015}]\label{lem:similar}
	Let $G$ be a graph with two disjoint cliques $A$ and $B$ with $A= a_1, a_2, \ldots, a_m$ and $B = b_1, b_2 , \ldots , b_m$ such that $N(a_i) \setminus A \subseteq N(b_i) \setminus B$ for all $1 \leq i \leq m$. Then $G$ is not $k$-vertex-critical.
\end{lemma}

While we are still focused on showing there are no vertex-critical graphs in ${\cal H}^{*}$, the following lemma will be required.

\begin{lemma}\label{lem:nocritinG2G3}
There are no vertex-critical graphs in $\mathcal{G}_2\cup \mathcal{G}_3$.
\end{lemma}
\begin{proof}
Suppose $G$ is a $k$-vertex-critical graph in $\mathcal{G}_2\cup\mathcal{G}_3$ for some $k$  {and recall that $Q_i$ denotes the set of vertices replacing $v_i$ as labelled in Figure~\ref{fig:G1G7}}. By Lemma~\ref{lem:module}, it follows that $G\in \mathcal{G}_2^{\ast}\cup\mathcal{G}^{\ast}_3$. If $|Q_6| \geq |Q_5|$, then $Q_6$ contains a clique $Q'_6$ with $|Q'_6| = |Q_5|$ vertices such that $Q'_6$ and $Q_5$ that satisfy the hypothesis of Lemma~\ref{lem:similar}, so $G$ is not $k$-vertex-critical, a contradiction. So, we have $|Q_6|<|Q_5|$.

Fix a $k$-colouring of $c:V(G)\to\{1,\ldots,k\}$ of $G$ such that $x\in Q_6$ is the only vertex with colour $k$. Since $\chi(G)=k$, all colours in $\{1,\ldots,k-1\}$ appear in the neighbourhood of $x$, so $c(Q_5)\subseteq c(Q_6)\cup c(Q_3)$. Note that $|c(Q_6)|=|Q_6|$ since $Q_6$ is a clique.

Since $|Q_6|< |Q_5|$, it follows that there is a colour $c_5\in c(Q_5)\cap c(Q_3)$. Let $x_5\in Q_5$ be such that $c(x_5)=c_5$. Further, if $c(Q_1)\subseteq c(Q_4)\cup c(Q_3)$, then $k=|Q_6|+|Q_4|+|Q_3|=\chi(G[Q_3\cup Q_4\cup Q_6])$, contradicting $G$ being $k$-vertex-critical. Therefore, there is a colour $c_1\in c(Q_1)\setminus (c(Q_4)\cup c(Q_3))$. Let $x_1\in Q_1$ be such that $c(x_1)=c_1$. Thus, $c':V(G)\to\{1,2,\ldots,k\}$ defined by $c'(v)=c(v)$ for all $v\in V(G)\setminus\{x_1,x_5\}$ and $c'(x_5)=c_1$ and $c'(x_1)=c_5$ is a $k$-colouring of $G$. Finally, we can now change the colour of $x$ to $c_1$ to obtain a $(k-1)$-colouring of $G$, a contradiction. Therefore $\chi(G-x)=\chi(G)$ for all $x\in Q_6$, and since $Q_6\neq \emptyset$, this contradicts $G$ being $k$-vertex-critical. Therefore, there are no $k$-vertex-critical graphs in $\mathcal{G}_2\cup\mathcal{G}_3$.
\end{proof}

\begin{lemma}[\cite{KarthickMaffraygemcogem}]\label{lem:chromaticnumberofC5exp}
If $G\in \mathcal{G}_{1}^{\ast}$ has order $n$, then $\chi(G)=\max\left(\omega(G),\left\lceil\frac{n}{2}\right\rceil\right)$.
\end{lemma}

\begin{theorem}\label{thm:nocritinH}
There are no critical graphs in ${\cal H}$.
\end{theorem}
\begin{proof}
By induction on $k$.  For $k=1,2,3$ the result is clear. For some $k>3$, suppose there are no $k'$-vertex-critical graphs in ${\cal H}$ for all $k'<k$ and suppose $G$ is a $k$-vertex-critical in ${\cal H}$. By Lemma~\ref{lem:module}, it follows that $G\in {\cal H}^{*}$. Let $T_1,T_2,\ldots,T_m$ be the components of $A_7$, $C=\cup_{i=1}^{5}A_i$, and $q=\omega(G)$. Note that $q<k$ since $G$ is $k$-vertex-critical and not complete.
From Lemma~\ref{lem:similar}, it follows that
\begin{eqnarray}
\omega(A_6)&<&|A_2|\label{eq:A6<A2}\\
\omega(A_6)&<&|A_5|.\label{eq:A6<A5}
\end{eqnarray}

\noindent It follows from (\ref{eq:A6<A2}) and (\ref{eq:A6<A5}) that no clique in $A_1\cup A_6$ is a $q$-clique.  {We now consider two cases}.\\

\noindent\textit{{Case 1:}} {For some $c_1,c_2\in C$ and $t_i\in T_i$, $S=\{c_1,c_2,t_1,t_2,\ldots, t_m\}$ is a good stable set in $G$.}\\

Then $G-S$ must be $(k-1)$-chromatic and $\omega(G-S)<q\le k-1$. If $G-S$ does not have an induced $C_5$, then since it is $(P_5,\text{ gem})$-free, it must be perfect by the Strong Perfect Graph Theorem~\cite{Chudnovsky2006} and therefore must contain $K_{k-1}$, contradicting $\omega(G-S)<k-1$.  Thus we may assume $G - S$ contains an induced $C_5$.  Since there are no $(k-1)$-vertex-critical graphs in ${\cal H}$ by the inductive hypothesis, $G-S$ must have a $(k-1)$-vertex-critical induced subgraph $F$ that is in $\mathcal{G}_i^{\ast}$ for some $i\in \{1,\ldots,10\}$ by Theorem~\ref{thm:P5gemfreestructure}.  By Lemma~\ref{lem:HisG4G10free}, we have $F\in \mathcal{G}_{1}^{\ast}\cup \mathcal{G}_{2}^{\ast}\cup \mathcal{G}_{3}^{\ast}$. From Lemma~\ref{lem:nocritinG2G3}, we must have $F\in \mathcal{G}_{1}^{\ast}$. Since $\omega(G-S)<k-1$, Lemma~\ref{lem:chromaticnumberofC5exp} gives $k-1=\left\lceil\frac{|V(F)|}{2}\right\rceil$. By Lemma~\ref{lem:C5sinH}, we must have $V(F)\subseteq C\setminus\{c_1,c_2\}$ or $V(F)\subseteq A_6\cup A_7\setminus\{t_1,t_2\ldots, t_m\}$.

If $V(F)\subseteq C\setminus\{c_1,c_2\}$, then $G[V(F)\cup\{c_1,c_2\}]\in \mathcal{G}_1^{\ast}$ and has order $|V(F)|+2$. Therefore, by Lemma~\ref{lem:chromaticnumberofC5exp}, we have $\chi(G[V(F)\cup\{c_1,c_2\}])=k$, which contradicts $G$ being $k$-vertex-critical. So, we may assume  $V(F)\subseteq A_6\cup A_7\setminus\{t_1,t_2\ldots, t_m\}$. If there is at most one $i\in\{1,\ldots,m\}$ such that $T_i\cap V(F)\neq\emptyset$, then $A_6\cap V(F)$ must induce a graph containing a $P_4$. But then $G[A_6\cup\{x\}]$ contains an induced gem for all $x\in A_1\cup A_3\cup A_4$, contradicting $G$ being gem-free. So, for some distinct $i,j\in \{1,\ldots, m\}$, we have $T_i\cap V(F)\neq \emptyset$ and $T_j\cap V(F)\neq\emptyset$. Therefore, $G[V(F)\cup\{t_i,t_j\}]\in \mathcal{G}_1^{\ast}$ and has order $|V(F)|+2$. Hence $\chi(G[V(F)\cup\{t_i,t_j\}])=k$ by Lemma~\ref{lem:chromaticnumberofC5exp}, which contradicts $G$ being $k$-vertex-critical.\\

\noindent\textit{ {Case 2:}}  {For all $c_1,c_2\in C$ and $t_i\in T_i$, $S=\{c_1,c_2,t_1,t_2,\ldots, t_m\}$ is not a good stable set in $G$.}\\

Then since no clique in $A_1\cup A_6$ is a $q$-clique, from the proof of \cite[Theorem 8]{Chudetal2020}, we have $|A_1|=|A_3|=|A_4|=a$ and $|A_2|=|A_5|=q-a$ (for example, if $A_1\cup A_2$ is not a $q$-clique, then the set $S$ defined in the next sentence is a good stable set). Let $S=\{x_3,x_5,t_1,t_2,\ldots, t_m\}$ for some $t_i\in T_i$, $x_3\in A_3$, and $x_5\in A_5$. If $q<k-1$, then since $G$ is $k$-vertex-critical and $S$ is a stable set, we have $\chi(G-S)=k-1$ and have $\omega(G-S)\le q< k-1$ and similar to above, $G-S$ must have a $(k-1)$-vertex-critical induced subgraph $F$, such that $F\in \mathcal{G}_1^{\ast}$ and either $G[F\cup\{x_3,x_5\}]$ or $G[F\cup \{t_i,t_j\}]$ are $k$-chromatic proper subgraphs of $G$ for some $i,j\in\{1,\ldots,m\}$, contradicting $G$ being $k$-vertex-critical.

So suppose $q=k-1$. If $a=1$, then $G-(A_6\cup A_7)$ is a $k$-vertex-critical graph, a contradiction. So suppose $a>1$. Since $G$ is $k$-vertex-critical, there is a  $k$-colouring $G$ such that $x_5\in A_5$ is the only vertex with colour $k$ and $k-1$ colours appear in the neighbourhood of $x_5$. Without loss of generality, we may assume  colours $a+2,a+3,\ldots, k-1,k$ appear in $A_5$. Since $A_5$ is complete to $A_4$, the $k$-colouring must use $a$ new colours on $A_4$, so we may assume that colours $2,3,\ldots, a,a+1$ appear in $A_4$. It follows that the vertices of $A_1$ are coloured with colours $1,2,\ldots, a$. Since $A_2$ is complete to $A_1$ and $x_5$ is the only vertex coloured $k$, the $(q-a)$-clique $A_2$ is coloured with colours $a+1,a+2,\ldots,k-1$. But now, since $A_3$ is complete to both $A_4$ and $A_2$ and  colour $k$ does not appear in $A_3$, colour $1$ is the only colour that can be used to colour the $a$-clique $A_3$. However, $a>1$, so no such $k$-colouring exists.   This contradicts $G$ being $k$-vertex-critical.

Therefore, there are no critical graphs in ${\cal H}$.
\end{proof}

\begin{theorem}\label{thm:finitemanycritp5gemfree}
There are only finitely many $k$-vertex-critical $(P_5,\text{ gem})$-free graphs for all $k$.
\end{theorem}
\begin{proof}
Let $G$ be a $k$-vertex-critical $(P_5,\text{ gem})$-free graph. By Lemma~\ref{lem:kcritinGastorHast}, we know that $G=K_k$, $G\in \mathcal{G}_i^{\ast}$, or $G\in{\cal H}^{*}$. From Theorem~\ref{thm:nocritinH}, we know that $G\not\in{\cal H}^{*}$. If $G\in\mathcal{G}_i^{\ast}$, then as in~\cite{CaiGoedgebeurHuang2021}, each $Q_j$ can only be a clique of order at most $k-2$ (otherwise, since there are no isolated vertices in any of the $G_i$'s, a $k$-clique will be a proper subgraph of $G$) $|V(G)|\le (k-2)^{|V(G_i|}\le (k-2)^{9}$, and therefore there are only finitely many such graphs.
\end{proof}

\section{Complete Classification for $k=6$ and $k=7$}\label{sec:6and7}
As outlined in the proof of Theorem~\ref{thm:finitemanycritp5gemfree}, every $k$-vertex-critical ($P_5$, gem)-free graph is $K_k$ or a clique expansion of $G_i$ for some $i\in \{1,4,5,6,7,8,9,10\}$. This allowed us to implement an exhaustive computer search for  all such graphs for $k\le 7$ in Sagemath. Our search yielded exactly $19$ $6$-vertex-critical $(P_5,\text{ gem})$-free graphs and $46$ $7$-vertex-critical $(P_5,\text{ gem})$-free graphs.
The code used to generate these graphs as well as the graphs in graph6 format are available at~\cite{graphfiles}. The idea behind our generation is quite simple: for each $G_i$ for $i\in\{1,4,5,6,7,8,9,10\}$ we construct all possible nonisomorphic clique expansions with clique number at most $k-1$ and perform a brute-force test to check if each graph is $k$-vertex-critical or not. To generate all nonisomorphic clique expansions with sufficiently low clique number, we first enumerate all possible $|V(G_i)|$-tuples where each element is between $1$ and $k-2$ and these will correspond to the replacing vertex $v_j$ of $G_i$ with a clique whose order is equal to the $j$-th element in the tuple. The tuples ones that would create cliques that are too large or would produce an isomorphic graph to a previous $|V(G_i)|$-tuple are removed from the list. It should be noted that both tuple generation and the vertex-critical test are computationally expensive\footnote{The computation for $k=7$ took approximately 72 hours on a machine with an AMD Ryzen 7 4800H Processor and 16GB of RAM.}. The complete lists are included below. For these lists, let $G_1(n_1,n_2,n_3,n_5)$ ($G_7(n_1,n_2,n_3,n_5,n_6,n_7)$) denote the clique expansion of $G_1$ ($G_7$) with vertex $v_i$ replaced by a clique of order $n_i$.\\

The $19$ $6$-vertex-critical $(P_5,\text{ gem})$-free graphs are given in the following list:\\
\tiny
\begin{multicols}{2}
\begin{enumerate}[nosep]
 \item $K_6$
 \item $G_1(4, 1, 4, 1, 1)$
 \item $G_1(4, 1, 3, 2, 1)$
 \item $G_1(3, 2, 3, 2, 1)$
 \item $G_1(3, 2, 2, 3, 1)$
 \item $G_1(3, 2, 2, 2, 2)$
 \item $G_7(1, 3, 1, 3, 1, 1, 4, 2)$
 \item $G_7(1, 2, 2, 3, 1, 1, 4, 2)$
 \item $G_7(1, 1, 3, 3, 1, 1, 4, 2)$
 \item $G_7(1, 2, 2, 2, 1, 2, 4, 2)$
 \item $G_7(2, 2, 1, 2, 2, 1, 3, 3)$
 \item $G_7(2, 1, 2, 2, 2, 1, 3, 3)$
 \item $G_7(1, 1, 3, 2, 2, 1, 3, 3)$
 \item $G_7(2, 1, 2, 2, 1, 2, 3, 3)$
 \item $G_7(1, 2, 2, 2, 1, 1, 4, 3)$
 \item $G_7(1, 1, 3, 2, 1, 1, 4, 3)$
 \item $G_7(1, 2, 2, 1, 1, 2, 4, 3)$
 \item $G_7(1, 1, 3, 1, 1, 1, 4, 4)$
 \item $G_7(1, 1, 2, 1, 1, 2, 4, 4)$
\end{enumerate}
\end{multicols}

\normalsize
\vspace{5mm}
The $46$ $7$-vertex-critical $(P_5,\text{ gem})$-free graphs are given in the following list:\\
\tiny
\begin{multicols}{3}
\begin{enumerate}[nosep]
 \item $K_7$
 \item $G_1(5, 1, 5, 1, 1)$
 \item $G_1(5, 1, 4, 2, 1)$
 \item $G_1(4, 2, 4, 2, 1)$
 \item $G_1(5, 1, 3, 3, 1)$
 \item $G_1(4, 2, 3, 3, 1)$
 \item $G_1(3, 3, 3, 3, 1)$
 \item $G_1(4, 2, 2, 4, 1)$
 \item $G_1(4, 2, 3, 2, 2)$
 \item $G_1(3, 3, 3, 2, 2)$
 \item $G_1(3, 3, 2, 3, 2)$
 \item $G_7(1, 4, 1, 4, 1, 1, 5, 2)$
 \item $G_7(1, 3, 2, 4, 1, 1, 5, 2)$
 \item $G_7(1, 2, 3, 4, 1, 1, 5, 2)$
 \item $G_7(1, 1, 4, 4, 1, 1, 5, 2)$
 \item $G_7(1, 3, 2, 3, 1, 2, 5, 2)$
 \item $G_7(1, 2, 3, 3, 1, 2, 5, 2)$
 \item $G_7(2, 3, 1, 3, 2, 1, 4, 3)$
 \item $G_7(2, 2, 2, 3, 2, 1, 4, 3)$
 \item $G_7(1, 3, 2, 3, 2, 1, 4, 3)$
 \item $G_7(2, 1, 3, 3, 2, 1, 4, 3)$
 \item $G_7(1, 2, 3, 3, 2, 1, 4, 3)$
 \item $G_7(1, 1, 4, 3, 2, 1, 4, 3)$
 \item $G_7(2, 2, 2, 3, 1, 2, 4, 3)$
 \item $G_7(2, 1, 3, 3, 1, 2, 4, 3)$
 \item $G_7(2, 2, 2, 2, 2, 2, 4, 3)$
 \item $G_7(1, 2, 3, 2, 2, 2, 4, 3)$
 \item $G_7(1, 3, 2, 3, 1, 1, 5, 3)$
 \item $G_7(1, 2, 3, 3, 1, 1, 5, 3)$
 \item $G_7(1, 1, 4, 3, 1, 1, 5, 3)$
 \item $G_7(1, 3, 2, 2, 1, 2, 5, 3)$
 \item $G_7(1, 2, 3, 2, 1, 2, 5, 3)$
 \item $G_7(1, 1, 3, 3, 1, 2, 5, 3)$
 \item $G_7(2, 2, 2, 2, 2, 1, 4, 4)$
 \item $G_7(2, 1, 3, 2, 2, 1, 4, 4)$
 \item $G_7(1, 1, 4, 2, 2, 1, 4, 4)$
 \item $G_7(2, 2, 2, 2, 1, 2, 4, 4)$
 \item $G_7(2, 1, 3, 2, 1, 2, 4, 4)$
 \item $G_7(1, 1, 3, 2, 2, 2, 4, 4)$
 \item $G_7(1, 2, 3, 2, 1, 1, 5, 4)$
 \item $G_7(1, 1, 4, 2, 1, 1, 5, 4)$
 \item $G_7(1, 2, 3, 1, 1, 2, 5, 4)$
 \item $G_7(1, 2, 2, 2, 1, 2, 5, 4)$
 \item $G_7(1, 1, 3, 2, 1, 2, 5, 4)$
 \item $G_7(1, 1, 4, 1, 1, 1, 5, 5)$
 \item $G_7(1, 1, 3, 1, 1, 2, 5, 5)$
\end{enumerate}
\end{multicols}
\normalsize

\section{Conclusion}\label{sec:conclusion}
We provided an alternate, stronger proof that there are only finitely many $k$-vertex-critical $(P_5,\text{ gem})$-free graphs for all $k$. We used our proof to completely classify all such graphs for $k\le 7$, extending the lists given in~\cite{CaiGoedgebeurHuang2021} for $k\le 5$. While our results give the most thorough structural description of vertex-critical $(P_5,\text{ gem}$)-free graphs, there are still a few open problems.

\begin{problem}\label{prob:8andup}
Classify all $k$-vertex-critical $(P_5,\text{ gem})$-free graphs for fixed $k\ge 8$.
\end{problem}

To our knowledge, complete lists of $k$-vertex-critical $(H_1,H_2)$-free graphs for $k\ge 8$ are only known for (gem, co-gem)-free graphs, where they are given for all $k\le 16$~\cite{AbuadasCameronHoangSawada2022}. Therefore, it is of interest to generate these complete lists for other families of graphs. Our methods from Section~\ref{sec:6and7} work in theory to find these classifications, however, limitations on computational resources prevent us from generating these lists for now. One approach to reduce the search space of our methods in Section~\ref{sec:6and7} is to show that $G_i$ has no $k$-vertex-critical clique expansions for some $1\le i\le 10$. In fact, for $k\le 7$, there are no clique expansions of any $G_i$ other than for $i=1$ or $i=7$ and we proved in Lemma~\ref{lem:nocritinG2G3} that there are none for $i=2$ and $i=3$ for any $k$. This leads to the following open problem:

\begin{problem}\label{prob:cliqueexpsofotherGis}
Determine if there are any $k$-vertex-critical clique expansions of $G_i$ for $i\neq 1$ and $i\neq 7$.
\end{problem}

If there are no such graphs for Open Problem~\ref{prob:cliqueexpsofotherGis}, then the $G_1$ clique expansions have already been classified from the recent work on (gem, co-gem)-free graphs~\cite{AbuadasCameronHoangSawada2022}, so to classify all $k$-vertex-critical $(P_5,\text{ gem})$-free graphs would require only searching all clique expansions of $G_7$.

\section*{Acknowledgements}
This work was supported by the Canadian Tri-Council Research Support Fund. Both authors were each supported by individual NSERC Discovery Grants.

\bibliographystyle{abbrv}