1. Introduction
The question of whether the Hawking evaporation [
1] violates unitarity and, therefore, results in the loss of information [
2], has remained unresolved since Hawking’s seminal discovery. The proposed solutions include black hole complementarity [
3], firewalls [
4,
5] (see, for example, [
6,
7], for a recent review and [
7,
8,
9] for a counterargument), soft hairs [
10], black hole remnants [
11], islands [
12,
13], replica wormholes [
14,
15], and instanton tunneling between multiple histories of Euclidean path integrals [
16]. So far, the investigations remain mostly theoretical since it is almost impossible to settle this paradox through direct astrophysical observations, as typical stellar-size black holes are cold and young; however, the solution to the paradox depends crucially on the end-stage of the black hole evaporation.
There have been proposals for laboratory investigations of the Hawking effect, including sound waves in moving fluids [
17,
18], electromagnetic waveguides [
19], traveling index of refraction in media [
20], ultra-short laser pulse filament [
21], Bose–Einstein condensates [
22], and electrons accelerated by intense lasers [
23]. It should be emphasized that the Chen–Tajima proposal [
23] differs from other concepts mentioned above in that it is based on the equivalence principle, which mimics the Hawking radiation of an
eternal (non-dynamical) black hole. Experimentally, Reference [
22] reported on the observation of a thermal spectrum of the Hawking radiation in the analog system and its entanglement. However, most of these are limited to verifying the thermal nature of the Hawking radiation.
It has long been recognized that accelerating mirrors can mimic black holes and emit Hawking-like thermal radiation [
24]. In 2017, Chen and Mourou proposed a scheme to physically realize a relativistic mirror by using a state-of-the-art high-intensity laser to impinge a plasma target with a decreasing density [
25,
26]. The proposal follows the same philosophy as Reference [
23], but differs in that it mimics the Hawking radiation from
gravitational collapse and, therefore, from a dynamical black hole, which is a more direct analogy to the original Hawking evaporation. It is also unique in that it does not rely on a certain fluid to mimic the curved spacetime around a black hole, but rather a more direct quantum field theoretical analogy between the spacetime geometry defined by a black hole and a flying mirror.
Based on this concept, an international AnaBHEL collaboration has been formed to carry out the Chen–Mourou scheme, which is the only experimental proposal of its kind in the world. Our ultimate scientific objectives are to detect analog Hawking radiation for the first time in history and through the measurement of the quantum entanglement between the Hawking particles and their vacuum fluctuating pair partner particles, to shed some light on the unresolved information loss paradox. From this perspective, the AnaBHEL experiment may be regarded as a
flying EPR (Einstein–Podolsky–Rosen) experiment [
27].
The concept of a flying plasma mirror was proposed by Bulanov et al. [
28,
29,
30,
31]. It provides an alternative approach to the free electron laser (FEL) in generating high-frequency coherent radiation. The flying plasma mirror approach provides a great prospect for future applications. A series of proof-of-principle experiments led by Kando at KPSI in Japan [
32,
33,
34,
35] has validated the concept. However, the mirror reflectivity (as a function of frequency) as well as other physical properties, such as the reflection angular distribution, etc., have not been characterized in those two experiments.
In this paper, we first review the physics of a flying mirror as an analog black hole. We then reveal the concept of accelerating relativistic plasma mirrors as analog black holes, with the attention paid to the aspects pertinent to the investigation of the Hawking radiation and the information loss paradox, including the laser–plasma dynamics that give rise to the acceleration of the plasma mirror, the reflectivity, the frequency shift of the reflected spectrum, and corrections due to the finite-size and semi-transparency effects of a realistic plasma mirror to the blackbody spectrum of the analog Hawking radiation based on an idealized, perfectly reflecting, point mirror. We then report on the progress of our R&D, i.e., of the key components in the AnaBHEL experiment, including those of the supersonic gas jet and the superconducting nanowire single-photon Hawking detector. We conclude by projecting our experimental outlook.
2. Flying Mirror as Analog Black Hole
Figure 1 depicts the analogy between the Hawking radiation of a real BH (left) and that of an accelerating mirror (right). The fact that accelerating mirrors can also address the information loss paradox was first suggested by Wilczek [
36]. As is well-known, the notion of black hole information loss is closely associated with quantum entanglement. In order to preserve the “black hole unitarity”, Wilczek argued that, based on the moving mirror model, in vacuum fluctuations, the
partner modes of the Hawking particles would be trapped by the horizon until the end of the evaporation, where they would be released and the initial pure state of a black hole would be recovered with essentially zero cost of energy. More recently, Hotta et al. [
37] argued that the released partner modes are simply indistinguishable from the zero-point vacuum fluctuations. On the other hand, there is also the notion that these partner modes would be released in a burst of energy, for example, in the Bardeen model [
38] (See
Figure 2).
One common
drawback in all analog black hole concepts involves setting up in a laboratory with flat spacetime (therefore, the standard quantum field theory is known to be valid); it is inevitable that any physical process, including analog black hole systems, must preserve the unitarity. Therefore, none of the proposed analog black holes can in principle
prove the loss of information even if that is indeed so. The real issue is, therefore, not so much about whether the unitarity is preserved, but more about
how it is preserved. That is, it is even more important to determine how the black hole information is retrieved. Does it follow the Page curve [
39], a modified Page curve where the Page time is significantly shifted towards the late time [
16], or alternative scenarios [
40]? The measurement of the entanglement between the Hawking particles and the partner particles as well as the evolution of the entanglement entropy [
41], should help to shed much light on the black hole information loss paradox. As pointed out by Chen and Yeom [
41], different scenarios of black hole evolution can be tested by different mirror trajectories [
42,
43].
4. Analog Hawking Temperature
There exists a wealth of literature on the vacuum fluctuating modes of quantum fields, their reflections from a flying mirror, and the analog “Hawking temperature” of such a flying mirror as an analog black hole [
49]. In general, such analog Hawking temperature depends on the actual mirror trajectory. According to Reference [
26],
where
is the speed of light in the plasma medium, which is position dependent. In our conception [
25], the plasma target thickness is supposed to be much larger than the characteristic scale of the density variation, i.e.,
. In this situation, it is safe to extend the integration to
(and
). Taking this approximation, we find
where
and
. This is identical to the Davies–Fulling trajectory, i.e., Equation (4.51) of Reference [
49],
where
are positive constants and
.
Transcribing the
coordinates to the
coordinates, where
and
, we see that only null rays with
can be reflected. All rays with
will pass undisturbed. The ray
, therefore, acts as an effective horizon [
49]. Following the standard recipe [
49], we obtain the Wightman function as
where
in the
limit. The constant factors in the argument of the log function in the above equation do not contribute to the nontrivial part of the physics. Note that in our notation,
t is the time when the ray hits the mirror. Let us denote the observation time and position by
T and
X. Then
. For large
t,
. This leads to
for a static mirror at
Integrating over
T and
, we then have, in the asymptotic limit of
,
This leads to the response function (of the particle detector) per unit of time with the form
where the analog Hawking temperature of the mirror measured by a stationary particle detector is
Here, is the Boltzmann constant. It is interesting to note that the analog Hawking temperature associated with our constant-plus-exponential-squared density profile depends strongly on the characteristic length D and only weakly on the plasma density (through ). This points to the possibility of employing gaseous instead of solid plasma targets, which would greatly simplify our proposed experiment.
5. Conceptual Design
The original experimental concept proposed by Chen and Mourou [
25] invoked a two-plasma-target approach, where the first plasma target converts an optical laser into an X-ray pulse through the flying plasma mirror mechanism. The converted X-ray pulse then impinges on a nano-thin-film that is fabricated with a graded density in different layers. This design has the advantage of having a solid-state density, providing a higher plasma frequency, which is proportional to the square root of the plasma density and, therefore, a higher density gradient for maximizing the Hawking temperature. On the other hand, the drawbacks of this concept are multiple. First, the typical conversion efficiency of flying plasma mirrors is ∼
, rendering it difficult for the converted X-ray pulse to remain in the nonlinear regime. Second, the solid plasma target would induce extra backgrounds, which are linearly proportional to the target density.
In 2020, Chen and Mourou proposed a second design concept [
26], where the conversion of optical laser to X-ray was no longer needed and, thus, the first plasma target was removed, and the nano-thin-film solid plasma target was replaced by a supersonic gas jet. This largely simplifies the design and the technical challenges.
Figure 5 shows a schematic conceptual design of the single-target, optical laser approach. The key components now reduce to a supersonic gas jet with a graded density profile and a superconducting nanowire single-photon Hawking detector, the R&D progress of which will be described in later sections.
In our design of the AnaBHEL experiment, we assume the driving laser has the frequency
and the wavelength
. For the plasma target, we set
in Equation (
2) so that
, and we assume
. The corresponding plasma frequency is
and the plasma wavelength
. Next, we design the plasma target density profile. Since our formula is not constrained by the adiabatic condition, we are allowed to choose a minute characteristic length
. Then we find
which corresponds to a characteristic Hawking radiation frequency
. Thus, the Hawking photons can propagate through and out of the plasma for detection.
6. Hawking Photon Yield
Among the proposed models, the physics of flying/moving mirrors is perhaps the one closest to that of real black holes, since in both cases the radiation originated from vacuum fluctuations. The essence of the Hawking radiation lies in the gravitational redshift of the field modes’ phases. Since the key is the phase shift, various analog models or experimental proposals attempt to generate the same phase shift as that of the Hawking radiation but now in flat spacetime, i.e., laboratory. Indeed, in the flying mirror model, the gravitational redshift is mimicked by the Doppler redshift.
Due to the spherically symmetric nature of typical black hole spacetimes, the spherical coordinate origin is effectively a perfectly reflecting point mirror and the corresponding Hawking radiation is expected to be emitted radially, hence the situation is effectively (1+1)-dimensional and, thus, most of the flying mirror literature only considers a real perfectly reflecting point mirror in (1+1)-dimensional flat spacetime. Nevertheless, in the laboratory, the spacetime is (1+3)-dimensional. In addition, our proposed relativistic flying mirror generated through laser–plasma interaction has a low reflectivity [
26] and a finite transverse/longitudinal size; therefore, it is necessary to take these practical effects into consideration to estimate the particle production yield.
The standard treatment in the flying mirror model [
24,
48,
49,
51] considers a real scalar field in (1+1)D flat spacetime subjected to a single, relativistic, time-dependent Dirichlet boundary condition in space to represent a relativistic perfectly reflecting point mirror. Since the boundary condition is externally provided, the breakdown of Poincaré invariance leads to the possibility of spontaneous particle creations following quantum field theory.
The generalization of this standard calculation to a flying plasma mirror with a finite reflectivity in
n-dimensional flat spacetime can be made by starting with the action functional [
52]:
where natural units are employed,
is the coupling constant with dimension of mass,
is the fine structure constant,
is the surface density of the electrons on the mirror, and
encodes the mirror’s trajectory
, longitudinal/transverse distribution
, and the Lorentz factor
.
Solving the equation of motion for
with the in-mode/out-mode boundary conditions in (1+1) dimensions, one finds (assuming the field to be in the in-vacuum state
with the mirror flying to the negative
x-direction) the created particles (due to the field mode reflected to the mirror’s right to have the frequency spectrum) [
53,
54,
55,
56]:
where
and
is the incident/emitted plane wave mode’s frequency,
is the mirror’s reflectivity,
, and
is the phase shift/ray-tracing function induced upon reflection off the receding mirror. From Equation (
18), one sees that for a given trajectory
, the spectrum would be different depending on the reflectivity.
A simple model that mimics the formation and evaporation of a Schwarzschild black hole is the collapse of a spherical null shell. In this scenario, the relevant ray-tracing function is
where
is the black hole’s surface gravity, and
is the past event horizon, which is conventionally set to zero. For field modes propagating in the vicinity of
(late time),
, and
(extreme gravitational/Doppler redshift), one obtains
for a perfectly reflecting point mirror, and
for a semi-transparent point mirror, where
is the analog Hawking temperature. In general, the accelerating mirror radiates along the entire worldline, but only those radiated in the late time are relevant to the analog Hawking radiation. In particular, the spectrum Equation (
20) for a perfectly reflecting point mirror is in exact accordance with the Hawking radiation emitted by a Schwarzschild black hole. Although a semi-transparent point mirror possesses a different spectrum due to the time-dependent and frequency-dependent reflectivity, it nevertheless has the same temperature as that of a perfectly reflecting point mirror.
As previously mentioned, practical considerations in the laboratory force us to work in (1+3)-dimensional spacetime and a mirror with some kind of longitudinal/transverse distribution. In the case of a semi-transparent mirror, it is possible to find the corresponding analytic spectrum through a perturbative approach. The result is
where [
56]
where
is the incident/emitted plane wave mode frequency, respectively, and
is the form factor due to the mirror’s longitudinal and transverse geometry, which is independent of the mirror’s motion and reflectivity.
According to particle-in-cell (PIC) simulations [
50], a mirror of square-root-Lorentzian density distribution and a finite transverse area can generate a good-quality mirror. Thus, we shall consider the case:
where
W is the half-width at half maximum of the square-root-Lorentzian distribution and
is the transverse area. In addition, according to the plasma density profile designed in Reference [
26], the mirror follows the trajectory:
where
are positive plasma mirror parameters and time
t is written as a function of the trajectory
. This trajectory is designed such that it approximates the black hole-relevant trajectory:
either (i) at the late-time
for any value of
b, or (ii) in a near-uniform plasma background
during the entire accelerating phase. In either case, the spectrum relevant for the analog Hawking radiation is
where
is the effective temperature,
, and
are complicated form factors due to the mirror’s transverse/longitudinal distributions given in Reference [
56]. Notice that the form factor
leads to diffraction, whereas
may enhance the production rate.
Using the PIC simulation parameter values:
eV,
eV (
),
eV,
eV
(1.5 nm),
eV
(50
), and
, the resulting analog Hawking temperature is
eV (369 K) in the far infrared regime and the number of produced analog Hawking particles per laser crossing is
where 0.27 and 0.02 correspond to the red and the blue areas in
Figure 6, respectively.
Assuming a petawatt-class laser, such as that in the Apollon Laser Facility in Saclay, France, which can provide 1 laser shot per minute and 8 h of operation time per day, a 20-day experiment with a 100% detector efficiency would give the total yield of events as
It should be reminded that this value is highly idealized. Fluctuations of the physical parameters, especially that of the characteristic length of the density gradient, D, which we have not yet measured, would impact the expected Hawking photon yield.
7. Supersonic Gas Jet
As estimated in Reference [
25], the gradient of the electron number density required for the experiment is ∼
, which is attainable with a supersonic gas jet. There are several methods proposed in the literature, such as a shock wave generated induced by a laser that propagates perpendicular to the gas jet [
57,
58], and a supersonic gas flow impinged by a thin blade [
59,
60]. The estimated gradients of the electron number densities reached by different groups in [
57,
59,
60] are summarized in
Table 1. It is clear that, in principle, both methods can provide gradients that satisfy our requirement. As our first attempt, we chose the latter method for its simplicity.
The supersonic gas jet can be realized by passing high-pressure gas through the de Laval nozzle, which is also known as the converging-diverging nozzle. The gas flow will reach sonic speed at the throat of the nozzle and then be accelerated in the diverging section to reach supersonic speed. Based on the design of the nozzle in [
61], we produce our own nozzle to generate supersonic gas flow.
Figure 7 shows the inner geometry and the image of the nozzle we built. The nozzle is connected to the tank of an air compressor that can provide air with pressure up to 8 atm. An electrically controlled valve is placed between the nozzle and the tank to control the flow.
There are several techniques to quantitatively characterize the density of a supersonic gas jet, including interferometry and shadowgraphy [
62,
63,
64,
65], tomography [
62,
66,
67], planar laser-induced fluorescence (PLIF) [
60,
68,
69], Schlieren optics [
70,
71] (more references can be found in [
71]). As the first step, we built a Schlieren imaging system in the lab for the jet characterization. Our Schlieren optics is equipped with a rainbow filter, which allows for the visualization of the gas jet as well as quantitative analysis of its refractive index.
Figure 8 demonstrates the schematic diagram of our system.
The principle behind the Schlieren optics is that the variation of the refractive index would diffract light. A rainbow filter that intercepts the diffracted light then provides information that would quantitatively determine the diffraction angle according to the color codes. The imaging system is calibrated with a plano-convex lens, whose refractive index is known. In this way, the map of the refractive index gradient, which is directly related to the gas density gradient, can be obtained.
Figure 9 shows the image using our Schlieren optics. The figure shows the supersonic jet produced by the nozzle. The so-called “shock diamonds” are clearly demonstrated, which is an indicator of the jet propagating with supersonic speed in the atmosphere.
The design of the nozzle is verified by comparing the shock diamond structure from the data with the computational fluid dynamic (CFD) simulation result. The 3D fluid simulation was performed with OpenFOAM code. In the simulation, a compressible Navier–Stokes flow solver, rhoCentralFoam [
72], is used to study the behavior of the supersonic jet.
With the conventional Abel inversion technique, the gradient of the refractive index was reconstructed and compared with the simulation result in
Figure 10. Line profiles at different horizontal positions,
y, relative to the axial center of the gas jet are shown in
Figure 11. We found the positions of several peaks in the data agree reasonably with simulation results. This implies the behavior of our self-made supersonic nozzle is as expected and our Schlieren optics can characterize the profile of the supersonic jet. Further improvement is ongoing to obtain results with higher accuracy.
8. Superconducting Nanowire Single-Photon Hawking Detector
Observing the Hawking photons is the main goal and one of the major challenges of the planned AnaBHEL experiment. There is probably no single technology that satisfies all requirements. The detector must be a single photon detector, with efficiency close to 100%. The desired Hawking photon sensitivity wavelength range should be from 10 m to 100 m. A second detector design is required for the forward-moving partner photon with sensitivity at the UV (1–100 nm). The low expected signal yield and the potentially large asynchronous thermal and plasma-induced backgrounds set stringent detector timing requirements (to picosecond level or better). Since within the data acquisition timing window accidental coincidences may still be present, single photon pair polarization measurement will be required in order to unambiguously tag the pair as Hawking and partner photons. In addition to the above requirements, the detector should have a very fast recovery to avoid photon pile-up, a very low dark current rate (DCR), and the ability to cover relatively large areas.
Superconducting nanowire single-photon detectors (SNSPDs) satisfy most of the above requirements [
73]. Thin superconducting films (∼10 nm) from materials such as NbN and WSi are sputtered on substrates. Subsequently, electron nanolithography is used to etch narrow wire structures (50–100 nm wide). The detector operates at a temperature below the Curie temperature
at an appropriate bias current that maximizes efficiency. Additional cavity structures are needed in order to bring the efficiency close to 100%.
The intrinsic time jitter of SNSPDs is ∼1 ps. Recently, time jitters using short straight nanowires and are found to be <3 ps for NbN [
74] and 4.8 ps for WSi wires [
75]. Thanks to their short reset time, these devices exhibit very high count rates at the level of hundreds of MHz. Although the expected Hawking photon yield is low, such a fast recovery detector reduces dramatically the probability of photon pileup (multiple counts in the same time window). The dark count rate (DCR) is extremely low at the level of one count for a period of hours, depending on the operating temperature and the bias current.
Typical SNSPD designs relevant to AnaBHEL are based on a superconducting nanowire patterned from a thin film of thickness between 5 and 10 nm. The most common nanowire design follows a meandering structure geometry. However, in our case, we need to consider specific structures that have sensitivity to polarization. SNSPDs are DC-biased with operation currents close to their critical currents so that efficiency is maximized. As discussed in [
73], the detection process is divided into the following steps: (I) Photon absorption; (II) Creation of quasiparticles and phonons combined with their diffusion; (III) Emergence of a non-superconducting nanowire segment; (IV) Redirection of the bias current in readout circuitry, leading to a voltage pulse; and (V) detector recovery.
During step (II), the impinging near-IR photon photo-excites an electron (the relaxation of which leads to the formation of a cloud of quasiparticles and phonons). An instability of the superconducting state emerges due to the quasiparticle cloud, which results in the reduction of the effective critical current density and a part of the nanowire experiences a transition to the non-superconducting state (III). The occurrence of a normal-conducting hot spot in the nanowire can lead to the detection of the photon event as the current flowing through the bias resistor (bias current) is re-directed. Due to internal Joule heating, the resistive domain of the nanowire keeps growing, which leads to increased resistance at the level of k. This significant non-zero resistivity causes the redirection of the bias current from the nanowire to the readout electronics (IV). Finally, the resistive domain is cooled down and the superconductivity is restored, bringing the nanowire back to its initial state (V).
Specific requirements of the AnaBHEL experiment photon sensors are summarized in
Table 2 (first row). Realistic operational parameters and performance for typical SNSPD materials are also presented.
In most applications, SNSPDs are coupled to fibers with a typical operation wavelength at the telecom window (1550 nm). AnaBHEL is an open-air experiment with a tight requirement of operation at mid to far infrared (
m) regime. As reported in [
76,
77], significant progress has been made for open-air longer wavelength operating SNSPDs. To achieve sensitivity for wavelengths longer than 10
m, materials of lower Curie temperatures must be used. WSi is an example of such material. However, further R&D on other materials is needed.
In addition to efficiency, successful detection of the Hawking and partner photons in AnaBHEL requires good detector acceptance in both the forward and backward parts of the experimental apparatus. A single-pixel SNSPD covers a very small active area of the order of
m
. To maximize photon acceptance, a
mm
pixel array would be preferred. This kilopixel array has already been produced [
78] and used in exoplanet transit spectroscopy in the mid-infrared regime.
Hawking Photon Sensor Fabrication and Characterization
In 2021, a R&D program was initiated in Taiwan to develop photon sensors for Hawking photon detection. Academia Sinica, NTU, and NCU groups are currently sharing equipment and laboratories for the fabrication and testing of prototype SNSPDs, the preferred technology for Hawking-photon sensors.
We have been producing NbN films of 10 nm thickness using the Academia Sinica magnetron sputtering machine shown in
Figure 12 (Kao Duen Technology, Model: KD-UHV, N-11L17). The films grown on two different substrates, MgO and distributes Bragg reflector (DBR), were used. The films were sputtered at UHV pressure of
Torr. Sample NbN films on a sample holder are shown in
Figure 13.
The superconducting transition properties of the NbN films have been determined using magnetic susceptibility measurements with a SQUID, as well as electric resistivity measurements. In the left side of
Figure 14, the MPMS3 SQUID magnetometer is used to measure the magnetic susceptibility of the NbN samples grown on MgO. On the right side of the same figure, the superconducting transition is shown as the material becomes diamagnetic. A
NbN sample was placed in the SQUID and its magnetic susceptibility was measured in the temperature range of 2–20K in steps as small as 0.1K per step as it approached the Curie temperature
.
Electric resistivity measurements were performed with the Triton 500 cryogenics system set up by the NTU-CCMS group, shown in
Figure 15 (left). A superconducting transition measurement for a NbN film sample is shown in
Figure 15 (right). Samples of 3 × 3 mm
sizes were prepared and glued on a sample holder with CMR-direct GE varnish. The sample was wire-bonded to readout pads on the sample holder using aluminum wires. The holder carried 20 readout pads, allowing us to perform more than the minimum requirement of 4 bonds. In this way, we ensured that we still had connectivity in case some bonds broke in very low operating temperatures. The resistivity was first measured at room temperature to check for possible oxidation or defects in the film growth process, and to test the connectivity of the wire bonds.
After the successful characterization of the NbN-sample superconducting properties, we proceed with the production of prototype nanowire sensors. The performance requirements for the Hawking photon sensors necessitate the use of the electron beam lithography (EBL) for the etching of nanowires from the NbN films. Currently, nanowire prototypes of different widths and lengths are under design. The baseline design using an autoCAD drawing of a 20 × 20
m
sensing area, with a nanowire with a width of 100 nm and a pitch of 100 nm, is shown in
Figure 16 and the zoom-in is shown in
Figure 17.
The SNSPD sensor prototypes are produced by the EBL, ELS-7000 ELIONIX, machine located in the Academia Sinica laboratories, as shown in
Figure 18 (Left). Mean while, our AnaBHEL Collaboration has purchased a Junsun Tech MGS-500 sputtering machine installed at the NEMS center of NTU
Figure 18 (Right), which will be fully utilized.
In order to maximize the single photon detection efficiency, various structures such as cavities or Bragg reflectors can be utilized. As part of the ongoing R&D, distributed Bragg reflectors (DBR) have been grown in the NTU MEMS facility. The measured reflectivity of a DBR for a sensor sensitive at 1550 nm is shown in
Figure 19. The good agreement with a finite-difference time-domain method (FDTD) simulation of the structure, gives us the confidence to proceed with new cavity designs optimal for longer wavelengths.
We are currently in the process of setting up a system test bench to characterize the Hawking sensors, using single photons at the infrared. The setup includes a SPAD commercial sensor for single photon calibration shown in
Figure 20. We plan to first verify the sensor operation at 1550 nm where most commercially available SNSPDs operate, as shown in
Figure 21. Finally, the sensors will be tested at longer wavelengths relevant to the AnaBHEL experiment shown in
Figure 22.
9. Experimental Backgrounds
The propagation of the high-intensity laser through a plasma target would necessarily induce background photons that would compete against the rare Hawking signals. The plasma electrons perturbed by the propagating laser would execute non-trivial motions and can therefore emit photons. In addition, they can interact with the electromagnetic fields induced by the laser and charged particles, and also with the plasma ions through scatterings.
The radiations induced from interactions between the electrons and the background ions can be categorized into Thomson/Compton scattering and Bremsstrahlung. These processes have long been well studied and the radiation so induced can be estimated when the electron trajectories are given.
There is also the possibility of radiation caused by electron acceleration. The analytic solution for plasma accelerating in the blowout regime of plasma wakefield excitations has been studied by Stupakov [
79], where it was shown that there are not only accelerated plasma inside the bubble but also charged particles that oscillate along the boundary of the plasma bubble. The work [
79] was for the case of the plasma wakefield accelerator (PWFA) [
45], but the method can also be applied to a laser wakefield accelerator (LWFA) [
44], which is the basis of our flying plasma mirror. Thus we also expect to have the same type of electron motions that are oscillating around the plasma bubble. These electrons in the plasma wakefields perform a
figure-8 motion in the plasma, and they can emit low-energy photons through synchrotron radiation. These photons are propagating in the direction parallel to the laser, which could affect the observation of the partner photons downstream. Therefore, we should study these electrons carefully.
In the following, we categorize the trajectories of the plasma electrons obtained from simulations by using a machine learning-based technique. We classify the electrons into several categories, according to their characteristic motions. After this classification, we are able to identify the leading radiation processes for the electrons and evaluate the radiation spectrum. We use SMILEI [
80] for particle-in-cell (PIC) simulations and python and the scikit-learn library [
81] for the clustering analysis.
9.1. Simulation Setup
The PIC simulations are in 2D and we refer to the coordinate as
x and
y. The simulation box size is 250
150
, i.e.,
, −75
75
, with
grids. A Gaussian laser with 800nm wavelength and
is applied at the boundary of
(left end of the simulation) and travels in the
x-positive direction. We place helium gas in the simulation box that can be ionized by the impinging laser. The helium density
is given by
where
and
. We do the simulation for 265 time steps, where each step is
femtoseconds in real-time.
9.2. Categorization of Electron Motions
Following the categorization technique introduced in [
82], we identify electron trajectories that would induce photons that dominate the background signals. The trajectory categorization introduced in [
82] is essentially clustering in momentum space using k-mean clustering method.
Let us denote the i-th particle’s trajectory by , where and are the x, y coordinate of the i-th particle at time t, respectively. The total time steps of the simulation are denoted as T. (In this case, T = 265). If we have N particles to track, then our data set will be . Let us denote the Fourier coefficient of and as and , respectively. The categorization will be done with the following steps.
- 1.
Restrict the tracked particle data to those that have been simulated for more than 380 femtoseconds.
- 2.
Prepare a data set,
where
and
are the mean of momentum of the
i-th particle in
x and
y direction, respectively,
is the mean of the
y-coordinate of the
i-th particle,
and
are the maximum and minimum of the acceleration in the
y-direction, and
is the corresponding frequency of the
t-th Fourier coefficient.
- 3.
Calculate k principal component values (PCVs) from the data set. This reduces the space of clustering from -dimensional vector space to k dimensional vector space.
- 4.
Perform k-mean clustering in the k-dimensional space, for a given number of clusters K.
Our choice of data set at step 2 is different from the one used in [
82], where we have additional value
and the information of its acceleration in the
y-direction. We added these since the longitudinal behavior is quite important for the experimental purpose, and indeed, by adding these we were able to separate the modes into more reasonable categories.
9.3. Classification Results
We have used , in the following, i.e., we classify the particles into 12 sets by using 30 PCs. Although we have classified them into 12 categories, since we have included the mean value of y coordinate, , into the data, we obtain pairs of categories that are almost symmetric along . In the following, we classify those two into the same category, since their physical processes are the same.
9.3.1. Wakefield Accelerated Electrons
The first kind is the electrons accelerated with the LWFA process. They are accelerated in the forward direction, to a highly relativistic regime
. These are shown in
Figure 23.
These electrons can radiate photons by interacting with the nuclei, i.e., through Thomson/Compton scattering or as Bremsstrahlung.
9.3.2. Snowplowed Electrons
Snowplowed electrons are the ones that are pushed forward by the laser’s ponderomotive potential and are clustered at the front of the laser pulse.
Figure 24 is a snapshot of the snowplowed electrons.
9.3.3. Backward Scattered Electrons
These electrons typically have
and are shown in
Figure 25.
They might contribute to the background radiation via Thomson/Compton scattering or Bremsstrahlung.
9.3.4. Slide-Away Electrons
There are certain fractions of plasma electrons that are pushed by the transverse plasma wakefields and propagate in the transverse direction. In practice, they would not affect the experiment since they are not moving toward the sensor, however, one would have to consider their hitting and reflection from the gas nozzle, which would induce background photon events.
Figure 26 is a snapshot of slide-away electrons. As pointed out previously, slide-away electrons slide toward the positive
y direction, and are classified into a different category through the process.
9.3.5. Transverse Oscillating Electrons
The last ones are the oscillating electrons. These are the electrons that are attracted by the Coulomb force of the plasma ion bubble and they oscillate around the laser trajectory in the traverse direction.
Figure 27 shows the density distribution, velocity, and acceleration in the
y-direction of these electrons.
This distribution has a tail that extends into the non-relativistic region. We expect that they would emit photons through synchrotron radiation, which can affect the identification of the Hawking photon signals.
9.3.6. Low-Frequency Soliton Radiations
We note that laser–plasma interaction also induces additional low-frequency background photons emitted by collective effects such as solitons, which are not included in the above discussion. Such near-plasma-frequency radiation was first pointed out by Bulanov [
83] and has been recorded experimentally [
35]. The radiation released by the solitons propagates essentially in the forward direction, where the seek-after
partner photon signals are expected to be in the EUV range. Thus such soliton-induced radiation signals may not render competing backgrounds to our experiment. Nevertheless, we will further investigate this collective soliton effect to determine whether some of such radiation might be reflected backward so as to confuse the Hawking photons, whose wavelengths are indeed close.
Our next step is to estimate the radiation with the corresponding process according to these categories, and compare the result with the radiation spectrum generated by the PIC simulation code and assess their impacts on the AnaBHEL experiment.
10. Strategy of AnaBHEL
We execute the AnaBHEL project based on the following strategy.
Stage-1
R&D of the key components, namely the superconducting nanowire single-photon Hawking detector and the supersonic gas jet with the designed density profile, are mainly carried out at the Leung Center for Cosmology and Particle Astrophysics (LeCosPA), National Taiwan University. These are going well, as reported in the previous sections.
Stage-2
Dynamics of the laser-induced plasma mirror trajectory and its correspondence with the plasma density profile. The first attempt was scheduled at Kansai Photon Science Institute (KPSI) in Kyoto, Japan, using its PW laser facility, in the summer of 2022. We expect that the iterative interplay between the gas jet design and the laser–plasma interaction data acquisition is indispensable.
Stage-3
The full-scale analog black hole experiment used to detect Hawking and partner photons will be pursued when the Hawking detector is fully developed and the plasma mirror trajectory is characterized. It is our desire that the Stage-3 experiment be carried out at the Apollon Laser Facility in Saclay, France.