Tunneling Time in Attosecond Experiments and Time Operator in Quantum Mechanics

Attosecond science is of a fundamental interest in physics. The measurement of the tunneling time in attosecond experiments, offers a fruitful opportunity to understand the role of time in quantum mechanics (QM). We discuss in this paper our tunneling time model in relation to two time operator definitions introduced by Bauer and Aharonov–Bohm. We found that both definitions can be generalized to the same type of time operator. Moreover, we found that the introduction of a phenomenological parameter by Bauer to fit the experimental data is unnecessary. The issue is resolved with our tunneling model by considering the correct barrier width, which avoids a misleading interpretation of the experimental data. Our analysis shows that the use of the so-called classical barrier width, to be precise, is incorrect.


Introduction
The paper concerns the time operator in relation to the definitions introduced by Bauer [1] and Aharonov-Bohm [2].Attosecond science (attosecond = 10 −18 s) concerns primarily electronic motion and energy transport on atomic scales.In previous work [3][4][5], we presented a tunneling model and a formula to calculate the tunneling time (T-time) by exploiting the time-energy uncertainty relation (TEUR), precisely that time and energy are a (Heisenberg) conjugate pair.Our tunneling time is in a good agreement with the attosecond (angular streaking) experiment for the He atom [3] with the experimental finding of Eckle et al. [6][7][8] and for the hydrogen atom [9] with the experimental finding of Sainadh et al. [10].Our model presents a real T-time picture or a delay time with respect to the ionization at atomic field strength F a (see below; compare Figure 1).It is also interesting for the tunneling theory in general, because of its relation to the height of the barrier [3,4].Indeed, the role of time has been controversial since the appearance of quantum mechanics (QM).The famous example is the Bohr-Einstein weighing photon box Gedanken experiment (BE-photon-box-GE) [11], in which a photon is allowed to escape from a box through a hole, which is closed and opened temporarily by a shutter.The period of time is determined by a clock, which is part of the box system.This means that the time is intrinsic and dynamically connected with the system under consideration.The total mass of the box before and after a photon passes is measured, which leads to uncertainty in the gravitational field, hence an uncertainty in time.Our tunneling time picture [3] shows an intriguing similarity to the BE-photon-box-GE (see for example [12], p. 132), where the former can be seen as a realization of the latter.
Concerning the time operator in QM, recently, Galapon [13][14][15] has shown that there is no a priori reason to exclude the existence of a self-adjoint time operator, canonically conjugate to a semibounded Hamiltonian, contrary to the (famous) claim of Pauli.The result is, as noted earlier by Garrison [16], for a canonically conjugate pair of operators of a Heisenberg type (i.e., uncertainty relation), that Pauli theorem does not apply; unlike a pair of operators that form a Weyl pair (or Weyl system.) Our simple tunneling model was introduced in [3] (see Figure 1,below).In this model, an electron can be ionized by a laser pulse with an electric field strength (hereafter, field strength) F. Ionization happens directly when the field strength is larger than a threshold called the atomic field strength F a = I 2 p /(4Z e f f ) [17,18], where I p is the ionization potential of the system (atom or molecule) and Z e f f is the effective nuclear charge in the single-active electron approximation.However, for field strengths F < F a , ionization can happen by a tunneling mechanism, through a barrier built by the effective potential due to the Coulomb potential of the nucleus and the electric field of the laser pulse.It can be expressed in a one-dimensional form: compare Figure 1.In this model, the tunneling process can be described solely by the ionization potential I p of the valence (the interacting) electron and the peak field strength F, which leads to the quantity δ z = I 2 p − 4Z e f f F, where F stands (throughout this work) for the peak electric field strength at maximum.In Figure 1 (for details, see [3]), the inner (entrance x e,− ) and outer (exit x e,+ ) point x e,± = (I p ± δ)/(2F), the barrier width d W = x e,+ − x e,− = δ z /F and the barrier height where the direct ionization starts, the barrier-suppression ionization.
In the (low-frequency) attosecond experiment, the laser field is comparable in strength to the electric field of the atom.Usually, intensities of ∼10 14 W•cm −2 are used.A key quantity is the Keldysh parameter [19], where ω 0 is the central circular frequency of the laser pulse and τ K denotes the Keldysh time.
In Equation ( 2) and hereafter, we adopt the atomic units (au), where the electron's mass and charge and the Planck constant are set to unity, h = m = e = 1.In Equation ( 2), at values γ K > 1, the dominant process is the multiphoton ionization (MPI).On the opposite side, i.e., for γ K < 1, the ionization (or field-ionization) happens by a tunneling process, which occurs for F < F a .This is the famous Keldysh result [19].In the tunneling process (for F < F a ), the electron does not have enough energy to ionize directly.The electron tunnels (tunnel ionizes) through a barrier made by the Coulomb potential and the field of the laser pulse.The electron escapes the barrier at the exit point x e,+ to the continuum, as shown in Figure 1 (a sketch for He atom).For our discussion later, we mention that in the experiment of [6], an elliptically-polarized laser pulse was used with ω 0 = 0.0619 au (λ = 736 nm) and with ellipticity = 0.87.The electric field strength was varied in the range F = 0.04 − 0.11 and for He atom I p = 0.90357 au.In the attosecond angular streaking experiment, one uses a close-to-circular polarized laser pulse; the direction of laser field ensures a unique relationship between the time at which the electron exits the tunnel and the direction of its momentum after the laser pulse.The measured momentum vector of the electron hence serves as the hand of a clock (attoclock), which indicates (the streaking angle) the time when the electron appeared from the tunnel in the laser field.It determines the time (streaking angle divided by the laser frequency) the electron has spent in the tunnel in the classically forbidden region inside the potential wall.This difference or delay is referred to as the tunneling delay time.It has to be mentioned that the nonadiabatic effects lead to large measurement error, reflected in large error bars.Compare [6], and see Section 3.3.The impact of deriving pulse parameters (shape, duration, amplitude) can, however, be controlled by choosing optimal pulse parameters, and the readout error can be minimized, as shown in [20,21], which is important due to the relation between tunneling and quantum information processing.x axes in a.u. -

The Tunneling Time
The tunneling process in the attosecond experiment can be visualized by a simple picture, as shown in Figure 1.This picture is related to the strong field approximation (SFA) or the Keldysh-Faisal-Reiss approximation [19,22,23]; for introductory reviews, see [24][25][26].The main idea of the SFA is that the tunneled electron escapes the barrier at the exit point x e,+ (see Figure 1), with zero kinetic energy.Precisely, the electron velocity along the (opposite) field direction is zero and negligible in the other directions.One neglects the effect of the atomic potential in the under-the-barrier region (see also [27].The electron wave function after x e,+ (the barrier exit) can be then approximated by a type of plane wave called the Volkov state.The result of SFA is of exponential accuracy (i.e., with an approximate pre-factor).It is, however, important to determine the final electron amplitude of the tunneling or ionization.This simple picture (our tunneling model), which is useful to explain the experiment, is strictly true only in the length gauge (see [28][29][30][31]).In a previous work [3], it was shown that the uncertainty in the energy, which is related to the height of the barrier h B (x m ), is quantitatively connected to the atomic potential energy at the exit point, ∆E ∼| V(x e ) |=|−Z e f f /x e |, for arbitrary strengths F ≤ F a .

The Tunneling Time
The tunneling process in the attosecond experiment can be visualized by a simple picture, as shown in Figure 1.This picture is related to the strong field approximation (SFA) or the Keldysh-Faisal-Reiss approximation [19,22,23]; for introductory reviews, see [24][25][26].The main idea of the SFA is that the tunneled electron escapes the barrier at the exit point x e,+ (see Figure 1), with zero kinetic energy.Precisely, the electron velocity along the (opposite) field direction is zero and negligible in the other directions.One neglects the effect of the atomic potential in the under-the-barrier region (see also [27].The electron wave function after x e,+ (the barrier exit) can be then approximated by a type of plane wave called the Volkov state.The result of SFA is of exponential accuracy (i.e., with an approximate pre-factor).It is, however, important to determine the final electron amplitude of the tunneling or ionization.This simple picture (our tunneling model), which is useful to explain the experiment, is strictly true only in the length gauge (see [28][29][30][31]).In a previous work [3], it was shown that the uncertainty in the energy, which is related to the height of the barrier h B (x m ), is quantitatively connected to the atomic potential energy at the exit point, ∆E ∼| V(x e ) |=|−Z e f f /x e |, for arbitrary strengths F ≤ F a .
By exploiting the TEUR ∆E • ∆T ≥ 1/2, we obtained what we call the symmetrical (or total) T-time: [3]: In this relation, δ z = I 2 p − 4Z e f f F for a single active electron model (for the hydrogen atom, Z e f f = 1); for details, see [3,4,9].The symmetrical T-time or the total time was obtained by a symmetrization procedure (similar to the Aharonov time operator [2,11]; see Section 3.2) from the unsymmetrized (unsymmetrical) T-time relation: Equation ( 3), besides its mathematical simplicity, aids in the conceptual reasoning [3,4].The physical reasoning is the following: the barrier itself causes a delaying time τ T,d , where: τ T,d presents the time delay with respect to the ionization at atomic field strength F a , where the barrier is absent (the barrier height, the barrier width d W = δ z /F and δ z are zero).It is the time duration for a particle to pass the barrier region (between x e ,− , x e,+ ) and escape at the exit point x e,+ to the continuum [3].The first term τ T,i in Equation ( 3) is the time needed to reach the entrance point x e,− from the initial point x i ; compare Figure 1.At the limit F → F a (δ z → 0), the total time becomes the ionization time For F > F a , the barrier-suppression ionization sets up [32,33].At the opposite side, i.e., for F → 0, δ z → I p and τ T,d → ∞, hence nothing happens (undisturbed electron), and the electron remains in its ground state, which shows that our model is consistent.For details, see [3][4][5].

Time Operator
In the early days of the QM, a TEUR (analogous to the position-momentum relation) and the existence of the time operator faced the well-known objection of Pauli.According to the Pauli theorem, the introduction of the time operator is basically forbidden, and the time t in QM must necessarily be considered as an ordinary number ('c' number).On the other hand, the famous example and its debate is the BE-photon-box-GE.The crucial point is that, to date, no general time operator has been found [34]; thus, a sort of time operator and the uncertainty relation are used dependent on the study cases [35].There is still a common opinion that time plays a role essentially different from the role of the position in QM (although, this is not in line with special relativity).Hilgevoord concluded in his work [36] that when looking to a time operator, a distinction must be made between universal time coordinate, the t, a c-number like a space coordinate and the dynamical time variable of a physical system situated in space-time; i.e., clocks.The search for a time operator has a long history [37][38][39].This led Busch [40,41] to classify three types of time in QM: external time (parametric or laboratory time), intrinsic or dynamical time and observable time.In our case study, the time, and hence the tunneling time, is intrinsic or dynamical type [4].

Bauer's Time Operator
Bauer [1,42] introduced a dynamical self-adjoint time operator in the framework of Dirac's relativistic quantum mechanics (DRQM).Bauer's time operator (BTO) is defined as the following, where α, β are the well-known Dirac matrices, c the speed of light and r the three-dimensional space vector.τ 0 represents in principle an internal property of the system, determined to be the de Broglie period h/(m 0 c 2 ) (h is the Planck constant, and m 0 is the rest mass of the particle).τ 0 is not important for our discussion, since it cancels for time intervals, or time differences.
The operator defined in Equation ( 6) has been shown to commute with the Dirac free particle Hamiltonian H D = c α • p + βm 0 c 2 .Bauer proved the Heisenberg commutation relation, analyzed the dynamical character of T and found for a free particle: where v gp is the group and v ph (Equation ( 9), see Bauer [42]) is the phase velocity of the particle.According to Bauer, T(t) is the internal time and t the external (laboratory) time, where T(t) < t.
In the limit, m 0 = 0, then v gp = c and T(t) = t.The parametric time t, according to Bauer, can be interpreted as the laboratory time, which is the time variable appearing in the time-dependent Schrödinger equation and characterizing the dynamical evolution of microscopical systems Bauer argued that DRQM allows the introduction of a dynamical time operator that is self-adjoint, unlike the parametric time entering in the Schrödinger equation, while Galapon [13][14][15] showed that in non-relativistic QM (NRQM), there is no reason to exclude the existence of a time operator for a Heisenberg pair, and consequently an observable of time and a dynamical time, as already mentioned.
For a non-relativistic particle, Bauer found [1] a relation between the parametric (external in Bauer's notation) time intervals and the dynamical (internal in Bauer's notation) intervals, The parametric time interval ∆t is enhanced relative to the internal time interval ∆T, which by the virtue of Equation ( 6) is related to the time the light takes to travel the same distance.On the other hand, for high ("relativistic") energies, one obtains: In other words, in the relativistic case, internal time intervals coincide with external time intervals, whereas in the non-relativistic case, the latter is enhanced relative to the internal time intervals.At this point, it is important to note that in the presence of a potential dependent only on position, e.g., Coulomb type potentials, hence, the commutation relation of the time operator is reduced to the commutation relation with the relativistic free particle of Dirac operator ĤD , since the latter is a linear function of relativistic momentum p. From Equation ( 6), one finds that Equation ( 11) is reduced to the position momentum commutation relation [r, p].For the tunneling in the attosecond experiment, Bauer uses the argument of Kullie [3] that the potential energy at the exit point defines the uncertainty in the energy, leads to the time of passage of the barrier or the time needed to cross through the exit point and represents a tunneling internal time of the system τ T .With his view of ∆T B τ T , Bauer obtained a relation for the laboratory time lapse to cross the barrier, where Γ = (v gp /c) −2 is called the enhancement factor, ξ is a phenomenological parameter (see below) and τ T the internal time interval.Compare Equation (10).

Properties of BTO
The BTO is interesting and satisfies the property and the conditions of an ordinary time operator in the relativistic framework of the QM.Despite this, there is an unexpected feature of the relation (12) given by Bauer, as a consequence of the time operator definition in (6).The laboratory time interval ∆t is connected to the internal time interval ∆T ≡ ∆T B by Equation ( 12).However, with Equations ( 6), ( 8) and ( 10) and similar to Equation ( 12): Consequently, one finds: = ξ c Bauer introduced the phenomenological parameter ξ < 1 so that Υ T in Equation ( 12) somehow fits the experimental data and the Feynman path integral (FPI) calculation, presented by Landsman [6]; compare Figure 1 of [1] (the same plot is given in Figure 3, Section 3.3).He concluded that ξ = 0.45 gives the best fit to the experiment.The result of Bauer fits the FPI result of Landsman [6] well; see Figure 1 in [1] (as in Figure 3, Section 3.3).However, there is no theoretical justification for the choice of ξ as a parameter or its value ξ = 0.45.This is a rather unexpected result, then without the parameter ξ (i.e., ξ = 1), one obtains v gb = 2c, which is a consequence of the definition in Equations ( 6) and (12), where the speed of light is used in the denominator.In other words, the value ξ = 0.5 gives v gb = c, as it should be in accordance with the definition in Equation (6).The idea of Bauer to set a universal internal time Equation ( 6) is reasonable; however, using it to measure external time intervals, i.e., the relation between internal time intervals and external time intervals Equation ( 12), leads to unexpected implications.We will see that our tunneling model and our definition, or the generalization of time operators of Bauer and Aharonov (see below), clarifies this implication and that there is no need to introduce a phenomenological parameter.

Time Operator of the Type Aharonov-Bauer
The recent interesting work of Bauer and the introduction of a time operator in the frame work of the DRQM, together with the well-known Aharonov time operator (Equation ( 19) below), have stimulated us to define a generalized form of the time operator(s) of the same types.We suggest that the time operator definition of Aharonov and Bauer can be extended and combined as the following.For one particle (in atomic units): where v ph , v gp are the phase and the group velocity, respectively.In the following, we look to a one-dimensional case with the radial coordinate r, i.e., we neglect the factor 1/(4π) used by Bauer of the three-dimensional case, where r = 4π r.The Dirac matrices α will be set to unity α = 1.Under the definition in Equation ( 16), the notation internal, external time operator is misleading.We denote tgp the dynamical time operator (dynamical TO) and tph the phase (or phase-velocity) time operator (phase TO), without any connection to the notation external, internal time classification of Busch [35].The relation between the dynamical and phase times follows immediately from Equation ( 16), using an interval ∆r = ∆t ph Because of the relativistic relation v ph v gp = c 2 for a matter particle with a mass m, E = mc 2 , v ph = E/p, one obtains tgp ≥ tph , where p is the momentum of the particle.For this reason, it is better to adopt the notations: tph the phase and tgp the dynamical time operator.For the light (photon) particle v gp = v ph = c, one obtains the BTO: For the non-relativistic one-dimensional case (α = 1), it is straightforward to obtain the well-known Aharonov-Bohm [2] time operator (ABTO) ( tgp = r/v gp ) or the quantum mechanical symmetric operator: For tph , the symmetrization has no meaning, because v ph is not an observable, unlike v gp .The commutation relations were verified in both cases, the non-relativistic Bohm-Aharonov operator [2,43] and the relativistic Bauer operator [1,42].For tph , because v ph is not an observable, the commutation relation is reduced to the known commutation relation [r, p].
The first consequence of our definition is the equivalence to the BTO, as given in Equation ( 18), and to the ABTO Equation (19).There is also an equivalence between BTO and ABTO; then, from Equation ( 19), for a light particle with v ph = v gp = c: and for a matter particle v ph • v gp = c 2 : as already found by Bauer.Compare Equations ( 7), ( 12) and (13).No approximation is used, but with ∆t gp instead of the parametric time ∆t and ∆t ph instead the internal time ∆T of Bauer's notations.Bauer obtained the factor ξ (1/2) −1 = 2 ξ in Equation ( 12) by going from relativistic to non-relativistic approximation (the factor (1/2) −1 ) and using a phenomenological parameter ξ to fit the experimental data.tph is the phase time, whereas Bauer refers to it as internal time T B (using v ph = c).In the Bauer case (notation), one has ∆T B ≤ ∆t B ; see Equations ( 10) and (12).Likewise, we have a relation between the dynamical and phase times, ∆t ph ≤ ∆t gp , Equation ( 21).However, one has T B = r/c ≥ t ph = r/v ph because v ph ≥ c.Further, we get the Bauer procedure for the three-dimensional case by the replacing ∆r/c with ∆r/v ph , which with ∆r = 4π∆r, gives Equation ( 13): The only difference between ( 13) and ( 22) is that Bauer used 2 ξ = 2 × 0.45 = 0.9 (compare Equations ( 12) and ( 13)), whereas it is equal 2 × 0.5, or exactly ξ = 1/2 in our case without any approximation.The difference is very small, and ξ = 0.45 does not fit perfectly to the experimental data of [6]; see below.Thus, referring to ∆T = ∆r/c as an internal time and introducing a phenomenological parameter ξ have no justification, unless one relates every time interval to its counterpart of a unique time interval of the light propagation, which is ∆r/c by replacing the phase time of a matter particle by v ph = c; that is straightforward, and a parameter ξ is redundant.
With his approximation, Bauer obtained Equation ( 12) or ( 13), whereas on the basis of our tunneling model, Equation ( 5) can be written, after a small manipulation, in the form (F is the field strength): where d W = δ z /F is the barrier width and x e,− = (I p − δ z )/(2F) [3]; compare Figure 1.In the following, to avoid a confusion, we refer in the general case to a barrier width as D BW ; Whereas in our model, we set D BW = d W , and for the classical barrier width we set D BW = d C .On the basis of numerical values from the experimental data in [6], with approximate barrier width d C (compare [6] with [1]) and with the values d C = 13 au (in atomic units), τ = 40 as (in attosecond), F ≈ 0.069 au, v gb = 13/40 (au/as) = 6.88/40 (Å/as), Bauer obtained [1], = 16.14 where D BW is a barrier width, Å is Angström length unity and au atomic units.The factor 44 is calculated by using the numerical data above, and the best fit to the experimental data according to Bauer is ξ = 0.45.One notes that Landsman [6] assumed a classical barrier width, i.e., D BW = d C = I p /F which is larger than d W (compare Figure 1).It is usually taken to be approximately valid and is adopted by Bauer.However, as we will see below, this leads to a confusion in the evaluation of the tunneling time against the barrier width; whereas from our model Equations ( 5) and ( 23), it follows that τ d = 3.6 (d W + x e,− ) as, using Z e f f = 1.6875 of Clementi [44] (for small barrier width), or τ d = 4.4 (d W + x e,− ) as, using Z e f f = 1.375 of Kullie (for large barrier width) (see [3]); however, no fitting procedure is used, and our τ d is in good agreement with the experimental data.One can imagine that BTO, Equation (6), presents a universal time scale or internal clock of a light particle, a photon, but its relation to the external time or clocks is then not presented by Equations ( 12) and ( 13); see below.We think that our definition is a generalized form, Equation ( 16), with a straightforward transition to both the BTO and the non-relativistic ABTO.

Experimental Affirmation
It is worthwhile to mention that the velocity with which the electron passes through the tunnel varies slowly with the barrier width.In our model and with some manipulation, one obtains the mean velocity as a function of d W : For Z e f f = 1.6875, respectively 1.375, we obtain the values v gp (d W ) < 4Ze e f f = 6.75, respectively 5.5.Using the experimental data at one point (see Equation ( 24)), Bauer extracted a value (2 × ξ/(v gb /c) 2 ) = ξ × 608.44 (compare Equation (22) or v C gp = (d C /τ)| F=0.069 = (13 au)/(40 as) = (13/1.66)au = 7.86 au, which was used as a fixed (independent of F) mean value in Equation ( 24) independent of the barrier width.It is sightly larger than our values 6.75 (5.5), which is caused by the use of the classical barrier width d C = I p /F = d W + 2x e,− (compare Equations ( 12) and ( 24)), and we think it is one of the reasons why Bauer needed to introduce a phenomenological parameter.
In Figure 2, we plot our result of tunneling time Υ K (d W ), Equation (23), against the barrier width d W (F) = δ z /F, where Z e f f = 1.6875 is used, together with Bauer's tunneling time formula Υ B (D BW = d W ), Equation (24), with ξ = 0.45 used by Bauer, with our value ξ = 0.5 and the values ξ = 0.3, 0.6 used by Bauer.In addition, the experimental tunneling time data (plotted against D BW = d W ) together with the error bars are shown.The elementary data or the experimental values of time of [6] (at the corresponding F values) were sent by Landsman, where d W was easily calculated from d W = δ z /F.As seen in Figure 2, Equation ( 23) (solid red line) shows the best fit with the experimental data.The dashed (pink) line of Equation (24) with D BW = d W and ξ = 0.5 is slightly below the red solid line (for small d W ), and the dashed dotted (blue) one of Equation ( 23) with D BW = d W and ξ = 0.45 is slightly below both.
One notices that we used Z e f f = 1.6875, which is suitable for small barrier widths, whereas Z e f f = 1.375 is better for large barrier widths, with which the lines will slightly shift towards higher values.The red line is then closer to the experimental values for large barrier width (compare Figure 5 in [3]), and the pink and blue lines stay below the red line in this region.
Figure 2 has a small difference from the figure shown in [6] (Figure 3d, FPI) and to Figure 1 from [1] (same as Equation (24), with ξ = 0.3, 0.45, 0.6).Because in both works, D BW = d C = I p /F = d W + 2x e,− was used, where the plotted lines are slightly above the experimental data, i.e., the lines showed less agreement than in our evaluation plotted in Figure 2. In other words, a parameter ξ is not needed, when one uses the correct barrier width d W = δ/F (compare Figure 1), and a better agreement with the experimental data is obtained.Bauer eq (24) D BW =dw, ζ=0.5 Bauer eq (24) D BW =d W ζ=0. 45 Bauer eq (24) D BW =d W ζ=0. 3 Bauer eq (24) d BW =d W ζ=0.6 with ξ = 0.3, 0.45, 0.5, 0.6.The experimental data are also shown, at the corresponding barrier width d W with the error bars.See [3].The experimental data were kindly sent by Landsman; see [6].
For clarity, we plot in Figure 3 a copy of Figure 1 of [1] (extracted data by using a web digitizer), i.e., the data of Bauer Υ B (D BW = d C ) for ξ = 0.45 (dashed, black) and also the FPI (solid, light blue), which was used by Bauer after it was extracted from Figure 3d of [6].Additionally, we plot Υ B (D BW = d C ) with ξ = 0.45 (dashed doted, magenta), which as expected, reproduce the line of Bauer (dashed, black); the tiny difference is only because we extracted the data of Figure 1 of Bauer [1] by a web digitizer.Furthermore Υ B (D BW = d C ) with ξ = 0.5 is plotted (dashed dotted dotted, blue), which lies slightly higher.We can reproduce the data of Bauer, which makes our conclusion reliable.Thus, we can see why Bauer found that ξ = 0.45 fits better to FPI; this is because the use of the approximate barrier width D BW = d C ; precisely speaking, the use of D BW = d C is incorrect.The small difference was not crucial for Landsman in the work [6] (Figure 3d), but Landsman noted that it is an approximate barrier width, unlike our model [3] (published later), where a correct barrier width D BW = d W = x e,+ − x e,− = δ z /F was obtained; compare Figure 1.Our conclusion is that, although the difference between d W and d C seems to be not crucial as thought by Landsman [6] (and many other authors, for example [10,45]; see also [9]), the use of a barrier width D BW = d C instead of d W leads to confusion and is misleading; especially when plotting the tunneling time data against the barrier.This led Bauer to introduce a parameter with the value ξ = 0.45, which in our view is unnecessary, regardless of how one understands Bauer's definition of "internal" time operator in Equation ( 6).Thus, our definition of a time operator Equation ( 16) is reasonable; it is a general form with a straightforward transitions to BTO Equation ( 6) [1,42] and ABTO [2].Eq (24),with D BW =d C , ζ=0.5 Eq (24), with D BW =d C ζ=0.45As a final note, we demonstrate the importance of our model to the tunneling theory in general, because it relates the tunneling time to the barrier height.The T-time found in Equation (3) can be also derived in a simple way, when we assume that the barrier height corresponds to a maximallysymmetric operator as the following.The barrier height h B (x m ) ≡ h M can be related to an (real) operator h B (x m ) and the uncertainty in the energy caused by the barrier: From this, we get ∆E ± = abs(−I p ± δ z ), i.e., we have to add (subtract) the internal energy of the system (the ionization potential).Hence, we get τ T,d , τ T,i and τ sym T by the virtue of the time-energy uncertainty relation, where we assumed that the time is intrinsic and has to be considered (a delay time) with respect to the ionization (time) at atomic field strength, i.e., with respect to the internal energy or the ionization potential I p .This suggests to consider h M as a perturbation (energy) operator, where the full operator of the system can be taken as H 0 + h ± M .In fact, one can argue that −I p must be taken to avoid the divergence of the time to infinity for F = F a , because it is physically incorrect, as δ z (F a ) = 0 ⇒ τ = 1/(2δ z ) → ∞, which in turn, can be seen as an initialization of the internal clock, i.e., the T-time is counted as a delay with respect to the ionization at F a (the limit of the subatomic field strength), after which no tunneling occurs, and the ionization is classically an allowed process, the barrier-suppression ionization.
Actually, we see immediately from

Conclusions
In this work, we showed that our tunneling time model for the tunneling in the attosecond angular streaking experiment enables us to discuss and understand the time and time operator in quantum mechanics.We found a generalized form of a time operator with straightforward transitions to the Bauer time operator, which was introduced in the framework of relativistic Dirac theory, and to the Aharonov-Bohm time operator, which was introduced in the non-relativistic quantum mechanics.Moreover, we found that the introduction of a phenomenological parameter as done by Bauer is unnecessary.The issue is resolved, where the confusion was caused by the use of the classical barrier width.It is clarified by using the correct barrier width, as found by our tunneling model.

Funding:
The work is supported by the Open Access Publishing Fund of the University Library Kassel.

Figure 1 .
Figure 1.(Color online) Graphic display of the potential and the effective potential curves, the two inner and outer points x e,± = (I p ± δ z )/2F and the barrier width d W = x e,+ − x e,− = δ z /F.δ z = I 2 p − 4Z e f f F, and I p is the ionization potential, Z e f f the effective nuclear charge and F the electric field strength of the laser pulse at maximum.x e,c = I p /F ≡ d C is the "classical exit" point, x m (F) = Z e f f /F is the position at the maximum of the barrier height h(x m ), x a = x m (F = F a ) and F a is the atomic field strength; see the text.The plot is for the He atom in the single-active-approximation model with Z e f f = 1.6875 and I p = 0.90357 au.Note, for atoms with different Z e f f and I p , the overall picture does not change.

Figure 1 .
Figure 1.(Color online) Graphic display of the potential and the effective potential curves, the two inner and outer points x e,± = (I p ± δ z )/2F and the barrier width d W = x e,+ − x e,− = δ z /F.δ z = I 2 p − 4Z e f f F, and I p is the ionization potential, Z e f f the effective nuclear charge and F the electric field strength of the laser pulse at maximum.x e,c = I p /F ≡ d C is the "classical exit" point, x m (F) = Z e f f /F is the position at the maximum of the barrier height h(x m ), x a = x m (F = F a ) and F a is the atomic field strength; see the text.The plot is for the He atom in the single-active-approximation model with Z e f f = 1.6875 and I p = 0.90357 au.Note, for atoms with different Z e f f and I p , the overall picture does not change.

Figure 2 .
Figure 2. (Color online) Graphic display.The x-axis is the barrier width d W in au and the vertical axis the tunneling time in as of formula Equation (23) Υ K (d W ), and of formula Equation (24) Υ B (D BW = d W )with ξ = 0.3, 0.45, 0.5, 0.6.The experimental data are also shown, at the corresponding barrier width d W with the error bars.See[3].The experimental data were kindly sent by Landsman; see[6].
as Expr.data extracted from fig 1 of ref[1] Data FPI from fig 1 of ref[1] Data of Bauer, ζ=0.45 from fig of ref[1]

Figure 3 .
Figure 3. (Color online) Graphic display.The horizontal axis (classical) is barrier width d C in au and the vertical axis the tunneling time in as.The figure reproduces Figure 1 of[1] (using a figure digitizer), the tunneling time of the Bauer formula for ξ = 0.45, the Feynman path integral (FPI) and the experimental data; note that these data are extracted from Figure1of[1]); compare Figure3dof[6].Together with the formula in Equation (24) Υ B (d C ) (same formula of[1]), for ξ = 0.45, ξ = 0.5, we reproduced the result of Bauer.
Figure1that the maximum of the effective potential curve is equal to −2 Z e f f F[3]; it goes lower with increasing field strength, where at F = F a , it becomes−2 Z e f f F = −I p or I 2 p − 4Z e f f F = 0. Hence, for F ≤ F a , one gets I 2 p − 4 Z e f f F = δ 2 z ≥ 0.Now, the distance (or the height) to the −I p (horizontal line) level (at maximum) is then |h| = δ 2 z .δ z presents the change of the potential energy in the barrier region (as it should be) ∆V(x)| x e,+ x e,− = −Z e f f /x e,+ − (−Z e f f /x e,− ) = δ z .Hence, the energy gap ∆E ± = abs(−I p ± |h|); from which we get the tunneling time τ T,d and τ T,i by virtue of the TEUR.