Observation of a phase space horizon with surface gravity water waves

In 1974, Stephen Hawking predicted that quantum effects in the proximity of a black hole lead to the emission of particles and black hole evaporation. At the very heart of this process lies a logarithmic phase singularity which leads to the Bose-Einstein statistics of Hawking radiation. An identical singularity appears in the elementary quantum system of the inverted harmonic oscillator. In this Letter we report the observation of the onset of this logarithmic phase singularity emerging at a horizon in phase space and giving rise to a Fermi-Dirac distribution. For this purpose, we utilize surface gravity water waves and freely propagate an appropriately tailored energy wave function of the inverted harmonic oscillator to reveal the phase space horizon and the intrinsic singularities. Due to the presence of an amplitude singularity in this system, the analogous quantities display a Fermi-Dirac rather than a Bose-Einstein distribution.

In 1974, Stephen Hawking predicted that quantum effects in the proximity of a black hole lead to the emission of particles and black hole evaporation.At the very heart of this process lies a logarithmic phase singularity which leads to the Bose-Einstein statistics of Hawking radiation.An identical singularity appears in the elementary quantum system of the inverted harmonic oscillator.In this Letter we report the observation of the onset of this logarithmic phase singularity emerging at a horizon in phase space and giving rise to a Fermi-Dirac distribution.For this purpose, we utilize surface gravity water waves and freely propagate an appropriately tailored energy wave function of the inverted harmonic oscillator to reveal the phase space horizon and the intrinsic singularities.Due to the presence of an amplitude singularity in this system, the analogous quantities display a Fermi-Dirac rather than a Bose-Einstein distribution.
When a massive star collapses, a black hole 1 is born.During the collapse, the matter is compressed to an infinitesimal small volume of infinite density leading to a singularity in spacetime, surrounded by a domain where gravity is strong enough to capture light.This area is bounded by an event horizon dividing spacetime into two disjunct regions.
Since light cannot escape, one might think that black holes are black 1 .However, Hawking 2,3 postulated that a black hole emits radiation with a spectrum governed by the Bose-Einstein distribution and similar to that of a black body.This phenomenon is a direct consequence [4][5][6] of quantum field theory and the curvature of spacetime at the event horizon of the black hole.Essential for the so-called Hawking radiation is a logarithmic phase singularity 7,8 in the mode functions of the quantized light field.Indeed, the characteristic Bose-Einstein distribution is a consequence 8 of the Fourier transform of this phase, as was already discussed back in 1974 in the seminal paper of Hawking 2 .
In order to study quantum effects of this type, access to a black hole is not mandatory.Many analog systems such as negative-frequency waves [9][10][11][12] , Bose-Einstein condensates 13 , optical fibers 14 , and shallow water waves 15 are experimentally accessible.While the focus of these works lies on the observation of effects similar to Hawking radiation, here we are interested in the measurement of its origin 8,16 , that is, of the logarithmic phase singularity in a mode function close to an event horizon.For this purpose, we exploit an analogy between a black hole and an inverted harmonic oscillator 17,18 .The energy eigenfunctions of this system also display 19 a horizon with a logarithmic phase singularity, however, now in phase space.Since this phenomenon only takes place in specific phase space variables, we use the free time evolution of these energy wave functions to bring out the effects of this phase space horizon most clearly.This technique also allows us to show experimentally that for the inverted harmonic oscillator, the energy distribution associated with this horizon is of the Fermi-Dirac rather than Bose-Einstein type.
In our article, we use the one-dimensional inverted harmonic oscillator as a system to study black hole physics.For this purpose, we describe its inherent properties such as the phase space horizons, the logarithmic phase singularity, as well as the Fermi-Dirac distribution.We demonstrate how a free propagation of its quantum-mechanical energy eigenfunctions enables an observation of these features in position space.We present an experiment based on a classical analog of a quantum system by employing surface gravity water waves 20,21 .Our measurements confirm the predicted horizon as well as the associated singularities and allow us to extract the Fermi-Dirac distribution.

Results and discussion
Phase space horizons, logarithmic phase singularity, and Fermi-Dirac distribution We start by reviewing the classical dynamics of a particle moving in an inverted parabolic potential.As shown in Fig. 1a (bottom), for energies E below the top of the barrier, that is, E < 0, the particle is reflected, whereas for energies above, that is, E > 0, it is transmitted.When the energy is equal to the maximum of the barrier, that is, E = 0, the particle approaches the top and cannot cross over to the other domain.
This particular dynamic gives rise to four disjunct regions in phase space, which are separated by the lines p = mωx and p = −mωx, as depicted in Fig. 1a (top).Here we have introduced the position x and the momentum p of the particle, while m and ω denote its mass and the steepness of the parabolic potential, respectively.The line p = mωx separates particles coming from the left from particles coming from the right, and thus corresponds to a horizon in phase space.
However, the line p = −mωx also represents a horizon in phase space.This fact stands out most clearly in Fig. 1(b), where we present the corresponding quantum picture 17 of the energy eigenfunction of the inverted harmonic oscillator with a positive energy eigenvalue ε, by depicting the corresponding Wigner function 19,22 Here N þ and D ν denote the normalization constant and the parabolic cylinder function 23 .Moreover, the scaling parameter κ ensures that the argument of D ν is dimensionless.For more details, we refer to the Methods section.
The Wigner function W þ ϕ;ε vanishes in a half-plane of phase space separated by the horizon p = −mωx, is constant along all classical trajectories in the other half-plane and, in particular, along the horizon, while displaying a dominant positive maximum (red domain) in the neighborhood of the trajectory corresponding to the eigenvalue ε.It oscillates between positive and negative values for trajectories (yellow and blue hyperbolic bands) determined by energies larger than the eigenvalue, but decays with oscillations for trajectories (light blue and yellow domains) governed by negative energies reaching into the classically non-accessible quadrant.At the bottom, we display the probabilities in the position variable x obtained by integration of the Wigner function over the momentum p at a given coordinate.
We emphasize that the energy eigenfunction in position representation defined by Eq. (1) does not give any indication of a singular behavior.However, as shown in the Methods section, a logarithmic phase singularity and an amplitude singularity manifest 19 themselves in a specific representation of the corresponding state.
To uncover these features, we now make use of the free time evolution 24 , which translates to a sheering in phase space as shown in the transition from Fig. 1b-d.The situation in Fig. 1c depicts the moment t þ c 1=ω in time, where the Wigner function is constrained to the right side of x = 0.In this case, the position density is restricted to positive values of x only, with an inverse square-root amplitude and a logarithmic phase singularity as expressed by the wave function where φ is a phase and Θ denotes the Heaviside step function.
The square-root amplitude singularity in this probability amplitude is a consequence of integrating the constant Wigner function at x = 0 over the momentum.For more details, we refer to the Methods section.At a later time, shown in Fig. 1d, the probability density returns to the domain x < 0, and covers again the complete coordinate axis.
Next, we recall that the Bose-Einstein statistics of the Hawking radiation is a consequence 8,19 of the Fourier transform of the logarithmic phase singularity.In order to obtain its analog for the inverted harmonic oscillator, we could make use of the backward propagation of the wave function ϕ þ ε ðxÞ until the critical time t À c ¼ À1=ω.Instead, we consider here the forward propagation of the energy eigenfunction determined by the complex conjugate of Eq. ( 1), which also allows us to extract the analog distribution for the inverted harmonic oscillator.Here the As demonstrated in the Methods section, when we propagate the wave function , we obtain the Fermi-Dirac distribution 25 FðεÞ ½1 þ expð2πεÞ À1 determining the transmission and reflection coefficients T(ε) = F(−ε) and R(ε) = F(ε) of the inverted harmonic oscillator from the expression For a given energy ε, the steepness of the position distribution associated with the simple pole at x = 0 determines R and T, and thus F. Indeed, positive x-values yield R(ε) while negative ones lead to T(ε).We find the dependence of R and T, and thus of F on ε, by propagating wave functions corresponding to different energies.

Experimental observation with surface gravity water waves
In order to observe the predicted phenomena 19 of a logarithmic phase singularity and a Fermi-Dirac distribution in the reflection and transmission coefficients of the inverted harmonic oscillator, we take advantage of the concept of analog experiments 14,26,27 and utilize the analogy between the propagation of quantum-mechanical waves and surface gravity water waves 20,21,28 .In the co-moving frame with group velocity c g , the water waves are governed by a wave equation that corresponds to the Schrödinger equation, where time and space are interchanged.This analogy enables us to transfer the wave properties of the quantum-mechanical problem to purely classical waves 29 .Moreover, in contrast to a quantum system, we can measure simultaneously their amplitude and phase 20 .For a more detailed discussion of the quantum analogy of surface gravity water waves, we refer to ref. 21 and the Methods section.
Figure 2 shows a schematic illustration of our experimental setup.A wave maker at one end of the water tank, which is 5 m long, 0.4 m wide, and 0.19 m deep, creates the initial wave packet that propagates along the water tank.Four capacitance type of wave gauges record the surface elevation, which corresponds to the complex envelope A jAj expðiφ A Þ. Hence, the imaginary part of A follows 29 from the Hilbert transform 30 of the real-valued elevation, see Methods section for more details.
In our experiments, we modulate the envelope of the surface gravity water wave to achieve the desired wave packet by the real function hðt; x ¼ 0Þ a 0 jAðtÞj cosðΩ 0 t þ φ A Þ, where a 0 is the amplitude of the carrier wave with frequency Ω 0 .In the Methods section, we present a comprehensive discussion on the generation of the truncated state A = A(t).Within this context, we have already made the necessary interchange between time t and position x, in accordance with the requirements of the water wave equation.
We measured the time-dependent elevations of the wave at different positions in the tank and stored them in a computer.The data recorded at the end of the canal was sent to the wave maker for a new excitation based on the previous measurement 31 .This process was repeated 4 times, enabling to achieve an effective propagation distance of 11 m in a 5 m long water tank.
In Fig. 3a, b, we compare and contrast the numerical simulation and the experimental results for the envelope of this free propagation.We prepare a surface gravity water wave corresponding to the wave packet, Eq. ( 1), with ε = 0.25 at x = 0.The free propagation of this wave packet leads to focusing at x = 10.4 m (red spot), indicating the amplitude singularity.From the raw data of the wave packets exemplified by (i) and (ii) and measured at the two positions indicated in (d) by arrows, we reconstruct the amplitude (c) and the phase (d) of the wave packet along the horizontal axis t = x/c g .Here, dots represent experimentally obtained values and the blue lines are numerical solutions of the truncated initial envelope with the same temporal truncation length as used in the experiments.The red curves represent analytical expressions outlined in the Methods section.
These results bring out most clearly the square-root singularity in the amplitude as well as the onset of logarithmic phase singularity.They not only manifest themselves along the spatial axis but also in the orthogonal direction, as shown in Fig. 3e, f, where the black curves represent the experimental values of amplitude and phase while the blue and red curves correspond to the simulation and analytical expressions, respectively.However, we note that in a truncated system, the expected growth does not exhibit full divergence.Rather, the amplitude peak tends to exhibit a smoother or flattened profile.Moreover, the phase singularity in such a system is influenced by truncation.In the context of an infinite wave packet, one would anticipate the logarithmic phase singularity to manifest as a clear divergence.However, the inherent finiteness of a truncated wave packet restricts this manifestation.Further details and methodologies associated with this observation can be found in the Methods section.
Our system also allows us to deduce the dependence of the reflection and transmission coefficients 32 on the energy reminiscent of the Fermi-Dirac distribution 25 F displayed in Fig. 4. For this purpose, we prepare the complex conjugate of the wave function defined by Eq. ( 1), and consider the propagation along the canal up to the position x = 10.4 m where the singularity emerges, and measure the transverse distribution of the amplitude envelope.For a fixed energy ε, we observe the amplitudes exemplified in Fig. 4a, c

Conclusions
In our article we have employed the similarity of the wave equations for quantum and surface gravity water waves to bring out the essence of Hawking radiation.Indeed, in our experiment, we observe the fingerprint of the logarithmic phase singularity emerging analogously at the event horizon of a black hole.For this purpose, we take advantage of the phase space horizons of the inverted harmonic oscillator.Here the logarithmic phase singularity manifests itself in the corresponding energy eigenfunctions and appears in phase space variables, that are not easily accessible.We emphasize that the observation of this effect does not require the presence of Illustration of a laboratory water wave tank measuring 5 m in length, 0.4 m in width, and 0.19 m in depth.The tank is constructed with a transparent glass sidewall and base, enclosed in an aluminum frame, and mounted on eight shock-absorbing legs for stability.A computer-controlled wave maker at one end generates water waves, while a wave energy absorbing beach at the opposite end minimizes residual reflections.Wave dynamics within the tank can be observed from all angles due to the tank's transparency.The instantaneous elevation of the water surface is monitored by four-wave gauges mounted on a movable bar, facilitating precise wave field analysis.
a potential barrier, but instead a free propagation of appropriately tailored initial wave packets.Although we are only able to prepare truncated wave packets our experimental results allow a clear identification of amplitude and phase singularities as appearing in an ideal system.
We point out that a square-root amplitude singularity is absent in the mode functions of the light field around the black hole.For this reason, we have observed a Fermi-Dirac rather than a Bose-Einstein distribution governing the transmission and reflection coefficients of the inverted harmonic oscillator.
These insights open pathways to another branch of experiments in black hole analogs that will propel our understanding of Hawking radiation.At the same time, many open questions emerge, and it suffices to mention only three: What would be an analog experiment to simulate the particle creation 3 in the ergosphere of a rotating black hole?Is there a deeper Bright and dark colors represent large and small values of the elevation, as indicated by the color bars on the right side.Our analysis is based on the raw data exemplified in (i) and (ii), yielding the amplitude (c) and the phase (d) of the water wave along the horizontal axis at t = x/c g .The experimental data (white dots), together with simulations (blue) and analytical predictions (red) given by Eq. ( 2), show both an amplitude and a logarithmic phase singularity at x = 10.4 m.In e, f, we display the transverse amplitude and phase at x = 10.4 m and compare the observation (black), simulation (blue), and analytical predictions (red).The horizon at t = x/c g separates a vanishing amplitude (t < x/c g ) from non-vanishing contributions (t > x/c g ), which also display an amplitude and a logarithmic phase singularity.
connection between the singularities and spin?Could black hole physics open a window toward the spin-statistics theorem 33,34 ?

Methods
In this section we briefly summarize elements of the quantum description 19 of the inverted harmonic oscillator crucial for our article, and review key properties of surface gravity water waves 20 .In addition, we provide an extensive overview of the experimental techniques which were used to perform the experiments.We start by re-deriving the amplitude and the logarithmic phase singularities in the energy wave functions of the inverted harmonic oscillator expressed in rotated phase space variables.Then, we recall the corresponding Wigner functions as well as the reflection and transmission coefficients.Moreover, we demonstrate that the singularities manifest themselves in space and time during the free propagation of these wave functions.We also address the limitations observed in truncated Weber wave packets.
We then turn to the discussion of the experimental details.In particular, we introduce the Schrödinger-like wave equation for surface gravity water waves and outline the conditions under which it is valid.Subsequently, we delve into the methods of extracting phase and amplitude data using the Hilbert transform, and present the primary features of our experimental setup, such as the water tank, the computer-controlled wave maker, and the use of capacitance-type wave gauges.The section concludes with a summary of the experimental protocol to generate and observe truncated Weber wave packets and their properties.

Energy wave functions
The quantum inverted harmonic oscillator of mass m and steepness ω is described by the Hamiltonian where the position operator x and the momentum operator p satisfy the familiar commutation relation ½x; p ¼ i_.
Next, we introduce the dimensionless operators and The operators ξ and η are intimately related to the familiar annihilation and creation operators â and ây of the standard harmonic oscillator by the transformation ω → iω and an overall phase factor.Hence, ξ and η are Hermitian operators whereas â and ây are non-Hermitian ones.
In terms of ξ and η, the Hamiltonian given by Eq. ( 4) takes the symmetric form Subsequently, the familiar eigenvalue equation determines the energy eigenstates |εi corresponding to the dimensionless energy eigenvalue ε.
Analogously, the ξ-representation Ψ ε (ξ) ≡ 〈ξ|ε〉 of an energy eigenstate |εi with regard to an eigenstate |ξi of the operator ξ given by Eq. ( 5) is Amplitude[mm] determined by the differential equation which yields the set of orthogonal solutions Expressed in the rotated phase space coordinates η and ξ, the horizons inherent in the stationary states of the inverted harmonic oscillator stand out most clearly.In fact, both sets of linearly independent solutions Φ ± ε and Ψ ± ε in ηand ξ-representation, respectively, have a similar functional dependence, differing only by the sign in front of ε.These energy eigenfunctions display a square-root amplitude singularity which is a consequence of the canonical commutation relation, and a logarithmic phase singularity at the horizons in phase space at η = 0 and ξ = 0, respectively, expressed by the Heaviside step function Θ.
We note that in ξ-representation each energy eigenstate |εi can be expressed as a superposition of the eigenfunctions Ψ þ ε and Ψ À ε .Analogously, in η-representation |εi can be decomposed into a superposition of the eigenfunctions Φ þ ε and Φ À ε .Moreover, we emphasize that the horizon and the logarithmic phase singularity in the energy eigenfunctions only become apparent in the coordinates ξ and η.
In fact, in position representation, the energy eigenfunctions and corresponding to Eqs. ( 10) and ( 12), respectively, do not display a logarithmic phase singularity nor a horizon.Instead, these functions are governed by a parabolic cylinder function 23 D ν (z) and a normalization factor Here, we have introduced the scaling factor Finally, we point out that the horizons in the inverted harmonic oscillator become evident in a phase space representation.Indeed, in complete analogy to the energy eigenfunctions Ψ ± ε and Φ ± ε , the Wigner functions display a horizon in phase space.Consequently, the Wigner function W þ Φ;ε corresponding to the wave function Φ þ ε , Eq. ( 10), vanishes on the half plane η < 0, while it depends only on the product ξη in the other half plane η > 0. Analogously, the Wigner function W þ Ψ;ε corresponding to the wave function Ψ þ ε , Eq. ( 12), vanishes on the half plane ξ < 0, while it depends only on the product ξη in the other halfplane ξ > 0.
Moreover, by making use of the c-number relations for ξ and η corresponding to Eqs. ( 5) and ( 6), we are able to obtain the Wigner functions and in terms of the position x and the momentum p.

Free propagation
Next, we consider the free propagation of an arbitrary energy eigenfunction of the inverted harmonic oscillator, expressed as a superposition of the orthogonal energy eigenfunctions ϕ ± ε , Eq. ( 13), with complex coefficients c ± that satisfy the normalization condition |c + | 2 + |c − | 2 = 1.Consequently, the time-evolved wave function φ = φ(x, t) is obtained as solution of the Schrödinger equation with the initial condition φ(x, 0) ≡ φ ε (x).
In the following we demonstrate that the function φ(x, t) reveals an amplitude as well as a phase singularity at the times t ± c ± 1=ω.Indeed, the solution of Eq. (23) given by the Fresnel transform takes the explicit form φðx; tÞ ¼ , being valid for any energy eigenfunction φ ε .For t > t þ c we instead obtain the expression φðx; tÞ ¼ which holds again for any energy eigenfunction φ ε .We point out that Eq. ( 25) is reminiscent of the one obtained 24 for the free time evolution of a symmetric energy eigenfunction of the standard harmonic oscillator.
In addition, at the center x = 0, the wave packet φð0; tÞ ¼ displays both an amplitude and a phase singularity at time t ¼ t ± c .So far, we have only studied the behavior of the wave function φ(x, t) for , that is, before and after the emergence of the singularity.In order to examine the wave function in the limit t !t ± c , we now analyze different options for the initial energy eigenfunction φ ε at t = 0. Indeed, we demonstrate that (i) for φ ε ¼ ϕ þ ε , Eq. ( 13), a horizon and the logarithmic phase singularity is revealed at 14), we are able to extract the transmission and reflection coefficients of the inverted harmonic oscillator at time t ¼ t þ c , which are reminiscent of the Fermi-Dirac distribution.

The horizon and the logarithmic phase singularity
In order to transfer the horizon and the logarithmic phase singularity of the inverted harmonic oscillator hidden in the variables η and ξ to the position coordinate, we consider the free time evolution of a particular energy eigenstate |εi, whose η-representation Φ þ ε ðηÞ is given by Eq. (10).For this purpose, we prepare at t = 0 the energy eigenfunction defined by Eq. ( 13), where the normalization constant N þ ðεÞ is given by Eq. ( 15).In the limit t !t þ c , we obtain according to Eqs. ( 25) and ( 26) the probability amplitude and consequently, the probability density Indeed, at the time t ¼ t þ c we are able to observe the logarithmic phase singularity as well as the horizon of the energy wave function Φ þ ε , Eq. ( 10), in the position variable x.Moreover, apart from a phase quadratic in x, the overall form of the propagated wave function given by Eq. ( 29) is reminiscent of the initial wave function Φ þ ε in the η-coordinate given by Eq. ( 10).We emphasize that a similar treatment for Ψ ± ε reveals a logarithmic phase singularity and a horizon in the position coordinate at t ¼ t À c , which is obtained by a free time evolution backward in time.

Influence of truncation parameter
We emphasize that in an idealized scenario, the infinite extension of the energy eigenfunctions of the inverted harmonic oscillator employed as an initial wave packet leads to singularities in both amplitude and phase during free propagation.However, in a realistic experimental setup, an infinite wave packet is unattainable and thus must be truncated.For the sake of simplicity, we employ for the truncation a rectangular window of width 2γ leading us to the initial wave function prepared at time t = 0.This truncation has consequences for the singularity appearing at time c in the free propagation ϕ þ ε ðx; t; γÞ of the initial wave function ϕ þ ε ðx; γÞ as shown in Fig. 5.The truncation of the wave packet imposes finite boundaries, which inherently limits the growth of the wave amplitude.In a theoretical infinite wave packet, the amplitude could grow without bound, leading to a singularity.When truncated, this growth is stymied, and instead of a true divergence, a smoothing or flattening of the amplitude peak may be observed.In order to bring this out most clearly, we show again in Fig. 6a the influence of the truncation parameter on the amplitude at t ¼ t þ c for the chosen parameters.
A similar behavior is observed for the phase singularity as indicated by Fig. 6b.In an infinite wave packet, the logarithmic phase singularity would manifest itself as a divergence, but the finite nature of a truncated wave packet constrains this behavior.The result is a tempered, less pronounced phase shift, which still follows a logarithmic trend but without exhibiting a clear singularity.

Fermi-Dirac distribution
Next, we show that the free time evolution not only reveals the inherent horizon and logarithmic phase singularity of the inverted harmonic oscillator, but also the characteristic transmission and reflection coefficients of the parabolic barrier.For this purpose, we make use of the energy wave functions Ψ ± ε , Eq. ( 12), which display a logarithmic phase singularity and a horizon in the ξ-representation.
According to Eq. ( 14), we thus prepare at time t = 0 the initial wave function where the normalization factor N À ðεÞ is given by Eq. ( 15).At time t ¼ t þ c , the time-evolved wave function is governed by the position-dependent functions with energy-dependent coefficients Consequently, the function ψ þ ε ðx; t þ c Þ consists of two distinct contributions in the domains x < 0 (−) and x > 0 (+), which are separated by a singularity at x = 0.
In terms of the probability density the respective regions in position are associated with either the transmission coefficient or the reflection coefficient of the parabolic scattering potential.Consequently, by propagating the energy eigenfunction ψ þ ε for different values of ε during the time t ¼ t þ c , we uncover the characteristic transmission and reflection coefficients 32 T and R of the inverted harmonic oscillator.
; γÞ of the propagated wave packet with truncation parameter γ as a function of the position x.For γ = ∞ (red), we obtain the ideal case displaying both singularities.For the truncation parameters γ = 5 (blue) and γ = 10 (black) the singularities are blurred and the respective values become finite.
At this point, it is worthwhile mentioning that T(ε) and R(ε) resemble a particular quantum statistics, namely the Fermi-Dirac 19 distribution.In fact, it is the presence of the amplitude singularity in addition to the logarithmic phase singularity in the energy eigenfunction of the inverted harmonic oscillator that results in the positive sign in the denominator of T and R. In contrast, disregarding the amplitude singularity and solely retaining the logarithmic phase singularity leads to a negative sign in the denominator, being reminiscent of the Bose-Einstein statistics.

Surface gravity water waves
In a frame moving with the group velocity c g , the evolution of the slowly varying complex-valued envelope A ≡ A(τ, ζ) of a surface gravity water wave follows from the wave equation reminiscent 29,35,36 of the Schrödinger equation ( 23) of a free particle, where A takes over the role of the wave function φ and the substitutions t → ζ, x → τ, ℏ → 1, m → 1/2, and i → − i have been applied in regard to the quantummechanical case.
The scaled dimensionless variables ζ and τ are related to the propagation coordinate x and the time t for the surface gravity water waves by ζ s 2 0 k 0 x and τ s 0 Ω 0 ðx=c g À tÞ.The carrier wave number k 0 and the angular carrier frequency Ω 0 satisfy the deep-water dispersion relation Ω 2 0 ¼ k 0 g with g being the gravitational acceleration, and define the group velocity c g ≡ Ω 0 /2k 0 .The parameter s 0 ≡ k 0 a 0 characterizing the wave steepness is assumed to be small (s 0 ≪ 1) in the linear regime.

Phase and amplitude from Hilbert transform
The time-dependent elevation hðx; tÞ ¼ a 0 Aðx; tÞ cos½Ω 0 t À k 0 x þ φ A ðx; tÞ a 0 u x ðtÞ ð40Þ of the water surface at any point x in the tank involves the envelope function, characterized by the amplitude A(x, t) and the phase φ A (x, t), which both vary slowly with respect to the carrier wave of frequency Ω 0 .The Hilbert transform 30 of the real-valued function u x (t) reduces in our case to v x ðtÞ ≈ Aðx; tÞ Using the Euler formula, we can define the complex function: with the total phase φ A tot ðx; tÞ Ω 0 t À k 0 x þ φ A ðx; tÞ and thereby obtain the amplitude and phase of the envelope function directly from the measured function u x (t) and its Hilbert transform v x (t).
To numerically compute the Hilbert transform from the measured surface elevation data, we use the Matlab toolbox function 'hilbert'.

Experimental facility
The experimental facility consists of three essential ingredients (i) a wave tank, (ii) a wave maker, and (iii) wave gauges.In this section, we summarize important details.
Surface gravity water wave tank.Our experiments were performed in a 5 m long, 0.4 m wide, and 0.19 m deep laboratory wave tank, as illustrated in Fig. 2. The wave tank is encompassed by an aluminum extrusion frame and supported by eight shock-absorbing legs.The sidewalls and base of the wave tank are made of transparent glass to permit flow visualization of the wave field and the observation of the waves from all angles.At each end of the test section, openings in the tank floor permit tank filling and discharging via particle filter.Prior to each experiment, the water surface was cleaned to remove any surface film that could influence the results.The use of transparent glass allowed the observation of the water's surface from above, and capacitance-type sensors were inserted into the test section at any distance from the inlet.The measuring equipment, power supplies, and sensors are supported by the instruments carriage, which is constructed of aluminum extrusions and affixed on a rail along the test section.The position of the carriage along the test section at the intended fetch and the input at the wave maker are the controlled experimental parameters.
In our experiments, we have used the carrier frequency Ω 0 = 15 rad/s and the initial amplitude a 0 = 2.0 mm.Moreover, k 0 satisfies the deep-water condition 21  Computer-controlled wave maker.A mechanical wedge-type wave maker (Linmot T01-72/420-1ph) is used to generate surface gravity water waves.It is composed of a wedge-shaped plate, a motor or actuator, a frame to hold the wedge, and a water tank.The wedge plate is mounted on the frame, and the frame is positioned above the water tank.The motor or actuator is connected to the wedge plate, which is positioned at the surface of the water and moves back and forth in a reciprocating motion.As the wedge plate moves, it creates a disturbance in the water, which generates waves that propagate outwards from the wedge.The size and frequency of the waves can be adjusted by changing the amplitude and frequency of the wedge's motion.
The initial envelope of the surface water gravity wave is determined by the quantum-mechanical energy eigenfunctions of the inverted harmonic oscillator.Therefore, the wave packet generated by the wave maker reads with the carrier frequency Ω 0 = 15 rad/s and the initial amplitude a 0 = 2 mm.Here A = A(t) and φ A = φ A (t) represent the amplitude and phase of the eigenfunctions obtained from Eqs. ( 13) and (14).For this purpose we make use of the substitutions as defined previously in order to translate quantummechanical variables to the formalism of surface gravity water waves.Then, the time-dependent elevation of the wave was measured using wave gauges at different positions in the tank and stored in a computer.This data was numerically truncated at the location of the next temporal slit, and was then sent to the wave maker for a new excitation based on the previous measurement.In this way, we were able to cascade several slits in the time domain and observe the effect of diffractive guiding.
Capacitance-type wave gauges.We utilize a capacitive wave gauge transducer consisting of a 0.3 mm thin tantalum wire coated with a uniform thin layer of tantalum oxide.Before the measurement of the surface wave height using wave gauges, the calibration of the wave sensor is performed in three steps: (i) we first set the vertical position of the wire wave gauge in such a way that the mean water level is approximately in the middle of the wire.(ii) We perform automatic calibration of the wave gauge using a lab-view-made routine by submerging the wave sensor at different depths and recording the mean voltage for 5 s for each depth.(iii) Finally, a polynomial fit is applied to the recorded data in order to find out the dependence between measured height H[mm] and the gauge voltage and verify visually the fitted calibration polynomial as shown in Fig. 7.We note that the analysis of the relationship between height and voltage in capacitor type wave gauges, the precision of measurements was evaluated by employing statistical error calculations based on five separate measurements for each data point.This approach ensured that the estimated errors reflect the variability inherent in the experimental setup and measurement process.Notably, the resulting error bars were found to be smaller than the symbols representing the data points on the plotted graph.

Experimental protocol
In this section, we outline the experimental protocol to generate the truncated Weber wave packets and observe the singularities appearing in Eq. (10).In particular, we discuss the propagation along the wave tank, and its recording and transformation into a full complex wave function.The iterative procedure is repeated until the observation of a singularity in both amplitude and phase.
We proceed in the following steps: (i) Generation of truncated Weber wave packet.The Weber wave packet defined by Eq. ( 1) is not normalizable.For the experiment, we had to truncate them to a finite duration in time.The temporal window of truncation was carefully chosen based on experimental constraints and theoretical predictions.In Fig. 8, we show the full wave packet (black curve) and the truncated Weber wave packet (red).The truncated Weber wave packet was generated by the computer-controlled wave maker.(ii) Propagation and recording.The so-generated truncated Weber wave packet was then allowed to propagate along the carefully controlled water tank.The height of the propagating wave packet, that is the real part, was recorded using the high-resolution wave gauges positioned at specific locations along the tank discussed before.(iii) Translation into complex wave function.Each recorded real part of the wave packet envelope was then translated into a full complex wave function using the Hilbert transform discussed above.(iv) Reconstructing and sending back to the wave maker.The wave packet was then propagated to the final point, x f , recorded, and recreated again at the wave maker.(v) Iterative procedure and singularity observation.This entire process was repeated iteratively-each time observing the properties of the propagated, transformed, and re-sent wave packet.The iteration was continued until we observed a singularity in both the amplitude and phase of the wave packet.The criteria for the observation of a singularity were based on a combination of theoretical predictions and experimental viability.Each data point is derived from an average of five separate measurements at each height to ensure accuracy and reliability.This experimental data undergoes an analysis through linear regression methods to establish the relationship between height and voltage accurately.To evaluate the precision of these measurements, statistical error calculations were employed, which factored in the variability inherent in the experimental setup and the measurement process.Notably, the calculated error bars, indicative of measurement precision, were found to be significantly smaller than the symbols used to represent the data points on the graph, underscoring the high degree of accuracy in the experimental findings.

Fig. 1 |
Fig. 1 | Emergence of a horizon in phase space and an amplitude singularity in the free time evolution of an energy eigenstate of the inverted harmonic oscillator.a Phase space trajectories of classical particles (top) approaching a parabolic barrier V(x) (bottom) from the left (right) are represented by the blue (green) and red (orange) hyperbolas.Depending on their energy E the particles are transmitted (E > 0) or reflected (E < 0).For E = 0, they define a horizon in phase space (black), which separates the different energy domains.b Phase space representation of a quantum particle (top) with energy ε corresponding to the energy wave function, Eq. (1), as a function of position x and momentum p.In the neighborhood of the classical trajectory corresponding to the energy E = ℏωε, the corresponding Wigner function W þ ϕ;ε ðx; pÞ has a dominant maximum (red).In contrast to a classical particle, the Wigner function covers two domains in phase space due to tunneling.The horizons confine the Wigner function to one half-plane in phase space.We show the corresponding position distribution jϕ þ ε ðxÞj 2 (bottom) resulting from integration over the momentum variable.The time evolution of such a Wigner function in the absence of any potential is governed by a sheering in phase space.c At a particular time t ¼ t þ c an amplitude singularity emerges in the position distribution while the horizon at x = 0 separates vanishing from non-vanishing parts.d At later times t > t þ for the energies ε = −0.4 and ε = 0.4.This procedure applied to different energies yields the experimental values (open circles) in Fig. 4b which agree well with the theoretical prediction (solid lines) of the Fermi-Dirac distribution.

Fig. 2 |
Fig.2| Experimental setup for generating wave packets of surface gravity water waves and measuring their propagation dynamics in both amplitude and phase.Illustration of a laboratory water wave tank measuring 5 m in length, 0.4 m in width, and 0.19 m in depth.The tank is constructed with a transparent glass sidewall and base, enclosed in an aluminum frame, and mounted on eight shock-absorbing legs for stability.A computer-controlled wave maker at one end generates water waves, while a wave energy absorbing beach at the opposite end minimizes residual reflections.Wave dynamics within the tank can be observed from all angles due to the tank's transparency.The instantaneous elevation of the water surface is monitored by four-wave gauges mounted on a movable bar, facilitating precise wave field analysis.

Fig. 3 |
Fig. 3 | Experimental observation of a horizon in phase space giving rise to amplitude and logarithmic phase singularities in surface gravity water waves.The free propagation of the energy eigenfunction of the inverted harmonic oscillator determined by Eq. (1) with ε = 0.25 displays in the simulation (a) as well as in the experiment (b) a singularity at t = x/c g and x = 10.4 m, indicated by the orange arrows.Bright and dark colors represent large and small values of the elevation, as indicated by the color bars on the right side.Our analysis is based on the raw data exemplified in (i) and (ii), yielding the amplitude (c) and the phase (d) of the water wave along the

Fig. 4 |
Fig.4| Experimental observation of the Fermi-Dirac distribution with surface gravity water waves.The free propagation of the energy eigenfunction of the inverted harmonic oscillator, determined by the complex conjugate of Eq. (1), results in an amplitude singularity at x = 10.4 m and t = x/c g .In a we display the amplitude singularity for an initial wave packet below the parabolic barrier (ε = −0.4).We present the observation (black) together with the simulation (blue) and the analytical prediction (red).The transmission and reflection coefficients T(ε) and R(ε) are extracted following Eq.(3) from the measured amplitude at x = 10.4 m by fits along the transverse coordinate for t < x/c g and t > x/c g , respectively.b We recover the dependence of T (green) and R (orange) on ε by varying the energy of the initial wave packet.Here the open circles represent the experimentally obtained values in comparison to the analytical expressions (solid lines), given by R(ε) = F(ε) and T(ε) = F(−ε) where FðεÞ ½1 þ expð2πεÞ À1 represents the Fermi-Dirac distribution.c For an initial wave packet above the parabolic barrier (ε = 0.4), we observe a mirror image of a at t = x/c g for the analytical prediction (red) and numerical simulation (blue) leading to the symmetry of the transmission T and reflection coefficient R with respect to the energy eigenvalue ε = 0.

Fig. 5 |
Fig. 5 | Influence of the truncation of the initial wave packet ϕ þ ε ðx; γÞ on the amplitude for three values of the truncation parameter γ.For finite values of γ, that is, a truncated initial Weber wave packet, illustrated here for a γ = 5 and b γ = 10, the singularity at x = 0 is smoothed out.Only for γ = ∞ in c we obtain a divergence at t ¼ t þ c and x = 0, indicating the amplitude singularity.

Fig. 6 |
Fig. 6 | Influence of the truncation of the initial wave function on the amplitude and phase singularity.At time t¼ t þ c , we display the a amplitude jϕ þ ε ðx; t ¼ t þ c ; γÞj 2 and b phase ϑ þ ϵ ðx; t ¼ t þ c; γÞ of the propagated wave packet with truncation parameter γ as a function of the position x.For γ = ∞ (red), we obtain the ideal case displaying both singularities.For the truncation parameters γ = 5 (blue) and γ = 10 (black) the singularities are blurred and the respective values become finite.
k 0 d > π, where d = 0.19 m denotes the depth of the tank and the corresponding steepness is s 0 < 0.04 guaranteeing the validity of the linear Schrödinger equation.More details on the experimental setup can be found in refs.20,36.

Fig. 7 |
Fig.7| Set of calibration curves for the four capacitance wave gauges.Comprehensive calibration curves for four distinct capacitance wave gauges, detailing measurements across six different heights (H) in millimeters (mm) and their respective recorded voltages (V) in volts.Each data point is derived from an average of five separate measurements at each height to ensure accuracy and reliability.This experimental data undergoes an analysis through linear regression methods to establish the relationship between height and voltage accurately.To evaluate the

Fig. 8 |
Fig.8| Illustration of the truncated initial wave packet for surface gravity water waves.Full initial envelope waveform based on Eq. (1) for ε = 0.25 (black solid lines) and the truncated one for a truncation window of 50 s (red solid lines).
Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material.If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.