# Exclusion Limits on a Scalar Decaying to Photons and Distinguishing Its Production Mechanisms.

1. IntroductionThere are numerous theoretical motivations to expect that the Standard Model (SM) is not the complete story and the scalar spectrum of a larger theory may be richer than to possess only one neutral scalar, the Higgs boson. From this expectation searches for new scalars are continuously being carried out at the Large Hadron Collider (LHC) in various channels. Along these directions, no confirmed hint has been found so far in any of these searches. Nevertheless, some anomalies in some of these searches have drawn significant attention in the high energy physics community, recently. Among them, the most famous one is the 750 GeV diphoton excess [1, 2] which created a lot of excitements in the community. Before the excess went away with more data, numerous attempts have been made to explain the excess (see [3] for a review and a long list of references). Another important excess was the diboson excess around 2 TeV resonance mass [4-6]which also later turned out to be a statistical fluctuation.

Searches for heavy scalars at the LHC are generally being carried out in the diphoton, diboson, or dijet resonance searches. The diphoton channel, among them, is particularly important as this channel provides a comparatively cleaner background. Higgs boson was first discovered in the diphoton channel at the LHC [7, 8]. TeV-scale scalars decaying into a diphoton system is one of the key predictions of many beyond the Standard Model (BSM) theories. Various possibilities have been extensively explored in the context of the 750 GeV diphoton excess (see the reference list of [3] for various models that predict TeV-scale diphoton resonances). To test these predictions, LHC run-II provides us with a great opportunity to observe a diphoton resonance of mass up to a few TeV. In this paper we particularly focus on the diphoton final state for these reasons.

If a particle decays to diphotons, it must either be spin-0 or spin-2 in nature as the spin-1 particles decay to on-shell diphotons is forbidden by the Landau-Yang theorem [9, 10]. A spin-2 particle or graviton couples universally to all matter fields through energy-momentum tensor. Various extra dimensional models like the ADD model [11] or the RS model [12] predict the existence of graviton. If a resonance in the diphoton system mediated by graviton is observed, one would expect resonances at the same mass in other possible channels also. Therefore, simultaneous studies in various channels might be more illuminating for the spin-2 particle. The current limits on the graviton mass are already quite high, around ~2-3 TeV [13, 14]. On the other hand, scalars of mass ~1 TeV decaying into diphotons, which is a typical signature of many models, are still allowed by the LHC data. These scalars can be produced at the LHC in various ways, namely. through the gg, qq, [gamma][gamma], [gamma]Z, WW, or ZZ fusions. In this paper we consider a model-independent effective field theory (EFT) approach where the scalar can be produced and decayed (two-body) in different possible ways as mentioned. But we only concentrate on the diphoton decay mode in this paper as stated earlier.

First, we derive the available parameter space for a scalar (produced in different ways) decays to diphotons using our EFT approach. These limits will be grossly model independent and can be used to set limits on other models wherever applicable. If a scalar resonance will actually be seen in future, the most obvious question that will arise is how the scalar is produced? A most common way to decipher the production mechanism of a heavy scalar is to look at various kinematic distributions especially various jet observables, which are important in this regard. This has been investigated to some extent in the literature in the context of the 750 GeV diphoton excess [15-19]. In this paper we revisit some of the jet observables and show their effectiveness in distinguishing different production modes. We, then, use a multivariate analysis (MVA) by combining many kinematic variables to distinguish different production modes more efficiently.

This paper is organized as follows: in Section 2, we employ an effective Lagrangian for the scalar [PHI]; in Section 3, we discuss about the decays and various production modes of [PHI] at the LHC and derive exclusion limits on the mass and couplings from the latest diphoton resonance search data. In the same section, we discuss how two different production modes of [PHI] can be distinguished using a MVA analysis. Finally, we conclude in Section 4.

2. Effective Lagrangian

We consider an EFT approach where a heavy scalar [PHI] interacts with the SM gauge bosons through the dimension-5 operators and with the SM quarks through the dimension-4 operators. Assuming [PHI] is a CP-even real scalar, we employ the following effective Lagrangian:

[mathematical expression not reproducible], (1)

where the field-strength tensors corresponding to gluon (g), photon ([gamma]), [W.sup.[+ or -]], and Z bosons are [G.sup.a.sub.[mu][nu]], [A.sub.[mu][nu]], [([W.sup.[+ or -]]).sub.[mu][nu], and [Z.sub.[mu][nu]], respectively, and their generic form is [V.sub.[mu][nu]] = [[partial derivative].sub.[mu]][V.sub.[nu]] - [[partial derivative].sub.[nu]][V.sub.[mu]] where V = {g, [gamma], [W.sup.[+ or -]], Z}. All the dimension-5 operators are suppressed by the new physics scale [LAMBDA]. In general, [LAMBDA] could be different for different operators, but we assume they are the same for all the operators. Note that only the [PHI][bar.q]q is a dimension-4 operator and we introduce the electroweak symmetry breaking scale v [approximately equal to] 246 GeV in association with [[kappa].sub.qq] to bring the scale [LAMBDA] in the interaction. The motivation behind this is that if [PHI][bar.q]q operators are effectively originated from new physics then effective couplings are expected to contain the imprint of the scale [LAMBDA] (possibly in the form v/[LAMBDA] or some power of this ratio). This way of parameterizing the [PHI][bar.q]q couplings also enables us to present exclusion limits on all the couplings in the [[kappa].sub.xy]/[LAMBDA] form. Here, we use the notation [[kappa].sub.xy] to denote a generic dimensionless coupling associated with the [PHI]xy vertex. The scalar [PHI] can, in general, couple differently with the different SM quarks. For simplicity, in this analysis, we assume a single coupling [[kappa].sub.qq], the same for all the SM quarks. Note that, in all interactions with the gauge bosons, the normalization factor is so chosen such that the corresponding Feynman rule takes the form

(i[[kappa].sub.xy]/[LAMBDA]) ([g.sup.[alpha][beta]][p.sub.1] x [p.sub.2] - [p.sup.[beta].sub.1][p.sup.[alpha].sub.2]), (2)

where [p.sub.1] and [p.sub.2] are the 4 momenta of two gauge bosons [V.sup.[alpha].sub.1] and [V.sup.[beta].sub.2], respectively, directed towards the vertex. The Feynman rule for the [PHI][bar.q]q interaction is i[[kappa].sub.qq]v/[LAMBDA].

In general, the new scalar [PHI] can mix with the 125 GeV scalar ([h.sub.125]) with a mixing angle [alpha]. This leads to the scaling of all the couplings of [h.sub.125] by a factor cos [alpha]. Although this would not change the branching ratios (BRs) of [h.sub.125], it would change the production cross section of [h.sub.125] by a factor [cos.sup.2][alpha]. Since all the measured signal strengths are pretty close to unity, this will make cos [alpha] close to one. That is why, in this paper, we have neglected any mixing between [PHI] and [h.sub.125] for simplicity.

3. Phenomenology

In addition to the SM Lagrangian, we implement the effective Lagrangian of [PHI] shown in (1) in FeynRules [20] to generate the Universal FeynRules Output [21] model files for the MadGraph [22] event generator. We use the MMHT14LO [23] parton distribution functions (PDFs) for event generation. This PDF set includes the photon PDF which has been computed following the approach described in [16, 24]. We use the factorization scale [[mu].sub.F] and the renormalization scale [[mu].sub.R] at [M.sub.[PHI]] in our analysis. Generated events are further showered and hadronized including multiple parton interactions by using Pythia8 [25]. We perform detector simulation using Delphes [26] which uses FastJet [27] for jet clustering. Jets are clustered using the anti-[k.sub.T] algorithm [28] with R = 0.4. We analyze the reconstructed objects by implementing ATLAS selection cuts [29], which we summarize in Section 3.3. For MVA, we use the adaptive Boosted Decision Tree (BDT) algorithm in the TMVA [30] framework.

3.1. Decays of [PHI]. From the Lagrangian in (1), we have the following two-body decay modes of [PHI], namely, [PHI] [right arrow] xy, where xy = {gg, qq, [gamma][gamma], [gamma]Z, WW, ZZ}. The partial widths for these decay modes are given by the following expressions:

[mathematical expression not reproducible], (3)

where V denotes the electroweak gauge bosons [W.sup.[+ or -]] and Z. There could be subdominant three-body decays of [PHI] possibly mediated through an off-shell gauge boson. If the intermediate gauge boson is massless, in case of gluons or photons, the three-body BRs are nonnegligible especially when [M.sub.[PHI]] is large [31]. In this analysis, we consider the two-body and three-body decays of [PHI] to obtain the total width where the three-body decay widths are computed numerically using MadGraph. Partial widths of three-body decay modes where an off-shell gauge boson goes to [W.sup.+][W.sup.-] pair grow very rapidly with increasing scalar mass. This is due to the contribution coming from the longitudinal polarization of W bosons. Therefore, in high mass region, BR for [PHI] [right arrow] [gamma][gamma] reduces substantially.

3.2. Production of [PHI] at the LHC. When all [[kappa].sub.xy] in (1) are nonzero, the scalar can be produced from the gg, qq, [gamma][gamma], [gamma]Z, WW, and ZZ fusions at the LHC. In Figure 1, we show the partonic cross sections of different production modes of [PHI] at the 13 TeV LHC for [[kappa].sub.xy] = 1 (taking one at a time) and [LAMBDA] = 1 TeV. In case of the production of [PHI] through the WW or ZZ fusions, initial W and Z come from the quark splitting. Therefore, [PHI] is produced in association with at least two jets for this case. Similarly, for the [gamma]Z initiated production, [PHI] is produced in association with at least one jet. Partonic cross sections are computed by applying the following generation level cuts on the jets (j) and photons ([gamma]) wherever applicable

[p.sub.T] (j), [p.sub.T] ([gamma]) > 20 GeV; [absolute value of ([eta] (j))] < 5, [absolute value of ([eta] ([gamma]))] < 2.5; [DELTA]R (jj), [DELTA]R ([gamma][gamma]), [DELTA]R (j[gamma]) > 0.4. (4)

Here, transverse momentum, pseudorapidity, and separation in the [eta]-[phi] plane are denoted by [p.sub.T], [eta] and [DELTA]R, respectively. These basic cuts are used to avoid any soft divergence present at the event generation level and stricter selection cuts are applied at the level of reconstructed event analysis after detector simulation. Note that all cross sections scale as [([[kappa].sub.xy]/[LAMBDA]).sup.2], and therefore, we present them by choosing [k.sub.xy] = 1 and [LAMBDA] = 1 TeV such that one can translate it easily for other values.

For all the six types of production of [PHI], we generate parton level events with up to two jets in the final state. These events are passed to Pythia8 [25] for showering and hadronization. This process may introduce double counting between the matrix element partons and the parton showers. To generate inclusive signal events by avoiding any double counting, we use the MLM matching [32] technique to match the matrix element partons with the parton shower. Inclusive signal events including up to two jets for the gg, qq, and [gamma][gamma] fusions are generated by combining the following processes:

[mathematical expression not reproducible], (5)

where we set the matching scale [Q.sub.cut] ~ 125 GeV. The curved connections above two photons signify that they come from the decay of [PHI]. To determine the appropriate [Q.sub.cut] for these production processes, we have done three important checks, namely, smooth transition in the differential jet-rate distributions between events with N and N + 1 jets; matched cross sections are within ~10% of the zero jet contribution and also do not vary much with the [Q.sub.cut] variation once we have chosen it properly. For the [gamma]Z, WW, or ZZ fusion productions, the initial W or Z come from the quark splitting and we have additional jets at the Born level process. Therefore, the WW and ZZ fusion events are generated only at the pp [right arrow] [PHI]jj level and no matching is required for these cases. But for the [gamma]Z fusion, we do use matching by combining the processes [gamma]p [right arrow] [PHI]j and pp [right arrow] [PHI]jj with [Q.sub.cut] ~ 30 GeV. The dominant SM [gamma][gamma] background (about 90% of the total) comes from the q[bar.q] [right arrow] [gamma][gamma] process. We generate this background by matching up to 2 jets with [Q.sub.cut] ~ 20 GeV.

3.3. Exclusion from the LHC Data. Diphoton resonance searches at the LHC using run-I and run-II data set strong upper limits (ULs) on [sigma] x BR of a spin-0 or spin-2 resonances [29, 33]. It should be noted that these searches are generally optimized for an s-channel resonance production through the gg fusion followed by its decay to two photons. If the resonance is not produced from the gg fusion, the selection cut efficiencies can vary depending on the different production mechanisms of the resonance. For a particular production mechanism, it can also vary significantly on the number of selected photons and jets. Therefore, in order to derive exclusion limits on the model parameters by recasting the limits on [sigma] x BR from an experiment, one has to properly take care of the selection cut efficiencies. This can be done properly by using the following relation [31]:

[N.sub.s] = [([sigma] x BR).sub.s] x [[epsilon].sub.s] x L = [summation over (i)] [([sigma] x BR).sub.i] x [[epsilon].sub.i] x L, (6)

where [N.sub.s] is the UL on the number of signal events, which can be written as the product of the signal cross section [([sigma] x BR).sub.s] (produced through a particular mechanism used in the analysis), the corresponding signal cut efficiency [[epsilon].sub.s], and the luminosity L. When different types of production mechanisms contribute to any experiment, [N.sub.s] can be expressed by the sum [[summation].sub.i][([sigma] x BR).sub.i] x [[epsilon].sub.i] x L. Here, i runs over all the contributing production mechanisms.

To see the change in efficiency for the different production mechanisms and also for the different resonance masses, we roughly employ the following event selection cuts used by the ATLAS collaboration for their spin-0 diphoton resonance search as listed below [29].

(1) Transverse energy of the two selected photons satisfies [E.sub.T]([[gamma].sub.1]) > 40 GeV and [E.sub.T]([[gamma].sub.2]) > 30 GeV and transverse momenta of selected jets satisfy [p.sub.T](j) > 25 GeV for [absolute value of ([eta](j))] < 2.5 and [p.sub.T](j) > 50 GeV for [absolute value of ([eta](j))] > 2.5. Here, [[gamma].sub.1] and [[gamma].sub.2] denote the highest-[p.sub.T] and second highest-[p.sub.T] photons, respectively.

(2) Pseudorapidity of the selected photons satisfies [absolute value of ([eta]([gamma]))] < 2.37 excluding the barrel-endcap region 1.37 < [absolute value of ([eta]([gamma]))] < 1.52 and jets [absolute value of ([eta](j))] < 4.4.

(3) Separation in the [eta]-[phi] plane between the two selected photons or any photon-jet or jet-jet pair satisfies [DELTA]R([[gamma].sub.1], [[gamma].sub.2]), [DELTA]R([[gamma].sub.1], j), [DELTA]R([[gamma].sub.2], j), [DELTA]R(j, j) > 0.4.

(4) Invariant mass of the two selected photons M([[gamma].sub.1], [[gamma].sub.2]) satisfies [E.sub.T]([[gamma].sub.1]) > 0.4M([[gamma].sub.1], [[gamma].sub.2]) and [E.sub.T]([[gamma].sub.2]) > 0.3M ([[gamma].sub.1], [[gamma].sub.2]).

In addition to the above set of cuts, we also apply default photon isolation cuts given in Delphes for the ATLAS detector. Jets with high-[eta]mainly come from the vector boson fusion topologies. We use a threshold [p.sub.T](j) > 50 GeV for [absolute value of ([eta](j))] > 2.5 for better sensitivity [34]. In Figure 2, we show cut efficiencies for the cuts listed above for different production modes as functions of [M.sub.[phi]]. The cut efficiency for ATLAS for their spin-0 resonance produced through the gg fusion is roughly about 62% [29] and we find very close agreement (around 60%) using our analysis codes. After validating our codes, we compute cut efficiencies for the other production modes for the selection cuts mentioned above and find that they do not vary much, only up to ~15% for different production modes. It is pointed out in the ATLAS paper [29] that the cut efficiencies for different production modes would not differ much for their signal criteria (fiducial region). As expected, in the high mass region, [M.sub.[PHI]] [greater than or equal to] 1 TeV, cut efficiencies become insensitive to the mass.

In our EFT approach, there are six free couplings [[kappa].sub.xy] that affect the production of [PHI]. But taking all [[kappa].sub.xy] nonzero at the same time will make the analysis very complicated. Therefore, for simplicity, we choose only one [[kappa].sub.xy] as nonzero at a time, in addition to nonzero [[kappa].sub.[gamma][gamma]], and show the two-dimensional (2D) exclusion regions (colored) in the [[kappa].sub.xy]/[LAMBDA]-[[kappa].sub.[gamma][gamma]]/[LAMBDA] plane for four benchmark masses, [M.sub.[PHI]] = 0.75, 1, 1.5, 2 TeV (presented in Figure 3). Only in Figure 3(a), we show the exclusion regions (colored) in the [M.sub.[PHI]]-[[kappa].sub.[gamma][gamma]]/[LAMBDA] plane assuming all [[kappa].sub.xy] are zero except [[kappa].sub.[gamma][gamma]]. To derive these limits, we recast 95% confidence level (CL)UL on the [sigma] x BR for the [gamma][gamma] spin-0 resonance search by the ATLAS collaboration at the 13 TeV with L = 36.7 [fb.sup.-1] [29]. This analysis is done for the resonance width [GAMMA] = 4 MeV. If the width of a particle is very small compared to its mass, one can safely use the narrow width approximation (NWA). In all our results we use the NWA ignoring any interference effect between the signal and the background.

The bumpy nature in the exclusion limit on [[kappa].sub.[gamma][gamma]]/[LAMBDA] in Figure 3(a) is due to nonsmooth nature of the observed UL on [sigma] x BR. The highest value of [[kappa].sub.[gamma][gamma]]/[LAMBDA] that is excluded is ~0.05 around [M.sub.[PHI]] ~ 1 TeV. By choosing a value for [[kappa].sub.[gamma][gamma]], one can translate this limit to [LAMBDA]. For instance, choosing [[kappa].sub.[gamma][gamma]] = 1 for [M.sub.[PHI]] ~ 1 TeV, one finds that [LAMBDA] [less than or equal to] 20 TeV is excluded. Basic shape of the exclusion regions in Figures 3(b) and 3(c) is different from the ones in Figures 3(d), 3(e), and 3(f). This is because the gg and the qq fusion productions dominate over the [gamma][gamma] fusion mode for [[kappa].sub.xy] = 1 as seen in Figure 1. On the other hand, cross sections for the WW, ZZ, and [gamma]Z fusion modes are smaller than the [gamma][gamma] mode for unity [[kappa].sub.xy]. One should also notice that exclusion regions in Figures 3(d) and 3(f) are slightly asymmetric around [[kappa].sub.[gamma][gamma]] = 0 axis. This is due to the interference effect between the [gamma][gamma] and the [gamma]Z or ZZ production modes. On the other hand, no interference is possible between the [gamma][gamma] and the gg, qq, or WW fusion modes. Note that exclusion limits become insensitive to [[kappa].sub.xy]/[LAMBDA] as we go to higher values. This is because the production cross section [sigma] varies as [[kappa].sup.2.sub.xy] and BR to diphoton B[R.sub.[gamma][gamma]] varies as ~[[kappa].sup.2.sub.[gamma][gamma]]/[[kappa].sup.2.sub.xy] for [[kappa].sub.xy] [much greater than] [[kappa].sub.[gamma][gamma]] region. This makes [sigma] x BR ~ [[kappa].sup.2.sub.[gamma][gamma]] for large [[kappa].sub.xy].

To include higher order effects, we use a constant next-to-leading order (NLO) K-factor of 2 for the gg fusion [35]. The NLO corrections to a heavy scalar produced from the b[bar.b] fusion are computed in [36] where it is found that the NLO K-factor is close to 1 for heavier masses. If the scalar is produced from the light quark fusions, one might expect slightly bigger K-factor. Here, we assume it to be 1 since it is not available in the literature. For the [gamma][gamma], [gamma]Z, WW, and ZZ we assume it to be 1.3 [37]. The actual values of the K-factors for different channels can be slightly different from the constant values we have used but they have very little effect on the exclusion limits.

In general, when [[kappa].sub.[gamma][gamma]] and any one [[kappa].sub.xy] are nonzero, the production cross section can be expressed as

[mathematical expression not reproducible], (7)

where [[sigma].sup.0.sub.I] is the interference contribution and [[sigma].sup.0.sub.xy] ([[sigma].sup.0.sub.[gamma][gamma]]) in the r.h.s. is the production cross sections through the xy ([gamma][gamma]) fusion (see Figure 1). These cross sections as functions of mass have been computed numerically by interpolating cross sections points in the mass range [M.sub.[PHI]] = 0.5-2.5 TeV. When more than one [[kappa].sub.xy] are nonzero, the combined selection cut efficiency, in general, depends on [M.sub.[PHI]] and [[kappa].sub.xy]. Including cut efficiencies in (7) (omitting the functional dependence on [M.sub.[PHI]]), we get

[sigma] x [epsilon] = [[LAMBDA].sup.-2] {[[kappa].sup.2.sub.xy][[sigma].sup.0.sub.xy] [[epsilon].sub.xy] + [[kappa].sup.2.sub.[gamma][gamma]][[sigma].sup.0.sub.[gamma][gamma]] [[epsilon].sub.[gamma][gamma]] + [[kappa].sub.xy][[kappa].sub.[gamma][gamma]] [[sigma].sup.0.sub.I] [[epsilon].sub.I]}, (8)

where [[epsilon].sub.xy], [[epsilon].sub.[gamma][gamma]] are the cut efficiencies for the pure xy and pure [gamma][gamma] fusion production modes, respectively, and they are functions of [M.sub.[PHI]] only, whereas the combined efficiency [epsilon] and [[epsilon].sub.I] associated with the interference term are functions of [M.sub.[PHI]], [[kappa].sub.xy], and [[kappa].sub.[gamma][gamma]]. We have seen that [[epsilon].sub.I] is mostly sensitive to [M.sub.[PHI]] but not to the couplings. Therefore, for simplicity we use [[epsilon].sub.I] = [[epsilon].sub.xy] ([M.sub.[PHI]]) for [[kappa].sup.2.sub.xy][[sigma].sup.0.sub.xy] > [[kappa].sup.2.sub.[gamma][gamma]][[sigma].sup.0.sub.[gamma][gamma]] region and [[epsilon].sub.I] = [[epsilon].sub.[gamma][gamma]]([M.sub.[PHI]]) for [[kappa].sup.2.sub.[gamma][gamma]][[sigma].sup.0.sub.[gamma][gamma]] > [[kappa].sup.2.sub.xy][[sigma].sup.0.sub.xy] region. Branching fraction in the [gamma][gamma] channel can be expressed as

B[R.sub.[gamma][gamma]] ([M.sub.[PHI]], [[kappa].sub.xy], [[kappa].sub.[gamma][gamma]]) = [[[kappa].sup.2.sub.[gamma][gamma]][[GAMMA].sub.[gamma][gamma]] ([M.sub.[PHI]])]/[[[kappa].sup.2.sub.[gamma][gamma]][[GAMMA].sub.[gamma][gamma]] ([M.sub.[PHI]]) + [[kappa].sup.2.sub.xy][[GAMMA].sub.xy] ([M.sub.[PHI]])], (9)

where [GAMMA]'s are known analytically from (3). Finally, we derive exclusion regions in Figure 3 by using (8) and (9) in (6).

In Figure 4, we show combined cut efficiency for the coupling assumption [[kappa].sub.gg], [[kappa].sub.[gamma][gamma]] [not equal to] 0 (all other [[kappa].sub.xy] are zero). This combined cut efficiency should lie between the two individual efficiencies [[epsilon].sub.[gamma][gamma]] ~ 55% and [[epsilon].sub.gg] ~ 60% according to its definition. One might also be interested to see the behavior of [sigma], [sigma] x BR and [sigma] x BR x [epsilon] for different coupling assumptions. In Figure 5, we show these three quantities in the 2D plane for the two cases: [[kappa].sub.gg], [[kappa].sub.[gamma][gamma]] [not equal to] 0 and [[kappa].sub.WW], [[kappa].sub.[gamma][gamma]] [not equal to] 0 (other [[kappa].sub.xy] are set to zero).

3.4. Distinguishing Different Production Modes. A common way to distinguish different production modes of a heavy scalar is to scrutinize various kinematic distributions especially the jet activities associated with the scalar. It was pointed out in [15-19] that the jet multiplicity ([N.sub.jet]) distribution could be very important in this regard. In Figure 6, we show the normalized [N.sub.jet] distributions for various production modes of the scalar and compare them with the SM prediction. These distributions are obtained assuming [M.sub.[PHI]] = 1 TeV at the 13 TeV LHC with 50 [fb.sup.-1] integrated luminosity with the diphoton invariant mass (M[gamma][gamma]) satisfying [absolute value of (M[gamma][gamma] - [M.sub.[PHI]])] < 150 GeV, in addition to the set of cuts defined earlier. Our jet selection cuts are [p.sub.T](j) > 25 GeV for [absolute value of ([eta](j))] < 2.5 and [p.sub.T](j) > 50 GeV for [absolute value of ([eta](j))] > 2.5. The dominant background contribution of about 90% comes from the SM q[bar.q] [right arrow] [gamma][gamma] process and, in this analysis, we only consider this as the background which we estimate from our simulation. The error bars associated with the background represent the statistical uncertainly only. In reality, various components of systematic uncertainties like the jet energy scale, jet energy resolution, and uncertainty in the luminosity must also be considered to obtain the total uncertainty [33]. But the systematic uncertainty becomes small compared to the statistical one when background distributions are obtained from data.

It is visibly clear that the different production modes display very different jet multiplicity distributions. The distributions for the [gamma][gamma] and the qq fusion modes peak at 0 jets but the peak for the qq mode is not as sharp as the [gamma][gamma] mode. Cross section for the 0-jet bin for the [gamma][gamma] mode is roughly about 60% of the total cross section. On the other hand, it is about 45% for the qq fusion case. The SM background [N.sub.jet] distribution also peaks at 0 jets but contains only 30% of the total cross section. The gg fusion shows a peak at 1 and 2 jets whereas the vector boson fusion production through the WW and ZZ fusions shows peak at 2 and 3 jets. The [gamma]Z fusion mode, on the other hand, shows a peak at 1 jet.

Different nature of the [N.sub.jet] distribution can be captured by the average jet multiplicity associated with the diphoton resonance. We compute the average jet multiplicity of different production modes and the background and report these numbers in Table 1. It is expected that if the scalar is produced through the [gamma][gamma] fusion, then the average jet multiplicity is lower compared to the gg or the qq fusions. This is because, in case of the [gamma][gamma] fusion, a hard jet in the final state can originate from the q [right arrow] q[gamma] splitting. However, this is suppressed compared to the leading order (LO) process with zero jet by the small probability of q [right arrow] q[gamma] branching and also by the smallness of [alpha]. On the other hand, colored particles in the initial state, that is, in case of the gg or the qq fusions, lead to higher jet multiplicity. The average jet multiplicity is greater than two for the vector boson fusion modes because one would anticipate to get at least two hard jets in most of the events since two initial V's come from the q [right arrow] q'V branching. For the [gamma]Z initial state, one expects at least one hard jet from the q [right arrow] qZ splitting.

3.5. Multivariate Analysis. In the previous subsection, we show as a demonstration that the jet multiplicity distributions of two different production modes (and also for the background) can be quite different. Apart from the [N.sub.jet] distribution, there are other kinematic variables which also show some differences in their shapes for different production modes. For example, in [15], the authors showed that various distributions like the scalar sum of transverse energy [H.sub.T], pseudorapidity ([eta]) of the selected photons and jets, central rapidity gap ([DELTA][eta]) between the jets, and the scalar show some visible differences for the gg and the [gamma][gamma] production modes.

A cut-based analysis which employs a set of rectangular cuts may not perform well to decipher the underlying production mechanism of the scalar. In order to effectively distinguish two different production modes, one can use various kinematic variables that show some (small) differences in their shapes simultaneously in a MVA whose output might show large differences in their shapes. If appropriate variables are chosen, a MVA is expected to perform better than a cut-based analysis. Generally, MVA techniques are used to separate signal from background. Here, we use a MVA technique (BDT) to distinguish two different production mechanisms more efficiently than a simple cut-based analysis. In particular, we use the adaptive BDT algorithm in the TMVA framework. We train the algorithm by tuning various parameters like the number of trees, minimum size of the node, and so forth for proper training of different production modes. Optimal values of these parameters are not fixed and they can differ for each analysis.

For MVA, we select events with at least one jet and construct twelve simple kinematic variables as shown in Table 2. This includes [N.sub.jet], [H.sub.T], [DELTA][eta] between [PHI] and leading jet, [eta] and [p.sub.T] of two selected photons, and the leading jet and the separation [DELTA]R in the [eta]-[phi] plane between the photons and the leading jet. These twelve variables are finalized from a bigger set of variables by looking at their discriminatory power and less correlation. In particular, the variables we use are not correlated more than ~40% for signal. But these correlations might be different for the background. Next to each variable in Table 2, we show their relative importance in the BDT response and these numbers are obtained from TMVA using the gg and the [gamma][gamma] production modes. Relative importance is a fraction (where all importance sums up to unity) which is used to identify the ranking of the variables in MVA. In other words, greater relative importance of a variable signifies that the variable is a better discriminator. For actual definition of relative importance, interested readers may look into the TMVA manual. From Table 2, we see that [N.sub.jet] is the best discriminator to differentiate the gg and the [gamma][gamma] production modes. Other variables like [DELTA]R([[gamma].sub.1], [j.sub.1]), [DELTA]R([[gamma].sub.2], [j.sub.1]), [p.sub.T]([j.sub.1]), [eta]([j.sub.1]), and [DELTA][eta]([PHI], [j.sub.1]) also act as good discriminators. Here, our main aim is to distinguish different production mechanisms of the scalar using a suitable MVA. Before arriving to this step, one might be interested to see comparisons of various kinematic distributions for the different signal modes with the background. In Appendix, we show distributions of some input variables for the signal and the background for the interested readers.

It should be remembered that relative importance or in other words the ranking of a variable might change for different production modes and also for different parameters like [M.sub.[PHI]], [square root of s], and so forth which can change the shape of the kinematic distributions. It is important to mention that this set of twelve variables used here may not be the optimal one. One can always improve the analysis with cleverer choices of variables.

In Figure 7, we show the BDT response by comparing two different production modes at a time. The WW and the ZZ fusion modes are very similar in nature and, therefore, it is extremely difficult to distinguish them. We do not consider the ZZ fusion further as it is very much identical to the WW mode. We show, by picking two production modes at a time, ten such possible BDT responses in Figure 7. These responses are substantially different for most of the combinations and therefore can be distinguished very efficiently. We observe that it is hard to distinguish the [gamma][gamma] and the [gamma]Z production modes as their BDT responses are not very different from each other. One should also notice that there are two peaks in the BDT response of the WW mode. This is because two types of different topologies, that is, the associated production and the vector boson fusion, contribute to the [PHI]jj final state. In case of the bimodal distributions like this, one can use two different BDTs that are trained for two different topologies to further improve the analysis. This type of advanced analysis is beyond the scope of this paper.

As a side remark, one should always be careful about overtraining while using the BDT algorithm (or any other algorithm which uses nonlinear cuts). This can happen without the proper choices of the algorithm specific tuning parameters. One can check whether a test sample is overtrained or not by using the Kolmogorov-Smirnov (KS) statistics. Generally, if KS probability lies within the range 0.1 to 0.9 this guarantees that the test sample is not overtrained. For this purpose, one uses two statistically independent samples, one for training and the other for testing.

4. Conclusions

Among the various resonance search channels at the LHC, the diphoton channel is particularly important as this channel provides a comparatively cleaner background. Generally, the diphoton resonance searches at the LHC assume that the resonance is produced from the gg fusion. Apart from the gg fusion production, many BSM theories predict TeV-scale scalars that decay to diphotons can dominantly be produced by other means, namely, through the quark-quark (qq) fusion or through the gauge boson fusions ([gamma][gamma], [gamma]Z, WW, and ZZ). In this paper we consider an effective field theory of a heavy scalar that decays to diphotons. In this model-independent approach, the scalar can be produced in all the possible types mentioned above. We derive the exclusion limits on the mass and the effective couplings of the scalar using the latest 13 TeV ATLAS diphoton resonance search data with L = 36.7 [fb.sup.-1]. While deriving the limits, we consider, for simplicity, only one effective coupling other than [[kappa].sub.[gamma][gamma]] (since we only focus on the diphoton final state) which is nonzero. We have properly taken care of the modified cut efficiencies while recasting the limits set by the ATLAS collaboration. We find that when the scalar is dominantly produced from the [gamma][gamma] fusion, the latest LHC diphoton resonance search data sets limit on the new physics scale [LAMBDA] [greater than or equal to] 20 TeV for the coupling [[kappa].sub.[gamma][gamma]] ~ 1 for [M.sub.[PHI]] ~ 1 TeV.

In future, if a scalar resonance is seen at the LHC in the diphoton channel, the immediate important issue one has to investigate is how the scalar is produced. Some preliminary analyses have already been done in the context of the 750 GeV resonance where it is shown that the jet multiplicity distributions can be very different for the different production modes. In this paper we revisit the issue and show that the average jet multiplicity and the [N.sub.jet] distribution can act as good discriminators. For better discrimination, we use a sophisticated multivariate analysis by combining twelve simple kinematic variables to distinguish one production mechanism from the other. Our analysis shows that one can identify different production mechanisms very efficiently at the LHC.

https://doi.org/10.1155/2018/9150617

Appendix

Distributions for Signal and Background

For the interested readers, we show here the distributions of some input variables for the signal and the background. The signals in Figures 8 and 9 are for the gg and the WW fusion production modes, respectively. These distributions are obtained by applying the selection cuts defined in Section 3.3. The BDT responses for these two production modes with background are presented in Figure 10. We observe that the signal and the background distributions are very different in nature and, therefore, one could use a MVA to isolate the signal from the background. After filtering out the signal events from the background, one can use our method to identify the underlying production mechanism. It is expected that the signal distributions deviate more and more from the background as we increase the resonance mass. Therefore, isolation of the signal from the background becomes easier for heavier resonances. One can, therefore, tune MVA for lower masses and use the same optimized analysis for higher masses, for simplicity.

Notice that there is a second bump around 500 GeV and in the range 300-500 GeV in the background [p.sub.T]([[gamma].sub.1]) and [p.sub.T]([[gamma].sub.2]) distributions, respectively. Similarly, there is a second bump in the background [H.sub.T] distribution around 1000 GeV (this is expected since [H.sub.T] is correlated with the transverse momenta of the photons). This unusual shape of these distributions also leads to the bimodal nature of the background BDT responses in Figure 10. The origin of these peculiar second bumps in the background distributions is due to the selection cuts [E.sub.T]([[gamma].sub.1]) > 0.4M([[gamma].sub.1], [[gamma].sub.2]) and ET([[gamma].sub.2]) > 0.3M([[gamma].sub.1], [[gamma].sub.2]) used to obtained these plots. We have confirmed that these bumps go away with the removal of the abovementioned correlated cuts.

Conflicts of Interest

The author declares that they have no conflicts of interest.

Acknowledgments

The author thanks Valery A. Khoze, Lucian A. Harland-Lang, Rikard Enberg, and Gunnar Ingelman for helpful discussions. This work is supported by the Swedish Research Council under Contracts 621-2011-5107 and 2015-04814 and by the Carl Trygger Foundation under Contract CTS-14:206.

References

[1] M. Aaboud, G. Aad, B. Abbott et al., "Search for resonances in diphoton events at [square root of s] = 13 TeV with the ATLAS detector," High Energy Physics, vol. 1609, 2016.

[2] V. Khachatryan et al., "Search for resonant production of high-mass photon pairs in proton-proton collisions at [square root of s] = 8 and 13 TeV," Physical Review Letters, vol. 117, no. 5, 2016.

[3] A. Strumia, "Interpreting the 750 GeV digamma excess: a review," High Energy Physics--Phenomenology, 2016, https:// arxiv.org/abs/1605.09401.

[4] G. Aad et al., "Search for high-mass diboson resonances with boson-tagged jets in proton-proton collisions at [square root of s] = 8 TeV with the ATLAS detector," Journal of High Energy Physics, vol. 1512, no. 055, 2015.

[5] G. Aad et al., "Search for resonant diboson production in the llq[bar.q] final state in pp collisions at [square root of s] = 8 TeV with the ATLAS detector," The European Physical Journal C, vol. 75, p. 69, 2015.

[6] G. Aad et al., "Search for production of WW/WZ resonances decaying to a lepton, neutrino and jets in pp collisions at[square root of s] = 8 TeV with the ATLAS detector," The European Physical Journal C, vol. 75, no. 5, p. 209, 2015, Erratum: [The European Physical Journal C vol. 75, pp. 370, 2015].

[7] G. Aad et al., "Observation of a new particle in the search for the Standard Model Higgs boson with the ATLAS detector at the LHC," Physics Letters B, vol. 716, pp. 1-29, 2012.

[8] S. Chatrchyan et al., "Observation of a new boson at a mass of 125 GeV with the CMS experiment at the LHC," Physics Letters B, vol. 716, p. 30, 2012.

[9] L. D. Landau, "On the angular momentum of a system of two photons," Doklady Akademii Nauk Seriya Fizicheskaya, vol. 60, no. 2, pp. 207-209, 1948.

[10] C. N. Yang, "Selection rules for the dematerialization of a particle into two photons," Physical Review A, vol. 77, no. 2, pp. 242-245, 1950.

[11] N. Arkani-Hamed, S. Dimopoulos, and G. Dvali, "The hierarchy problem and new dimensions at a millimeter," Physics Letters B, vol. 429, no. 3-4, pp. 263-272, 1998.

[12] L. Randall and R. Sundrum, "Large mass hierarchy from a small extra dimension," Physical Review Letters, vol. 83, no. 17, pp. 3370-3373, 1999.

[13] M. Aaboud et al., "Search for new phenomena in dijet events using 37 [fb.sup.-1] of pp collision data collected at [square root of s] = 13 TeV with the ATLAS detector," Physical Review D, vol. 96, no. 5, Article ID 052004, 2017.

[14] CMS Collaboration, "Searches for dijet resonances in pp collisions at [square root of s] = 13 TeV using data collected in 2016," CMS-PASEXO-16-056.

[15] C. Csaki, J. Hubisz, S. Lombardo, and J. Terning, "Gluon versus photon production of a 750 GeV diphoton resonance," Physical Review D, vol. 93, no. 9, 2016.

[16] L. A. Harland-Lang, V. A. Khoze, and M. G. Ryskin, "The production of a diphoton resonance via photon-photon fusion," Journal of High Energy Physics, vol. 1603, no. 182, 2016.

[17] M. A. Ebert, S. Liebler, I. Moult et al., "Exploiting jet binning to identify the initial state of high-mass resonances," Physical Review D, vol. 94, no. 5, 2016.

[18] M. Dalchenko, B. Dutta, Y. Gao, T. Ghosh, and T. Kamon, "Exploring the jet multiplicity in the 750 GeV diphoton excess," Physics Letters B, vol. 761, pp. 77-80, 2016.

[19] B. Fuks, D. W. Kang, S. C. Park, and M. Seo, "Investigating the jet activity accompanying the production at the LHC of a massive scalar particle decaying into photons," Physics Letters B, vol. 761, pp. 344-349, 2016.

[20] A. Alloul, N. D. Christensen, C. Degrande, C. Duhr, and B. Fuks, "FeynRules 2.0--A complete toolbox for tree-level phenomenology," Computer Physics Communications, vol. 185, no. 8, pp. 2250-2300, 2014.

[21] C. Degrande, C. Duhr, B. Fuks, D. Grellscheid, O. Mattelaer, and T. Reiter, "UFO--the universal FeynRules output," Computer Physics Communications, vol. 183, no. 6, pp. 1201-1214, 2012.

[22] J. Alwall, R. Frederix, S. Frixione et al., "The automated computation of tree-level and next-to-leading order differential cross sections, and their matching to parton shower simulations," Journal of High Energy Physics, vol. 1407, no. 079, 2014.

[23] L. A. Harland-Lang, A. D. Martin, P. Motylinski, and R. S. Thorne, "Parton distributions in the LHC era: MMHT 2014 PDFs," The European Physical Journal C, vol. 75, no. 5, p. 204, 2015.

[24] L. A. Harland-Lang, V. A. Khoze, and M. G. Ryskin, "The photon PDF in events with rapidity gaps," The European Physical Journal C, vol. 76, no. 5, p. 255, 2016.

[25] T. Sjostrand, S. Mrenna, and P. Z. Skands, "A brief introduction to PYTHIA8.1," Computer Physics Communications, vol. 178, no. 11, pp. 852-867, 2008.

[26] J. de Favereau et al., "DELPHES 3: a modular framework for fast simulation of a generic collider experiment," Journal of High Energy Physics, vol. 1402, no. 057, 2014.

[27] M. Cacciari, G. P. Salam, and G. Soyez, "FastJet user manual," The European Physical Journal C, vol. 72, no. 3, article 1896, pp. 1-54, 2012.

[28] M. Cacciari, G. P. Salam, and G. Soyez, "The Anti-k(t) jet clustering algorithm," Journal of High Energy Physics, vol. 0804, 2008.

[29] L. Bartoszek et al., "Search for new phenomena in high-mass diphoton final states using 37 [fb.sup.-1] of proton-proton collisions collected at [square root of s] = 13 TeV with the ATLAS detector," Physics Letters B, vol. 775, pp. 105-125, 2014.

[30] A. Hoecker, P. Speckmayer, J. Stelzer et al., "TMVA--Toolkit for Multivariate Data Analysis," PoS ACAT, vol. 040, 2007, https://arxiv.org/abs/physics/0703039.

[31] U. Danielsson, R. Enberg, G. Ingelman, and T. Mandal, "Heavy photophilic scalar at the LHC from a varying electromagnetic coupling," Nuclear Physics B, vol. 919, pp. 569-582, 2017.

[32] M. L. Mangano, M. Moretti, F. Piccinini, and M. Treccani, "Matching matrix elements and shower evolution for top-pair production in hadronic collisions," Journal of High Energy Physics, vol. 701, no. 013, 2007.

[33] V. Khachatryan, A. M. Sirunyan, A. Tumasyan et al., "Search for high-mass diphoton resonances in proton-proton collisions at 13 TeV and combination with 8 TeV search," Physics Letters B, vol. 767, no. 1, article 147, 2017.

[34] A. Kruse, A. S. Cornell, M. Kumar, B. Mellado, and X. Ruan, "Probing the Higgs boson via vector boson fusion with single jet tagging at the LHC," Physical Review D, vol. 91, no. 5, Article ID 053009, 2015.

[35] D. de Florian, C. Grojean, F. Maltoni et al., "Handbook of LHC higgs cross sections: 4. deciphering the nature of the higgs sector," LHC Higgs Cross Section Working Group, p. 868, 2016, https://arxiv.org/abs/1610.07922.

[36] C. Balazs, H. J. He, and C. P. Yuan, "QCD corrections to scalar production via heavy quark fusion at hadron colliders," Physical Review D, vol. 60, Article ID 114001, 1999.

[37] K. Arnold, M. Bahr, G. Bozzi et al., "VBFNLO: A parton level Monte Carlo for processes with electroweak bosons," Computer Physics Communications, vol. 180, pp. 1661-1670, 2009.

Tanumoy Mandal (iD)

Department of Physics and Astronomy, Uppsala University, Box 516, 751 20 Uppsala, Sweden

Correspondence should be addressed to Tanumoy Mandal; tanumoy.mandal@physics.uu.se

Received 7 July 2017; Revised 21 November 2017; Accepted 11 December 2017; Published 9 January 2018

Academic Editor: Luca Stanco

Caption: Figure 1: Partonic cross sections of various production modes of [PHI] as functions of [M.sub.[PHI]] computed at [[mu].sub.R] = [[mu].sub.F] = [M.sub.[PHI]] at the 13 TeV LHC. Here, [[sigma].sup.0.sub.[PHI]xy] denotes the cross section of [PHI] produced through the xy fusion with [[kappa].sub.xy] = 1 and [LAMBDA] = 1 TeV. Initial g, q, and [gamma] come from the PDFs of proton whereas initialWand Z come from initial quarks. These cross sections are computed by applying some basic generation level cuts as defined in (4).

Caption: Figure 2: Cut efficiencies for different production modes as functions of [M.sub.[phi]] for the ATLAS selection cuts as defined in the text. The ATLAS cut efficiency (blue dots) is for a spin-0 resonance produced through the gg fusion [29].

Caption: Figure 3: Exclusion regions using the observed upper limits (with 95% confidence level) on [sigma] x BR for the spin-0 diphoton resonance search by ATLAS [29] ([square root of s] = 13 TeV, L = 36.7 [fb.sup.-1]). (a) Exclusions in the [M.sub.[PHI]]-[[kappa].sub.[gamma][gamma]]/[LAMBDA] plane setting all other [[kappa].sub.xy] = 0. (b)-(f) Exclusions in the [[kappa].sub.xy]/[LAMBDA]-[[kappa].sub.[gamma][gamma]]/[LAMBDA] plane while setting all other [[kappa].sub.xy] to zero for [M.sub.[PHI]] = 0.75, 1, 1.5, 2 TeV.

Caption: Figure 4: Combined cut efficiency in the [[kappa].sub.gg]/[LAMBDA]-[[kappa].sub.[gamma][gamma]]/[LAMBDA] plane (all other [[kappa].sub.xy] are zero) for [M.sub.[PHI]] = 1 TeV for the ATLAS selection cuts as defined in the text.

Caption: Figure 5: 2D plots of [sigma], [sigma] x BR, and [sigma] x BR x [epsilon] for [M.sub.[PHI]] = 1 TeV for two coupling assumptions: [[kappa].sub.gg], [[kappa].sub.[gamma][gamma]] [not equal to] 0 (a) and [[kappa].sub.WW], [[kappa].sub.[gamma][gamma]] [not equal to] 0 (b) while setting other [[kappa].sub.xy] to zero.

Caption: Figure 6: Jet multiplicity ([N.sub.jet]) distributions of various production modes of the scalar [PHI] for [M.sub.[PHI]] = 1 TeV at the 13 TeV LHC. The [N.sub.jet] distribution for the dominant SM background, that is, for the pp [right arrow] [gamma][gamma] process, is obtained from our simulation. The uncertainty associated with this distribution is statistical uncertainty only which is computed for L = 50 [fb.sup.-1].

Caption: Figure 7: BDT responses of two different production modes at a time. These are obtained using a MVA with twelve input variables as listed in Table 2 for [M.sub.[PHI]] = 1 TeV at the 13 TeV LHC. Diphoton events are selected by applying [absolute value of ([M.sub.[gamma][gamma]] - [M.sub.[PHI]])] < 150 GeV cut in addition to the ATLAS selection cuts as listed in the text.

Caption: Figure 8: Normalized distributions of a few sample input variables for the signal (gg fusion with [M.sub.[PHI]] = 1 TeV) and background at the LHC ([square root of s] = 13 TeV). These distributions are obtained by applying the selection cuts defined in Section 3.3.

Caption: Figure 9: Normalized distributions of a few sample input variables for the signal (WW fusion with [M.sub.[PHI]] = 1 TeV) and background at the LHC ([square root of s] = 13 TeV). These distributions are obtained by applying the selection cuts defined in Section 3.3.

Caption: Figure 10: The BDT response for the signal ((a) gg fusion and (b)WW fusion) and the background for [M.sub.[PHI]] = 1 TeV at the [square root of s] = 13 TeV LHC.

Table 1: Average jet multiplicity for different production modes of [PHI] at the 13 TeV LHC with L = 50 [fb.sup.-1]. Average jet multiplicity is defined by the sum, [[summation].sub.i][(BH).sub.i] [N.sub.i] where B[H.sub.i] represents the bin height of the ith bin of the normalized [N.sub.jet] distribution and [N.sub.i] = i - 1 is the number of jets associated with the ith bin. gg [right qq [right [gamma][gamma] [gamma]Z [right arrow] [PHI] arrow] [PHI] [right arrow][PHI] arrow] [PHI] 1.94 0.92 0.61 1.75 WW [right ZZ [right arrow] [PHI] arrow] [PHI] [gamma][gamma] SM 2.85 2.87 1.42 Table 2: Input variables used for MVA to separate gg and [gamma][gamma] production modes and their relative importance. Variable Importance [N.sub.jet] 1.8 x [10.sup.-1] [H.sub.T] 4.5 x [10.sup.-2] [DELTA][eta]([PHI], [j.sub.1]) 8.1 x [10.sup.-2] [p.sub.T]([[gamma].sub.1]) 4.9 x [10.sup.-2] [p.sub.T]([[gamma].sub.2]) 6.0 x [10.sup.-2] [p.sub.T]([j.sub.1]) 9.1 x [10.sup.-2] [eta]([[gamma].sub.1]) 7.5 x [10.sup.-2] [eta]([[gamma].sub.2]) 6.2 x [10.sup.-2] [eta]([j.sub.1]) 9.4 x [10.sup.-2] [DELTA]R([[gamma].sub.1], [[gamma].sub.2]) 6.7 x [10.sup.-2] [DELTA]R([[gamma].sub.1], [j.sub.1]) 9.9 x [10.sup.-2] [DELTA]R([[gamma].sub.2], [j.sub.1]) 1.0 x [10.sup.-1]

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Mandal, Tanumoy |

Publication: | Advances in High Energy Physics |

Article Type: | Report |

Geographic Code: | 1USA |

Date: | Jan 1, 2018 |

Words: | 8594 |

Previous Article: | New Possibilities of Hybrid Texture of Neutrino Mass Matrix. |

Next Article: | Spread of Entanglement in Non-Relativistic Theories. |

Topics: |