Home Search Collections Journals About Contact us My IOPscience Quantum non-barking dogs This content has been downloaded from IOPscience. Please scroll down to see the full text. 2014 New J. Phys. 16 063026 (http://iopscience.iop.org/1367-2630/16/6/063026) View the table of contents for this issue, or go to the journal homepage for more Download details: IP Address: 149.169.148.202 This content was downloaded on 29/08/2014 at 17:16 Please note that terms and conditions apply. Quantum non-barking dogs Sara Imari Walker1,2,3, Paul C W Davies1, Prasant Samantray4 and Yakir Aharonov5,6 1 Beyond Center for Fundamental Concepts in Science, Arizona State University, Tempe AZ 85287-1504, USA 2 School of Earth and Space Exploration, Arizona State University, Tempe AZ 85287-1504, USA 3 Blue Marble Space Institute of Science, Seattle WA, USA 4 International Centre for Theoretical Sciences, IISc Campus, Bengaluru 560012, India 5 Department of Physics, Computational Science, and Engineering, Schmid College of Science, Chapman University, Orange CA 92866, USA 6 School of Physics and Astronomy, Tel Aviv University, Tel Aviv, Israel E-mail: sara.i.walker@asu.edu and paul.davies@asu.edu Received 29 January 2014, revised 27 April 2014 Accepted for publication 13 May 2014 Published 13 June 2014 New Journal of Physics 16 (2014) 063026 doi:10.1088/1367-2630/16/6/063026 Abstract Quantum weak measurements with states both pre- and post-selected offer a window into a hitherto neglected sector of quantum mechanics. A class of such systems involves time dependent evolution with transitions possible. In this paper we explore two very simple systems in this class. The first is a toy model representing the decay of an excited atom. The second is the tunneling of a particle through a barrier. The post-selection criteria are chosen as follows: at the final time, the atom remains in its initial excited state for the first example and the particle remains behind the barrier for the second. We then ask what weak values are predicted in the physical environment of the atom (to which no net energy has been transferred) and in the region beyond the barrier (to which the particle has not tunneled). Thus, just as the dog that didnʼt bark in Arthur Conan Doyleʼs story Silver Blaze gave Sherlock Holmes meaningful information about the dogʼs non-canine environment, here we probe whether the particle that has not decayed or has not tunneled can provide measurable information about physical changes in the environment. Previous work suggests that very large weak values might arise in these regions for long durations between pre- and Content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI. New Journal of Physics 16 (2014) 063026 1367-2630/14/063026+14$33.00 © 2014 IOP Publishing Ltd and Deutsche Physikalische Gesellschaft S Imari Walker et al New J. Phys. 16 (2014) 063026 post-selection times. Our calculations reveal some distinct differences between the two model systems. Keywords: quantum mechanics, weak measurement, post-selection 1. Introduction Given that quantum mechanics is several decades old, it is remarkable that a significant sector of the theory lay unexplored until recently. To glimpse this sector it is merely necessary to consider the expectation value of an observable A in a quantum state Ψ , assumed for the moment to be a stationary state, and note that by inserting a complete set of states it may be decomposed as follows: ⎡ ϕ ∣ A∣Ψ i Ψ ∣ A∣ Ψ = ∑ ∣ Ψ ∣ ϕi ∣ ⎢ ⎢⎣ ϕ ∣ Ψ i i 2 ⎤ ⎥, ⎥⎦ (1) where ϕi are eigenstates of a different observable B (i.e., B does not commute with A). The first term of the summand is recognized as the probability that, on measurement of B, the system will be found in state ϕi . The second term (in the square braces) is called the weak value of A (one value for each specific ϕi ). For an individual system the weak value has little meaning, but it does contain non-trivial statistical information about the world when one considers a large ensemble of identical systems, each prepared in state Ψ . Under those circumstances, the weak value can be shown to be precisely the mean value of an ensemble of measurements of B in which the coupling between the system and the measurement device is made arbitrarily weak (and hence non-disturbing) [1, 2]. Note that weak values are not eigenvalues. Rather, they are statistical averages of weak measurements. They may take values outside the spectrum of eigenvalues [3]. They may not even be real numbers; the real and imaginary parts have separate physical interpretation [4]. Weak values can and have been measured, and are the subject of considerable theoretical and experimental interest (see e.g., [1, 5–12]). The subject of weak values becomes of greater interest when combined with post-selection of states. That is possible because, for a large enough ensemble of identical systems in identically prepared initial states (pre-selection), there will always be a sub-ensemble of systems which are also found, on subsequent measurement of B, to be in any given eigenstate state ϕi . A generic expression for weak values may be written, schematically, wA = out A in out in , (2) for observable A, with the system prepared in state ∣ in , and post-selected for state ∣ out . In this paper we consider weak values with post-selection of states for systems that are intrinsically time-dependent; that is, they evolve unitarily in time away from a stationary state, so that the time dependence is not merely the result of measurement. Under these circumstances the weak value is given by 2 S Imari Walker et al New J. Phys. 16 (2014) 063026 wA ( t ) = out ∣ U † ( t − tout ) AU ( t − tin ) ∣ in out ∣ U † ( t − tout ) U ( t − tin ) ∣ in , (3) where U(t) is the unitary evolution operator for the system. We consider two examples in this paper: the decay of an excited atom, and the tunneling of a particle through a barrier. The postselection criteria are chosen as follows: at the final time, the atom remains in its initial excited state for the first example and the particle remains behind the barrier for the second. We then ask what weak values are predicted in the physical environment of the atom (to which no net energy has been transferred) and in the region beyond the barrier (to which the particle has not tunneled). Thus, just as the dog that didnʼt bark in Arthur Conan Doyleʼs story Silver Blaze gave Sherlock Holmes meaningful information about the dogʼs non-canine environment [13], here we demonstrate that, for each respective case, the particle that has not decayed or the particle that has not tunneled provide measurable information about physical changes in the environment (readers of a later generation may prefer a Rolling Stoneʼs analogy: ‘I hear the telephone that hasnʼt rung.’). 2. Decay of an excited atom Consider an atom prepared in an excited state at initial time t = ti . Over time it will decay. The decay may be described using first order perturbation theory, which predicts that the expectation value of the projection operator onto the excited state will have the well-known time dependence e−γ ( t − ti ). In the context of weak measurements, it is then possible to ask the following question. Suppose the atom is inspected at some time t f , and found to have definitely decayed. What result would be obtained for a weak measurement of the projection operator made at time t, in the interval ti < t < t f , with the state of the atom at t f post-selected to be decayed (understood in the context of measurements averaged over large ensemble of identical systems)? This problem was studied by one of us (PD) using an exactly solvable model of an initially excited two-level reference atom coupled equally to all members of a large bath of similar two-level atoms, all prepared initially in their ground states [9]. The weak value of the projection operator onto the excited state of the reference atom as a function of time could then be evaluated for various choices of post-selection for the bath atom states. In the special case that the post-selected bath state corresponds to an excited atom with energy levels that coincide with the reference atom, the result was found to be real, and simple: ⎡ 1 − e−γ ( t f − t ) ⎤ ⎥, w = e−γ ( t − ti ) ⎢ ⎢⎣ 1 − e−γ ( t f − ti ) ⎥⎦ (4) which reduces to the standard exponential decay law when t f → ∞. The same model can also be used to solve the complementary question of post-selecting the atom to have definitely not decayed at time t f , and that is the case we wish to address here. The weak value of the projection operator onto the excited state of the reference atom at time t in that case is w = 1. However, our interest here lies with the bath atoms. It might be supposed that if the reference atom, henceforth labeled 0, is both pre-and post-selected to be in its excited state, then weak measurements of the state of the bath in the interval ti < t < t f would inevitably 3 S Imari Walker et al New J. Phys. 16 (2014) 063026 yield values corresponding to the ground states of every bath atom, since no energy will have flowed from the excited atom 0 to the bath at the conclusion of the interval ⎡⎣ ti, t f ⎤⎦. Intriguingly, this supposition is incorrect. To demonstrate this, we calculate the weak values of projection operators onto the excited states of the bath atoms, in the case that atom 0 is both pre- and post-selected to be in its excited state. For simplicity, consider a population of N bath atoms with upper energy levels equispaced and given by En − E0 = nE − N ⩽ n ⩽ N, (5) i.e., the excited states are distributed symmetrically about the excited state of atom 0. Let the weak value for the projection operator onto the excited state of atom n be denoted wn, and the bra vector for the initial state of the total system be denoted as ( 1, 0, 0, 0, …), the first entry corresponding to atom 0 in its excited state and the remaining entries to the n bath atoms in their ground states. The projection operator Pn onto the excited state of atom n will then be, in this notation, a square matrix with all elements 0 except the entry for row n, column n, which will be 1. The Schrödinger equation for this system is a set of coupled differential equations a˙ 0 = − i ∑Han e−inΔEt n a˙n = − iHa 0 einΔEt , (6) where an is the probability amplitude that the atom labeled by n is in the excited state. We set  = 1 for convenience, and choose H to be real for simplicity. In the limit that 2 N → ∞ , ΔE → 0 , H → 0, and HΔEπ → γ (where γ is defined as the decay constant), the above set of equations can be solved exactly using Laplace transforms. The evolution operator (which in this case is a 2N + 1 × 2N + 1 matrix) can be written as: U00 ( t ) = e−γ t −iE0 t Un0 ( t ) = He−iEn t e−γ t −iEn t − 1 γ − inΔE U0n ( t ) = He−iEn t e−γ t −iEn t − 1 . γ + inΔE (7) The elements Unm are not required for what follows. It may be readily verified that the above operator satisfies the unitarity constraint UU † = 1 (for the elements given), and the evolution condition U t f − t U ( t − ti ) = U t f − ti . The weak values of interest are given by ( ) ( ) T [ 1, 0, 0, 0 ,..] U ( t f − t ) Pn U ( t − ti ) [ 1, 0, 0, 0 ..] wn = , (8) T [ 1, 0, 0, 0 ,..] U ( t f − ti ) [ 1, 0, 0, 0 ..] T where [ 1, 0, 0, 0 ,..] is the transpose of the column vector [ 1, 0, 0, 0 ,..], and use has been ( ) ( ) made of the relation U † t − t f = U t f − t . The matrix multiplications are straightforward, 4 S Imari Walker et al New J. Phys. 16 (2014) 063026 and using equation (7) we find, in the aforementioned limits, ⎛ ⎞ γ t −t H2 −γ t − t −inΔE ( t f − ti ) wn = ⎜ 2 ⎟ e ( f i) e ( f i) + e ⎝ γ + n 2ΔE ⎠ ( ) − e−γ ( t − ti )− inΔE( t f − t ) − e−γ ( t f − t )− inΔE( t − ti ) . (9) Equation (9) gives the value that the nth atom, weakly measured, is in an excited state, subject to our specification of pre- and post-selected states. Note that w vanishes (by construction) at t f and ti . Inspection of equation (9) reveals that the weak values of individual bath atoms are generally non-zero in the interval ⎡⎣ t , t ⎤⎦. However, it is readily verified, in the i f limit ΔE → 0 where the sums may be performed explicitly, that ∑wn = 0, (10) n as it must, by unitarity. Thus, on average, the bath atoms are undisturbed during the interval ⎡⎣ t , t ⎤⎦. This can occur only if some w are negative. The fact that an operator (the projection n i f operator onto a bath atomʼs excited state) can be negative seems counter-intuitive to those used to thinking of quantum states in terms of eigenvalues, but it is well-established that such ‘weird’ weak values are commonplace for various post-selections [1, 3]. Another distinctive property of wn is that it increases exponentially with time. This is seen most strikingly in the case that we again choose n = 0, i.e., the bath atoms have the same energy levels as the reference atom, when the expressions are all real: w0 = H2 ( γ 2 1 + e γ ( t f − ti ) − e γ ( t f − t ) − e γ ( t − ti ) ) , (11) The interpretation of this result is as follows. The probability of finding the chosen postselected state—atom 0 undecayed—becomes exponentially small at times in excess of the normal half-life of the atom. That is, an exponentially smaller fraction of N in the ensemble of identical systems will be found on measurement at time t f to be in the state with atom 0 still ( ) excited. Thus U t f − ti ∣i , being the initial state unitarily evolved to the final time t f , will be almost orthogonal to the post-selected final state ∣f : ( ) ( ) f ∣ U t f − ti ∣i = f ∣ U † t − t f U ( t − ti ) ∣ i ≪ 1, (12) Inspection of equation (3) then indicates that w will be very large (positive or negative). The longer the interval t f − ti the smaller equation (12) becomes, and the larger the values of wn become, exponentially. Nevertheless, because wn can be both positive and negative according to the value of n, equation (10) still applies. What can we conclude from this calculation? The result demonstrates that if an excited atom is found after a period of time to have not decayed, this does not mean the electromagnetic field in the vicinity of the atom is undisturbed. The use of weak measurements can reveal activity in the field. This activity will average to zero, but individual weak values of, say, the energy density of the electromagnetic field, will be non-zero and will in fact grow exponentially 5 S Imari Walker et al New J. Phys. 16 (2014) 063026 large (positive and negative) as the post-selection time becomes much longer than the half-life of the excited state. Thus, the atom that doesnʼt decay gives measurable information about physical changes in its environment. We now turn to a second example that demonstrates the same phenomenon, but with some importantly different features. 3. Particle tunneling through a barrier As a second example, let us consider the quantum tunneling of a particle trapped inside a potential well. We model the system by the potential:  2κ δ ( x ); x > − 2L V ( x) = μ = ∞ ; x ⩽ − 2L , (13) where κ > 0, i.e., we model the system with a potential barrier centered at x = 0. The problem is set up as follows. A particle, modeled by a wave packet, is moving to the right and is prepared such that at time t = 0 it is positioned at the center of the well at x = −L . Eventually it encounters the delta function potential and is partially reflected and partially transmitted. The reflected component moves leftward until it encounters the infinite wall at x = −2L , from which it is totally reflected, to repeat its journey in the direction of the delta function potential, from which a second reflection and transmission process occurs. This back and forth motion is repeated a large number of times. Each reflection at the delta function diminishes the amplitude of the reflected component and creates one additional transmitted packet, also with diminished amplitude. If we wait for a sufficiently long time, we will find that the wave packet has almost completely tunneled through the potential. As with the previous case of the atom that did not decay, utilizing weak-measurements permits us to probe interesting features of the system if the particle, at some later time T, is found to not have tunneled through the barrier. We are interested in weak values at an intermediate time 0 < t < T in the region outside the well, subject to the post-selection that the particle is still confined to the well at time T. For simplicity, we model the trapped particle by the Gaussian wave packet 2 Φ ( x, t ) = b 2 b + i t μ e ( x − x 0 − vt ) ik 0 x − x 0 − vt − 2 2 b 2 + imt ( ) ( ), (14) where μ is the particle mass, and x0 is the position of maximal amplitude Φ ( x, t ) 2 at t = 0. To avoid the complications associated with spreading of the wave packet, we assume the following conditions: (i) the mass μ is very large such that b 2 ≫ μt in the range of time that we are interested in, (ii) the initial width b ≪ L , and k (iii) the dominant wave number k 0 is very large such that v = μ0 is finite. 6 S Imari Walker et al New J. Phys. 16 (2014) 063026 With the above conditions, the wave packet can be approximated as 2 Φ ( x, t ) = e ( ) ik 0 x − x 0 − vt − 2 ( x − x 0 − vt ) 2b 2 , (15) which suffices for the purposes of our calculation. Pre-selecting our wave-packet to be located at x = −L at t = 0 moving to the right yields: 2 Φpre ( x, 0) = e ik 0( x + L ) − (x+L) 2b 2 . (16) Likewise, we post-select the state at a later time T to be again located at x = −L and moving to the left: 2 Φpost ( x, T ) = e −ik 0( x + L ) − (x+L) 2b 2 . (17) We have chosen the post-selected state to be the same as the pre-selected state, except for the fact that the post-selected wave is traveling in the opposite direction. These choices of preand post-elections are chosen for convenience without any loss of generality. In order to compute the weak value, one requires an exact time-dependent solution for the potential in equation (13), which for an incident Gaussian wave packet is intractable. However, one can derive the reflection (ρ) and transmission (τ) coefficients for the time-independent Schrödinger equation in the presence of a delta function barrier alone (with no reflecting wall) [14–17], which yields: ρ= − iκ , k 0 + iκ τ= k0 . k 0 + iκ (18) Therefore, to approximate the time-dependent solution without introducing much error, it can be safely assumed that every time a narrow Gaussian wave packet is incident upon a delta barrier, the reflected component (assumed to be a Gaussian) has its amplitude reduced by a factor ρ. The rest of the wave packet is transmitted with a factor τ times the amplitude of the incident wave packet. To add the effects of an infinite wall to the existing delta barrier is non-trivial. For our purpose, it suffices to assume that the wall is at a large distance from the barrier, i.e., Lb ≫ 1 and acts as a mirror. Also b 2 ≫ mt , and the typical time scale will be given by t ∼ L v from dimensional grounds. However, since mv = k 0, we require that bk 0 ≫ L ≫ 1. b (19) With these approximations, we can construct the forward time evolution of the preselected wave packet (16) as 2 2 ⎛ ⎞ ( x + ( 4N + 1) L− vt ) vt − ( x − ( 4N − 3) L + vt ) 4 3 ik x N L − − − + ik 0 ( x + ( 4N + 1) L − vt ) − ( ) 0 ( ) 2 2 ⎟⎟ , 2b 2 2b 2 Φpre ( x, t ) = ρ ⎜⎜ e −e ⎝ ⎠ N 7 (20) S Imari Walker et al New J. Phys. 16 (2014) 063026 for the Nth reflected component in the region −2L < x < 0, and 2 N Φpre ( x, t ) = ∑τρ n−1 ×e ( ) ik 0 x + ( 4n − 3) L − vt − 2 ( x + ( 4n − 3) L− vt ) 2b , 2 (21) n=1 for the transmitted component(s) in the region 0 < x < ∞, where N is the total number of interactions with the delta barrier for the forward evolving wave-packet. As required, the solution satisfies the boundary condition at the wall Φ ( −2L , t ) = 0. It can be seen that with increasing time, the amplitude of the oscillating wave packet decreases by successive powers of the factor ρ. A succession of N such interactions with the barrier creates a train of N transmitted wave packets in the region x > 0 as intuitively expected. Following similar reasoning, the backward time evolution of the post-selected wave packet can be expressed as 2 ⎛ ⎛ v ( t − T ) ⎞ ( x + ( 4 M + 1) L + v ( t − T ) ) ⎟⎟ − −ik 0 ⎜⎜ x + ( 4 M + 1) L + 2 2 ⎠ 2b Φpost ( x, t ) = ρ ⎜ e ⎝ ⎜ ⎝ 2 ⎛ v ( t − T ) ⎞ ( x − ( 4 M − 3) L − v ( t − T ) ) ⎞ ⎟⎟ − ik 0 ⎜⎜ x − ( 4 M − 3) L − 2 ⎠ ⎟, 2b 2 −e ⎝ ⎟ ⎠ M (22) for the Mth reflected component in the region −2L < x < 0, and M Φpost ( x, t ) = ∑ τρ m− 1 m=1 ⎛ v ( t − T ) ⎞ x + ( 4 m − 3) L + v ( t − T ) ⎟− −ik 0 ⎜⎜ x + ( 4 m − 3) L + 2 ⎟⎠ 2b 2 ⎝ ( ×e 2 ) , (23) for the transmitted component(s) in the region 0 < x < ∞, where M denotes the total number of interactions with the delta barrier for the backward evolving wave-packet, and T = tout − tin. Recall that weak values at time t for time-dependent systems are computed at a time t by evolving the pre-selected states forward in time from tin to t and evolving post-selected states backward in time from tout to t (see equation (3)). The weak value at position xk may be calculated using the projection operator P = ∣ x − xk x − xk ∣ , (24) where xk > 0 corresponds to a projection taken outside of the well, and xk < 0 to a projection taken inside the well. Using this projection operator, the weak value becomes: * ( xk , t ) Φpre ( xk , t ) Φpost . (25) w ( xk , t ) = * ( x, t ) Φpre ( x, t ) dx ∫ Φpost It is easy to see that integrating this expression over all space yields, ∫ wdx k = 1, (26) 8 S Imari Walker et al New J. Phys. 16 (2014) 063026 Figure 1. Space-time diagram showing all sweet spots outside the well for post-selected times T = 6L v (point A only), and T = 14L v (points A, B, C, and D). The dotted line follows the space–time trajectory of the reflected and transmitted components of the preselected wave packet, while the dashed line follows that of the post-selected packet. The impenetrable wall is located at x = −2L and the delta function barrier at x = 0. which provides a useful check for the accuracy of the weak value calculations in the discussion below. 3.1. The general solution Substantial weak values will arise at space-time regions where the overlap between the forward and backward evolving wave functions is significant. Inspection of equations (21) and (23) reveals that the transmitted components of these forward and backward evolving wave-packets overlap at certain specific points, or ‘sweet spots’ in the region x > 0. Sweet spots in the spacetime diagram are denoted by their coordinates (n, m), with n and m respectively tracking the number of interactions of the pre- and post-selected transmitted wave packets with the delta barrier at the time of transmission (i.e., the contributions to the sums in equations (21) and (23)). Figure 1 shows the space-time diagram of these sweet spots for the specific case of post-selected times T = 6L v and T = 14L v . Similar space-time diagrams may be drawn for any postselected time T satisfying the condition T = ( 4i + 2) L v, (27) where i = 1, 2, 3, 4… is an integer. The resulting spacetime diagram is symmetric, and the calculations simplify considerably without compromising important qualitative features. We therefore use post-selected times consistent with equation (27) throughout the remainder of this paper. From equation (25), the general solution for the weak value at any sweet spot ( n, m ) with post-selected time consistent with equation (27) is: 9 S Imari Walker et al New J. Phys. 16 (2014) 063026 ⎛ κ ⎞n + N + m − M k ( k 0 + 2iκ ) ⎜ ⎟ ⎝ ik 0 − κ ⎠ 2 0 w ( n, m, N , M ) = ⎛ ⎛ κ ⎞ 2N ⎞ b 2k 02 ⎜ π bκ ⎜ − k 0 + e ( k 0 + 2iκ ) − 2iκ ⎜ ⎟ ⎟⎟ − ik κ ⎠ ⎠ ⎝ 0 ⎝ ( 2 ) (b k +L ( −2i +( 2+2i) n−( 2+2i) N −( 2−2i) m+( 2−2i) M )+ix) (b k +i (L (−2+(2+2i)n−(2+2i)N +(2−2i)m−(2−2i)M )+x)) 2 0 2 0 ×e . b2 (28) The above expression is valid for any N < M , subject to the constraints n ⩽ N and m ⩽ M . Thus, equation (28) is valid for any sweet spot in the lower half of the space–time diagram where the number of forward reflections of the pre-selected packet is less than the number of backward reflections of the post-selected packet (by symmetry an analogous expression applies for the upper half plane with ( n, N ) ↔ ( m, M ) and M < N ). We will first consider this result for the specific post-selection times T = 6L v and T = 14L v before returning to consider the behavior of the weak value for arbitrary post-selection time T. 3.1.1. Post-selection at T ¼ 6L=v . In the atomic model discussed in section 2, the condition in equation (10) ensured that the weak values in the set of bath atoms averaged to zero, consistent with the reference atom being in its excited state at the post-selected time. Similarly in the present model, we expect the weak values outside the well, though individually non-zero (and potentially very large), also to average to zero. This is indeed the case, as may be shown explicitly in the simplest case where the wave-packet interacts only once with the delta function barrier before post-selection at T = 6L v . For this case, N = M = 1 and the weak measurement is taken at t = 3L v (corresponding to point A in figure 1). Using equation (28) the weak value is 2 ( b k +i ( x−2L) ) 2 k 02 e w ( x, N = 1, T = 6L v) = ( π b k 02 + 0 b2 ( 2 2 ) ) , x > 0. (29) − 1 + eb k 0 κ 2 This solution oscillates rapidly between positive and negative values as a function of x for both the real and imaginary parts, with decreasing magnitude as one moves away from the central sweet spot location at x = 2L as shown in figure 2. Using equations (20), (21), (22) and (23) we obtain: 2 ( b k −i ( 2L+x) ) 2 0 w ( x, N = 1, T = 6L v) = − e b2 ( (−1+e π b k 02 + ( 2ik 0 ( 2L + x ) 2 2 2 )κ ) ) 2 , x < 0, (30) − 1 + eb k 0 κ 2 for the reflected component inside the potential well. It is readily verified that ∞ ∞ ∫ w ( x )inside dx + ∫ w ( x )outside dx = 1, as expected from equation (26). −2L 0 3.1.2. Post-selection at T ¼ 14L=v . A more interesting case is presented for larger postselected times T, where the number of interactions with the delta barrier is N > 1. Consider for 10 S Imari Walker et al New J. Phys. 16 (2014) 063026 Figure 2. Shown are the (A) real (blue) and (B) imaginary (pink) components of the weak value outside the well for post-selection at T = 6L v . Constants are set as b = 1; μ = 1 000 ; κ = 1 000; k 0 = 5 000; and L = 100. example, the case for post-selected time T = 14L v shown in figure 1. equation (28) may be used to calculate the weak value at the specific space-time points B (n = m = 2) and D (n = m = 1), where t = 7L v and N = M = 2. The weak values are: w( x, n = 2, T = 14L v) B = 2 ( b k −2iL+ix) 2 0 e − b2 k 02 κ 2 ( π b k 04 + 2ik 03 κ − 2k 02 κ 2 − ( 2 2 ) ) (31) − 1 + eb k 0 κ 4 and w( x, n = 1, T = 14L v)D = 2 ( b k −6iL+ix) 2 0 e b2 ( k 02 ( k 0 + iκ ) π b k 04 + 2ik 03 κ − 2k 02 κ 2 − ( 2 2 2 ) ) (32) − 1 + eb k 0 κ 4 for points B and D, respectively. These expressions are plotted in figure 3 and show the characteristic oscillatory behavior. Comparing figure 3(A) and (B) demonstrates that the weak value increases as one moves away from the delta function barrier (as we discuss in the subsequent section, the ratio of equations (31) and (32) taken at the peaks at x = 2L and x = 6L , respectively, is consistent with equation (34)). 3.1.3. Behavior of the weak value for arbitrary post-selection time T. To gain insight into the behavior of the weak value as a function of the time of post-selection, we evaluate equation (28) for the interior set of points where n = N and m = M (i.e., the time series of sweet spots centered at x = 2L ). Thus, the results of section 3.1.1 represent a special case of those presented here for N = M = 1. For general N and M: 11 S Imari Walker et al New J. Phys. 16 (2014) 063026 Figure 3. Shown are the weak values for post-selection at T = 14L v calculated for (A) N = M = 2, centered at x = 2L corresponding to point B in figure 1 and (B) N = M = 1, centered at x = 6L , corresponding to point D in figure 1. Constants are set as b=1; μ = 1 000 ; κ = 1 000; k 0 = 5 000; and L = 100. w (N, M, N, M ) = ⎛ κ ⎞ 2N k ( k 0 + 2iκ ) ⎜ ⎟ ⎝ ik 0 − κ ⎠ 2 ( b k −2iL+ix) 2 0 e b 2 0 2 ⎛ ⎛ κ ⎞ 2N ⎞ 2⎜ b 2k 02 π bκ ⎜ − k 0 + e ( k 0 + 2iκ ) − 2iκ ⎜ ⎟ ⎟⎟ − ik κ ⎝ ⎠ ⎠ 0 ⎝ ( . (33) ) This solution is independent of number of reflections of the backward evolving wave packet M, and thus independent of the post-selection time T. For any arbitrarily chosen postselected time (consistent with equation (27)), we arrive at the curious result that the weak value at fixed space-time point x and t does not depend on how far into the future we choose the postselected state. Thus, for example, a weak measurement centered at point A in figure 1 (with N = 1) taken at t = 3L v0 will yield a weak value given by equation (29) for any post-selected time T consistent with equation (27). Likewise due to the symmetry of the problem, the same weak value will be observed for a measurement again centered at x = 2L with t = T − 3L v0 (i.e., with M = 1) for any T. Therefore, the weak value can only increase as we move away from the delta-function barrier (e.g. for fixed post-selection time T). For fixed N and M as one moves away from the delta barrier, comparing the weak value at any two neighboring sweet spots yields: w ( n, m, N , M ) w ( n − 1, m − 1, N , M ) = κ2 ( ik 0 + κ )2 , (34) for 1 < n ⩽ N and 1 < m ⩽ M (this of course is not unexpected since on each successive interaction with the delta barrier the transmitted component of the wave packet will diminish in amplitude by a factor τρ). Thus, along the central axis N = M the largest weak value is located at n = m = 1 with an amplification factor ( N−1 κ2 2 ( ik 0 + κ ) ) x = 2L . 12 over the interior most point at centered at S Imari Walker et al New J. Phys. 16 (2014) 063026 4. Discussion When combined with post-selection of states, quantum weak measurements can expose a hidden, measurable, and hitherto neglected, sector of (standard) quantum mechanics. In this paper we have explored one aspect of this sector, namely, time-dependent systems that undergo a transition. We discussed two idealized models: the decay of an excited atom into a large bath of unexcited atoms, and the tunneling of a particle, represented by a Gaussian wave packet, through a thin barrier. The system is prepared at time ti with the reference atom excited and the bath atoms in their ground states, and the packet concentrated behind the barrier, respectively. We chose the final states to correspond to a situation where the system had not undergone the transition, that is, at time t f the atom is determined to be excited in the first example, and the particle is found to be behind the barrier in the second example. It might be supposed that since the system at the end of the experiment is in essentially the same state as it was at the beginning, nothing much of interest could be said about the external region (bath atoms and region outside the barrier, respectively). This is, however, incorrect. Weak measurements in the interval ⎡⎣ ti, t f ⎤⎦ can uncover cryptic activity there. The systems considered have a natural half-life for the transition to occur. The probability that it has failed to occur by time t f falls sharply when t f − ti exceeds this half-life. Nevertheless, given a sufficiently large ensemble of identically prepared systems, there will always be a sub-ensemble that satisfies both the pre- and post-selection criteria. The strength of the cryptic activity in the environs rises with time, i.e., the lower the probability of finding a given system in the post-selected state, the bigger the disturbance revealed by the weak measurements. An important feature of the weak values concerned is that they must average to zero, a condition that follows from unitarity. Thus we find weak values of opposite sign that do indeed sum to zero. Although the net disturbance in the external region is zero, individual weak measurements can be very large, and grow larger with time. The negative weak values found in these examples are strange. We consider projection operators of the bath atoms onto their excited states. Such projection operators always have eigenvalues ⩾ 0. Yet weak values routinely lie outside the spectrum of eigenvalues: many other examples may be found in the literature (see e.g., [1–4]). Similar comments apply to our second example. An unexpected difference emerged between our two examples. In the case of the excited atom, the bath atoms displayed weak values that rose exponentially with time: the longer the atom was detained in its excited state, the bigger the weak values grew among the bath atoms. Their exponential rise mirrors the exponential nature of the decay (or the exponentially small probability of finding members of the ensemble still un-decayed at late times). In the tunneling example, however, there was steady growth, but no exponential surge in weak values. The probability of finding the particle trapped behind the barrier decreases as a power-law in the number of reflections N. One critical difference between the two models, which may explain the results, is that, in the one-dimensional tunneling case, there is a finite probability that the particle will tunnel through the barrier and then tunnel back again. (A Gaussian wave packet moving to the right will still contain some plane-wave components representing left-moving particles.) This probability will be independent of how far the tunneled particle has traveled to the right. On the other hand, for an atom decaying into a large bath, the probability that decay will be followed by a transfer of energy back from the bath to the atom approaches zero in the limit of an infinite bath. 13 S Imari Walker et al New J. Phys. 16 (2014) 063026 Although our examples are idealized and simplified, we contend that they capture an important and quite general feature of time-dependent quantum systems. We make no attempt here to discuss the practicality of measuring the weak values in the systems’ environment, but we are confident that realistic quantum systems that reveal this aspect can be found. Acknowledgements The authors wish to thank Lev Vaidman and Jeff Tollaksen for constructive conversations. This work has been supported in part by the Israel Science Foundation Grant No. 1125/10. References [1] Aharonov Y and Rohrlich D 2005 Quantum Paradoxes (Weinheim: Wiley-VCH) [2] Aharonov Y and Vaidman L 1990 Properties of a quantum system during the time interval between two measurements Phys. Rev. A 41 11 [3] Aharonov Y, Casher A, Albert D and Vaidman L 1987 Surprising quantum effects Phys. Lett. A 124 199 [4] Aharonov Y, Albert D and Vaidman L 1988 How the result of a measurement of a component of the spin of a spin- particle can turn out to be 100 Phys. Lett. A 60 1351 [5] Aharonov Y and Rohrlich D 1990 Towards a two vector formulation of quantum mechanics Quantum Coherence 221 ed J Anandan (Singapore: World Scientific) [6] Steinberg A M 1995 How much time does a tunneling particle spend in the barrier region? Phys. Rev. Lett. 74 2405 [7] Ferreira P J S G, Kempf A and Reis M J C S 2007 Construction of Aharonov and Berryʼs superoscillations J. Phys. A: Math. Theor. 40 5141 [8] Tollaksen J 2007 Probing contextuality with pre-and postselection J. Phys.: Conf. Ser. 70 012014 [9] Davies P C W 2009 Time-dependent quantum weak values: decay law for postselected states Phys. Rev. A 79 032103 [10] Tollaksen J, Aharonov Y, Casher A, Kaufherr T and Nussinov S 2010 Quantum interference experiments, modular variables and weak measurements New. J. Phys. 12 013023 [11] Aharonov Y, Popescu S and Tollaksen J 2010 A time-symmetric formulation of quantum mechanics Phys. Today 63 27 [12] Dressel J, Malik M, Miatto F M, Jordan A N and Boyd R W 2014 Understanding quantum weak values: basics and applications Rev. Mod. Phys. 86 307–16 [13] Doyle A C 1987 Silver Blaze and Other Stories (Portsmouth, NH: Heinemann) [14] Cox T and Lekner J 2008 Reflection and non-reflection of particle wavepackets Eur. J. Phys. 29 671–9 [15] Andreata M A and Dodonov V V 2004 Tunnelling of narrow Gaussian packets through delta potentials J. Phys. A: Math. Gen 37 2423–38 [16] van Dijk W and Nogami Y 2002 Analytical approach to the wave function of a decaying quantum system Phys. Rev. C 65 024608 [17] Yearsley J M 2008 The propagator for the step potential and delta function potential using the path decomposition expansion J. Phys. A: Math. Theor. 41 285301 14