\section{Introduction}

The study of electromagnetic wave interactions with time-varying media dates back to more than half a century ago, when several pioneering works \cite{morgenthaler1958velocity, simon1960action, cassedy1963dispersion, cassedy1965waves, holberg1966parametric, cassedy1967dispersion, auld1968signal, felsen1970wave, fante1971transmission, chu1972wave} revealed the rich physics arising due to such interactions. Although in the subsequent decades several other works have further explored the implications of wave propagation in dynamic media \cite{yablonovitch1973spectral, stepanov1976dielectric, kalluri1992frequency, stepanov1993waves}, it was not until recently when the field gained considerable momentum \cite{caloz2019spacetime} partly owing to the improved fabrication and characterization techniques at various frequency regimes \cite{shaltout2019spatiotemporal, taravati2020space, reshef2019nonlinear}. Consequently, it was soon discovered that time-dependent materials can enable many interesting and novel phenomena, such as interband photonic transitions \cite{winn1999interband,yu2009complete}, nonreciprocal systems \cite{sounas2017non, williamson2020integrated}, efficient wavelength conversion \cite{notomi2006wavelength, preble2007changing, zhou2020broadband, khurgin2020adiabatic}, temporal aiming \cite{pacheco2020temporal}, Fresnel drag \cite{huidobro2019fresnel}, negative extinction \cite{shcherbakov2019time}, entangled photon generation \cite{cirone1997photon, sloan2021casimir}, parametric oscillations \cite{lee2021parametric}, synthetic dimensions \cite{yuan2021tutorial}, topological phase transition \cite{ozawa2019topological}, nonreciprocal gain \cite{koutserimpas2018nonreciprocal}, among others \cite{caloz2019spacetime2}. Several recent publications also further investigated the fundamental aspects of wave propagation in time-varying media, such as Refs. \cite{shvartsburg2005optics, biancalana2007dynamics, xiao2011spectral, bakunov2014reflection, plansinis2015temporal, plansinis2016spectral}, thus providing a better understanding of the physical implications of time-varying media.

Apart from enabling novel effects, time-varying materials and components may be also used to improve a given device to provide better functionality in terms of operational bandwidth, size, efficiency of the relevant process, etc. \cite{caloz2019spacetime2}. In some cases, such improvements may even exceed what is normally allowed by well-established physical bounds, as conventional physical bounds are typically derived under the assumption of linearity and time-invariance (LTI systems). Indeed, as mentioned above, time-variance in the system properties may enable various functionalities not present in a conventional setting, such as parametric gain \cite{koutserimpas2018nonreciprocal} or spectral changes \cite{winn1999interband, notomi2006wavelength, preble2007changing, zhou2020broadband, khurgin2020adiabatic,hayran2021spectral}, which are generally not accounted for when deriving conventional bounds for LTI systems, and can therefore be exploited to bypass such physical limits. Consequently, time-varying materials and components provide an exceptional opportunity to realize superior devices not restricted by various conventional limits.

\section{Applications}

In this section we provide a concise overview of the most relevant applications of time-varying electromagnetic systems in the context of overcoming various physical bounds of conventional LTI systems. In this regard, this section is divided into various subsections, each of which focuses on a specific application or research topic, with the goal to systematically guide the reader through this emerging research area. We will not discuss some major research directions that also seek to overcome various issues present in conventional systems and broaden their functionality, such as nonreciprocity \cite{sounas2017non, williamson2020integrated} or synthetic dimensions \cite{yuan2021tutorial} based on time modulation, since these topics are already covered in detail in the literature. Instead, we will focus on topics of particular interest to the IEEE Antennas and Propagation community, including how time-varying systems impact the fundamental limitations of antennas \cite{hansen1981fundamental}, broadband impedance matching networks \cite{fano1950theoretical}, and radar absorbers \cite{rozanov2000ultimate}.

\subsection{Impedance Matching}
\label{section_impedance}

Efficient power transfer from a source to a load with minimal reflections is arguable one of the most important tasks in electromagnetic engineering. Thus, when designing impedance matching networks to minimize reflections, it becomes essential to identify the limitations and the optimal configurations of the system. In this regard, a matching-bandwidth limit was introduced for linear, time-invariant and passive systems by Bode in 1945 \cite{bode1945network}, which was later generalized by Fano for arbitrary load impedances \cite{fano1950theoretical}. The resulting Bode-Fano limit provides a theoretical upper bound on the bandwidth over which a certain reflection reduction can be attained, independent of the complexity of the employed matching network (only assumed to be reactive), and it is therefore a very powerful tool to assess the maximum bandwidth of various electromagnetic systems, such as antennas \cite{gustafsson2006bandwidth} or optical invisibility cloaks \cite{monticone2016invisibility}. To overcome this limitation one strategy that has been explored in the literature is to break the underlying assumption of passivity. Indeed, several proposals utilizing active ``non-Foster'' elements have shown that the bandwidth of small antennas \cite{sussman2009non} and cloaking devices \cite{hrabar2010towards, chen2013broadening, soric2014wideband, chen2019active} can be extended beyond what is possible with a passive system. However, active elements typically increase the noise level of the system \cite{jacob2016gain} and can lead to potential instabilities (unbounded oscillations in time) \cite{chen2019active, ugarte2012stability, abdelrahman2021physical}, thereby disrupting the desired functionality.

An alternative option to go beyond the Bode-Fano limit is to break the assumption of time-invariance, i.e., allow temporal variation in the system properties. In this context, Ref. \cite{wang2007time2} employed time-varying ``switched'' matching networks to achieve a better bandwidth performance compared to time-invariant networks. More recently, Ref. \cite{shlivinski2018beyond} provided another insightful method, which involves switching the system parameters while a broadband signal (a pulse) is within the matching network and before it reaches the load, as shown in Figs. \ref{impedance}(a) and \ref{impedance}(b). Typically, the temporal switching would lead to reflected waves due to the temporal discontinuity \cite{mendoncca2002time}, however it was shown that the transmission efficiency can still be maximized by optimizing several parameters, including the impedance of the matching transmission line and the phase velocity of the pulse within the matching line before and after the switching event \cite{shlivinski2018beyond}. Since temporal switching increases the number of free parameters that can be tuned compared to the time-invariant case, one may achieve optimal performances that are otherwise out of reach in a conventional time-invariant setting, especially for large load-source impedance contrast values (see Fig. \ref{impedance}(c)) \cite{shlivinski2018beyond}. This proposal was further extended from hard (abrupt) to ``soft'' temporal switching in Ref. \cite{hadad2020soft}, which could be more feasible to implement in practice. An adverse byproduct of the temporal switching operation is the spectral compression or broadening of the pulse, which necessitates the use of additional analog or digital devices to restore the original pulse spectrum and, thus, complicates the overall matching system. Moreover, this method is feasible only for short-time signals (as quasi-monochromatic signals would require impractically long matching lines), and requires knowledge that the pulse is present within the matching line (which implies the presence of some fast detection mechanism able to trigger the switching event). Even with these limitations, this approach show a realistic path toward ultra-broadband matching networks beyond the Bode-Fano constraints.

\begin{figure}[t]
\centering
\includegraphics[width=1.0\linewidth]{impedance.pdf}
\caption{(a) A temporally switched transmission line can be employed to match a source impedance $R\textsubscript{G}$ to a load impedance $R\textsubscript{L}$ beyond the Bode-Fano limit. (b) Such a switched matching line can be realized with parallel varactor diodes and banks of series inductors, whose properties are changed abruptly in time at $t = t\textsubscript{s}$. (c) The maximum transmission efficiency for various $R\textsubscript{L}/R\textsubscript{G}$ values shows that the temporally switched case (solid blue) provides a better performance compared to the non-switched time-invariant case (dashed red line) especially for large source-load contrast values. (d),(e) A transmission line terminated with a time-varying inductance $L(t)$. (d) With a suitable time modulation, the inductor can accumulate energy without any reflections. (e) By time-reversing the temporal modulation in (d) the energy trapped by the inductor can be released. (f) Time modulation profile of the inductance comprising both the trapping (before $t=2.5$ s) and releasing regimes (after $t=2.5$ s). (g) Temporal variation of the incident signal and the reflected signal (obtained theoretically and through simulations) at the load position for a load inductance temporally modulated as in panel(f). Panels (a-c) are reproduced with permission from Ref. \cite{shlivinski2018beyond}, American Physical Society. Panels (d-g) are reproduced with permission from Ref. \cite{mirmoosa2019time}, American Physical Society.}
\label{impedance}
\end{figure}

For the quasi-monochromatic regime, on the other hand, a transmission line where the load reactance itself is varied in time has been used to impedance match a time-harmonic source a purely reactive load and to accumulate and release energy on demand (see Fig. \ref{impedance}(d-g)) \cite{mirmoosa2019time}. In this context, it was noted that a time-modulated reactance can emulate a resistance and, can, therefore act as a ``virtual'' absorber \cite{mirmoosa2019time}. An intuitive picture explaining this result was also provided, showing that temporally modulating the reactance with such a tangent function is equivalent to a short-circuited transmission line with the shorted end moving away from the source with the same velocity as the phase velocity of the incident signal \cite{mirmoosa2019time}. Therefore, the incident wave does not reflect back since it never ``reaches'' the reflecting termination. Consequently, energy from a time-harmonic source can be trapped with no reflections and released on demand (see Figs. \ref{impedance}(f) and \ref{impedance}(g)).

Lastly, another interesting application of temporal modulation in the context of impedance matching is related to the problem of capturing and trapping an incident electromagnetic wave into a \emph{lossless} cavity, with no reflections. This can be achieved by temporally modulating the radiation leakage from the cavity, such that it perfectly cancels the direct reflection from the cavity at all time instants \cite{sounas2020virtual}. In this regard, we have also shown that by temporally modulating the shell of a compact core-shell scatterer (acting as an open cavity supporting a nonradiating mode), an incident broadband pulse can be captured within the cavity with no reflection, effectively forcing a broadband signal into a narrowband resonator with ideal efficiency (this may also be interpreted as spectrum compression through temporal modulation with the additional constraint of zero reflections). \cite{hayran2021capturing} See Supplementary Material for further details.

\subsection{Antenna Performance Enhancement}

Antenna miniaturization is another major research area within the applied electromagnetism community with many potential benefits for various areas, such as communication systems, sensor networks, implanted devices, and navigation systems. While it is generally possible to reduce the size of an antenna down to subwavelength scales for a given design frequency, this typically comes at the price of reduced bandwidth and radiation efficiency. In this regard, the trade-off between size and performance for electrically small antennas was investigated more than half a century ago in several seminal papers by Wheeler \cite{wheeler1947fundamental}, Chu \cite{chu1948physical}, and later by Harrington \cite{harrington1960effect}, resulting in, as it is known today, the Chu-Harrington limit \cite{yarman2008design}. The limit as derived by Chu is based on an equivalent circuit model for the modal impedance of the lowest-order spherical wave radiated by an antenna enclosed within a spherical volume. Regardless of the antenna structure and LTI material, this equivalent circuit for spherical waves becomes more reactive as the size of the spherical volume is reduced, leading to a higher Q factor and a smaller bandwidth (consistent with the fact that, for example, a short wire dipole has a strongly reactive input impedance). These analyses were later refined by other researchers (e.g. \cite{hansen1981fundamental, mclean1996re, gustafsson2007physical}), further underlying the fact that for a linear time-invariant passive antenna one cannot simultaneously optimize size, bandwidth and efficiency.

Following Chu's work, antenna engineers realized very early the need for active (i.e., non-passive) designs, nonlinear materials, and time-varying components to increase the communication bandwidth of wireless systems based on small antennas. In this context, for example, nonlinear ferrite materials were exploited to realize time-varying inductors, which were used to overcome the limitations of small antennas in VLF transmitters by shifting the carrier frequency \cite{jacob1954keying}. Subsequent works (see, for instance, Refs. \cite{wolff1957high, galejs1963switching, yao2004radiating, wang2007time, wang2010analysis}) further demonstrated the unique opportunities of employing time-varying ``switched'' elements (which affect the transient response, and therefore the bandwidth, of radiation) to design antennas with improved characteristics.

More recently, with the renewed surge of interest in temporal modulation in electromagnetic systems, research into time-varying antenna designs have been reinvigorated. For example, Ref. \cite{manteghi2009antenna} confirmed that the antenna bandwidth can be increased significantly through impedance modulation. Another more recent work employed ``Floquet'' impedance matching \cite{li2019beyond} (see Fig. \ref{chu}) to overcome the Chu-Harrington limit by noting that this limitation is fundamentally rooted within the more general Bode-Fano limit (see Section \ref{section_impedance}) \cite{lopez2013historical}. Indeed, as recognized by Chu, a more accurate definition of the fractional bandwidth of the antenna is in terms of the tolerable reflection coefficient over the band (and not simply the inverse of the Q factor), and therefore the antenna bandwidth is limited by the Bode-Fano limit applied to Chu's equivalent circuit for spherical waves. As discussed in \cite{li2019beyond} , this limitation can be surpassed by periodically modulating the matching network components at twice the frequency of the antenna resonance frequency to impart parametric gain \cite{lee2021parametric} into the system. A realistic simulation of the proposed concept was also demonstrated based on a small loop antenna (see Figs. \ref{chu}(b)), which includes periodically modulated variable capacitances (varactors). Through numerical calculations (see Fig. \ref{chu}(c)), it was shown that the bandwidth of the transducer gain can be significantly enhanced owing to the temporal modulation compared to the time-invariant case, and can even surpass the Chu-Harrington limit, which is a rather remarkable result given the fact that parametric phenomena are typically narrowband. A possible drawback of this Floquet scheme is that the parametric gain may potentially lead to instabilities. However, unlike antenna systems based on active non-Foster elements \cite{zhu2012broad,church2014uhf,shi2019improved}, the dispersion of the parametric gain can be easily controlled through the modulation properties to bring the system into the stable regime. Moreover, it was also pointed out that deliberately inducing an impedance mismatch by modifying the source impedance can act as a negative feedback to help ensure stability and also increase the operational bandwidth \cite{li2019beyond}. Indeed, a stable parametric enhancement effect was later experimentally demonstrated in Ref. \cite{mekawy2021parametric} in a fabricated time-varying loop antenna operating at RF frequencies (see Fig. \ref{chu}(d)). The parametric enhancement due to the time-modulated varactor was verified through reflection $|S_{11}|$ measurements (where $|S_{11}|>1$ was observed over a relatively broad bandwidth), which showed good agreement with effective circuit model calculations, as shown in Fig. \ref{chu}(e). The measured data was then used to evaluate the radiated power from the antenna, which was shown to surpass the limitations of linear time-invariant passive antennas \cite{mekawy2021parametric}. In this context, other recent studies also investigated the use of time-varying parametric amplifiers with the aim of broadening the antenna bandwidth \cite{loghmannia2020low, loghmannia2021broadband}. This was achieved by first trading efficiency with bandwidth (in line with the Bode-Fano limit), and subsequently using parametric up-converters based on time-varying reactive elements to increase the antenna gain. It was also noted that such parametric amplifiers have various advantages, including inherent stability and low-noise operation \cite{loghmannia2020low, loghmannia2021broadband}.

\begin{figure}[t]
\centering%
\includegraphics[width=0.9\linewidth]{chu.pdf}
\caption{(a) A ``Floquet'' matching network with periodically time-varying elements can be employed to overcome the Chu-Harrington limit by imparting parametric gain to the system. (b) An electrically small loop antenna (at the design frequency of 300 MHz) matched with periodically-driven variable capacitors (varactors). (c) Transducer gain profiles obtained through numerical simulations show that the bandwidth of the time-modulated antenna (green dotted line) can exceed the maximum bandwidth dictated by the Chu-Harrington limit (red solid line). The time-invariant antenna (blue dashed line), on the other hand, is bounded by the Chu-Harrington limit. The incident frequency is normalized to the antenna design frequency. $\Delta C$ denotes the capacitance modulation amplitude, while $C_0 $ ($\approx$3 pF) is the static capacitance. (d) Photograph of the experimentally realized time-varying loop antenna (left) and the effective circuit model (right). (e) Experimentally measured reflection spectrum superimposed with effective circuit model calculations. Panels (a-c) are reproduced with permission from Ref. \cite{li2019beyond}, American Physical Society. Panels (d-e) are reproduced with permission from Ref. \cite{mekawy2021parametric}, American Physical Society.}
\label{chu}
\end{figure}

\subsection{Electromagnetic Absorption}

\begin{figure}[t]
\centering%
\includegraphics[width=1.0\linewidth]{absorption_sc.pdf}
\caption{(a) An electromagnetic absorber with a time-switched permittivity can be employed to overcome the Rozanov limit for broadband electromagnetic absorption (shown by the green dashed lines). Different reflection spectra are superimposed for incident Gaussian pulses with different FWHM. (b) Alternatively, one may temporally modulate the complex permittivity (both its real and imaginary parts simultaneously) such that the modulation spectrum becomes one-sided (or ``spectrally casual''), which implies that back-reflections due to the temporal modulation are prohibited over a broad bandwidth. (c) Such a complex modulation may potentially increase the bandwidth of absorbers, as in the time-switched case, but without inducing any reflection within or outside the operational bandwidth (as seen in the figure, the field is fully transmitted and absorbed, with no reflections, at any frequency). Panel (a) is reproduced with permission from Ref. \cite{li2021temporal}, Optical Society of America. Panels (b) and (c) are reproduced with permission from Ref. \cite{hayran2021spectral}, Optical Society of America.}
\label{absorption}
\end{figure}

Wave absorption plays a crucial role in various electromagnetic systems, ranging from solar energy harvesting to scattering minimization for radar or anechoic chamber applications. Such systems are designed to have highly absorbing properties within the bandwidth of interest. A key limitation of conventional electromagnetic absorbers originates from the trade-off between the absorption bandwidth and the thickness of the absorbing material. Within this context, using the analytical properties of the reflection coefficient, Rozanov \cite{rozanov2000ultimate} showed that there exists a ultimate thickness to bandwidth ratio $\Delta\lambda/d$ for passive time-invariant absorbers backed by a perfect conductor. Specifically, for the case of a nonmagnetic absorber, the upper limit for this ratio becomes $\Delta\lambda/d<16/|\ln\rho_0|$, where $\rho_0$ is the lowest reflectance within the operational bandwidth. Clearly, the absorption bandwidth can therefore be enlarged by increasing the thickness of the absorber or the tolerable reflection coefficient. Prominent examples of broadband absorbers include cascaded \cite{he2015ultra, yang2021visible} or adiabatically tapered \cite{cui2012ultrabroadband, liang2013metamaterial, hayran2017stopped} structures which typically have thicknesses on the order of several wavelengths. However, for many applications, such as radar absorbing or solar energy harvesting, a thin absorber with a large bandwidth is highly desirable due to practical and economic reasons, which motivates the large interest in overcoming the ``Rozanov bound''. In this direction, it was shown in Ref. \cite{li2021temporal} that temporally switching the relative permittivity of the absorbing material can be used to extend the bandwidth of thin absorbers beyond what is possible with time-invariant structures (see Fig. \ref{absorption}(a)), similar to the use of temporally switched transmission lines for broadband impedance matching as discussed in Section \ref{section_impedance}. An important drawback of these approaches is that, while reflections are reduced within the bandwidth of the incident wave, significant reflections may be induced outside this frequency range as the temporal switching can lead to frequency conversion. While this may not pose a problem for certain narrowband systems, it would become an issue in the context of, for example, broadband radar and stealth technology, as the target may become even more easily detectable due to the increased reflections outside the original signal bandwidth, which could then be picked up by sufficiently broadband or multi-band detectors.  To overcome this problem, a possible strategy is to modulate both the real and imaginary part of the permittivity in such a way that the modulation spectrum becomes ``unilateral'' (or ``spectrally causal''), as we theoretically demonstrated in a recent publication \cite{hayran2021spectral}. Such an approach ensures that back-reflections due to the temporal modulations are prohibited over a very wide frequency range and, therefore, the incident wave is perfectly absorbed with no reflections over wide bandwidths (see Fig. \ref{absorption}(b),(c)). Along similar lines, a recent work theoretically demonstrated, through Green's function analysis and numerical optimizations, that a suitable switching of both the conductivity and permittivity leads to a reduction of reflection over the whole frequency spectrum \cite{firestein2022absorption}. As an alternative strategy, a recent experimental work \cite{yang2022broadband} also demonstrated broadband absorption by creating an energy-trap through switched electronic components triggered by the pulse entering the absorbing region.
While the strategies described above work for broadband pulses and are based on non-periodic temporal modulations (e.g., ``switching''), perfect absorption for quasi-monochromatic signals has also been realized by employing periodic temporal modulations \cite{suwunnarat2019dynamically}. Specifically, it was shown that perfect absorption can be achieved even with infinitesimal losses owing to coherent multichannel illumination with Floquet time-periodic driving schemes. A single-channel coherent perfect absorber based on a time-Floquet engineered reflector without any physical mirror (ground plane) was also demonstrated, which might be useful in certain application scenarios where metallic or Bragg mirrors are not available.

\subsection{Other Applications}

Our list of applications which can benefit from temporal modulation is certainly not exhaustive. In the Supplementary Material, we provide several other application examples that might be of interest to readers also from different backgrounds, such as optics and photonics.

\section{Outlook}
Although the study of time-varying electromagnetic systems dates back more than half a century, recent years have seen a surge of renewed interest into such dynamic schemes, largely motivated by the need to realize novel electromagnetic systems that are superior to conventional LTI ones for various applications, from high-speed full-duplex communication systems, to broadband stealth technology, to high-efficiency energy harvesting and energy transfer. However, despite the promising recent progress in this field, we believe several challenges should be taken into consideration and carefully addressed before such systems can be adopted into practical use.

Since a time-varying system requires external power to operate, such active systems might become susceptible to instabilities (unbounded oscillations) that can disrupt the desired functionality. Therefore, especially for systems involving parametric modulations, a careful analysis is essential to assess and ensure stability \cite{li2019beyond}. Moreover, the additional energy pumped into the system \cite{secondo2020absorptive} should be carefully assessed to determine the figure of merits (such as power efficiency) of the design in a proper manner \cite{jayathurathnage2021time} and also to assess whether the power required by the temporal modulation is beyond practical reach and/or may damage the system.

Frequency dispersion is another fundamental property of materials arising due to causality. Hence, any realistic system involves, to some degree, frequency-variations in the material constitutive parameters. However, especially for numerical studies involving time-varying structures, such dispersion effects are often ignored and the time-variation is assumed to be independent of frequency for simplicity \cite{yu2009complete}. Therefore, any realistic implementation of such proposals will necessarily experience performance deterioration due to frequency dispersion of the material properties. We note that several recent works \cite{mirmoosa2020dipole, solis2021functional, bakunov2021constitutive} have studied material dispersion effects in the context of temporal modulations. We believe further investigations in this area are important to enhance our understanding of time-varying materials and may lead to improved designs with more realistic properties.

For schemes requiring synchronization between the incident wave and the temporal variation, another challenge is that fast detectors and switching mechanisms are essential to trigger the modulation in the presence of the incident wave. Such additional components might complicate the design and reduce its practical applicability. Hence, more research is required to integrate such components into time-varying systems in a practically useful manner.

Dissipation, unless desired to realize absorbers, is another major challenge that needs to be overcome to realize certain practical applications with time-varying components. For instance, epsilon-near-zero materials have recently become the subject of interest in the context of time modulation at optical frequencies, as they can be efficiently controlled through external optical pumps \cite{reshef2019nonlinear}. However, such materials still suffer from high losses especially around their epsilon-near-zero wavelengths. Hence, the search for better epsilon-near-zero materials with lower losses is an important area of research (note that, theoretically, nothing prevents the existence of an epsilon-near-zero material that is dissipationless, albeit frequency dispersive, within a certain frequency window). An alternative strategy could be to exploit (rather than to eliminate or compensate) the losses in such materials to realize time-dependent loss profiles that accompany the time-varying refractive index profiles. In this context, we recently provided a general strategy to exploit the synergy between these two time-modulation profiles \cite{hayran2021spectral}. Investigations along these directions may also enhance our understanding of wave dynamics in complex (``non-Hermitian'') time-varying media.

Despite these challenges, we believe that the use of time-varying materials to overcome the limitations of conventional LTI electromagnetic systems is an exciting and promising research area that might lead to novel and superior devices in the future as more research is conducted. Spatially-varying engineered structures, such as frequency selective surfaces, diffraction gratings, photonic crystals, metamaterials, etc., have revolutionized the field of electromagnetics in recent decades. Along the same vein, it would be safe to say that the new insight and the new tools provided by the field of time-varying electromagnetic systems will enable further advances and possible revolutions in the forthcoming decades for electromagnetic research.

\section*{Acknowledgements}
The authors acknowledge support from the Air Force Office of Scientific Research with Grant No. FA9550-19-1-0043 through Dr. Arje Nachman.

\bibliographystyle{IEEEtran}


\section{Introduction}

Heavy quarks (charm and beauty) have masses much larger than the characteristic energy scale of QCD interactions, $\Lambda_\mathrm{QCD}$. Due to this, they are typically produced in hard scattering processes with large $Q^2$, meaning that their production can be described well using perturbative QCD (pQCD) calculations. Typically the production of heavy-flavour hadrons is factorised into three main components: the parton distribution functions (PDFs) that describe the momentum distributions of quarks in the colliding hadrons; the partonic interaction cross sections to produce the heavy quarks; and the fragmentation functions, which describe the hadronisation of quarks to specific species. The non-perturbative fragmentation functions are determined from experiment in \ensuremath{\mathrm{e}^+\mathrm{e}^-}\xspace and \ensuremath{\mathrm{e}^-\mathrm{p}}\xspace collisions, and are assumed to universally define the hadronisation ratios of charm hadron species in other collision systems. Hadron-to-hadron production ratios are particularly sensitive to changes in the hadronisation mechanisms, since the contributions from the PDFs and scattering cross sections cancel out, making charm baryon measurements an important tool in testing the assumptions of the factorisation approach.

Recent measurements by the ALICE Collaboration of the production of D mesons showed that the strange and non-strange meson-to-meson yield ratios $\ensuremath{\mathrm{D}_\mathrm{s}^{+}}\xspace/\ensuremath{\mathrm{D}^{0}}\xspace$ and $\ensuremath{\mathrm{D}^{+}}\xspace/\ensuremath{\mathrm{D}^{0}}\xspace$ are flat as a function of \ensuremath{p_\mathrm{T}}\xspace~\cite{alice:Dmespp5}, meaning that there are no momentum-dependent modifications to charm meson hadronisation. The ratios are also in line with the expected values from previous \ensuremath{\mathrm{e}^+\mathrm{e}^-}\xspace measurements. In contrast, the \ensuremath{\Lambda_\mathrm{c}^{+}}\xspace/\ensuremath{\mathrm{D}^{0}}\xspace yield ratios in pp and \textrm{p--Pb}\xspace collisions at \ensuremath{\sqrt{s_\mathrm{NN}} = 5.02~\TeV}\xspace are significantly higher than the \ensuremath{\mathrm{e}^+\mathrm{e}^-}\xspace measurements for $\ensuremath{p_\mathrm{T}}\xspace<10\,\gevc$ and exhibit a decreasing trend towards high transverse momentum~\cite{alice:LcpppPb5}. In addition, MC calculations tuned to hadronisation ratios from \ensuremath{\mathrm{e}^+\mathrm{e}^-}\xspace collisions, such as PYTHIA with the Monash tune~\cite{monash}, are unable to reproduce the measurements in pp collisions. Instead, additional hadronisation mechanisms must be considered in order to describe the data. This implies that the previously assumed universality of fragmentation as the main hadronisation process in the factorisation approach does not hold true for baryon production. In order to fully characterise these effects in hadronic colliders, it is important to make precise measurements of further charm baryon states with the widest possible momentum coverage, as well as measuring differentially in terms of other observables such as the event multiplicity.

\section{Results}

The baryon-to-meson yield ratio $\ensuremath{\Lambda_\mathrm{c}^{+}}\xspace/\ensuremath{\mathrm{D}^{0}}\xspace$ is presented down to $\ensuremath{p_\mathrm{T}}\xspace=0$ in \textrm{p--Pb}\xspace collisions at \ensuremath{\sqrt{s_\mathrm{NN}} = 5.02~\TeV}\xspace, and for the first time in pp collisions at \ensuremath{\sqrt{s} = 5.02~\TeV}\xspace, in \figref{fig:lclow} (left). In both collision systems, the measurements in the interval $0<\ensuremath{p_\mathrm{T}}\xspace<1\,\gevc$, indicated by the open circles, were performed in the channel $\ensuremath{\Lambda_\mathrm{c}^{+}}\xspace\to\mathrm{pK}^0_\mathrm{S}$, with candidates reconstructed from their decay tracks using the KFParticle package~\cite{kfparticle} and machine learning selections applied using the XGBoost gradient boosting algorithm~\cite{xgboost}. The measurements for $\pt>1\,\gevc$ are from combined measurements of the decay channels $\ensuremath{\Lambda_\mathrm{c}^{+}}\xspace\to\mathrm{pK}^0_\mathrm{S}$ and $\ensuremath{\Lambda_\mathrm{c}^{+}}\xspace\to\mathrm{pK}\pi$ in both collision systems~\cite{alice:LcpppPb5}. A non-flat distribution is evident in both systems as a function of \ensuremath{p_\mathrm{T}}\xspace, and there is a hint of a hardening in the \ensuremath{p_\mathrm{T}}\xspace spectrum in \textrm{p--Pb}\xspace collisions with respect to pp. In the right panel, the pp measurements are compared with model calculations. Models tuned to hadronisation ratios from \ensuremath{\mathrm{e}^+\mathrm{e}^-}\xspace collisions, such as PYTHIA with the Monash tune~\cite{monash} and HERWIG~\cite{herwig} are unable to describe either the shape or the magnitude of the \ensuremath{\Lc/\Dz}\xspace yield ratio in pp collisions. Models including additional hadronisation effects provide a better description of the results, for instance the Catania model~\cite{catania}, which considers quark coalescence in addition to fragmentation; PYTHIA calculations with enhanced colour reconnection beyond the leading order~\cite{pythCR}; and the SH + RQM model~\cite{rqm}, which is a statistical hadronisation approach including feed-down from yet-unmeasured resonant charm baryon states.

\begin{figure}[h!tb] \centering
\begin{minipage}{0.42\linewidth}
\centerline{\includegraphics[width=0.9\linewidth]{LcDppVspPb_5TeV-lowpT_0}}
\end{minipage}
\begin{minipage}{0.42\linewidth}
\centerline{\includegraphics[width=0.9\linewidth]{LcDpp5vsModels-lowpT_0}}
\end{minipage}
\caption[]{ALICE measurements of the \ensuremath{\Lc/\Dz}\xspace baryon-to-meson yield ratio at \ensuremath{\sqrt{s_\mathrm{NN}} = 5.02~\TeV}\xspace. Left: Comparison of pp and \textrm{p--Pb}\xspace collisions~\cite{alice:LcpppPb5}, including preliminary measurement of the interval $0<\ensuremath{p_\mathrm{T}}\xspace<1\,\gevc$. Right: Comparison between \ensuremath{\Lc/\Dz}\xspace in pp collisions with models~\cite{monash,herwig,alice:LcpppPb5,pythCR,catania,rqm}.}
\label{fig:lclow}
\end{figure}

The \ensuremath{\Sigma_\mathrm{c}^{0,++}}\xspace baryon states are potential contributors to the prompt cross section of \ensuremath{\Lambda_\mathrm{c}^{+}}\xspace baryons, as they decay strongly and are indistinguishable from prompt production. These states are expected to be suppressed with respect to \ensuremath{\Lambda_\mathrm{c}^{+}}\xspace due to the relative rarity of spin-1 diquark states produced in collisions. Measurements of \ensuremath{\Sigma_\mathrm{c}^{0,++}}\xspace production are shown in~\figref{fig:sigc}, in particular the middle panel which shows the $\ensuremath{\Sigma_\mathrm{c}^{0,++}}\xspace/\ensuremath{\mathrm{D}^{0}}\xspace$ production cross section, and the right panel which shows the feed-down contribution from \ensuremath{\Sigma_\mathrm{c}^{0,++}}\xspace to \ensuremath{\Lambda_\mathrm{c}^{+}}\xspace~\cite{alice:Sigc}. A significant increase, by a factor of approximately 10, is seen in the $\ensuremath{\Sigma_\mathrm{c}^{0,++}}\xspace/\ensuremath{\mathrm{D}^{0}}\xspace$ ratio with respect to the Monash tune of PYTHIA, while models with coalescence effects and statistical hadronisation with additional baryon states describe the data well. The fraction $\ensuremath{\Lambda_\mathrm{c}^{+}}\xspace(\leftarrow\ensuremath{\Sigma_\mathrm{c}^{0,++}}\xspace)\times3/2$ is also significant, with a \ensuremath{p_\mathrm{T}}\xspace-integrated value of $0.38\pm0.06\mathrm{(stat.)}\pm0.06\mathrm{(syst)}$. This implies that a sizable proportion of the measured prompt \ensuremath{\Lambda_\mathrm{c}^{+}}\xspace originates from the decays of heavier resonant states, and notably the PYTHIA predictions with enhanced colour reconnection appear to overestimate this contribution.

\begin{figure}[h!tb] \centering
\includegraphics[width=0.82\linewidth]{plot_3pads_new_allModes_LineWidth2models_logXLcD0_newCatania_21_noBRsepareted}
\caption[]{Measurements of (left) \ensuremath{\Lc/\Dz}\xspace, (middle) $\ensuremath{\Sigma_\mathrm{c}^{0,++}}\xspace/\ensuremath{\mathrm{D}^{0}}\xspace$, and (right) contribution of \ensuremath{\Sigma_\mathrm{c}^{0,++}}\xspace feed-down to the prompt \ensuremath{\Lambda_\mathrm{c}^{+}}\xspace cross section in pp collisions at \ensuremath{\sqrt{s} = 13~\TeV}\xspace~\cite{alice:Sigc}.}
\label{fig:sigc}
\end{figure}

The potential role of strangeness in the hadronisation processes may be considered by measuring \ensuremath{\Xi_\mathrm{c}^{0,+}}\xspace and \ensuremath{\Omega_\mathrm{c}^{0}}\xspace baryons. Both have been measured in pp collisions at \ensuremath{\sqrt{s} = 13~\TeV}\xspace, and as an example the $\ensuremath{\Omega_\mathrm{c}^{0}}\xspace/\ensuremath{\Xi_\mathrm{c}^{0}}\xspace$ ratio is shown in the left panel of~\figref{fig:Oc}. As the branching ratio of the decay channel $\ensuremath{\Omega_\mathrm{c}^{0}}\xspace\to\Omega^-\pi^+$ is not yet experimentally known, the \ensuremath{\Omega_\mathrm{c}^{0}}\xspace cross section is not corrected for the branching ratio. This baryon-to-meson yield ratio is underpredicted by two orders of magnitude by PYTHIA calculations using the Monash tune, and by one order of magnitude for PYTHIA including beyond-leading-order colour reconnection. The Catania model in this case comes closest to describing the data, but while the non-strange baryons were described well by the coalesence + fragmentation case, for $\ensuremath{\Omega_\mathrm{c}^{0}}\xspace/\ensuremath{\Xi_\mathrm{c}^{0}}\xspace$ adding an extra contribution of feed-down from resonant baryon states better describes the data.

Differential measurements of baryon and meson production in terms of the charged-particle multiplicity are shown in~\figref{fig:LcDmult} (right). While a significant modification of the \ensuremath{p_\mathrm{T}}\xspace distribution of the \ensuremath{\Lc/\Dz}\xspace yield ratio has been observed between low- and high-multiplicity pp collisions at \ensuremath{\sqrt{s} = 13~\TeV}\xspace~\cite{alice:LcvsMult}, the \ensuremath{p_\mathrm{T}}\xspace-integrated \ensuremath{\Lc/\Dz}\xspace yield ratio spanning three orders of magnitude in multiplicity show no significant dependence on the event multiplicity. This further implies that the collision system dependence of the overall \ensuremath{\Lc/\Dz}\xspace yield ratio between leptonic and hadronic collisions is not simply related to the number of charged particles produced, but rather additional hadronisation processes on top of vacuum fragmentation must be considered for hadronic collisions.

\begin{figure}[h!tb]\centering
\begin{minipage}{0.44\linewidth}
\centerline{\includegraphics[width=0.9\linewidth]{brxomegac0toxic0_pp13tev_wmodel_0}}
\end{minipage}
\begin{minipage}{0.44\linewidth}
\centerline{\includegraphics[width=0.85\linewidth]{LcOverD0_ptint_acrosssystems_wSTAR_fullPythia_0}}
\end{minipage}
\caption[]{Left: $\mathrm{BR}(\ensuremath{\Omega_\mathrm{c}^{0}}\xspace\to\Omega^-\pi^+)\times\sigma(\ensuremath{\Omega_\mathrm{c}^{0}}\xspace)/\sigma(\ensuremath{\Xi_\mathrm{c}^{0}}\xspace)$ baryon-to-meson yield ratio in pp collisions at \ensuremath{\sqrt{s} = 13~\TeV}\xspace. Right: \ensuremath{p_\mathrm{T}}\xspace-integrated \ensuremath{\Lc/\Dz}\xspace yield ratio as a function of charged-particle multiplicity, comparing results from pp, \textrm{p--Pb}\xspace and \textrm{Pb--Pb}\xspace collisions in ALICE with Au--Au collisions at STAR.}
\label{fig:Oc}
\label{fig:LcDmult}
\end{figure}

The measurement of all ground-state charm hadron species in a broad \ensuremath{p_\mathrm{T}}\xspace range allows the total charm production cross section to be derived at mid-rapidity with minimal dependence on models, as shown in \figref{fig:crosssec} (left) as a function of collision energy. The total $\mathrm{c\bar{c}}$ cross section was computed for the first time in \textrm{p--Pb}\xspace collisions at \ensuremath{\sqrt{s_\mathrm{NN}} = 5.02~\TeV}\xspace and shown as the open blue circle, along with previous ALICE measurements~\cite{alice:ppfrag} in pp collisions at $\ensuremath{\sqrt{s}}\xspace=2.76, 5.02,$ and $7\,\ensuremath{\mathrm{TeV}}\xspace$ shown as solid points. The results are compared with pQCD calculations under the FONLL and NNLO schemes, and the ALICE measurements are consistently seen to lie at the upper edge of the theory uncertainty bands. The cross section can also be split into the individual contributions from hadron species to give the relative hadronisation fractions, $f(\mathrm{c}\to\mathrm{H}_\mathrm{c})$, including the first measurement of the $f(\ensuremath{\Xi_\mathrm{c}^{0}}\xspace)$. The measurements in pp collisions~\cite{alice:ppfrag} and preliminary measurements in \textrm{p--Pb}\xspace collisions are shown in the right panel of~\figref{fig:crosssec}, in comparison with measurements in \ensuremath{\mathrm{e}^+\mathrm{e}^-}\xspace and \ensuremath{\mathrm{e}^-\mathrm{p}}\xspace collisions. The two hadronic collision systems are consistent with one another, but a significant enhancement of \ensuremath{\Lambda_\mathrm{c}^{+}}\xspace and depletion of \ensuremath{\mathrm{D}^{0}}\xspace production are seen with respect to leptonic collisions.

\begin{figure}[h!tb]\centering
\begin{minipage}{0.43\linewidth}
\centerline{\includegraphics[width=0.9\linewidth]{cc_Vs_energy_withpPb_0}}
\end{minipage}
\begin{minipage}{0.43\linewidth}
\centerline{\includegraphics[width=0.9\linewidth]{FF5tev_withpPb}}
\end{minipage}
\caption[]{Left: Total charm production cross section at mid-rapidity in pp and p--Pb collisions as a function of collision energy. Right: Relative hadronisation fractions of ground-state charm hadron species measured by ALICE in pp and p--Pb collisions, compared with \ensuremath{\mathrm{e}^+\mathrm{e}^-}\xspace and \ensuremath{\mathrm{e}^-\mathrm{p}}\xspace collisions.}
\label{fig:crosssec}
\end{figure}

Looking to the upcoming Run 3 of the LHC and beyond to ALICE 3, the higher rates of data taking and upgrades to detector systems will allow the precision and \ensuremath{p_\mathrm{T}}\xspace coverage of charm baryon measurements to be greatly improved as well as experimental access to new species, e.g. multi-charm baryons such as $\Xi_\mathrm{cc}^{++}$.

\section*{References}