## Second-order nonlocal shifts of scattered wave-packets: What can be measured by Goos-H¨ anchen and Imbert-Fedorov effects ?

Klaus Morawetz 1 2 , ,

1 M¨ unster University of Applied Sciences, Stegerwaldstrasse 39, 48565 Steinfurt, Germany and 2 International Institute of Physics- UFRN, Campus Universit´rio Lagoa nova, 59078-970 Natal, Brazil a

The scattering of wavepackets with arbitrary energy dispersion on surfaces has been analyzed. Expanding up to second order in scattering shifts, it is found that besides the known Goos-H¨nchen a or Imbert-Fedorov spatial offset, as well as the Wigner delay time, new momentum and frequency shifts appear. Furthermore, the width of the scattered wave packet becomes modified as well, which can lead to a shrinking of pulses by multiple scattering. For a model of dielectric material characterized by a longitudinal and transverse dielectric function the shifts are calculated analytically. From the Goos-H¨ anchen and Imbert-Fedorov shifts one can access the longitudinal and transversal dielectric function. Perfectly aligned crystal symmetry axes with respect to scattering beam shows no Imbert-Fedorov effect. It is found that the Goos-H¨nchen and Imbert-Fedorov effect are absent a for homogeneous materials. Oppositely it is found that the Wigner delay time and the shrinking of the temporal pulse width allows to access the dielectric function independent on the beam geometry.

## I. INTRODUCTION AND RESULTS

Since the experiments of Goos and H¨ anchen, it has been known that a wave packet or a beam of light suffers a nonlocal shift when reflecting from surfaces [1]. The shift in the focal plane is called the Goos-H¨nchen effect, a and the shift out of the plane is called the Imbert-Fedorov effect [2]. For an overview, see [3, 4]. Both deviations from geometrical optics predictions can be distinguished and depend on the shape of the incident beam, its polarization, and the material composition of the reflecting surface [5]. Treatments consider beam shifts for pairs of plane waves [6] or show that a classical spinning photon yields an 'exotic particle' on a curved surface [7]. Experimentally, it was shown that the degree of spatial coherence influences the angular beam shifts, while the spatial beam shifts are unaffected [8]. The effect of beams with orbital angular momentum was clarified in [9, 10]. Connected with these effects are the spin separations of light for femtosecond laser pulses due to the spin-Hall effect [11]. As special forms of beams rotating elliptical Gaussian [12] or Airy beams [13] are used. Airy vortex and Airy vector beams are created by modulation of dynamic and geometric phases [14]. Circular Airy vortex beams can be created in the terahertz regime [15]. Large spatial shifts of a reflected Airy beam on the surface of hyperbolic crystals were found [16]. Reflection and transmission of an Airy beam impinging on a dielectric surface has been investigated in [17]. Other applications of these shifts consider graphene [18, 19] or optical vortex beams [20]. The effect of an independent quantum degree of freedom on the barycenter of a diffraction-free light beam was calculated in [21].

The theoretical basis of these shifts is the energy and/or momentum dispersion of the wavepacket or beam. This results in nonlocal shifts when wavepackets are scattered from a surface [22]. Quite frequently the treatment averages over the Fresnel coefficients [12, 23] or uses the transmission coefficients [24, 25]. There exist various other schemes to describe such shifts ranging from averaging over the center of mass of field energy density [26] to averaging over the incoming and outgoing fields [27], analogously to quantum expectation values [28]. A complete quantum kinetic theory, including nonlocal shifts, can be found in [29-31]. All quantum effects of scattering can be recast in a set of nonlocal spatial, temporal, momentum, and energy shifts. This results in contributions to the thermodynamic variables due to binary correlations [32].

In this letter, we consider the various offsets by expanding up to second order in various dispersions. The formalism can be applied to any surface, be it metals, molecules, dielectric, or other materials with any dispersion ω k ( ) of beams. We assume simply that the scattering is described by a proper scattering amplitude and will show that an incident wave packet scattered at a surface becomes modified by six effects. These effects can be expressed as nonlocal shifts in terms of derivatives of the scattering amplitude

$$f \left ( \frac { \mathbf p } { p }, \omega \right ) = | f | e ^ { i \Psi }.$$

One needs simply derivatives with respect to energy

$$\delta = \varphi + i \Delta = \partial _ { \omega } \ln f = \partial _ { \omega } \ln | f | + i \partial _ { \omega } \Psi \quad \ ( 2 )$$

and the vector shifts

$$\underset { \text{ee of} } { \text{ense} } \quad \delta = \varphi _ { r } + i \Delta _ { r } = \frac { 1 } { k } \partial _ { \theta, \phi } \ln f = \frac { 1 } { k } \partial _ { \theta, \phi } \ln | f | + \frac { i } { k } \partial _ { \theta, \phi } \Psi \left ( 3 \right )$$

with the orbital derivatives

$$\partial _ { \theta, \phi } = - \mathrm e _ { \theta } \partial _ { \theta } - \frac { \mathrm e _ { \phi } } { \sin \theta } \partial _ { \phi } = \mathbf k \partial _ { k } - k \partial _ { k }. \quad \ \ ( 4 )$$

As a result, the scattered wavepacket with velocity v k = ∂ ω k k

1. becomes delayed by the Wigner delay time ∆.

- 2. obtains a spatial offset ∆ r , which includes the Goos-H¨ anchen

$$\Delta ^ { G H } = \Delta _ { x } + \tan \vartheta _ { 0 } \Delta _ { z } \quad \quad ( 5 )$$

and Imbert-Fedorov effects

$$\Delta ^ { I F } = \Delta _ { y }$$

for light if the incident angle to the z-axes is ϑ 0 .

- 3. suffers a shift of momentum

$$\mathbf K = k - \sigma ^ { 2 } ( \varphi _ { r } - \mathbf v _ { \mathbf k } \varphi )$$

with the momentum width of the wavepacket σ .

- 4. experiences a change in the temporal width of the packet by

$$\bar { \sigma _ { t } } ^ { 2 } = \sigma _ { t } ^ { 2 } - \partial _ { \omega } \varphi.$$

- 5. shows a frequency shift of

$$- \frac { \varphi } { \bar { \sigma } _ { t } ^ { 2 } }.$$

- 6. gets a modified momentum width

$$\frac { 1 } { \bar { \sigma } ^ { 2 } } = \frac { 1 } { \sigma ^ { 2 } } + \frac { \varphi _ { \theta } - v _ { k } \varphi } { k } + i \left ( \frac { \Delta _ { \theta } - v _ { k } \Delta } { k } + \Gamma t \right ) \quad ( 7 ) \quad \.$$

with the second-order derivative of dispersion Γ = ∂ k 2 ω k .

For elastic impurity scattering, the modulus | f | is constant, and effects 3, 4, and 5 are not present. Effects 1 and 2 are independent on the form of the beam or wavepacket.

In the following chapter we give a derivation of these results and calculate the shifts for a model of light scattering at a dielectric material in chapter III.

## II. SCATTERING OF A WAVE PACKET ON A SURFACE

We consider a three-dimensional incident wave packet

$$\Psi _ { \text{in} } ( \mathbf r, \mathbf k, t ) & = \int \frac { d ^ { 3 } p d \omega } { ( 2 \pi ) ^ { 4 } } \text{e} ^ { i \mathbf p \cdot \mathbf r - \frac { ( p - k ) ^ { 2 } } { 2 \sigma ^ { 2 } } - i \omega t - \frac { \sigma _ { 1 } ^ { 2 } } { 2 } ( \omega - \omega _ { p } ) ^ { 2 } } \quad ( 8 ) \quad \text{term} \\ \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots$$

with any disperion ω p which scatters at a surface with the scattering amplitude f ( p p , ω )

$$\Psi _ { \text{out} } ( { \mathbf r }, { \mathbf k }, t ) = & \int \frac { d ^ { 3 } p d \omega } { ( 2 \pi ) ^ { 4 } } \text{e} ^ { i \widetilde { \mathbf p } \cdot { \mathbf r } - \frac { ( \mathbf p - { \mathbf k } ) ^ { 2 } } { 2 \sigma ^ { 2 } } - i \omega t - \frac { \sigma ^ { 2 } } { 2 } ( \omega - \omega _ { p } ) ^ { 2 } } f \left ( \frac { \mathbf p } { p }, \omega \right ) \quad \cdot$$

where the outgoing momenta ¯ p interchanges the sign of the z-component according to the reflection at a surface z = 0. Since we have ¯ p · r = ¯ r · p we work with ¯ r in the following. Since the wavepacket is sharply peaked around p ≈ k and ω ≈ w p we expand in two steps. First we expand the energy up to second order

$$f \left ( \frac { \mathbf p } { p }, \omega \right ) & = f \left ( \frac { \mathbf p } { p }, \omega _ { p } \right ) \left [ 1 + \frac { \partial _ { \omega _ { p } } f } { f } ( \omega - \omega _ { p } ) + \frac { \partial _ { \omega _ { p } } ^ { 2 } f } { 2 f } ( \omega - \omega _ { p } ) ^ { 2 } \right ] \\ \cdot \vartheta _ { 0 }. & = f \left ( \frac { \mathbf p } { p }, \omega _ { p } \right ) \, e ^ { \delta ( p ) ( \omega - \omega _ { p } ) + \frac { \partial _ { \omega \delta } } { 2 } ( \omega - \omega _ { p } ) ^ { 2 } } \quad ( 1 0 ) \\ & \quad \text{where we rewrote the Tavlor expansion as exponential}$$

where we rewrote the Taylor expansion as exponential leading to the energy derivative of the shift δ for the second order term. The appearing shifts have to be expanded around p ≈ k to provide

$$\L f \text{ the } \quad \delta ( p ) \ & = \ \delta ( k ) + ( \mathbf p - \mathbf k ) \cdot \partial _ { k } \delta \\ & = \ \delta ( k ) + ( \mathbf p - \mathbf k ) \cdot \mathbf v _ { k } \partial _ { \omega _ { k } } \delta + \frac { 1 } { k } ( \mathbf p - \mathbf k ) \cdot \partial _ { \theta, \phi } \delta \\ \partial _ { \omega _ { p } } \delta ( p ) \ & = \ \partial _ { \omega _ { k } } \delta ( k )$$

up to second order. Here we used

$$\partial _ { k } = \mathrm e _ { k } \partial _ { k } + \frac { 1 } { k } \partial _ { \theta, \phi } = \mathbf v _ { k } \partial _ { \omega _ { k } } + \frac { 1 } { k } \partial _ { \theta, \phi }. \quad ( 1 2 )$$

In the second step we expand the scattering amplitude

$$f \left ( \frac { \mathbf p } { p }, \omega _ { p } \right ) = f _ { k } \mathrm e ^ { \frac { ( p - k ) \cdot \partial _ { k } f } { f } + ( p - k ) \cdot \partial _ { k } \frac { ( p - k ) \cdot \partial _ { k } f } { 2 f } }. \quad ( 1 3 )$$

With the shifts introduced in (2) and (3) we have

$$\frac { \partial _ { k } f } { f } = \mathbf v _ { k } \delta - \delta \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \ \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \ }$$

and a helpful relation between the second derivatives

$$v _ { k } \partial _ { \theta, \phi } \delta = - \delta _ { r } - v _ { k } \delta ^ { \prime } _ { r }.$$

In the following we abbreviate δ ′ = ∂ ω k δ . Now we calculate

$$& \mathbf A \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
 \quad \text{set} \quad & \quad \left ( \mathbf p \cdot \mathbf e _ { k } \partial _ { k } + \frac { 1 } { k } \mathbf p \cdot \partial _ { \theta, \phi } \right ) \left ( \mathbf p \cdot \mathbf v _ { k } \delta + \frac { 1 } { k } \mathbf p \cdot \frac { \partial _ { \theta, \phi } f } { f } \right ) \quad ( 1 6 ) \\ & \quad \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
 \ \quad \text{$m$} \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

term by term. The first one leads to

$$p \cdot e _ { k } \partial _ { k } ( \mathbf p \cdot \mathbf v _ { k } \delta ) = ( \mathbf p \cdot e _ { k } ) ^ { 2 } \Gamma \delta + ( \mathbf p \cdot \mathbf v _ { k } ) ^ { 2 } \delta ^ { \prime } \quad ( 1 7 )$$

with Γ = ∂ k 2 ω k . The second one

$$\begin{array} { c } \text{with $1-\upsilon_{k}\ \omega_{k}$.                                                                                                                                                                                                         \quad                                                                                                                                                                                                       
                                                                                                                                                                                                        }                                                                                                                                                                                                        < \quad \quad \quad$$

Here we used k ( ∂ θ,φ ) ( i e k ) j = ( e θ ) ( i e θ ) j +( e φ i ) ( e φ j ) and the relation (15) in the first step. The third term reads

$$\mathbf p \cdot \mathbf e _ { k } \partial _ { k } \left ( \frac { \mathbf p \cdot \partial _ { \theta, \phi } f } { k f } \right ) = \frac { \mathbf p \cdot \mathbf v _ { k } } { k v _ { k } } \mathbf p \cdot \delta _ { r } - ( \mathbf p \cdot \mathbf v _ { k } ) ( \mathbf p \cdot \delta _ { r } ^ { \prime } ). \ ( 1 9 ) \quad \ \text{c}$$

The only problematic term is the last one. Due to (3) the spatial shifts have only two orbital components δ r = e θ δ rθ + e φ δ rφ and we obtain

$$& \frac { 1 } { k } \mathbf p \cdot \partial _ { \theta, \phi } \left ( \frac { 1 } { k } \mathbf p \cdot \frac { \partial _ { \theta, \phi } f } { f } \right ) = \\ & - \frac { 1 } { k } \left ( p _ { \theta } ^ { 2 } \partial _ { \theta } \delta _ { r \theta } + 2 p _ { \theta } p _ { \phi } \frac { \partial _ { \theta } \delta _ { r \phi } } { \sin \theta } + p _ { \phi } ^ { 2 } \frac { \partial _ { \phi } \delta _ { r \phi } } { \sin \theta ^ { 2 } } \right ) \quad & \text{W} \\ & = - \frac { 1 } { k } \left ( p ^ { 2 } - \frac { ( \mathbf p \cdot \mathbf v _ { k } ) ^ { 2 } } { v _ { k } ^ { 2 } } \right ) \partial _ { \theta } \delta _ { r \theta } \quad & ( 2 0 ) \quad & \text{inar} \quad \\ \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \end{cases}$$

since we can rotate the coordinates of p -integration such that p φ = 0 and can express p θ 2 = p 2 -( p v · k ) 2 v 2 k .

Collecting all four terms together we obtain

$$& p \cdot \partial _ { k } \frac { \mathbf p \cdot \partial _ { k } k f } { f } = \begin{array} { c } \underset { \begin{array} { c } \text{$i$zatio$} \\ \text{the w} \\ \text{obtai$} \\ \text{$5$} \end{array} } \\ & ( \mathbf p \cdot \mathbf v _ { k } ) ^ { 2 } \left [ \delta \left ( \Gamma - \frac { v _ { k } } { k } \right ) \frac { \delta } { v _ { k } ^ { 2 } } + \delta ^ { \prime } \right ] \\ & + \frac { v _ { k } } { k } p ^ { 2 } \delta - 2 ( \mathbf p \cdot \mathbf v _ { k } ) ( \mathbf p \cdot \delta _ { r } ) - \frac { 1 } { k } p _ { \theta } ^ { 2 } \partial _ { \theta } \delta _ { r \theta }. \quad ( \text{21)} \end{array}$$

Now we can calculate the ω -and p -integration in (9). For that purpose we shift ω -ω p → ω and p -k → p . This produces an exponential factor iω t p which we expand up to second order as well

$$\omega _ { p } = \omega _ { k } + ( { \mathbf p } - { \mathbf k } ) \cdot { \mathbf v } _ { k } + ( p - k ) ^ { 2 } \frac { \Gamma } { 2 }. \quad \ ( 2 2 )$$

The factor for the ω integration reads then

$$\exp \left \{ i \left ( \frac { \mathbf p \cdot \delta _ { r } } { k v _ { k } } - i \delta ^ { \prime } \mathbf p \cdot \mathbf v _ { k } + i \mathbf p \cdot \delta ^ { \prime } _ { r } - t - i \delta } \right ) \omega - \frac { \sigma _ { t } ^ { 2 } - \delta ^ { \prime } } { 2 } \omega ^ { 2 } \right \} \quad \text{proced} \quad \text{be two} \quad \text{in the}$$

which Gaussian integral is readily integrated. Separating the time-dependence from the p -dependence we obtain

$$\Psi _ { \text{out} } = & \quad f \left ( \frac { \mathbf k } { k }, \omega _ { k } \right ) \frac { \mathrm e ^ { - i \omega _ { k } t - \frac { ( t + i \delta ) ^ { 2 } } { 2 ( \sigma _ { t } ^ { 2 } - \delta ^ { \prime } ) } } } { \sqrt { 2 \pi ( \sigma _ { t } ^ { 2 } - \delta ^ { \prime } ) } } & \text{ with} \\ & \quad \times \int \frac { d ^ { 3 } p } { ( 2 \pi ) ^ { 3 } } \mathrm e ^ { i \mathbf p \cdot \tilde { r } + ( \mathbf p \cdot \nu _ { k } ) ( \mathbf p \cdot \mathbf c ) + \frac { ( \mathbf p \cdot \mathbf d ) ^ { 2 } } { 2 ( \sigma _ { t } ^ { 2 } - \delta ^ { \prime } ) } + b p ^ { 2 } }. ( 2 4 ) & \overline { \mathfrak { r } }$$

Here the abbreviations are

$$\i d \quad & \text{Here the abbreviations are} \\ b \, & = \, - \frac { 1 } { 2 \sigma ^ { 2 } } + \frac { v _ { k } } { 2 k } \delta - \frac { \partial _ { \theta } \delta _ { r \theta } } { 2 k } - i \frac { \Gamma } { 2 } t \\ c \, & = \, \frac { v _ { k } } { 2 } \left [ \delta ^ { \prime } \left ( 1 + \frac { \delta ^ { \prime } } { \bar { \sigma } _ { t } ^ { 2 } } \right ) + \delta \left ( \frac { \Gamma } { v _ { k } ^ { 2 } } - \frac { 1 } { v _ { k } } \right ) + \frac { \partial _ { \theta } \delta r _ { \theta } } { k v _ { k } ^ { 2 } } \right ] \\ & \quad - \frac { \delta ^ { \prime } \delta _ { r } } { k v _ { k } \bar { \sigma } _ { t } ^ { 2 } } - \delta _ { r } ^ { \prime } \left ( 1 + \frac { \delta ^ { \prime } } { \bar { \sigma } _ { t } ^ { 2 } } \right ) \\ d \, & = \, \delta _ { r } ^ { \prime } + \frac { \delta _ { r } } { k v _ { k } } \\ \tilde { r } \, & = \, \bar { r } - \nu _ { k } t - i \nu _ { k } \delta + i \delta _ { r } - ( t + i \delta ) \frac { \nu _ { k } \delta ^ { \prime } - \delta _ { r } ^ { \prime } } { \bar { \sigma } _ { t } ^ { 2 } }. \quad ( 2 5 ) \\ & \quad \text{We see that the width of the time-dependent pulse be-} \quad \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \end{$$

We see that the width of the time-dependent pulse becomes modified by δ ′ . Working it out for real and imaginary parts according to (2) it becomes

$$\text{such} \quad & \frac { e ^ { - i \omega _ { k } t - \frac { ( t + \delta ) ^ { 2 } } { 2 ( \sigma _ { t } ^ { 2 } - \delta ^ { \prime } ) } } } { \sqrt { 2 \pi ( \sigma _ { t } ^ { 2 } - \delta ^ { \prime } ) } } = \frac { e ^ { - i ( \omega _ { k } - \frac { \omega } { \partial _ { t } ^ { 2 } } ) t - \frac { ( t - \Delta ) ^ { 2 } } { 2 \bar { \sigma } _ { t } ^ { 2 } } + i \frac { \Delta ^ { \prime } } { 2 \sigma _ { t } ^ { 2 } } } } { \sqrt { 2 \pi \bar { \sigma } _ { t } ^ { 2 } } } + o \left ( \frac { \delta ^ { 2 } } { \sigma _ { t } ^ { 2 } }, \frac { \delta ^ { \prime } } { \sigma _ { t } ^ { 2 } } \right ).$$

Besides an overall complex phase shift due to renormalization we see that the time evolution is delayed by ∆, the width is modified by ¯ σ 2 t = σ 2 t -ϕ ′ , and the frequency obtains a shift ϕ/σ ¯ 2 t which establishes the results 1,4, and 5.

To calculate (24) we can facilitate the algebra restricting to linear orders in the shifts in the sense of (26). This omits the ( pd ) 2 term in (24) and we have

$$\int & \frac { d ^ { 3 } p } { ( 2 \pi ) ^ { 3 } } e ^ { i \mathfrak { p } \cdot \tilde { r } + ( \mathfrak { p } \cdot \mathfrak { A } \cdot \mathfrak { p } ) } = \frac { e ^ { \frac { 1 } { 4 } \tilde { r } \cdot \mathfrak { A } ^ { - 1 } \cdot \tilde { r } } } { ( 2 \pi ) ^ { 3 / 2 } \sqrt { 2 | \mathbf A | } } \quad \ ( 2 7 ) \\ \text{with the matrix}$$

with the matrix

$$\mathbf A = b \mathbf I + \mathbf v _ { k } \otimes \mathbf c.$$

The inverse and the determinant can be found

$$\mathbf A ^ { - 1 } \ & = \ \frac { 1 } { b } - \frac { \mathbf v _ { k } \otimes \mathbf c } { b ( b + \mathbf v _ { k } \cdot \mathbf c ) } \\ | \mathbf A | \ & = \ b ^ { 2 } ( b + \mathbf c \cdot \mathbf v _ { k } ) \\ \cdot \quad \cdot \quad \cdot \quad \frac { \mathbf v _ { k } } { \mathbf v _ { k } } \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdot \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdclots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \ddots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \ddots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \vdots \quad \end{cases}$$

proved by inspection. The scalar products with (25) can be worked out somewhat tediously and with linear orders in the shift according to (26) we can finally write

$$\text{rating} \quad & \Psi _ { \text{out} } = \quad f \left ( \frac { \text{k} } { k }, \omega _ { k } \right ) \frac { 1 } { \sqrt { 2 \pi ( \sigma _ { t } ^ { 2 } - \varphi ^ { \prime } ) } } \text{e} ^ { - i ( \omega _ { K } - \frac { \varphi } { \sigma _ { t } ^ { 2 } } ) t - \frac { ( t - \Delta ) ^ { 2 } } { 2 ( \sigma _ { t } ^ { 2 } - \varphi ^ { \prime } ) } } \\ & \quad \times \frac { \bar { \sigma } ^ { 3 } } { ( 2 \pi ) ^ { 3 / 2 } } \text{e} ^ { i \mathbf K \cdot \bar { r } - \frac { \sigma ^ { 2 } } { 2 } ( \bar { r } - \Delta _ { r } - \mathbf v _ { k } ( t - \Delta ) ) ^ { 2 } } \quad ( 3 0 )$$

with the modified momentum width

$$\frac { 1 } { \bar { \sigma } ^ { 2 } } = \frac { 1 } { \sigma ^ { 2 } } + \frac { \varphi _ { \theta } - v _ { k } \varphi } { k } + i \left ( \frac { \Delta _ { \theta } - v _ { k } \Delta } { k } + \Gamma t \right ). \ \ ( 3 1 )$$

The latter one can be separated into real and imaginary parts in the same way as done for the time width. The second-order dispersion Γ of (22) leads to a broadening of the pulse with increasing time as it is well known. This width of the spatial wave-packet becomes modified by real and imaginary parts of the scattering shifts which yields the effect 6. As new effects 3 and 5 we see that the momentum and frequency are shifted according to

$$\text{K} \ & = \text{ k} - \sigma ^ { 2 } ( \varphi _ { r } - \varphi \mathbf v _ { k } ) \quad & \text{den} \\ \omega _ { K } \ & = \text{ } \omega _ { k } - \sigma ^ { 2 } ( \varphi _ { r } - \varphi \mathbf v _ { k } ) \cdot \mathbf v _ { k }. \quad & ( 3 2 ) \\. \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot$$

That the latter one appears exactly as the first-order expansion of the frequency in terms of the shifted momentum is a check of internal consistency of the expansion. As an additional justification, one sees that the delay time appears consistently in the time of the wavepacket and in the traveling spatial wavepacket. So far we have derived the modification of wave packets valid for an arbitrary dispersion w k = w k ( ). In the next chapter we will calculate these shifts for light scattering w = ck and a model of scattering on a dielectric material.

## III. MODEL OF LIGHT SCATTERING ON DIELECTRIC MATERIALS

We consider a dielectric material which is characterized by a longitudinal /epsilon1 l ( ω, k ) = /epsilon1 z and transverse /epsilon1 t ( ω, k ) = /epsilon1 x = /epsilon1 y dielectric function such that the dielectric tensor reads

$$\varepsilon = \alpha _ { \phi _ { 1 } } ^ { T } \alpha _ { \theta _ { 1 } } ^ { T } \text{diag} \{ \epsilon _ { x }, \epsilon _ { y }, \epsilon _ { z } \} \alpha _ { \theta _ { 1 } } \alpha _ { \phi _ { 1 } } \text{ \quad \ \ ( 3 3 ) \quad,$$

where θ 1 is the rotation around y -axes and φ 1 the rotation around the z -axes out of the symmetry axes of the crystal. We keep for general case three different /epsilon1 for the different directions. A scattering of electromagnetic waves can be described by the scattering amplitude [33] up to momentum and energy-independent constants

$$f \sim \mathbf e \cdot \varepsilon \cdot \mathbf e _ { a } \, \quad \, \quad \, ( 3 4 ) \quad =$$

Here the direction of the incoming wave with wavevector k is chosen as e = (sin θ, 0 cos , θ ) and the direction e a of the scattered wave has an opposite z -direction. We will calculate any further shifts by derivatives of this expression and will average about the Gaussian profile of the beam. This integration is to be performed about the upper half plane ( z &gt; 0) as possible incoming directions. We can average about the total space if we extend the expression (34) by e a → e such that finally we will calculate

$$\delta = \langle \partial \ln ( \mathbf e \cdot \varepsilon \cdot \mathbf e ) \rangle _ { G a u \i } \quad \quad ( 3 5 ) \quad \mathbf F \quad \mathbf e \quad \mathbf e$$

First we consider the spatial shifts (3) by performing the angular derivatives (4) with (34). We consider dielectric functions /epsilon1 l,t ( ω , k k ) = /epsilon1 l,t ( ω k ) with ω k = ck which means that the derivatives on the dielectric functions itself vanishes

$$\overset { \text{$7$ of} } { \text{his} } \quad \frac { 1 } { k } \partial _ { \theta, \phi \epsilon } = ( \text{e} \partial _ { k } - \partial _ { k } ) \epsilon ( \omega _ { k } ) = c \partial _ { \omega } \epsilon \left ( \frac { \text{k} } { k } - \text{e} \right ) = 0. \ ( 3 6 ) \\ \underset { \text{th} } { \text{th} } \quad \text{$2$} \quad \epsilon \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots$$

Therefore only the derivatives with respect to the direction factors e in (34) matters. This would be different for dielectric functions with explicit momentum dependencies. Since k = k e we use

$$\omega _ { 1 } \right ) & \quad \left ( \mathbf e _ { \theta } \partial _ { \theta } + \frac { \mathbf e _ { \phi } } { \sin \theta } \partial _ { \phi } \right ) \mathbf b \cdot \mathbf e \, = \, ( \mathbf e _ { \theta } \cdot \mathbf b ) \mathbf e _ { \theta } + ( \mathbf e _ { \phi } \cdot \mathbf b ) \mathbf e _ { \phi } \\ \mathbf e r \\ \rangle & \quad \text{with $b$} \, - \, \mathbf e _ { 1 } \sim \mathbf e _ { 2 } \, \text{that} \, \mathbf a \, \text{matrix} \, \text{domain} \, \text{and} \, \text{and} \, \text{a}$$

with b = e · ε to obtain a material-dependent and a material-independent part

$$\delta = - \frac { k \cdot \varepsilon + \varepsilon \cdot k } { k \cdot \varepsilon \cdot k } + 2 \frac { k } { k ^ { 2 } } = \delta _ { 1 } + \delta _ { 2 }. \quad \ ( 3 8 ) \\ \intertext { since we consider a Gaussian beam around the mean mo-}$$

Since we consider a Gaussian beam around the mean momentum

$$\mathbf k _ { 0 } = k _ { 0 } ( \sin \vartheta _ { 0 }, 0, \cos \vartheta _ { 0 } )$$

and a momentum spreading of σ we have to average about this Gaussian beam. In the second part of (38) we scale the momenta x 0 = k / 0 √ 2 σ and ¯ = k k/ √ 2 σ and obtain

$$\langle$$

$$\text{level} \, \text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text	$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text
$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text#\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{	$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\font{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\font$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text.$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text}$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text<$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text'$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$
$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\test{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\test$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text1$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$$$$

Here we used the complex error function erfi( x ) = -i erf( ix ) and the Dawson integral D x ( ) = exp( -x 2 ) ∫ ∞ 0 dy exp( y 2 ). This second part of the spatial shift is independent on the material and purely real. It follows the beam direction (39) and gives a contribution to the Goos-H¨ anchen effect exclusively by the beam shape and no contribution to the Imbert-Fedorov effect. As a check we see that the limit of vanishing momentum width σ → 0 appears correctly which means x 0 →∞ corresponding to averaging about the δ ( k -k 0 ) function.

For the first and material-dependent part of (38) we perform a partial integration

$$\stackrel { \log } { _ { \mathbb { m } { h } } } } _ { \mathbb { m } { h } } ^ { \log } \int \frac { d ^ { 3 } k } { ( 2 \pi \sigma ^ { 2 } ) ^ { 3 / 2 } } e ^ { - \frac { ( k - k _ { 0 } ) ^ { 2 } } { 2 \sigma ^ { 2 } } } \ln ( k \cdot \varepsilon \cdot k ) \quad ( 4 1 )$$

and understand the derivative not concerning ε according to (36). For the logarithm we use the trick ln c = ∫ ∞ 0 dξ [exp( -ξ ) -exp( -cξ )] /ξ such that we obtain

$$\langle \delta _ { 1 } \rangle = \partial _ { \mathbf k _ { 0 } } e ^ { - \frac { k _ { 0 } ^ { 2 } } { 2 \sigma ^ { 2 } } } \underset { 0 } { \overset { \infty } { \int } } \frac { d \xi } { \xi } \int \frac { d ^ { 3 } k } { ( 2 \pi \sigma ^ { 2 } ) ^ { 3 / 2 } } e ^ { \mathbf k \cdot \mathbf p - \frac { 1 } { 2 } \mathbf k \cdot \mathbf A \cdot \mathbf k - \frac { k _ { 0 } ^ { 2 } } { 2 \sigma ^ { 2 } } } ( 4 2 ) \quad \varphi _ { r } ^ { \mathbf c }$$

with p = k 0 /σ 2 , A = bI +2 ξε and b = 1 /σ 2 . We further approximate ε k ( ) ≈ ε k ( 0 ) and can apply the formula (27). For that purpose we need

$$\mathbf A ^ { - 1 } \, = \, \frac { 1 } { b } I$$

$$& + \alpha _ { \phi _ { 1 } } ^ { T } \alpha _ { \theta _ { 1 } } ^ { T } \text{diag} \{ \frac { - 2 \xi \epsilon _ { x } } { 2 \xi \epsilon _ { x } + b }, \frac { - 2 \xi \epsilon _ { y } } { 2 \xi \epsilon _ { y } + b }, \frac { - 2 \xi \epsilon _ { z } } { 2 \xi \epsilon _ { z } + b } \} \alpha _ { \theta _ { 1 } } \alpha _ { \phi _ { 1 } } \\ | A | \ = \ & ( 2 \xi \epsilon _ { x } + b ) ( 2 \xi \epsilon _ { y } + b ) ( 2 \xi \epsilon _ { z } + b ) & ( 4 3 ) \quad \text{FIG. 1.} \\ \text{which is easily seen from } ( 3 3 ) \ \text{With the substitution } \quad \text{transv}$$

which is easily seen from (33). With the substitution 2 σ ξ 2 = z and the abbreviations

$$a & = \sin \vartheta _ { 0 } \cos \theta _ { 1 } \cos \phi _ { 1 } - \cos \vartheta _ { 0 } \sin \theta _ { 1 } \quad \cdot \\ b & = - \sin \vartheta _ { 0 } \sin \theta _ { 1 } \\ c & = \sin \vartheta _ { 0 } \sin \theta _ { 1 } \cos \phi _ { 1 } + \cos \vartheta _ { 0 } \cos \theta _ { 1 } \quad \ \ ( 4 4 )$$

we obtain

$$\langle \delta _ { 1 } \rangle \ = \ \frac { 1 } { k _ { 0 } } \left [ \mathbf e _ { 1 } \partial _ { c } + \mathbf e _ { \theta _ { 1 } } \partial _ { a } + \mathbf e _ { \phi _ { 1 } } \partial _ { b } b \right ] I ( a, b, c ) \ \ ( 4 5 ) \quad - \mathbf 0.$$

with

$$\dots \dots & & - \infty _ { 0 } \\ I ( a, b, c ) = \underset { 0 } { \overset { \infty } { \int } } & e ^ { - x _ { 0 } ^ { 2 } \left [ \frac { z a ^ { 2 } \epsilon _ { x } } { 1 + \epsilon _ { x } z } + \frac { z b ^ { 2 } \epsilon _ { y } } { 1 + \epsilon _ { y } z } + \frac { z c ^ { 2 } \epsilon _ { z } } { 1 + \epsilon _ { z } z } \right ] } } \\ & \text{(4)} & \text{(4)} & \text{$-0.i$} \\ & \text{the spherical unit vectors} \, \text{e} _ { 1 }, \text{e} _ { \theta _ { 1 } }, \text{e} _ { \phi _ { 1 } }, \text{ and } x _ { 0 } = \frac { k _ { 0 } } { \sigma \sqrt { 2 } }.$$

the spherical unit vectors e 1 , e θ 1 , e φ 1 , and x 0 = k 0 σ √ 2 .

From the prefactors (44) we see already that for symmetrically oriented probes with θ 1 = 0 as well as perpendicular beams ϑ 0 = 0 the Imbert-Fedorov shift disappears.

Just for illustrative purpose let us discuss the case of homogeneous materials /epsilon1 = /epsilon1 l = /epsilon1 t . Then we have a 2 + b 2 + c 2 = 1 and the prefactor in (45) becomes

$$\mathbf e _ { 1 } c + \mathbf e _ { \theta _ { 1 } } a + \mathbf e _ { \phi _ { 1 } } b = \frac { \mathbf k _ { 0 } } { k _ { 0 } } = ( \sin \vartheta _ { 0 }, 0, \cos \vartheta _ { 0 } ) \quad ( 4 7 ) \quad \text{FIG.}$$

which means that Imbert-Fedorov effect is absent. Scaling z/epsilon1 → /epsilon1 the first part becomes independent on the material

$$\maxerra \\ \langle \delta _ { r } ^ { 1 } \rangle \ & = \ - \frac { 2 \mathbf k _ { 0 } } { k _ { 0 } ^ { 2 } } x _ { 0 } ^ { 2 } \int \frac { \infty } { z ( 1 + z ) ^ { 5 / 2 } } e ^ { - x _ { 0 } ^ { 2 } \frac { \xi } { 1 + z } } \quad \text{and} \quad \text{for h} \\ & = \ - \frac { 2 \mathbf k _ { 0 } } { k _ { 0 } ^ { 2 } } \left [ 1 - \frac { D ( x _ { 0 } ) } { x _ { 0 } } \right ] \quad \text{(48)} \quad \text{Goo} \\ \text{which exactly compensates the second part } ( 40 ) \text{ such that} \quad \text{(\begin{matrix} \Delta ^ { \mathcal { C } } \\ \Delta ^ { 1 } \end{matrix} }$$

which exactly compensates the second part (40) such that the total spatial shift vanishes

$$\langle \delta _ { r } \rangle = 0 \quad \quad \ \ ( 4 9 ) \quad w i$$

FIG. 1. The Goos-H¨nchen shift versus width of wavenumbers a for a plasma frequency /planckover2pi1 ω p = 5 8eV of Au and the ratio of . transverse to longitudinal relaxation times. The transverse relaxation time was chosen as τ = 20fs.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00375v4-figure-9bd89ee7-0f01-437b-aafb-6d5df86b0266.jpg)

FIG. 2. The Wigner delay time (above) and the real part (below) for at width of σ = /planckover2pi1 / 30nm and various relaxation times.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00375v4-figure-803a89c0-5fb5-4a1d-ad38-abb2358899c9.jpg)

and no Goos-H¨ ahnchen or Imbert-Fedorov shift appears for homogeneous dielectric materials.

For demonstration we assume ¯ = ( /epsilon1 /epsilon1 l -1) / /epsilon1 ( t -1). Then we can directly calculated the Gaussian-averaged Goos-H¨ anchen (5) and Imbert-Fedorov shifts (6)

$$\mathfrak { h } \, \text{th that} \quad \begin{pmatrix} \Delta ^ { \text{GH} } \\ \Delta ^ { \text{IF} } \end{pmatrix} = \left \langle \begin{pmatrix} \sin \theta \cos ^ { 2 } \theta \cos \phi - \sin ^ { 3 } \theta \\ \sin \theta \cos ^ { 2 } \theta \sin \phi \end{pmatrix} \frac { 2 ( \bar { \epsilon } - 1 ) } { k ( \sin ^ { 2 } \theta + \bar { \epsilon } \cos ^ { 2 } \theta ) } \right \rangle$$

which shows that the Imbert-Fedorov shift vanishes due

to averaging about sin φ and the Goos-H¨nchen shift bea comes with the notation of (40), x = cos θ , and t = x 2

$$\Delta ^ { \text{GH} } = - \frac { \sqrt { 2 } ( \bar { \epsilon } - 1 ) } { \sqrt { \pi } } e ^ { - x _ { 0 } ^ { 2 } } \int _ { 0 } ^ { \infty } \tilde { d } \tilde { k } e ^ { - \tilde { k } ^ { 2 } } \int _ { - 1 } ^ { 1 } d x \frac { ( 1 - x ^ { 2 } ) ^ { 3 / 2 } e ^ { x ^ { 2 } x _ { 0 } ^ { 2 } } } { 1 + ( \bar { \epsilon } - 1 ) x ^ { 2 } } \quad \text{(} \\ = - \frac { 2 ( \bar { \epsilon } - 1 ) } { k _ { 0 } } x _ { 0 } ^ { 2 } e ^ { - x _ { 0 } ^ { 2 } } \int _ { 0 } ^ { 1 } d t \frac { ( 1 - t ) ^ { 3 / 2 } e ^ { t x _ { 0 } ^ { 2 } } } { 1 + ( \bar { \epsilon } - 1 ) t }. \quad \text{(} \\ = \text{(} \\ \text{deriva} }$$

$$x _ { 0 } ^ { 2 } \int _ { 0 } ^ { 1 } d t \frac { ( 1 - t ) ^ { 3 / 2 } \mathrm e ^ { t \colon } } { 1 + ( \bar { \epsilon } - 1 ) }$$

FIG. 3. The ratio of the reduced temporal width effect 4 vs. wavelength width (above) for λ = 770nm and vs. wavelength (below) for ∆ λ = 100nm.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00375v4-figure-520c0642-27a0-40d6-8c6f-edca1063dd29.jpg)

We choose as an example a simple Drude model in the following

$$\epsilon _ { j } = 1 - \frac { \omega _ { p } ^ { 2 } } { \omega ( \omega + \frac { i } { \tau _ { j } } ) } \quad \quad ( 5 2 ) \quad R e$$

with a plasma frequency for Au of /planckover2pi1 ω p = 5 8eV [34] . and different relaxation times in longitudinal j = l and transverse j = t direction. The result is plotted in figure 1. We assume different relaxation times in longitudinal and transverse direction. It vanishes quadratically with x 0 = k /σ 0 √ 2 and in the homogeneous limit τ t = τ l . For τ l &gt; τ t the sign of the shift is changed.

Next we discuss the time shifts (2) which are calculated analogously

$$\langle \delta \rangle = \varphi + i \Delta = \partial _ { \omega } I ( a, b, c ) = \langle \partial _ { \omega } \ln \mathbf e \cdot \epsilon I \cdot \mathbf e \rangle = \frac { \partial _ { \omega } \epsilon } { \epsilon } \quad \begin{array} { c } \mathbf i \, \text{in} \, p \\ \text{quenc} \\ \text{enhar} \end{array}$$

where the last expression is valid for homogeneous materials. It shows that the geometry of the Gaussian beam drops out due to the logarithmic derivative of scattering amplitude (34) for homogeneous dielectric materials. Therefore measuring the time delay ∆ (effect 1) and the shrinking of the temporal width of the pulse ∂ ω ϕ (effect 4) one has the possibility to access the dielectric function of the material directly.

The behaviour of the time shifts an be seen by the derivatives of the dielectric function (52)

$$\langle \vec { \ }, \quad \text{$i$ ne behavior of the time snirts an be seen by $ne$} \\ \varphi = \text{Re} \, \frac { \partial _ { \omega } \epsilon } { \epsilon } = \frac { \omega _ { p } ^ { 2 } \tau ^ { 2 } \left ( 2 \tau ^ { 2 } \omega ^ { 2 } \left ( \omega ^ { 2 } - \omega _ { p } ^ { 2 } \right ) - \omega _ { p } ^ { 2 } \right ) } { \omega \left ( \tau ^ { 2 } \omega ^ { 2 } + 1 \right ) \left ( \tau ^ { 2 } \left ( \omega ^ { 2 } - \omega _ { p } ^ { 2 } \right ) ^ { 2 } + \omega ^ { 2 } \right ) } \\ \approx \frac { 2 \omega _ { p } ^ { 2 } } { \omega ^ { 3 } - \omega \omega _ { p } ^ { 2 } } + \omega _ { p } ^ { 2 } \frac { 4 \omega ^ { 4 } - 3 \omega ^ { 2 } \omega _ { p } ^ { 2 } + \omega _ { p } ^ { 4 } } { \tau ^ { 2 } \omega ^ { 3 } \left ( \omega _ { p } ^ { 2 } - \omega _ { 2 } \right ) ^ { 3 } } + o \left ( \frac { 1 } { \tau } \right ) ^ { 3 } \\ \Delta = \text{Im} \, \frac { \partial _ { \omega } \epsilon } { \epsilon } = \frac { \omega _ { p } ^ { 2 } \tau ^ { 3 } \left ( \omega _ { p } ^ { 2 } - 3 \omega ^ { 2 } \right ) - \omega _ { p } ^ { 2 } \tau } { \left ( \tau ^ { 2 } \omega ^ { 2 } + 1 \right ) \left ( \tau ^ { 2 } \left ( \omega ^ { 2 } - \omega _ { p } ^ { 2 } \right ) ^ { 2 } + \omega ^ { 2 } \right ) } \\ \approx \omega _ { p } ^ { 2 } \frac { \omega _ { p } ^ { 2 } - 3 \omega ^ { 2 } } { \tau \omega ^ { 2 } \left ( \omega ^ { 2 } - \omega _ { p } ^ { 2 } \right ) ^ { 2 } } + o \left ( \frac { 1 } { \tau } \right ) ^ { 3 } \\ \text{ presented in figure 2. One sees that Wigner's delay time} \\ \text{changes the sign for} \\ \text{$i$,} 2 / \text{ $i$ } \, \ \end{cases}$$

presented in figure 2. One sees that Wigner's delay time changes the sign for

$$\omega _ { \Delta } ^ { 2 } & = \frac { \omega _ { p } ^ { 2 } } { 3 } \left ( \omega _ { p } \tau - \frac { 1 } { \omega _ { p } \tau } \right ) \\ \omega _ { \Delta } \quad \omega _ { \Delta } \quad \dots \quad \dots \quad \dots$$

which means for τ &gt; 1. The real part changes sign at

$$\omega _ { \varphi } ^ { 2 } = \frac { \omega _ { p } ^ { 2 } } { 2 } \left [ 1 + \sqrt { \left ( 1 + \frac { 2 } { \omega _ { p } ^ { 2 } \tau ^ { 2 } } \right ) } \right ]$$

for any τ . For large relaxation times the shifts develop a pole at the plasma frequency.

The shrinking of the temporal width is given by the frequency derivative of the real part ∂ ω ϕ according to effect 4 and is shown in figure 3. It leads to a maximum at the plasma frequency for large relaxation times. Further it shows a zero at certain wavelength or frequencies

$$\begin{array} { c } \text{$i$ shows a zero at certain wavelengin or frequencies} \\ ( 5 2 ) & \text{Re} \, \partial _ { \omega } \left ( \frac { \partial _ { \omega } \epsilon } { \epsilon } \right ) \approx \omega _ { p } ^ { 2 } \frac { 2 \omega _ { p } ^ { 2 } - 6 \omega ^ { 2 } } { \omega ^ { 2 } \left ( \omega ^ { 2 } - \omega _ { p } ^ { 2 } \right ) ^ { 2 } } \\ \text{$V$ [34]$} \\ l \, \text{and} & + \frac { 2 0 \omega ^ { 6 } - 1 7 \omega ^ { 4 } \omega _ { p } ^ { 2 } + 1 2 \omega ^ { 2 } \omega _ { p } ^ { 4 } - 3 \omega _ { p } ^ { 6 } } { \omega _ { p } ^ { 2 } \tau ^ { 2 } \omega ^ { 4 } \left ( \omega ^ { 2 } - \omega _ { p } ^ { 2 } \right ) ^ { 4 } }. + o \left ( \frac { 1 } { \tau } \right ) ^ { 3 }. \\ \text{in fig-} \\ \text{ngitu-} \\ \text{tically} & \text{One sees that the real part $\varphi(\omega)$ has an extreme value} \end{array}$$

One sees that the real part ϕ ω ( ) has an extreme value at ω = ω / p √ 3 corresponding to a zero at the shift of temporal width ∂ ω ϕ for large τ . We see in figure 3 that the effect 4 of the shrinking of the temporal width of the wave packet is dependent on the wavelength and width. Though is seems to be very small it should be noted that it is proportional to the square of the inverse plasma frequency. A reduced plasma frequency leads to a quadratic enhancement of effect 4.

## IV. SUMMARY

By expanding the scattered wavepacket up to second order, we obtain six different effects that can be expressed by scattering shifts. These shifts are derivatives of the scattering amplitude with respect to energy and momentum. Their imaginary parts, as derivatives of the scattering phase, determine the Wigner delay time and the spatial displacements known as Goos-H¨nchen or Imberta Fedorov effects in optics. The real parts, as derivatives of the modulus of the scattering amplitude, yield energy and momentum shifts as well as a shrinking of the pulse width.

The shifts are calculated analytically for a materials described by longitudinal and transverse dielectric func-

- [1] F. Goos and H. H¨nchen, Ein neuer und fundamentaler a Versuch zur Totalreflexion, Annalen der Physik 436 , 333 (1947).
- [2] F. I. Fedorov, To the theory of total reflection, Journal of Optics 15 , 014002 (2013).
- [3] K. Y. Bliokh and A. Aiello, Goos-H¨nchen and Imberta Fedorov beam shifts: An overview, J. Opt. 15 , 014001 (2013).
- [4] K. Y. Bliokh and F. Nori, Transverse and longitudinal angular momenta of light, Physics Reports 592 , 1 (2015), transverse and longitudinal angular momenta of light.
- [5] A. Aiello, Goos-H¨nchen and Imbert-Fedorov shifts: a a novel perspective, New Journal of Physics 14 , 013058 (2012).
- [6] M. R. Dennis and J. B. G¨tte, Beam shifts for pairs of o plane waves, Journal of Optics 15 , 014015 (2013).
- [7] C. Duval, P. A. Horv´thy, and P. M. Zhang, Transverse a shifts in paraxial spinoptics, Journal of Optics 15 , 014005 (2013).
- [8] W. L¨ offler, A. Aiello, and J. P. Woerdman, Spatial coherence and optical beam shifts, Phys. Rev. Lett. 109 , 213901 (2012).
- [9] W. L¨ offler, N. Hermosa, A. Aiello, and J. P. Woerdman, Total internal reflection of orbital angular momentum beams, Journal of Optics 15 , 014012 (2013).
- [10] M. Merano, N. Hermosa, J. P. Woerdman, and A. Aiello, How orbital angular momentum affects beam shifts in optical reflection, Phys. Rev. A 82 , 023817 (2010).
- [11] Y. Qin, Y. Li, J. Ren, Q. Wen, J. Zhang, Y.-F. Xiao, H. Yang, and Q. Gong, Spin separations of light at the air-glass interface for femtosecond laser pulses, Journal of Optics 15 , 014006 (2013).
- [12] Q. Liu, W. Zhen, M. Gao, and D. Deng, Goos-H¨nchen a and Imbert-Fedorov shifts for the rotating elliptical Gaussian beams, Results in Physics 18 , 103297 (2020).
- [13] M. Ornigotti, Goos-H¨nchen and Imbert-Fedorov shifts a for airy beams, Opt. Lett. 43 , 1411 (2018).
- [14] J. Zhou, Y. Liu, Y. Ke, H. Luo, and S. Wen, Generation of airy vortex and airy vector beams based on the modulation of dynamic and geometric phases, Opt. Lett. 40 , 3193 (2015).
- [15] C. Liu, J. Liu, L. Niu, and et al, Terahertz circular airy vortex beams, Sci Rep 7 , 3891 (2017).

tion. It is found that the Imbert-Fedorov effect is absent for homogeneous materials or perfectly aligned crystal axes to the scattering plane. For homogeneous material the Wigner time delay and the shrinking of the temporal width as the frequency derivative of the Wigner delay time can directly access the dielectric function of the material. This could be a suggestion for experiments.

## ACKNOWLEDGMENTS

Evgeny Gurevich is thanked for stimulating discussions and his convincing arguments that these effects can be observed.

- [16] H.-Y. Song, Z.-X. Chen, Y.-B. Li, S.-P. Hao, Q. Zhang, S. Zhou, S. fang Fu, and X.-Z. Wang, Large spatial shifts of a reflected airy beam on the surface of hyperbolic crystals, J. Opt. Soc. Am. B 40 , 1240 (2023).
- [17] X. Yang, T. Qu, Z. Wu, H. Li, L. Bai, L. Gong, and Z. Li, Characteristics of an airy beam at a dielectric interface, Optics &amp; Laser Technology 156 , 108607 (2022).
- [18] S. Grosche, A. Szameit, and M. Ornigotti, Spatial GoosH¨ anchen shift in photonic graphene, Phys. Rev. A 94 , 063831 (2016).
- [19] N. A. F. Zambale, J. L. B. Sagisi, and N. P. Hermosa, Goos-H¨ anchen shifts due to graphene when intraband conductivity dominates, Optics Communications 433 , 25 (2019).
- [20] X. Guo, L. Zhang, X. Zhang, and B. Shen, Deflection of a reflected intense spatiotemporal optical vortex beam, Opt. Lett. 48 , 1610 (2023).
- [21] S.-Y. Yang and C.-F. Li, Properties of the barycenter of a diffraction-free light beam, Journal of Optics 15 , 014016 (2013).
- [22] R. F. Gragg, The total reflection of a compact wave group: Long-range transmission in a waveguide, American Journal of Physics 56 , 1092 (1988).
- [23] W. Zhen and D. Deng, Goos-H¨ anchen and Imbert-Fedorov shifts in temporally dispersive attenuative materials, Journal of Physics D: Applied Physics 53 , 255104 (2020).
- [24] O. J. S. Santana and L. E. E. de Araujo, Goos-H¨nchen a and Imbert-Federov shifts of vortex beams near critical incidence, J. Opt. Soc. Am. B 38 , 300 (2021).
- [25] M. Mazanov and K. Y. Bliokh, Wigner time delays and Goos-H¨ anchen shifts of 2d quantum vortices scattered by potential barriers, Journal of Physics A: Mathematical and Theoretical 55 , 404005 (2022).
- [26] M. Ornigotti and A. Aiello, Goos-H¨ anchen and Imbert-Fedorov shifts for bounded wavepackets of light, Journal of Optics 15 , 014004 (2013).
- [27] A. Nieminen, A. Marini, and M. Ornigotti, Goos-H¨ anchen and Imbert-Fedorov shifts for epsilonnear-zero materials, Journal of Optics 22 , 035601 (2020).
- [28] F. T¨ oppel, M. Ornigotti, and A. Aiello, Goos-H¨nchen a and Imbert-Fedorov shifts from a quantum-mechanical

perspective, New Journal of Physics 15 , 113059 (2013).

- [29] V. ˇ piˇka, P. Lipavsk´, and K. Morawetz, Nonlocal corS c y rections to boltzmann equation for dense fermi systems, Phys. Lett. A 240 , 160 (1998).
- [30] K. Morawetz, P. Lipavsk´, and V. ˇ piˇka, Retarded very S c sus time-nonlocal quantum kinetic equations, Ann. of Phys. 294 , 135 (2001).
- [31] K. Morawetz, Interacting systems far from equilibrium -quantum kinetic theory (Oxford University Press, Oxford, 2017).
- [32] K. Morawetz, Nonequilibrium thermodynamics with binary quantum correlations, Phys. Rev. E 96 , 032106
- (2017), critics: P. Lipavsky, Phys. Rev. E 97, 066103 (2018) and reply: K. Morawetz, arXiv:1806.11324.
- [33] B. J. Berne and R. Pecora, Dynamic light scattering (Dover, Mineola, New York, 2000).
- [34] D. Langley, R. A. C. Jr., L. A. Starman, and S. Rogers, Optical metamaterials for photonics applications, in Adaptive Coded Aperture Imaging, Non-Imaging, and Unconventional Imaging Sensor Systems , Vol. 7468, edited by D. P. Casasent, S. Rogers, J. J. Dolne, T. J. Karr, and V. L. Gamiz, International Society for Optics and Photonics (SPIE, 2009) p. 74680H.
K. Morawetz
Second-order nonlocal shifts of scattered wave-packets: What can be measured by Goos-Hänchen and Imbert-Fedorov effects ?
The scattering of wavepackets with arbitrary energy dispersion on surfaces
has been analyzed. Expanding up to second order in scattering shifts, it is
found that besides the known Goos-H\"anchen or Imbert-Fedorov spatial offset,
as well as the Wigner delay time, new momentum and frequency shifts appear.
Furthermore, the width of the scattered wave packet becomes modified as well,
which can lead to a shrinking of pulses by multiple scattering. For a model of
dielectric material characterized by a longitudinal and transverse dielectric
function the shifts are calculated analytically. From the Goos-H\"anchen and
Imbert-Fedorov shifts one can access the longitudinal and transversal
dielectric function. Perfectly aligned crystal symmetry axes with respect to
scattering beam shows no Imbert-Fedorov effect. It is found that the
Goos-H\"anchen and Imbert-Fedorov effect are absent for homogeneous materials.
Oppositely it is found that the Wigner delay time and the shrinking of the
temporal pulse width allows to access the dielectric function independent on
the beam geometry.