Non-zero temperature study of spin 1/2 charmed baryons using lattice gauge theory

We study the behaviour of spin $1/2$ charmed baryons as the temperature increases. We make use of anisotropic lattice QCD simulations with $N_f = 2 + 1$ dynamical flavours. After determining the positive and negative parity ground state masses at the lowest temperature, we investigate the effect of rising temperature using ratios of thermal lattice correlators with both so-called reconstructed correlators and with simple model correlators. This avoids difficulties associated with non-zero temperature fitting or spectral reconstruction. We find that temperature effects are prominent throughout the hadronic phase for all negative parity channels considered and for some positive parity channels. Subsequently and where possible, we determine the masses of the ground states as a function of temperature. Finally we consider the effect of chiral symmetry restoration and extract an estimate of the pseudocritical temperature from singly charmed baryonic correlators.


Introduction
This work focuses on the fate of spin 1/2 charmed baryons as the temperature increases and the confining hadronic medium at low temperature smoothly transitions to a quark-gluon plasma (QGP), as predicted by Quantum Chromodynamics (QCD) [1].While the light degrees of freedom deconfine at the crossover [2,3], there is evidence that heavier hadrons comprised of charm and bottom quarks can survive as bound states in the QGP.Since the seminal work by Matsui and Satz [4], this has been studied extensively for charmonium (cc) and bottomonium (b b) states, as it provides important insights into the length scales in the QGP [5][6][7].For results obtained with lattice QCD, see e.g. the comprehensive reviews [8,9] and references therein.
Here we consider baryons containing charm quarks, building upon previous work studying the behaviour of light and strange baryons at non-zero temperature [10][11][12].The focus in this earlier work was parity doubling, the emergence of a degeneracy between positive-and negativeparity correlators and ground state masses, which is a signal for chiral symmetry restoration.Due to the heav-ier mass of the charm quark, one expects a reduction in the level of parity doubling observed when compared to the light and strange baryons.Instead, baryons containing multiple charm quarks may be more similar to heavyquark mesons (quarkonia) and therefore not immediately dissolve in the QGP.Hence, it will be interesting to analyse the structure of charmed baryons across a range of temperatures.
Heavy-light charmed D and D s mesons have been studied using thermal lattice QCD in Refs.[13,14].In particular, Ref. [14] considered the response of D (s) mesons to an increasing temperature without the explicit need for non-zero temperature fits and spectral reconstruction.The outcome of this analysis was subsequently used to inform the applicability of standard fits to the correlation function used to determine the ground state mass.The analysis studied ratios of lattice correlators to singlestate model correlators to determine when the spectral content of the correlator was different from that of the zero-temperature correlator.Here we extend this treatment to the baryonic sector for the first time.
Motivated and directed by the outcome of the analysis of correlator ratios across temperatures, we fit correlators arXiv:2308.12207v2[hep-lat] 20 Mar 2024 with standard exponential fits when ratios suggest that thermal effects are small.Hence these fits assume that the hadron is a (narrow) bound state.Extensive use is made of model averaging methods [15,16] in order to produce a reliable extraction of the mass.It is noted that the behaviour of charmed baryon masses at non-zero temperature may be of phenomenological interest [17].
To make use of the full range of temperatures available, we also compare positive and negative parity channels using the correlation functions directly [11,12,18], to reveal the effect of chiral symmetry restoration in the QGP.Parity doubling is not expected due to the large charm quark mass; nevertheless, for singly charmed baryons a clear change of behaviour is observed.
The paper is organised as follows.In Section 2 we summarise the FASTSUM Generation 2L ensembles, correlator methods, and the fitting techniques used.We consider ratios of the lattice correlator with the "reconstructed" correlator and with a single-state model correlator in Section 3. In Section 4 we determine masses across a range of temperatures provided the analysis of the previous section suggests this is appropriate.The fit-independent parity doubling ratio is presented in Section 5. Conclusions and future work are discussed at the end.

Lattice Correlators
In the following subsections we discuss the gauge field ensembles, operators, and fitting procedure used to determine ground state masses at the lowest temperature.

Ensembles
We make use of the thermal ensembles of the FASTSUM collaboration [18], with 2 + 1 flavours of Wilson fermions on anisotropic lattices.The renormalised anisotropy is ξ ≡ a s /a τ = 3.453 (6) [18,20].The lattice action follows that of the Hadron Spectrum Collaboration [21] and is a Symanzikimproved [22,23] anisotropic gauge action with tree-level mean-field coefficients and a mean-field-improved Wilsonclover [24,25] fermion action with stout-smeared links [26].Full details of the action and parameter values can be Table 1.FASTSUM Generation 2L ensembles used in this work.The lattice size is 32 3 × Nτ , with temperature T = 1/ (aτ Nτ ).The spatial lattice spacing is as = 0.11208 (31) fm, renormalised anisotropy ξ = as/aτ = 3.453 (6) and the pion mass mπ = 239(1) MeV [19].We use ∼ 1000 configurations and eight (random) sources for a total of ∼ 8000 measurements at each temperature.The estimate for Tpc comes from an analysis of the renormalised chiral condensate and equals Tpc = 167(2)(1) MeV [14,18].Full details of these ensembles may be found in Refs.[14,18].found in Ref. [18].We use the "Generation 2L" ensembles which have a pion mass of m π = 239(1) MeV (in the previous "Generation 2" ensembles the pion mass was m π = 384(4) MeV [27][28][29]).While this is still heavier than physical, it represents an important step towards the physical regime.The strange quark has been tuned to its physical value via the tuning of the light and strange pseudoscalar masses [27,30,31].The ensembles are generated using a fixed-scale approach, such that the temperature is varied by changing N τ , as T = 1/ (a τ N τ ).A summary of the ensembles is given in Table 1.There are five ensembles below the pseudocritical temperature T pc = 167(2)(1) MeV, one close to T pc and five above T pc .The estimate for T pc comes from an analysis of the renormalised chiral condensate [18].Note that here we have used the updated lattice spacing of Ref. [19], which has been implemented in our analysis in Ref. [14].

Operators
The J = 1/2 singly and doubly charmed baryons, with quark content, are grouped according to the underlying flavour symmetry of QCD, namely and for doubly charmed baryons, following the decomposition of SU (4) as in the PDG 2020 [32].
The baryon correlators we consider are of the form where O = O † γ 4 and α, α ′ are Dirac indices.The operators used follow Refs.[33,34].Denoting the three quarks as q 1 , q 2 and q 3 (from left to right), the operator for the case q 1 = q 2 (Ξ cc , Ω c and Ω cc ) is given by For the non-degenerate case, Σ c and Ξ ′ c belong to the SU (3) 6 flavour multiplet, with operators where q 1 = u, q 3 = c and q 2 = d/s respectively.Λ c and Ξ c enjoy SU (3) 3 flavour symmetry, with operators where q 1 = u, q 3 = c and q 2 = d/s respectively.The Euclidean gamma matrices satisfy γ and C is the charge conjugation matrix [35].
The contractions of these operators produce terms proportional to where S 1,2,3 is a quark propagator of flavour q 1,2,3 and Γ so/si is the gamma matrix in the source/sink operator respectively.These are the principal building blocks of the baryon correlation functions.The positive and negative parity projectors are [33,35] and we denote the projected correlation functions at vanishing spatial momentum as [10][11][12] G ± (τ ) = tr P ± G(τ ).
These are related as [11] G implying that the forward-(backward-) propagating states of G + (τ ) are states with positive (negative) parity.Gaussian smearing is applied to the source and sink using [36] where η is the bare (delta function) source, κ and n determine the amount of smearing, H is the spatial hopping part of the Dirac operator and C norm is an appropriate normalisation.Here we used κ = 5.5, n = 100 at all temperatures.The root-mean-square radius of this profile is ∼ 6.8 lattice sites.These parameters were chosen such that the positive parity nucleon ground state at the lowest temperature displays good ground state isolation, as in Refs.[11,12].

Mass extraction
A systematic extraction of hadron masses from lattice QCD is an area of active development, with methods becoming increasingly sophisticated [15,16,[37][38][39][40][41][42][43].Broadly these can be divided into methods which improve the underlying correlation function [37,38,41] and those which aim to more reliably extract the mass from a given correlation function.As mentioned above, we use gauge invariant Gaussian source and sink smearing in order to improve the overlap of the operator with the (zero-temperature) ground state.We also adopt the model averaging methods of both Refs.[16,43] and Ref. [15].By comparing these two distinct methods we ensure a robust determination of the mass.
We start at the lowest temperature (T = 47 MeV, N τ = 128) and consider fit functions of the form where G(τ ) = G + (τ ) /G + (0) and M ± n is the n th positive (negative) parity state.The number of exponentials N is allowed to vary from one to three.The final number of exponentials is set by examining the Gaussian Bayes factor and stepping back one exponential when it ceases to change as more exponentials are added [39,[44][45][46].This is equivalent to stopping when the data can no longer support additional exponential terms.We illustrate some of our findings below using the positive and negative parity Ξ cc ground states as examples.
Typically the correlator G (τ ) is considered from a starting time τ min up until an end time τ max .This is due to excited state contamination at early times and the onset of noise of late times (on the positive parity side and trivially reversed for the negative parity side).The choice of "fit window" can have a substantial effect on the extracted mass; this is particularly evident for the effective mass [35] fits shown in Figure 1 for the positive parity Ξ cc ground state at the lowest temperature.Here we have fit the effective mass, with δ τ = 2 across an interval τ min ≤ τ ≤ τ max , with varying τ min .Specifically, we used τ min ≥ 0.1 N τ = 13 and τ max ≤ N τ /2 − 0.05 N τ = 57 (note that in this section we use "temporal lattice units" for both Euclidean time and quoted mass values).This range was chosen to remove early time slices dominated by excited state effects and "noisy" time slices near the middle of the time extent.At small values of τ max clear excited state contamination is manifest in the increased value from the fit.As τ min and τ max are increased, excited state contamination is decreased as the ground state becomes more dominant, but the statistical uncertainty also increases.While it may be possible to extract a robust ground state mass from the effective mass for N τ = 128, in general this will not be possible at higher temperatures, due to a reduction in the number of temporal points.Exponential fits electing a single fit window ignore the information from other fit windows.It is therefore beneficial to use model averaging methods to determine weighted fit parameters.The first method, introduced in Ref. [16] uses a weight w(1) where p f is the p-value of the fit f , δM f is the uncertainty in the fit parameter of fit f and there are N fits to be averaged.This method penalises both poor fits and unconstraining fits [42].
Fig. 2. Mass of the positive parity Ξcc ground state at the lowest temperature, extracted from multi-exponential fits as in Eq. ( 12).For each τmin, τmax increases from left to right in each pane.The top 30% highest weighted fits from Eq. ( 15) are shown, as are the corresponding model averaged values from Eqs. ( 14) and ( 15).X-axis tick labels are not shown where they would overlap.
This is in contrast to the second method [15] which uses an exponential of the modified Akaike information criterion [15,47] where N cut is the number of data points not fit to, k is the number of fit parameters and χ 2 aug is the augmented chi-squared [44].Agreement between the two methods provides confidence in the extracted fit parameter.
A large number of roughly equivalent fit windows is shown in Figure 2. The exponential fits have range τ min = 1 ≤ τ ≤ τ max ≤ N τ /2 − 0.05 N τ .In particular, in Figure 2 only the top 30% highest weighted fits using the weight of Eq. ( 15) are shown.
To exclude outlying fits, a further cut is made.The top 30% of weighted fits using the modified Akaike information criterion weights of Eq. ( 15) are separated and re-averaged using both methods.This is demonstrated in Figure 3, in which the effective masses for both the positive and negative parity Ξ cc ground state are shown and compared to the model averaged results from the multiexponential fits after this cut has been applied.The resulting fits are in good agreement for both the sectors giving us assurance that the fitted masses are robustly determined.
In practice, to produce a single number for a mass, a single averaged value is chosen by examination of the weights, such that the averaged value clearly represents the fits.This can be seen in Figure 4 where the averaged value encompasses the majority of the fit weights (which sum to 1).The difference to the other three methods of Eqs. ( 14) and (15) with and without the top 30% fit cut is then added as a systematic uncertainty.An unrepresentative case would have the weight dominated by a single, Table 2. Positive and negative parity ground state masses M ± at the lowest temperature, in units of aτ and in GeV.The systematic error due to the choice of model averaging method has been added in quadrature to the error from the selected model averaging procedure.The final two columns show the Particle Data Group values [48], when known, averaged over charge partners with systematic uncertainties added in quadrature.
Fig. 3. Comparison of the effective mass for the positive and negative parity Ξcc ground state at the lowest temperature with the results of the two model averaging methods of Eqs. ( 14) and (15).
outlying fit.If this is the case, that method is not used when determining the systematic uncertainty.

Spectrum
We are now in a position to present the spectrum of J = 1/2 charmed baryons with positive and negative parity at the lowest temperature, in Figure 5 and Table 2.While exact agreement with the experimental results is not expected due to our heavier-than-physical light quarks and lack of continuum and infinite volume extrapolation, there is encouraging similarity in the observed masses in both the positive and negative parity channels.Better agreement is observed for hadrons with more strange and charm quarks; this can be explained by the tuning of these heavier quarks to their physical masses [27,30,31].We remind the reader that we do not attempt precision zerotemperature spectroscopy, see e.g.Refs.[30,[49][50][51][52][53][54][55][56][57][58][59][60][61][62] instead, but that our interest is in the response to heating up the system, to which we turn in the next section.15) are also shown, as is the resulting model averaged fit value.

Temperature Effects
To investigate the effect of increasing the temperature of the hadronic medium, we first attempt to assess this using only information contained in the correlators, i.e., without performing fits at finite temperature.We use two complementary approaches; first we perform an appropriate resummation of the correlator at the lowest temperature to account for the difference in temporal extent at finite temperature.This leads to the so-called "reconstructed" correlators [63], extended to the fermionic case here.Subsequently, we consider a minimal model of the thermal correlator using information only from the well-defined ground state fits at the lowest temperature and take appropriate ratios, following Ref.[14].In both these approaches, the key idea is that we attempt to isolate the effect of changing the temperature on the spectrum, while removing the effect of reducing the temporal extent as a trivial kinematic construction on the Euclidean lattice.

Reconstructed Correlator
The spectral relation for baryons reads, see e.g.Eqs.(2.16, 2.17, 2.27) of Ref. [11], where K F is the fermionic kernel Here G (τ ; T ) and ρ (ω) can have Dirac indices, or be projected with parity operators, see Ref. [11].This kernel can be compared with the building block for bosonic (mesonic) correlators Often this is multiplied with exp(ω/2T ) in the numerator and the denominator, to arrive at the usual sinh(ω/2T ) in the denominator, but this is not needed here.
To relate a correlator at a higher temperature T to one at a lower temperature T 0 by resummation, we use a simple identity which follows from factorisation.We switch to lattice units such that T = 1/N τ , T 0 = 1/N 0 , and N 0 /N τ = m is an (odd) integer.In the bosonic case where m is any integer.The fermion case is for m an odd integer.The bosonic kernel can hence be written as for integer m, and the fermionic kernel as for integer odd m.Inserting Eq. ( 22) into the spectral relation ( 16) for a fermionic correlator at temperature T = 1/N τ relates this correlator to a resummation of one at a lower temperature T 0 = 1/N 0 = 1/(mN τ ), assuming that the spectral content content is unchanged.Thus yields the reconstructed correlator for fermions, Switching back to denoting the temperatures with T and T 0 , it allows us to compare the actual correlator at temperature T , G(τ ; T ), with the expected correlator if the spectral content is unchanged from a lower reference temperature T 0 , namely G rec (τ ; T, T 0 ), via the ratio For bosons, this technique was used in e.g.Refs.[13,63].
For fermions, we are not aware of a previous application.As m = N 0 /N τ = T /T 0 must be an odd integer for fermions, the lattice sizes which can be used are limited in principle.This can be resolved by "padding" the correlator with the value at the minimum.No qualitative difference is found by padding with zero instead.We consider only adding points to extend the correlator as this can be done at all temperatures.In principle one could also consider removing points from the reference correlator at the lower temperature T 0 = 1/N 0 .No substantial difference is observed in our tests.

Model Correlator
Alternatively, we may construct simple model correlators at a higher temperature T by employing the ground state masses determined at the lowest temperature T 0 in the previous section, built on the assumption that spectral content has not changed and using the baryon spectral relation (16) for the zero-momentum projected positiveparity correlator.We write for the spectral function where M ± 0 is the positive (negative) parity mass at temperature T 0 and A ± are the corresponding amplitudes.The model correlator at temperature T is then The main assumption here is that the width of the state extracted is negligible.We can now take ratios of the actual and the model correlator, as well as double ratios, to potentially eliminate the effects of excited states and other features at early Euclidean times.This approach has been studied in some detail in Ref. [14] for the case of D (s) mesons.Note that in all ratios we may consider separately the positive parity correlator G + (τ ) and the negative parity correlator G − (τ ) = −G + (1/T − τ ).Moreover, we always divide the correlators by their value at τ = 0, to set a consistent scale.For further analysis on the reconstructed and double ratios, see Ref. [64].

Results
We now present some results for ratios of correlators, using both the reconstructed correlator and the model correlator in the denominator.We focus here on the Σ c (udc) channel; Appendix A contains additional results in the Ξ cc (ccu) and Ω cc (ccs) channels.These are representative; we have of course analysed all channels.We start with the ratio using the reconstructed correlators, see Eq. ( 24).In Figure 6 we show the result in the Σ c (udc) channel, for positive parity (left) and negative parity (right), at low temperatures (top) and at higher temperatures (bottom).Note that T = 127 MeV is shown in both sets of figures.In the positive parity case, we note that the ratio is close to 1 at the lower temperatures; the bending upwards around τ /a τ = N τ /2 is expected due to the appearance of states with opposite parity.At the higher temperatures, the deviation becomes more significant around T ∼ T pc and is of the order of 10% at T = 190 MeV.Note the different vertical scales.We conclude that in the positive parity sector thermal effects are at the percent level up to T ∼ T pc .In the negative parity case on the other hand, there are strong temperature effects at all temperatures, even deep in the hadronic phase.
The reconstructed correlator gives insight into possible changes into the underlying spectral functions as the temperature increases.A closer look at the fate of the ground state can be obtained by taking the ratio with the model correlator, see Eq. ( 27).This analysis relies on the assumption that at the lowest temperature the width of the state extracted is negligible and that it is possible to distinguish the ground state from the excited states.As the temperature increases the temporal extent of the lattice decreases and so the ground state may not be distinguishable from the excited states, even if the spectral information is unchanged.If the ground states are precisely determined and indeed independent of temperature, the ratio Eq. ( 27) will be independent of τ for large enough Euclidean time, understood as the region where the ground state dominates, but where the "bending" of the correlator near N τ /2 is not yet encountered.As the model contains only the ground state and the correlator G (τ ; T ) has excited states present, this ratio is not expected to be completely independent of τ .This is demonstrated in Figure 7, again in the Σ c (udc) channel.On the positive parity side, we do indeed observe a horizontal plateau for T ≲ 127 MeV, indicating that the ground state is well described by a single exponential.The negative parity sector is less clear due to the skewness of the correlation function causing fewer time slices to be available before the mixing of forward and backward states becomes evident.This could be improved by using methods to reduce the amount of excited state contamination in the correlator.The single ratio is similar to an effective mass in that it shows when excited state effects cannot be ignored.The presence of plateaus in the single ratio is indicative that the narrow peak ansatz for the spectral function is suitable for exponential fits as in Eq. (12).At our lowest temperature, T = 47 MeV, excellent plateau behaviour is seen for all channels investigated, providing further support for the results for the ground states masses presented in Table 2.
From the plots with single ratios, it is clear that also the structure at earlier Euclidean times is approximately temperature independent, in particular at lower temperatures and in the positive parity sector, see Figure 7 (top left).We can factor out the common behaviour by considering the double ratio of Eq. ( 28), involving temperatures T and T 0 .The results are shown in Fig. 8.We note that the magnitude of deviation from 1 is comparable between the double ratio and the ratio with the reconstructed correlator.Quantitatively, they show different behaviour: as mentioned, the double ratio only uses the information on the ground state at T = T 0 , assuming that it is well described by a narrow peak, whereas the reconstructed correlator uses the full spectral function at T 0 .
To determine where exponential fits can be reasonably performed to extract the ground state mass in the thermal case, we will require plateau-like behaviour with the double ratio not exceeding a difference of ∼ 10% from one.This is illustrated in Fig. 8, where evidence can be seen   of a change in behaviour around the pseudocritical temperature in the positive parity channel.The double ratio is close to one until T = 169 MeV ∼ T pc .The negative parity sector is much noisier -note the different scale -and deviates from one at much lower temperatures, removing the justification of using exponential fits for the negative parity state above T = 127 MeV in this channel.
When the negative parity sector has clear evidence of a change in spectral content, we do not report a negative parity mass obtained by fitting to an exponential Ansatz, as it is not justified to do so.As mentioned, Appendix A contains additional results in the Ξ cc (ccu) and Ω cc (ccs) channels.

Masses at non-zero temperature
We now extend the analysis detailed in Section 2.3 to nonzero temperature.We use the analysis and methods of the previous sections as a guide to when the spectral content has changed such that the exponential Ansatz of Eq. ( 12) is no longer appropriate.Consequently we show mass results for only a subset of the temperatures available to us.The fitting method used is identical to those in Section 2.3; the ratio analysis of the previous section provides guidance up to which temperature such an analysis can be applied.The singly charmed baryon spectrum is presented in Figure 9.All masses have been normalised with the positive parity ground state mass at the lowest temperature, m + (T 0 ), in the corresponding channel.While the uncertainties increase as temperature does -reflective of fewer temporal points and a potential shift away from a zerowidth state -the pattern suggests that in the positive parity sector the masses increase with temperature and that this starts relatively deep in the confining phase.On the negative parity side, acceptable results are obtained, even though the negative parity signal is shorter lived than the positive parity, as illustrated in e.g. Figure 3.Here the increase in uncertainties is such that it is difficult to determine if masses of negative parity states are changing in a systematic way.We note that the ratio analysis of Section 3 did indeed suggest that the negative parity sector is more strongly affected by temperature effects and as such we do not report masses at as high a temperature for the negative parity states as compared to the positive parity ones.
Masses of the doubly charmed baryons Ξ cc (ccu) and Ω cc (ccs) are shown in Figure 10.Temperature effects are much less evident here, in both parity sectors, with the positive parity Ξ cc (ccu) mass remaining approximately constant well into the deconfining phase, up to T ∼ 190 MeV.To provide support for this finding we show the single and double ratios for this baryon in Figures 13 and  15 in Appendix A. The positive parity double ratio in Figure 15 is particularly striking, with a value close to one at all temperatures shown here, indicating the apparent absence of thermal effects in this channel.It would be interesting to further investigate the difference in thermal effects for doubly charmed baryons compared to singly charmed baryons using effective potential models, similar to the well-studied meson (quarkonium) cases; preliminary work in this direction has been carried out [65].
While there is a possibility that new continuum states emerge in each of the channels considered at or near each respective ground state and conspire to keep the (double) ratio mostly unchanged, we consider this unlikely to occur in each channel.The simpler solution is a modification of the ground state masses.At a single volume, or without spectral function reconstruction methods, it is not possible for us to eliminate the possibility of new continuum states.

Parity doubling
The restoration of chiral symmetry is a signature for the formation of a QGP [2,3,66,67].We investigate chiral symmetry restoration via the parity doubling phenomenon [11], in which positive and negative parity correlation functions become degenerate when chiral symmetry is unbroken.While correlators for light and strange baryons exhibit a clear parity doubling signal [11,12], this is not expected here as the charm quark mass strongly and explicitly breaks chiral symmetry.Nevertheless, the light quarks contained within the charmed baryons are susceptible to chiral symmetry, resulting in a restoration signal as we now demonstrate.We define the baryon R parameter [10-12, 18, 68] from the summed ratio of positive and negative correlators, Here G + (τ ) is the positive parity correlator, G + (1/T − τ ) = −G − (τ ) is the negative parity correlator -see Eq. ( 10), σ R (τ n ) denotes the statistical error for R (τ n ) and we consider a sum over time slices τ n /a τ ∈ [n 0 , N τ /2 − 1] at all temperatures.We choose n 0 = 4 to suppress lattice and excited state artefacts at small τ n [10].Shifting n 0 slightly does not have a qualitative effect on the results and Gaussian smearing at the source and sink further suppresses the effect of excited states.Lattice artefacts are due to e.g. the chiral symmetry violating Wilson term at larger energy scales.As suggested by Eq. ( 29), the R parameter describes the difference between the positive and negative parity correlators.When chiral symmetry is unbroken, the correlators are degenerate and R = 0. Broken chiral symmetry produces R ̸ = 0.In the limit where the correlators are dominated by their respective ground states and the positive parity mass is much smaller than the negative parity mass, R ≈ 1.We hence expect that R is close to one in the hadronic phase.Since the charm mass breaks chiral symmetry at all temperatures considered, we do not expect R to go to zero.
Figure 11 shows the R parameter for the singly and doubly charmed J = 1 2 baryons we considered.We observe that the values are indeed close to one at the lowest temperature, as expected.As the temperature increases, R decreases for singly charmed baryons, indicating the sensitivity of the constituent light quarks to chiral symmetry restoration.Above the crossover temperature, the R values are distinct from zero but continue to decrease.It is expected that at very high temperature the effect of the charm quark mass will vanish eventually, as m c /T → 0. For doubly charmed baryons, the R parameter is much less The lines connecting the data points in Figure 11 are cubic splines fitted to the points.For singly charmed baryons, these splines have an inflection point indicated by the vertical lines and the inflection point temperatures are listed in Table 3.We note that within a few MeV the inflection point temperatures agree with the pseudocritical temperature obtained from the renormalised chiral condensate.This approximate agreement was observed earlier for the case of light and strange baryons, see Refs.[12,18].It is somewhat remarkable to note that even for singly charmed baryons the constituent light and strange quarks are sufficiently sensitive to chiral symmetry restoration to induce inflection points at these temperature values.On the other hand, the R parameters for doubly charmed baryons do not show an inflection point in this temperature range.

Conclusions
In this work we presented a detailed study of the spin 1/2 charmed baryons throughout the hadronic and into the deconfining phase.After determining the spectrum at the lowest temperature, the temperature dependence in each channel was investigated using ratios of thermal correlators with model or reconstructed correlators built on spectral content from the lowest temperature.This permits us to identify where thermal effects become important.It was found that the negative parity sector is sensitive to thermal effects at a much lower temperature than the positive parity sector.For doubly charmed, positive parity baryons we observed an approximate temperature independence of correlator ratios up to T ∼ 190 MeV, the highest temperature we studied.The analysis of these ratios is robust as it does not rely on a particular Ansatz for spectral content at finite temperature.
Based upon the analysis of correlator ratios, we justify the use of multi-exponential fits to determine baryon masses in the cases where the spectral content appears to be largely unchanged.This allows us to extract positive parity masses at temperatures up to T ∼ 190 MeV, albeit with a large uncertainty for singly charmed baryons, while the negative parity masses can only be determined up to T ∼ 127 MeV.The masses of the positive parity ground states of doubly charmed baryons are approximately independent of temperature up to T ∼ 190 MeV.
To study the effect of chiral symmetry restoration for charmed baryons, we investigated the R parameter which encodes the difference between positive and negative parity correlators as a function of temperature.While parity doubling is not manifest -as expected due to the large charm quark mass -a crossover effect is observed for singly charmed baryons.Interestingly, the temperatures of the inflection points of this quantity are close to the pseudocritical temperature obtained from the renormalised chiral condensate, indicating that the constituent light and strange quarks are sufficiently sensitive to chiral symmetry restoration to induce this effect also in charmed correlators.
Future studies could examine the spectrum using either a more sophisticated operator basis or a lattice with more temporal points [69].Both these methods are under investigation.Alternatively one could use spectral function reconstruction methods to directly examine the mass and width of the ground states as a function of temperature to determine when thermal effects become significant.
The analysis in this work makes ample use of the Python packages gvar [72] and lsqfit [46].Additional data analysis tools included matplotlib [73,74] and numpy [75].The dataset and scripts used for this paper can be found at [76].Fig. 12. Ratio with the reconstructed correlator as in Figure 6, in the Ωcc (ccs) channel.

Authors' Contributions
-Bignell: Data analysis, plot generation and primary manuscript production.-Jäger: Data production and preliminary data analysis.
-Aarts: Initial proposal [77], development of the reconstructed correlator for baryons and draft of the manuscript.-Allton, Anwar, Burns, Skullerud: Contributions to analysis and physics interpretations of results, and draft of the manuscript.

A Ratios in other channels
For comparison, we present the ratios of thermal correlators with reconstructed and model correlators in the Ξ cc (ccu) and the Ω cc (ccs) channels in Figs.12-16.Details are as in Section 3.

Fig. 4 .
Fig.4.Highest weighted fit results for the positive and negative parity Ξcc ground state at the lowest temperature as a function of fit window.Both sides of the correlator are fit symmetrically and so e.g."20-57" fits points[20,57] and[71, 108]   simultaneously.The weights w(2) f of Eq. (15) are also shown, as is the resulting model averaged fit value.

Fig. 5 .
Fig. 5. Positive (diamonds ♢) and negative (stars ⋆) parity ground state masses M ± at the lowest temperature.Inner error bars represent the statistical uncertainty and the outer bars incorporate a systematic error from the choice of averaging method.The grey bars represent experimental results from the Particle Data Group, when known [48].

Fig. 9 .
Fig.9.Ground state masses of singly charmed spin 1/2 baryons, normalised with the positive parity ground state mass at the lowest temperature, as a function of temperature.Filled (open) symbols are used for positive (negative) parity states.The inner error bar represents the statistical uncertainty and the outer incorporates the systematic from the choice of averaging method.Horizontal dashed lines show the result from the lowest temperature.The uncertainty of the lowest temperature reflects the relative uncertainty of the lowest temperature mass.Masses are shown only when the analysis of Section 3 suggests the mass can be extracted using an exponential Ansatz.

Fig. 11 .
Fig. 11.Crossover behaviour of the R parameter of Eq. (30) for J = 1 2 baryons.The lines connecting the data points are cubic splines; the vertical lines indicate the inflection points for singly charmed baryons.No inflection point is found for the doubly-charmed baryons.

Table 3 .
(30)erature of the inflection points (in MeV) of the R parameter(30)for singly charmed baryons presented in Figure11, categorised by strangeness (S) and charm (C).Uncertainties are statistical only.For reference, the pseudocritical temperature from the renormalised chiral condensate is Tpc = 167(2)(1) MeV.
G.A., C.A., R.B. and T.J.B. are grateful for support via STFC grant ST/T000813/1.M.N.A. acknowledges support from The Royal Society Newton International Fellowship.This work used the DiRAC Extreme Scaling service at the University of Edinburgh, operated by the Edinburgh Parallel Computing Centre and the DiRAC Data Intensive service operated by the University of Leicester IT Services on behalf of the STFC DiRAC HPC Facility (www.dirac.ac.uk).This equipment was funded by BEIS capital funding via STFC capital grants ST/R00238X/1, ST/K000373/1 and ST/R002363/1 and STFC DiRAC Operations grants ST/R001006/1 and ST/R001014/1.DiRAC is part of the UK National e-Infrastructure.We acknowledge the support of the Swansea Academy for Advanced Computing, the Supercomputing Wales project, which is part-funded by the European Regional Development Fund (ERDF) via Welsh Government, and the University of Southern Denmark and ICHEC, Ireland for use of computing facilities.This work was performed using PRACE resources at Cineca (Italy), CEA (France) and Stuttgart (Germany) via grants 2015133079, 2018194714, 2019214714 and 2020214714.