Heavy Flavour Physics and CP Violation at LHCb: a Ten-Year Review

Abstract: Heavy flavour physics provides excellent opportunities to indirectly search for new physics at very high energy scales and to study hadron properties for deep understanding of the strong interaction. The LHCb experiment has been playing a leading role in the study of heavy flavour physics since the start of the LHC operations about ten years ago, and made a range of high-precision measurements and unexpected discoveries, which may have far-reaching implications on the field of particle physics. This review highlights a selection of the most influential physics results on CP violation, rare decays, and heavy flavour production and spectroscopy obtained by LHCb using the data collected during the first two operation periods of the LHC. The upgrade plan of LHCb and the physics prospects are also briefly discussed.

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

With the discovery of the Higgs boson at the Large Hadron Collider (LHC) in 2012~\cite{Bruening:782076,ATLAS:2012yve,CMS:2012qbp}, all fundamental particles expected in the Standard Model (SM) of particle physics are found and the SM is finally completed.  The SM has achieved tremendous success in explaining experimental results in high-energy physics.
However, there are still many key questions that are  not answered in the SM, such as the mechanism to generate the matter-antimatter asymmetry in the Universe, the origin of the three generations of fermions and their mixing,
and the nature of dark matter and dark energy.
It is commonly believed that new physics (NP) beyond the SM should exist at or above the $\kern -0.25em\aunit{Te\kern -0.1em V}\xspace$ energy scale. Flavour physics can provide a unique approach to indirectly probe NP at energy scales far above  \aunit{Te\kern -0.1em V}\xspace via precision study of charge-parity  (CP) violation and rare phenomena, complementary to the direct search for new particles and interactions at the energy frontier.
Flavour physics also serves as a natural laboratory to test quantum chromodynamics (QCD), the theory of the strong interaction, via measurements of hadron production and spectroscopy.
The Large Hadron Collider beauty (LHCb) experiment has been playing a leading role in the study of heavy-flavour physics since the start of the LHC, and has made a series of discoveries and improvements in CP violation, rare decays, and hadron production and spectroscopy.
This review aims to present a selection of the high-impact physics results on the above subjects from the LHCb experiment, and to briefly discuss the future prospects.
Due to the limited space, not all interesting results can be covered here.
For a complete list of all LHCb physics results, please refer to the official LHCb summary~\cite{LHCbpapers}.

The LHCb detector~\cite{LHCb-DP-2008-001} is optimised for the study of the decays of heavy-flavour hadrons, \mbox{\itshape i.e.}\xspace hadrons containing heavy quarks ({\ensuremath{\Pb}}\xspace or {\ensuremath{\Pc}}\xspace quarks, often collectively referred to as $Q$).
In proton-proton ($pp$) collisions at LHC energies, the {\ensuremath{\bquark\bquarkbar}}\xspace pairs are produced dominantly through the gluon fusion process $gg\ensuremath{\rightarrow}\xspace {\ensuremath{\bquark\bquarkbar}}\xspace$.
Due to the large Lorentz boost along the proton beam in the laboratory frame, the {\ensuremath{\Pb}}\xspace and {\ensuremath{\overline \bquark}}\xspace quarks
generated in a pair are highly correlated in their momentum directions,
either both in the forward region or both in the backward region in the majority of  cases.
In order to take advantage of this characteristic of the {\ensuremath{\bquark\bquarkbar}}\xspace pair at the LHC,
 the LHCb detector is designed to have a forward geometry
as  shown in Fig.~\ref{fig:01:01:LHCbDetector} to cover the forward region of $pp$ collisions.
\begin{figure}[!b]
 \begin{center}
  \includegraphics[width=0.70\textwidth]{figs/introduction/Fig1LHCbDetector.pdf}
 \end{center}
 \caption{
Layout of the LHCb detector~\cite{LHCb-DP-2008-001}.
 }
 \label{fig:01:01:LHCbDetector}
\end{figure}
The flavour of the  {\ensuremath{\Pb}}\xspace-hadron
under study can usually  be tagged by the other {\ensuremath{\Pb}}\xspace-hadron that is also inside the LHCb acceptance.
This enhances the potential of the LHCb experiment in the study of CP violation and mixing with ${\ensuremath{\B^0}}\xspace$ and ${\ensuremath{\B^0_\squark}}\xspace$ decays.

Excellent vertex and momentum resolution, particle identification \mbox{\itshape etc.}\xspace are key ingredients for  flavour physics measurements.
A silicon vertex locator (VELO) surrounding
the $pp$
collision region is used to precisely determine  the primary interaction vertices (PVs) and the displaced secondary/tertiary vertices (SVs/TVs) formed by the decay products of heavy-flavour hadrons.
The VELO system offers decay time measurements with a typical resolution better than $50\aunit{fs}$,  thus allows the LHCb experiment to make precision measurement of hadron lifetimes and resolve the  fast ${\ensuremath{\B^0_\squark}}\xspace$-${\ensuremath{\Bbar{}^0_\squark}}\xspace$ oscillation, which has a period of about 350\aunit{fs}.
The clear separation of SVs from PVs also allows for substantial suppression of the combinatorial background, which is extremely high in $pp$ collisions.
In addition to the VELO, the LHCb tracking system includes also four trackers, namely the TT and T1-T3 stations, located  upstream and downstream of the dipole magnet, respectively.
Together with a magnet that has a bending power of $4\,\mathrm{Tm}$, the tracking system provides precise measurements of the momenta of charged particles. The momentum resolution ($\Delta p/p$) is typically 0.5\% for
low momentum tracks and 1.0\% for track momentum up to $200\ensuremath{\aunit{Ge\kern -0.1em V\!/}c}\xspace$.
The mass resolution for ${\ensuremath{\Pb}}\xspace$-hadrons can be as good as $8\ensuremath{\aunit{Me\kern -0.1em V\!/}c^2}\xspace$, precise enough to distinguish decays of ${\ensuremath{\B^0}}\xspace$ and ${\ensuremath{\B^0_\squark}}\xspace$ mesons to the same final state.

There are two ring-imaging Cherenkov detectors (RICH1 and RICH2)  used to identify charged hadrons in the momentum range $p\in(2,100)\ensuremath{\aunit{Ge\kern -0.1em V\!/}c}\xspace$.
The  RICH detectors are very powerful at suppressing misidentification background for {\ensuremath{\Pb}}\xspace-hadron  decays to final states containing
charged kaons, pions or protons.
The muon system provides excellent muon identification and is essential for reconstruction of decays with muons in the final states.
The electromagnetic calorimeter is used  for photon and electron reconstruction and identification. Together with the hadronic calorimeter, it also provides information for event trigger.

The trigger system is crucial for the success of the LHCb experiment.
The hardware trigger at the first level
reduces the data rate from $40\ensuremath{\aunit{MHz}}\xspace$ down to $1\ensuremath{\aunit{MHz}}\xspace$, at
which point the flexible software-based trigger takes over to further reduce the rate to around $12\ensuremath{\aunit{kHz}}\xspace$ for offline processing and  analysis.
The ability to sustain such large rates enables the LHCb experiment to trigger with high efficiency on decay processes across a wide range of final states
and to provide large data samples for study of exclusive rare decay modes as well as for inclusive data mining.
Of particular importance is the LHCb muon trigger, which allows events containing one muon or two muons to be selected with greater than 95\% efficiency. This is ideal for studying $b$-hadron decays to
${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace$ or ${\ensuremath{\Ppsi{(2S)}}}\xspace$ mesons that further decay to ${\ensuremath{\Pmu^+}}\xspace\mun$ pairs and semileptonic decays into muons.
Meanwhile, there are also highly efficient triggers that are purely based on the multi-body topology of the final state hadrons. These triggers are important for  reconstructing decays of heavy-quark hadrons to final states without muons.

The LHCb experiment  collected a data sample corresponding to an integrated luminosity of  3\ensuremath{\fb^{-1}}\xspace in $pp$ collisions  at centre-of-mass energies $\ensuremath{\protect\sqrt{s}}\xspace=7$ and $8\aunit{Te\kern -0.1em V}\xspace$ from 2011 to 2012 (Run 1), and another sample of 6\ensuremath{\fb^{-1}}\xspace at $\ensuremath{\protect\sqrt{s}}\xspace=13\aunit{Te\kern -0.1em V}\xspace$ from 2015 to 2018 (Run 2).
Results discussed in this review are based on either full Run 1 and Run 2 data samples or a subset.
Since the cross-sections for $c$- and $b$-quark production in  $pp$ collisions at  13\aunit{Te\kern -0.1em V}\xspace  are about twice of the cross-sections at 7 and 8\aunit{Te\kern -0.1em V}\xspace, and  the trigger scheme for Run 2 has also been improved compared with Run 1, the number of recorded charm and beauty decays available for physics analysis is more than four times higher in the Run 2 data  than in Run 1 data.
In total,
more than $10^{11}$ {\ensuremath{\Pb}}\xspace-hadrons and $10^{12}$ {\ensuremath{\Pc}}\xspace-hadrons have been produced within the LHCb detector acceptance. The typical trigger and selection efficiencies are of the order of $10^{-3}$ to $10^{-2}$
for  decays only to  charged particles
and  $10^{-4}$ to $10^{-3}$ for decays to final states  involving photons, ${\ensuremath{\pion^0}}\xspace$, ${\ensuremath{\PLambda}}\xspace^0$ or ${\ensuremath{\kaon^0_{\mathrm{S}}}}\xspace$ particles.
The enormous $b$- and $c$-hadron samples
 form the basis of
precision measurements of CP violation, exploration of rare decays, and searches for new hadrons.

This review is structured as follows, with each section
covering a different subject.
Recent results of heavy-flavour production and spectroscopy are shown in Section~2. For heavy-flavour production, recent results of associated production and the studies in heavy-ion collisions are shown; for spectroscopy, the results of conventional hadrons and exotic hadrons are summarised.
Section~3 discusses rare $B$-meson decays. The results on purely leptonic $B$-meson decays, semileptonic $b\ensuremath{\rightarrow}\xspace{\ensuremath{\Ps}}\xspace{\ensuremath{\ell^+}}\xspace{\ensuremath{\ell^-}}\xspace$ decays, and radiative $b\ensuremath{\rightarrow}\xspace{\ensuremath{\Ps}}\xspace\gamma$ decays are shown.\footnote{Charge conjugation is implied throughout unless stated otherwise.}
The observed anomalies in angular distributions and  lepton flavour universality tests are discussed in details.
Section~4 presents the latest results on CP violation in the beauty sector.
Emphasis is put on the progresses that have been made for precision test of the Cabibbo-Kobayashi-Maskawa (CKM) mechanism, such as significant improvement in the determination of the parameters $\gamma$, $\beta$, $\beta_s$, ${\ensuremath{V_{\uquark\bquark}^{\phantom{\ast}}}}\xspace$, ${\ensuremath{V_{\cquark\bquark}^{\phantom{\ast}}}}\xspace$, and $\Delta m_{s/d}$.
Section~5 provides recent results of CP violation in the charm sector, including those for charm mixing and the observation of CP violation in ${\ensuremath{\D^0}}\xspace$ decays.
The final section provides a picture of the LHCb upgrade, as well as a brief summary of the content in this review. The main goals and modifications to the detector in Upgrade I and Upgrade II are introduced, and the prospects of some key measurements are presented.

\clearpage
\section{Heavy-flavour production and spectroscopy}
\label{sec:production}

Rich information on QCD dynamics can be deciphered from measurements of heavy-flavour production ~\cite{Cacciari:2012ny,Andronic:2015wma,Zhang:2021any,Chen:2021obo,Chen:2021tmf,Wang:2020jjq,Chen:2020ajc,Yang:2019lls,He:2019qqr,Butenschoen:2011yh,Lansberg:2016deg,Zhang:2014ybe,Ma:2012hh,Feng:2018ukp,Wang:2014vsa,Tang:2014tga,Sun:2018yam,Li:2017qmj,Han:2014jya,Han:2014kxa,Zhang:2020atv,Liu:2019iml,Ma:2017xno,Sun:2014gca,Ma:2010jj,Ma:2010yw,Li:2009zu,Ma:2010vd,Zhang:2009ym,Shao:2012iz,Chang:2003cr,Niu:2018tvo,He:2017had,Zhang:2021wjj,Chen:2019ykv,Hu:2021gdg}
and studies of heavy-flavour spectroscopy~\cite{Brambilla:2019esw,Chen:2016qju,Guo:2017jvc,Swanson:2006st,Olsen:2017bmm,Liu:2019zoy,Chen:2016spr,Liu:2013waa,Guo:2019twa,Yuan:2018inv,Liu:2009qhy,Dong:2021bvy,Shi:2018rhk,Ali:2017jda,Esposito:2014rxa,Maiani:2014aja,Peng:2020xrf,Yang:2020atz,Richard:2021nvn,Karliner:2017qhf,Richard:2016eis,Wu:2021tzo}.
The mass of a heavy quark, which is much larger than the nonperturbative QCD scale $\Lambda_\mathrm{QCD}$,  provides an energy scale that allows for perturbative calculation of heavy-quark production.
The production of $Q\overline{Q}$ pairs
is predominantly  in the initial stage of the collision, thus can be used to probe properties of the colliding system and the possibly created QCD medium~\cite{Wicks:2005gt,Zhou:2014kka}.
The presence of heavy quark(s) also provides practical benefits for theoretical and experimental studies of spectroscopy.
Heavy quarks are approximately nonrelativistic in hadrons, which makes it possible to simplify theoretical calculations.
The large mass and the weak decay of heavy-flavour hadrons offer essential features, such as decay products with high transverse momentum, \ensuremath{p_{\mathrm{T}}}\xspace, and vertices displaced from the PVs, which can be exploited to reject the huge QCD background at hadron colliders.

As discussed in Section~\ref{sec:introduction},
the excellent performance of the LHCb detector, owing to the dedicated design for heavy-flavour hadrons, enables LHCb to make great achievements in the study of heavy-flavour production and spectroscopy, \mbox{\itshape e.g.}\xspace the observation of pentaquark states and the doubly charmed baryon ${\ensuremath{\PXi}}\xspace_{cc}^{++}$.
In this section, relevant results from the LHCb experiment are reviewed, with a focus on recent developments.

\subsection{Production}
A summary of LHCb production measurements for open heavy-flavour hadrons, heavy quarkonia and pairs of heavy-flavour hadrons at LHCb~\cite{LHCb-CONF-2019-004,LHCb-PAPER-2011-003,LHCb-PAPER-2011-013,LHCb-PAPER-2011-019,LHCb-PAPER-2011-030,LHCb-PAPER-2011-034,LHCb-PAPER-2011-036,LHCb-PAPER-2011-043,LHCb-PAPER-2011-045,LHCb-PAPER-2012-003,LHCb-PAPER-2012-028,LHCb-PAPER-2012-039,LHCb-PAPER-2012-041,LHCb-PAPER-2013-004,LHCb-PAPER-2013-028,LHCb-PAPER-2013-052,LHCb-PAPER-2013-066,LHCb-PAPER-2014-004,LHCb-PAPER-2014-015,LHCb-PAPER-2014-029,LHCb-PAPER-2014-031,LHCb-PAPER-2014-040,LHCb-PAPER-2014-050,LHCb-PAPER-2015-016,LHCb-PAPER-2015-032,LHCb-PAPER-2015-037,LHCb-PAPER-2015-041,LHCb-PAPER-2015-045,LHCb-PAPER-2015-046,LHCb-PAPER-2015-058,LHCb-PAPER-2016-042,LHCb-PAPER-2016-057,LHCb-PAPER-2016-064,LHCb-PAPER-2017-014,LHCb-PAPER-2017-015,LHCb-PAPER-2017-037,LHCb-PAPER-2018-002,LHCb-PAPER-2018-021,LHCb-PAPER-2018-035,LHCb-PAPER-2018-047,LHCb-PAPER-2018-048,LHCb-PAPER-2018-049,LHCb-PAPER-2019-024,LHCb-PAPER-2019-033,LHCb-PAPER-2019-035,LHCb-PAPER-2020-010,LHCb-PAPER-2020-023,LHCb-PAPER-2020-046,LHCb-PAPER-2020-048} are listed in Table~\ref{sec:production:tab:productionHF}.
Inclusive hadroproduction of open heavy-flavour hadrons ($H_Q$) factorises into
three components in perturbative QCD (pQCD) calculations: the parton distribution function (PDF) in the two initial projectiles $f_{i,j}$, the parton level cross-section $\sigma_{ij\ensuremath{\rightarrow}\xspace Q+X}$ for a heavy-quark $Q$ production, and
the heavy-quark fragmentation function $D_{Q\ensuremath{\rightarrow}\xspace H_Q}$~\cite{Pumplin:2002vw}.  Differential cross-section is expressed as
\begin{equation}
    \ensuremath{\mathrm{d}} \sigma_{AB\ensuremath{\rightarrow}\xspace H_Q} = \sum_{i,j}(f_i^A\otimes f_j^B )\otimes\ensuremath{\mathrm{d}}\sigma_{ij\ensuremath{\rightarrow}\xspace Q+X}\otimes D_{Q\ensuremath{\rightarrow}\xspace H_Q},\label{sec:prod:heavyflavour}
\end{equation}
where the indices $i,j$ run over all possible parton species, and at LHC energies heavy-flavour production is dominated by gluons.
The PDF and fragmentation function include nonperturbative effects, and can be determined from a global fit of available
data~\cite{Pumplin:2002vw,Andronic:2015wma}.
The results of open charm and beauty production at LHCb are consistent with pQCD models, for example the calculation based on
fixed-order plus next-to-leading logs (FONLL)~\cite{Cacciari:1998it}.
It turns out that the LHCb results on charm and beauty cross-sections have a better precision than theoretical calculations~\cite{LHCb-PAPER-2015-041,LHCb-PAPER-2017-037,LHCb-PAPER-2016-042}, and can be used to reduce the
uncertainties on gluon PDF, in particular in the small Bjorken-$x$ region, $x\lsim 10^{-5}$~\cite{Gauld:2016kpd}.

\begin{table}[!tb]
\caption{LHCb measurements of production cross-sections (or ratios) for various heavy hadrons in $pp$ and $p$Pb collisions at different centre-of-mass energies.}
{\small
    \begin{tabular}{m{0.06\textwidth}|m{0.14\textwidth}|m{0.23\textwidth}|m{0.23\textwidth}|m{0.19\textwidth}}
        \hline
        \multirow{2}{*}{System}  &\multicolumn{4}{c}{$\sqrt{s}_{(NN)}$ ($\kern -0.25em\aunit{Te\kern -0.1em V}\xspace$)}\\
        \cline{2-5}
        & \hspace{0.5cm}$5\,(2.76)$ &  \hspace{1.5cm}$7$ & $\hspace{1.2cm}8\,(8.16)$ & \hspace{1.5cm}$13$\\
        \hline
            $pp$  &
           \mbox{${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace$\cite{LHCb-PAPER-2012-039}},
           \mbox{$\PUpsilon$\cite{LHCb-PAPER-2013-066}},
           \mbox{$D$\cite{LHCb-PAPER-2016-042}}
            &
           \mbox{${\ensuremath{\Peta_\cquark}}\xspace$\cite{LHCb-PAPER-2014-029}},
           \mbox{${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace$\cite{LHCb-PAPER-2011-003}},
           \mbox{$\chi_{c1}(3872)$\cite{LHCb-PAPER-2011-034}},
           \mbox{${\ensuremath{\Pchi_\cquark}}\xspace/{\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace$\cite{LHCb-PAPER-2011-030}},
            \mbox{${\ensuremath{\Pchi_{\cquark 2}}}\xspace/{\ensuremath{\Pchi_{\cquark 1}}}\xspace$\cite{LHCb-PAPER-2011-019,LHCb-PAPER-2013-028}},
           \mbox{${\ensuremath{\Ppsi{(2S)}}}\xspace$\cite{LHCb-PAPER-2011-045,LHCb-PAPER-2018-049}},
           \mbox{$\PUpsilon$\cite{LHCb-PAPER-2011-036,LHCb-PAPER-2015-045}},
           \mbox{${\ensuremath{\Pchi_{\bquark 2}}}\xspace/{\ensuremath{\Pchi_{\bquark 1}}}\xspace$\cite{LHCb-PAPER-2014-040,LHCb-PAPER-2014-031}},
           \mbox{$D,{\ensuremath{\Lz^+_\cquark}}\xspace$\cite{LHCb-PAPER-2012-041}},
           \mbox{$B$\cite{LHCb-PAPER-2011-043,LHCb-PAPER-2013-004,LHCb-PAPER-2017-037,LHCb-PAPER-2020-046}},
           \mbox{${\ensuremath{\Lz^0_\bquark}}\xspace$\cite{LHCb-PAPER-2014-004,LHCb-PAPER-2015-032}},
           \mbox{${\ensuremath{\Xires_\bquark}}\xspace$\cite{LHCb-PAPER-2018-047}},
           \mbox{${\ensuremath{\B_\cquark^+}}\xspace$\cite{LHCb-PAPER-2012-028}},
           \mbox{${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace\jpsi$\cite{LHCb-PAPER-2011-013}},
           \mbox{$\PUpsilon D$\cite{LHCb-PAPER-2015-046}},
           \mbox{${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace D, DD, D{\ensuremath{\Dbar{}^0}}\xspace$\cite{LHCb-PAPER-2012-003}}
            &
           \mbox{${\ensuremath{\Peta_\cquark}}\xspace$\cite{LHCb-PAPER-2014-029}},
           \mbox{${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace$\cite{LHCb-PAPER-2015-016}},
           \mbox{${\ensuremath{\Pchi_{\bquark 2}}}\xspace/{\ensuremath{\Pchi_{\bquark 1}}}\xspace$\cite{LHCb-PAPER-2014-040,LHCb-PAPER-2014-031}},
           \mbox{$\chi_{c1}(3872)/{\ensuremath{\Ppsi{(2S)}}}\xspace$\cite{LHCb-PAPER-2020-023}},
           \mbox{$\PUpsilon$\cite{LHCb-PAPER-2015-016,LHCb-PAPER-2015-045}},
           \mbox{${\ensuremath{\B^0_\squark}}\xspace/{\ensuremath{\B^0}}\xspace$\cite{LHCb-PAPER-2020-046}},
           \mbox{${\ensuremath{\Lz^0_\bquark}}\xspace$\cite{LHCb-PAPER-2015-032}},
           \mbox{${\ensuremath{\Xires_\bquark}}\xspace$\cite{LHCb-PAPER-2018-047}},
           \mbox{$\PUpsilon D$\cite{LHCb-PAPER-2015-046}},
           \mbox{${\ensuremath{\B_\cquark^+}}\xspace$\cite{LHCb-PAPER-2014-050}}
            &
           \mbox{${\ensuremath{\Peta_\cquark}}\xspace$\cite{LHCb-PAPER-2019-024}},
           \mbox{${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace$\cite{LHCb-PAPER-2015-037,LHCb-PAPER-2016-064}},
           \mbox{${\ensuremath{\Ppsi{(2S)}}}\xspace$\cite{LHCb-PAPER-2018-049}},
           \mbox{$\PUpsilon$\cite{LHCb-PAPER-2018-002}},
           \mbox{$D$\cite{LHCb-PAPER-2015-041}},
           \mbox{${\ensuremath{\Xires^{++}_{\cquark\cquark}}}\xspace$\cite{LHCb-PAPER-2019-035}},
           \mbox{${\ensuremath{\Bu}}\xspace$\cite{LHCb-PAPER-2017-037}},
           \mbox{${\ensuremath{\Xires_\bquark}}\xspace$\cite{LHCb-PAPER-2018-047}},
           \mbox{${\ensuremath{\B_\cquark^+}}\xspace$\cite{LHCb-PAPER-2019-033}},
           \mbox{${\ensuremath{\B^0_\squark}}\xspace/{\ensuremath{\B^0}}\xspace$\cite{LHCb-PAPER-2020-046}},
           \mbox{${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace\jpsi$\cite{LHCb-PAPER-2016-057}}
            \\
        \hline
            $p$Pb &
            \mbox{${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace$\cite{LHCb-PAPER-2013-052}},
            \mbox{${\ensuremath{\Ppsi{(2S)}}}\xspace$\cite{LHCb-PAPER-2015-058}},
            \mbox{$\PUpsilon$\cite{LHCb-PAPER-2014-015}},
            \mbox{${\ensuremath{\D^0}}\xspace$\cite{LHCb-PAPER-2017-015}},
            \mbox{${\ensuremath{\Lz^+_\cquark}}\xspace$\cite{LHCb-PAPER-2018-021}}
            &                 &
           \mbox{${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace$\cite{LHCb-PAPER-2017-014}},
           \mbox{$\PUpsilon$\cite{LHCb-PAPER-2018-035}},
           \mbox{${\ensuremath{\D^0}}\xspace$\cite{LHCb-CONF-2019-004}},
           \mbox{${\ensuremath{\Bu}}\xspace, {\ensuremath{\B^0}}\xspace, {\ensuremath{\Lz^0_\bquark}}\xspace$\cite{LHCb-PAPER-2018-048}},
           \mbox{${\ensuremath{\Pchi_{\cquark 2}}}\xspace/{\ensuremath{\Pchi_{\cquark 1}}}\xspace$\cite{LHCb-PAPER-2020-048}},
           \mbox{$DD,D{\ensuremath{\Dbar{}^0}}\xspace$\cite{LHCb-PAPER-2020-010}}
            &\\
        \hline
\end{tabular}
\label{sec:production:tab:productionHF}
}
\end{table}

\begin{figure}[tb]
\begin{center}
    \includegraphics[width=0.37\textwidth]{figs/production/JpsiProduction.pdf}
     \includegraphics[width=0.58\textwidth]{figs/production/JpsiPolarization.pdf}
    \caption{(Left) Differential cross-section of ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace$ in $pp$ collisions at $\ensuremath{\protect\sqrt{s}}\xspace=13\aunit{Te\kern -0.1em V}\xspace$, taken from Ref.~\cite{LHCb-PAPER-2015-037}. (Right) ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace$ polarisation in $pp$ collisions at $\ensuremath{\protect\sqrt{s}}\xspace=7\aunit{Te\kern -0.1em V}\xspace$, taken from Ref.~\cite{LHCb-PAPER-2013-008}. The measurements are compared with NRQCD~\cite{Chao:2012iv,Butenschoen:2012qh,Gong:2012ug} and colour singlet calculations~\cite{Butenschoen:2012qh}.}
\label{sec:prod:fig:JpsiProduction}
\end{center}
\end{figure}

For quarkonium production,  assumptions have to be put on
how heavy-quark pairs, produced with various possible colour, spin and parity configurations,  transform into specific colourless
quarkonia~\cite{Ma:2014mri,Zhang:2013igb,Brambilla:2010cs,Brambilla:2004jw}. Cross-section measurements favour calculations using the nonrelativistic QCD (NRQCD)
framework~\cite{Bodwin:1994jh}, as shown on the right of Fig.~\ref{sec:prod:fig:JpsiProduction}, for ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace$ production in $pp$ collisions at $\ensuremath{\protect\sqrt{s}}\xspace=13\aunit{Te\kern -0.1em V}\xspace$. The NRQCD framework
introduces long distance matrix elements (LDMEs) as model parameters~\cite{Ma:2016exq} to account for transition probabilities from heavy-quark pairs to quarkonia. LDMEs are
assumed to be independent of quarkonium production environments and kinematics, and are fixed by matching predicted $\ensuremath{p_{\mathrm{T}}}\xspace$ spectrum to data.
The polarisation of heavy quarkonia is another observable sensitive to the $Q\overline{Q}$ production mechanism and LDMEs.
Inconsistencies are observed between LHCb data and theoretical calculations on the $\psi$ and $\PUpsilon$
polarisations~\cite{LHCb-PAPER-2017-028,LHCb-PAPER-2013-067,LHCb-PAPER-2013-008}. Only a level of 10\% level or smaller polarisation
is observed in the LHCb acceptance, in contrast to a dominant transverse polarisation predicted by NRQCD~\cite{Chao:2012iv,Gong:2012ug,Butenschoen:2012qh}. The measurement of ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace$ polarisation is shown on the right of Fig.~\ref{sec:prod:fig:JpsiProduction}. Even though the discrepancy can be reduced by tuning the LDMEs, a coherent description of production
cross-section and polarisation is still a difficult theoretical problem~\cite{Gong:2012ug,Shao:2014yta,Han:2014jya,Feng:2018ukp,Shao:2014yfa,Butenschoen:2012px,Ma:2018qvc,Shao:2014fca,Shao:2012fs}. If only the $Q\overline{Q}$ state that has the same quantum number as the final quarkonium is considered, the NRQCD framework reduces to the colour singlet model, which underestimate production cross-sections~\cite{Butenschoen:2012qh} and disagrees with data on $\psi$ and $\PUpsilon$ polarisations~\cite{LHCb-PAPER-2013-008,LHCb-PAPER-2013-067,LHCb-PAPER-2017-028}.

\subsubsection{Associated production}
Recent LHCb production measurements focus on associated production of multiple heavy flavours and quantities
probing properties of QCD matter.
Associated heavy-flavour production provides an approach to study the multiple parton interactions (MPIs).
MPIs are sensitive to correlations between partons in space, momentum, flavour, colour, spin \mbox{\itshape etc.}\xspace inside  the colliding projectiles~\cite{Chapon:2020heu,Kom:2011bd,Shao:2019qob,He:2007te,Lansberg:2013qka,Shao:2020kgj}.
Usually in an MPI process each parton-scattering is assumed to be independent from each other.
Under this assumption, the cross-section for associated production of $AB$ through a double-parton-scattering (DPS) process is related to the single inclusive production of $A$ and $B$ as~\cite{Chapon:2020heu}
\begin{equation}
\sigma^{AB}=\kappa\frac{\sigma^A\sigma^B}{\sigma_\mathrm{eff}},\label{sec:produ:eq:dps}
\end{equation}
where $\kappa$ is a symmetry factor with $\kappa=1$ if $A\ne
B$ and the effective cross-section $\sigma_\mathrm{eff}$ is assumed to be  universal.
Heavy-quark fragmentations  in MPIs  are implied to be identical to that in inclusive production defined in Eq.~\ref{sec:produ:eq:dps}.
In particular, the kinematics of $A$ and $B$ is uncorrelated and each of them is similar to that in single particle inclusive production. Besides MPI, the single parton scattering (SPS) is also able to generate associated production, but in SPS the final-state kinematics is correlated. This difference between DPS and SPS is used to identify DPS. Studies of DPS include measurements of the $\sigma_\mathrm{eff}$ parameter and tests of its universality for different states, and investigations of kinematic correlations between $A$ and $B$. One example of correlation variables is the relative azimuthal angle  $\Delta\phi$ between $A$ and $B$. For DPS production $\Delta\phi$ distribution is approximately flat, while in SPS events a concentration at $\Delta\phi\sim0$ or $\pi$ is expected.

Measurements of associated production in $pp$ collisions are made at LHCb for two open charm hadrons~\cite{LHCb-PAPER-2012-003}, heavy quarkonia plus open charm~\cite{LHCb-PAPER-2015-046,LHCb-PAPER-2012-003}, and double ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace$ mesons~\cite{LHCb-PAPER-2016-057}.
To subtract the SPS contribution from data, one usually relies on theoretical inputs for cross-sections of SPS or fits to data using templates of correlation variables built for both SPS and DPS production.
Note that theoretical uncertainties are still much larger than
experimental ones for these measurements.

For some associated production, SPS is estimated or assumed to be negligible, resulting in a sample of approximately pure DPS events.
In this case the $\sigma_\mathrm{eff}$ parameter measured using Eq.~\ref{sec:produ:eq:dps} is around $15\aunit{mb}\xspace$ for ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace D$ and $\PUpsilon D$ production, independent of the $D$ species and
collision centre-of-mass energies~\cite{LHCb-PAPER-2012-003,LHCb-PAPER-2015-046}.
The results are similar to those extracted using multi-jet production at Tevatron~\cite{CDF:1997yfa}.
However, the values obtained using same-sign $DD$ pairs are around $20\aunit{mb}\xspace$,
and that for ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace\jpsi$ pairs is $7\aunit{mb}\xspace$~\cite{LHCb-PAPER-2016-057}.
The former are consistently higher than the value of $15\aunit{mb}\xspace$~\cite{LHCb-PAPER-2012-003}, while the latter is significantly lower.
Higher values of $\sigma_\mathrm{eff}$ for ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace\jpsi$ production are obtained if  a fraction of SPS is subtracted. The SPS fraction is estimated to be between $20\%$ and $40\%$ depending on the choice of control
variables and input templates for SPS and DPS distributions~\cite{LHCb-PAPER-2016-057}.
The smaller $\sigma_\mathrm{eff}$ measurement for ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace\jpsi$ pairs confirms previous
observations by D0, CMS and ATLAS experiments for quarkonium pairs~\cite{ATLAS:2016ydt,D0:2015dyx,Lansberg:2014swa}.  For correlation variables, as shown on the left of Fig.~\ref{sec:prod:fig:dps}, the $\Delta\phi$ distributions of ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace D$ events are reasonable flat~\cite{LHCb-PAPER-2012-003}, consistent with the DPS production in which ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace$ and $D$ kinematics is uncorrelated. This observation is a sign
of dominant or pure DPS contribution for ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace D$ samples. For same-sign $DD$ production, the correlation variables also favour DPS dominance~\cite{LHCb-PAPER-2012-003}.

The $\ensuremath{p_{\mathrm{T}}}\xspace$ distribution of each hadron in the pair production is also studied, which is expected to be similar to that in single inclusive production for DPS events. For $\PUpsilon D$ samples, both the $\PUpsilon$ and the $D$ meson have a  $\ensuremath{p_{\mathrm{T}}}\xspace$ distribution similar  to that in single
inclusive production~\cite{LHCb-PAPER-2015-046}. The same conclusion holds for $D$ \ensuremath{p_{\mathrm{T}}}\xspace distribution in ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace D$ events. However, the $\ensuremath{p_{\mathrm{T}}}\xspace$ of ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace$  mesons in ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace D$ events is
significantly harder than that in inclusive production, indicated by the right plot of Fig.~\ref{sec:prod:fig:dps}. For same-sign $DD$ production, the $D$ $\ensuremath{p_{\mathrm{T}}}\xspace$ distributions are also significantly harder than
those in single inclusive $D$ production, but are similar to those in opposite-sign $D{\ensuremath{\Dbar{}^0}}\xspace$ samples~\cite{LHCb-PAPER-2012-003}. This
result is inconsistent with the estimation  using correlation variables, which hints at a dominant DPS (SPS) contribution in the same-sign (opposite-sign) pair production.  Note that for associated
production of $DD$ pairs, the $D$ $\ensuremath{p_{\mathrm{T}}}\xspace$ distribution is similar for different $D$ species, indicating that charm hadron fragmentations are not modified, so that the unexpected $\ensuremath{p_{\mathrm{T}}}\xspace$ distribution is not due to the fragmentation process. A detailed theoretical calculation on ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace D$ production was performed recently to understand the problem, but a solid conclusion is not available yet~\cite{Shao:2020kgj}.

\begin{figure}[!tb]
\begin{center}
    \includegraphics[width=0.4\textwidth]{figs/production/DoubleDpp.pdf}
    \includegraphics[width=0.4\textwidth]{figs/production/DoubleDJpsipp.pdf}
    \caption{(Left) Distribution of relative azimuthal angle, $\Delta\phi$, between ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace$ and $D$ mesons in pair production, showing flat behaviour; (right) the $\ensuremath{p_{\mathrm{T}}}\xspace$ distribution of ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace$ mesons in ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace D$ pair production, compared with that in inclusive production (shown in black open circle)~\cite{LHCb-PAPER-2011-003}. Figures are taken from Ref.~\cite{LHCb-PAPER-2012-003}. }
\label{sec:prod:fig:dps}
\end{center}
\end{figure}

\subsubsection{Production in $p$Pb collisions}
Charm pair production is also studied in $p$Pb collisions of a center-of-mass energy per nucleon pair $\ensuremath{\protect\sqrt{s_{\scriptscriptstyle\text{NN}}}}\xspace=8.16\aunit{Te\kern -0.1em V}\xspace$~\cite{LHCb-PAPER-2020-010}. The DPS cross-section in $p$Pb collisions is expected to
scale with three times of the Pb mass number ($A_\mathrm{Pb}=208$) compared with that in $pp$ data at the same $\ensuremath{\protect\sqrt{s_{\scriptscriptstyle\text{NN}}}}\xspace$, rather than scale with a simple factor of $A_\mathrm{Pb}$, when nuclear matter effects are not considered~\cite{Baranov:2011ch}. The $A$-scaling is relevant for
SPS production in the absence of nuclear matter effects~\cite{Baranov:2011ch}, so that DPS production has a factor of three enhancement compared with SPS. The  cross-section ratio between same-sign $DD$ and
opposite-sign $D{\ensuremath{\offsetoverline{\PD}}}\xspace$ signals is measured to be around three times of that in $pp$ collisions~\cite{LHCb-PAPER-2020-010}. The result is in favour of the expected factor-three enhancement. The parameter $\sigma_\mathrm{eff}$ is measured with ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace D$ and same-sign $DD$ production as shown on the left of Fig.~\ref{sec:prod:fig:heavyIons} for the positive rapidity region, which corresponds to the Pb beam direction (high Bjorken-$x$ of Pb nucleus), and the negative rapidity region, which corresponds to the $p$ beam direction (low Bjorken-$x$ of Pb nucleus). The measurements show that, similar to the results in $pp$ data, the $\sigma_\mathrm{eff}$ parameter for ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace D$ production
is about 30\% smaller than that for
same-sign $DD$ production. Besides, the results in negative rapidity
hint at smaller values than those in positive rapidity for both ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace D$ and $DD$ pairs. It may be
a sign of universality violation of $\sigma_\mathrm{eff}$. This will be explored with better precision in Run 3 heavy-ion collision data, where ten times more
luminosity is expected to be collected.

\begin{figure}[!tb]
\begin{center}
    \includegraphics[width=0.49\textwidth]{figs/production/DPSpA.pdf}
    \includegraphics[width=0.4\textwidth]{figs/production/UpsilonpA.pdf}
    \caption{(Left) The parameter $\sigma_\mathrm{eff}$ measured in $p$Pb collisions at $\ensuremath{\protect\sqrt{s_{\scriptscriptstyle\text{NN}}}}\xspace=8.16\aunit{Te\kern -0.1em V}\xspace$ using (green) ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace{\ensuremath{\D^0}}\xspace$ and (brown) ${\ensuremath{\D^0}}\xspace\Dz$ production for (positive rapidity) $p$ beam and (negative rapidity) Pb beam direction, taken from~\cite{LHCb-PAPER-2020-010}. The prediction shown in shaded area is from Ref.~\cite{dEnterria:2012jam}. (Right) The ratio of cross-sections between $\PUpsilon(3S)$ and $\PUpsilon(1S)$ over that in $pp$ data at $\ensuremath{\protect\sqrt{s_{\scriptscriptstyle\text{NN}}}}\xspace=8\aunit{Te\kern -0.1em V}\xspace$, taken from Ref.~\cite{LHCb-PAPER-2018-035}, is compared with the comovers model~\cite{Ferreiro:2018wbd}. The rapidities are defined in the rest frame of two colliding nucleons with respect to proton beam direction.}
\label{sec:prod:fig:heavyIons}
\end{center}
\end{figure}

Besides the enhancement of DPS production, heavy nuclear collisions have many more new phenomena compared with $pp$ collisions, collectively called nuclear matter
effects. Presence of nuclear matter effects in $p$Pb collisions could modify the PDF, or reduce parton energies or  dissociate heavy quarkonia, which can be probed using heavy-flavour production in $p$Pb data compared with the $pp$ cross-section scaled by $A_\mathrm{Pb}$~\cite{Gavin:1991qk,Armesto:2006ph,Arleo:2012rs}. Measurements of ${\ensuremath{\D^0}}\xspace$ and
${\ensuremath{\Bu}}\xspace$ production in $p$Pb data suggest heavy-quark production in the $p$ beam direction is significantly
suppressed compared with the $A_\mathrm{Pb}$ scaling, by about 30\%, while the production in Pb beam direction
approximately scales with $A_\mathrm{Pb}$~\cite{LHCb-PAPER-2017-015,LHCb-PAPER-2018-048}. The results are
consistent with modifications of the gluon PDF in a Pb nucleus compared with that in a free nucleon. The LHCb measurements are found to be able to reduce the
gluon PDF uncertainties in the Pb nucleus by about a factor of three compared with the commonly used nuclear PDF sets~\cite{Kusina:2017gkz}. A new precise measurement of ${\ensuremath{\D^0}}\xspace$ production in $p$Pb shows that
the magnitude of the ${\ensuremath{\D^0}}\xspace$ suppression in $p$ beam
direction reduces significantly at high $\ensuremath{p_{\mathrm{T}}}\xspace({\ensuremath{\D^0}}\xspace)$ and the suppression
seems to vanish at $\ensuremath{p_{\mathrm{T}}}\xspace({\ensuremath{\D^0}}\xspace)>8\ensuremath{\aunit{Ge\kern -0.1em V\!/}c}\xspace$~\cite{LHCb-CONF-2019-004}.
However according to predictions using the nuclear PDF, the suppression is
about 30\%, independent of $\ensuremath{p_{\mathrm{T}}}\xspace({\ensuremath{\D^0}}\xspace)$. The observed suppression trend for ${\ensuremath{\D^0}}\xspace$ may be caused by the parton energy loss effect which alters heavy-flavour kinematic distribution, and its impact is reduced at high $\ensuremath{p_{\mathrm{T}}}\xspace$~\cite{Arleo:2021bpv}.

Measurement of ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace$ production in
$p$Pb data shows a similar trend of suppression compared with open heavy-flavour hadrons~\cite{LHCb-PAPER-2013-052,LHCb-PAPER-2017-014}, suggesting that they suffer from a common influence by nuclear matter effects. However, the result for ${\ensuremath{\Ppsi{(2S)}}}\xspace$ in $p$Pb data
suggests a stronger suppression compared with ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace$, in particularly in the Pb beam direction~\cite{LHCb-PAPER-2015-058}. Similarly,
the $\PUpsilon(3S)$ meson is measured to be more suppressed compared with $\PUpsilon(1S)$~\cite{LHCb-PAPER-2018-035}, as shown on the right of Fig.~\ref{sec:prod:fig:heavyIons}. The stronger suppression for excited quarkonium can't be explained using the nuclear PDF
modification or the parton energy loss effect. A model introducing final state interactions between a heavy quarkonium and comoving particles
is able to explain data~\cite{Ferreiro:2018wbd}. The comovers effect is stronger in events of higher occupancy and for particles with
larger sizes, such that it is more
pronounced in Pb beam direction and for excited states. The comovers mechanism also exists in $pp$ collisions, and is probed using heavy quarkonium production.
The cross-section ratio between prompt $\chi_{c1}(3872)$ and ${\ensuremath{\Ppsi{(2S)}}}\xspace$ mesons is measured to decrease with the increase of the number of reconstructed tracks in
the vertex detector~\cite{LHCb-PAPER-2020-023}. It suggests that the $\chi_{c1}(3872)$ state has a larger size or a smaller binding energy compared  with the ${\ensuremath{\Ppsi{(2S)}}}\xspace$ meson, and is consistent with a component of ${\ensuremath{\D^{*0}}}\xspace{\ensuremath{\Dbar{}^0}}\xspace+{\ensuremath{\Dbar{}^{*0}}}\xspace{\ensuremath{\D^0}}\xspace$ hadron molecule in the $\chi_{c1}(3872)$ wave function~\cite{Braaten:2020iqw,Esposito:2020ywk}.
The same measurement in $p$Pb data feasible in LHC Run 3 period would be very important to confirm the results.

In the near future, new data will provide enough statistics for associated production of triple heavy flavours and heavy quarkonium pairs beyond ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace\jpsi$, which help to further understand
MPI and heavy quarkonium production mechanism~\cite{Chapon:2020heu}. Concerning studies of nuclear matter effects using heavy-flavour production, a rich program is foreseen in Run 3, including measurements probing the DPS enhancement, the modification of nuclear PDF and heavy-quark
fragmentation and the heavy quarkonium dissociation mechanism.

\subsection{Spectroscopy}
The strong interaction confines quarks (and/or gluons) to form various colour-singlet hadrons that are accessible
experimentally. This confinement phenomenon is nonperturbative and is not understood yet from the current QCD theory. In
analogy to photon spectroscopy in atomic physics, hadron spectroscopy provides a way to understand dynamics of
QCD at low energy.  Hadrons composed of a quark and an antiquark are called mesons, those of three quarks are called baryons, and those
composed of more than three quarks are usually referred to as exotic hadrons. Existence of exotic hadrons have been predicted since the birth of the quark model and their properties are reexamined by refined theoretical approaches in the past decades~\cite{Gell-Mann:1964ewy,Zweig:570209,Zhu:2003ba,Sun:2011uh,Liu:2009ei,Liu:2008fh,Zhu:2007wz,Chen:2019asm,Liu:2007bf}. The past years witnessed the observations of a plethora of new conventional and
exotic hadrons containing heavy quarks and LHCb is one of the leading players in the field.
Figure~\ref{sec:prod:fig:hadrons} displays the 62 new hadrons that are discovered by the LHC experiments, and most of
them by LHCb.

\begin{figure}[!htpb]
\begin{center}
    \includegraphics[width=0.9\textwidth]{figs/production/LHChadrons_s.png}
    \caption{New hadrons that are discovered by the LHC experiments, taken from
    \cite{lhcHadron}.}
\label{sec:prod:fig:hadrons}
\end{center}
\end{figure}

\subsubsection{Conventional hadrons}

\paragraph{Classification of heavy baryons} Following the heavy-quark symmetry, baryons with a heavy quark $Q$ are organised into multiplets according to quantum configurations of the two light quarks~\cite{Liu:2007fg,Yu:2021zvl,Wang:2021bmz,Matsuki:2017tqf,Wang:2017kfr,Mao:2017wbz,Cheng:2017ove,Yang:2020zrh,Chen:2018orb,Wang:2017goq,Wei:2016jyk,Lu:2016gev,He:2021xrh,Zhang:2008pm,Zhang:2008iz,Wei:2010zza,Guo:2008ns}.
The total wave function including flavour ($F$), spin ($s_{qq'}$) and orbital angular momentum ($l_{qq'}$) must be symmetric for the two light quarks $qq'$ to form an antisymmetric state
together with their antisymmetric colour configuration. Baryons with $l_{qq'}=s_{qq'}=0$ (antisymmetric in spin space) have a spin-parity of  $J^P=1/2^{+}$, and are grouped into a multiplet of three flavour-antisymmetric states for each heavy quark $Q$. While baryons with $s_{qq'}=1, l_{qq'}=0$
(being symmetric)  have
$J^P=1/2^+$ or $J^P=3/2^+$,  and form two multiplets of six flavour-symmetric states for the two $J^P$. These three different multiplets are shown in the bottom row of Fig.~\ref{sec:prod:fig:HbB}. Orbital and radial excitation can happen inside the two light quarks
($\rho$-mode) or between $Q$ and the $qq'$ system ($\lambda$-mode). The parity of a baryon is determined to be  $P=(-1)^{l_\rho+l_\lambda}$, where $l_\rho=l_{qq'}$, and $l_\lambda$ is the orbital angular momentum between $Q$ and $qq'$.
Beauty and charm baryons with $l_\lambda=l_{qq'}=s_{qq'}=0$ decay weakly and have been well established, but a chart of their excited states are far from being
complete. As a matter of fact only a few low lying states are observed, in particularly for beauty baryons as can be seen in Fig.~\ref{sec:prod:fig:HbB}.
Up to date no sign of $\rho$-mode states have been identified experimentally, probably because
they are too wide (hundreds of $\kern -0.125em\aunit{Me\kern -0.1em V}\xspace$) to be resolved from underlying background.

\paragraph{Charm baryons} In the invariant mass spectrum of ${\ensuremath{\PXi}}\xspace^+_c{\ensuremath{\kaon^-}}\xspace$ hadrons shown on the left of Fig.~\ref{sec:prod:fig:charmBaryons}, LHCb observed five states whose quark contents are considered to be $css$:  ${\ensuremath{\POmega}}\xspace_c(3000)^0$, ${\ensuremath{\POmega}}\xspace_c(3050)^0$,
${\ensuremath{\POmega}}\xspace_c(3066)^0$, ${\ensuremath{\POmega}}\xspace_c(3090)^0$ and ${\ensuremath{\POmega}}\xspace_c(3119)^0$~\cite{LHCb-PAPER-2017-002}.
All these states have narrow widths,  below $10\aunit{Me\kern -0.1em V}\xspace$, and their mass differences are only tens of $\kern -0.125em\aunit{Me\kern -0.1em V}\xspace$.
The first four states are
confirmed by Belle in ${\ensuremath{\Pe^+\Pe^-}}\xspace$ collisions~\cite{Belle:2017ext} and by LHCb in the exclusive ${\ensuremath{\POmega}}\xspace_b^{-}\ensuremath{\rightarrow}\xspace {\ensuremath{\PXi}}\xspace^+_c{\ensuremath{\kaon^-}}\xspace{\ensuremath{\pion^-}}\xspace$ decay~\cite{LHCb-PAPER-2021-012}. The spin assignments of the first four states favour $1/2, 3/2, 3/2, 5/2$, consistent with the expectations for $P$-wave $\lambda$ excitation~\cite{LHCb-PAPER-2021-012}. A determination of their parities will help to make firm conclusions. According to phenomenological models~\cite{Karliner:2017kfm,Wang:2017hej,Xu:2019kkt,Chen:2017gnu,Wang:2017vnc,Chen:2017sci,Yang:2017rpg}, one of the five $1P$ states with $l_\rho=0,l_\lambda=1$ in the mass region of observed states is missing, and the ${\ensuremath{\POmega}}\xspace_c(3119)^0$ state may be a $2S$ or $D$-wave baryon. In the high mass region of LHCb data, a hint of a wide state ${\ensuremath{\POmega}}\xspace_c(3188)^0$ is present,  to be confirmed in future analysis with additional statistics. It is noted some of these states are considered to be exotic states of quark constituents $cssu{\ensuremath{\overline \uquark}}\xspace$ rather than conventional $css$ baryons~\cite{Wang:2021cku,Wang:2018alb,Chen:2017xat,Wang:2017smo,Wang:2017xam}.

Similarly, excited ${\ensuremath{\PXi}}\xspace^0_c$ states are searched for by LHCb~\cite{LHCb-PAPER-2020-004} in the ${\ensuremath{\Lz^+_\cquark}}\xspace{\ensuremath{\kaon^-}}\xspace$ invariant mass spectrum shown on the right of Fig.~\ref{sec:prod:fig:charmBaryons}. A
new state ${\ensuremath{\Xires_\cquark}}\xspace(2965)^0$ is observed, and the state ${\ensuremath{\Xires_\cquark}}\xspace(2930)^0$ claimed by the Belle experiment~\cite{Belle:2017jrt} now splits into two structures, ${\ensuremath{\Xires_\cquark}}\xspace(2923)^0$ and
${\ensuremath{\Xires_\cquark}}\xspace(2939)^0$. The widths of these three states are determined to be around $10\aunit{Me\kern -0.1em V}\xspace$. These states and previously known ${\ensuremath{\Xires_\cquark}}\xspace(2790)^0$ and
${\ensuremath{\Xires_\cquark}}\xspace(2815)^0$ lie in the mass region $1P$  excitation~\cite{Ebert:2007nw,Roberts:2007ni,Migura:2006ep,Yang:2021lce,Chen:2021eyk}. There are in total seven $1P$ ${\ensuremath{\PXi}}\xspace^0_c$ states of the $\lambda$ excitation. At least two of these $1P$ states are still missing;
Besides, a complete and solid matching of these observed states to predicted spectrum is not resolved yet~\cite{Nieves:2019jhp}.

\begin{figure}[!tb]
\begin{center}
    \includegraphics[width=0.9\textwidth]{figs/production/HbPart.pdf}
    \caption{Experimentally identified bottom baryons, grouped according to spin and flavour symmetry of the light quarks as well as the  excitation between the light quark system and $b$ quark. The assignments of quantum numbers are based on measured properties and relevant predictions. Note that for recently discovered states their $J^P$ need confirmation and the  assignments of excited ${\ensuremath{\POmega}}\xspace_b$ are not certain, though they likely have $J^P=1^-$. In the figure, 1$P$ excitation of flavour symmetric bottom baryons are gathered into a single multiplet (shown in purple). Only $\lambda$ mode radial and orbital angular momentum excitation are considered. Missing states are marked as question marks in the figure. The notation $s_{qq'}$ indicates the total spin of the light quark system.}
\label{sec:prod:fig:HbB}
\end{center}
\end{figure}

\begin{figure}[!tb]
\begin{center}
    \includegraphics[width=0.41\textwidth]{figs/production/Omegac.pdf}
    \includegraphics[width=0.52\textwidth]{figs/production/Xicstar.pdf}
    \caption{(Left) Invariant mass distribution of ${\ensuremath{\PXi}}\xspace^+_c{\ensuremath{\kaon^-}}\xspace$, overlaid with the fit including five narrow ${\ensuremath{\POmega}}\xspace_c$ states, taken from Ref.~\cite{LHCb-PAPER-2017-002}. (Right) Invariant mass distribution of ${\ensuremath{\Lz^+_\cquark}}\xspace{\ensuremath{\kaon^-}}\xspace$, overlaid with the fit including three narrow ${\ensuremath{\Xires_\cquark}}\xspace$ states (red, brown and purple dashed lines) and partially reconstructed components, taken from Ref.~\cite{LHCb-PAPER-2020-004}.}
\label{sec:prod:fig:charmBaryons}
\end{center}
\end{figure}

\paragraph{Bottom baryons}
A summary of all bottom baryons with clear experiment evidences are shown in Fig.~\ref{sec:prod:fig:HbB}, organised in multiplets of $qq'$ flavour symmetry and $\lambda$-mode excitation when possible. The spin-parity quantum numbers for most of these states are not measured, so we rely on theoretical calculations as a guidance to make the classification. Actually, there is not always a consensus on the $J^P$ of each state, in particularly for ${\ensuremath{\POmega}}\xspace_b$ baryons~\cite{Xu:2020ofp,Karliner:2020fqe,Xiao:2020gjo,Yang:2020zjl,Xiao:2020oif,Wang:2020pri,Liang:2020dxr,Chen:2020mpy,Liang:2020hbo,Lu:2019rtg,Chen:2018vuc,Wang:2020vwl,Mao:2015gya,Lu:2014ina,Yang:2018lzg,Huang:2018bed}.

In total, five states have been reported in the ${\ensuremath{\Lz^0_\bquark}}\xspace{\ensuremath{\pion^+}}\xspace{\ensuremath{\pion^-}}\xspace$ mass spectrum: ${\ensuremath{\Lz_\bquark}}\xspace(5912)^0$ and $
{\ensuremath{\Lz_\bquark}}\xspace(5920)^0$ with widths below $1 \aunit{Me\kern -0.1em V}\xspace$~\cite{LHCb-PAPER-2012-012}, ${\ensuremath{\Lz_\bquark}}\xspace(6146)^0$ and ${\ensuremath{\Lz_\bquark}}\xspace(6152)^0$ with
widths of about $2\aunit{Me\kern -0.1em V}\xspace$~\cite{LHCb-PAPER-2019-025}, and ${\ensuremath{\Lz_\bquark}}\xspace(6072)^0$ with a width around
$70\aunit{Me\kern -0.1em V}\xspace$~\cite{LHCb-PAPER-2019-045}. Their masses match two $1P$, two $1D$ and $2S$ ${\ensuremath{\Lz^0_\bquark}}\xspace$ $\lambda$-mode excitation respectively, though other assignments are also discussed~\cite{Chen:2019ywy,Yang:2019cvw,Wang:2019uaj,Wang:2018fjm,Mao:2020jln}. It is useful to note that intermediate
${\ensuremath{\PSigma}}\xspace^{(*)\pm}_b(\ensuremath{\rightarrow}\xspace{\ensuremath{\Lz^0_\bquark}}\xspace\pi^{\pm})$ states are found to be present in the ${\ensuremath{\Lz_\bquark}}\xspace^{*(*)0}\ensuremath{\rightarrow}\xspace{\ensuremath{\Lz^0_\bquark}}\xspace{\ensuremath{\pion^+}}\xspace{\ensuremath{\pion^-}}\xspace$ decays.

The ground
${\ensuremath{\PSigma}}\xspace^{\pm}_b$ and ${\ensuremath{\PSigma}}\xspace^{*\pm}_b$ states were first detected in the ${\ensuremath{\Lz^0_\bquark}}\xspace\pi^{\pm}$ mass spectrum  by CDF~\cite{CDF:2007oeq}.
In the same final state, two new ones, ${\ensuremath{\Sigmares_\bquark}}\xspace(6097)^+$ and ${\ensuremath{\Sigmares_\bquark}}\xspace(6097)^-$, are observed by LHCb~\cite{LHCb-PAPER-2018-032}, whose widths are about $30\aunit{Me\kern -0.1em V}\xspace$. These two new states
belong to the $P$-wave family, and many more of them are still missing, like charm baryons.

In analogy, excited ${\ensuremath{\Xires_\bquark}}\xspace$ states are searched for by the
LHC experiments in the ${\ensuremath{\PXi}}\xspace^{0}_{\ensuremath{\Pb}}\xspace\pi^{-}/{\ensuremath{\PXi}}\xspace^{-}_{\ensuremath{\Pb}}\xspace\pi^{+}$ spectra. New states close to ${\ensuremath{\Xires_\bquark}}\xspace\pi$ mass thresholds are observed, which include the low lying ${\ensuremath{\PXi}}\xspace^{*0}_{\ensuremath{\Pb}}\xspace$ baryon discovered
by CMS~\cite{CMS:2012frl} and ${\ensuremath{\PXi}}\xspace^{'-}_{\ensuremath{\Pb}}\xspace$, ${\ensuremath{\PXi}}\xspace^{*-}_{\ensuremath{\Pb}}\xspace$ states discovered by LHCb~\cite{LHCb-PAPER-2014-061}.
The states ${\ensuremath{\PXi}}\xspace^{'-}_{\ensuremath{\Pb}}\xspace$ and ${\ensuremath{\PSigma}}\xspace^{'\pm}_b$ belong to the  flavour symmetric multiplet with $s_{qq'}=1, J^P=1/2^+$, while ${\ensuremath{\PXi}}\xspace^{*0}_{\ensuremath{\Pb}}\xspace$, ${\ensuremath{\PXi}}\xspace^{*-}_{\ensuremath{\Pb}}\xspace$ and ${\ensuremath{\PSigma}}\xspace^{*\pm}_b$ belong to the flavour symmetric  multiplet with $s_{qq'}=1, J^P=3/2^+$. Going to higher mass region, a state ${\ensuremath{\Xires_\bquark}}\xspace(6227)^-$, with a width around $20 \aunit{Me\kern -0.1em V}\xspace$,  is found in both ${\ensuremath{\Xires^0_\bquark}}\xspace{\ensuremath{\pion^-}}\xspace$ and
${\ensuremath{\Lz^0_\bquark}}\xspace{\ensuremath{\kaon^-}}\xspace$ final states~\cite{LHCb-PAPER-2018-013}, and its partner ${\ensuremath{\Xires_\bquark}}\xspace(6227)^0$  is found
in the ${\ensuremath{\Xires^-_\bquark}}\xspace{\ensuremath{\pion^+}}\xspace$ mass spectrum~\cite{LHCb-PAPER-2020-032}. They can be matched to $P$-wave
states or a mixture of several $P$-wave states with masses close to $6227\ensuremath{\aunit{Me\kern -0.1em V\!/}c^2}\xspace$. Very recently, two new states ${\ensuremath{\Xires_\bquark}}\xspace(6327)^0$ and
${\ensuremath{\Xires_\bquark}}\xspace(6333)^0$, with widths below $2\aunit{Me\kern -0.1em V}\xspace$, are found in the ${\ensuremath{\Lz^0_\bquark}}\xspace{\ensuremath{\kaon^-}}\xspace{\ensuremath{\pion^+}}\xspace$ mass spectrum~\cite{LHCb-PAPER-2021-025}, consistent with the
$1D$ excitation of the ${\ensuremath{\Xires^0_\bquark}}\xspace$ baryon. These two states may also be present in the ${\ensuremath{\Xires^0_\bquark}}\xspace{\ensuremath{\pion^+}}\xspace{\ensuremath{\pion^-}}\xspace$ sample as well, demanding a future investigation of this decay mode by LHCb. In fact, in the ${\ensuremath{\Xires^-_\bquark}}\xspace{\ensuremath{\pion^+}}\xspace{\ensuremath{\pion^-}}\xspace$ spectrum, a ${\ensuremath{\Xires_\bquark}}\xspace(6100)^-$ state is observed by CMS~\cite{CMS:2021rvl}, consistent with the $1P$ excitation of the flavour antisymmetric ${\ensuremath{\Xires^-_\bquark}}\xspace$ state and has $J^P=3/2^-$. Apparently, the other $1P$ state with $J^P=1/2^-$ and a mass around $6100\ensuremath{\aunit{Me\kern -0.1em V\!/}c^2}\xspace$ is missing. No states with a higher mass, for example flavour partners of ${\ensuremath{\Xires_\bquark}}\xspace(6327)^0$ and ${\ensuremath{\Xires_\bquark}}\xspace(6333)^0$,  are reported in CMS data, which may be explained by lower production rate for these states.

Excited ${\ensuremath{\Omegares^-_\bquark}}\xspace$ states are searched for in the ${\ensuremath{\Xires^0_\bquark}}\xspace{\ensuremath{\kaon^-}}\xspace$ mass spectrum~\cite{LHCb-PAPER-2019-042}. Four narrow (width $<5\aunit{Me\kern -0.1em V}\xspace$) peaking
structures are identified with two of them having high significance, named as ${\ensuremath{\POmega}}\xspace_b(6340)^-$ and
${\ensuremath{\POmega}}\xspace_b(6350)^-$ respectively. These states lie in the mass region of  $P$-wave excitation. More statistics in Run 3 will  allow for a further investigation of these states.

\begin{figure}[!tb]
\begin{center}
    \includegraphics[width=0.5\textwidth]{figs/production/Bc2S.pdf}
    \includegraphics[width=0.43\textwidth]{figs/production/XiccObs.pdf}
    \caption{(Left) Invariant mass distribution of $\Delta M\equiv m_{{\ensuremath{\B_\cquark^+}}\xspace{\ensuremath{\pion^+}}\xspace{\ensuremath{\pion^-}}\xspace}-m_{{\ensuremath{\B_\cquark^+}}\xspace}$, overlaid with the fit including two narrow structures, corresponding ${\ensuremath{\PB}}\xspace_c^{*}(2S)^+$ and ${\ensuremath{\PB}}\xspace_c(2S)^+$ respectively, taken from Ref.~\cite{LHCb-PAPER-2019-007}. (Right) Distribution of the reconstructed invariant mass $m_\mathrm{cand}({\ensuremath{\Xires^{++}_{\cquark\cquark}}}\xspace)\equiv m({\ensuremath{\Lz^+_\cquark}}\xspace{\ensuremath{\kaon^-}}\xspace{\ensuremath{\pion^+}}\xspace) -m_\mathrm{cand}({\ensuremath{\Lz^+_\cquark}}\xspace)+m_\mathrm{PDG}({\ensuremath{\Lz^+_\cquark}}\xspace)$, overlaid with the fit projections, taken from Ref.~\cite{LHCb-PAPER-2017-018}.}
\label{sec:prod:fig:BcXicc}
\end{center}
\end{figure}

\paragraph{Doubly heavy hadrons} LHCb has observed a few new hadrons with two heavy quarks. In the final states of ${\ensuremath{\D^0}}\xspace{\ensuremath{\Dbar{}^0}}\xspace$ and
${\ensuremath{\D^+}}\xspace{\ensuremath{\D^-}}\xspace$ mesons produced promptly, a new particle $X(3842)$,  with a width of about $3\aunit{Me\kern -0.1em V}\xspace$, is discovered~\cite{LHCb-PAPER-2019-005}. It is consistent with the spin-3 conventional charmonium  $\psi_3({}^3D_3)$ with $J^{PC}=3^{--}$~\cite{Barnes:2005pb}.

As shown on the left of Fig.~\ref{sec:prod:fig:BcXicc}, two narrow structures are detected in the ${\ensuremath{\B_\cquark^+}}\xspace{\ensuremath{\pion^+}}\xspace{\ensuremath{\pion^-}}\xspace$ invariant mass spectrum of LHCb data~\cite{LHCb-PAPER-2019-007}. The left one corresponds to the ${\ensuremath{\PB}}\xspace_c^{*}(2^3S_1)^+\ensuremath{\rightarrow}\xspace {\ensuremath{\PB}}\xspace_c^{*}(1^3S_1)^+{\ensuremath{\pion^+}}\xspace{\ensuremath{\pion^-}}\xspace$ decay with the
photon in  decay of ${\ensuremath{\PB}}\xspace_c^{*}(1^3S_1)^+\ensuremath{\rightarrow}\xspace{\ensuremath{\B_\cquark^+}}\xspace\gamma$ not detected. The peak on the right is consistent with the ${\ensuremath{\PB}}\xspace_c(2^1S_0)^+\ensuremath{\rightarrow}\xspace
{\ensuremath{\B_\cquark^+}}\xspace{\ensuremath{\pion^+}}\xspace{\ensuremath{\pion^-}}\xspace$ decay. Almost in parallel, these two states are also independently observed by CMS~\cite{CMS:2019uhm}. As the ${\ensuremath{\B_\cquark^+}}\xspace$ meson is  composed of two heavy quarks, its excitation mass spectroscopy can be calculated using models similar to those applied to heavy quarkonia~\cite{Gupta:1995ps,Chen:1992fq,Ding:2021dwh,Chen:2020ecu,Chang:2019wpt,Chang:2019eob} despite their hadroproduction mechanisms are very different~\cite{Chang:2003cq,Chang:2005hq}.

LHCb opens a new era in studies of doubly heavy baryons by observing
the ${\ensuremath{\Xires^{++}_{\cquark\cquark}}}\xspace$ baryon in the ${\ensuremath{\Xires^{++}_{\cquark\cquark}}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\Lz^+_\cquark}}\xspace{\ensuremath{\kaon^-}}\xspace{\ensuremath{\pion^+}}\xspace\pip$ mass spectrum~\cite{LHCb-PAPER-2017-018}, as shown on the right of Fig.~\ref{sec:prod:fig:BcXicc}. This discovery decay mode is predicted to have a relatively large branching fraction~\cite{Yu:2019lxw,Hu:2019bqj}. The ${\ensuremath{\Xires^{++}_{\cquark\cquark}}}\xspace$ state is later confirmed in
${\ensuremath{\Xires^{++}_{\cquark\cquark}}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\Xires^+_\cquark}}\xspace{\ensuremath{\pion^+}}\xspace$ decays~\cite{LHCb-PAPER-2018-026}. Its lifetime is measured to be about
$0.25\ensuremath{\aunit{ps}}\xspace$~\cite{LHCb-PAPER-2019-037} and its mass is  precisely determined to be
$3621.55\pm0.38\ensuremath{\aunit{Me\kern -0.1em V\!/}c^2}\xspace$~\cite{LHCb-PAPER-2019-037}. Its $SU(3)$ partners,  ${\ensuremath{\Xires^+_{\cquark\cquark}}}\xspace$ and ${\ensuremath{\Omegares^+_{\cquark\cquark}}}\xspace$,  are  also searched for by LHCb, but with no
sign of observation yet~\cite{LHCb-PAPER-2019-029,LHCb-PAPER-2021-011,LHCb-PAPER-2021-019}. Theoretically, an important question is to understand the mass spectrum of the doubly charm  baryons and related systems, which depends on the binding energy between the two heavy quarks~\cite{Tong:2021raz,Li:2020uok,Wang:2021rjk,Liu:2019vtx,Dias:2019klk,Yu:2019yfr,Li:2019ekr,Chen:2017sbg,Wang:2017qvg,Guo:2021yws,Qin:2021dqo,Gao:2020ogo,Weng:2018mmf}. The success discovery of ${\ensuremath{\Xires^{++}_{\cquark\cquark}}}\xspace$ has triggered wide theoretical work to understand the properties of baryons with more than a heavy quark~\cite{Han:2021gkl,Li:2021rfj,Han:2021azw,Wang:2020avt,Xiao:2017dly,Shi:2019hbf,Cheng:2019sxr,Zhang:2018llc,Jiang:2018oak,Zhao:2018zcb,Wang:2018utj,Cui:2017udv,Shi:2017dto,Meng:2017dni,Geng:2017mxn,Lu:2017meb,Xiao:2017udy,Li:2017pxa,Wang:2017azm,Wang:2017mqp,Yu:2017zst,Chen:2013upa,Wang:2017dcq}. In addtion, the ${\ensuremath{\Xires^{++}_{\cquark\cquark}}}\xspace$ baryon mass  is used to study the stability of tetraquark states with $QQ$ contents~\cite{Cheng:2020wxa,Lu:2020rog} with the assumption that $QQ$ form a heavy diquark~\cite{Shi:2020qde,An:2019idk,Meng:2018zbl,Zhou:2018bkn,Chen:2016ont}.

The beauty hadron decay is an ideal place to study excited charm hadrons, since
$b\ensuremath{\rightarrow}\xspace c$ is  the dominant transition of the $b$ quark. For example, in $B^0\ensuremath{\rightarrow}\xspace {\ensuremath{\D^-}}\xspace{\ensuremath{\D^+}}\xspace{\ensuremath{\kaon^+}}\xspace{\ensuremath{\pion^-}}\xspace$ decays, LHCb observes a new excited
$D_s^+$ meson, $D_s(2590)^+$~\cite{LHCb-PAPER-2020-034}. Its width is measured to about $100\aunit{Me\kern -0.1em V}\xspace$ and spin-parity determined to be $J^P=0^-$. This state is consistent
with the radial excited state $D_s(2^1S_0)^+$. The  mass and width of the new $D_s(2^1S_0)^+$ state may help to understand the excitation spectrum of ${\ensuremath{\D^+_\squark}}\xspace$ mesons which are found to be not fully consistent with the quark model predictions~\cite{Luo:2021dvj,Ni:2021pce,Xie:2021dwe,Yang:2021tvc,Wang:2021orp,Wang:2012wk}.

\subsubsection{Exotic hadrons}

Many theoretical efforts have been placed to understand how quarks are combined to form a multi-body system~\cite{Dong:2021rpi,Ding:2021igr,Chen:2021crg,Yang:2020fou}. At the same time, more and more new states are observed experimentally which cannot fit into the conventional hadron spectra~\cite{Olsen:2017bmm}. As a result, the study on exotic hadrons has been a hot topic for the past decade. New results from LHCb are discussed below.

\paragraph{Tetraquark states}  LHCb provides essential information
for the understanding of previously known tetraquark states. For example LHCb determined the  quantum number of the $X(3872)$ state first reported by Belle~\cite{Belle:2003nnu} to
be $J^{PC}=1^{++}$ through a full amplitude analysis~\cite{LHCb-PAPER-2015-015}. LHCb also precisely measured the mass of the $X(3872)$ state  (referred to as $\chi_{c1}(3872)$ in Ref.~\cite{PDG2020}) to be~\cite{LHCb-PAPER-2020-008,LHCb-PAPER-2020-009},  $m_{\chi_{c1}(3872)}-m_{{\ensuremath{\Ppsi{(2S)}}}\xspace}=185.49\pm0.06\pm0.03\ensuremath{\aunit{Me\kern -0.1em V\!/}c^2}\xspace$. Its Breit-Wigner (BW) width is determined by LHCb to be~\cite{LHCb-PAPER-2020-008,LHCb-PAPER-2020-009},
$\Gamma_{\chi_{c1}(3872)}^\mathrm{BW}=0.96^{+0.19}_{-0.18}\pm0.21\aunit{Me\kern -0.1em V}\xspace$. The $\chi_{c1}(3872)$ state is the mostly studied exotic candidate, and its exotic behaviours compared with conventional charmonia include a extremely  narrow width, isospin breaking decays and a mass close to the $D^*D$ threshold. However we are still not sure whether it is a compact tetraquark state, a ${\ensuremath{\D^{*0}}}\xspace{\ensuremath{\Dbar{}^0}}\xspace$ hadron molecule, a mixture of ${\ensuremath{\D^{*0}}}\xspace{\ensuremath{\Dbar{}^0}}\xspace$ molecule with a ${\ensuremath{\Pchi_{\cquark 1}}}\xspace(2P)$ charmonium component or just caused by kinematic rescattering effect~\cite{Zhang:2020mpi,Meng:2021kmi,He:2021ihy,Maiani:2004vq,Wang:2013vex,Chen:2010ze,Li:2004sta,Guo:2014iya,Guo:2013zbw,Meng:2007cx,Guo:2016bjq,Luo:2017eub,Kang:2016jxw,Yuan:2005dr,Cui:2006mp,Meng:2005er,Zhang:2009vs,Zhang:2009bv,Meng:2013gga,Ding:2009vj,Chen:2016otp,Li:2012cs,Liu:2019stu,Chen:2013pya,Dong:2021juy,Guo:2014taa,Guo:2019qcn,Zhao:2014gqa,Liu:2008qb,Ke:2012gm,Geng:2017hxc,Meng:2014ota,Chen:2017vai,Shen:2010ky,Zhang:2007xa,Wang:2012cp,Wang:2015rcz,Yang:2017prf,Liang:2021fzr,Sun:2017wgf,Yu:2011wb,Cao:2020dlz,Qiao:2013xca}.

The exotic candidate $Z_c(4430)^-$ was first observed by Belle in the ${\ensuremath{\Ppsi{(2S)}}}\xspace{\ensuremath{\pion^-}}\xspace$ mass spectrum in  ${\ensuremath{\B^0}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\Ppsi{(2S)}}}\xspace{\ensuremath{\kaon^+}}\xspace{\ensuremath{\pion^-}}\xspace$ decays~\cite{Belle:2007hrb}. An
amplitude analysis of the ${\ensuremath{\B^0}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\Ppsi{(2S)}}}\xspace{\ensuremath{\kaon^+}}\xspace{\ensuremath{\pion^-}}\xspace$ decay is performed at LHCb, confirming the existence of the $Z_c(4430)^-$ state and determining it to be consistent
with a Breit-Wigner resonance with $J^P=1^+$~\cite{LHCb-PAPER-2014-014}. The quark contents of $Z_c(4430)^-$, $c{\ensuremath{\overline \cquark}}\xspace{\ensuremath{\overline \uquark}}\xspace d$,  are the same as the $Z_c(3900)^+$ state observed by the BESIII experiment~\cite{BESIII:2013ris}. Being charged, they are definitely not consistent with conventional charmonia and many phenomenological calculations are performed to explain their internal structure and properties~\cite{Wu:2021ezz,Wang:2021aql,Zhang:2020rqr,Xiao:2019spy,He:2017mbh,Chen:2015ata,Chen:2015fsa,Liu:2014eka,Ma:2014zva,Wang:2014vha,Ma:2014zua,Ke:2013gia,Zhao:2014qva,Gong:2016hlt,Gui:2018rvv,Wang:2019tlw,He:2014nxa,Liu:2009wb}.

The ${\ensuremath{\Bu}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace\phi{\ensuremath{\kaon^+}}\xspace$ decay is a zoo of exotic hadrons. In 2009, a narrow state $X(4140)$ was reported by CDF in the
${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace\phi$ mass spectrum of the ${\ensuremath{\Bu}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace\phi{\ensuremath{\kaon^+}}\xspace$ decay~\cite{CDF:2009jgo,CDF:2011pep}, and is later confirmed by CMS~\cite{CMS:2013jru}. The quark
contents of the $X(4140)$ state is likely to be ${\ensuremath{\Pc}}\xspace{\ensuremath{\overline \cquark}}\xspace{\ensuremath{\Ps}}\xspace{\ensuremath{\overline \squark}}\xspace$, consistent with an exotic hadron~\cite{Hu:2010fg,Meng:2019fan}, even though excited conventional
charmonia may have the chance to decay into ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace\phi$ too. In the amplitude analysis by LHCb using Run 1 data, four exotic candidates $X(4140)$, $X(4274)$, $X(4500)$ and $X(4700)$ are observed~\cite{LHCb-PAPER-2016-018}. Currently, these $X$ states are considered to be either hadron molecules or compact tetraquark states or high-mass conventional charmonia in various calculations~\cite{Cao:2019wwt,Hao:2019fjg,Wang:2020dya,Chen:2016ncs,Lu:2016cwr,Wu:2016gas,Liu:2016onn,Wang:2016ujn,Wang:2014gwa,Zhang:2009st,Liu:2009iw,Chen:2016oma,Wang:2009wk,Wang:2017mrt,He:2013oma,Ferretti:2020civ,Deng:2019dbg,Wang:2021gml,Shi:2021jyr,Liu:2021xje,Ge:2021sdq,Yang:2021sue,Yang:2020nrt,Wang:2020rcx,Chen:2013wca,Jin:2020yjn}.
The LHCb analysis is updated recently with a sample that has six times more statistics, in which three more $X$ states are reported~\cite{LHCb-PAPER-2020-044}. In addition,  two $Z_{cs}^+$ structures are observed in the ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace{\ensuremath{\kaon^+}}\xspace$
mass spectrum. The mass spectra and fit projections are shown in Fig.~\ref{sec:prod:fig:JpsiPhiK} and the properties of these exotic candidates are summarised in Table~\ref{sec:prod:tab:B2JpsiPhiK}.  The $Z_{cs}^+$ states mark the first observation of  exotic hadrons
with an $s$ quark through beauty decays. It is noted that another $Z_{cs}^+$ state is reported  by BESIII in the final state of
$D_s^-D^{*0}+D_s^{*-}{\ensuremath{\D^0}}\xspace$ pairs~\cite{BESIII:2020qkh}, with a mass and width different from those observed by LHCb. There should be more states of $c{\ensuremath{\overline \cquark}}\xspace u{\ensuremath{\overline \squark}}\xspace$ quark contents, and their discovery will definitely help to understand the internal structure of strange tetraquark states~\cite{Cao:2021ton,Duan:2021pll}.

\begin{figure}[!tb]
\begin{center}
    \includegraphics[width=0.8\textwidth]{figs/production/JpsiPhiK.pdf}
    \caption{Distributions of (left) $\phi{\ensuremath{\kaon^+}}\xspace$, (middle)${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace\phi$  and (right) ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace{\ensuremath{\kaon^+}}\xspace$ invariant masses for (black data points) the ${\ensuremath{\Bu}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace\phi{\ensuremath{\kaon^+}}\xspace$ candidates  compared with  (red solid lines) the fit results using  (top row) the final model and (bottom row) the model of LHCb Run 1 analysis~\cite{LHCb-PAPER-2016-018}. Figure taken from Ref.~\cite{LHCb-PAPER-2020-044}.}
\label{sec:prod:fig:JpsiPhiK}
\end{center}
\end{figure}

\begin{table}[!hbt]
\begin{center}
    \caption{Spin-parity, significance, masses, widths, and fit fractions of exotic candidates observed in the ${\ensuremath{\Bu}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace\phi{\ensuremath{\kaon^+}}\xspace$ decay
    by LHCb~\cite{LHCb-PAPER-2020-044}.}
    \label{sec:prod:tab:B2JpsiPhiK}
\vskip -0.2cm
\def1.2{1.2}
\begin{tabular}{ccr@{\:$\pm$\:}lr@{\:$\pm$\:}lr@{\:$\pm$\:}l}
\hline
Contribution &Significance\,[$\times \sigma$]&            \multicolumn{2}{c}{$M_0$\,[$\kern -0.125em \aunit{Me\kern -0.1em V}\xspace$]}   & \multicolumn{2}{c}{$\Gamma_0$\,[$\kern -0.125em\aunit{Me\kern -0.1em V}\xspace$]}      & \multicolumn{2}{c}{FF\,[\%]}          \\
\hline \hline
$X(2^-)$ & & \multicolumn{2}{c}{} & \multicolumn{2}{c}{} &  \multicolumn{2}{c}{}  \\
$X(4150)$            & 4.8 (8.7) & $4146$&$18\pm33$ & $135$&$28\,_{-\,30}^{+\,59}$   & $2.0$&$0.5\,_{-\,1.0}^{+\,0.8}$ \\
\hline
$X(1^-)$ & & \multicolumn{2}{c}{} & \multicolumn{2}{c}{} &  \multicolumn{2}{c}{}  \\
$X(4630)$            & 5.5 (5.7) & $4626$&$ 16\,_{-\,110}^{+\,\phantom{0}18}$ & $174$&$27\,_{-\,\phantom{0}73}^{+\,134}$  & $2.6$&$0.5\,_{-\,1.5}^{+\,2.9}$ \\
\hline

All $X(0^+)$ & &\multicolumn{2}{c}{} &\multicolumn{2}{c}{} & $20$&$5\,_{-\,\phantom{0}7}^{+\,14}$  \\
$X(4500)$           & 20 (20) & $4474$&$3 \pm3 $ & $77$&$6\,_{-\,\phantom{0}8}^{+\,10}$ &  $5.6$&$0.7\,_{-\,0.6}^{+\,2.4}$ \\
$X(4700)$           & 17 (18) & $4694$&$ 4 \,_{-\,\phantom{0}3}^{+\,16}$ & $87$&$8\,_{-\,\phantom{0}6}^{+\,16}$ &  $8.9$&$1.2\,_{-\,1.4}^{+\,4.9}$ \\
\hline
All $X(1^+)$ & &\multicolumn{2}{c}{} &\multicolumn{2}{c}{} & $26$&$3\,_{-\,10}^{+\,\phantom{0}8}$ \\
$X(4140)$           & 13 (16) & $4118$&$ 11 \,_{-\,36}^{+\,19}$ & $162$&$21\,_{-\,49}^{+\,24}$  & $17$&$3\,_{-\,6}^{+\,19}$ \\
$X(4274)$           & 18 (18) & $4294$&$4 \,_{-\,6}^{+\,3}$     & $53$ &$5\pm5$                 & $2.8$&$0.5\,_{-\,0.4}^{+\,0.8}$ \\
$X(4685)$           & 15 (15) & $4684$&$ 7\,_{-\,16}^{+\,13}$   & $126$&$15 \,_{-\,41}^{+\,37}$ & $7.2$&$1.0 \,_{-\,2.0}^{+\,4.0}$ \\
\hline\hline
All $Z_{cs}(1^+)$ & & \multicolumn{2}{c}{}& \multicolumn{2}{c}{}& $25$&$5\,_{-\,12}^{+\,11}$ \\
$Z_{cs}(4000)$          & 15 (16) & $4003$&$ 6\,_{-\,14}^{+\,\phantom{0}4}$  & $131$&$15\pm26$  & $9.4$&$2.1\pm3.4$ \\
$Z_{cs}(4220)$          & 5.9 (8.4) & $4216$&$24\,_{-\,30}^{+\,43}$ & $233$&$52\,_{-\,73}^{+\,97}$ & $10$&$4\,_{-\,\phantom{0}7}^{+\,10}$ \\
\hline\hline
\end{tabular}
\end{center}
\end{table}

Exotic hadrons are also searched for in open charm final states using fully reconstructed beauty hadron decays. A Dalitz analysis of the ${\ensuremath{\Bu}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\D^+}}\xspace{\ensuremath{\D^-}}\xspace{\ensuremath{\kaon^+}}\xspace$ decay is performed by LHCb~\cite{LHCb-PAPER-2020-024,LHCb-PAPER-2020-025}, and two exotic states, $X_0(2900)$ and $X_1(2900)$, are required to have a good
fit to the ${\ensuremath{\D^-}}\xspace{\ensuremath{\kaon^+}}\xspace$ invariant mass spectrum. The  ${\ensuremath{\D^-}}\xspace{\ensuremath{\kaon^+}}\xspace$  invariant-mass distribution and the fit projections are shown on the left of Fig.~\ref{sec:prod:fig:X2900DD}. Their spin-parities are measured to be $J^P=0^+$ and $1^-$, and widths to be about $50\aunit{Me\kern -0.1em V}\xspace$ and $100\aunit{Me\kern -0.1em V}\xspace$ respectively. The quark contents of these $X$ states are
${\ensuremath{\overline \cquark}}\xspace {\ensuremath{\overline \squark}}\xspace du$.  It is noted that these two $X$ states
contribute up to $35\%$ of the total ${\ensuremath{\Bu}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\D^+}}\xspace{\ensuremath{\D^-}}\xspace{\ensuremath{\kaon^+}}\xspace$ decay branching fraction, a magnitude similar to those of conventional charmonia in the decay.
A final state rescattering effect is considered to explain  such a large branching fraction~\cite{Chen:2020eyu}. The observation of $X$ states in the ${\ensuremath{\Bu}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\D^+}}\xspace{\ensuremath{\D^-}}\xspace{\ensuremath{\kaon^+}}\xspace$ decay  opens a new avenue for studies of exotic hadrons composed of four different quark flavours~\cite{Lu:2020qmp,Liu:2020nil,Hu:2020mxp,Tan:2020cpu,Dong:2020rgs,Dai:2015bcc,Cheng:2020nho,Wang:2021lwy,Wang:2020xyc,Huang:2020ptc,He:2020jna,Chen:2020aos,Zhang:2020oze,Liu:2020orv}. Actually there are more than two dozens of $B\ensuremath{\rightarrow}\xspace D D K(\pi)$  decays, and also a few similar decays for $b$-baryons, and  it is promising that more $X$ states will be observed in these decays.

Prompt production of di-charm hadrons has been suggested to search for exotic states containing multiple charm quarks~\cite{Yang:2009zzp}.
In the invariant mass spectrum of ${\ensuremath{\D^0}}\xspace\Dz{\ensuremath{\pion^+}}\xspace$, shown on the right of Fig.~\ref{sec:prod:fig:X2900DD}, a new structure is observed close to the $D^{*+}{\ensuremath{\D^0}}\xspace$ mass
threshold~\cite{LHCb-PAPER-2021-031, LHCb-PAPER-2021-032}. The structure is measured to be consistent with the ground state of a $T_{cc}^+$ isoscalar
tetraquark,  with $J^P=1^+$ and quark contents $cc{\ensuremath{\overline \uquark}}\xspace{\ensuremath{\overline \dquark}}\xspace$. Its Breit-Wigner mass is measured to be
$-273\pm61\pm5^{+11}_{-14}\ensuremath{\aunit{ke\kern -0.1em V\!/}c^2}\xspace$ below $m_{D^{*+}}+m_{\ensuremath{\D^0}}\xspace$, and  it BW width is $\Gamma_\mathrm{BW}=410\pm165\pm43^{+18}_{-38}\aunit{ke\kern -0.1em V}\xspace$.
The same state also appears in the ${\ensuremath{\D^0}}\xspace\Dz$ and ${\ensuremath{\D^0}}\xspace{\ensuremath{\D^+}}\xspace$ mass spectra, with a ${\ensuremath{\pion^+}}\xspace$, ${\ensuremath{\pion^0}}\xspace$ or
$\gamma$ in the $T_{cc}^+\ensuremath{\rightarrow}\xspace {\ensuremath{\D^0}}\xspace\Dz{\ensuremath{\pion^+}}\xspace$, $T_{cc}^+\ensuremath{\rightarrow}\xspace{\ensuremath{\D^0}}\xspace{\ensuremath{\D^+}}\xspace{\ensuremath{\pion^0}}\xspace$ or $T_{cc}^+\ensuremath{\rightarrow}\xspace{\ensuremath{\D^0}}\xspace{\ensuremath{\D^+}}\xspace\gamma$ decays undetected respectively.
Dedicated studies of the $T_{cc}^+$ resonance lineshape are performed using a unitarised Breit-Wigner distribution, considering $T_{cc}^+$ decays in ${\ensuremath{\D^0}}\xspace\Dz{\ensuremath{\pion^+}}\xspace, {\ensuremath{\D^0}}\xspace{\ensuremath{\D^+}}\xspace{\ensuremath{\pion^0}}\xspace$ and ${\ensuremath{\D^0}}\xspace\Dz\gamma$ final states. The pole mass of the resonance in this advanced model is measured to be $-360\pm4^{+4}_{-0}\aunit{ke\kern -0.1em V}\xspace$
below $m_{\ensuremath{\D^*}}\xspace+m_{\ensuremath{\D^0}}\xspace$, and the pole width is $\Gamma_\mathrm{pole}=-48\pm2^{+0}_{-14}\aunit{ke\kern -0.1em V}\xspace$. This extremely narrow width has caused many theoretical interests~\cite{Ling:2021bir,Meng:2021jnw,Dai:2021wxi,Li:2021zbw,Guo:2021mja,Chen:2021kad}.  The $T_{cc}^+$ state is the first observed
tetraquark candidates with two heavy quarks of the same flavour. Many theoretical models have been applied to explain the existence and structure of such a tetraquark state with a large fraction of them favouring a $D^{*+}{\ensuremath{\D^0}}\xspace$ hadron molecule interpretation~\cite{Xing:2021yid,An:2021vwi,Deng:2018kly,Wang:2017dtg,Karliner:2017qjm,Eichten:2017ffp,Zhu:2016arf,Yang:2019itm,Yan:2018gik,Tan:2020ldi,Kong:2021ohg,Hsiao:2021tyq,Chen:2021erj,Dong:2021lkh,He:2021oyy,Majarshin:2021hex,Wang:2021ajy,Du:2021zzh,Ling:2021asz,Baru:2021ldu,Ren:2021dsi,Wang:2017uld,Chen:2021htr}.

\begin{figure}[!tb]
\begin{center}
    \includegraphics[width=0.47\textwidth]{figs/production/X2900.pdf}
    \includegraphics[width=0.4\textwidth]{figs/production/Tcc.pdf}
    \caption{(Left)  Invariant-mass distribution of ${\ensuremath{\D^-}}\xspace{\ensuremath{\kaon^+}}\xspace$  in ${\ensuremath{\Bu}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\D^+}}\xspace{\ensuremath{\D^-}}\xspace{\ensuremath{\kaon^+}}\xspace$ decays in data and the fit projections, taken from Ref.~\cite{LHCb-PAPER-2020-025}. The yellow dot-dashed and red dotted line correspond to $X_0(2900)$ and $X_1(2900)$ respectively, and more details can be found in Ref.~\cite{LHCb-PAPER-2020-025}.
    (Right) Distribution of ${\ensuremath{\D^0}}\xspace\Dz{\ensuremath{\pion^+}}\xspace$ invariant mass overlaid with fit projections.  The contribution of the non-${\ensuremath{\D^0}}\xspace$ background has been statistically subtracted. Vertical dashed lines show the $D^{*+}{\ensuremath{\D^0}}\xspace$ and ${\ensuremath{\D^{*0}}}\xspace{\ensuremath{\D^+}}\xspace$ mass thresholds respectively for a comparison. Figure taken from Ref.~\cite{LHCb-PAPER-2021-032}.
    }
\label{sec:prod:fig:X2900DD}
\end{center}
\end{figure}

\begin{figure}[!tb]
\begin{center}
    \includegraphics[width=0.58\textwidth]{figs/production/DiJpsi.pdf}
    \includegraphics[width=0.41\textwidth]{figs/production/PcStates.pdf}
    \caption{
    (Left) Invariant mass distribution of the combination of two ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace$ mesons overlaid with fit projections,  taken from Ref.~\cite{LHCb-PAPER-2020-011}. The fit model doesn't contain interference between any component.
    (Right) Invariant mass distribution of ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace p$ fitted with three BW functions plus polynomial background, taken from Ref.~\cite{LHCb-PAPER-2020-011}. The vertical lines mark the ${\ensuremath{\PSigma}}\xspace^+_c{\ensuremath{\Dbar{}^0}}\xspace$  and the ${\ensuremath{\PSigma}}\xspace^+_c {\ensuremath{\Dbar{}^{*0}}}\xspace$ mass thresholds respectively.
    }
\label{sec:prod:figs1}
\end{center}
\end{figure}

In addition to two open charm final states,  di-${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace$ mass spectrum of prompt production is also studied using full LHCb data~\cite{LHCb-PAPER-2020-011}. Two peaking structures are observed in the mass range
$6.2<m_{{\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace\jpsi}<7.4\ensuremath{\aunit{Ge\kern -0.1em V\!/}c^2}\xspace$, where fully charmed tetraquark states are predicted~\cite{Iwasaki:1976cn,Chao:1980dv,PhysRevD.25.2370,PhysRevD.29.426,Badalian:1985es,
Berezhnoy:2011xn,Wu:2016vtq,Karliner:2016zzc,
Barnea:2006sd,Debastiani:2017msn,Liu:2019zuc,
Chen:2016jxd,
Wang:2019rdo,Bedolla:2019zwg,Lloyd:2003yc,
Chen:2020lgj,Wang:2018poa,
Anwar:2017toa,Esposito:2018cwh,
Becchi:2020mjz,Bai:2016int,Richard:2017vry,Vega-Morales:2017pmm,Chen:2019vrj,Berezhnoy:2012tu}. The first
structure (referred to as the threshold peak) covers the range between $6.2$ and $6.6\ensuremath{\aunit{Ge\kern -0.1em V\!/}c^2}\xspace$ close to the di-${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace$ mass threshold, and the
other one sits at $6.9\ensuremath{\aunit{Ge\kern -0.1em V\!/}c^2}\xspace$, as shown on the left of Fig.~\ref{sec:prod:figs1}. The di-${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace$ mass spectrum is modelled with a combination of
BW functions for the two peaking structures and empirical smooth functions for SPS, DPS production of nonpeaking background. When no interference between BW and SPS is applied, the
threshold structure can be described by two BW functions and  the one at $6.9\ensuremath{\aunit{Ge\kern -0.1em V\!/}c^2}\xspace$ is well described by a BW.
The narrow structure, denoted as $X(6900)$, is measured to have a mass of $m_{X(6900)}=6905\pm11\pm7\ensuremath{\aunit{Me\kern -0.1em V\!/}c^2}\xspace$ and a width of
$\Gamma_{X(6900)}=80\pm19\pm33\aunit{Me\kern -0.1em V}\xspace$.
Interpretations of the threshold peak using feeddown decays from excited quarkonium pairs are also possible. For various fit models without any interference, the dip around $6.8\ensuremath{\aunit{Ge\kern -0.1em V\!/}c^2}\xspace$ can't be well described. Advance fit studies are performed introducing interference between SPS and resonant
structures. In one such fit, two BW functions are considered: a broad BW interfering with SPS used to describe the
threshold structure, and a stand-alone narrow one used to model the $6.9\ensuremath{\aunit{Ge\kern -0.1em V\!/}c^2}\xspace$ peak. This new model could fit well the overall spectrum, and the
broad structure is now measured to have a mass around $6.7\ensuremath{\aunit{Ge\kern -0.1em V\!/}c^2}\xspace$ and a width of about $0.3\aunit{Ge\kern -0.1em V}\xspace$, while the $X(6900)$ structure has
a mass consistent with the no-interference fit model, but its width becomes about twice larger. As the fit results for the broad
structure is not stable in different models, its nature is not fully resolved and more data are needed to provide
better information.
The $X(6900)$ state is consistent with a genuine fully charmed tetraquark, however, in some models possible origins due to rescatterings of multiple charmonia or dibaryon molecules \mbox{\itshape etc.}\xspace are also discussed~\cite{Chao:2020dml,Weng:2020jao,An:2020vku,Lu:2020cns,Lee:2011rka,Jin:2020jfc,Lu:2017dvm,Meng:2017fwb,Wang:2020dlo,Tang:2019nwv,Wan:2020fsk,Yang:2020wkh,Li:2019uch,Liu:2020tqy,Zhao:2020zjh,Yan:2021glh,Liu:2021rtn,Wang:2021hdk,Li:2021ygk,Ke:2021iyh,Huang:2020dci,An:2020jix,Gong:2020bmg,Zhu:2020snb,Zhang:2020vpz,Cao:2020gul,Guo:2020pvt,Zhu:2020xni,Zhang:2020xtb,Faustov:2020qfm,Feng:2020riv,Ma:2020kwb,Dong:2020nwy,Karliner:2020dta,Wang:2020wrp,Yang:2020rih,Deng:2020iqw}. Other fully heavy tetraquarks, such as
${\ensuremath{\Pb}}\xspace{\ensuremath{\overline \bquark}}\xspace{\ensuremath{\Pc}}\xspace{\ensuremath{\overline \cquark}}\xspace$ and ${\ensuremath{\Pb}}\xspace{\ensuremath{\overline \bquark}}\xspace{\ensuremath{\Pb}}\xspace{\ensuremath{\overline \bquark}}\xspace$,
can be searched for in $\PUpsilon{\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace$ and $\PUpsilon\PUpsilon$ or similar final states, however current data are estimated to have small sensitivities for these states, demanding the increased luminosity in LHCb upgrades.

\paragraph{Pentaquark states} Following the successful discoveries of tetraquarks with $Q\overline{Q}$ contents, pentaquark states with $Q\overline{Q}$ were predicted, in the form of either meson-baryon hadron molecules or compact five-quark hadrons~\cite{Yang:2011wz,Zhu:2015bba,Li:2012bt,Liu:2018bkx}. Since 2015, several pentaquark candidates are found in LHCb, as summarised in Table~\ref{sec:prod:tab:Pc}.
The first observation of them is made in ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace p$ final states in ${\ensuremath{\Lz^0_\bquark}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace p{\ensuremath{\kaon^-}}\xspace$ decays through
an amplitude analysis~\cite{LHCb-PAPER-2015-029}. Two states were reported, $P_c(4380)^+$ and $P_c(4450)^+$, and their
evidence is also found in the ${\ensuremath{\Lz^0_\bquark}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace p {\ensuremath{\pion^-}}\xspace$ decays with a similar amplitude study~\cite{LHCb-PAPER-2016-015}.  A model independent moment analysis of the
${\ensuremath{\Lz^0_\bquark}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace p{\ensuremath{\kaon^-}}\xspace$ decay concludes that
contributions of ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace p$ exotics are essential, since only allowing ${\ensuremath{\PLambda}}\xspace^*\ensuremath{\rightarrow}\xspace p{\ensuremath{\kaon^-}}\xspace$ resonances in the decay are not sufficient to describe
data~\cite{LHCb-PAPER-2016-009}.  With full LHCb data, an amplitude analysis of ${\ensuremath{\Lz^0_\bquark}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace p{\ensuremath{\kaon^-}}\xspace$ decays becomes computationally very difficult. On the other hand, benefiting from the high statistics, one dimensional ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace p$ mass spectrum is
investigated to look for narrow pentaquark states~\cite{LHCb-PAPER-2019-014}. In this new analysis, the $P_c(4450)^+$ structure is found to consist of two narrow overlapping peaks
$P_c(4440)^+$ and $P_c(4457)^+$, and a new structure $P_c(4312)^+$ is observed. It is noted that the $P_c(4312)^+$ and $P_c(4457)^+$
states are close to the ${\ensuremath{\PSigma}}\xspace^+_c{\ensuremath{\Dbar{}^0}}\xspace$ and ${\ensuremath{\PSigma}}\xspace^+_c {\ensuremath{\Dbar{}^{*0}}}\xspace$ mass
thresholds respectively, as shown on the right of Fig.~\ref{sec:prod:figs1}, making them ideal candidates of meson-baryon
molecules.
Recently, LHCb reported the evidence of a new pentaquark state, $P_c(4337)^+$, in the ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace p({\ensuremath{\overline \proton}}\xspace)$ mass spectrum of
${\ensuremath{\B^0_\squark}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace p{\ensuremath{\overline \proton}}\xspace$ decays~\cite{LHCb-PAPER-2021-018}. This possible state is different from those observed in ${\ensuremath{\Lz^0_\bquark}}\xspace$ decays, making beauty meson decays a new place to search for pentaquarks.

Pentaquark candidates with strangeness are predicted in the ${\ensuremath{\Xires^-_\bquark}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace {\ensuremath{\PLambda}}\xspace{\ensuremath{\kaon^-}}\xspace$ decay~\cite{Wang:2019nvm,Chen:2015sxa}, which is an analogy of the ${\ensuremath{\Lz^0_\bquark}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace p{\ensuremath{\kaon^-}}\xspace$ channel by replacing the $d$ quark by the $s$ quark. An amplitude analysis of the ${\ensuremath{\Xires^-_\bquark}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace {\ensuremath{\PLambda}}\xspace{\ensuremath{\kaon^-}}\xspace$ decay is performed using full LHCb data, resulting in the evidence of a
new state, $P_{cs}(4459)^0$, in the ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace {\ensuremath{\PLambda}}\xspace$ mass spectrum, with quark contents, ${\ensuremath{\Pc}}\xspace{\ensuremath{\overline \cquark}}\xspace{\ensuremath{\Pu}}\xspace{\ensuremath{\Pd}}\xspace{\ensuremath{\Ps}}\xspace$. Its mass is measured to be
$4458.8\pm2.9^{+4.7}_{-1.1}\ensuremath{\aunit{Me\kern -0.1em V\!/}c^2}\xspace$  and width to be $17.3\pm6.5^{+8.0}_{-5.7}\aunit{Me\kern -0.1em V}\xspace$~\cite{LHCb-PAPER-2020-039}. According
to the prediction in Ref.~\cite{Wang:2019nvm}, there are two ${\ensuremath{\Xires_\cquark}}\xspace{\ensuremath{\D^*}}\xspace$ hadron molecules with masses within a few $\ensuremath{\aunit{Me\kern -0.1em V\!/}c^2}\xspace$ around the observed
structure. If the $P_{cs}(4459)^0$ structure is found to be composed of two two nearby states with LHCb upgrade data, it will be a strong proof of the molecular interpretation of such states.

Observations of pentaquark states with heavy-quark contents have triggered many studies from theorists, with the purpose of understanding their nature. In general, these states  are considered to either be compact tetraquarks, hadronic molecules or just bumps due to kinematic rescatterings~\cite{Wang:2019nwt,Du:2019pij,Meng:2019ilv,Weng:2019ynv,Xu:2019zme,Ozdem:2017exj,Huang:2015uda,Du:2021bgb,Chen:2021cfl,Hu:2021nvs,Yalikun:2021bfm,Lu:2021irg,Du:2021fmf,Zhu:2021lhd,Nakamura:2021dix,Shi:2021wyt,Phumphan:2021tta,Nakamura:2021qvy,Xiao:2021rgp,Chen:2020kco,Wang:2020eep,Peng:2020gwk,Xu:2020gjl,Giachino:2020rkj,Chen:2019thk,Wang:2019ato,Pimikov:2019dyr,Wang:2019hyc,Yamaguchi:2019seo,Cheng:2019obk,Wang:2019got,Fernandez-Ramirez:2019koa,Xiao:2019mvs,Guo:2019kdc,He:2019ify,Guo:2019fdo,Chen:2019bip,Liu:2018nse,Ferretti:2018ojb,Hiyama:2018ukv,Qin:2018dqp,Richard:2017una,He:2016pfa,Eides:2017xnt,Yamaguchi:2017zmn,Dong:2017gaw,Lin:2017mtz,Chen:2016ryt,Guo:2016bkl,Santopinto:2016pkp,Shen:2016tzq,Shimizu:2016rrd,Lu:2016nnt,Oset:2016lyh,Chen:2016heh,Wang:2015ixb,Yang:2015bmv,Burns:2015dwa,Scoccola:2015nia,Wang:2015epa,Ghosh:2015xqp,Li:2015gta,Anisovich:2015cia,Mikhasenko:2015vca,Lebed:2015tna,Liu:2015fea,He:2015cea,Maiani:2015vwa,Guo:2015umn,Mironov:2015ica,Roca:2015dva,Chen:2015moa,Chen:2015loa,Karliner:2015ina}.
Measurements of their $J^P$ and production properties, and finding their flavour partners will help to shed light on the problem.
However, it is likely that debates on the nature of exotic hadrons will continue before we have a complete and coherent theory to explain all of them. In the past years, phenomenology models on the hadron spectroscopy evolved quickly and some patterns have been revealed, for example, dynamics close to two-hadron mass thresholds~\cite{Dong:2020hxe}. Besides, lattice QCD simulations have improved in computational performances and application scopes substantially over the years and will play more and more important roles in our understanding of low-energy QCD~\cite{Chiu:2006hd,Guo:2013nja,Bi:2015ifa,Liu:2016kbb,Leskovec:2014gxa,Gayer:2021xzv,Cheung:2020mql}. Hopefully, one day we can  predict low energy QCD phenomena as precise as other parts of the SM.

\begin{table}[!hbt]
\begin{center}
    \caption{Detection decay channels, experimental significance, masses and widths of pentaquark states reported by LHCb.}
    \label{sec:prod:tab:Pc}
    \begin{tabular}{c|cccc}
\hline
        State           & Decays     & Significance [$\sigma$] & Mass [$\ensuremath{\aunit{Me\kern -0.1em V\!/}c^2}\xspace$] & Width [$\aunit{Me\kern -0.1em V}\xspace$]\\
        \hline
        $P_c(4312)^+$   & ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace p$  & $7.3\sigma$             & $4311.9\pm0.7^{+6.8}_{-0.6}$   & $\phantom{0}9.8\pm2.7^{+\phantom{0}3.7}_{-\phantom{0}4.5}$\\
        $P_c(4440)^+$   & ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace p$  & $5.4\sigma$             & $4440.3\pm1.3^{+4.1}_{-4.7}$   & $20.6\pm4.9^{+\phantom{0}8.7}_{-10.1}$\\
        $P_c(4457)^+$   & ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace p$  & $5.4\sigma$             & $4457.3\pm0.6^{+4.1}_{-1.7}$   & $0.53\pm2.0^{+\phantom{0}5.7}_{-\phantom{0}1.9}$\\

        $P_c(4337)^+$   & ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace p$  & $3.1\sigma$             & $\phantom{0}4337\phantom{0}^{+\phantom{0}7}_{-\phantom{0}4}\phantom{0}\pm2$&
        $\phantom{0}\phantom{0}29^{+\phantom{0}26}_{-\phantom{0}12}\pm14$\\
\hline
        $P_{cs}(4459)^0$   & ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace {\ensuremath{\PLambda}}\xspace$  & $3.1\sigma$        & $4458.8\pm2.9^{+4.7}_{-1.1}$   & $17.3\pm6.5^{+\phantom{0}8.0}_{-\phantom{0}5.7}$\\
        \hline
\hline
\end{tabular}
\end{center}
\end{table}

\clearpage

\section{Rare beauty hadron decays}
\label{sec:raredecay}

By rare decays we mainly refer to   flavour-changing-neutral-current (FCNC) processes, which are  highly suppressed in the SM by the Glashow-Iliopoulos-Maiani (GIM)  mechanism~\cite{Glashow:1970gm}.
However, the rare decays could receive significant contributions from new particles or new interactions beyond the SM.
Precision measurements  of their properties play a special role in search of  physics beyond the SM.

The LHCb collaboration has given priority to the study of FCNC $b\ensuremath{\rightarrow}\xspace s$ transitions, focusing on theoretically clean observables such as decay rates of purely leptonic $B$-meson decays, angular coefficients  in $b\ensuremath{\rightarrow}\xspace s \ensuremath{\ell^+ \ell^-}\xspace$ decays, and ratio  of decay rates  between $b\ensuremath{\rightarrow}\xspace s \ensuremath{\ell^+ \ell^-}\xspace$ processes with different lepton flavours.
Analyses of $pp$ collision data  collected in the Run 1 and Run 2 periods have already led to some very important and interesting findings, including but not limited to the first observation of the purely leptonic decay ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace$, anomalous angular distribution in the decay  ${\ensuremath{\B^0}}\xspace \ensuremath{\rightarrow}\xspace {\ensuremath{\kaon^{*0}}}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace$,  and  indication of lepton flavour universality violation in ${\ensuremath{\Bu}}\xspace \ensuremath{\rightarrow}\xspace {\ensuremath{\kaon^+}}\xspace \ensuremath{\ell^+ \ell^-}\xspace$ decays.

\subsection{Effective field theory for \texorpdfstring{$b\ensuremath{\rightarrow}\xspace s$}{TEXT} transitions}

The effective Hamiltonian describing the quark-level $b\ensuremath{\rightarrow}\xspace s$ transitions is given by~\cite{Buras:1994dj, Buras:1998raa, Buchalla:1995vs, Altmannshofer:2008dz,Kruger:2005ep,Descotes-Genon:2011nqe,Lunghi:2006hc}
\begin{equation}
   H_{\rm eff}(b\ensuremath{\rightarrow}\xspace s) =-{\frac{4G_F}{\sqrt{2}}} {\ensuremath{V_{\tquark\bquark}^{\phantom{\ast}}}}\xspace{\ensuremath{V_{\tquark\squark}^\ast}}\xspace \sum_{i=1}^{10} C_i \Ope{i},
   \label{eq:hamil}
\end{equation}
with $\mathcal{O}_{i}$ denoting the local operators in the SM and  $C_i$ indicating the corresponding Wilson coefficients.
Of particular interest is the electromagnetic dipole operator corresponding to penguin diagrams mediated by photons,
\begin{equation}
   \Ope{7} = {\dfrac{e}{ 16 \pi^2}} m_b ({\ensuremath{\overline \squark}}\xspace \sigma_{\mu\nu} P_R b)  F^{\mu\nu},
   \label{eq:O7}
\end{equation}
and  semileptonic operators corresponding to loop diagrams mediated by $Z^0$ or ${\ensuremath{\PW^\pm}}\xspace$ bosons,
\begin{equation}
   \Ope{9} ={\dfrac{e^2}{16 \pi^2}} ({\ensuremath{\overline \squark}}\xspace \gamma_{\mu} P_L b) (\overline{\ell} \gamma^{\mu} \ell), \,\,
   \Ope{10} ={\dfrac{e^2}{16 \pi^2}} ({\ensuremath{\overline \squark}}\xspace \gamma_{\mu} P_L b) (\overline{\ell} \gamma^{\mu} \gamma_5 \ell),
   \label{eq:O9O10}
\end{equation}
where $P_L=(1-\gamma_5)/2$ and $ P_R=(1+\gamma_5)/2$.
Contribution from  physics  beyond the SM can either alter the values of the Wilson coefficients and/or give rise to new operators that are absent or highly suppressed in the SM, such as the scalar and pseudo-scalar operators
\begin{equation}
   \Ope{S} ={\dfrac{e^2}{16 \pi^2}} ({\ensuremath{\overline \squark}}\xspace P_L b) (\overline{ell}  \ell), \,\,
   \Ope{P} ={\dfrac{e^2}{16 \pi^2}} ({\ensuremath{\overline \squark}}\xspace P_L b) (\overline{\ell} \gamma_5 \ell)\;,
   \label{eq:OsOp}
\end{equation}
and the  chirality-flipped operators $\mathcal{O}_{7,8,9,10,S,P}'$,  which are obtained by changing $P_{L(R)}$ to $P_{R(L)}$ in $\mathcal{O}_{7,8,9,10,S,P}$.

The Wilson coefficients $C^{(\prime)}_7$ can be probed in radiative $b$-hadron decays such as ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace \phi \gamma$ and ${\ensuremath{\Lz^0_\bquark}}\xspace \ensuremath{\rightarrow}\xspace {\ensuremath{\PLambda}}\xspace \gamma$, $C^{(\prime)}_{9,10}$ probed in semileptonic decays such as $B^0\ensuremath{\rightarrow}\xspace {\ensuremath{\kaon^+}}\xspace\ensuremath{\ell^+ \ell^-}\xspace$ and $B^+\ensuremath{\rightarrow}\xspace {\ensuremath{\kaon^{*0}}}\xspace  \ensuremath{\ell^+ \ell^-}\xspace$, and $C^{(\prime)}_{S,P}$ probed in purely leptonic decays such as $B^0_s \ensuremath{\rightarrow}\xspace \ensuremath{\ell^+ \ell^-}\xspace$.
Recent LHCb results on $b\ensuremath{\rightarrow}\xspace s$ transitions are summarised in the remaining of this section.

\subsection{Purely leptonic \texorpdfstring{$B$ meson decays}{TEXT}}
The decays ${\ensuremath{\B_{(\squark)}^0}}\xspace\ensuremath{\rightarrow}\xspace \ensuremath{\ell^+ \ell^-}\xspace (\ell=\mu,\,e)$ are  among the most interesting probes of new physics.
They are theoretically clean and are expected to be extremely rare  in the SM due to  helicity suppression in addition to the FCNC loop suppression.
Their branching fractions in the SM are precisely predicted to be~\cite{Beneke:2019slt}
\begin{equation}
\begin{aligned}
{\ensuremath{\mathcal{B}}}\xspace({\ensuremath{\B^0_\squark}}\xspace\ensuremath{\rightarrow}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace) &= (3.66\pm 0.14 )\times 10^{-9}, \\
{\ensuremath{\mathcal{B}}}\xspace({\ensuremath{\B^0}}\xspace\ensuremath{\rightarrow}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace) &= (1.03\pm 0.05 )\times 10^{-10}, \\
{\ensuremath{\mathcal{B}}}\xspace({\ensuremath{\B^0_\squark}}\xspace\ensuremath{\rightarrow}\xspace {\ensuremath{\Pe^+\Pe^-}}\xspace) &= (8.60\pm 0.36 )\times 10^{-14},  \\
{\ensuremath{\mathcal{B}}}\xspace({\ensuremath{\B^0}}\xspace\ensuremath{\rightarrow}\xspace {\ensuremath{\Pe^+\Pe^-}}\xspace) &= (2.41\pm 0.13 )\times 10^{-15}.
\end{aligned}
\label{eq:b2llSM}
\end{equation}
Note the  ${\ensuremath{\B^0}}\xspace \ensuremath{\rightarrow}\xspace \ensuremath{\ell^+ \ell^-}\xspace$ decays proceed via $b\ensuremath{\rightarrow}\xspace d$ transitions, thus are further suppressed with respect to the ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace \ensuremath{\ell^+ \ell^-}\xspace$ decays
by a factor of $|{\ensuremath{V_{\tquark\dquark}^{\phantom{\ast}}}}\xspace/{\ensuremath{V_{\tquark\squark}^{\phantom{\ast}}}}\xspace|^2\sim \lambda^2$.
The decay rates of ${\ensuremath{\B_{(\squark)}^0}}\xspace\ensuremath{\rightarrow}\xspace \ensuremath{\ell^+ \ell^-}\xspace$ processes are highly sensitive to  (pseudo-)scalar interactions beyond the SM.

A joint analysis of data from the LHCb and CMS experiments  collected in Run~1 led to the observation of the ${\ensuremath{\B^0_\squark}}\xspace\ensuremath{\rightarrow}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace$ decay with a significance exceeding six standard deviations, and determined the branching fraction to be ${\ensuremath{\mathcal{B}}}\xspace({\ensuremath{\B^0_\squark}}\xspace\ensuremath{\rightarrow}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace)=(2.8^{+0.7}_{-0.6})\times 10^{-9}$~\cite{LHCb-PAPER-2014-049}.
This result was later updated by LHCb~\cite{LHCB-PAPER-2017-001} and CMS~\cite{CMS:2019bbr}  by adding the 2016 data.
The ATLAS collaboration  reported an evidence
for the decay ${\ensuremath{\B^0_\squark}}\xspace\ensuremath{\rightarrow}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace$ with a significance of $4.6\sigma$  using the  data collected between 2011 and 2016.
However, no significant signal for the decay ${\ensuremath{\B^0}}\xspace\ensuremath{\rightarrow}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace$ has been found by any experiment yet.
A combination of the results from ATLAS, CMS and LHCb gives the branching fraction ${\ensuremath{\mathcal{B}}}\xspace({\ensuremath{\B^0_\squark}}\xspace\ensuremath{\rightarrow}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace)=(2.69^{+0.37}_{-0.35})\times 10^{-9}$, and sets an upper limit of ${\ensuremath{\mathcal{B}}}\xspace({\ensuremath{\B^0}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\Pmu^+\Pmu^-}}\xspace)<1.9\times 10^{-10}$ at 95\% confidence level~\cite{LHCb-CONF-2020-002}.
Figure~\ref{fig:b2mumu_comb} shows the constraints on ${\ensuremath{\mathcal{B}}}\xspace({\ensuremath{\B^0_\squark}}\xspace\ensuremath{\rightarrow}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace)$ and ${\ensuremath{\mathcal{B}}}\xspace({\ensuremath{\B^0_\squark}}\xspace\ensuremath{\rightarrow}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace)$ from  the three experiments and the combined results, which  are compatible with the SM predictions~\cite{Beneke:2019slt}  within 2.1 standard deviations. The difference is mainly driven by the ATLAS results,  which have an optimal solution  outside the physical region and are slightly in tension with the SM predictions.

Very recently, LHCb reported updated results on the ${\ensuremath{\B^0_\squark}}\xspace\ensuremath{\rightarrow}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace$ and ${\ensuremath{\B^0}}\xspace\ensuremath{\rightarrow}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace$ decays  using all data collected in Run 1 and Run 2.
The results are ${\ensuremath{\mathcal{B}}}\xspace({\ensuremath{\B^0_\squark}}\xspace\ensuremath{\rightarrow}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace)=(3.09^{+0.46+0.15}_{-0.43-0.11})\times 10^{-9}$,
${\ensuremath{\mathcal{B}}}\xspace({\ensuremath{\B^0}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\Pmu^+\Pmu^-}}\xspace)<2.6\times 10^{-10}$ at 95\% confidence level~\cite{LHCb-PAPER-2021-007,LHCb-PAPER-2021-008}, which are in good agreement with the SM expectations, as shown in Fig.~\ref{fig:b2mumu2021}.

In addition to the decay rate, other interesting quantities, such as effective lifetime and CP asymmetry, can also be measured to search for possible non-SM contribution to the decay ${\ensuremath{\B^0_\squark}}\xspace\ensuremath{\rightarrow}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace$~\cite{DeBruyn:2012wk}.
Pioneering studies  of its  effective lifetime   using all   data collected in Run 1 and Run 2 have been performed at LHCb, leading to the result of  $\tau^{\rm eff}({\ensuremath{\B^0_\squark}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\Pmu^+\Pmu^-}}\xspace) =2.07\pm 0.29 \pm 0.03$ \ensuremath{\aunit{ps}}\xspace~\cite{LHCb-PAPER-2021-007,LHCb-PAPER-2021-008}.
As a long-term goal, the effective lifetime and time-dependent CP violation of the decay ${\ensuremath{\B^0_\squark}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\Pmu^+\Pmu^-}}\xspace$ will be fully exploited in LHCb upgrade II, which aims to accumulate a $pp$ collision data sample of 300\ensuremath{\fb^{-1}}\xspace~\cite{LHCb-PII-Physics}.

\begin{figure}[!tb]
\begin{center}
\includegraphics[width=0.45\textwidth]{./figs/raredecay/b2mumu_comb1.pdf}
\includegraphics[width=0.45\textwidth]{./figs/raredecay/b2mumu_comb2.pdf}
\caption{
(Left) Likelihood contours in the ${\ensuremath{\mathcal{B}}}\xspace({\ensuremath{\B^0_\squark}}\xspace\ensuremath{\rightarrow}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace)$ - ${\ensuremath{\mathcal{B}}}\xspace({\ensuremath{\B^0_\squark}}\xspace\ensuremath{\rightarrow}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace)$ plane corresponding to $-2\Delta L =2.3,\, 11.6$ and $11.8$, for ATLAS, CMS and LHCb experiments and  the combination.
(Right) Likelihood contours corresponding to $-2\Delta L =2.3,\, 11.6,\, 11.8,\, 19.3$ and $30.2$. The data sets used were collected from 2011 to 2016.
Figures are extracted from Refs.~\cite{LHCb-CONF-2020-002}.
}
\label{fig:b2mumu_comb}
\end{center}
\end{figure}

\begin{figure}
\begin{center}
\includegraphics[width=0.45\textwidth]{./figs/raredecay/B2mumu_9invfb_mass.pdf}
\includegraphics[width=0.45\textwidth]{./figs/raredecay/B2mumu_9invfb_br.pdf}
\caption{ (Left) Invariant mass distribution of the selected ${\ensuremath{\B_{(\squark)}^0}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\Pmu^+\Pmu^-}}\xspace$ candidates with output of the used multivariate classifier above 0.5, superimposed with the fit result. (Right) Confidence intervals in the plane of the branching fractions of ${\ensuremath{\B^0_\squark}}\xspace\ensuremath{\rightarrow}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace$ and ${\ensuremath{\B^0}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\Pmu^+\Pmu^-}}\xspace$.
The data sets used for the blue solid  (brown dotted)  contours were collected from 2011 to 2018 (2016). Figures are extracted from Refs.~\cite{LHCb-PAPER-2021-007,LHCb-PAPER-2021-008}. }
\label{fig:b2mumu2021}
\end{center}
\end{figure}

The decays  ${\ensuremath{\B^0_\squark}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\Pe^+\Pe^-}}\xspace$ are even more rare than ${\ensuremath{\B^0_\squark}}\xspace\ensuremath{\rightarrow}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace$ and can provide powerful tests of lepton flavour universality.
To date, the searches performed by the CDF and LHCb experiments have found no evidence for either ${\ensuremath{\B^0_\squark}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\Pe^+\Pe^-}}\xspace$ or ${\ensuremath{\B^0}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\Pe^+\Pe^-}}\xspace$.
The most stringent upper limits  on their branching fractions are
${\ensuremath{\mathcal{B}}}\xspace({\ensuremath{\B^0_\squark}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\Pe^+\Pe^-}}\xspace)<11.2\times 10^{-9}$ and
${\ensuremath{\mathcal{B}}}\xspace({\ensuremath{\B^0}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\Pe^+\Pe^-}}\xspace)<3.0\times 10^{-9}$ at 95\% confidence level
set by LHCb~\cite{LHCb-PAPER-2020-001}.

In addition to decays to dileptons, the  LHCb experiment  has also searched for decays of neutral $B$ mesons to four leptons.
A recent search using  the full Run~1 and Run~2 data sample found no hint of such decays, and  upper limits at 90\% confidence level for the nonresonant decays are determined to be ${\ensuremath{\mathcal{B}}}\xspace({\ensuremath{\B^0_\squark}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\Pmu^+\Pmu^-}}\xspace\mumu)<8.6\times10^{-10}$ and  ${\ensuremath{\mathcal{B}}}\xspace({\ensuremath{\B^0}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\Pmu^+\Pmu^-}}\xspace\mumu)<1.8\times10^{-10}$~\cite{LHCb-PAPER-2021-039}.
More stringent limits are set for decays involving the ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace$ resonance or a promptly decaying intermediate scalar particle with a mass of $1\ensuremath{\aunit{Ge\kern -0.1em V\!/}c^2}\xspace$.

\subsection{Semileptonic \texorpdfstring{$b\ensuremath{\rightarrow}\xspace s \ensuremath{\ell^+ \ell^-}\xspace$}{TEXT} decays}
Semileptonic ${\ensuremath{\Pb}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\Ps}}\xspace\ensuremath{\ell^+ \ell^-}\xspace$ decays provide valuable insight into possible non-SM contributions that affect  the Wilson coefficients $C_9$ and $C_{10}$ of the electromagnetic operators.
The presence of hadrons in the final state makes the search for new physics in semileptonic decays more complicated than that in purely leptonic decays.
The challenges in hadronic form-factor calculations lead to significant uncertainties in the SM predictions of their decay rates.
Fortunately, a number of clean observables that are less affected by the form factors than the total decay rates have been identified, including  some special observables in angular distributions and observables for lepton universality test. A comprehensive study of these observables in a series of  $b\ensuremath{\rightarrow}\xspace s \ensuremath{\ell^+ \ell^-}\xspace$ processes by the LHCb collaboration has led to several anomalous results,  which collectively and coherently indicate possible non-SM  contribution to $b\ensuremath{\rightarrow}\xspace s \ensuremath{\ell^+ \ell^-}\xspace$ transitions.

\subsubsection{Differential decay rates with respect to $q^2$}
\label{dBdq2}

The differential branching fraction $d{\mathcal B}/dq^2$ can be measured in intervals of $q^2$, the invariant mass squared of the lepton pair, and compared with   SM predictions.
Calculations of form factors are needed for making the SM predictions.
Such calculations are challenging, and require different treatments depending on the $q^2$ regions.
Light-cone sum rule calculations~\cite{Altmannshofer:2014rta, Bharucha:2015bzk, Altmannshofer:2015sma}  and lattice QCD calculations~\cite{Horgan:2013hoa, Horgan:2015vla} are often  used to determine the form factors in low- and high-$q^2$ regions, respectively.

After measuring the branching fraction of the decay ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace \phi {\ensuremath{\Pmu^+\Pmu^-}}\xspace$ to be about $3\sigma$ below the SM expectation value~\cite{LHCb-PAPER-2015-023}, the LHCb experiment further studied its differential branching fraction  as a function of $q^2$. The top left plot in Fig.~\ref{fig:b2smumu_br} shows the latest results of  $d{\ensuremath{\mathcal{B}}}\xspace({\ensuremath{\B^0_\squark}}\xspace\ensuremath{\rightarrow}\xspace \phi {\ensuremath{\Pmu^+\Pmu^-}}\xspace)/dq^2$ obtained using Run 1 data~\cite{LHCb-PAPER-2021-014}, where the ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace$ and ${\ensuremath{\Ppsi{(2S)}}}\xspace$ regions are excluded.
A puzzle appears in the  range $1<q^2<6\gevgevcccc$.
 The  branching fraction integrated over this range is measured to be $(2.88\pm 0.22 ) \times 10^{-8}$.  Currently, the most precise SM prediction for this range is   $(5.37\pm 0.66 ) \times 10^{-8}$, obtained from a combination of light-cone sum rule and lattice QCD calculations. A discrepancy of  3.6  standard deviations is observed. Similar patterns are also seen in the LHCb measurements of differential branching fractions in the decays ${\ensuremath{\Lz^0_\bquark}}\xspace\ensuremath{\rightarrow}\xspace {\ensuremath{\PLambda}}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace$~\cite{LHCb-Paper-2015-009},
  ${\ensuremath{\B^0}}\xspace \ensuremath{\rightarrow}\xspace {\ensuremath{\kaon^{*0}}}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace$~\cite{LHCb-Paper-2016-012},
${\ensuremath{\Bu}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\kaon^+}}\xspace{\ensuremath{\Pmu^+\Pmu^-}}\xspace$,
${\ensuremath{\B^0}}\xspace \ensuremath{\rightarrow}\xspace {\ensuremath{\kaon^0_{\mathrm{S}}}}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace$ and
${\ensuremath{\Bu}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\kaon^{*+}}}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace$~\cite{LHCb-Paper-2014-006}, as shown in Fig.~\ref{fig:b2smumu_br}.

\begin{figure}[!tb]
\begin{center}
\includegraphics[width=0.45\textwidth]{./figs/raredecay/phimumu_f2.pdf}
\includegraphics[width=0.45\textwidth]{./figs/raredecay/Lb2Lmdmumu_br.pdf}
\includegraphics[width=0.45\textwidth]{./figs/raredecay/Bd2kstmumu_br.pdf}
\includegraphics[width=0.45\textwidth]{./figs/raredecay/Bplus2Kstmumu_br.pdf}
\includegraphics[width=0.45\textwidth]{./figs/raredecay/Bd2Kmumu_br.pdf}
\includegraphics[width=0.45\textwidth]{./figs/raredecay/Bplus2Kmumu_br.pdf}
\caption{
Top left: measured $d{\ensuremath{\mathcal{B}}}\xspace({\ensuremath{\B^0_\squark}}\xspace\ensuremath{\rightarrow}\xspace \phi {\ensuremath{\Pmu^+\Pmu^-}}\xspace)/dq^2$~\cite{LHCb-PAPER-2021-014} overlaid with the SM predictions based on light-cone sum rule calculations~\cite{Altmannshofer:2014rta, Bharucha:2015bzk, Altmannshofer:2015sma}  at low $q^2$ and lattice QCD calculations~\cite{Horgan:2013hoa, Horgan:2015vla} at high $q^2$;
top right: measured $d{\ensuremath{\mathcal{B}}}\xspace({\ensuremath{\Lz^0_\bquark}}\xspace \ensuremath{\rightarrow}\xspace {\ensuremath{\PLambda}}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace)/dq^2$~\cite{LHCb-Paper-2015-009} compared with the SM predictions using form factors from lattice QCD calculations~\cite{Detmold:2012vy};
middle left: measured $d{\ensuremath{\mathcal{B}}}\xspace({\ensuremath{\B^0}}\xspace \ensuremath{\rightarrow}\xspace {\ensuremath{\kaon^{*0}}}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace)/dq^2$~\cite{LHCb-Paper-2016-012} compared with the SM predictions using form factors from line-cone sum rule and lattice QCD calculations~\cite{Bharucha:2015bzk,Horgan:2013hoa};
middle right, bottom left and bottom right:
measured $d{\ensuremath{\mathcal{B}}}\xspace({\ensuremath{\Bu}}\xspace \ensuremath{\rightarrow}\xspace {\ensuremath{\kaon^{*+}}}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace)/dq^2$,
$d{\ensuremath{\mathcal{B}}}\xspace({\ensuremath{\B^0}}\xspace \ensuremath{\rightarrow}\xspace {\ensuremath{\kaon^0_{\mathrm{S}}}}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace)/dq^2$ and
$d{\ensuremath{\mathcal{B}}}\xspace({\ensuremath{\Bu}}\xspace \ensuremath{\rightarrow}\xspace {\ensuremath{\kaon^+}}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace)/dq^2$~\cite{LHCb-Paper-2014-006}
compared with the SM predictions using form factors from line-cone sum rule and lattice QCD calculations~\cite{Bobeth:2011gi, Bobeth:2011nj}.
}
\label{fig:b2smumu_br}
\end{center}
\end{figure}

\subsubsection{Angular distributions}

Angular distributions in $b\ensuremath{\rightarrow}\xspace s \ensuremath{\ell^+ \ell^-}\xspace$ decays contain rich information about  interference  between the SM and non-SM contributions that may not  be accessible via  decay rates integrated over angular variables. A set of  $q^2$-dependent angular coefficients   can be extracted from the angular distributions and used as probes for new physics, which are complementary to  branching fractions and $d{\ensuremath{\mathcal{B}}}\xspace/dq^2$.
Based on these  coefficients,  we can define some theoretically clean observables with reduced dependency on form factors.

Of particular interest is the angular distribution of the $\mbox{{\ensuremath{\B^0}}\xspace\ensuremath{\rightarrow}\xspace {\ensuremath{\kaon^{*0}}}\xspace(\ensuremath{\rightarrow}\xspace {\ensuremath{\kaon^+}}\xspace{\ensuremath{\pion^-}}\xspace) {\ensuremath{\Pmu^+\Pmu^-}}\xspace}$ decay , which  has been extensively studied by \mbox{BaBar}\xspace~\cite{BaBar:2006tnv}, \mbox{Belle}\xspace~\cite{Belle:2016fev}, ATLAS~\cite{ATLAS:2018gqc}, CMS~\cite{CMS:2017rzx} and LHCb~\cite{LHCb-PAPER-2011-020,LHCb-PAPER-2013-019, LHCb-PAPER-2015-051, LHCb-PAPER-2020-002}.
Followinng the definitions in Ref.~\cite{Altmannshofer:2008dz},
the CP-averaged angular distribution of the decay ${\ensuremath{\B^0}}\xspace\ensuremath{\rightarrow}\xspace {\ensuremath{\kaon^{*0}}}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace$ with ${\ensuremath{\kaon^{*0}}}\xspace\ensuremath{\rightarrow}\xspace {\ensuremath{\kaon^+}}\xspace{\ensuremath{\pion^-}}\xspace$ is given by
\begin{equation}
\begin{aligned}
    \frac{1}{d(\Gamma + \bar{\Gamma})/dq^2}
    \frac{d^4(\Gamma +\bar{\Gamma})}{dq^2 d\vec{\Omega}}
    =  \frac{9}{32\pi}&  \left[ \tfrac{3}{4}(1-F_{\rm L})\sin^2\theta_K
    +F_{\rm L}\cos^2\theta_K\right. \\
    & +  \tfrac{1}{4} (1-F_{\rm L})\sin^2\theta_K \cos2\theta_\ell
    -F_{\rm L}\cos^2\theta_K \cos 2\theta_\ell \\
    &+ S_4 \sin 2\theta_K \sin 2\theta_\ell \cos\phi + S_5 \sin 2\theta_K \sin \theta_\ell \sin \phi \\
    &+ \tfrac{4}{3} A_{\rm FB} \sin^2\theta_K \cos \theta_\ell + S_7 \sin 2\theta_K \sin\theta_\ell \sin \phi \\
    &+ \left. S_8 \sin 2\theta_K \sin 2 \theta_\ell \sin \phi + S_9 \sin^2\theta_K \sin^2\theta_\ell \sin 2\phi \right],
\end{aligned}
\label{eq:kstarmumu}
\end{equation}
where $\vec{\Omega}=(\cos\theta_K, \cos \theta_\ell, \phi)$, $\theta_K$ is the angle between the directions of the ${\ensuremath{\kaon^+}}\xspace$ (${\ensuremath{\kaon^-}}\xspace$) and  ${\ensuremath{\B^0}}\xspace$ (${\ensuremath{\Bbar{}^0}}\xspace$) in the rest frame of the ${\ensuremath{\kaon^{*0}}}\xspace$ (${\ensuremath{\Kbar{}^{*0}}}\xspace$) system, $\theta_\ell$ is the angle between the direction of the  ${\ensuremath{\Pmu^+}}\xspace$ ($\mun$) and the opposite direction of the ${\ensuremath{\B^0}}\xspace$ (${\ensuremath{\Bbar{}^0}}\xspace$) in the rest frame of the ${\ensuremath{\Pmu^+\Pmu^-}}\xspace$ system, $\phi$ is the angle between the plane defined by the muon pair and the plane defined by the kaon and pion in the ${\ensuremath{\B^0}}\xspace$ (${\ensuremath{\Bbar{}^0}}\xspace$) rest frame.
Eight observables  can be extracted, including  the fraction of the longitudinal polarisation of the ${\ensuremath{\kaon^{*0}}}\xspace$ meson ($F_{\rm L}$),
the forward-backward asymmetry of the ${\ensuremath{\Pmu^+\Pmu^-}}\xspace$ system ($A_{\rm FB}$),  and six other angular coefficients ($S_i$, $i=3,4,5,7,8,9$).
Using the $S_i$ coefficients,
new  observables less sensitive to form factor uncertainties are defined, such as  $P^{\prime}_i =S_i/\sqrt{F_{\rm L}(1-F_{\rm L})}$ $(i=3,4,5)$.

The latest LHCb results on  angular analysis of the ${\ensuremath{\B^0}}\xspace\ensuremath{\rightarrow}\xspace {\ensuremath{\kaon^{*0}}}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace$ decay are obtained  using data collected in 2011, 2012 and 2016~\cite{ LHCb-PAPER-2020-002}. The majority of the angular  observables
are  consistent with the SM predictions~\cite{Altmannshofer:2014rta, Descotes-Genon:2014uoa} obtained from a combination of  light-cone sum rule calculations~\cite{Bharucha:2015bzk} for low-$q^2$ regions and lattice QCD calculations~\cite{Horgan:2013hoa, Horgan:2015vla}  for high-$q^2$ regions.
A clear exception is seen with the robust observable  $P^{\prime}_5$ defined using $S_5$, as shown in  Fig.~\ref{fig:p5pr_puzzle}.
The measured values of $P^{\prime}_5$ in the intervals $4.0<q^2<6.0\gevgevcccc$ and $6.0<q^2<8.0\gevgevcccc$ are found to be higher than the SM predictions~\cite{Descotes-Genon:2014uoa} by 2.5 and 2.9 standard deviations, respectively.
These results confirm the  discrepancy in $P^{\prime}_5$ observed in an earlier LHCb analysis with Run 1 data~\cite{LHCb-PAPER-2015-051}.  According to model-independent fits using the FLAVIO software package~\cite{Straub:2018kue},  the overall tension with the SM  is increased  from 3.0 to 3.3 standard deviations.
The fits reveal that the current measurements of the angular observables in ${\ensuremath{\B^0}}\xspace\ensuremath{\rightarrow}\xspace {\ensuremath{\kaon^{*0}}}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace$  can be accommodated by shifting the real part of the Wilson coefficient $C_9$ from its SM value by $0.99^{+0.25}_{-0.21}$~\cite{ LHCb-PAPER-2020-002}.

\begin{figure}[!tb]
\begin{center}
\includegraphics[width=0.45\textwidth]{./figs/raredecay/Kstmumu_FL.pdf}
\includegraphics[width=0.45\textwidth]{./figs/raredecay/Kstmumu_AFB.pdf}
\includegraphics[width=0.45\textwidth]{./figs/raredecay/Kstmumu_S5.pdf}
\includegraphics[width=0.45\textwidth]{./figs/raredecay/Kstmumu_p5pr.pdf}
\caption{ Results of the observables $F_{\rm L}$, $A_{\rm FB}$, $S_5$ and $P^{\prime}_5$ in bins of $q^2$, compared with the SM predictions.
Figures are extracted from Ref.~\cite{LHCb-PAPER-2020-002}. }
\label{fig:p5pr_puzzle}
\end{center}
\end{figure}

Recently, LHCb reported  results of  angular analysis of the $B^+\ensuremath{\rightarrow}\xspace K^{*+}(\ensuremath{\rightarrow}\xspace {\ensuremath{\kaon^0_{\mathrm{S}}}}\xspace{\ensuremath{\pion^+}}\xspace) {\ensuremath{\Pmu^+\Pmu^-}}\xspace$ decay using Run~1 and Run~2 data\cite{LHCb-PAPER-2020-041}.
A trend of deviations from the SM predictions in~$P^{\prime}_5$, similar to that in the isospin partner decay $B^0\ensuremath{\rightarrow}\xspace {\ensuremath{\kaon^{*0}}}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace$, is shown in the left of Fig.~\ref{fig:Bp2Kstmumu}.
Meanwhile,  a large  discrepancy in the measurement of  $P_2=\frac{2}{3}A_{\rm FB}/(1-F_{\rm L})$  has also been observed in the  $6.0<q^2<8.0\gevgevcccc$ region, where the measurement deviates from its SM prediction~\cite{Altmannshofer:2014rta, Descotes-Genon:2014uoa} by 3.0 standard deviations (Fig.~\ref{fig:Bp2Kstmumu} right ).

\begin{figure}
\begin{center}
\includegraphics[width=0.45\textwidth]{./figs/raredecay/Bp2kstmumu_p5pr.pdf}
\includegraphics[width=0.45\textwidth]{./figs/raredecay/Bp2kstmumu_p2.pdf}
\caption{ Results for the observables  $P^{\prime}_5$ and $P_2$ in $B^+\ensuremath{\rightarrow}\xspace K^{*+} {\ensuremath{\Pmu^+\Pmu^-}}\xspace$,
 compared with SM predictions. Figures are extracted from Ref.~\cite{LHCb-PAPER-2020-041}. }
\label{fig:Bp2Kstmumu}
\end{center}
\end{figure}

An untagged time-integrated angular analysis of the decay ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace \phi {\ensuremath{\Pmu^+\Pmu^-}}\xspace$ has  been performed by the LHCb collaboration using data collected in 2011 to 2012 and 2016 to 2018~\cite{LHCb-PAPER-2021-022}.
In this channel, the same particles ($\mu^+$ and ${\ensuremath{\kaon^+}}\xspace$) are used to define the angular variables for both ${\ensuremath{\B^0_\squark}}\xspace$ and ${\ensuremath{\Bbar{}^0_\squark}}\xspace$ decays, since the final state is not self-tagging as in the ${\ensuremath{\B^0}}\xspace\ensuremath{\rightarrow}\xspace {\ensuremath{\kaon^{*0}}}\xspace{\ensuremath{\Pmu^+\Pmu^-}}\xspace$ case.
With this convention, the coefficients of the terms in the CP-averaged  time-integrated angular distributions corresponding to the interference between CP-even ($0$ or $||$) and CP-odd ($\perp$ or $S$) amplitudes  are the CP asymmetries $A^{{\ensuremath{C\!P}}\xspace}_{\rm FB}$ and $A_{5,8,9}$, rather than the CP-average observables $A_{\rm FB}$ and $S_{5,8,9}$ in Eq.\ref{eq:kstarmumu}.
These asymmetries can arise from either direct CP violation or  nonzero effective mixing phase,  with the latter contribution suppressed by the small value of $\Delta\Gamma_s/\Gamma_s$.
They are  predicted to be close to  zero in the SM but has some sensitivity to new physics contributions~\cite{Bobeth:2008ij}.
The measurements of $A^{CP}_{\rm FB}$ and $A_{5,8,9}$ in intervals of $q^2$ are shown in Fig.~\ref{fig:Bs2phimumu}, which are consistent with CP invariance.
Much more information on  CP violation can be obtained from time-dependent angular analysis of tagged ${\ensuremath{\B^0_\squark}}\xspace\ensuremath{\rightarrow}\xspace \phi {\ensuremath{\Pmu^+\Pmu^-}}\xspace$ decays~\cite{Descotes-Genon:2015hea}, which may become feasible with the huge amount of  data  that will be collected with the upgraded LHCb detector in the coming data-taking periods.

\begin{figure}[!tb]
\begin{center}
\includegraphics[width=0.45\textwidth]{./figs/raredecay/phimumu_1.pdf}
\includegraphics[width=0.45\textwidth]{./figs/raredecay/phimumu_2.pdf}
\includegraphics[width=0.45\textwidth]{./figs/raredecay/phimumu_3.pdf}
\includegraphics[width=0.45\textwidth]{./figs/raredecay/phimumu_4.pdf}
\caption{ CP asymmetries $A^{CP}_{\rm FB}$ and $A_{5,8,9}$ in intervals of $q^2$ in the decay ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace \phi {\ensuremath{\Pmu^+\Pmu^-}}\xspace$ measured using Run 1 and Run 2 data (black) in Run 1 data only (grey). Figures are extracted from Ref.~\cite{LHCb-PAPER-2021-022}. }
\label{fig:Bs2phimumu}
\end{center}
\end{figure}

Angular coefficients in  the $B^+ \ensuremath{\rightarrow}\xspace K^+ {\ensuremath{\Pmu^+\Pmu^-}}\xspace$ decay
have  been measured by LHCb~\cite{LHCb-PAPER-2014-007} and CMS~\cite{CMS:2018qih} using Run 1 data.
Angular observables in ${\ensuremath{\Lz^0_\bquark}}\xspace \ensuremath{\rightarrow}\xspace p{\ensuremath{\kaon^-}}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace$ decays have been determined by LHCb~\cite{LHCb-PAPER-2018-029} from a moment analysis using data collected between 2011 and 2016.
These results are  consistent with the SM predictions but limited by statistical uncertainties.

\subsubsection{Lepton flavour  universality tests}
In the SM, the couplings of  the three generations of leptons to the electroweak gauge bosons $Z^0$ and ${\ensuremath{\PW^\pm}}\xspace$
are assumed to be identical. This is known as lepton flavour  universality (LFU). Under this assumption,  processes involving the three flavours of charged leptons, $e$, $\mu$ and $\tau$, have equal  rates up to corrections caused by different lepton masses, which can be trivially taken into account. Contributions of new particles or new interactions may lead to violation of LFU, particularly in FCNC processes such as $b\ensuremath{\rightarrow}\xspace s \ensuremath{\ell^+ \ell^-}\xspace$ decays.
Stringent tests of LFU can be performed by measuring the ratio of the branching fractions between $B\ensuremath{\rightarrow}\xspace X {\ensuremath{\Pmu^+\Pmu^-}}\xspace$ and $B\ensuremath{\rightarrow}\xspace X {\ensuremath{\Pe^+\Pe^-}}\xspace$ decays~\cite{Hiller:2003js,Bordone:2016gaq,Isidori:2020acz} outside the charmonium regions in the dilepton mass spectrum, with $X$ indicating the hadron(s) in the decays.
The ratio is denoted by
\begin{equation}
    R_{X} \equiv\frac{{\ensuremath{\mathcal{B}}}\xspace(B\ensuremath{\rightarrow}\xspace X {\ensuremath{\Pmu^+\Pmu^-}}\xspace)}
                {{\ensuremath{\mathcal{B}}}\xspace(B \ensuremath{\rightarrow}\xspace X {\ensuremath{\Pe^+\Pe^-}}\xspace)}\,,
\end{equation}
In the SM, $R_X$ is expected to be very close to unity with negligible theoretical uncertainty,
due to the small and precisely known difference between the  muon and electron masses.
On the experimental side, reconstruction of electrons is challenging due to the Bremsstrahlung radiation. The decays $B\ensuremath{\rightarrow}\xspace X {\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace(\ensuremath{\rightarrow}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace)$ and $B\ensuremath{\rightarrow}\xspace X {\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace(\ensuremath{\rightarrow}\xspace {\ensuremath{\Pe^+\Pe^-}}\xspace)$ are used as control channels for cancellation of systematic uncertainties associated with electron reconstruction.
Experimentally, the ratio $R_X$ is measured using a double-ratio technique following the equation
\begin{equation}
    R_X =\frac{{\ensuremath{\mathcal{B}}}\xspace(B\ensuremath{\rightarrow}\xspace X {\ensuremath{\Pmu^+\Pmu^-}}\xspace)/{\ensuremath{\mathcal{B}}}\xspace(B\ensuremath{\rightarrow}\xspace X {\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace(\ensuremath{\rightarrow}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace))} {{\ensuremath{\mathcal{B}}}\xspace(B \ensuremath{\rightarrow}\xspace X {\ensuremath{\Pe^+\Pe^-}}\xspace)/{\ensuremath{\mathcal{B}}}\xspace(B\ensuremath{\rightarrow}\xspace X {\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace (\ensuremath{\rightarrow}\xspace e^+e^-))}\;,
\end{equation}
where ${\ensuremath{\mathcal{B}}}\xspace(B\ensuremath{\rightarrow}\xspace X {\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace(\ensuremath{\rightarrow}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace))/{\ensuremath{\mathcal{B}}}\xspace(B\ensuremath{\rightarrow}\xspace X {\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace(\ensuremath{\rightarrow}\xspace {\ensuremath{\Pe^+\Pe^-}}\xspace))$ is known to be very close to unity~\cite{PDG2020}.

The LHCb collaboration  previously  measured  $R_K=0.846^{+0.060}_{-0.054}\aunit{(stat)}\xspace{}^{+0.016}_{-0.014}\aunit{(syst)}\xspace$
in ${\ensuremath{\Bu}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\kaon^+}}\xspace\ensuremath{\ell^+ \ell^-}\xspace\; (\ell=e,\,\mu)$ decays in the dilepton
mass-squared range $1.1<q^2<6.0\gevgevcccc$
using Run 1 and part of Run 2 data~\cite{LHCb-PAPER-2019-009}.
The result was below the SM expectation~\cite{Straub:2018kue} by 2.5 standard deviations.
Recently, LHCb  updated the $R_K$ measurement
using the full Run 1 and Run 2 sample~\cite{LHCb-PAPER-2021-004}.
The mass distributions of the ${\ensuremath{\Bu}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\kaon^+}}\xspace\ensuremath{\ell^+ \ell^-}\xspace$ candidates are shown in Fig.~\ref{fig:RK1}.
The $R_K$ value is measured to be
\begin{equation*}
R_K(1.1<q^2<6.0\gevgevcccc)=
0.846^{+0.042}_{-0.039}\aunit{(stat)}\xspace{}^{+0.013}_{-0.012}\aunit{(syst)}\xspace\;,
\end{equation*}
 which deviates from the SM prediction,  $1.00\pm 0.01$~\cite{Descotes-Genon:2015uva, Bobeth:2007dw, Bordone:2016gaq,van_dyk_danny_2016_159680, Straub:2018kue},
 by 3.1 standard deviations, giving the first evidence for LFU violation in $B^+ \ensuremath{\rightarrow}\xspace K^+ \ensuremath{\ell^+ \ell^-}\xspace$ decays.
A comparison of the LHCb $R_K$ result with the values measured by  \mbox{BaBar}\xspace~\cite{BaBar:2012mrf} and \mbox{Belle}\xspace~\cite{Belle:2009zue} is shown in Fig.~\ref{fig:RK2}.

Tests of LFU have also been performed
in other $b\ensuremath{\rightarrow}\xspace s \ensuremath{\ell^+ \ell^-}\xspace\;(\ell=e,\mu)$ decays.
Based on the Run 1 data sample, the LHCb collaboration has determined
the ratios of branching fractions of
$B^0\ensuremath{\rightarrow}\xspace {\ensuremath{\kaon^{*0}}}\xspace \ensuremath{\ell^+ \ell^-}\xspace\; (\ell=e,\mu)$  decays
in two regions of dilepton mass-squared below the ${\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace$ resonance to be~\cite{LHCb-PAPER-2017-013}
\begin{equation*}
\begin{aligned}
R_{{\ensuremath{\kaon^{*0}}}\xspace}(0.045<q^2<1.1\gevgevcccc) & = 0.66^{+0.11}_{-0.07} \aunit{(stat)}\xspace \pm 0.03 \aunit{(syst)}\xspace\;, \\
R_{{\ensuremath{\kaon^{*0}}}\xspace}(1.1<q^2<6.0\gevgevcccc) & =  0.69^{+0.11}_{-0.07} \aunit{(stat)}\xspace\pm 0.05 \aunit{(syst)}\xspace\;.
\end{aligned}
\end{equation*}
These results are in tension with the SM predictions~\cite{Bordone:2016gaq,Descotes-Genon:2015uva,Capdevila:2016ivx,Capdevila:2017ert,Serra:2016ivr,van_dyk_danny_2016_159680, Straub:2018kue,Bharucha:2015bzk,Jager:2014rwa} at the level
of $2.1-2.3$ and $2.4-2.5$ standard deviations, respectively.
LHCb has also measured
\begin{equation*}
R_{pK}(0.1<q^2<6.0\gevgevcccc) =0.86^{+0.14}_{-0.11} \aunit{(stat)}\xspace\pm 0.05 \aunit{(syst)}\xspace
\end{equation*}
in ${\ensuremath{\Lz^0_\bquark}}\xspace\ensuremath{\rightarrow}\xspace p{\ensuremath{\kaon^-}}\xspace\ensuremath{\ell^+ \ell^-}\xspace\;(\ell=e,\mu)$ decays
using Run 1 and part of Run 2 data~\cite{LHCb-PAPER-2019-040}.

Very recently, LHCb reported the observation of
the decays ${\ensuremath{\B^0}}\xspace\ensuremath{\rightarrow}\xspace {\ensuremath{\kaon^0_{\mathrm{S}}}}\xspace {\ensuremath{\Pe^+\Pe^-}}\xspace$ and  ${\ensuremath{\Bu}}\xspace\ensuremath{\rightarrow}\xspace {\ensuremath{\kaon^{*+}}}\xspace {\ensuremath{\Pe^+\Pe^-}}\xspace$ and
the measurements of  LFU observables $R_{{\ensuremath{\kaon^0_{\mathrm{S}}}}\xspace}$ and $R_{{\ensuremath{\kaon^{*+}}}\xspace}$
 using the full Run 1 and Run 2 data samples~\cite{LHCb-PAPER-2021-038}.
 The obtained results of $R_{{\ensuremath{\kaon^0_{\mathrm{S}}}}\xspace}$ and $R_{{\ensuremath{\kaon^{*+}}}\xspace}$  are
\begin{equation*}
    \begin{aligned}
    R_{{\ensuremath{\kaon^0_{\mathrm{S}}}}\xspace}(1.1<q^2<6.0\gevgevcccc) & = 0.66^{+0.20}_{-0.14} \aunit{(stat)}\xspace{}^{+0.02}_{-0.04} \aunit{(syst)}\xspace\;, \\
R_{{\ensuremath{\kaon^{*+}}}\xspace}(0.045<q^2<6.0\gevgevcccc) & = 0.70^{+0.11}_{-0.07} \aunit{(stat)}\xspace{}^{+0.03}_{-0.04} \aunit{(syst)}\xspace\;,
    \end{aligned}
\end{equation*}
which are lower than but consistent with the SM predictions at 1.5 and 1.4 standard deviations, respectively.

In summary,  some hints of LFU violation are observed  in the low-$q^2$ region of a series of $b\ensuremath{\rightarrow}\xspace s \ensuremath{\ell^+ \ell^-}\xspace$ decays.
The branching fractions of  the muon modes are systematically lower than those of the electron modes  in the low-$q^2$ region, though the  measurements are currently with limited precision. This  is consistent with the finding that
 the differential decay rates in  several $b \ensuremath{\rightarrow}\xspace s {\ensuremath{\Pmu^+\Pmu^-}}\xspace$ decays are below SM predictions  in the low-$q^2$ region.
 Meanwhile, there are also   puzzles in the angular distributions of $B\ensuremath{\rightarrow}\xspace {\ensuremath{\kaon^*}}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace$ decays in the low-$q^2$ region.
 These anomalies  point to potential contributions of new physics in  $b \ensuremath{\rightarrow}\xspace s \ensuremath{\ell^+ \ell^-}\xspace$ transitions~\cite{Alguero:2019ptt,Ciuchini:2020gvn, Aebischer:2019mlg,Alguero:2021anc,Altmannshofer:2021qrr,Geng:2021nhg,Cornella:2021sby,Isidori:2021vtc,Isidori:2021tzd,Hurth:2021nsi}. Particularly, new physics scenarios that mainly affect  $b \ensuremath{\rightarrow}\xspace s {\ensuremath{\Pmu^+\Pmu^-}}\xspace$ transitions are preferred, according to
  global analysis in the framework of SM effective field theory~\cite{Aebischer:2019mlg,Alguero:2021anc,Altmannshofer:2021qrr,Geng:2021nhg,Isidori:2021tzd,Hurth:2021nsi}.  Improved understanding of the  anomalies can  be achieved using the huge amount of data that will be recorded in the Run 3 data taking period  and beyond.
  This will not only significantly increase the precision of the benchmark measurements in $b \ensuremath{\rightarrow}\xspace s$ transitions, but also provide great opportunities to explore new observables and new rare decay modes that are currently inaccessible, such as time-dependent observables in ${\ensuremath{\B^0}}\xspace \ensuremath{\rightarrow}\xspace {\ensuremath{\kaon^0_{\mathrm{S}}}}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace$~\cite{Descotes-Genon:2020tnz,Kosnik:2021owl} and ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace \phi {\ensuremath{\Pmu^+\Pmu^-}}\xspace$~\cite{Descotes-Genon:2015hea}   decays,
 and   lepton universality ratios and angular observables in heavily suppressed  $b \ensuremath{\rightarrow}\xspace d \ensuremath{\ell^+ \ell^-}\xspace$ transitions~\cite{Bordone:2021olx,Rusov:2019ixr,Soni:2020bvu,Kindra:2018ayz}.

\begin{figure}[!tb]
\begin{center}
\includegraphics[width=0.45\textwidth]{./figs/raredecay/RK2021_f1a.pdf}
\includegraphics[width=0.45\textwidth]{./figs/raredecay/RK2021_f1b.pdf}
\caption{ Invariant mass distributions of the selected (left) ${\ensuremath{\Bu}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\kaon^+}}\xspace {\ensuremath{\Pe^+\Pe^-}}\xspace$  and (right) ${\ensuremath{\Bu}}\xspace\ensuremath{\rightarrow}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace$ candidates, superimposed by the fit results. Figures are extracted from Ref.~\cite{LHCb-PAPER-2021-004}.
}
\label{fig:RK1}
\end{center}
\end{figure}

\begin{figure}[!tb]
\begin{center}
\includegraphics[width=0.45\textwidth]{./figs/raredecay/RK2021_f2.pdf}
\caption{ Comparison of the LHCb $R_K$ result with the measurements by B-factories~\cite{BaBar:2012mrf,Belle:2009zue}.
The vertical line indicates the SM expected value. Figures are extracted from Ref.~\cite{LHCb-PAPER-2021-004}.
}
\label{fig:RK2}
\end{center}
\end{figure}

\subsection{Radiative \texorpdfstring{$b\ensuremath{\rightarrow}\xspace s \gamma$}{TEXT} decays}

The effective Hamiltonian for  $b\ensuremath{\rightarrow}\xspace s \gamma$ transitions can be approximately written as
\begin{equation}
   H_{\rm eff}(b\ensuremath{\rightarrow}\xspace s \gamma) =-{\dfrac{4G_F}{\sqrt{2}}} V_{tb}V^*_{ts} ( C_7 \Ope{7} + C^{\prime}_7 \mathcal{O}^{\prime}_7),
   \label{eq:hamil2}
\end{equation}
where only  the leading operator $\mathcal{O}_{7}$ and its chirality-flipped counterpart  $\mathcal{O}^{\prime}_7$ are included. In the SM, the   coefficient $C^{\prime}_7$ is given by $C^{\prime}_7=\frac{m_s}{m_b} C_7$ due to the chiral $V-A$ structure of the weak interaction, where $m_s$ ($m_b$) indicates the mass of the $s$ ($b$) quark.
Consequently, the photons emitted in radiative $b$-hadron (${\ensuremath{\overline \bquark}}\xspace$-hadron) decays
 are predominantly left-handed (right-handed).
 Amplitudes with  right-handed photons, $A_R$, are suppressed by the ratio $m_s/m_b$  compared with those with left-handed photons, $A_L$, but could be enhanced in new physics scenarios with right-handed charged current, such as supersymmetric grand unified theories and  left-right symmetric models~\cite{Atwood:1997zr, Everett:2001yy,Grinstein:2004uu, Becirevic:2012dx,Kou:2013gna,Haba:2015gwa,Paul:2016urs}.

Rich information on photon polarisation can be obtained from time-dependent analysis of  $B_q^0 (\bar{B}^0_q) \ensuremath{\rightarrow}\xspace f_{\rm CP} \gamma$ decays~\cite{Atwood:1997zr, Atwood:2004jj, Muheim:2008vu}, where $f_{\rm CP}$ is a CP eigenstate with eigenvalue $\eta$. The time-dependent decay rates summing over  left-handed and right-handed photons are expressed as
\begin{equation}
  P(t) =P_0 e^{{-\Gamma_q t}}
  \left[\cosh(\Delta\Gamma_q t/2)-A^{\Delta} \sinh(\Delta\Gamma_q t/2) +\xi  C\cos(\Delta m_qt)-\xi S \sin(\Delta m_q t) \right]\;,
   \label{eq:b2sg_rate}
\end{equation}
where $\xi$ takes the value of $+1$ ($-1$) for an initial $B^0_q$ ($\bar{B}^0_q$) meson.
The coefficient $C$ quantifies CP violation in the decay. This type of CP violation  has been constrained to be small in radiative $B$ meson decays by \mbox{BaBar}\xspace, \mbox{Belle}\xspace and \mbox{LHCb}\xspace~\cite{BaBar:2012fqh,Belle:2017hum,LHCb-PAPER-2012-019}.
Assume no CP violation in the decay for simplicity, it is convenient to write
\begin{equation}
\begin{aligned}
&\mathrm{Favoured:}\,\,\,\phantom{0}    A(\bar{B}^0_q \ensuremath{\rightarrow}\xspace f_{\rm CP} \gamma_{L}) = a_{L} e^{i \delta_{L}} e^{i \phi_{L}} \,\;\Rightarrow\;
  A(B^0_q \ensuremath{\rightarrow}\xspace f_{\rm CP} \gamma_{R}) =\eta \, a_{L} e^{i \delta_{L}} e^{-i\phi_{L}}\,,\\
&\mathrm{Suppressed:\,}   A(\bar{B}^0_q \ensuremath{\rightarrow}\xspace f_{\rm CP} \gamma_{R}) = a_{R} e^{i \delta_{R}} e^{i \phi_{R}}\; \Rightarrow\;A(B^0_q \ensuremath{\rightarrow}\xspace f_{\rm CP} \gamma_{L}) =\eta \,  a_{R} e^{i \delta_{R}} e^{-i\phi_{R}}\,,
\end{aligned}
\label{eq:}
\end{equation}
where $a_{L(R)}$, $\delta_{L(R)}$ and $\phi_{L(R)}$ are the size, strong phase and weak phase of \mbox{$A(\bar{B}^0_q \ensuremath{\rightarrow}\xspace f_{\rm CP} \gamma_{L(R)})$}, respectively, and $\eta$ is the CP parity of $f_{\rm CP}$.
The terms $\sin(\Delta m_q)$ and $\cosh(\Delta \Gamma_q t/2)$ in Eq.~\ref{eq:b2sg_rate}
arise from  interference of the amplitudes of direct decay, $A(\bar{B^0_q} \ensuremath{\rightarrow}\xspace f_{\rm CP} \gamma_L)$ or $A(B^0_q \ensuremath{\rightarrow}\xspace f_{\rm CP} \gamma_R)$, and the decay via $B^0_q$-$\bar{B}^0_q$ mixing, $p/q A(B^0_q \ensuremath{\rightarrow}\xspace f_{\rm CP} \gamma_L)$ or $q/p A(\bar{B}^0_q \ensuremath{\rightarrow}\xspace f_{\rm CP} \gamma_R)$.
The mixing-induced observables $S$ and $A^{\Delta}$ are  given by~\cite{ Muheim:2008vu}
\begin{equation}
\begin{aligned}
 S &\approx  \frac{2\eta\, r}{1+r^2} \cos (\delta_L-\delta_R)  \sin (\phi_q-\phi_L-\phi_R) \,,\\
 A^{\Delta } &\approx  \frac{2\eta\, r}{1+r^2} \cos (\delta_L-\delta_R)  \cos (\phi_q-\phi_L-\phi_R) \;,
\end{aligned}
\end{equation}
where  $r\equiv|a_R/a_L|\approx|C^{\prime}_7/C_7|$, $\phi_q$ is the $B^0_q$-$\bar{B}^0_q$ mixing phase.
 The values of $S$ and $A^{\Delta}$ are expected to be small  in the SM due to the suppression by the ratio $r\approx m_s/m_b$.
 Since   $S$ and $A^{\Delta}$ are approximately linearly dependent on  $r$, they are sensitive to even a small increase of right-handed photons.

 Currently, the observables $S$ and $A^{\Delta}$ are only weakly constrained.
 The ${\ensuremath{\B^0}}\xspace \ensuremath{\rightarrow}\xspace {\ensuremath{\kaon^0_{\mathrm{S}}}}\xspace {\ensuremath{\pion^0}}\xspace \gamma$ decay is a golden channel to study photon polarisation at $B$ factories.
The mixing-induced CP asymmetry in this channel has been measured to be
 $S_{{\ensuremath{\kaon^0_{\mathrm{S}}}}\xspace {\ensuremath{\pion^0}}\xspace \gamma} = -0.10\pm 0.31 \aunit{(stat)}\xspace \pm 0.07 \aunit{(syst)}\xspace$~\cite{Belle:2006pxp} and
 $S_{{\ensuremath{\kaon^0_{\mathrm{S}}}}\xspace {\ensuremath{\pion^0}}\xspace \gamma} = -0.78\pm 0.59 \aunit{(stat)}\xspace \pm 0.09 \aunit{(syst)}\xspace$~\cite{BaBar:2008okc} by the \mbox{Belle}\xspace and \mbox{BaBar}\xspace collaborations, respectively, both consistent with the SM expectation value of roughly $\frac{m_s}{m_b} \sin 2\beta  $.
 The coefficient of the $\sinh{\Delta\Gamma_d t/2}$ term,
 $A^{\Delta}$,
 is inaccessible  in  ${\ensuremath{\B^0}}\xspace$ decays  due to the tiny value of the ${\ensuremath{\B^0}}\xspace$ width difference, $\Delta\Gamma_d$.

Reconstruction of  ${\ensuremath{\B^0}}\xspace \ensuremath{\rightarrow}\xspace {\ensuremath{\kaon^0_{\mathrm{S}}}}\xspace {\ensuremath{\pion^0}}\xspace \gamma$ decays is challenging at the LHCb  experiment. Alternatively, LHCb can measure mixing-induced CP violation in $B^0\ensuremath{\rightarrow}\xspace {\ensuremath{\kaon^0_{\mathrm{S}}}}\xspace {\ensuremath{\pion^+}}\xspace {\ensuremath{\pion^-}}\xspace \gamma$ decays through a time-dependent amplitude analysis~\cite{LHCb-PII-Physics}. A more promising channel to probe right-handed NP
is the  decay ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace \phi \gamma$. Both  $S$ and   $A^{\Delta}$ can be measured in this channel and they are predicted to be close to zero in the SM~\cite{Muheim:2008vu}:
\begin{equation*}
S_{\phi \gamma} ({\rm SM}) = 0.000 \pm 0.002\;,\;
A^{\Delta}_{\phi \gamma } ({\rm SM}) =0.047\pm 0.039\;.
\end{equation*}
Using data collected in Run 1, the LHCb collaboration studied the tagged time-dependent decay rates of ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace \phi \gamma$, which are shown in Fig.~\ref{fig:phig}. The mixing-induced observables are measured to be~\cite{LHCb-PAPER-2019-015}
\begin{equation*}
 S_{\phi \gamma } =
0.43 \pm 0.30 \aunit{(stat)}\xspace \pm 0.11 \aunit{(syst)}\xspace\;,\;
A^{\Delta}_{\phi \gamma}=-0.67 ^{+0.37}_{-0.41} \aunit{(stat)}\xspace \pm 0.17 \aunit{(syst)}\xspace\;,
\end{equation*}
 which are  in agreement with the SM expectations.

\begin{figure}[!tb]
\begin{center}
\includegraphics[width=0.45\textwidth]{./figs/raredecay/phigamma1.pdf}
\includegraphics[width=0.45\textwidth]{./figs/raredecay/phigamma2.pdf}
\caption{ Decay-time distributions of the tagged (left) ${\ensuremath{\B^0_\squark}}\xspace\ensuremath{\rightarrow}\xspace \phi \gamma$ and (right) ${\ensuremath{\Bbar{}^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace \phi \gamma$ candidates, superimposed
by the fit projections. Plots are taken from Ref.~\cite{LHCb-PAPER-2019-015}.
}
\label{fig:phig}
\end{center}
\end{figure}

The  ${\ensuremath{\B^0}}\xspace\ensuremath{\rightarrow}\xspace {\ensuremath{\kaon^{*0}}}\xspace {\ensuremath{\Pe^+\Pe^-}}\xspace$ decay in the low-$q^2$ region  offers a powerful probe of right-handed new physics.
In the vicinity of the  photon pole, the decay amplitudes are dominated by contributions from the electromagnetic Wilson coefficients $C^{(\prime)}_7$.
An  angular analysis  can be performed in
a similar way as that in the ${\ensuremath{\B^0}}\xspace \ensuremath{\rightarrow}\xspace {\ensuremath{\kaon^{*0}}}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace$ case.
In the angular distribution, there are two terms arising from interference of the left-handed and right-handed decay amplitudes that  are proportional to $C_7$ and $C^{\prime}_7$, respectively.
For small values of $r$ ($\equiv|C_7^\prime/C_7)|$),
the coefficients of these two terms, denoted
 $A^{(2)}_{\rm T}$ and $A^{\rm Im}_{\rm T}$,
 are approximately expressed as~\cite{Becirevic:2011bp,Paul:2016urs}

\begin{equation}

   A^{(2)}_{\rm T}\approx r\,\cos (\phi_L-\phi_R)\;,\;
   A^{\rm Im}_{\rm T}\approx r\,\sin (\phi_L-\phi_R)\;,
\end{equation}
where $\phi_{L}$ and $\phi_R$ represent the phases of $C_7$ and $C^{\prime}_7$, respectively.
Like the  mixing-induced observables in ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace \phi \gamma$, $A^{(2)}_{\rm T}$ and $A^{\rm Im}_{\rm T}$ depend approximately linearly on $r$,
thus can provide high sensitivity  to right-handed currents in the small $r$ region.
Using data collected in Run 1 and Run 2, the LHCb collaboration has measured    $A^{(2)}_{\rm T}$ and $A^{\rm Im}_{\rm T}$  to be~\cite{LHCb-PAPER-2020-020}
\begin{equation*}
 A^{(2)}_{\rm T} =
0.11 \pm 0.10 \aunit{(stat)}\xspace \pm 0.02 \aunit{(syst)}\xspace\;,\;
A^{\rm Im}_{\rm T}=0.02 \pm 0.10 \aunit{(stat)}\xspace \pm 0.01 \aunit{(syst)}\xspace\;.
\end{equation*}
These results are compatible with the following SM predictions calculated using the FLAVIO software package~\cite{Straub:2018kue}:
\begin{equation*}
 A^{(2)}_{\rm T} ({\rm SM}) =
0.033 \pm 0.020 \;,\;
A^{\rm Im}_{\rm T} ({\rm SM})=-0.00012\pm0.00034\;,
\end{equation*}
and provide the most stringent constraint on the $b\ensuremath{\rightarrow}\xspace s \gamma$ photon polarisation.

Photon polarisation in $b \ensuremath{\rightarrow}\xspace s \gamma$ transitions
can also be probed by exploiting the angular correlations
in radiative decays of $b$ baryons or charged $b$ mesons.
 Since current detection technology cannot distinguish left-handed and right-handed photons, the final states with both left-handed and right-handed photons  are summed together. The left-handed  amplitude $A_L$ and right-handed amplitude $A_R$ add incoherently in the form of  $|A_L|^2+|A_R|^2$, without any interference.
 In certain  cases, the angular distributions allow for determining a parity violation parameter, $A_{\rm parity}$, which is proportional to the photon polarisation~\cite{Gronau:2002rz, Kou:2010kn},
 \begin{equation}
A_{\rm parity} \propto  \lambda_{\gamma}\equiv \frac{|A_L|^2-|A_R|^2}{|A_L|^2+|A_R|^2}\approx \frac{1-r^2}{1+r^2}\;,
\end{equation}
 This approach is powerful  in probing  large right-handed currents  but has limited sensitivity to any small  right-handed component.

 The LHCb collaboration observed a significantly non-zero up-down asymmetry of  the photons in $B^-\ensuremath{\rightarrow}\xspace K^-{\ensuremath{\pion^+}}\xspace {\ensuremath{\pion^-}}\xspace \gamma$ decays~\cite{LHCb-PAPER-2014-001} with respect to the plane defined by the three final-state hadrons in their rest frame.
 This observation demonstrates that the photons are indeed polarized.
 However, it is difficult to translate the measured asymmetry into a constraint on the polarisation parameter $\lambda_{\gamma}$, due to currently limited knowledge of the intermediate  resonances involved in this process.
 The baryonic decay  ${\ensuremath{\Lz^0_\bquark}}\xspace\ensuremath{\rightarrow}\xspace {\ensuremath{\PLambda}}\xspace \gamma$, observed by  the LHCb experiment using data collected in 2016~\cite{LHCb-PAPER-2019-010},
 provides a more convenient way  to measure the photon polarisation  in $b \ensuremath{\rightarrow}\xspace s \gamma$ transitions~\cite{Gremm:1995nx, Mannel:1997xy, Hiller:2001zj}.
 The angular distribution of this process is given by the differential rate
\begin{equation}
    \frac{d\Gamma}{d\cos\theta_p} \propto 1-\alpha_{{\ensuremath{\PLambda}}\xspace} \lambda_{\gamma} \cos \theta_p \;,
\end{equation}
where $\theta_p$  is the helicity angle of the proton in the ${\ensuremath{\PLambda}}\xspace$ rest frame with respect to the opposite direction of the photon, $\alpha_{{\ensuremath{\PLambda}}\xspace}$ is the   decay parameter of the weak process ${\ensuremath{\PLambda}}\xspace \ensuremath{\rightarrow}\xspace p {\ensuremath{\pion^-}}\xspace$. The photon polarisation parameter has recently been measured to be
$\lambda_{\gamma}=0.82^{+0.17}_{-0.26}{}^{+0.04}_{-0.13}$~\cite{LHCb-PAPER-2021-030} by the LHCb experiment using all data from Run 2 and the average of the decay parameter values of ${\ensuremath{\PLambda}}\xspace$  and
${\ensuremath{\offsetoverline{\PLambda}}}\xspace$ measured by BESIII, $\alpha_{{\ensuremath{\PLambda}}\xspace}=0.754\pm 0.004$~\cite{BESIII:2018cnd}. This result is in agreement with the SM prediction.

\clearpage

\section{CP violation in beauty and CKM parameters}
\label{sec:beautyCPV}

CP violation is a necessary condition to explain the matter dominated universe.

 While the SM with the CKM mechanism can account for the current experimental results on CP violation, it fails to explain  the cosmological matter–antimatter imbalance. Searching for new sources of CP violation  is one of the primary goals of flavour physics. This can be done by overconstraining the CKM matrix using measurements of the matrix  elements in many different processes.

The decays of $b$-hadrons provide a number of key measurements to access the five  CKM matrix elements  related to the $b$ or $t$ quark.
Taking advantages of the intense source of b-hadrons at the LHC and a detector designed to probe CP violation in heavy-flavour decays,  the LHCb experiment  has been the leading experiment in the field of $B$ physics  in the past ten years, and achieved  some of the most precise measurements of CP violation and mixing of $B$ mesons. Particularly, the  precision of the CKM angle $\gamma$ is now approaching that of the indirect determination;
the CP violation parameter $\phi_s$ and mixing parameter  ${\ensuremath{\Delta m_{{\ensuremath{\Ps}}\xspace}}}\xspace$ of the {\ensuremath{\B^0_\squark}}\xspace system, which are key observables for NP searches, have been pinned down with unprecedented precision.

This section describes the  key measurements of the CKM elements in the beauty sector by the LHCb experiment using the data taken in Run 1 and Run 2. The parameters $\gamma$, $\phi_s$ and ${\ensuremath{\Delta m_{{\ensuremath{\Ps}}\xspace}}}\xspace$, which have received the most significant improvements, are discussed in detail.
Other observables related to CKM global fit will also be discussed briefly while many other interesting topics, such as  CP violation in $b$-baryon decays, are not mentioned.

\subsection{CKM angle \texorpdfstring{$\gamma$}{TEXT}}
The angle $\gamma$, defined as $\arg[-(V_{ud}V_{ub}^*)/(V_{cd}V_{cb}^*)]$, is one of the key observables related to the CKM matrix. As can be seen from Fig.~\ref{fig:rhoeta_large}, one of the main limitations of global constraints comes from the angle $\gamma$. To be noted, the fits have already included recent $\gamma$ measurements from LHCb which improves the sensitivity on $\gamma$ from $14^{\circ}$, established at the era of $B$-factories, to around $5^{\circ}$. Future improvements on the sensitivities on $\gamma$ can be foreseen with the upgrade of LHCb and running of the Belle II experiment. In the following sections, we  briefly overview the main developments on $\gamma$ measurements in the past several years from the LHCb experiment.

The direct determination of the angle $\gamma$ is obtained through interference between $b \ensuremath{\rightarrow}\xspace c$ and $b \ensuremath{\rightarrow}\xspace u$ tree-level processes, where new physics hardly enters~\cite{Brod:2014bfa}.
The hadronic parameters of the system are all determined from experimental data and related theoretical uncertainty is negligible~\cite{Brod:2013sga}. The direct $\gamma$ measurements thus serves as a key input for SM predictions, which can be compared with other NP sensitive measurements to search for physics beyond the SM.

Several methods have been proposed to measure the angle $\gamma$, based on the types of $D$ decays. In this paper, when not specified, $D$ means a mixture state of $D^0$ and $\overline{D}^0$. The GLW method~\cite{Gronau:1991dp,Gronau:1990ra} refers to those decays with $D$ into a CP eigenstate or multi-body $D$ decays which can be effectively considered as a CP eigenstate using a CP-even fraction $F_+$.
The ADS method~\cite{Atwood:1996ci} refers to two-body $D$ decays or multi-body $D$ decays where the detailed structures over phase space are considered by introducing a global coherent factor $R_D$ and an effective strong phase $\delta_D$. In this case, the interference happens between $b\ensuremath{\rightarrow}\xspace c$ transition with $D$ decaying into doubly-Cabibbo-suppressed final states while $b\ensuremath{\rightarrow}\xspace u$ transition with $D$ decaying into Cabibbo favored final state.  The BPGGSZ method~\cite{Giri:2003ty,Bondar:2005ki,Bondar:2008hh} refers to $D$ decaying into multi-body final state where the phase space is binned to make full use of the statistic power of the decay. In addition, the angle $\gamma$ can also be extracted from time-dependent CP violation measurements of ${\ensuremath{\B^0_\squark}}\xspace$ decays. The measurements from the LHCb experiments with these methods are discussed in detail in the following section.

\subsubsection{GLW and ADS measurements}\label{sec:GLWADS}
The GLW and ADS channels are usually considered together when measured due to similar final states and experimental treatments. The GLW and ADS measurements have been performed in many decay channels to obtain the best sensitivity on $\gamma$. The full list of the measured channels can be found in Ref.~\cite{LHCb-PAPER-2021-033}.

The decay rate of the GLW method is
\begin{equation}
\label{EQ_GLW}
    \Gamma(B^- \ensuremath{\rightarrow}\xspace D[\ensuremath{\rightarrow}\xspace f_{\textrm{GLW}}]h^{-}) \propto 1+r_B^2  + 2\kappa_B r_B(2F_+ -1) \cos(\delta_B - \gamma),
\end{equation}
where $F_+ = 1 (0)$ means a pure CP-even (CP-odd) state while decays with $F_+ =0.5$ does not have any sensitivity on $\gamma$. The value of $F_+$ can be determined from quantum-coherent data from BESIII and CLEO-c experiments in a model-independent way. For example, the decays $B^- \ensuremath{\rightarrow}\xspace DK^-$, $D \ensuremath{\rightarrow}\xspace \pi^+ \pi^- \pi^0$ and $D \ensuremath{\rightarrow}\xspace K^+ K^- \pi^0$ have been used to measure the angle $\gamma$ in LHCb using the measured $F_+^{\pi^+\pi^-\pi^0}$ and $F_+^{K^+K^-\pi^0}$ of $0.973 \pm 0.017$ and $0.73 \pm 0.06$, respectively~\cite{Malde:2015mha}. One can see that almost full sensitivity can be achieved in the $\pi^+\pi^-\pi^0$ channel without considering structures over the phase space while further binning the phase space of $K^+K^-\pi^0$ will help get more sensitivity due to small $F_+^{K^+K^-\pi^0}$. In fact, the amplitude analyses of $D^0 \ensuremath{\rightarrow}\xspace \pi^+ \pi^- \pi^0$ and $D^0 \ensuremath{\rightarrow}\xspace K^+ K^- \pi^0$  channels have both been performed~\cite{BaBar:2007dro,CLEO:2005uoz,BaBar:2007soq,CLEO:2006tjv}. The relative strong phase between $D^0 \ensuremath{\rightarrow}\xspace \rho^+ \pi^-$ and $D^0 \ensuremath{\rightarrow}\xspace \rho^- \pi^+$ is $(-2.0 \pm 0.8)^{\circ}$ while that between $D^0 \ensuremath{\rightarrow}\xspace K^{*}(892)^+K^-$ and $D^0 \ensuremath{\rightarrow}\xspace K^{*}(892)^-K^+$ is $(-37.0 \pm 2.9)^{\circ}$, this is consistent with the fact $F_+^{\pi^+\pi^-\pi^0}$ is close to 1 and is much larger than $F_+^{K^+K^-\pi^0}$.
The parameters $r_B$ and $\delta_B$ are the amplitude ratio and phase difference between $b \ensuremath{\rightarrow}\xspace u$ and $b \ensuremath{\rightarrow}\xspace c$ processes, $\kappa_B$ is the coherent factor to take into account the sensitivity lost due to contamination of other contributions when the bachelor $h$ is a broad resonant structure, e.g. $K^*(892)$. The coherent factor $\kappa_B$ is obtained based on amplitude models used to describe multi-body $B$ decays. In the formula, only $B^-$ is written, however, it also applies for $B^0$ decays. The sensitivity of $\gamma$ is directly linked to  the size of $r_B$ and $\kappa_B$, larger $r_B \times \kappa_B$ gives better sensitivity on $\gamma$. The CP measurements from the GLW channels are only sensitive to $\cos(\delta_B - \gamma)$ which has four-fold ambiguity on determination of angle $\gamma$.

In the GLW modes, as $D$ decays into CP eigenstates, in case of no CP violation in $D$ decays, the decay amplitudes and their phases are the same for $D^0$ and $\overline{D}^0$. However, LHCb has discovered direct CP violation in $D^0 \ensuremath{\rightarrow}\xspace K^+ K^-$ and $D^0 \ensuremath{\rightarrow}\xspace \pi^+ \pi^-$ by looking at the difference of CP violation between the two decay channels~\cite{LHCb-PAPER-2019-006}. This indicates small difference of amplitudes and phases between the two channels, which affects $\gamma$ determination. Studies~\cite{Wang:2012ie} show that the effect on $\gamma$ determination is smaller than 0.5$^{\circ}$. The CP violation measured between $D^0 \ensuremath{\rightarrow}\xspace K^+ K^-$ and $D^0 \ensuremath{\rightarrow}\xspace \pi^+ \pi^-$ has been considered in the LHCb $\gamma$ combinations shown later.

Based on Eq.~\ref{EQ_GLW}, the angle $\gamma$ can be accessed by measuring the relative decay rate difference for $B^+$ and $B^-$ mesons. An example of the decay rate difference can be seen directly from the raw yields of the invariant mass distributions of the decays $B^{\pm} \ensuremath{\rightarrow}\xspace D(\ensuremath{\rightarrow}\xspace K^+K^-)K^{\pm}$ as shown in Fig.~\ref{fig:GLW}. The plots are from the latest LHCb measurements using $9\ensuremath{\fb^{-1}}\xspace$ data~\cite{LHCb-PAPER-2020-036}. By further considering the production asymmetry of $B^+$ and $B^-$ and detection efficiency difference between $K^+$ and $K^-$, the size of the CP violation of the decay can be determined.

\begin{figure}
\begin{center}
\includegraphics[width=0.45\textwidth]{./figs/beautyCPV/GLW_Bp.pdf}
\includegraphics[width=0.45\textwidth]{./figs/beautyCPV/GLW_Bm.pdf}
\caption{ The invariant mass distributions of $(K^+K^-)_DK^+$ (left) and $(K^+K^-)_DK^-$ (right) from \cite{LHCb-PAPER-2020-036}.
}
\label{fig:GLW}
\end{center}
\end{figure}

The decay rate of the ADS channel has a similar form of
\begin{equation}
 \Gamma(B^- \ensuremath{\rightarrow}\xspace D[\ensuremath{\rightarrow}\xspace f_{\textrm{ADS}}]h^-) \propto r_D^2+r_B^2  + 2\kappa_B r_B r_D R_D \cos(\delta_D + \delta_B - \gamma),
\end{equation}
where $r_D$ and $\delta_D$ are the average amplitude ratio and phase difference between doubly Cabibbo-suppressed and Cabibbo-favored $D$ decays. The value of $r_D$ is at similar magnitude as $r_B$ and thus leads to larger CP violation. However, the statistics of the ADS channel is also suppressed. This can be seen from the measurements done by the LHCb experiment using $B^{\pm} \ensuremath{\rightarrow}\xspace (K^{\mp}\pi^{\pm})K^{\pm}$ decays as show in Fig.~\ref{fig:ADS}~\cite{LHCb-PAPER-2020-036}.
\begin{figure}
\begin{center}
\includegraphics[width=0.45\textwidth]{./figs/beautyCPV/ADS_Bp.pdf}
\includegraphics[width=0.45\textwidth]{./figs/beautyCPV/ADS_Bm.pdf}
\caption{ The invariant mass distributions of $(K^+\pi^-)_D K^+$ (left) and $(\pi^+K^-)_D K^-$ (right) from \cite{LHCb-PAPER-2020-036}.
}
\label{fig:ADS}
\end{center}
\end{figure}
The measurements also have four-fold ambiguity, however, as $\delta_D$ is not zero, combining with GLW mode gives two-fold ambiguity on determination of angle $\gamma$.

The coherent factor $R_D$ equals to one for two-body final states and is less than one for multi-body final states to take into account dilutions due to different resonant contributions from Cabibbo-favored and doubly-Cabibbo-suppressed $D$ decays.
The values of $r_D$, $R_D$ and $\delta_D$ can also be determined using the quantum-coherent data from BESIII and CLEO-c experiments. For example, the coherent factor of ${\ensuremath{\D^0}}\xspace({\ensuremath{\Dbar{}^0}}\xspace) \ensuremath{\rightarrow}\xspace K^+\pi^+ \pi^- \pi^-$, $R_{K3\pi}$ is measured by the two experiments to be $0.43^{+0.17}_{-0.13}$~\cite{Evans:2016tlp}, which is significantly lower than one. Further binning the phase space according to the variation of strong phase determined from amplitude analysis can improve the sensitivity on $\gamma$~\cite{Evans:2019wza}.

The GLW and ADS measurements have been performed by the LHCb experiment with $D$ decaying mainly into charged final states. However, final states with neutral particles are also studied, either by reconstructing the neutral particles $\pi^0$ or $\gamma$, or in a partially reconstructed method, where $\pi^0$ and $\gamma$ from $D^{*0}$ are not reconstructed~\cite{LHCb-PAPER-2020-036}. In the second case, the reconstruction efficiency is much higher, however, the sensitivity of $\gamma$ is limited by the background contributions, especially those from $B^0 \ensuremath{\rightarrow}\xspace D^{*-}K^+$ decays which has similar line-shapes as signal.

\subsubsection{Measurements using the GGSZ method}
As has discussed in the above section, for multi-body $D$ decays, if CP even fraction $F_+$ for self-conjugated decay or coherent factor $R$ for semi-flavour tagged decay is significantly smaller than one, further sensitivity can be achieved by considering the variation of $r_D$ and $\delta_D$ over the phase space. One can model the $D$ decays with an amplitude model which provides the information of strong phase $\delta_D$ over the phase space. However, this may suffer from large systematic uncertainties due to modelling of amplitude distributions.

Alternative method is to bin the phase space and the effective $r_D^{i}$ and $\delta_D^{i}$ (or $c_i$ and $s_i$) in bin $i$ are defined as \begin{equation}
    c_i + i s_i \equiv R_D^{i} e^{i\delta_D^{i}}  \equiv \frac{\int_{\textrm{bin i}} A_f(p) \overline{A}_f(p)^* dp}{\sqrt{\int_{\textrm{bin i}} |A_f(p)|^2 dp} \sqrt{\int_{\textrm{bin i}} |\overline{A}_f(p)|^2dp}}.
\end{equation}
The GGSZ mode refers to decays into self-conjugated final states. BPGGSZ is used when referring to the analysis method using a binned strategy. In the BPGGSZ analysis, $c_i$ and $s_i$ are used instead of effective $R_D$ and $\delta_D$ for the benefit of better statistical performance.
The values of $c_i$ and $s_i$ can be determined using quantum coherent data collected by BESIII and CLEO-c experiments, where a mixture of $D^0$ and $\overline{D}^0$ can be achieved from $\psi(3770) \ensuremath{\rightarrow}\xspace D^0 \overline{D}^0$. As $\psi(3770)$ is a parity-odd state, the $D^0$ and $\overline{D}^0$ are in quantum-correlated state of $(|D^0 \overline{D}^0> - |\overline{D}^0 D^0>)/\sqrt{2}$.
Using the measured $c_i$ and $s_i$ values, the angle $\gamma$ and strong parameters $r_B$ and $\delta_B$ can be extracted. The method has been applied to $B^- \ensuremath{\rightarrow}\xspace DK^-$, $D \ensuremath{\rightarrow}\xspace {\ensuremath{\kaon^0_{\mathrm{S}}}}\xspace \pi^+ \pi^- $ and $D \ensuremath{\rightarrow}\xspace {\ensuremath{\kaon^0_{\mathrm{S}}}}\xspace K^+ K^-$ decays, where 16 bins and 4 bins are used according to the statistics of the decays, respectively. In order to optimise the sensitivity, binning schemes are chosen according to the strong phase variation over the Dalitz plot and possible background contamination and efficiency effects. Symmetry between ${\ensuremath{\kaon^0_{\mathrm{S}}}}\xspace \pi^+ ({\ensuremath{\kaon^0_{\mathrm{S}}}}\xspace K^+)$ and  ${\ensuremath{\kaon^0_{\mathrm{S}}}}\xspace \pi^- ({\ensuremath{\kaon^0_{\mathrm{S}}}}\xspace K^-)$ is also used to increase the sensitivity on $c_i$ and $s_i$. The binning scheme optimised for the ${\ensuremath{\kaon^0_{\mathrm{S}}}}\xspace \pi^+ \pi^-$ decays is shown in Fig.~\ref{fig:GGSZ}.

Each bin of the BPGGSZ method can offer constrains to the angle $\gamma$ and to the strong parameters $r_B$ and $\delta_B$. The number of measurements are much more than the number of unknown parameters, thus global production asymmetry and detection asymmetry can be treated as fit variables and related systematic uncertainties  are reduced. While for the ADS and GLW methods, production asymmetry and detection asymmetry have to be considered using control channels.
The measured yield difference in each bin is shown in Fig.~\ref{fig:GGSZ} for the ${\ensuremath{\kaon^0_{\mathrm{S}}}}\xspace\pi^+\pi^-$ and ${\ensuremath{\kaon^0_{\mathrm{S}}}}\xspace K^+K^-$ decays. Clear CP violation can be found. By combining the statistical power of these bins, the angle $\gamma$ is measured to be $(68.7^{+5.2}_{-5.1})^{\circ}$~\cite{LHCb-PAPER-2018-017}.
\begin{figure}
\begin{center}
\includegraphics[width=0.45\textwidth]{./figs/beautyCPV/Binning.pdf}
\includegraphics[width=0.45\textwidth]{./figs/beautyCPV/CPV_Bin.pdf}
\caption{
(Left) Binning scheme of $D \ensuremath{\rightarrow}\xspace {\ensuremath{\kaon^0_{\mathrm{S}}}}\xspace \pi^+ \pi^-$ decays used for $\gamma$ measurements; (right) Observed raw CP asymmetries and the predicted values using obtained CP parameters in different bins~\cite{LHCb-PAPER-2018-017}.
}
\label{fig:GGSZ}
\end{center}
\end{figure}

\subsubsection{Multi-body $B$ decays}
Similar to multi-body $D$ decays, multi-body $B$ decays can also be used where $r_B$ and $\delta_B$ now is a function of $B$ decay Dalitz plot. However, in this case, one can only use a model to describe different resonant contributions, where some can only be obtained through $b \ensuremath{\rightarrow}\xspace c$ process and some can be obtained through both $b \ensuremath{\rightarrow}\xspace c$ and $b \ensuremath{\rightarrow}\xspace u$ processes and the interference between them gives sensitivity to the angle $\gamma$. The measurements have been performed by the LHCb collaboration in $B^0 \ensuremath{\rightarrow}\xspace D K^+ \pi^-$ decays with $D \ensuremath{\rightarrow}\xspace K^+K^- (\pi^+\pi^-)$~\cite{LHCb-PAPER-2015-059} using 3 fb$^{-1}$ data. However, due to limited statistics, the sensitivity on $\gamma$ is still low. Further measurements with all the data collected by the LHCb experiment will be very interesting.

\subsubsection{Time-dependent ${\ensuremath{\B^0_\squark}}\xspace$ decays}
The angle $\gamma$ can be measured through time-dependent ${\ensuremath{\B^0_\squark}}\xspace$ and $B^0$ decays where the weak phases extracted are $(\gamma - 2\beta_s)$ and $(\gamma + 2\beta)$, respectively. In hadron collider experiments like LHCb, the golden channels are ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace D_s^{\pm} K^{\mp}$, ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace D_s^{\pm} K^{\mp} \pi^+\pi^-$ and ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace D \phi$ decays. The time-dependent analyses have been performed for the first two channels~\cite{LHCb-PAPER-2020-030,LHCb-PAPER-2017-047}, while only branching fraction has been measured for the ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace \overline{D}^0 \phi$ decay~\cite{LHCb-PAPER-2013-035,LHCb-PAPER-2018-015}.

As ${\ensuremath{\B^0_\squark}}\xspace$ mixing is involved, a time-dependent analysis is needed to extract CP parameters. The time-dependent decay rate of the ${\ensuremath{\B^0_\squark}}\xspace$ decay into a final state $f$ is given by
\begin{equation}
    \frac{d\Gamma_{{\ensuremath{\B^0_\squark}}\xspace\ensuremath{\rightarrow}\xspace f}(t)}{d t} \propto e^{-\Gamma_s t} [\cosh(\frac{\Delta \Gamma_s t}{2}) + A_f^{\Delta \Gamma} \sinh(\frac{\Delta \Gamma_s t}{2}) + C_f \cos(\Delta m_s t) - S_f\sin(\Delta m_s t)],
\end{equation}
where $\Delta \Gamma_s = \Gamma_{B_L} - \Gamma_{B_H}$ and $\Delta m_s = m_{B_H} - m_{B_L}$ are the decay-width and mass differences between the light ($B_L$) and heavy ($B_H$) ${\ensuremath{\B^0_\squark}}\xspace$ mass eigenstates and $\Gamma_s$ is the average ${\ensuremath{\B^0_\squark}}\xspace$ decay width. For the decays to the CP-conjugated final states, the CP violation parameters $C_f$, $S_f$ and $A_f^{\Delta \Gamma}$ are replaced with $C_{\overline{f}}$, $S_{\overline{f}}$ and
$A_{\overline{f}}^{\Delta \Gamma}$. These CP violation parameters are related with $\gamma$ through
\begin{eqnarray}\label{Eq:TD_Gamma}
    C_f &=& -C_{\overline{f}} =  \frac{1-r_B^2}{1+r_B^2}, \\\nonumber
    S_f &=& \frac{2r_B \sin(\delta_B-(\gamma-2\beta_s))}{1+ r_B^2}, S_{\overline{f}} = \frac{-2r_B \sin(\delta_B+(\gamma-2\beta_s))}{1+ r_B^2}, \\\nonumber
    A_f^{\Delta \Gamma} &=& \frac{-2r_B\cos(\delta_B - (\gamma-2\beta_s))}{1+ r_B^2}, A_{\overline{f}}^{\Delta \Gamma} = \frac{-2r_B\cos(\delta_B + (\gamma-2\beta_s))}{1+ r_B^2}
\end{eqnarray}
following the definitions in Ref.~\cite{LHCb-PAPER-2014-038}.
In the above formula, we have assumed no CP violation in either the mixing and in the decay amplitude.

The time-dependent CP violation has been measured using 3.0 fb$^{-1}$ and 9.0 fb$^{-1}$ data for ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace D_s^{\pm} K^{\mp}$ and ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace D_s^{\pm} K^{\mp} \pi^+ \pi^-$, respectively. The measured values of the angle $\gamma$, using the world-average value of $-2\beta_s$, are $\gamma = (128^{+17}_{-22})^{\circ}$ and $\gamma = (44 \pm 12)^{\circ}$, respectively (modulo 180$^{\circ}$). The ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace D_s^{\pm} K^{\mp} \pi^+ \pi^-$ is complicated due to multiple bachelor particles, and an amplitude analysis is needed. However, as discussed before, multi-body $B$ decays effectively introduce a dilution factor $\kappa$, as Eq.~\ref{Eq:TD_Gamma} have five constraints while together with $\kappa$, there are four unknown variables. The measurement could also be performed in a model-independent way, which leads to $\gamma = (44^{+20}_{-13})^{\circ}$. The sensitivity is worse than model-dependent results as expected.

Sensitivity studies with ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace D \phi$ using a time-integrated method has also been performed~\cite{Ao:2020cwh}, the expected statistical sensitivity of $\gamma$ is about $(8-19)^{\circ}$ using the \mbox{9\ensuremath{\fb^{-1}}\xspace} $pp$ collision data collected by the LHCb experiment. It points out that additional sensitivity on $\gamma$ can be achieved using the longitudinal polarized part of ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace D^* \phi$ decays with a partially reconstructed technique~\cite{LHCb-PAPER-2018-015}. Besides, further sensitivities on $\gamma$ can also be obtained by using other quasi-two-body decays in ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace D K^+K^-$ decays~\cite{Wang:2011zw}.

\subsubsection{Combination on $\gamma$}
\label{sec:04:01:05gammaCombination}
A combination of the parameter $\gamma$ was recently performed using all available results of $\gamma$ measurements from the LHCb experiment~\cite{LHCb-PAPER-2021-033}.
The determination of $\gamma$ relies on the inputs from charm decays; on the other hand, precise $\gamma$ measurements and strong parameters of $B$ decays can offer valuable constraints on the charm parameters, which, in turn, help constraining the mixing parameters of ${\ensuremath{\D^0}}\xspace$ mesons.
Therefore, the measurements that are sensitive to the charm mixing parameters are also used in the combination. The CKM angle $\gamma$ and charm mixing parameters are simultaneously determined with significant improvements.
Here we focus on the $\gamma$ combination results, leaving those on the charm mixing parameters to be discussed later in Section~\ref{sec:05:01:CharmMixing}.

\begin{figure}
\begin{center}
\includegraphics[width=0.6\textwidth]{./figs/beautyCPV/Comb_B.pdf}
\caption{$\gamma$ combinations from different $B$ hadrons~\cite{LHCb-PAPER-2021-033}.
}
\label{fig:Comb_B}
\end{center}
\end{figure}
Figure~\ref{fig:Comb_B} shows the $\gamma$ contributions from different $B$ decays.
One can immediately see that the sensitivity on $\gamma$ mainly comes from $B^+$ decays. The sensitivity on $\gamma$ is driven totally by statistics, yields of those channels from $B^0$ and ${\ensuremath{\B^0_\squark}}\xspace$ decays are much smaller than those from $B^+$ decays, thus offers much poorer sensitivity. On the other hand, one can also see that the central values of $\gamma$ determination from $B^0$ and ${\ensuremath{\B^0_\squark}}\xspace$ are around 20$^{\circ}$ higher, which motivates further measurements with the two $B$ mesons to check overall consistency between different $B$ mesons. Special attention may be needed on the ${\ensuremath{\B^0_\squark}}\xspace$ meson, where the ${\ensuremath{\B^0_\squark}}\xspace$ and ${\ensuremath{\B^0_\squark}}\xspace$ mixing is involved and new physics contributions can easily enter inside.

The 1~$\sigma$ contour of the constraint on $\gamma$ and strong parameters of $B^+$ decays are shown in Fig.~\ref{fig:Comb_rB_deltaB}. Ambiguities on determination of $\gamma$ from two-body $D$ decays can be seen as has been discussed in sec.~\ref{sec:GLWADS}, however, they constrain the parameter space into a very narrow region and together with the unique determination of $\gamma$ from the BPGGSZ channel, $\gamma$ can be precisely obtained.
The combined $\gamma$ is found to be $\gamma = (65.4^{+3.8}_{-4.2})^{\circ}$ and is the most precise determination from a single experiment.
\begin{figure}
\begin{center}
\includegraphics[width=0.45\textwidth]{./figs/beautyCPV/Comb_rB.pdf}
\includegraphics[width=0.45\textwidth]{./figs/beautyCPV/Comb_deltaB.pdf}
\caption{ $\gamma$ combinations from all LHCb measurements, 2D contours of $r_B$ vs $\gamma$ (left) and $\delta_B$ vs $\gamma$ have been shown at 68.3\% confidence level~\cite{LHCb-PAPER-2021-033}.
}
\label{fig:Comb_rB_deltaB}
\end{center}
\end{figure}

\subsection{CKM angle \texorpdfstring{$\beta$}{TEXT}}
The angle $\beta$, defined as arg$(-V_{cd}V_{cb}^*/V_{td}V_{tb}^*)$, is approximately the phase of $V_{td}^*$ in the Wolfenstein parameterisation~\cite{Wolfenstein:1983yz}.
 It enters the  decay time distributions  of {\ensuremath{\B^0}}\xspace  and ${\ensuremath{\Bbar{}^0}}\xspace$ meson decays   due to $B^0-{\ensuremath{\Bbar{}^0}}\xspace$ oscillation.
 The effective value of
 $\sin(2\beta)$, which could have been altered by NP contributions in $B^0-{\ensuremath{\Bbar{}^0}}\xspace$  mixing , can be extracted from the time-dependent CP asymmetries of  {\ensuremath{\B^0}}\xspace  decays  via $b \ensuremath{\rightarrow}\xspace {\ensuremath{\cquark\cquarkbar}}\xspace s$ transitions  following the relation
\begin{equation}
    A_{\textrm{CP}}(t) = \frac{\Gamma({\ensuremath{\Bbar{}^0}}\xspace(t)\ensuremath{\rightarrow}\xspace f_{\rm CP}) - \Gamma(B^0(t)\ensuremath{\rightarrow}\xspace f_{\rm CP})}
    {\Gamma({\ensuremath{\Bbar{}^0}}\xspace(t)\ensuremath{\rightarrow}\xspace f_{\rm CP}) + \Gamma(B^0(t)\ensuremath{\rightarrow}\xspace f_{\rm CP})} \approx \eta_{f} \sin(2\beta)\sin(\Delta mt)\;,
\end{equation}
where $f_{\rm CP}$ is a CP eigenstate with eigenvalue $\eta_f$,
and the approximation assumes no  CP violation in the mixing or decay. Any significant deviation of the measured $\sin(2\beta)$ value from the indirect determination of $\sin(2\beta)$ through a global CKM fit excluding  $\sin(2\beta)$  measurements is a clear sign of NP.

The ${\ensuremath{\cquark\cquarkbar}}\xspace$ pair could appear either
in a charmonium meson or in two charmed mesons
in the final state.
The precision of the $\beta$  is mainly driven by {\ensuremath{\B^0}}\xspace decays to charmonium final-states due to their large decay rates and distinct characteristics for identification.
The world-average of $\sin(2\beta)$ of all the charmonium measurements is  sin$(2\beta)=0.699 \pm 0.017$~\cite{HFLAV18}. The LHCb experiment has performed measurements of $\sin(2\beta)$
in the decays ${\ensuremath{\B^0}}\xspace \ensuremath{\rightarrow}\xspace {\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace(\ensuremath{\rightarrow}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace) K^0_S$~\cite{LHCb:2015ups}, ${\ensuremath{\B^0}}\xspace \ensuremath{\rightarrow}\xspace {\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace(\ensuremath{\rightarrow}\xspace e^+e^-) K^0_S$  and  ${\ensuremath{\B^0}}\xspace \ensuremath{\rightarrow}\xspace \psi(2S)(\ensuremath{\rightarrow}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace) K^0_S$~\cite{LHCb:2017mpa}. The combined  value is $\sin(2\beta)=0.760\pm 0.034$,
the precision of which is already comparable to that of the  \mbox{BaBar}\xspace result  $\sin(2\beta)=0.09\pm 0.03 \pm 0.01$~\cite{BaBar:2009byl} and the Belle result $\sin(2\beta)=0.67\pm 0.02 \pm 0.01$ from Belle~\cite{Belle:2012paq},
An improvement  by a factor of two is expected from measurements including  LHCb Run 2 data.

The presence of small penguin contributions in $b \ensuremath{\rightarrow}\xspace {\ensuremath{\cquark\cquarkbar}}\xspace s$ processes  may shift the measured values of $\sin(2\beta)$ by up to few percents~\cite{Ciuchini:2005mg}. The decays via tree-level $b \ensuremath{\rightarrow}\xspace {\ensuremath{\Pc}}\xspace{\ensuremath{\overline \uquark}}\xspace{\ensuremath{\Pd}}\xspace$ transitions, though having smaller signal yields due to the small branching fractions of $D$ decays, are free of the  penguin effects and thus theoretically clean.
A recent joint analysis  of the decay
$\overline{B}^0 \ensuremath{\rightarrow}\xspace Dh^0$ with  $D\ensuremath{\rightarrow}\xspace {\ensuremath{\kaon^0_{\mathrm{S}}}}\xspace h^+ h^-$ by  the \mbox{BaBar}\xspace and Belle  experiments
measured $\sin(2\beta)=0.80 \pm 0.14\pm 0.04\pm 0.03$ and $\cos(2\beta)=0.91 \pm 0.22 \pm 0.09 \pm 0.07$, which ruled out the other solution of $\beta$ at 7.3$\sigma$~\cite{BaBar:2018agf,BaBar:2018cka}.
Analysis of this decay is challenging at LHCb due to the presence of  ${\ensuremath{\kaon^0_{\mathrm{S}}}}\xspace$ and $\pi^0$ mesons in the final state. On the other hand,
the decay $\overline{B}^0 \ensuremath{\rightarrow}\xspace D\pi^+\pi^-$ followed by $D\ensuremath{\rightarrow}\xspace K^+K^- (\pi^+\pi^-)$ only involves charged particle  and thus is ideal for LHCb to pursue.

\subsection{CKM angle \texorpdfstring{$\beta_s$}{TEXT}}
 The angle $\beta_s$, defined as $-$arg$(-V_{cb}V_{cs}^*/V_{tb}V_{ts}^*)$, is approximately the phase of $V_{ts}$ in the Wolfenstein parameterisation.
 The effective value of $-2\beta_s$  can be measured in the time-dependent CP asymmetries of ${\ensuremath{\B^0_\squark}}\xspace$ decays to  CP eigenstates via $b \ensuremath{\rightarrow}\xspace {\ensuremath{\cquark\cquarkbar}}\xspace s$ transitions, and is  denoted $\phi_s^{{\ensuremath{\cquark\cquarkbar}}\xspace s}$.
In contract to the angle $\beta$, $\beta_s$ is very small. The SM prediction for $\phi_s $ is $\phi_s^{\rm SM} = -2\beta_s=-0.03696 \pm 0.0004$~\cite{CKMfitter2015}, which is
subject to small corrections due to the neglected penguin contributions in  $b \ensuremath{\rightarrow}\xspace {\ensuremath{\cquark\cquarkbar}}\xspace s$ decays.
Presence of new particles  in ${\ensuremath{\B^0_\squark}}\xspace$-${\ensuremath{\Bbar{}^0_\squark}}\xspace$ mixing diagrams may have a sizeable effect on  $\phi_s$, making it an sensitive probe of physics beyond the SM.

The  LHCb experiment has performed measurements of $\phi_s^{{\ensuremath{\cquark\cquarkbar}}\xspace s}$ in the decays ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace {\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace(\ensuremath{\rightarrow}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace) \phi$~\cite{LHCb-PAPER-2019-013}, ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace {\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace(\ensuremath{\rightarrow}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace) \pi^+ \pi^-$~\cite{LHCb-PAPER-2019-003}, ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace {\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace(\ensuremath{\rightarrow}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace) K^+ K^-$ with $m(K^+K^-)>1.05\ensuremath{\aunit{Ge\kern -0.1em V\!/}c^2}\xspace$~\cite{LHCb-PAPER-2017-008}, ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace {\ensuremath{\Ppsi{(2S)}}}\xspace(\ensuremath{\rightarrow}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace)\phi$~\cite{LHCb-PAPER-2016-027}, ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace D_s^+ D_s^-$~\cite{LHCb-PAPER-2014-051} and more recently in ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace {\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace(\ensuremath{\rightarrow}\xspace {\ensuremath{\Pe^+\Pe^-}}\xspace) \phi$~\cite{LHCb-PAPER-2020-042}. A combination of the measurements  in ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace {\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace(\ensuremath{\rightarrow}\xspace \mu^+\mu^-) \phi$ and  ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace {\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace(\ensuremath{\rightarrow}\xspace \mu^+\mu^-) \pi^+ \pi^-$ obtained using data taken during 2011-2016 gives $\phi_s^{{\ensuremath{\cquark\cquarkbar}}\xspace s}= -0.042 \pm 0.025$ rad.

The average of all LHCb measurements of  $\phi_s^{{\ensuremath{\cquark\cquarkbar}}\xspace s}$
is compared with the ATLAS~\cite{ATLAS:2020lbz} and  CMS~\cite{CMS:2020efq} results in Fig.~\ref{fig:PhisDGs}, where 2-dimensional contours in the plane of $\phi_s^{{\ensuremath{\cquark\cquarkbar}}\xspace s}$ versus the {\ensuremath{\B^0_\squark}}\xspace decay width difference ($\Delta \Gamma_s$) are displayed at 68\% confidence level~\cite{HFLAV18}.
Note ATLAS and CMS have only performed measurements in ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace {\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace(\ensuremath{\rightarrow}\xspace {\ensuremath{\Pmu^+\Pmu^-}}\xspace) \phi$ decays, due to constraints from their trigger systems,.
The current world-average value is $\phi_s^{{\ensuremath{\cquark\cquarkbar}}\xspace s}= -0.050 \pm 0.019$~rad~\cite{HFLAV18}.
In Fig~\ref{fig:PhisDGs}, one can see a good agreement in the $\phi_s^{{\ensuremath{\cquark\cquarkbar}}\xspace s}$ measurements
from different measurements. However, some tension is observed for  $\Delta \Gamma_s$ and $\Gamma_s$. Factors of 2.5 and 1.77 have been applied to scale up the uncertainties of $\Gamma_s$ and $\Delta \Gamma_s$ in the combination. Further investigations by the relevant experiments are needed to solve this problem.
With the uncertainty of $\phi_s^{{\ensuremath{\cquark\cquarkbar}}\xspace s}$
well below its SM value,
the study of CP violation in {\ensuremath{\B^0_\squark}}\xspace decays enter an era of precision test, and control of the penguin pollution in $\phi_s^{{\ensuremath{\cquark\cquarkbar}}\xspace s}$ using data-driven methods is essential for identification of NP signals in  ${\ensuremath{\B^0_\squark}}\xspace-{\ensuremath{\Bbar{}^0_\squark}}\xspace$ mixing~\cite{Liu:2013nea, Faller:2008zc, Nagahiro:2008cv, DeBruyn:2010hh, Jung:2012mp, DeBruyn:2014oga, Frings:2015eva, Barel:2020jvf}.
In this regard, the LHCb experiment has measured CP violation in the penguin-enhanced $b\ensuremath{\rightarrow}\xspace {\ensuremath{\cquark\cquarkbar}}\xspace d$ decays ${\ensuremath{\B^0}}\xspace \ensuremath{\rightarrow}\xspace {\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace \rho^0$~\cite{LHCb-PAPER-2014-058} and ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace {\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace {\ensuremath{\Kbar{}^{*0}}}\xspace$~\cite{LHCb-PAPER-2015-034} decays. The measurements  are used to estimate the penguin shift of $\phi_s^{{\ensuremath{\cquark\cquarkbar}}\xspace s}$ measured in ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace {\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace \phi$,  assuming SU(3) flavour symmetry. The shift is found to be compatible  with zero~\cite{LHCb-PAPER-2015-034}, with a uncertain well below the statistical uncertainties  of the current $\phi_s^{{\ensuremath{\cquark\cquarkbar}}\xspace s}$ measurements.

\begin{figure}
\centering
\includegraphics[width=0.45\textwidth]{./figs/beautyCPV/PDG2021_Phis_vs_DGs.pdf}
\includegraphics[width=0.46\textwidth]{./figs/beautyCPV/PDG2021_Gs_vs_DGs.pdf}
\caption{ Combination  of $\phi_s^{{\ensuremath{\cquark\cquarkbar}}\xspace s}$ measurements by HFLAV~\cite{HFLAV18}, where 2D contours of (left) $\phi_s^{{\ensuremath{\cquark\cquarkbar}}\xspace s}$ vs $\Delta\Gamma_s$ and (right) $\Gamma_s$ vs $\Delta \Gamma_s$ are displayed at 68\% confidence level.
}\label{fig:PhisDGs}
\end{figure}

Similar CP-violating phases can also be measured in  decays of {\ensuremath{\B^0_\squark}}\xspace mesons via
$b \ensuremath{\rightarrow}\xspace {\ensuremath{\squark\squarkbar}}\xspace s$ and $b \ensuremath{\rightarrow}\xspace {\ensuremath{\dquark\dquarkbar}}\xspace s$ transitions,  denoted $\phi_s^{{\ensuremath{\squark\squarkbar}}\xspace s}$ and $\phi_s^{{\ensuremath{\dquark\dquarkbar}}\xspace s}$, respectively.
Since  these decays are dominated by penguin diagrams with internal top quarks,
the phases $\phi_s^{{\ensuremath{\squark\squarkbar}}\xspace s}$ and $\phi_s^{{\ensuremath{\dquark\dquarkbar}}\xspace s}$ receive  contributions from the decay amplitudes that cancels out the contribution of $-2\beta$  from the ${\ensuremath{\B^0_\squark}}\xspace - {\ensuremath{\Bbar{}^0_\squark}}\xspace$ mixing, resulting to  vanishing predictions in the SM. Measurements of these quantities can  probe NP in these FCNC decays.
 LHCb has measured  $\phi_s^{{\ensuremath{\squark\squarkbar}}\xspace s}$ and $\phi_s^{{\ensuremath{\dquark\dquarkbar}}\xspace s}$ in the ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace \phi \phi$ and ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace {\ensuremath{\kaon^*}}\xspace(892){\ensuremath{\kaon^*}}\xspace(892)$ decays using data taken during 2011-2016, and the results are
  $\phi_s^{{\ensuremath{\squark\squarkbar}}\xspace s} = -0.073 \pm 0.115 \pm 0.027$~rad~\cite{LHCb-PAPER-2019-019} and $\phi_s^{{\ensuremath{\dquark\dquarkbar}}\xspace s} = -0.10 \pm 0.13 \pm 0.14$~rad~\cite{LHCb-PAPER-2017-048}, respectively, where the first uncertainties are statistical and the second systematic.

\subsection{CKM elements \texorpdfstring{$V_{ub}$}{TEXT} and \texorpdfstring{$V_{cb}$}{TEXT}}
The amplitudes of the CKM matrix elements $V_{ub}$ and $V_{cb}$ are measured through semi-leptonic transitions of $b \ensuremath{\rightarrow}\xspace u \ell\nu$ and $b \ensuremath{\rightarrow}\xspace c \ell \nu$. They have been extensively studied previously in $B$-factories using the so-called exclusive and inclusive methods, which infer to whether a specific decay channel is used or not. The two approaches suffer from different theoretical and experimental uncertainties and offer important cross-checks between each other. In the inclusive measurements, the Heavy Quark Expansion (HQE) is used as $\Lambda_{\textrm{QCD}}/m_b$ is small and the Operator Product Expansion (OPE) calculates non-perturbative contributions involved.
For the exclusive measurements, parameterisation of the form factor of the corresponding decay is needed, where inputs are obtained from light-cone sum rules (LCSR) or from lattice QCD (LQCD). However, tensions have been found between the inclusive and exclusive results as can be seen in Fig.~\ref{fig:VubVcb}~\cite{CKMfitter2005}.
\begin{figure}
\centering
\includegraphics[width=0.8\textwidth]{./figs/beautyCPV/VubVcb.pdf}
\caption{ Global fits of $|V_{ub}|$ and $|V_{cb}|$ from the CKMfitter group~\cite{CKMfitter2005} where measurements using different methods have been displayed.
}\label{fig:VubVcb}
\end{figure}
Efforts from both experimental and theoretical sides are needed to understand the discrepancy.

Unlike $B$-factories, where full kinematics can be obtained, the LHCb can not obtain kinematic information of neutrinos from energy-momentum conservation due to many tracks produced in  $pp$ collisions. In addition, background contributions from other $b$ and $c$ hadrons, and also huge combinatorial backgrounds randomly combined from tracks other than signal makes the analyses of semi-leptonic decays very complicated. Despite of these difficulties, the LHCb experiment has successfully measured the ratio of $|V_{ub}|$ and $|V_{cb}|$ using ${\ensuremath{\Lz^0_\bquark}}\xspace \ensuremath{\rightarrow}\xspace p \mun {\ensuremath{\neu_\mu}}\xspace$ ($q^2 > 15\gevgevcccc$) and ${\ensuremath{\Lz^0_\bquark}}\xspace \ensuremath{\rightarrow}\xspace {\ensuremath{\Lz^+_\cquark}}\xspace \mu^- {\ensuremath{\neu_\mu}}\xspace$ ($q^2 > 7\gevgevcccc$)~\cite{LHCb-PAPER-2015-013}. The usage of the control channel ${\ensuremath{\Lz^0_\bquark}}\xspace \ensuremath{\rightarrow}\xspace {\ensuremath{\Lz^+_\cquark}}\xspace \mu^- {\ensuremath{\neu_\mu}}\xspace$ not only cancels out common systematic uncertainties between the two channels, but also offers a global scale needed to determine $|V_{ub}|$ from branching fractions. Using the updated branching fraction measurement of ${\ensuremath{\Lz^+_\cquark}}\xspace \ensuremath{\rightarrow}\xspace pK^-\pi^+$, benefiting from the ${\ensuremath{\Lz^+_\cquark}}\xspace{\ensuremath{\Lbar{}^-_\cquark}}\xspace$ data collected by the BESIII experiment~\cite{BESIII:2015bjk}, the ratio is determined to be $|V_{ub}|/|V_{cb}| = 0.079 \pm 0.009$.

Using the same approach, the ratio of $|V_{ub}|$ and $|V_{cb}|$ is also determined using ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace K^- \mu^+ \nu_{\mu}$ and ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace D_s^- \mu^+ {\ensuremath{\neu_\mu}}\xspace$~\cite{LHCb-PAPER-2020-038}. In the measurement, two $q^2$ regions are used, \mbox{$q^2 > 7\gevgevcccc$} and $q^2 < 7\gevgevcccc$ where the form factors are obtained from LQCD and LCSR, respectively. However, the ratios of $|V_{ub}|$ and $|V_{cb}|$ obtained from the two methods differ significantly, $|V_{ub}|/|V_{cb}|_{q^2 < 7\gevgevcccc} = 0.0607 \pm 0.0015 \pm 0.0013 \pm 0.0008 \pm 0.0030$ and $|V_{ub}|/|V_{cb}|_{q^2 > 7\gevgevcccc} = 0.0946 \pm 0.0030^{+0.0024}_{-0.0025} \pm 0.00013 \pm 0.0068$, where the first uncertainty of each result is statistical, the second systematic, the third due to $D_s^-$ branching fraction, and the last one from the form factor.  The discrepancy between the two
results clearly indicates that more efforts from theoretical side are needed to resolve the tension on $|V_{ub}|$ and $|V_{cb}|$ measurements.

In addition to the ratio between the two CKM matrix elements, the LHCb experiment is also exploring its potential in determine the $|V_{cb}|$ alone using ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace D_s^{(*)-}\mu^+ \nu_{\mu}$~\cite{LHCb-PAPER-2019-041}. The branching fraction of ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace D_s^{(*)-}\mu^+ \nu_{\mu}$ is needed to set the global scale for the determination of $|V_{cb}|$. This is obtained using the control channel $B^0 \ensuremath{\rightarrow}\xspace D^{(*)-}\mu \nu_{\mu}$, where the ratios of
\begin{eqnarray}
    R &=& \frac{{\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace D_s^{-}\mu^+ \nu_{\mu}}{B^0 \ensuremath{\rightarrow}\xspace D^{-}\mu \nu_{\mu}}, \\
    R^* &=& \frac{{\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace D_s^{*-}\mu^+ \nu_{\mu}}{B^0 \ensuremath{\rightarrow}\xspace D^{*-}\mu \nu_{\mu}},
\end{eqnarray}
are determined.
The data from LHCb offers $q^2$ dependence needed to extract $|V_{cb}|$ together with non-perturbative inputs.
Using form factor parameterisations from Caprini, Lellouch and Neubert~\cite{Caprini:1997mu} or from Boyd, Grinstein and Lebed~\cite{Boyd:1994tt,Boyd:1997kz}, the measured values of $|V_{cb}|$ are $(41.4 \pm 0.6 \pm 0.9 \pm 1.2) \times 10^{-3}$ and $(42.3 \pm 0.8 \pm 0.9 \pm 1.2) \times 10^{-3}$, respectively, where the first uncertainty is statistical, the second systematic, and the last one due to external inputs such as theoretical inputs on form factors, branching fractions of $D_s^-$ or $D^-$ decays, ${\ensuremath{\B^0_\squark}}\xspace$ lifetime etc.

\subsection{\texorpdfstring{$\Delta m_d$}{TEXT} and \texorpdfstring{$\Delta m_s$}{TEXT}}
The parameters $\Delta m_d$ and $\Delta m_s$
denote the mass differences between the heavy and light mass eigenstate of the ${\ensuremath{\B^0}}\xspace$ and {\ensuremath{\B^0_\squark}}\xspace systems, and define
 the oscillation frequencies of $B^0$ mixing and ${\ensuremath{\B^0_\squark}}\xspace$ mixing, respectively.
Currently, The  most precise determination of $\Delta m_d$ comes from the LHCb measurements in semileptonic decays with a $D^-$ or $D^{*-}$ meson using 3 fb$^{-1}$ of data. Combining the results obtained in the two decay modes yields $\Delta m_d = (0.5050 \pm 0.0021 \pm 0.0010)$ ps$^{-1}$~\cite{LHCb-PAPER-2015-031}, where the first uncertainty is statistical and the second systematic.
The world-average is
$\Delta m_d = 0.5065 \pm 0.0019$ ps$^{-1}$~\cite{HFLAV18}.
The determination of  $\Delta m_s$ is also led by the LHCb experiment. Combining the recent
 measurements in the decays  ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace D_s^- \pi^+$~\cite{LHCb-PAPER-2021-005} and ${\ensuremath{\B^0_\squark}}\xspace \ensuremath{\rightarrow}\xspace D_s^- \pi^+ \pi^- \pi^+$~\cite{LHCb-PAPER-2020-030} and earlier measurements  yields   $\Delta m_s =17.7656 \pm 0.0057$ ps$^{-1}$.

However, the constraints on the CKM matrix elements provided by the $\Delta m_d$ and $\Delta m_s$ measurements rely on  the decay constants and  Bag parameters of the the $B^0$ and ${\ensuremath{\B^0_\squark}}\xspace$ mesons, which are obtained from Lattice QCD calculations~\cite{FlavourLatticeAveragingGroup:2019iem}. The precision of these hadronic parameters
is much worse than the experimental precision, thus limits the constraining  power of $\Delta m_d$ and $\Delta m_s$ on the CKM  global fit. Further improvements in Lattice QCD calculations
are eagerly awaited.

\subsection{Global fit}

 The four parameters of the
 CKM  matrix are measured  in different processes and some of the key observables have been discussed in the above sections.  A global fit is needed to get the best sensitivity and to probe NP effects.
The $\chi^2$ value of the fit provides a measure of
the overall consistency between the different measurements,
while  the pull value for each measurement quantifies the difference between the measured value and  value predicted by the fit results.
Clues for new physics can be identified from large $\chi^2$ or pull values.

The constraining of the CKM matrix is usually illustrated in  complex planes using  triangles defined using unitarity relations of the CKM matrix. The most commonly quoted CKM triangle corresponds to the relation ${\ensuremath{V_{\uquark\dquark}^{\phantom{\ast}}}}\xspace{\ensuremath{V_{\uquark\bquark}^\ast}}\xspace + {\ensuremath{V_{\cquark\dquark}^{\phantom{\ast}}}}\xspace{\ensuremath{V_{\cquark\bquark}^\ast}}\xspace + {\ensuremath{V_{\tquark\dquark}^{\phantom{\ast}}}}\xspace{\ensuremath{V_{\tquark\bquark}^\ast}}\xspace =0$, where $V_{ij}$ is the CKM matrix element
between the quarks of the flavours $i$ and $j$.
The results of the state-of-the-art global fit  are shown in Fig.~\ref{fig:rhoeta_large} provided by the CKMfitter group~\cite{CKMfitter2015}.
\begin{figure}
\begin{center}
\includegraphics[width=0.5\textwidth]{./figs/beautyCPV/rhoeta_large.pdf}
\caption{Global fit results of the CKM matrix from different measurements, provided by the CKMfitter group~\cite{CKMfitter2005}. Different colours and labels on the plot indicate constraints from measurements of different observables.
}
\label{fig:rhoeta_large}
\end{center}
\end{figure}

Within current precision, different measurements cross on a single point and give an overall consistent picture. However,  the argument of matter-antimatter asymmetry suggests  that CP violation  from sources beyond the CKM matrix may break the consistency. One of the main efforts of the LHCb experiment is to search for such a discrepancy by further improving measurement precision in the beauty and charm sectors.

In the following decade, both the LHCb  and Belle II experiments will accumulate much more data to further constrain the CKM matrix. With data collected till 2015, either experiment will be able to reduce the uncertainty of $\gamma$ to around $1.5\ensuremath{^{\circ}}\xspace$, and  further improve it to $0.3\ensuremath{^{\circ}}\xspace$ after 2030s. Other CKM angles and matrix elements will also be significantly improved.
Details on the future outlook
can be found in section~\ref{sec:upgrade-prospect}. Together with improvements of other measurements and lattice calculations, NP may be observed from inconsistency between different measurements.

\clearpage
\section{Charm mixing and CP violation}
\label{sec:charmCPV}

Charm physics covers the studies of hadrons containing charm quarks.
CP violation in the charm sector is expected to be incredibly small in the SM, of the order $\mathcal{O}(10^{-3})$ or less~\cite{PDG2020}.
However, the presence of new physics may enhance the amount of CP violation,
which can be probed using  the enormously large sample of charmed hadrons at LHCb. Particularly, the study of mixing and CP violation of neutral $D$ mesons can provide unique probes of  NP in FCNC transitions in the up-type quark sector, complementary to the study of  mixing and CP violation in neutral $B$ and $K$ mesons, which are sensitive to NP in FCNC transitions of down-type quarks. For this reason,
this section mainly focuses on  results in  mixing and CP violation of {\ensuremath{\D^0}}\xspace mesons from the LHCb experiment.

\subsection{Neutral {\ensuremath{\PD}}\xspace meson mixing}
\label{sec:05:01:CharmMixing}
Similar to  neutral ${\ensuremath{\kaon^0}}\xspace$ and ${\ensuremath{\B_{(\squark)}^0}}\xspace$ mesons, the neutral charmed meson, ${\ensuremath{\D^0}}\xspace$,
can oscillate to its antiparticle partner, ${\ensuremath{\Dbar{}^0}}\xspace$, via the short-distance ${\ensuremath{\PW^\pm}}\xspace$ exchange or long-distance rescattering diagrams, as shown in Fig.~\ref{fig:charm:00}.
This phenomenon of oscillation or mixing can be characterised by
the normalised (dimensionless) mixing parameters $x$ and $y$, defined as
\begin{equation}
x\equiv \frac{\Delta M}{\Gamma},\quad y\equiv \frac{\Delta \Gamma}{\Gamma},
\end{equation}
where $\Delta M$ ($\Delta\Gamma$) is the mass (decay width) difference of the heavy and light mass eigenstates, and $\Gamma$ is the average decay width.
Unlike in the case of the ${\ensuremath{\kaon^0}}\xspace$ or ${\ensuremath{\B_{(\squark)}^0}}\xspace$ system, both  $x$ and $y$ in the {\ensuremath{\D^0}}\xspace system are significantly smaller than unity,
thus  very large data samples are required to observe {\ensuremath{\D^0}}\xspace mixing and determine the  tiny values of $x$ and $y$.

\begin{figure}[!b]
\begin{center}
\includegraphics[width=0.49\textwidth]{./figs/charmCPV/FeynmanBox.pdf}
\includegraphics[width=0.49\textwidth]{./figs/charmCPV/FeynmanBox2.pdf}\\
\includegraphics[width=0.49\textwidth]{./figs/charmCPV/FeynmanLongQCD.pdf}
\end{center}
\caption{Oscillation of ${\ensuremath{\D^0}}\xspace$ mesons via (top) ${\ensuremath{\PW^\pm}}\xspace$ exchanges or (bottom) rescattering effect.
}
\label{fig:charm:00}
\end{figure}

Evidence of ${\ensuremath{\D^0}}\xspace$-${\ensuremath{\Dbar{}^0}}\xspace$ mixing was first reported by \mbox{BaBar}\xspace~\cite{BaBar:2007kib} and \mbox{Belle}\xspace~\cite{BELLE:2007dgh} in 2007, and later also seen by \mbox{CDF}\xspace~\cite{CDF:2007bdz} in 2008.
Subsequent measurements by \mbox{BaBar}\xspace~\cite{BaBar:2008xkf,BaBar:2009kmh} with different {\ensuremath{\D^0}}\xspace decay channels provided more evidences of the mixing.
The combination of these measurements confirmed the existence of charm mixing with a significance more than $5\sigma$.
The first observation of ${\ensuremath{\D^0}}\xspace$-${\ensuremath{\Dbar{}^0}}\xspace$ mixing in a single measurement was achieved by LHCb~\cite{LHCb-PAPER-2012-038} in 2012 by using the data taken in 2011 to study the time-dependent ratio of ${\ensuremath{\D^0}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\kaon^+}}\xspace{\ensuremath{\pion^-}}\xspace$ (doubly Cabibbo-suppressed, DCS) to ${\ensuremath{\D^0}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\kaon^-}}\xspace{\ensuremath{\pion^+}}\xspace$ (Cabibbo favoured, CF) decay rates.
The {\ensuremath{\D^0}}\xspace candidates are selected from the  ${\ensuremath{\D^{*+}}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\D^0}}\xspace{\ensuremath{\pion^+}}\xspace$ decays, where the charge of the pion directly from each {\ensuremath{\D^{*+}}}\xspace decay is used to determine the {\ensuremath{\D^0}}\xspace flavour at its production time.

The ${\ensuremath{\D^{*+}}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\D^0}}\xspace(\ensuremath{\rightarrow}\xspace{\ensuremath{\kaon^-}}\xspace{\ensuremath{\pion^+}}\xspace){\ensuremath{\pion^+}}\xspace$ process, referred to as the right-sign (RS) process,
is dominated by a CF decay, contaminated with a small contribution from the {\ensuremath{\D^0}}\xspace-{\ensuremath{\Dbar{}^0}}\xspace mixing followed by the DCS decay;
the ${\ensuremath{\D^{*+}}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\D^0}}\xspace(\ensuremath{\rightarrow}\xspace{\ensuremath{\kaon^+}}\xspace{\ensuremath{\pion^-}}\xspace){\ensuremath{\pion^+}}\xspace$ process, referred to as wrong-sign (WS) process,
includes contributions from both the DCS decay and the {\ensuremath{\D^0}}\xspace-{\ensuremath{\Dbar{}^0}}\xspace mixing followed by the CF decay.
Under the assumption of small mixing and negligible CP violation, the time-dependent ratio of the WS to the RS decay rates, $R$, is given by~\cite{Bianco:2003vb}
\begin{equation*}
R(t)
\approx R_D + \sqrt{R_D}\,\ensuremath{y^{\prime}}\xspace\,\dfrac{t}{\tau}+\dfrac{{\ensuremath{x^{\prime}}\xspace}^2+{\ensuremath{y^{\prime}}\xspace}^2}{4}\left(\dfrac{t}{\tau}\right)^2,
\end{equation*}
where $t/\tau$ is the decay time normalised to the average ${\ensuremath{\D^0}}\xspace$ lifetime, $R_D$ is the ratio between the DCS and CF decay rates, and $\ensuremath{x^{\prime}}\xspace$ and $\ensuremath{y^{\prime}}\xspace$ are the mixing parameters `rotated' by the strong phase difference $\delta$ between the DCS and CF amplitudes:
$x'=x\cos \delta+y\sin\delta$ and $y'=y\cos \delta-x\sin\delta$.
The time evolution of the ratio $R$ is shown in Fig.~\ref{fig:charm:03}.
Further studies with larger data samples have  also  been performed by LHCb~\cite{LHCb-PAPER-2012-038,LHCb-PAPER-2013-053,LHCb-PAPER-2016-033,LHCb-PAPER-2017-046} and the results are summarized in Table~\ref{sec:charm:tab:1}.
\begin{figure}[!tb]
\begin{center}
\includegraphics[width=0.5\textwidth]{./figs/charmCPV/charm-Fig1.pdf}
\caption{
Ratio of the ${\ensuremath{\D^0}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\kaon^+}}\xspace\pi$ to ${\ensuremath{\D^0}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\kaon^-}}\xspace{\ensuremath{\pion^+}}\xspace$ decay rates, $R$, as a function of the {\ensuremath{\D^0}}\xspace decay time normalised by the average ${\ensuremath{\D^0}}\xspace$ lifetime~\cite{LHCb-PAPER-2012-038}.
The solid (dashed) line indicates the projection of the mixing  (no-mixing) fits.
}
\label{fig:charm:03}
\end{center}
\end{figure}

\begin{table}[!t]
\begin{center}
    \caption{ Summary of the charm mixing parameters measured by LHCb  using  ${\ensuremath{\D^0}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\kaon^+}}\xspace{\ensuremath{\pion^-}}\xspace$ decays.
    Results under the assumptions of CP invariance and CP violation are given.
    }
    \label{sec:charm:tab:1}
\vskip -0.2cm
{\small
\begin{tabular}{lcccc}
\hline
 & \multicolumn{2}{c}{no CPV} &  \multicolumn{2}{c}{CPV allowed}\\
Data sample & $x'^2$ ($\times 10^{-3}$)  & $y'$ ($\times 10^{-3}$) &   $x'^2$ ($\times 10^{-3}$) & $y'$ ($\times 10^{-3}$) \\
\hline
$1.0\ensuremath{\fb^{-1}}\xspace$, $D^*$\ tag~\cite{LHCb-PAPER-2012-038} & $-0.09\pm 0.13$ & $7.2\pm 2.4$  & -  &  - \\

 \multirow{2}*{$3.0\ensuremath{\fb^{-1}}\xspace$, $D^*$\ tag ~\cite{LHCb-PAPER-2013-053}}&   \multirow{2}*{$0.055\pm 0.049$} &  \multirow{2}*{ $2.8\pm 1.0$}  & ${\ensuremath{\D^0}}\xspace$: $\phantom{0} 0.049\pm 0.070$  & $5.1\pm 1.4$   \\

 & & & ${\ensuremath{\Dbar{}^0}}\xspace$: $\phantom{0} 0.060\pm 0.068$ & $4.5\pm 1.4$  \\

 \multirow{2}*{$3.0\ensuremath{\fb^{-1}}\xspace$, $B$\ tag ~\cite{LHCb-PAPER-2016-033}}&   \multirow{2}*{$0.028\pm 0.310$} &   \multirow{2}*{$4.6\pm 3.7$} & ${\ensuremath{\D^0}}\xspace$: $-0.019\pm 0.447$ & $5.81\pm 5.26$  \\

 & & &${\ensuremath{\Dbar{}^0}}\xspace$: $\phantom{0} 0.079\pm 0.433$ & $3.32\pm 5.23$  \\

 \multirow{2}*{$5.0\ensuremath{\fb^{-1}}\xspace$, $D^*$\ tag ~\cite{LHCb-PAPER-2017-046}}&   \multirow{2}*{$0.039\pm 0.027$} &  \multirow{2}*{$5.28\pm 0.52$}  & ${\ensuremath{\D^0}}\xspace$: $\phantom{0} 0.061\pm 0.037$  & $5.01\pm 0.74$  \\

 & & &${\ensuremath{\Dbar{}^0}}\xspace$: $\phantom{0} 0.016\pm 0.039$ & $5.54\pm 0.74$ \\
\hline
\end{tabular}
}
\end{center}
\end{table}

The measurements in the ${\ensuremath{\D^0}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\kaon^-}}\xspace{\ensuremath{\pion^+}}\xspace$ decay are sensitive to the normalized decay-width difference $y$ and the sum $x^2+y^2$ (under the assumption of negligible CP violation), but not to the sign of the normalized mass difference $x$. One approach to solve this problem is to study the Dalitz distributions of  three-body decays.
The `golden channel' at LHCb for such studies is the decay ${\ensuremath{\D^0}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\kaon^0_{\mathrm{S}}}}\xspace{\ensuremath{\pion^+}}\xspace{\ensuremath{\pion^-}}\xspace$, where the decay to the  ${\ensuremath{\kaon^0_{\mathrm{S}}}}\xspace{\ensuremath{\pion^+}}\xspace{\ensuremath{\pion^-}}\xspace$  final state  proceed mainly via the following three  processes with different intermediate  resonances:
1) the ${\ensuremath{\kaon^0_{\mathrm{S}}}}\xspace{\ensuremath{\rhomeson^0}}\xspace$ process with ${\ensuremath{\rhomeson^0}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\pion^+}}\xspace{\ensuremath{\pion^-}}\xspace$, which is common for both ${\ensuremath{\D^0}}\xspace$ and ${\ensuremath{\Dbar{}^0}}\xspace$ mesons;
2) the ${\ensuremath{\kaon^{*-}}}\xspace{\ensuremath{\pion^+}}\xspace$ process with ${\ensuremath{\kaon^{*-}}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\kaon^0_{\mathrm{S}}}}\xspace{\ensuremath{\pion^-}}\xspace$, which is a CF decay;
and 3) the ${\ensuremath{\kaon^{*+}}}\xspace{\ensuremath{\pion^-}}\xspace$ process with ${\ensuremath{\kaon^{*+}}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\kaon^0_{\mathrm{S}}}}\xspace{\ensuremath{\pion^+}}\xspace$, which is either a DCS decay or  {\ensuremath{\D^0}}\xspace-{\ensuremath{\Dbar{}^0}}\xspace oscillation followed by a CF decay.

In the Dalitz phase space, the DCS and CF decay amplitudes of the ${\ensuremath{\D^0}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\kaon^0_{\mathrm{S}}}}\xspace{\ensuremath{\pion^+}}\xspace{\ensuremath{\pion^-}}\xspace$ decay populate the same space and interfere. Therefore, the parameters $x$ and $y$ can be determined by measuring the strong phase difference between the contributing amplitudes in an amplitude analysis, or by importing the average strong-phase difference in  regions of phase space obtained by  ${\ensuremath{\Pe^+\Pe^-}}\xspace$ experiments operating at the energy of the $\psi(3770)$ resonance.
The latter approach is employed in several LHCb measurements~\cite{LHCb-PAPER-2015-042,LHCb-PAPER-2019-001,LHCb-PAPER-2021-009}, and the most recent measurement~\cite{LHCb-PAPER-2021-009} led to the first observation of a nonzero mass difference between the two mass eigenstates in the ${\ensuremath{\D^0}}\xspace$-{\ensuremath{\Dbar{}^0}}\xspace system.
The results of these measurements using  ${\ensuremath{\D^0}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\kaon^0_{\mathrm{S}}}}\xspace{\ensuremath{\pion^+}}\xspace{\ensuremath{\pion^-}}\xspace$ decays are summarised in Table~\ref{sec:charm:tab:2}.
\begin{table}[!tb]
\begin{center}
    \caption{
    Summary of the charm mixing parameters measured by LHCb using  ${\ensuremath{\D^0}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\kaon^0_{\mathrm{S}}}}\xspace{\ensuremath{\pion^+}}\xspace{\ensuremath{\pion^-}}\xspace$ decays.
    }
    \label{sec:charm:tab:2}
\vskip -0.2cm
\begin{tabular}{ccc}
\hline
 & \multicolumn{2}{c}{ CP-averaged parameters} \\
Data sample & $x$ ($\times 10^{-3}$)  & $y$ ($\times 10^{-3}$) \\
\hline
$1.0\ensuremath{\fb^{-1}}\xspace$, $D^*$\ tag~\cite{LHCb-PAPER-2015-042} & $-8.6\pm5.3\pm1.7$  & $0.3\pm4.6\pm1.3$  \\
$3.0\ensuremath{\fb^{-1}}\xspace$, $B$\ tag~\cite{LHCb-PAPER-2019-001} & $2.7\pm1.6\pm0.4$ & $7.4\pm3.6\pm1.1$ \\
$5.4\ensuremath{\fb^{-1}}\xspace$, $D^*$\ tag~\cite{LHCb-PAPER-2021-009} & $3.97 \pm 0.46 \pm 0.29$ & $4.59 \pm 1.20 \pm 0.85$ \\
\hline
 &  \multicolumn{2}{c}{ CP-violating parameters}\\
Data sample &  $\Delta x$ ($\times 10^{-3}$) & $\Delta y$ ($\times 10^{-3}$)\\
\hline
$3.0\ensuremath{\fb^{-1}}\xspace$, $B$\ tag~\cite{LHCb-PAPER-2019-001} & $-0.53\pm0.70\pm0.22$ & $0.6\pm1.6\pm0.3$ \\
$5.4\ensuremath{\fb^{-1}}\xspace$, $D^*$\ tag~\cite{LHCb-PAPER-2021-009} & $-0.27 \pm 0.18 \pm 0.01$ & $0.20 \pm 0.36 \pm 0.13$ \\
\hline
\end{tabular}
\end{center}
\end{table}

As mentioned in Section~\ref{sec:04:01:05gammaCombination}, different from the past LHCb $\gamma$ combinations, the recent combination exploited the LHCb measurements that are sensitive to the CKM angle $\gamma$ and to the charm mixing parameters, and the $\gamma$ angle and charm mixing parameters are simultaneously determined~\cite{LHCb-PAPER-2021-033}.
The motivation for the simultaneous combination is as follows:
\begin{itemize}
\item The $\gamma$ angle and the strong phase difference between the interfering $B$ decays are now so precisely constrained by the large $B$-meson samples that the strong phase difference, $\delta_{D}^{K\pi}$, between the decays ${\ensuremath{\D^0}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\kaon^-}}\xspace{\ensuremath{\pion^+}}\xspace$ and ${\ensuremath{\Dbar{}^0}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\kaon^-}}\xspace{\ensuremath{\pion^+}}\xspace$ can achieve a precision of about a factor of  two better than the previous world average~\cite{HFLAV18}.
This improvement can then used to improve the precision of the charm mixing parameters $x$ and $y$.
\item Due to non-negligible effects originating from charm-meson mixing, a simultaneous combination is needed to obtain an unbiased determination of the $\gamma$ angle and the charm mixing parameters $x$ and $y$.
\end{itemize}

In the charm sector, the inputs used in the combination are obtained from the time-dependent measurements of ${\ensuremath{\D^0}}\xspace\ensuremath{\rightarrow}\xspace h^+ h^-$, ${\ensuremath{\D^0}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\kaon^+}}\xspace{\ensuremath{\pion^-}}\xspace$, ${\ensuremath{\D^0}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\kaon^\pm}}\xspace{\ensuremath{\pion^\mp}}\xspace{\ensuremath{\pion^+}}\xspace{\ensuremath{\pion^-}}\xspace$, and ${\ensuremath{\D^0}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\kaon^0_{\mathrm{S}}}}\xspace{\ensuremath{\pion^+}}\xspace{\ensuremath{\pion^-}}\xspace$  decays performed by LHCb~\cite{LHCb-PAPER-2019-006,LHCb-PAPER-2014-013,LHCb-PAPER-2015-055,LHCb-PAPER-2018-038,LHCb-PAPER-2014-069,LHCb-PAPER-2016-063,LHCb-PAPER-2019-032,LHCb-PAPER-2020-045,LHCb-PAPER-2016-033,LHCb-PAPER-2017-046,LHCb-PAPER-2015-057,LHCb-PAPER-2019-001,LHCb-PAPER-2021-009,LHCb-PAPER-2015-042}.
Figure~\ref{fig:charm:05:01:combinedXY} shows the two-dimensional profile likelihood contours in the $x$-$y$ plane.
The values of $x$ and $y$, determined in the simultaneous combination, are found to be
\begin{equation*}
x=(0.400{}_{-0.053}^{+0.052})\%,\quad
y=(0.630{}_{-0.030}^{+0.033})\%.
\end{equation*}
These results provide the most precise determinations of the parameters $x$ and $y$.
Particularly, the precision of $y$ is improved by a factor of two with respect to the current world average~\cite{HFLAV18}.
\begin{figure}[!tb]
\begin{center}
\includegraphics[width=0.5\textwidth]{./figs/charmCPV/charm-Fig-combineXY.pdf}
\caption{
Profile likelihood contours (from 1 to 5 standard deviations) of the charm mixing parameters $x$ and $y$.
The solid (brown) contours are determined from the simultaneous combination~\cite{LHCb-PAPER-2021-033}, while the dashed (blue) indicate the current world average from Ref.~\cite{HFLAV18}.
}
\label{fig:charm:05:01:combinedXY}
\end{center}
\end{figure}

\subsection{CP violation}

\subsubsection{Time-integrated CP violation}

The time-integrated CP asymmetry, $\mathcal{A}_{CP}$,  in the  decay $D\ensuremath{\rightarrow}\xspace f$ is dominated by the direct CP asymmetry.
Its measurement follows the formula
\begin{equation}
    \mathcal{A}_{CP}=\mathcal{A}_{raw}-\mathcal{A}_{prod}-\mathcal{A}_{det},
\end{equation}
where $\mathcal{A}_{prod}$ denotes the
 meson production asymmetry between the $c$-hadron and its antiparticle,
$\mathcal{A}_{det}$ represents the detection asymmetry, and $\mathcal{A}_{raw}$ is the raw asymmetry between the yields of $D\ensuremath{\rightarrow}\xspace f$
and  $\overline{D}\ensuremath{\rightarrow}\xspace\overline{f}$ decays.
Often the difference of  CP asymmetries between two different decay processes  are measured, which is defined as
\begin{equation}
 \Delta \mathcal{A}_{CP}\equiv \mathcal{A}_{CP}(D\ensuremath{\rightarrow}\xspace f_1) -\mathcal{A}_{CP}(D\ensuremath{\rightarrow}\xspace f_2),
\end{equation}
 where $f_1$ and $f_2$ are two different final states with similar typologies. The effects of production asymmetry and indirect CP asymmetry as well as part of the detection asymmetries  on $\Delta \mathcal{A}_{CP}$
are largely cancelled.

Two-body decays of ${\ensuremath{\D^0}}\xspace$ mesons are particularly interesting due to their super large sample sizes,
 which are crucial for probing the tiny CP violating effects.
The first observation of CP violation in the charm sector was reported by LHCb in 2019 using the Run 2 data~\cite{LHCb-PAPER-2019-006}.
The difference of $\mathcal{A}_{CP}$ between the ${\ensuremath{\D^0}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\kaon^+}}\xspace{\ensuremath{\kaon^-}}\xspace$ and ${\ensuremath{\D^0}}\xspace\ensuremath{\rightarrow}\xspace{\ensuremath{\pion^+}}\xspace{\ensuremath{\pion^-}}\xspace$ decays, $\Delta \mathcal{A}_{CP}\equiv  \mathcal{A}_{CP}({\ensuremath{\kaon^+}}\xspace{\ensuremath{\kaon^-}}\xspace)-\mathcal{A}_{CP}({\ensuremath{\pion^+}}\xspace{\ensuremath{\pion^-}}\xspace)$, was measured
with a deviation from zero corresponding to a significance of $5.3$ standard deviations.
Table~\ref{sec:charm:tab:3} summarises the $\Delta \mathcal{A}_{CP}$ results of  a series of LHCb measurements~\cite{ LHCb-PAPER-2011-023,LHCb-PAPER-2013-003,LHCb-PAPER-2014-013,LHCb-PAPER-2015-055,LHCb-PAPER-2019-006}.
\begin{table}[!tb]
\begin{center}
    \caption{
    Summary of the $\Delta \mathcal{A}_{CP}$ results from the LHCb measurements in the charm sector.
    }
    \label{sec:charm:tab:3}
\vskip -0.2cm
\begin{tabular}{lc}
\hline
Data sample &  $\Delta \mathcal{A}_{CP} $ ($\times 10^{-3}$) \\
\hline
$0.62\ensuremath{\fb^{-1}}\xspace$, $D^*$\ tag~\cite{LHCb-PAPER-2011-023} &  $-8.2\pm 4.1\pm0.6$\\
$1.0\ensuremath{\fb^{-1}}\xspace$, $B$\ tag~\cite{LHCb-PAPER-2013-003} &  $4.9\pm3.0\pm1.4$\\
$3.0\ensuremath{\fb^{-1}}\xspace$, $B$\ tag~\cite{LHCb-PAPER-2014-013} &  $1.4\pm1.6\pm0.8$\\
$3.0\ensuremath{\fb^{-1}}\xspace$, $D^*$\ tag~\cite{LHCb-PAPER-2015-055} & $-1.0\pm0.8\pm0.3$ \\
$5.9\ensuremath{\fb^{-1}}\xspace$, $B$ or $D^*$\ tag~\cite{LHCb-PAPER-2019-006} & $-1.54\pm0.29$ \\
\hline
\end{tabular}
\end{center}
\end{table}

While  multi-body charm decays often have limited sample sizes, they can  provide excellent opportunities for CP violation measurements. The presence of intermediate resonances can lead to large variation of the strong phase difference between the interfering amplitudes, which can lead to sizeable local CP asymmetries.
Several techniques to search for CP violation in multi-body charm decays are exploited by LHCb, including amplitude analysis, the binned $\ensuremath{\chi^2}\xspace$ technique, and an unbinned technique called the energy test.
For the latter two methods, model-dependent analyses are eventually required to pin down the source in case significant CP violation were observed.

The binned $\ensuremath{\chi^2}\xspace$ technique computes  the distribution of local asymmetries and compare it with a normal distribution to judge if CP violation were observed.
An example of binned $\ensuremath{\chi^2}\xspace$ distribution in a Dalitz plot is shown in Fig.~\ref{fig:charm:Scp}.
This method relies on the optimal choice of the binning scheme.
Wide bins across resonances can lead to the cancellation of real CP asymmetries within a bin.

\begin{figure}
\begin{center}
\includegraphics[width=\textwidth]{./figs/charmCPV/charm-Xic.pdf}
\caption{(left) Distributions of local per-bin asymmetry significance $S_{CP}$ and (right) corresponding one-dimensional distributions obtained by the binned $\ensuremath{\chi^2}\xspace$ method for ${\ensuremath{\PXi}}\xspace_c^+\ensuremath{\rightarrow}\xspace p {\ensuremath{\kaon^-}}\xspace {\ensuremath{\pion^+}}\xspace$ decays. Figures taken from Ref.~\cite{LHCb-PAPER-2019-026}.
}
\label{fig:charm:Scp}
\end{center}
\end{figure}

The LHCb collaboration has developed a novel unbinned method, energy test~\cite{Williams:2011cd, Parkes:2016yie},  to perform model-independent
search for CP violation in many-body decays. With this method,  a test statistic, $T$,  is defined.  For a given data sample,
   a $p$-value for the hypothesis of CP invariance is assigned by   comparing the observed value of $T$
to the distribution of $T$ obtained from many random permutations of the data.
This method  has been applied to search for  CP violation in decays of charm mesons and beauty baryons. As an example, Fig.~\ref{fig:charm:ET} shows the global  test statistic compared with the distribution of the statistic  from many random permutations, and
the Dalitz plot distribution of significance of local test statistics
in ${\ensuremath{\D^0}}\xspace \ensuremath{\rightarrow}\xspace {\ensuremath{\pion^-}}\xspace{\ensuremath{\pion^+}}\xspace{\ensuremath{\pion^0}}\xspace$ decays.
Despite the many efforts made by LHCb and the significant improvements in the measurement precision, no  evidence of CP violation in multi-body charm decays  has ever been found  to date.
Table~\ref{sec:charm:tab:4} summarises the searches for direct CP violation in phase space of charm decays by LHCb.

\begin{figure}
\begin{center}
\includegraphics[width=0.49\textwidth]{./figs/charmCPV/charm-3pi2.pdf}
\includegraphics[width=0.49\textwidth]{./figs/charmCPV/charm-3pi1.png}
\caption{ (a) The global  test statistic compared with the distribution of the statistic calculated from a large number of random permutations obtained from the energy test applied in the CP violation searches in ${\ensuremath{\D^0}}\xspace\rightarrow {\ensuremath{\pion^-}}\xspace{\ensuremath{\pion^+}}\xspace\pi^0$ decays. (b) Dalitz plot distribution of significance of local test statistics.  Figures are taken from Ref.~\cite{LHCb-PAPER-2014-054}.
}
\label{fig:charm:ET}
\end{center}
\end{figure}

\begin{table}[!hbt]
\begin{center}
    \caption{ Summary of LHCb direct CP violation searches in phase space of charm decays.}
    \label{sec:charm:tab:4}
\vskip -0.2cm
\begin{tabular}{lcc}
\hline

Decay channel & Data sample & Method  \\
\hline
$D^+\rightarrow {\ensuremath{\kaon^-}}\xspace{\ensuremath{\kaon^+}}\xspace{\ensuremath{\pion^+}}\xspace$~\cite{LHCb-PAPER-2011-017} & $35\ensuremath{\pb^{-1}}\xspace$ & binned $\chi^2$  \\
${\ensuremath{\D^0}}\xspace\rightarrow {\ensuremath{\kaon^-}}\xspace{\ensuremath{\kaon^+}}\xspace{\ensuremath{\pion^-}}\xspace{\ensuremath{\pion^+}}\xspace$~\cite{LHCb-PAPER-2013-041} & $1.0\ensuremath{\fb^{-1}}\xspace$, $D^*$\ tag & binned $\chi^2$ \\
${\ensuremath{\D^0}}\xspace\rightarrow {\ensuremath{\pion^-}}\xspace{\ensuremath{\pion^+}}\xspace{\ensuremath{\pion^-}}\xspace{\ensuremath{\pion^+}}\xspace$~\cite{LHCb-PAPER-2013-041} & $1.0\ensuremath{\fb^{-1}}\xspace$, $D^*$\ tag & binned $\chi^2$ \\
$D^+\rightarrow {\ensuremath{\pion^-}}\xspace{\ensuremath{\pion^+}}\xspace\pip$~\cite{LHCb-PAPER-2013-057} &  $1.0\ensuremath{\fb^{-1}}\xspace$ & binned $\chi^2$ \\
${\ensuremath{\D^0}}\xspace\rightarrow {\ensuremath{\kaon^-}}\xspace{\ensuremath{\kaon^+}}\xspace{\ensuremath{\pion^-}}\xspace{\ensuremath{\pion^+}}\xspace$~\cite{LHCb-PAPER-2014-046} & $3.0\ensuremath{\fb^{-1}}\xspace$, $B$ tag & binned $\chi^2$ \\
${\ensuremath{\D^0}}\xspace\rightarrow {\ensuremath{\pion^-}}\xspace{\ensuremath{\pion^+}}\xspace\pi^0$~\cite{LHCb-PAPER-2014-054} & $2.0\ensuremath{\fb^{-1}}\xspace$, $D^*$\ tag& energy test \\
${\ensuremath{\D^0}}\xspace\rightarrow {\ensuremath{\pion^-}}\xspace{\ensuremath{\pion^+}}\xspace{\ensuremath{\pion^-}}\xspace{\ensuremath{\pion^+}}\xspace$~\cite{LHCb-PAPER-2016-044} & $3.0\ensuremath{\fb^{-1}}\xspace$, $D^*$\ tag& energy test \\
${\ensuremath{\D^0}}\xspace\rightarrow {\ensuremath{\kaon^-}}\xspace{\ensuremath{\kaon^+}}\xspace{\ensuremath{\pion^-}}\xspace{\ensuremath{\pion^+}}\xspace$~\cite{LHCb-PAPER-2018-041} & $3.0\ensuremath{\fb^{-1}}\xspace$, $B$ tag& amplitude analysis \\
$\Xi_c^+\rightarrow p {\ensuremath{\kaon^-}}\xspace {\ensuremath{\pion^+}}\xspace$~\cite{LHCb-PAPER-2019-026} & $3.0\ensuremath{\fb^{-1}}\xspace$&  binned $\chi^2$ \\
\hline
\end{tabular}
\end{center}
\end{table}

 \subsubsection{ $\ensuremath{y_{\CP}}\xspace$ and $\ensuremath{A_{\Gamma}}\xspace$ measurements}

The amplitudes of
the direct decay  of {\ensuremath{\D^0}}\xspace to a CP eigenstate and the decay after mixing can interfere and lead to indirect CP asymmetry. Its contribution to the time-integrated CP asymmetry is denoted $\mathcal{A}_{CP}^{ind}$.

Due to the  ${\ensuremath{\D^0}}\xspace-{\ensuremath{\Dbar{}^0}}\xspace$ mixing, the effective decay width of ${\ensuremath{\D^0}}\xspace$ decays to a CP-even final state (e.g. $f$=${\ensuremath{\kaon^+}}\xspace{\ensuremath{\kaon^-}}\xspace$ or ${\ensuremath{\pion^+}}\xspace{\ensuremath{\pion^-}}\xspace$),  $\Gamma_{CP+}$, differs from the average decay width, $\Gamma$.
We can define  the  parameter $\ensuremath{y_{\CP}}\xspace\equiv \Gamma_{CP+}/ \Gamma-1$ to represent the amount of CP-violation in mixing. The quantity $\ensuremath{y_{\CP}}\xspace$ is related to $x$ and $y$, $|q/p|$, and $\phi\equiv arg(q\overline{A}/pA)$,
$$
\ensuremath{y_{\CP}}\xspace\approx\dfrac{1}{2}\left(\left|\dfrac{q}{p}\right|+\left|\dfrac{p}{q}\right|\right)y\cos\phi-\dfrac{1}{2}\left(\left|\dfrac{q}{p}\right|-\left|\dfrac{p}{q}\right|\right)x\sin\phi\;.
$$
Only if CP is conserved, $\ensuremath{y_{\CP}}\xspace$ is equal to $y$.
The decay rate asymmetry is defined as
\begin{align}
    \ensuremath{A_{\Gamma}}\xspace\approx\dfrac{1}{2}\left(\left|\dfrac{q}{p}\right|-\left|\dfrac{p}{q}\right|\right)y\cos\phi-\dfrac{1}{2}\left(\left|\dfrac{q}{p}\right|+\left|\dfrac{p}{q}\right|\right)x\sin\phi \nonumber\;.
\end{align}
The asymmetry \ensuremath{A_{\Gamma}}\xspace is related to the
indirect CP asymmetry $\mathcal{A}_{CP}^{ind}$  through
$\mathcal{A}_{CP}^{ind}=-\ensuremath{A_{\Gamma}}\xspace$.

The quantities $\ensuremath{y_{\CP}}\xspace$ and $\ensuremath{A_{\Gamma}}\xspace$ can be determined by measuring the ratio of the effective lifetimes of ${\ensuremath{\D^0}}\xspace$ and $\bar{{\ensuremath{\D^0}}\xspace}$  decays to the same CP eigenstate:
\begin{align}
  \ensuremath{y_{\CP}}\xspace&=\dfrac{2\tau({\ensuremath{\D^0}}\xspace\rightarrow f_{\not{C}\not{P}})}{\tau({\ensuremath{\Dbar{}^0}}\xspace\rightarrow f_{CP})+\tau({\ensuremath{\D^0}}\xspace\rightarrow f_{CP})}-1, \\
\ensuremath{A_{\Gamma}}\xspace&=\dfrac{\tau({\ensuremath{\Dbar{}^0}}\xspace\rightarrow f_{CP})-\tau({\ensuremath{\D^0}}\xspace\rightarrow f_{CP})}{\tau({\ensuremath{\Dbar{}^0}}\xspace\rightarrow f_{CP})+\tau({\ensuremath{\D^0}}\xspace\rightarrow f_{CP})},
\end{align}
where   $f_{\not{C}\not{P}}$ denotes a non-CP-eigenstate, such as ${\ensuremath{\kaon^-}}\xspace{\ensuremath{\pion^+}}\xspace$.

In recent years, LHCb has preformed several measurements of $\ensuremath{y_{\CP}}\xspace$ and $\ensuremath{A_{\Gamma}}\xspace$, which are summarized in
Table~\ref{sec:charm:tab:ycp}~\cite{LHCb-PAPER-2011-032,LHCb-PAPER-2013-054,LHCb-PAPER-2014-069,LHCb-PAPER-2018-038,LHCb-PAPER-2016-063,LHCb-PAPER-2019-032}.
None of these measurements  shows any indication of CP violation in ${\ensuremath{\D^0}}\xspace$-${\ensuremath{\Dbar{}^0}}\xspace$ mixing or in the interference between mixing and decay. Figure~\ref{fig:charm:average} compares the $\ensuremath{y_{\CP}}\xspace$ and $\ensuremath{A_{\Gamma}}\xspace$ measurements performed by different experiments, and the  averages provided by the Heavy Flavour Averaging Group~\cite{HFLAV18} . The world averages are dominated by the measurements  by the LHCb experiment.

\begin{table}[!t]
\begin{center}
    \caption{ Summary of LHCb $\ensuremath{y_{\CP}}\xspace$ and $\ensuremath{A_{\Gamma}}\xspace$ measurements.}
    \label{sec:charm:tab:ycp}
\vskip -0.2cm
\begin{tabular}{lccc}
\hline
Data sample & Final state(s) & $\ensuremath{y_{\CP}}\xspace$ (\%)  & $\ensuremath{A_{\Gamma}}\xspace$ ($\times 10^{-3}$) \\
\hline
$29\ensuremath{\pb^{-1}}\xspace$, $D^*$\ tag~\cite{LHCb-PAPER-2011-032} & ${\ensuremath{\kaon^+}}\xspace {\ensuremath{\kaon^-}}\xspace$  & $0.55\pm0.63\pm0.41$ & $-5.9\pm 5.9\pm 2.1$ \\
$1.0\ensuremath{\fb^{-1}}\xspace$, $D^*$\ tag~\cite{LHCb-PAPER-2013-054} &  ${\ensuremath{\pion^+}}\xspace {\ensuremath{\pion^-}}\xspace$  & - & $0.33\pm1.06\pm 0.14$ \\
$1.0\ensuremath{\fb^{-1}}\xspace$, $D^*$\ tag~\cite{LHCb-PAPER-2013-054} &${\ensuremath{\kaon^+}}\xspace {\ensuremath{\kaon^-}}\xspace$  &  - & $-0.35\pm0.62\pm0.12$ \\
$3.0\ensuremath{\fb^{-1}}\xspace$, $B$\ tag~\cite{LHCb-PAPER-2014-069} &  ${\ensuremath{\pion^+}}\xspace {\ensuremath{\pion^-}}\xspace$ & - & $-0.92\pm2.6^{+0.25}_{-0.33}$\\
$3.0\ensuremath{\fb^{-1}}\xspace$, $B$\ tag~\cite{LHCb-PAPER-2014-069} & ${\ensuremath{\kaon^+}}\xspace {\ensuremath{\kaon^-}}\xspace$  & - & $-1.34\pm 0.77^{+0.26}_{-0.34}$\\
$3.0\ensuremath{\fb^{-1}}\xspace$, $B$\ tag~\cite{LHCb-PAPER-2014-069} & ${\ensuremath{\pion^+}}\xspace {\ensuremath{\pion^-}}\xspace$ \& ${\ensuremath{\kaon^+}}\xspace {\ensuremath{\kaon^-}}\xspace$  & - & $-1.25\pm0.73$  \\
$3.0\ensuremath{\fb^{-1}}\xspace$, $B$\ tag~\cite{LHCb-PAPER-2018-038} & ${\ensuremath{\pion^+}}\xspace {\ensuremath{\pion^-}}\xspace$ \& ${\ensuremath{\kaon^+}}\xspace {\ensuremath{\kaon^-}}\xspace$  & $0.57\pm 0.13\pm 0.09$ & -  \\
$3.0\ensuremath{\fb^{-1}}\xspace$, $D^*$\ tag~\cite{LHCb-PAPER-2016-063} &  ${\ensuremath{\pion^+}}\xspace {\ensuremath{\pion^-}}\xspace$ & - & $0.46\pm 0,58\pm 0.12$\\
$3.0\ensuremath{\fb^{-1}}\xspace$, $D^*$\ tag~\cite{LHCb-PAPER-2016-063} & ${\ensuremath{\kaon^+}}\xspace {\ensuremath{\kaon^-}}\xspace$  & - & $-0.30\pm 0.32\pm 0.10$\\
$3.0\ensuremath{\fb^{-1}}\xspace$, $D^*$\ tag~\cite{LHCb-PAPER-2016-063} & ${\ensuremath{\pion^+}}\xspace {\ensuremath{\pion^-}}\xspace$ \& ${\ensuremath{\kaon^+}}\xspace {\ensuremath{\kaon^-}}\xspace$  & - & $-0.13\pm 2.0\pm 0.7$  \\
$5.4\ensuremath{\fb^{-1}}\xspace$, $B$\ tag~\cite{LHCb-PAPER-2019-032} &   ${\ensuremath{\pion^+}}\xspace {\ensuremath{\pion^-}}\xspace$ & - & $0.22\pm 0.70\pm 0.08$   \\
$5.4\ensuremath{\fb^{-1}}\xspace$, $B$\ tag~\cite{LHCb-PAPER-2019-032} &  ${\ensuremath{\kaon^+}}\xspace {\ensuremath{\kaon^-}}\xspace$ & - & $-0.43\pm 0.36\pm 0.05$ \\
\hline
\end{tabular}
\end{center}
\end{table}

\begin{figure}
\begin{center}
\includegraphics[width=0.49\textwidth]{./figs/charmCPV/ycp_16jul20.png}\includegraphics[width=0.49\textwidth]{./figs/charmCPV/a_gamma_2jun21.png}
\caption{The world averages of  $\ensuremath{y_{\CP}}\xspace$ and $\ensuremath{A_{\Gamma}}\xspace$. Figures taken from Ref.~\cite{HFLAV18}
}
\label{fig:charm:average}
\end{center}
\end{figure}

\clearpage

\section{Prospects and summary}
\label{sec:upgrade}
\newcommand{\,\mathrm{cm}^{-2}\mathrm{s}^{-1}}{\,\mathrm{cm}^{-2}\mathrm{s}^{-1}}

\subsection{Upgrade plan of LHCb}
\label{sec:upgrade-plan}
The physics output discussed in this review shows that LHCb has successfully
deepened our understanding of flavour physics with experimental data taken
up to the year of 2018. Most of the the key flavour observables are measured to
an unprecedented precision, yet it is generally true that the uncertainties
are still dominated by statistical fluctuation. To further increase the availability
of high-quality collision data, the LHCb detector is currently under a major
upgrade~\cite{LHCb-TDR-012}, known as Upgrade Ia or simply Upgrade~I.
The installation will complete by the end of the second Long Shutdown (LS2)
of the LHC, and the upgraded detector will start taking data in  2022
with an instantaneous luminosity of $2 \times 10^{33} \,\mathrm{cm}^{-2}\mathrm{s}^{-1}$,
five times the value achieved so far.
Fig.~\ref{fig:upgrade-timeline} shows the plan for LHCb operation after upgrade.
LHCb aims to accumulate an integrated luminosity of approximately $23 \ensuremath{\fb^{-1}}\xspace$
by the end of Run 3 around 2025, and a total of $50 \ensuremath{\fb^{-1}}\xspace$ by the end of Run 4.
Note that during the Long Shutdown 3 (LS3) between Run 3 and 4, intensive work
will be done on the machine configuration to prepare the High Luminosity Large Hadron Collider (HL-LHC)~\cite{Aberle:2749422}.
During this period consolidation work (Upgrade~Ib) will be carried out at LHCb
with only minor change on the detector configuration or performance.
In order to fully exploit the HL-LHC potential in flavour physics, the collaboration
plans another major upgrade, Upgrade II~\cite{LHCb-PII-EoI}, to enable
the detector to operate at luminosity as high as $1.5 \times 10^{34} \,\mathrm{cm}^{-2}\mathrm{s}^{-1}$.
This will allow for an integrated luminosity of $\sim 300\ensuremath{\fb^{-1}}\xspace$ to be achieved
in the lifetime of the (HL-)LHC.

\begin{figure}[!b]
\begin{center}
\includegraphics[width=0.7\textwidth]{./figs/upgrade/upgrade_lumi_by_year.png}
\caption{Instantaneous and integrated luminosities of LHCb as functions of year.}
\label{fig:upgrade-timeline}
\end{center}
\end{figure}

To put the LHCb prospects in context, other players in the field of
flavour physics study has to be mentioned. The Belle II
experiment~\cite{Hashimoto:2004sm}, $B$-factory at the superKEKB~\cite{Ohnishi:2013fma},
has started to take collision data of ${\ensuremath{\Pe^+\Pe^-}}\xspace\ensuremath{\rightarrow}\xspace{\Y4S}\ensuremath{\rightarrow}\xspace {\ensuremath{\PB}}\xspace{\ensuremath{\offsetoverline{\PB}}}\xspace$ since 2018,
aiming to collect an integrated luminosity of $50 \ensuremath{\ab^{-1}}\xspace$ by 2025.
Belle II and LHCb are expected to be competitive in pushing the
measurement precision, though with very different systematic uncertainties,
of a number of key flavour observables such as the CKM unitarity angles,
the Wilson coefficients and the tests on lepton flavour universality.
Given the difference in production mechanism and corresponding detector features,
their complementarity should be more appreciated.
With the beam energy constraints Belle II will be generally good at treating
final states with neutral particles ($\gamma$, ${\ensuremath{\pion^0}}\xspace$) or missing energy (neutrinos);
It has unique advantage in studying the $\tau$ leptons
through ${\ensuremath{\Pe^+\Pe^-}}\xspace \ensuremath{\rightarrow}\xspace {\ensuremath{\Ptau^+\Ptau^-}}\xspace$ process.
With the huge cross section of heavy hadrons, LHCb will have larger yield in
most of the final states especially those with all-charged tracks.
The beauty hadrons produced at high-energy $pp$ collisions are highly boosted,
hence their decay vertices could be well separated from the primary vertices.
With an excellent vertex detector LHCb fully exploits this phenomenon
to suppress background for most $B$ signatures to a extremely low level.
In addition to $B$ mesons, the studies at LHCb extend to all type of heavy hadrons
like ${\ensuremath{\B^0_\squark}}\xspace$, ${\ensuremath{\B_\cquark^+}}\xspace$ and beauty baryons.
At the HL-LHC era ATLAS and CMS will keep pushing the high-energy frontier
by studying the Higgs and searching for signature of new physics beyond
the Standard Model, meanwhile the yield of beauty hadrons will be so considerable
that measurement of a few key flavour observables can be made precisely,
especially those with a pair of muons in their final state.
They are designed to perform well at high $\ensuremath{p_{\mathrm{T}}}\xspace$ and central rapidity,
perfectly complementing the kinematic range of LHCb.

\subsubsection{Upgrade I}
\label{sec:upgrade-UI}
The ongoing Upgrade I aims to increase the instantaneous luminosity of LHCb
from the current value of $4 \times 10^{32}\,\mathrm{cm}^{-2}\mathrm{s}^{-1}$ to $2 \times 10^{33}\,\mathrm{cm}^{-2}\mathrm{s}^{-1}$.
Na\"{i}vely it would mean a five-fold increase in all signal yield, which will be
roughly true for final states with muons. Actually the gain for all-hadronic
final state will be more than that due to a major change in the trigger
system~\cite{LHCb-TDR-016}.
The hardware trigger L0, which reduces the data rate from \mbox{40 MHz} to
\mbox{1 MHz}, will be completely removed after Upgrade I,
allowing a more flexible full software trigger.
Generally the trigger efficiencies for all-hadronic final states
are expected to be doubled, as taken in simulation study of upgrade
performances, however this number could vary depending on the individual channel.

The increased pile-up causes much higher combinatorial background and
more challenging track reconstruction, therefore the tracking
systems~\cite{LHCb-TDR-013,LHCb-TDR-015} have been completely redesigned
with higher granularity and better radiation tolerance so as to provide
uncompromised tracking performance at higher pile-up.
Components of the particle identification systems~\cite{LHCb-TDR-014}
will be reused as much as possible, yet the readout electronics
will be replaced in accordance with the 40 MHz readout rate.
As a result of higher luminosity, improved trigger rate and larger number of
output channels, the data volume to be treated in real time and offline
will be substantially higher, hence new software infrastructure
and computing models have been developed
correspondingly~\cite{LHCb-TDR-017,LHCb-TDR-018,LHCb-TDR-021} to ensure
physics data to be processed and stored in a timely manner.
A new subsystem will also be installed to enhance the detector's
capability in fixed target and heavy-ion studies without
disturbing the main physics program~\cite{LHCb-TDR-020}.
The simulation study shows that detector performance after upgrade
will be at least as good as before, with improvement at some areas.

\subsubsection{Upgrade II}
\label{sec:upgrade-UII}
By the end of Run 4 LHCb will have accumulated 50 $\ensuremath{\fb^{-1}}\xspace$ $pp$ collision data,
with many subdetectors reaching end of lifetime operation
at the same condition beyond the point would be less attractive.
To fully exploit the HL-LHC potential in flavour physics study,
the collaboration proposed Upgrade II towards an integrated luminosity of
$300 \ensuremath{\fb^{-1}}\xspace$~\cite{LHCb-PII-EoI}.
The physics cases with the luminosity an order of magnitude higher than
before HL-LHC time have been studied extensively by the collaboration
and summarised in a document in 2018~\cite{LHCb-PII-Physics}.
A few benchmarks will be discussed below.
Note that the HL-LHC baseline design assumes LHCb running condition
to be the same as in Run 3, the HL-LHC experts recently released a report on
the upgrade feasibility from the machine side~\cite{Efthymiopoulos:2319258}
showing that possible solution of operating at a luminosity of
$1.5 \times 10^{34}\,\mathrm{cm}^{-2}\mathrm{s}^{-1}$ will allow the target of $300 \ensuremath{\fb^{-1}}\xspace$ to be met.

The seven-fold increase of luminosity will again impose more technical
challenges for the experiment.
The expected number of interactions per crossing is around 40, twenty times
of current situation (or a hundred times of the LHCb initial design).
Fast timing resolution will be required in most subsystems to fight against
the combinatorial backgrounds caused by the pile-up.
Finer granularity in all tracking detectors is compulsory
under much higher multiplicity.
Radiation hardness will be more of concern especially for areas close
to the beampipe.
A daunting amount of 200 Tb of data will be produced every second, and has to
be reduced by four orders of magnitude before stored permanently.
New subsystems are being proposed in order to extend geometrical acceptance
for low-momentum tracks, and to improve $K/\pi$ separation at lower momentum.
A lot of development activities have been launched driven by these requirements,
while exploiting new technologies in detector and computing.
A framework Technical Design Report summarising these activities
is being prepared by the collaboration and will be released soon.

\subsection{Physics prospects}
\label{sec:upgrade-prospect}
Before the HL-LHC or by the end of Run 3, LHCb will have taken $23 \ensuremath{\fb^{-1}}\xspace$ data,
drastically reducing the statistical uncertainties for most of channels
compared with current measurements. The expected projections are studied
in detail~\cite{LHCb-TDR-012} and updated with inputs from experiences in
Run I~\cite{LHCb-PAPER-2012-031,Gerson:1604468}.
The physics opportunities in Upgrade II with $300 \ensuremath{\fb^{-1}}\xspace$
have also been studied~\cite{LHCb-PII-Physics}, which concludes
that the energy scale probed by flavour observables will be doubled
with respect to pre-HL-LHC era.
The sensitivity of a selection of key flavour observables after LHCb upgrades are listed in
Table~\ref{tab:upgrade-key-observables} and illustrated in
Fig.~\ref{fig:upgrade-key-sensitivity}, mostly from Ref.~\cite{LHCb-PII-Physics}
with minor updates when available. Note that Belle II will have completed
data taking when LHCb collects $23 \ensuremath{\fb^{-1}}\xspace$ data. Expected projection from Belle II,
ATLAS and CMS are listed for comparison when applicable.
A few highlights will be briefly mentioned here.

\begin{table}
  \begin{center}
    \caption{Sensitivity of selected key flavour observables for LHCb, ATLAS and CMS, taken from Ref.~\cite{LHCb-PII-Physics} with updates when available, and  Belle II sensitivities   from Ref.~\cite{Belle-II:2018jsg}. }
    \label{tab:upgrade-key-observables}
    \resizebox{\textwidth}{!}{%
    \begin{tabular}{lccccc}
    \hline
               & LHCb    &  LHCb       & LHCb         & Belle II & ATLAS \\
    Observable & current & ($23\ensuremath{\fb^{-1}}\xspace$) & ($300\ensuremath{\fb^{-1}}\xspace$) & ($50\ensuremath{\ab^{-1}}\xspace$) & \& CMS \\
    \hline
    \textbf{CKM tests} & & & & & \\
    $\gamma$ (all modes) & $4^{\circ}$~\cite{LHCb-CONF-2020-003,LHCb-PAPER-2021-033}
    & $1.5^{\circ}$ & $0.35^{\circ}$ & $1.5^{\circ}$ & $-$ \\
    $\gamma$ ($B_s^0 \ensuremath{\rightarrow}\xspace D_s^+ K^-$) & $(_{-22}^{+17})^{\circ}$ & $4^{\circ}$ & $1^{\circ}$
    & $-$ & $-$ \\
    $\sin 2\beta$ & 0.04~\cite{LHCb-PAPER-2017-029} & 0.011
    & 0.003 & 0.005 & $-$ \\
    \multirow{2}{*}{$\phi_s$ ($B_s^0 \ensuremath{\rightarrow}\xspace {\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace \phi$)} &
    \multirow{2}{*}{49 mrad~\cite{LHCb-PAPER-2014-059}} & \multirow{2}{*}{14 mrad} &
    \multirow{2}{*}{4 mrad} & \multirow{2}{*}{-} & 22 mrad~\cite{TheATLAScollaboration:2013ppf}\\
    & & & & & 5-6 mrad~\cite{CMS:2018wpq} \\
    $\phi_s$ ($B_s^0 \ensuremath{\rightarrow}\xspace D_s^+D_s^-$) & 170 mrad~\cite{LHCb-PAPER-2014-051} & 35 mrad
    & 9 mrad & $-$ & $-$ \\
    $\phi_s^{{\ensuremath{\squark\squarkbar}}\xspace{}s}$ ($B_s^0 \ensuremath{\rightarrow}\xspace \phi\phi$) & 154 mrad~\cite{LHCb-PAPER-2014-026} & 39 mrad
    & 11 mrad & $-$ & feasible~\cite{CMS:2017gvo} \\
    $a_{sl}^{s}$ & $33 \times 10^{-4}$~\cite{LHCb-PAPER-2016-013} & $10 \times 10^{-4}$
    & $3 \times 10^{-4}$ & $-$ & $-$ \\
    $|V_{ub}|/|V_{cb}|$ & 6\%~\cite{LHCb-PAPER-2015-013} & 3\% & 1\% & 1\% & $-$ \\
    \hline
    \textbf{Charm} & & & & & \\
    $\Delta {\ensuremath{{\mathcal{A}}^{{\ensuremath{C\!P}}\xspace}}}\xspace$ & $2.9 \times 10^{-4}$~\cite{LHCb-PAPER-2019-006} & $1.7 \times 10^{-4}$
    & $3.0 \times 10^{-5}$ & $5.4 \times 10^{-4}$ & $-$ \\
    $\ensuremath{A_{\Gamma}}\xspace$ & $1.3 \times 10^{-4}$~\cite{LHCb-PAPER-2020-045} & $4.2 \times 10^{-5}$
    & $1.0 \times 10^{-5}$ & $3.5 \times 10^{-4}$ & $-$ \\
    \hline
    \textbf{${\ensuremath{\B_{(\squark)}^0}}\xspace \ensuremath{\rightarrow}\xspace \mu^+ \mu^-$} & & & & & \\
    $\frac{{\ensuremath{\mathcal{B}}}\xspace(B^0 \ensuremath{\rightarrow}\xspace \mu^+ \mu^-)}{{\ensuremath{\mathcal{B}}}\xspace(B_s^0 \ensuremath{\rightarrow}\xspace \mu^+ \mu^-)}$
    & 71\%~\cite{LHCb-PAPER-2021-007,LHCb-PAPER-2021-008} & 34\% & 10\%
    & $-$ & 21\%~\cite{CMS:2015iha,CMS:2018uxm} \\
    $\tau_{B_s^0 \ensuremath{\rightarrow}\xspace \mu^+ \mu^-}$ & 14\%~\cite{LHCb-PAPER-2021-007,LHCb-PAPER-2021-008}
    & 8\% & 2\% & $-$ & $-$ \\
    \hline
    \textbf{EW penguins} & & & & & \\
    $R_K$ ($B^+ \ensuremath{\rightarrow}\xspace K^+ {\ensuremath{\ell^+}}\xspace{\ensuremath{\ell^-}}\xspace$) & 0.044~\cite{LHCb-PAPER-2021-004} & 0.025 & 0.007
    & 0.036 & $-$ \\
    $R_{K^{*}}$ ($B^0 \ensuremath{\rightarrow}\xspace {\ensuremath{\kaon^{*0}}}\xspace{\ensuremath{\ell^+}}\xspace{\ensuremath{\ell^-}}\xspace$) & 0.10~\cite{LHCb-PAPER-2017-013} & 0.031 & 0.008
    & 0.032 & $-$ \\
    \textbf{LFU tests} & & & & & \\
    $R_{D^*}$ (${\ensuremath{\B^0}}\xspace \ensuremath{\rightarrow}\xspace {\ensuremath{\D^{*-}}}\xspace {\ensuremath{\ell^+}}\xspace \nu$) & 0.026~\cite{LHCb-PAPER-2015-025,LHCb-PAPER-2017-027}
    & 0.007 & 0.002 & 0.005 & $-$ \\
    $R_{{\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace}$ ($B_c^+ \ensuremath{\rightarrow}\xspace {\ensuremath{{\PJ\mskip -3mu/\mskip -2mu\Ppsi}}}\xspace{\ensuremath{\ell^+}}\xspace \nu$) & 0.24~\cite{LHCb-PAPER-2017-035}
    & 0.07 & 0.02 & $-$ & $-$ \\
    \hline
  \end{tabular}
  }
\end{center}
\end{table}

\begin{figure}
\begin{center}
 \includegraphics[width=0.8\textwidth]{./figs/upgrade/key_sensitivity_CPV_RD_LFU.png}
\caption{Sensitivity to probe key CP violating variables, rare decay and lepton flavour universality tests expected from LHCb upgrades. Anticipated results from Belle II, ATLAS or CMS are listed when available. The figure is taken from Ref.~\cite{LHCb-PII-Physics}.}
\label{fig:upgrade-key-sensitivity}
\end{center}
\end{figure}

\textbf{Rare decays:} The decay of $B^0 \ensuremath{\rightarrow}\xspace \mu^+\mu^-$ is not very far from being observed
with imminent Upgrade I data, and its branching fraction with respect to $B_s^0 \ensuremath{\rightarrow}\xspace \mu^+ \mu^-$
will be measured with 10\% uncertainty giving a powerful test of minimal flavour violation.
A wide range of studies will be performed in $b \ensuremath{\rightarrow}\xspace s {\ensuremath{\ell^+}}\xspace{\ensuremath{\ell^-}}\xspace$ or $b \ensuremath{\rightarrow}\xspace u {\ensuremath{\ell^+}}\xspace{\ensuremath{\ell^-}}\xspace$ decays with
improved precision, so the current hint of discrepancy in $R_{K^{(*)}}$ with SM predictions
will be confirmed or excluded with confidence. A series of tests on lepton flavour
universality can be carried out in $b \ensuremath{\rightarrow}\xspace c {\ensuremath{\ell^-}}\xspace \nu$ decays. The precision of $R_{D^*}$
will reach per mille level.

\textbf{CKM tests:}  The CKM unitarity triangle will be so precisely determined in the future
that discrepancies between various measurement caused by physics beyond Standard Model will
be extremely difficult to hide, as shown in Fig.~\ref{fig:upgrade-CKM-angles}.
The angle $\gamma$, currently still the least well-known,
will be determined with an uncertainty of $1.5^{\circ}$ after Run 3, similar
as the precision expected from Belle II; The uncertainty will be further reduced to
$0.35^{\circ}$ after Upgrade II.
The expected precision on $B_s$ weak mixing angle $\phi_s$ will be pushed to a few mrad,
the same level as the current precision determined indirectly from CKM fit using
tree-level measurements.

\begin{figure}
\begin{center}
  \includegraphics[width=0.48\textwidth]{./figs/upgrade/CKM_UI.png}
  \includegraphics[width=0.48\textwidth]{./figs/upgrade/CKM_UII.png}
\caption{LHCb constraints to the unitarity triangle with anticipated improvement from (left) Upgrade I and (right) Upgrade II. The figure is taken from Ref.~\cite{LHCb-PII-Physics}.}
\label{fig:upgrade-CKM-angles}
\end{center}
\end{figure}

\textbf{CP violation in charm:} After Upgrade II LHCb will be able to probe CP violation
in charm with a sensitivity of $\mathcal{O}(10^{-5})$, the only future facility promising
in observing indirect CP violation in charm which is predicted to be $\mathcal{O}(10^{-4})$
or less in the SM.

\textbf{Heavy-flavour spectroscopy:}
LHCb has demonstrated its capability as a general purpose
detector in the forward region, and it will be better equipped in this respect after Upgrades.
With data sample over an order of magnitude large than currently available, LHCb will
systematically update our knowledge on heavy hadron spectroscopy, especially those labelled
as exotic nowadays such as pentaquark and tetraquark states.

\textbf{Beyond flavour physics:}
In the high $\ensuremath{p_{\mathrm{T}}}\xspace$ range LHCb could also make contribution complementary to ATLAS
and CMS, for instance to the precision determination of the effective weak mixing angle
$\sin\theta_W^2$ and the $W$ mass. Measurements of top pair and gauge boson production
at LHCb are also crucial to study the poorly known gluon parton distribution functions
at high-$x$ range. This is an important study in the QCD, which help to understand the
ubiquitous background for any new high-mass states in ATLAS or CMS.
LHCb will also push sensitivity in search for dark-photon and long-lived particles
predicted in several NP scenarios.

\subsection{Summary}
\label{sec:upgrade-summary}
This manuscript briefly reviews the recent experimental highlights using data collected
with the LHCb detector. They are not only from the flavour physics benchmarks that the experiment
was designed for, but also include unexpected discoveries revealing LHCb's
capability as a general purpose detector in the forward region:
\begin{itemize}
    \item A large variety of new particles are discovered, either filling gaps in conventional heavy hadron spectroscopy or establishing new types of their own, like pentaquarks or tetraquarks;
    \item Some processes predicted to be extremely rare in the SM are observed, such as the $B_s^0 \ensuremath{\rightarrow}\xspace \mu^+ \mu^-$ decay. Precise measurements are performed on semileptonic and radiative $b$-decays, where NP at high energy scale can be probed with promising sensitivity, and tensions with the SM are found in a series of studies;
    \item The CKM parameters are determined precisely using multiple approaches in a wide range of final states. The angle $\gamma$, which was the least known in the unitarity triangle, has been determined with an unprecedented precision of $4^{\circ}$;
    \item Heavy flavour study in the charm sector witnesses a couple of milestones, such as the observation of ${\ensuremath{\D^0}}\xspace-{\ensuremath{\Dbar{}^0}}\xspace$ mixing in a single measurement and observation of non-zero mass difference between the two mass eigenstates in the ${\ensuremath{\D^0}}\xspace-{\ensuremath{\Dbar{}^0}}\xspace$ system. Precision in probing the CP violation in charm keeps pushing forward.
\end{itemize}
This is by far not a complete list~\cite{LHCbpapers}, and many  interesting topics studied at LHCb
are not covered  due to limited space. With Upgrade I approaching completion, the LHCb detector
will resume operation at higher luminosity next year while preparing for the future Upgrade II.
New exciting physics results are expected, which will continue to shape the landscape of heavy flavour physics and beyond.

\clearpage
\section*{Acknowledgements}

This work is partially supported by
the National Key Research and Development Program of China under grant No. 2017YFA0402100,
National Natural Science Foundation of China (NSFC) under grant Nos. 11435003, 11575091, 11575094, 11925504, 11975015, 12175245, 12175005, 11705209, 11961141015, 12061141007,
Chinese Academy of Sciences, the Fundamental Research Funds for the Central Universities,
 Peking University Funds for the New Faculty Startup program.
We thank Franz Muheim and Niels Tuning for suggestions in improving the draft.

High Energy Physics - Experiment (hep-ex); High Energy Physics - Phenomenology (hep-ph)