Entanglement and Bell inequalities violation in H → ZZ with anomalous coupling

,


I. INTRODUCTION
Violation of Bell inequalities is one of the most striking properties of quantum theory.Such a violation has been observed in a variety of physical systems like e.g.pairs of photons [1][2][3], ions [4], electrons [5], superconducting currents [6] or solid state systems [7].
Recently, the possibility of observing quantum entanglement and violation of Bell-type inequalities in high energy physics has been put forward.In particular, scattering processes [8,9], systems of top quarks [10][11][12][13][14][15], B 0 B0 mesons [16] and pairs of vector bosons arising from Higgs particle decay [17][18][19][20][21][22] were proposed in this context.This last possibility for the first time was suggested by Alan Barr in [17].Barr considered there the violation of Clauser-Horn-Shimony-Holt (CHSH) and Collins-Gisin-Linden-Massar-Popescu (CGLMP) inequalities in a system of W W bosons arising in the decay of Higgs particle.In [18] the possible violation of CHSH, Mermin and CGLMP inequalities for a boson-antiboson system in an overall scalar state was discussed.In this paper the most general scalar state of two boson system was considered.If the bosons originate from the Higgs decay, then one of the components of such a general scalar state corresponds to an anomalous coupling of H with the daughter bosons (we explain this point in the present paper in Sec.II).Such a general scalar state of two vector bosons in the context of quantum correlations for the first time was discussed in [23] while the correlations of relativistic vector bosons in [24].The authors of [19] analyzed entanglement and violation of CGLMP inequality in the system of two Z bosons produced in the decay of a Higgs particle assuming the Standard Model interaction of H with the ZZ pair.In [21] entanglement of W bosons produced in H → W W → lνlν channel was considered.The paper [20] explores the possibility of using quantum state tomography methods to determine a density matrix of massive particles produced in weak decays.In [22] entanglement and Bell inequality violation in boson pairs arising in the Standard Model processes H → W W * , H → ZZ * , pp → W W , pp → W Z, and pp → ZZ is considered.
In this paper we discuss entanglement and violation of the CGLMP inequality for a ZZ system produced in Higgs decay.We consider anomalous (beyond the Standard Model) structure of the vertex describing interaction of a Higgs particle with two daughter bosons but limit ourselves to the case of a scalar Higgs.Anomalous coupling parameters in the amplitude describing the interaction between H and ZZ bosons are strongly constrained by measurements of Higgs properties performed at the LHC [25], we discuss this point in detail in next section, below Eq. ( 9).However, our analysis is based on the most general Lorentz-covariant, CPT conserving amplitude (Eq.( 8)) describing coupling of a (pseudo)scalar particle with two spin-1 particles with different masses.Thus, our considerations can be applied to each of such processes.The Higgs decay can be treated as an exemplary process of this kind.
It is also worth to notice the very recent papers [26,27].The authors of the former paper propose to use quantum tomography techniques to bound anomalous coupling in H → W W and H → ZZ decays while the authors of the latter one use entanglement to probe new physics in diboson production.

II. STATE OF THE TWO-BOSON SYSTEM ARISING FROM THE HIGGS DECAY
Our first goal is to construct a quantum state of two bosons arising in the process Let us denote by M the Higgs mass and by k, m 1 and p, m 2 the four-momenta and invariant masses of the in general off-shell Z bosons produced in the decay (1).In the actual decay of the Higgs particle into a pair of Z bosons typically one of them is nearly on-shell and the other one off-shell.Nevertheless, for the sake of generality, we consider Z bosons with arbitrary invariant masses.We decided to work in this framework as it covers all possible scenarios of a general on-shell scalar decaying into vector bosons.On the other hand, when we average the ZZ state over kinematical configurations we use the probability distribution P c (m 1 , m 2 ) (Eq. ( 32)) giving the probability that H decays into bosons with masses m 1 and m 2 .For the actual process (1) this probability density is peaked at m i = m Z (i = 1 or 2).For further remarks on this point see the paragraph above Eq.( 56).Moreover, we treat of-shell particles like ordinary on-shell particles with reduced masses.Similar approach has been applied in previous quantum-information-related studies [19,22] as well as more phenomenologically-oriented papers like [28,29].We will perform our computations in the center of mass (CM) frame.In this frame we denote energies of Z bosons as ω 1 and ω 2 , consequently Using similar notation as in [18,23,24], a general scalar state of two vector bosons with arbitrary masses can be written as where and |(k, λ); (p, σ)⟩ denotes the two-boson state, one boson with the four-momentum k and spin projection along z axis λ, second one with the four-momentum p and spin projection σ.The basis two-particle states fulfill the following orthogonality condition (for k ̸ = p) [30]: The explicit form of amplitude e µ λ (q) for the fourmomentum q = (q 0 , q) with q 0 2 − q 2 = m 2 reads [24] e(q) = [e µ σ (q)] = and These amplitudes fulfill standard transversality condition e µ σ (q)q µ = 0.
To find an interpretation of the parameter c introduced in Eqs.(2,3) let us notice that the two boson state can be computed using the structure of the vertex describing interaction of the Higgs particle with two daughter vector bosons.Following e.g.[28,29] the amplitude corresponding to the most general Lorentz-invariant, CPT conserving coupling of the (pseudo)scalar particle with two vector bosons can be cast in the following form: where λ, σ are spin projections of the final states, v 1 , v 2 , v 3 are three real coupling constants, and ε αβµν is a completely antisymmetric Levi-Civita tensor.
The Standard Model interaction corresponds to v 1 = 1, v 2 = v 3 = 0. On the other hand, v 3 ̸ = 0 implies that Higgs boson contains a pseudo-scalar component and indicates the possibility of CP violation.For the moment, let us consider the case v 3 = 0, v 1 ̸ = 0, and v 2 free.Comparing (2,3) with (8) and taking into account the transversality condition (7) we can relate the parameter c with the coupling constants Let us mention that the assumption v 1 = 0, apart from being unphysical (we know that the Higgs has v 1 ̸ = 0), together with v 3 = 0 leads to a separable state.Therefore, from now on we limit ourselves to the case v 3 = 0, v 1 ̸ = 0, v 2 free, that is we assume that the Higgs boson is a scalar but we admit a beyond Standard Model coupling v 2 ̸ = 0. We would like to stress here that there exist experimental bounds on anomalous couplings v 2 and v 3 .Strong bound comes from the measurements of Higgs boson particles performed at the LHC by the CMS Collaboration [25].The CMS Collaboration paper uses a different parametrization of the amplitude describing the interaction between H and two daughter bosons, instead of v 1 , v 2 they use parameters a ZZ 1 and a ZZ 2 -see Eq. ( 2) in [25].Comparing (2) in [25] and our Eq.( 8) we obtain the following relation between the parameterizations: (with the same proportionality constant).Thus, follow-are given in Table 7 of [25] To estimate the maximal value of 2(kp)/m 2 Z we use Eq.(A5) and assume that one of the Z bosons is onshell and the invariant mass of the off-shell Z boson is equal to zero.Inserting the measured values for the Higgs mass M = 125.25 GeV and Z mass m Z = 91.19 GeV [31] we finally obtain the following bound for experimentally admissible values of c in the process H → ZZ: Ref. [26] suggests that even stronger bound could be obtained using the tomography of the two-boson density matrix.Therefore, in the actual process H → ZZ the range of the parameter c which is not excluded by experimental data is rather narrow.However, as we mentioned in Introduction, the process H → ZZ can be treated as a model for the most general case of a decay of a (pseudo)scalar particle into two spin-1 particles with different masses.That is why in the following part of the paper we do not restrict values of c to the interval (−c max HZZ , c max HZZ ).The following comment is also in order here: when we consider a decay of a scalar particle into two gauge bosons, gauge invariance requires that v 2 coupling appears in the combination v 2 (k µ p ν − η µν (kp)) which is equivalent to c = −1.Thus, the cases c = 0 and c = −1 are indeed special as it was emphasized in [18].
The normalization of the scalar state defined in Eq. ( 2) is the following: with where, for further convenience, we have introduced the parameter κ.Using formulas (A3-A7) we find that in the CM frame κ depends only on masses M , m 1 , m 2 and the parameter c: where The range of possible values of κ depends on the value of c and is the following: Nevertheless, further theoretical constraints must be taken into account to give the physically allowed range for c.In particular, perturbative unitary (see [32] for a recent review applied to Higgs physics) imposes bounds over the values of the anomalous coupling v 2 .Namely, based on ZZ → ZZ scatterings, numerical bounds have been obtained for the H → ZZ anomalous couplings [33].Comparing our amplitude (Eq.( 8)) with Eq. ( 9) from [33] we get the following relation between our parametrization and parametrization used in [33]: Thus, following ( 9) Now, we use Eq.(A5) assuming that one of the Z bosons is on-shell and the invariant mass of the off-shell Z boson is equal to zero and apply Eqs.(11,12,25) from [33] to obtain which gives Therefore, taking into account that according to Table I in [33]: 97, we notice that c has a pole in the allowed range for a ZZH 2 and hence c ∈ (−∞, ∞).That is, we find that the requirement of perturbative unitarity does not limit accessible values of c in the process H → ZZ.Now, denoting in the CM frame we can write the ZZ scalar state (2) as where |k| are given by Eqs.(A6,A7,A4), respectively, and (V is defined in Eq. ( 6)).
Without loss of generality we take n along z axis, n = (0, 0, 1) and simplify our notation: With such a choice we have and the explicit form of the ZZ state in this case reads It is worth taking note that the form of the above state agrees with the most general form of the ZZ state arising in the H → ZZ process conserving parity-compare Eq. ( 7) from [19].Moreover, when we put c = 0 in Eq. ( 16) the state (31) coincides with the ZZ state discussed in [19] where only the Standard Model vertex has been considered.

A. Averaging over kinematical configurations
In the realistic case when the state is determined on the basis of data obtained from various kinematical configurations, one has to average over these configuration.Thus, in such a case we receive a mixed state which for a given c reads (32) where 31)) and P c (m 1 , m 2 ) is a normalized probability distribution.To determine the form of P c (m 1 , m 2 ), we assume that each of the Z bosons produced in the process (1) decays subsequently into massless fermion-antifermion pair, i.e. we consider the process The normalized differential cross section of the decay where Γ i is the decay matrix of Z i → f + i f − i , Ω i is the solid angle related to the final particle f i [19,34].Now, inserting (32) into (34), integrating with respect to solid angles, using the property that Γ i dΩ i = 4π 3 I, and differentiating with respect to m 1 , m 2 we obtain where by ρ(m 1 , m 2 , c) we have denoted non-normalized density matrix ρ(m 1 , m 2 , c), i.e., Next, Eq. ( 8) from [29] in our notation can be written as [35] 1 where N is a normalization factor independent of m 1 , m 2 and c; while the functions λ and D are defined in Eqs.(A1,A2).Comparing ( 35) with (37) (and taking into account ( 36)) we finally obtain the probability distribution For a given value of c the normalization factor N can be determined numerically.In [19] the probability distribution P c=0 (m 1 , m 2 ) has been obtained with the help of Monte Carlo simulation.The results coincide with those computed from (38).Therefore, using (16,32,38), the density matrix averaged over kinematical configurations can be written as 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 a 0 −d 0 a 0 0 0 0 0 0 0 0 0 0 0 0 0 −d 0 b 0 −d 0 0 0 0 0 0 0 0 0 0 0 0 0 a 0 −d 0 a 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 where for better visibility we have framed the non-zero matrix elements: and we have introduced the following notation Note that sometimes it is relevant for phenomenological purposes to implement cuts on the possible values of the boson masses (for example to remove part of the background of a certain scattering process).This feature is easily implemented theoretically via the integrals defining B(n).For instance, when considering a lower cut in the off-shell mass of the vector boson, m i ≥ m cut i , one just needs to modify the lower bound in the integral over the corresponding mass: (47)

III. BELL INEQUALITIES AND ENTANGLEMENT
Now, we are at a position to discuss the violation of Bell inequalities in a system of two ZZ bosons.Various Bell inequalities have been designed for detecting departures from local realism by quantum mechanical systems [36], the most popular one being the CHSH inequality [37].For a system consisting of two d-dimensional subsystems the optimal Bell inequality was formulated in [38,39] and is known as the CGLMP inequality.For two qubits it reduces to the CHSH inequality.We consider here two spin-1 particles therefore we present the CGLMP inequality for d = 3.We assume that Alice (Bob) can perform two possible measurements A 1 or A 2 (B 1 or B 2 ) on her (his) subsystem, respectively.Each of these measurements can have three outcomes: 0,1,2.Let P (A i = B j + k) denotes the probability that the outcomes A i and B j differ by k modulo 3, i.e., P (A i = B j + k) = l=2 l=0 P (A i = l, B j = l + k mod 3), and let us define the following quantity The CGLMP inequality has the form We assume that Alice can perform measurements on one of the Z bosons, Bob on the second one.In principle, what they can measure are spin projections on given directions.A few remarks are in order here.In [18] we discussed broader the problem of choice of the proper spin observable.We have advocated there the Newton-Wigner spin operator.Under our assumptions, when Alice applies this spin operator to the basis vector |λ, σ⟩, this action can be written as the action of (S λ ′ λ ⊗ I) on |λ ′ ⟩ ⊗ |σ⟩, where S i are standard spin-1 matrices (and analogously for Bob).Therefore, from now on we take Alice (Bob) observables as A ⊗ I (I ⊗ B) and identify |λ⟩ ⊗ |σ⟩ ≡ |λ, σ⟩.
A. Probabilities of spin projection measurements for a particular configuration In general, when discussing the violation of CGLMP inequality by the state ρ ZZ (c) (Eq.( 39)) we do not restrict our attention to spin projection measurements only.However, for a particular configuration, the probabilities P λσ that Alice and Bob receive λ and σ when they measure spin projections along directions a and b in the pure state (27), respectively, can be calculated explicitly.We present them below For the case m 1 = m 2 the above probabilities coincide with the probabilities found in [18].
It is worth noticing that if Alice and Bob are allowed to use only spin projections as observables then, in principle, the violation of the Bell inequality would be suboptimal as we are not covering the whole space of possible observables.

B. Bell inequalities in a general ZZ state
Now, we want to answer the question whether the state (39) violates the CGLMP inequality.In general, for a given state ρ there does not exist a simple way to find optimal observables A 1 , A 2 , B 1 , B 2 , i.e., such observables for which the value of I 3 is maximal in the state ρ.Usually, optimal observables are looked for with the help of a certain optimization procedure.In [19] such a procedure was proposed in the considered there Standard Model coupling case, i.e. for c = 0.In this procedure one modifies the well known optimal choice of observables for the maximally entangled state.For completeness we describe the details of this procedure in B 1.This procedure works very well for the case c = 0 and for c close to that value.However, for higher values of |c| this procedure gives the observables which do not violate the CGLMP inequality.
Thus, we have considered also a different optimization procedure.This procedure is inspired by the proof of Theorem 2 in [40].The details of this approach we described in B 2.
In Fig. 1 we present the maximal value of I 3 as a function of c obtained with the help of both mentioned above optimization strategies.This plot shows that the state ρ ZZ (c) can violate the CGLMP inequality (39) for all values of c.For c ∈ (c Z − , c Z + ), where c Z − = −1.3749,c Z + = 1.6690, the optimization procedure proposed in [19] gives higher violation of CGLMP inequality than the procedure proposed in B 2. For other values of c the situation is opposite.The highest value of I 3 we obtained is equal to 2.9047, it is attained for c Z max = −0.8536.Regarding the values |c| < c max HZZ , larger violation of order [2.5, 2.8] is obtained when the optimization strategy from B 1 is implemented, while violations of order [2.2, 2.3] are attained for the optimization presented in B 2.

C. Entanglement of a general ZZ state
To evaluate entanglement of the state (39) we use the computable entanglement measure called logarithmic negativity [41,42]  curve was obtained with the help of the optimization procedure described in B 1 while the curve with the help of the procedure from B 2. Vertical dotted lines delimit the range (−c max HZZ , c max HZZ ) (with c max HZZ = 0.26-compare Eq. ( 13)) of the parameter c admissible by experimental data for the process H → ZZ. where is the trace norm of a matrix A and T B denotes partial transposition with respect to the second subsystem.The trace norm of a matrix A is equal to the sum of all the singular values of A; when A is hermitian then ||A|| 1 is equal to the sum of absolute values of all eigenvalues of A.
If a state ρ is separable then the logarithmic negativity of ρ is equal to zero.Thus, E N (ρ) > 0 indicates that the state ρ is entangled.In Fig. 2 we have plotted the logarithmic negativity of the state (39) with a Z , b Z , and d Z given in Eqs.(45,46,47).We see that the state is entangled for all values of c, the maximal value of the logarithmic negativity equal to 0.9964 is attained for c = −0.7371.It is worth noticing that the state with the highest entanglement corresponds to c = −0.7371while the state with the highest violation of the CGLMP inequality corresponds to c = −0.8536,i.e., these states are different.This observation is consistent with the general property of CGLMP inequality [43].
Concerning the logarithmic negativity for |c| < c max HZZ , it this is considerably far apart from 0, within a range ∼ [0.85, 0.95].This indicates a high grade of entanglement in any ZZ pair stemming from Higgs decays.

D. Impact of a background
Because the reconstruction of ρ ZZ (c) in a collider experiment is done via quantum tomography methods [19,20,44], the presence of systematic and statistical errors as well as the existence of a small background in  39), E N (ρ ZZ (c)) as a function of c.We have inserted the measured values for the Higgs mass, Z mass and Z decay width, i.e., we put a, b, d given in Eqs.(45,46,47).The cusp occurs for the value of c for which the parameter d Z in the density matrix vanishes and is caused by the change in monotonicity of the trace norm of the partially transposed matrix.Vertical dotted lines delimit the range (−c max HZZ , c max HZZ ) (with c max HZZ = 0.26-compare Eq. ( 13)) of the parameter c admissible by experimental data for the process 2 processes lead to a modification of its exact form.Following the discussion given in [19], we will focus on analyzing H → ZZ → e + e − µ + µ − , since it constitutes one of the cleanest channels to be explored at the LHC.In principle, the two Z bosons are cleanly identified.Due to the fact that one of them is nearly onshell, it gives two leptons whose invariant mass is close to m Z , while the remaining two leptons have a much lower invariant mass.We have labeled the (very close to real) Z boson with largest invariant mass as Z 1 and its four-momentum could be reconstructed from its decay products l + 1 l − 1 .On the other hand, the off-shell Z boson is labeled as Z 2 and its momentum is determined summing up the momenta of its decay products l + 2 l − 2 .Concerning both the systematic and statistical errors of the tomography procedure, it was stated in [19] that the statistical one dominates with respect to the systematics and detector resolution.An estimation of the former was computed in this previous paper and the results show that, even with these errors, a violation of Bell inequalities for the ρ ZZ (0) state could be probed at the 4.5σ level in the HL-LHC.Regarding the background of the process, the main one comes from the electroweak one pp → ZZ/Zγ → 4l, being this one about 4 times smaller at the Higgs peak [45].Nevertheless, as claimed in [19], a background subtraction will be necessary before computing the entanglement observable and evaluating the CGLMP inequality.In general, the nonnegligible background will slightly contribute to the sta-tistical uncertainty of the measurements.Moreover, as proposed in [19], the larger the invariant mass m 2 of the off-shell Z boson, the more entangled the ρ ZZ (c) state.Therefore, requiring a lower cut on m 2 leads to an interplay between increasing the entanglement (hence the violation of the CGLMP inequality) and decreasing the statistics (thus increasing the uncertainty in the measurements).In Figs. 3, 4, 5 we have studied the theoretical dependence of I 3 and E N on c, once the cuts m 2 ≥ 0, 10, 20, 30 GeV are implemented.
Returning for a moment to our discussion of "offshellness" of Z bosons produced in the decay (1), let us notice that I  3 and  4) are identical as I (1) 3 and I (2) 3 plotted in Fig. 1 where we allowed arbitrary masses m 1 and m 2 (of course constrained by the four-momentum conservation).The same holds for E N plotted in Figs. 5 (blue line) and 2. Of course this coincidence is not accidental, it results from the fact that the probability distribution Finally, in order to estimate the overall allowed uncertainty in both the entanglement and the Bell inequality violation, we have estimated the noise resistance of the CGLMP inequality violation with respect to the white noise.To this aim we considered the state ρ ZZ (c) mixed with the identity operator, i.e. the state Now, the noise resistance we define as a minimal value of λ, λ min , for which the state (56) violates the CGLMP inequality.Inserting the state (56) into the CGLMP inequality (B1) and taking into account that Tr(O Bell ) = 0 we obtain We obtained the maximal value of Tr(ρ ZZ (c)O Bell ) = I 3 with the help of two different optimization procedures ( B 1 and B 2) and denoted as I 3 , respectively (compare Fig. 1).Therefore, we have and this value we have plotted in Fig. 6.The plots show that for values of c close to 0 (which are the expected ones, due to the present bounds in anomalous couplings for the HZZ vertex [46]), one can stand up to almost a 20% of noise and still attain a violation of the CGLMP inequality and hence an entangled state.(45,46,47).We have applied the optimization procedure described in B 1 and assumed the cuts m 2 ≥ 0, 10, 20, 30 GeV are implemented.Vertical dotted lines delimit the range (−c max HZZ , c max HZZ ) (with c max HZZ = 0.26-compare Eq. ( 13)) of the parameter c admissible by experimental data for the process H → ZZ.  (45,46,47).We have applied the optimization procedure described in B 2 and assumed the cuts m 2 ≥ 0, 10, 20, 30 GeV are implemented.Vertical dotted lines delimit the range (−c max HZZ , c max HZZ ) (with c max HZZ = 0.26-compare Eq. ( 13)) of the parameter c admissible by experimental data for the process H → ZZ.

IV. CONCLUSIONS
In conclusions, we have analyzed entanglement and Bell inequality violation in a system of two Z bosons produced in Higgs decay.We consider beyond the Standard Model structure of the vertex describing interaction of H with daughter bosons.The amplitude corresponding to the most general Lorentz-invariant, CPT conserving coupling of a (pseudo)scalar particle with two vector bosons depends on three coupling constants v 1 , v 2 , v 3 , and is explicitly given in Eq. ( 8).The Standard Model interaction corresponds to v 1 = 1, v 2 = v 3 = 0 while   (45,46,47).We have assumed the cuts m 2 ≥ 0, 10, 20, 30 GeV are implemented.Vertical dotted lines delimit the range (−c max HZZ , c max HZZ ) (with c max HZZ = 0.26-compare Eq. ( 13)) of the parameter c admissible by experimental data for the process H → ZZ. v 3 ̸ = 0 implies Higgs boson with a pseudo-scalar component and indicates the possibility of CP violation.In this paper we have considered the case v 3 = 0, v 1 ̸ = 0, v 2 free, i.e. we have assumed the scalar Higgs boson but admitted a beyond Standard Model coupling v 2 ̸ = 0.In such a case, the state of produced bosons, beyond four-momenta and spins, can be characterized by a single parameter c which, up to normalization is equal to v 2 /v 1 (compare Eq. ( 9)).Under such assumptions, in the center-of-mass frame, we have determined the most general state of ZZ boson pair for a particular event H → ZZ.Next, we have considered a more realistic case when data are collected from different kinematical configurations.In such a sit-uation a ZZ state can be calculated by averaging over those configurations with respect to a proper probability distribution function (PDF).Thus, assuming further that each Z boson decays into fermion-antifermion pair, we have derived the corresponding PDF and computed the ZZ boson density matrix.Finally, we have shown that this matrix is entangled and violates the CGLMP inequality for all values of coupling (i.e. for all values of c) including the range admissible by experimental data [25].The procedure to check this is completely general and can be applied for any other decay of a scalar particle into vector bosons, with their corresponding decay to fermions, once the PDF of the latter decay is known.Moreover, our preliminary studies show that the inclusion of a CP-odd anomalous coupling should not qualitatively change the results derived.However, in this case the optimization strategy is more involved and work is still in progress.
Summarizing, this work settles a constructive way of probing the entanglement and violation of Bell inequalities of any vector boson pair coming from a spin-0 particle, independently of the value of the couplings in hand (as long as the interactions among particles are CPT and Lorentz invariant).This feature is of a complete novelty in the literature and states the highly non-trivial fact that non locality (and hence entanglement) of vector bosons in these kinds of processes is theoretical ensured in any phenomenological model (it could have been the case in which although the state is entangled, it does not violate any Bell inequality).Thus, the only limitation in checking the quantum behavior of these processes comes from the experimental side.In particular, this work gives the theoretical framework to test the quantum nature of processes in a great variety of phenomenological models, covering for instance models with extended scalar sectors as well as axion like particles (ALP) models where we let the ALP to interact with vector bosons.permutation matrices: where O is the 3 × 3 null matrix and C is the 3 × 3 cyclic permutation matrix Each U from (B2) can be taken as an element of SU (3) group, this group has 8 parameters.Thus, to perform the full optimization of O Bell for a given state one should check the 8 4 dimensional parameter space.

Strategy 1
To simplify this task, in [19] the following approach was applied.It is known what is the form of the optimal Bell operator for the maximally entangled state where For the maximally entangled state Ω M E = 1 √ 3 I, thus this state is invariant on the action U ⊗ U * .Therefore, the optimal Bell observable for the state |ψ scalar ZZ ⟩| κ=1 , instead of the form (B7) can be written in an equivalent form where U is an arbitrary unitary matrix.The value of I 3 with the Bell operator given in (B10) in the state |ψ scalar ZZ ⟩| κ=1 is the same for all unitary matrices U .
For κ ̸ = 1 the Bell operator (B7) is no longer an optimal one.Moreover, different choices of U in (B10) lead to different values of I 3 in this case.Thus, one can look for an optimal choice taking the Bell operator in the form (B10) and optimizing over all U matrices.This can be simplified further by observing that the state (O A ⊗ I)|ψ scalar ZZ ⟩(m 1 , m 2 , c), according to (30,B9) is represented by the matrix To obtain Fig. 1, for each value of c we performed the optimization using the above parametrization.The optimal choice of R, θ, ϕ 1 , and ϕ 2 depends on the value of c, this change of parametrization is responsible for the cusp appearing in the plot.

Strategy 2
In this case we define a matrix and assume that observables used by Alice and Bob are represented by the following unitary matrices: (B17) Next, we are looking for such value of t which gives the highest violation of the CGLMP inequality.It appears that the optimal value is t = − π 4 .We used this value to plot Fig. 1.

FIG. 1 :
FIG. 1: In this figure we present the maximal value of I 3 in the state (39) as a function of c.We have inserted the measured values for the Higgs mass, Z mass and Z decay width, i.e., we put a, b, d given in Eqs.(45,46,47).The I (1) 3

FIG. 2 :
FIG. 2:In this figure we present logarithmic negativity of the state (39), E N (ρ ZZ (c)) as a function of c.We have inserted the measured values for the Higgs mass, Z mass and Z decay width, i.e., we put a, b, d given in Eqs.(45,46,47).The cusp occurs for the value of c for which the parameter d Z in the density matrix vanishes and is caused by the change in monotonicity of the trace norm of the partially transposed matrix.Vertical dotted lines delimit the range (−c max HZZ , c max HZZ ) (with c max HZZ = 0.26-compare Eq. (13)) of the parameter c admissible by experimental data for the process H → ZZ.
Actually, the resistance to noise increases with the invariant mass of the off-shell Z boson, reaching for the cut m 2 ≥ 30 GeV a resistance of 20% for c ∈ [−3, 3] and a resistance of almost 30% for c ∈ [−2, 2].

m 2
FIG.3:In this figure we present the maximal value of I 3 in the state (39) as a function of c.We have inserted the measured values for the Higgs mass, Z mass and Z decay width, i.e., we put a, b, d given in Eqs.(45,46,47).We have applied the optimization procedure described in B 1 and assumed the cuts m 2 ≥ 0, 10, 20, 30 GeV are implemented.Vertical dotted lines delimit the range (−c max HZZ , c max HZZ ) (with c max HZZ = 0.26-compare Eq. (13)) of the parameter c admissible by experimental data for the process H → ZZ.

m 2
FIG.4:In this figure we present the maximal value of I 3 in the state (39) as a function of c.We have inserted the measured values for the Higgs mass, Z mass and Z decay width, i.e., we put a, b, d given in Eqs.(45,46,47).We have applied the optimization procedure described in B 2 and assumed the cuts m 2 ≥ 0, 10, 20, 30 GeV are implemented.Vertical dotted lines delimit the range (−c max HZZ , c max HZZ ) (with c max HZZ = 0.26-compare Eq. (13)) of the parameter c admissible by experimental data for the process H → ZZ.
FIG. 5: In this figure we present logarithmic negativity of the state (39), E N (ρ ZZ (c)) as a function of c.We have inserted the measured values for the Higgs mass, Z mass and Z decay width, i.e., we put a, b, d given in Eqs.(45,46,47).We have assumed the cuts m 2 ≥ 0, 10, 20, 30 GeV are implemented.Vertical dotted lines delimit the range (−c max HZZ , c max HZZ ) (with c max HZZ = 0.26-compare Eq. (13)) of the parameter c admissible by experimental data for the process H → ZZ.
FIG.6:In this figure we present λ min (58), as a function of c.We have inserted the measured values for the Higgs mass, Z mass and Z decay width, i.e., we put b, d given in Eqs.(45,46,47).We have assumed the cuts m 2 ≥ 0, 10, 20, 30 GeV are implemented.Vertical dotted lines delimit the range (−c max HZZ , c max HZZ ) (with c max HZZ = 0.26-compare Eq. (13)) of the parameter c admissible by experimental data for the process H → ZZ.