Arxiv:2011.07560V4 [Quant-Ph] 9 Jul 2021 Aharonov Et Al

KEK-TH-2272

Weak Value Ampliﬁcation in High Energy Physics: A Case Study for Precision Measurement of CP Violation in B Meson Decays

Satoshi Higashino,1, ∗ Yuichiro Mori,2, † Yosuke Takubo,3, 4, ‡ Takeo Higuchi,5, § Akimasa Ishikawa,3, 4, 6, ¶ and Izumi Tsutsui2, 3, 4, ∗∗ 1Department of Physics, Kobe University, Hyogo 657-8501, Japan. 2Department of Physics, the University of Tokyo, Tokyo 113-0033, Japan. 3Institute of Particle and Nuclear Studies, High Energy Accelerator Research Organization (KEK), Ibaraki 305-0801, Japan. 4The Graduate University for Advanced Studies (SOKENDAI), Hayama 240-0193, Japan. 5Kavli Institute for the Physics and Mathematics of the Universe (WPI), the University of Tokyo, Kashiwa 277-8583, Japan. 6International Center for Elementary Particle Physics, the University of Tokyo, Tokyo 113-0033, Japan. (Dated: July 12, 2021) The technique of weak value amplification, proposed by Aharonov et al. in 1988, has been applied for various fields of physics for the purpose of precision measurement, which is made possible by exploiting the freedom of ‘postselection’ specifying actively the final state in the physical process. Here we report the feasibility of utilizing the technique of weak value amplification in high energy particle physics, especially in measuring the CP -violating parameters in B meson decays, where the effective lifetime of the decay mode is expected to be prolonged statistically due to the postselection. Our analysis shows that, when adopted in the Belle II experiment at the SuperKEKB collider, the effective lifetime may be prolonged up to 2.6 times, and that the measurement precision of the CP -violating parameters will also be improved by its effect.

I. INTRODUCTION treats them equally. On account of their active nature, the choices of the two states, |ψi and |φi, are referred to In physical sciences, measurement is arguably the most as preselection and postselection, respectively. It is then ˆ fundamental practice to find something new in a given expected that, given an observable A, one may consider theoretical framework. In quantum physics, which ap- the weak value plies to systems of diverse scales ranging from elementary hφ|Aˆ|ψi particles to condensed matter and even to the early uni- Aw = (1) hφ|ψi verse, the framework is provided by Hermitian operators of observables serving as basic tools to make theoretical as a tangible physical quantity, and that this weak value predictions under initially prepared quantum states. It can be observed by the procedure called weak measure- has been firmly believed that the basis of these predic- ment in which one implements the postselection at the tions must be laid on the eigenvalues of the operators end of the measurement in addtion to the preselection or their statistical average, and the measurements are made at the begining. Over the last two decades, this designed to compare these numbers with the outcomes has been confirmed in various experiments, yielding ma- obtained in the actual experiments. Concomitantly, our jor ramifications in two respects, one is philosophical and very notion of physical reality has been tightly tied to the other practical. the eigenvalues of these operators. On the philosophical side, the weak value offers novel An initiative to reconsider this conventional wisdom interpretations of physical reality in paradoxical phenom- was put forward in 1988 when the novel notion of phys- ena such as the double slit experiment involving the wave- ical quantity, known as weak value, was proposed by particle duality [2, 3], the decoupling of spin and position arXiv:2011.07560v4 [quant-ph] 9 Jul 2021 Aharonov et al. in a time-symmetric formulation of quan- degrees of freedom of neutrons [4] or the three box para- tum theory [1]. This formulation is an alternative but dox where probabilities may become negative [5] along equivalent in content to the standard one, and the idea with the related issue of the photon path [6, 7]. On the behind is quite simple: one just replaces the notion of practical side, which is directly related to the topic of ‘state’ by ‘process’, that is, instead of choosing only the the present paper, the weak measurement can be used initial state |ψi one also chooses the final state |φi and for precision measurement of parameters that are other- wise difficult to detect. This intriguing possibility arises thanks to the freedom in implementing the postselection in relation to the preselection. To be more explicit, the ∗ [email protected] weak value A can possibly be amplified by choosing † yuichiro@post.kek.jp w ‡ [email protected] the two states |ψi and |φi such that the denominator of § [email protected] Aw is rendered small while the numerator is kept finite. ¶ [email protected] This possibility, which was already suggested in [1], has ∗∗ [email protected] actually been utilized in recent years for a number of 2 purposes, including the ultrasensitive detection of beam the difference of decay probabilities of the initial flavor deflection and the observation of the spin Hall effect of eigenstates, the effective lifetime of the mass eigenstates light [8, 9]. A more recent study reports that the sen- can be expected to be prolonged by a proper postselec- sitivity has been further improved by several orders of tion of the final state |φi. This may also improve the sen- magnitude [10]. sitivity in the measurement of CP -violating parameters, Despite its potential universality, however, the applica- which should be studied together with other contributing tion of weak measurement appears to have been limited factors that entail the postselection. mostly to optical systems so far, with the exception of The Belle II experiment [18] is operated at the Su- neutrons [4] and atoms [11]. One may therefore wonder perKEKB [19] electron-positron collider, where one of the if the weak measurement can also be put to use in high objectives is to measure the CP -violating parameters in energy particle physics whose energy scale is far from rare decays of the B mesons, with the primary aim to those previously applied. In fact, we notice that, even observe new physics beyond the Standard Model (SM). before its application in quantum optics, the weak mea- In this experiment, a huge amount of B0B¯0 pairs will surement has already been in use for long time in particle be produced from the resonance decay of Υ(4S), which physics, given that the initial and final states in particle admits a big advantage in performing the required post- scattering may be regarded precisely as the states cho- selection of the final states under a large reduction of sen for the preselection and postselection. On the other statistics. hand, we also notice that the freedom of postselection In this paper, we shall argue that the freedom of posts- has never been utilized for the particular purpose of am- election may be exploited in the B0 → K∗0γ decay chan- plifying the prospected signals before. The reason for nel if we can choose the polarization of the photon to this may be ascribed to the fact that the possible range meet a certain consistency condition. As a result, we of final states in high energy experiments appears to be find that the effective lifetime of the B meson can be dictated by nature, not by experimenters, leaving no free- prolonged up to 2.6 times, and that the precision mea- dom of postselection to be adjusted. Besides, it has not surement of the CP -violating parameters can also be im- been clear under what context the weak value becomes proved, under the assumption of the realistic signal yield relevant in high energy physics. and the background contamination as well as the sys- One of the aims of the present paper is to point out tematic errors used in the Belle [14] and Belle II experi- that in general there do exist such a context in systems ment [20]. involving decaying/oscillatory phenomena among differ- This paper is organized as follows. In Section II we ent species of particles. One is then allowed to consider first discuss how the weak value arises in the context the weak value of the Hamiltonian (or energy), and if of particle decay when the final states are chosen, and the freedom of postselection can be exploited to achieve thereby see the reason why the effective lifetime of the the amplification of the weak value, one finds an effective particle can be prolonged by the weak measurement. We prolongation of lifetime of the decaying particle. Obvi- then present in Section III a theoretical model of the ously, such systems are not uncommon in atomic and B mesons, where their effective lifetime is compared be- high energy physics, with the B meson system being a tween the two cases with and without the postselection. prime example. Indeed, one can recognize, in retrospect, The realization of the choice at the decayed stage, which that the weak measurement has implicitly been realized is a new aspect of weak measurements, is discussed for 0 ∗0 at the B factories with various choices of the decay modes the decay channel B → K γ. Our simulation study as final states (see, e.g., [12–14]). Another aim of this pa- in the experiment at the SuperKEKB collider is given in per is to furnish a case study on the consequence of weak Section IV, where we argue that, if the present method of measurement when the freedom of postselection is fully weak measurement is adopted in the Belle II experiment, available in the B meson system used as an archetype. we can expect the aforementioned significant prolongation of the effective lifetime along with the increase in More specifically, here we discuss how the lifetime of the precision. Finally, Section V is devoted to our con- the B meson can be prolonged effectively by the weak d clusion and discussions on remaining issues and future value amplification, and also argue the prospect of im- prospects in the line of this research. proving the sensitivity in measuring CP -violating parameters. Throughout the paper we call Bd meson as B meson for brevity, and consider only the CP violation induced by flavor mixing [15–17]. As the mass eigenstates II. WEAK VALUE IN PARTICLE DECAY of the B meson (|BLi and |BH i) are inseparable due to the magnitude of the decay width comparable with their Before we delve into the discussion of applying the mass difference, one can expect that the effective lifetime weak measurement in high energy particle physics, we of the B meson provides a suitable testing ground for describe the basic reason why the weak value arises in the feasibility study of the weak measurement in parti- analyzing the particle decay in general. cle physics. In fact, since the mass eigenstates |BLi and To this end, we first note that the weak value Aw given |BH i are both superpositions of the flavor eigenstates in (1) is in general a complex number, and that it may |B0i and |B¯0i, and since the CP violation manifests in take any value one wishes, provided that the pair of states 3

|ψi and |φi can be chosen freely. This condition is in fact where Aw is the weak value (1) of the normalized operator the presumption of the time-symmetric formulation [1] Aˆ in (3). Plugging this expression into (2), we obtain where the process |ψi → |φi is specified by selecting both of the states |ψi and |φi independently and freely, rather Pr(∆t | ψ → φ) ' e−Γ∆t| hφ|ψi |2e2gIm[Aw]. (7) than letting the latter determined from the former by causal time development. For this reason, the states |ψi This shows that, as long as the linear approximation is and |φi are called preselected state and postselected state, valid, the probability distribution of the decay may be respectively. affected by the weak value Aw (its imaginary part). No- In the context of particle decay, one may expect that tably, by choosing the combination of the preselection the weak value amplification can be used to prolong effec- and the postselection properly, the weak value Aw can tively the decay time of the initial particle by considering be altered freely; for instance, Im[Aw] may be rendered the Hamiltonian operator Hˆ for the observable Aˆ. That large by choosing hφ|ψi small as stated in Section I. Even this is indeed the case can be seen as follows. Let |ψi be if the coupling constant g cannot be regarded small, the the preselected state of the initial particle and |φi be the system may allow one to evaluate the transition ampli- postselected state of the decayed particle(s). The proba- tude in full order to see if the selections affect the out- bility of transition between them during the time period come. Later, we shall see explicitly that the decay of the ∆t will then be given by B meson falls into that particular case.

Pr(∆t | ψ → φ) = |hφ|e−i∆tHˆ |ψi |2. (2) III. THEORETICAL BASIS Suppose that the Hamiltonian Hˆ is bounded as is usually the case for considering transitions among a finite set In this section, we present our theoretical basis for ap- of energy levels. Noting that Hˆ is not Hermitian (or plying the weak measurement for the detection of CP self-adjoint) to take account of the decaying processes, violation in the B meson system. In our case, the Hamil- we let m and Γ be the averages of the real part and tonian Hˆ takes the role of the observable Aˆ, and we shall the imaginary part of the eigenvalues of Hˆ , respectively. see that the imaginary part of the weak value Aw is re- To extract only the effects of interference between the lated to the effective lifetime of the B meson. The aim of components, we may define the ‘normalized’ operator of the weak measurement is then to find a source of CP vio- time development by lation in the conditional time distribution when our state is |B0i at the initial time and is an arbitrary chosen state 2 i |B i at the time of the decay (see (27)). Aˆ = Hˆ − m − Γ , (3) decay ∆m 2 where ∆m is the range of the real part, i.e., the difference A. Time evolution between the largest and the smallest eigenvalues in the real part. In particular, if there are only two states which The state |ψi which describes a neutral B meson is in enter in the transition, then the operator Aˆ in (3) may general a superposition of the flavor eigenstates |B0i and admit the form of the familiar Pauli matrices, especially |B¯0i which are the CP conjugate to each other, that is, ˆ 0 0 0 0 A = σ3 in the energy diagonal basis. |B¯ i = CP|B i and |B i = CP|B¯ i under the CP trans- Substituting Hˆ with the normalized one Aˆ in the tran- formation CP for which (CP)2 = I. Since mesons are sition amplitude, we find unstable particles, the time evolution is phenomenologi- cally described by a non-Hermitian Hamiltonian Hˆ . The −i∆tHˆ − 1 Γ∆t −im∆t −igAˆ hφ|e |ψi = e 2 e hφ|e |ψi , (4) eigenstates of the Hamiltonian are given by where i Hˆ |B i = m − Γ |B i, (8) L L 2 L L 1 1 ∆m g = ∆m∆t = Γ∆t (5) i 2 2 Γ Hˆ |B i = m − Γ |B i, (9) H H 2 H H is the effective coupling constant. Now, from (4) we realize that the effective range of the time period ∆t may where mL, mH , ΓL, ΓH are real positive numbers satisfy- be restricted by the order O(Γ∆t) = 1, and hence if ing O(∆m/Γ) 1, then g can be regarded as a small pa- m < m . (10) rameter O(g) 1. This allows us to employ the linear L H approximation, These eigenstates |BLi and |BH i can be written as −igAˆ ˆ 0 0 hφ|e |ψi ' hφ|I − igA|ψi |BLi = p|B i + q|B¯ i, (11)

−igAw 0 0 = hφ|ψi (1 − igAw) ' hφ|ψi e , (6) |BH i = p|B i − q|B¯ i, (12) 4 where p, q are complex numbers fulfilling |p|2 + |q|2 = 1. confirmed from (8) and (9) by expressing the Hamilto- It should be noted that these eigenstates are not orthog- nian in terms of the parameters p, q in the orthogonal onal, and hence the Hamiltonian Hˆ cannot be diagonal- basis of the flavor eigenstates |B0i and |B¯0i as ized. i Using (11) and (12), we find Hˆ (p, q) = (m − Γ) |B0ihB0| + |B¯0ihB¯0| 2 1 i p q |B0i = (|B i + |B i) , (13) − (∆m − ∆Γ) |B0ihB¯0| + |B¯0ihB0| , 2p L H 2 2q 2p 1 (17) |B¯0i = (|B i − |B i) . (14) 2q L H using In high energy experiments, we decide whether the initial mL + mH ΓL + ΓH state is the flavor eigenstate |B0i or its CP conjugate m = , Γ = , (18) 2 2 |B¯0i through a process called the flavor tagging. In this paper, the time at which this tagging is performed is and referred to as ‘initial’, and the time lapse from the initial time will be denoted by ∆t. When the meson is identified ∆m = mH − mL, ∆Γ = ΓH − ΓL. (19) as B0 at the initial time, it evolves into the state One then immediatety observes that the CP conjuga- ˆ tion of the Hamiltonian results in the interchange of |B0(∆t)i = e−i∆tH |B0i (15) the parameters, CPHˆ (p, q) CP = Hˆ (q, p). It then follows that, given the transition amplitude T (p, q) = at the time of the decay. Similarly, when the particle is −i∆tHˆ (p,q) identified as B¯0, we have hα|e |βi, the corresponding amplitude of the CP conjugate has the same property, 0 −i∆tHˆ 0 |B¯ (∆t)i = e |B¯ i. (16) ˆ ˆ T¯(p, q) = hα¯|e−i∆tH(p,q)|β¯i = hα|CPe−i∆tH(p,q)CP|βi ˆ In our discussion, for the preselected state |ψi we pri- = hα|e−i∆t CPH(p,q) CP |βi = T (q, p), (20) marily choose |ψi = |B0i except when |B¯0i is considered for comparison. Our choice for the postselected state |φi as claimed. will be given later in Subsection C. B. Time distribution without postselection For our later convenience, we note here that, as long as the dynamical development is concerned, the evolu- Due to the instability of the particle, the norm of the tion (16) of the CP conjugate B¯0 is obtained immediately particle’s state is not always unity, as can be seen explic- from (15) by interchanging p ↔ q. This can be readily itly as

† hB0(∆t)|B0(∆t)i = hB0|ei∆tHˆ e−i∆tHˆ |B0i

1 im ∆t− ΓL ∆t im ∆t− ΓH ∆t −im ∆t− ΓL ∆t −im ∆t− ΓH ∆t = e L 2 hB | + e H 2 hB | · e L 2 |B i + e H 2 |B i 4|p|2 L H L H 1 e−Γ∆t = e−ΓL∆t + e−ΓH ∆t + ei∆m∆thB |B i + e−i∆m∆thB |B i 4|p|2 4|p|2 L H H L 1 |p|2 − |q|2 = e−ΓL∆t + e−ΓH ∆t + e−Γ∆t cos (∆m∆t) . (21) 4|p|2 2|p|2

Similarly, when the initial state is |B¯0i, we ﬁnd

1 |p|2 − |q|2 hB¯0(∆t)|B¯0(∆t)i = e−ΓL∆t + e−ΓH ∆t − e−Γ∆t cos (∆m∆t) . (22) 4|q|2 2|q|2

To clarify the meaning of (21), we insert the identity into the left hand side of (21) to obtain hB0(∆t)|B0(∆t)i = |hB0|B0(∆t)i|2 + |hB¯0|B0(∆t)i|2. I = |B0i hB0| + |B¯0ihB¯0| (23) (24) 5

This shows that (21) provides the probability of remain- Since the actual values of ΓL and ΓH are both close√ to ing either as B0 or B¯0 without decaying until the time ∆t Γ, it is assuring to observe that, for |p| ≈ |q| ≈ 1/ 2, 0 0 0 after the particle starts as B at the initial time. Note we have τeff (B ) ≈ 1/Γ, which is the lifetime of the B that the result (22) is also obtained from (21) by the meson. interchange p ↔ q as we noted before. Taking the normalization condition of probability into account, we observe that the time distribution P (∆t) ex- C. Time distribution with postselection pressing the probability density of decay at time ∆t reads Now we wish to implement the postselection and see how it affects the time distribution obtained before with- d 0 0 0 0 P (∆t) = − 0 hB (∆t )|B (∆t )i out the postselection. The postselected state |φi can in d∆t 0 ∆t =∆t general be written as a superposition of the two states, 1 0 0 −ΓL∆t −ΓH ∆t |B i and |B¯ i. As we argue later in Subsection E, such = 2 ΓLe + ΓH e 4|p| a postselected state may be realized as a state |Bdecayi |p|2 − |q|2 from a particular decay mode f measured in the exper- + e−Γ∆t∆m sin (∆m∆t) 2|p|2 iment. This will be our postselected state, namely, our postselected state |φi = |B i is given by |p|2 − |q|2 decay + e−Γ∆tΓ cos (∆m∆t) . (25) 2 0 0 2|p| |Bdecayi = r|B i + s|B¯ i, (27) where r, s are complex numbers satisfying the normaliza- We note in passing that, if we implement the postselection condition, tion further, we will obtain the conditional probability Pr(∆t|ψ → φ) in (2). The result P (∆t) in (25) will then |r|2 + |s|2 = 1. (28) be derived, up to normalization, by summing it up for all |φi states. These two parameters r and s are assumed to be chosen The effective lifetime of this particle can then be eval- freely by the measurement process of the decay mode f. uated as In addition, when investigating CP violation, the CP conjugate of the state in (27), Z ∞ 0 0 0 0 0 0 τeff (B ) = d∆t ∆t P (∆t ) |B¯decayi = s|B i + r|B¯ i, (29) 0 Z ∞ d should also be considered. From now on, ∆t will be = − d∆t0 ∆t0 hB0(τ)|B0(τ)i 0 dτ τ=∆t0 considered to be the time of the decay, which implies 0 0 0 that |B (∆t)i and |B¯ (∆t)i in (15) and (16) are now the 0 0 0 0 0 ∆t =∞ = + ∆t hB (∆t )|B (∆t )i ∆t0=0 states realized by dynamical evolution at the time of the Z ∞ 0 0 0 0 0 decay. + d∆t hB (∆t )|B (∆t )i Before evaluating the time distribution with postselec- 0 tion, we first obtain the probability of finding the post- 1 1 1 |p|2 − |q|2 Γ selected state |B i specified in (27) when the state = 2 + + 2 2 2 . decay 4|p| ΓL ΓH 2|p| Γ + ∆m under consideration is the B meson state |B0(∆t)i given (26) in (15). The result is

0 2 ∗ 0 0 ∗ ¯0 0 2 hBdecay|B (∆t)i = r hB |B (∆t)i + s hB |B (∆t)i ∗ ∗ ∗ ∗ 2 r qs − ΓL ∆t −i∆tm r qs − ΓH ∆t −i∆tm = + e 2 e L + − e 2 e H . (30) 2 2p 2 2p

Introducing the relative phase parameters ϕ and θ we ﬁnd through 0 2 hBdecay|B (∆t)i |r|2 |q|2|s|2 e−ΓL∆t + e−ΓH ∆t = + 2 2|p|2 2 p |p| r |r| = eiϕ, = eiθ, (31) |r|2 |q|2|s|2 q |q| s |s| + − e−Γ∆t cos (∆m∆t) 2 2|p|2 |r||q||s| e−ΓL∆t − e−ΓH ∆t + cos (θ − ϕ) |p| 2 |r||q||s| − sin (θ − ϕ) e−Γ∆t sin (∆m∆t) . (32) |p| 6

¯ ¯0 2 The corresponding CP conjugate hBdecay|B (∆t)i is by then obtained simply by putting |p| ↔ |q| and ϕ ↔ −ϕ ¯ ¯0 2 in (32). 0 hBdecay|B (∆t)i P (∆t|B¯ → B¯decay) = . R ∞ 0 ¯ ¯0 0 2 d∆t hBdecay|B (∆t )i It is notable that either the second or the fourth term 0 (34) in (32) is nonvanishing for any combinations of the two parameters (r, s) specifying the postselection except for Using (32), we obtain the denominator of (33) as the case |p| = |q|. These terms have oscillatory property which can be controlled by the parameters |r| and θ. Z ∞ 0 0 0 2 d∆t hBdecay|B (∆t )i The time distribution with postselection is then de- 0 |r|2 |q|2|s|2 1 1 ﬁned by = + 2 + 2 2|p| 2ΓL 2ΓH |r|2 |q|2|s|2 Γ + − 2 2|p|2 Γ2 + ∆m2 0 2 0 hBdecay|B (∆t)i |r||q||s| 1 1 P (∆t|B → Bdecay) = R ∞ 0 0 0 2 + cos (θ − ϕ) − 0 d∆t |hBdecay|B (∆t )i| |p| 2ΓL 2ΓH (33) |r||q||s| ∆m − sin (θ − ϕ) . (35) |p| Γ2 + ∆m2

In particular, when the decay widths of the two mass which is the conditional time distribution with B0 at the eigenstates of the meson are considered to be equivalent, initial time and Bdecay at the time of the decay. Similarly, |ΓL −ΓH |/Γ 1 (we have |ΓL −ΓH |/Γ < 0.01 according 0 its CP conjugate (B¯ → B¯decay) version is also deﬁned to [21]), the time distribution (33) admits the form,

√ |r|2 |q|2(1−|r|2) |r|2 |q|2(1−|r|2) |q||r| 1−|r|2 + 2 + − 2 cos (∆m∆t) − sin (θ − ϕ) sin (∆m∆t) 0 −Γ∆t 2 2|p| 2 2|p| |p| P (∆t|B → Bdecay) = Γe √ . |r|2 |q|2(1−|r|2) |r|2 |q|2(1−|r|2) Γ2 |q||r| 1−|r|2 Γ∆m 2 + 2|p|2 + 2 − 2|p|2 Γ2+∆m2 − |p| sin (θ − ϕ) Γ2+∆m2 (36)

0 0 Figure 1 shows P (∆t|B → Bdecay) for various |r| time with the condition B → Bdecay, is obtained by values in the case of |p| = |q|, ∆m = 0.506 (ps)−1, 1/Γ = 1.519 ps, and the postselection parameter θ is ﬁxed to be 0. The time distribution is clearly shifted to right as |r| goes to zero, indicating the increase of the eﬀective lifetime.

D. Weak value ampliﬁcation

Z ∞ 0 0 0 0 0 τeff (B → Bdecay) = d∆t ∆t P (∆t |B → Bdecay), Having seen that different time distributions can be ob- 0 tained depending on whether the postselection was per- R ∞ 0 0 0 0 2 d∆t ∆t hBdecay|B (∆t )i formed or not, we now evaluate the effective lifetime of = 0 . R ∞ d∆t0 |hB |B0(∆t0)i|2 the B meson with the postselection. Actually, the mean 0 decay value of the time distribution, which is the effective life- (37) 7

The denominator of (37) is given by (35). The numerator is evaluated as

Z ∞ 0 0 0 0 2 d∆t ∆t hBdecay|B (∆t )i 0 |r|2 |q|2|s|2 1 1 = + 2 2 + 2 2 2|p| 2ΓL 2ΓH |r|2 |q|2|s|2 Γ2 − ∆m2 + − 2 2|p|2 (Γ2 + ∆m2)2 |r||q||s| 1 1 + cos (θ − ϕ) 2 − 2 |p| 2ΓL 2ΓH |r||q||s| 2Γ∆m − sin (θ − ϕ) . (38) |p| (Γ2 + ∆m2)2

|q|2|s|2 1 |q|2|s|2 Γ2−∆m2 |q||s| 4Γ∆m 1 + 2 2 2 + 1 − 2 2 2 2 2 + sin (θ − ϕ) 2 2 2 0 |p| |r| Γ |p| |r| (Γ +∆m ) |p||r| (Γ +∆m ) τeﬀ (B → Bdecay) = , |q|2|s|2 1 |q|2|s|2 Γ |q||s| 2∆m 1 + |p|2|r|2 Γ + 1 − |p|2|r|2 Γ2+∆m2 + |p||r| sin (θ − ϕ) Γ2+∆m2

2 2 1 + |A |2 1 + 1 − |A |2 Γ −∆m + 4Im [A ] Γ∆m w Γ2 w (Γ2+∆m2)2 w (Γ2+∆m2)2 = . (39) 2 1 2 Γ ∆m 1 + |Aw| Γ + 1 − |Aw| Γ2+∆m2 + 2Im [Aw] Γ2+∆m2

Here, Aw is the weak value given in (1), which now takes This amplification effect can also be seen in the detec- the form, tion of CP violation. For simplicity, we assume |p| = |q| and ∆m/Γ 1 so that the difference between the de- |q||s| |q||s| nominator of (36) and its counterpart for the CP conju- A = cos(θ − ϕ) + i sin(θ − ϕ), (40) w |p||r| |p||r| gate is small enough to be ignored. Then, the difference between their time distributions is found as upon using |ψi = |B0i, |φi = |B i and Aˆ defined decay 0 ¯0 ¯ in (3) together with (17). P (∆t|B → Bdecay) − P (∆t|B → Bdecay) Note that we have obtained the result (39) without re- p1 − |r|2 ∝ − [sin(θ − ϕ) − sin(θ + ϕ)] e−Γ∆t sin(∆m∆t) sorting to the first-order approximation for ∆m/Γ. In |r| fact, the full calculation is required for us here since p1 − |r|2 ∆m/Γ ' 0.77 for the B meson system we are consider- = 2 cos θ sin ϕe−Γ∆t sin(∆m∆t). (42) ing [21]. For B mesons, Figure 2 shows that this effective |r| lifetime can be 2.6 times larger than the lifetime τ(B0). Returning to the general case, we note that when This result indicates that the parameters |r| and θ in- ∆m/Γ 1 holds, one may expand the effective life- troduced in our postselection can be used to amplify the time (39) in terms of ∆m/Γ to obtain difference in the conditional time distribution.

0 1 1 ∆m 2 τeff (B → Bdecay) = + Im[Aw] + O (∆m/Γ) . E. Process of postselection and related issues Γ Γ Γ (41) In order to perform the above amplification, we need The extension of the effective lifetime is proportional to to give an operation to choose the state |Bdecayi. This the imaginary part of the weak value Aw in the first order will be performed by identifying the type of particles pro- of ∆m/Γ, which can also be confirmed directly from (7). duced at the B meson decay. 8

We thus see that choosing the decay mode |fi characterized by the parameters r and s via the amplitude Af and A¯f amounts to choosing the state |Bdecayi, i.e., the postselection. Our problem now boils down to ﬁnding how to choose the decay mode |fi with desired r and s. To demonstrate that this may in principle be possible, we recall that the Standard Model predicts that photons produced in the decay b → sγ are primarily left- handed [22]. For simplicity, in this paper, we assume that the involved photons are always left-handed, which implies that the photons produced from the decay of anti-b quarks are always right-handed. It then follows that the combinations of the decayed particles are

0 ∗0 0 ∗0 |B i ↔ |K i|γRi, |B¯ i ↔ |K¯ i|γLi, (49) where |γ i and |γ i represent the left-handed photon 0 L R FIG. 2. The ratio between the lifetime of B meson τ(B ) and the right-handed photon, respectively. Also, if we and the effective lifetime (39). The effective lifetime can be 0 can assume, in effect, that this process preserves the CP extended about 2.6 times larger than τ(B ) in this case when symmetry, [CP, Sˆ] = 0, then we have |r| ∼ 0.2 and sin(θ − ϕ) ∼ −1. 0 ∗0 Sˆ|B i = c |K i|γRi + ··· , (50) ¯0 ∗0 To be explicit, let |fi be the state of decayed parti- Sˆ|B i = c |K¯ i|γLi + ··· , (51) cles which we identify from the measurement. The state of the B meson |B0(∆t)i at the time of decay is given with a common constant c, where ··· stands for other by (15), and the transition probability from the state particles generated from the decay. |B0(∆t)i to the state |fi is given by | hf|Sˆ|B0(∆t)i |2 us- This alludes us to consider the decay mode |fi by measuring the states of the corresponding particles separately ing the unitary operator Sˆ corresponding to the S-matrix. in the decay, i.e., we now have Expanding the state |B0(∆t)i in terms of the basis states |B0i and |B¯0i as ∗ |fi = |Kf i|γf i, (52) 0 0 0 |B (∆t)i = a(∆t)|B i + b(∆t)|B¯ i, (43) ∗ where |Kf i and |γf i are the states of the respective par- ∗ we find that the probability becomes ticles determined by the measurement. Expanding |Kf i in terms of the basis states as | hf|Sˆ|B0(∆t)i |2 = |a(∆t) hf|Sˆ|B0i + b(∆t) hf|Sˆ|B¯0i |2, 2 |K∗i = ξ |K∗0i + ξ |K¯ ∗0i, (53) = |a(∆t)Af + b(∆t)A¯f | , (44) f 1 2 where and similarly |γf i as ˆ 0 ¯ ˆ ¯0 Af = hf|S|B i , Af = hf|S|B i , (45) |γf i = η1|γRi + η2|γLi, (54) are the corresponding decay amplitudes from the basis and using (50), (51) together with the projection on the states. Comparing (44) with (30), we observe that if the space of B meson, we obtain condition ∗ ∗ ˆ ∗ ∗ ∗0 ∗ ∗ ¯ ∗0 r Af hKf |hγf |S = ξ1 η1 hK |hγR| + ξ2 η1 hK |hγR| = , (46) B meson s∗ A¯ f ∗ ∗ ∗0 ∗ ∗ ¯ ∗0 ˆ + ξ1 η2 hK |hγL| + ξ2 η2 hK |hγL| S is satisfied, the probability (44) becomes identical to the probability (30) as a function of r and s up to normal- (|B0ihB0| + |B¯0ihB¯0|) ization, ∗ ∗ ¯0 ∗ ∗ 0 = c (ξ2 η2 hB | + ξ1 η1 hB |). (55) ˆ 0 2 0 2 | hf|S|B (∆t)i | ∝ |hBdecay|B (∆t)i| . (47) From (55), we find This relation allows us to find the state hB | corre- decay ∗ ˆ 0 ∗ ∗ sponding to the state hf| Sˆ by applying the projection Af = hKf | hγf |S|B i = c ξ1 η1 , (56) 0 0 ¯0 ¯0 ¯ ∗ ˆ ¯0 ∗ ∗ operator |B i hB | + |B i hB | onto the space of B me- Af = hKf | hγf |S|B i = c ξ2 η2 . (57) son as This implies that the consistency condition (46) is satis- ˆ 0 0 ¯0 ¯0 hBdecay| ∝ hf| S(|B i hB | + |B i hB |), fied if we choose the expanding parameters as 0 0 = Af hB | + A¯f hB¯ | , ξ1η1 r ∗ 0 ∗ ¯0 = . (58) ∝ r hB | + s hB | . (48) ξ2η2 s 9

photonic systems is also available for high-energy photons, a more realistic method of choosing the polarized state (61) may be provided by measurement of angular distribution of decayed particles, which we shall mention in some detail in Section V. To sum up, we have learned that our desired postselection (27) on the B meson may be accomplished by choosing a corresponding decay mode appropriately (see Figure 3). Incidentally, we note that a similar argument applies to the CP conjugate state, |B¯0(∆t)i = a(∆t)|B¯0i + b(∆t)|B0i, (62) reaching essentially the same result (47) with |B0(∆t)i replaced by |B¯0(∆t)i. Note also that K∗ and γ produced by the decay of the B meson evolve over time after the production. While the state of the photon does not change as long as it ﬂies in FIG. 3. After the decay of the B meson, the selection is per- ∗ ∗ ∗ the vacuum, the K state may change during the ﬂight. formed for each of K and γ. For the K , it is measured Therefore, an additional operation is required to specify whether it is |K∗0i given by (59) or |K∗0 i, which is the or- 1 1⊥ the state (59) taking into account this time evolution. thogonal state. Similarly, for the γ, we measure whether it is

|γX i given by (61) or its orthogonal state |γX⊥ i. In the case of low-energy photons (visible light) often used in weak measurements, the beam splitter is commonly used for the post- IV. SIMULATION STUDY selection, while in the case of high-energy photons, a proper measurement of angular distribution of decayed particles may A. Introduction be considered alternatively. We shall now move on to discuss the weak measurement of the B0 − B¯0 mixing using the 0 ∗0 0 0 We thus see that, if we can freely select ξ1, ξ2, η1, and η2, B → K γ → KSπ γ decay channel from for the decay mode |fi such that (58) is fulfilled, we can the experimental point of view. Since our forgo- freely select r and s for the postselected state |Bdecayi, ing observation shows that the postselection after the 0 ∗0 0 0 that is, the postselection can be made freely. B → K γ → KSπ γ decay can in effect be The parameters ξ1, ξ2, η1, η2 are not uniquely deter- replaced by the postselection before the decay, we fo- mined from the condition (58), and one possible choice cus on the latter in the following discussions. Here- ∗ ∗0 is provided by choosing |Kf i = |K1 i, where after, we shall restrict ourselves to the type of flavor tagging carried out by identifying one of the species in the (C = −1) entangled pair of B mesons as adopted ∗0 1 ∗0 ∗0 |K i = √ |K i + |K¯ i . (59) in the Belle II experiment. This requires us to take 1 2 into account negative ∆t as well as positive ones, since

∗0 the decay of one meson in the entangled pair may have This CP eigenstate |K1 i is experimentally distinguished 0 0 already occurred when the tagging of the other meson as it admits the decay mode into KSπ , which is not the is made. To accommodate this extension in the prob- case with its orthogonal state 0 ability density function P (∆t|B → Bdecay) in (36), we consider, in addition to the transition probability 1 0 2 |K∗0 i = √ |K∗0i − |K¯ ∗0i . (60) |hBdecay|B (∆t)i| in (30) for ∆t > 0, the correspond- 1⊥ 2 0 2 0 i∆tHˆ 2 ing one |hB |Bdecay(−∆t)i| = |hB |e |Bdecayi| for ˆ ˆ † ∗0 ∆t < 0. This implies H → H when we go over from Thus, the state |K1 i can be determined experimentally ∆t > 0 to ∆t < 0, and for our case where the Hamilto- without ambiguity (with no component of |K∗0 i), which 1⊥ nian admits the form (17), this amounts to the parameter ∗ ∗0 implies that the postselection |Kf i = |K1 i can be im- replacement (p, q, Γ, ∆Γ) → (q∗, p∗, −Γ, −∆Γ). plemented relatively easily. With this choice, one easily 0 In the probability density function P (∆t|B → Bdecay) finds from (58) that the photon state (54) is now deter- in (36), the effect of CP violation in mixing arises from mined as |γf i = |γX i where the difference between p and q. In the SM, the difference of their absolute scalars 1−|q|/|p| is expected to be small |γX i = r|γRi + s|γLi. (61) O(10−3) [21] in the B0 − B¯0 system. Then, assuming |q|/|p| = 1 for simplicity, we find from (36) and its CP Although it would be unrealistic to assume that the stan- conjugate partner that the probability density functions dard method of using the beam splitter for low-energy reduce to 10

− |∆t| e τ h i P (∆t|B0 → B ) = 1 + (2|r|2 − 1) cos (∆m∆t) − 2|r|p1 − |r|2 sin (θ − ϕ) sin (∆m∆t) , (63) phys decay 2N − |∆t| e τ h i P (∆t|B¯0 → B¯ ) = 1 + (2|r|2 − 1) cos (∆m∆t) − 2|r|p1 − |r|2 sin (θ + ϕ) sin (∆m∆t) , (64) phys decay 2N

where τ = 1/Γ is the B0 lifetime [21] and C. Event simulation 2|r|2 − 1 N = τ 1 + (65) Improvement of sensitivity to the CP violation with 1 + (τ∆m)2 weak measurement is evaluated in the B0 → K∗0γ de- is a normalization factor. The formula (64) is obtained cay as a benchmark by using pseudo-experiments. The ∗0 0 from (63) by putting ϕ → −ϕ, which is a consequence K → KSγ decay mode is selected in addition to fulfill of the interchange p ↔ q. The extension to ∆t < 0 in the consistency condition as discussed in Section III. E. this simplified case is achieved just by putting τ → −τ The number of events is assumed to be 550 × 108 BB¯ according to the note in the previous paragraph, result- in accordance with the expected integrated luminosity of ing in the exponential factor in (63) and (64). We note the Belle II experiment, where BB¯ contains both B0B¯0 that in these formulae the only measurable parameter is and B+B−. The branching fractions of Υ(4S) → B0B¯0 0 ∗0 ϕ which corresponds to 2φ1 in the SM (see Section 12 and B → K γ decays are considered as 0.49 and −5 0 in [21]), where φ1 is one of the interior angles of the uni- 4.2 × 10 , respectively [21]. Since the decay mode of KS 0 + − 0 ∗0 tarity triangle. If an extra phase arises from some new is selected to be KS → π π to reconstruct B → K γ ∗ ∗0 0 0 physics in B → K γ penguin diagram, the measured ϕ events via K → KSπ decay in the Belle experiment, will deviate from the 2φ1. only the same decay mode is selected in this study. The ∗0 0 0 0 + − branching fractions of K → KSπ and KS → π π are 0.17 and 0.69, respectively. B. Belle II experiment The efficiency to identify the flavor of B mesons, B0 0 or B¯ , decaying into ftag is assumed to be 0.136, where The Belle II experiment is currently operated at High the wrong tagging fraction is 0.02 [23]. The reconstruc- Energy Accelerator Research Organization (KEK) in tion and selection efficiency will also be referred to the Japan, where one of the purposes is to measure the CP results in the Belle experiment, which is estimated to be violation in rare decay processes of B mesons. The ex- 18.2% [24]. Finally, the signal yield is expected to be periment will collect 550 × 108 BB¯ pairs at 50 ab−1, 3.2 × 103. that are produced from resonance decays of Υ(4S) at the SuperKEKB. The SuperKEKB is an asymmetric energy electron-positron collider with 7 GeV and 4 GeV of electron and position beam energy, respectively. Υ(4S) is Lorentz-boosted with βγ = 0.28 along the direction of the electron beam. The Belle II detector consists of several sub-detectors surrounding the electron-positron collision point, which are designed to realize precise measurements of particle trajectories, production and decay vertices and particle identification. ∗0 The B meson decaying to the frec final state, K γ in this study, is reconstructed from the decay particles detected by the Belle II detector. The remaining particles in the detector, assumed as the particles from the other B meson decaying to the ftag final state, are used to tag the flavor of the B meson decayed to the frec state, where the tagging time ttag is the instant when the decay of the B meson into ftag occurred. The probability den- FIG. 4. Time distribution of the signal and background, ¯ 0 sity functions of B and B mesons, given by (36) and its Pphys(∆t|B → Bdecay) and Pbkg(∆t), respectively. The CP conjugate partner, describe the time variation of the solid lines are the signal distributions for various values of decays as a function of ∆t which can be determined from the postselection parameter, |r|, and the dashed line is the background. the displacement between frec and ftag vertices along the beam axis, thanks to the asymmetric electron-positron beam energy. Taking all these factors into account, we learn that 11 the generation of simulated events follows the probability density function,

Z ∞ 0 0 h 0 0 0 0 0 i P(∆t|B → Bdecay) = d(∆t ) fphysPphys(∆t |B → Bdecay)Rphys(∆t − ∆t ) + (1 − fphys)Pbkg(∆t )Rbkg(∆t − ∆t ) , −∞ (66)

0 and its CP conjugate partner P(∆t|B¯ → B¯decay) in discriminate the signal from backgrounds because the the corresponding form, both of which are characterized separation of them is more significant as |r| goes to zero. by the same fraction of the signal fphys = 0.66 taken The time resolution functions for the signal Rphys and from [24], together with Pphys and Pbkg which are the dis- background Rbkg are evaluated separately and convo- tributions of the signal equivalent to (63) and the back- luted to each distribution, while the same profile is used ground, respectively. The main background sources are both for Rphys and Rbkg. The time resolution depends the light quark pair production process (e+e− → qq¯ on the position resolution of decay vertices which should with q = u, d, s, c) and the e+e− → BB¯ process be considered to reproduce the timing response in the where the B meson decays into final states different actual experiment. Figure 5 shows the residual of ∆tsim from the signal leading to misidentification. The pro- and ∆ttrue, where ∆tsim is the simulated ∆t where the file of the background distribution is empirically de- position resolution of decay vertices is taken into account, termined from the results in the Belle experiment as based on the Belle experiment [24], and ∆ttrue is the true 0 0 1 −|∆t |/τbkg Pbkg(∆t ) = e , where τbkg is set to value of ∆t without the detector response. τbkg be 0.896 ps. Figure 4 shows the time distribution of To obtain the resolution function of ∆t, Figure 5 is 0 fitted by an empirical function, called ‘double-sided crys- the signal, Pphys(∆t|B → Bdecay) and backgrounds, P (∆t). The postselection can be a powerful tool to tal ball function’. The double-sided crystal ball function bkg consists of a combination of a Gaussian and power law tails, defined as

2  − X e 2 , for − αL ≤ X ≤ αH ,  2 −n  αL h i L  − αL nL e 2 ( − αL − X) , for X < −αL, R(x) = N × nL αL (67) 2 −n  αH h i H  − αH nH e 2 ( − αH + X) , for X > αH , nH αH

with the shorthand, this study, so that dependence of sensitivity to ϕ only on x − µ these parameters can be investigated. The method how X = , (68) to determine these postselection parameters is discussed σ in Section V. where N is the normalization factor, µ and σ are the mean and width of the Gaussian, αH and αL are the positions of transition from the Gaussian to the power D. Signal extraction low tails on the higher or lower sides, and nH and nL are the exponents of the high and low tails, respectively. The resolution functions are well modeled by fitting In order to evaluate the measurement accuracy of the the simulated data under the parameters shown there. CP parameter ϕ, an unbinned maximum likelihood fit is The postselection parameters |r| and θ in (63) and (64) performed for the distributions obtained in the pseudo- represent the mixing ratio and the phase difference be- experiments. The probability density function given tween |B0i and |B¯0i states, respectively. They are the in (66) depends on the fitting parameter ϕ as well as sys- parameters of the postselection that have to be deter- tematic uncertainties coming from the three sources: the mined by event selections in the experiment. Although background estimation, the timing scale and the timing their determination is an important technical subject in resolution, which correspond to fphys, µ and σ appeared the actual experiment, for simplicity we assume to have in (66) with (68). A set of systematic uncertainties is some means to determine them as we wish. The selec- parameterized by the so-called nuisance parameters θnuis tion efficiencies of |r| and θ are both assumed as 100% in that follow the standard Gaussian distribution G(θnuis) 12

world average value of φ1 = 22.2 ± 0.7 degrees [27]. The unbinned maximum likelihood fit is performed to each pseudo-experiment. Figure 6 shows the ∆t distributions when |r| and θ are selected to be 0.5 and 0.0 degree, respectively, which minimize the total uncertainty of measured ϕ. The fitted functions are shown in the plot, which are normalized to the number of events in each distribution. The shift of the peak position between B0 → K∗0γ and B¯0 → K∗0γ shows the effect of the CP violation in the B0−B¯0 mixing with the CP parameter ϕ.

FIG. 5. Distribution of the (∆tsim and ∆ttrue) diﬀerences (residuals). The points are simulated. The red line is the ﬁtted function which is empirically modeled by the double- sided crystal ball function with parameters shown above.

(see Section 40 in [21]). Then, the likelihood function L is defined by the product, L = P(∆t, θnuis , θnuis, θnuis) G(θnuis ) G(θnuis) G(θnuis), fphys µ σ fphys µ σ (69) where these G(θnuis) provide constraints on the system- FIG. 6. Distribution of ∆t for B0 and B¯0 initial states. The atic uncertainties [25]. postselection parameters |r| and θ are selected to be 0.5 and The unbinned maximum likelihood fit is performed si- 0.0 degrees, respectively, which minimize the total uncertainty 0 multaneously for the distributions of B → Bdecay and of measured ϕ. The red and blue lines represent the best fit 0 B¯0 → B¯ . Systematic uncertainties estimated in the probability density functions, P(∆t | B → Bdecay) and decay ¯0 ¯ Belle experiment [24, 26] are used in this study. The P(∆t | B → Bdecay), respectively. dominant systematic uncertainties derive from the background estimation, the timing scale and the resolution of uncertainties in the position of the reconstructed vertices. Figure 7 shows the total uncertainty of measured ϕ Only these uncertainties are considered in the unbinned as a function of |r| and θ which are used as the input maximum likelihood fit. values. This plot indicates that the measurement precision is improved by 20% with |r| = 0.5 compared to the case in which the final state is the CP eigenstate with E. Analysis result |r| = √1 ' 0.7. Note that the precision is virtually in- 2 sensitive to θ. The systematic uncertainty is in the range A pseudo-experiment is carried out and repeated a of 0.2 - 1.1 degrees while the statistical uncertainty varies thousand times to estimate the measurement precision from 3.4 to 9.7 degrees, depending on the postselection of ϕ per particular set of the post selection parameters parameters, indicating the dominance of the statistical |r| and θ, in which events are randomly generated, fol- uncertainty in our analysis. lowing the probability density functions, (66) and its CP The effective lifetime can be prolonged by selecting conjugate partner. These parameters are scanned in the small |r| in the postselection as shown in Figure 2. The range of 0.1 ≤ |r| ≤ 0.9 with the intervals of 0.1 and level of contamination of background peaked at ∆t = 0 similarly in the range of −180 ≤ θ ≤ 180 degrees with tends to diminish as the effective lifetime becomes longer, the intervals of 36 degrees. The initial value of ϕ is set i.e., the uncertainty on background estimation can be re- to 44.4 degrees which corresponds to ϕ = 2φ1 with the duced for smaller |r| up to ∼ 0.5. The uncertainty of ϕ, 13 however, starts to grow for |r| smaller than ∼ 0.5, because 1.5 the term containing sin(θ − ϕ) in (63) and (64) becomes ψ 0 [rad] J/ (nS)K

ϕ negligible then. This eﬀect can be seen in Figure 7. The 1 ρ+π - - ﬁtting method is validated with the number of pseudo- θ η 0 experiments we have performed, in which the results are K 0.5 π +π - consistent with the input value within the statistical un- ρ+ρ- certainty. 0 χ K 0 c1 S 0 −0.5 D* CPh J/ ψ (1S)ρ0 - −1 D* +D* + - 0 K K K S −1.5 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 r

FIG. 8. Distribution of 10 major postselections performed in the previous measurements of CP violation in B decays. The values are estimated from the data in PDG [21], pp.63-64, under the identiﬁcation C = 2|r|2 − 1, S = −2|r|p1 − |r|2 sin(θ − ϕ) obtained from our formulae (63), (64) in Section IV, for parameters r, θ of postselections. The range in the choices of |r| is seen√ to be restricted to a nar- row band around |r| = 0.7(≈ 1/ 2). For a ﬁxed value of ϕ, the distribution along the vertical axis (plotted in the unit of FIG. 7. Total uncertainty of the ϕ measurement (δϕtotal) as radian) amounts to the distribution of θ. a function of |r| and θ used as the input values.

fect in the coefficient of sin(∆m∆t). However, with the conventional method, we cannot separate the absolute amplitude and the phase in the new physics. The sepa- V. CONCLUSION AND DISCUSSIONS ration may become possible with the weak measurement if it is performed under different sets of parameters (r, θ), We have shown that the effect of postselection is clearly elucidating the details of the new physics we have found. seen in the time distribution of the B meson decay, and In high energy physics, due to the short coherence that the postselection has the effect of amplifying both length compelled by the large momentum involved in the effective lifetime and the CP violation. In our sim- the experiment, one normally supposes that approaches ulation study, which assumed that the measurement is which exploit superpositions of states such as the weak performed under the parameters of the Belle II experi- value amplification will not be available. However, our ment, the effective lifetime is found to be prolonged by analysis shows that the coherence time, rather than the 2.6 times, and the accuracy of CP -violating parameters coherence length, is a deciding factor for the amplifica- can also be improved when it is realized in the decay tion to work, as we have demonstrated in our analysis for mode we considered. Although these results may not be the B0 − B¯0 system. Similar positive outcomes will be as great as one might wish to see, we believe that our case expected generally for systems which admit transitions study presented in this paper certainly indicates the po- characterized by the average decay width Γ and the dif- tential use of weak measurement for high energy physics, ference of the energy ∆m. especially with respect to the exploitation of the freedom The success of the weak value amplification hinges on of postselection. Our attempt provides a first systematic the capability of exploiting the freedom of postselection. study for such exploitation, given that the freedom of In high energy experiments, this must be realized at the postselection has not been utilized comprehensively for level of final decay modes which are actually measured. amplifying the prospected signals before. This can be In our case, this is achieved by considering measurements seen clearly in Figure 8, where postselections realized by on the particles K∗ and γ generated from the B decay, previous measurements form only a limited class in the in which our proper choice of the final states satisfies the possible range leaving a vast room of freedom unexplored. consistency condition to reproduce the transition proba- Our case study may lead to more interesting results bility evaluated directly with the postselected states. As in view of the new physics related to the decay mode such, the amount of data to be used for amplification may B → K∗γ. If the effect signalling the new physics be considerably reduced rendering the statistical accu- associated with the right-handed current is sizable to racy poorer. It is thus desirable to look for other decay B → K∗γ, we should be able to measure the finite ef- modes, fulfilling the same consistency condition to im- 14

FIG. 9. Distribution of the polar angle of charged leptons from the photon conversion for |r| = 0.8 (left) and 0.5 (right). The z-axis is normalized to the number of the expected events in the Belle II experiment. The inset shows the same distribution zooming up the range of θ±.

prove the statistics and eventually the accuracy as well. ∼45% by selecting θ− < 2 mrad. Once |r| is identified, We also remark that, in order to improve the sensitiv- the procedure described in [30] then furnishes a possible ity of our weak measurement further, those contributions means to determine the other postselection parameter θ. from the direct and indirect CP violation, which we have Although the principle to determine the photon circular neglected for simplicity in our analysis, should also be polarization has been discussed already long time ago, taken into account. its experimental demonstration has not been done so far. For actual experimental implementation, one has to The Belle II experiment has a good transverse momen- devise a workable procedure to achieve the postselection tum (pT) resolution of 0.3-0.6% for charged particles with with a given set of parameters, |r| and θ, specifying the pT = 0.5-1.0 GeV/c [31], which will result in the angu- superposition of |B0i and |B¯0i. The circular polarization lar resolution better than 1.2%. This may be enough to of the photon is considered as a key element for the post- identify the photon circular polarization with the method selection for B → K∗γ in our study. Identification of the discussed above. photon polarization may be realized by utilizing the fact Furthermore, the technique of the developing machine that angular distribution of charged lepton pairs from learning may also be useful to find the optimal |r| and θ the photon conversion depends on it [28, 29]. This may and thereby ascertain the benefit of the weak measure- allow us to select one of the postselection parameter |r|, ment systematically. although the probability of photon conversion is typically Finally, we wish to argue that, in order to explore a low: approximately 6% at Belle II. In addition to pho- new horizon beyond the conventional means of exper- 0 ∗ + − ton conversion events, contribution from B → K l l imentation of high energy physics, it is important to +1.2 −7 channel with (9.9−1.1) × 10 of the measured branching actively utilize the properties of elementary particles as ratio [21] can be added. Figure 9 shows the distributions genuine quantum states. We hope that the present study of the polar angle θ± of the positive and negative charged encourages those who share the goal in the same direc- leptons created from the photon conversion with respect tion. to photon momentum in the lab frame for |r| = 0.8 (left) and 0.5 (right), assuming that photon energy is 2.6 GeV. The fact that the approximation E±θ±/ml ∼ 1 maxi- mizes sharply the amplitude of photon conversion [30] is used in this calculation, where E± is energy of the posi- ACKNOWLEDGMENTS tive and negative charged leptons, and ml is their mass. Figure 10 shows the probability density functions of |r| and θ±. The dependence of angular distribution on |r| We would like to thank Takanori Hara, Yutaka will allow us to implement the postselection for |r| by se- Shikano, Masataka Iinuma, Satoshi Iso, Kazutaka Sumi- lecting a certain region of the distribution. For example, sawa and Katsuo Tokushuku for useful discussions. This from the events distributed at low θ+ and θ− for small work was supported by JSPS KAKENHI Grant Number and large |r|, respectively, the efficiency is estimated as 20H01906. 15

FIG. 10. Probability density function of |r| and the polar angle of charged leptons from the photon conversion: θ+ (left) and θ− (right).

[1] Y. Aharonov, D. Z. Albert, and L. Vaidman, How the Phys. J. C 74, 3026 (2014), arXiv:1406.6311 [hep-ex]. result of a measurement of a component of the spin of a [13] I. I. Bigi and A. I. Sanda, CP Violation, 2nd ed. (Cam- spin-1/2 particle can turn out to be 100, Phys. Rev. Lett. bridge University Press, 2009). 60, 1351 (1988). [14] J. Brodzicka et al. (Belle Collaboration), Physics achieve- [2] S. Kocsis, B. Braverman, S. Ravets, M. J. Stevens, R. P. ments from the Belle experiment, Prog. Theor. Exp. Mirin, L. K. Shalm, and A. M. Steinberg, Observing the Phys. 2012, 04D001 (2012), arXiv:1212.5342 [hep-ex]. Average Trajectories of Single Photons in a Two-Slit In- [15] A. B. Carter and A. I. Sanda, CP Nonconservation in terferometer, Science 332, 1170 (2011). Cascade Decays of B Mesons, Phys. Rev. Lett. 45, 952 [3] T. Mori and I. Tsutsui, Weak value and the wave-particle (1980). Duality, Quantum Stud.: Math. Found. 2, 371 (2015), [16] A. B. Carter and A. I. Sanda, CP violation in B-meson arXiv:1410.0787 [quant-ph]. decays, Phys. Rev. D 23, 1567 (1981). [4] T. Denkmayr, H. Geppert, S. Sponar, H. Lemmel, [17] I. Bigi and A. Sanda, Notes on the observability of CP A. Matzkin, J. Tollaksen, and Y. Hasegawa, Observation violations in B Decays, Nucl. Phys. B 193, 85 (1981). of a quantum Cheshire Cat in a matter-wave interferom- [18] T. Abe et al., Belle II Technical Design Report, Tech. eter experiment, Nat. Commun. 5, 4492 EP (2014). Rep. (2010) comments: Edited by: Z. Doleˇzaland S. [5] Y. Aharonov and L. Vaidman, Complete description of a Uno, arXiv:1011.0352. quantum system at a given time, J. Phys. A: Math. Gen. [19] K. Akai, K. Furukawa, and H. Koiso (SuperKEKB), Su- 24, 2315 (1991). perKEKB Collider, Nucl. Instrum. Meth. A 907, 188 [6] A. Danan, D. Farfurnik, S. Bar-Ad, and L. Vaidman, (2018), arXiv:1809.01958 [physics.acc-ph]. Asking Photons Where They Have Been, Phys. Rev. Lett. [20] E. Kou et al., The Belle II Physics Book, Prog. Theor. 111, 240402 (2013). Exp. Phys. 2019, 123C01 (2019), arXiv:1808.10567 [hep- [7] L. Vaidman and I. Tsutsui, When Photons Are Lying ex]. about Where They Have Been, Entropy 20, 538 (2018). [21] P. A. Zyla et al. (Particle Data Group), Review of Par- [8] O. Hosten and P. Kwiat, Observation of the Spin Hall ticle Physics, Prog. Theor. Exp. Phys. 2020, 083C01 Effect of Light via Weak Measurements, Science 319, (2020). 787 (2008). [22] B. Grinstein, Y. Grossman, Z. Ligeti, and D. Pirjol, Pho- [9] P. B. Dixon, D. J. Starling, A. N. Jordan, and J. C. ton polarization in B → Xγ in the standard model, Phys. Howell, Ultrasensitive Beam Deflection Measurement via Rev. D 71, 011504 (2005). Interferometric Weak Value Amplification, Phys. Rev. [23] H. Kakuno et al. (Belle Collaboration), Neutral B flavor Lett. 102, 173601 (2009). tagging for the measurement of mixing-induced CP vio- [10] J. Huang, Y. Li, C. Fang, H. Li, and G. Zeng, Toward lation at Belle, Nucl. Instrum. Meth. A 533, 516 (2004), ultrahigh sensitivity in weak-value amplification, Phys. arXiv:hep-ex/0403022. Rev. A 100, 012109 (2019). [24] T. Horiguchi et al. (Belle Collaboration), Evidence for [11] I. Shomroni, O. Bechler, S. Rosenblum, and B. Dayan, Isospin Violation and Measurement of CP Asymmetries Demonstration of Weak Measurement Based on Atomic in B → K∗(892)γ, Phys. Rev. Lett. 119, 191802 (2017), Spontaneous Emission, Phys. Rev. Lett. 111, 023604 arXiv:1707.00394 [hep-ex]. (2013). [25] G. Cowan, K. Cranmer, E. Gross, and O. Vitells, Asymp- [12] A. J. Bevan et al., The Physics of the B Factories, Eur. totic formulae for likelihood-based tests of new physics, 16

Eur. Phys. J. C 71, 1554 (2011), [Erratum: Eur. Phys. (1965-1970) 40, 13 (1964). J. C 73, 2501 (2013)], arXiv:1007.1727 [physics.data-an]. [29] H. Olsen and L. C. Maximon, Photon and electron polar- [26] Y. Ushiroda et al. (Belle Collaboration), Time-dependent ization in high-energy bremsstrahlung and pair produc- 0 0 0 CP asymmetries in B → KS π γ transitions, Phys. Rev. tion with screening, Phys. Rev. 114, 887 (1959). D 74, 111104 (2006), arXiv:hep-ex/0608017. [30] F. Bishara and D. J. Robinson, Probing the photon po- [27] Y. S. Amhis et al. (HFLAV), Averages of b-hadron, c- larization in B → K∗γ with conversion, J. High Energ. hadron, and τ-lepton properties as of 2018, Tech. Rep. Phys. 09, 013, arXiv:1505.00376 [hep-ph]. (2019) arXiv:1909.12524 [hep-ex]. [31] I. Adachi et al. (sBelle Design Group), sBelle Design [28] H. Kolbenstvedt and H. Olsen, Circular photon polariza- Study Report, Tech. Rep. (2008) arXiv:0810.4084 [hep- tion detection by pair production, Il Nuovo Cimento A ex].