2.1. Liquid target chamber setup for neutron generation

A 1-kHz Ti:sapphire laser was incident on a sheet of room-temperature, free-flowing deuterium oxide ( ${\mathrm{D}}_2\mathrm{O}$ ). On-target intensities upwards of $5\times {10}^{18}$   $\mathrm{W}/{\mathrm{cm}}^2$ were reached with 8 mJ of energy in 40 fs and a 1.65 μm (full width at half maximum (FWHM)) spot size, focused with an f/1.0 gold-coated off-axis parabolic mirror (OAP, Aperture Optical Science). The laser has a central wavelength of 780 nm and is incident on the target s-polarized with an angle of incidence of approximately 45°.

The target is a sub-micrometer-thick liquid flowing sheet formed from two intersecting 25 μm diameter ${\mathrm{D}}_2\mathrm{O}$ cylindrical jets. As discussed by George et al. ^[ Reference George, Morrison, Feister, Ngirmang, Smith, Klim, Snyder, Austin, Erbsen, Frische, Nees, Orban, Chowdhury and Roquemore^28], this sub-micrometer scale target is extremely stable and can operate at a kHz repetition rate or above. Two pumps, one for each jet, push the ${\mathrm{D}}_2\mathrm{O}$ through at a rate of 1 mL/minute. The approximately 65 mL size of the smaller of the two pumps thus limits the duration of an experiment to around an hour. Without this constraint, an experiment could run indefinitely. With roughly tens of nano-liters of ${\mathrm{D}}_2\mathrm{O}$ ionized per shot and the ability to recycle what remains, the target material costs roughly 2 US dollars per minute of run-time.

A second laser, temporally locked to the main laser, is used to image the interaction region in a pump–probe scheme. The energy in this frequency-shifted probe (80 μJ) is significantly lower than that of the main laser pulse. The probe has a 420 nm central wavelength and 80 fs pulse duration^[ Reference Feister, Nees, Morrison, Frische, Orban, Chowdhury and Roquemore^31]. A simplified chamber diagram can be seen in Figure 1 and a 3D rendering is shown in Figure 2. The probe beam passes through the target and the microscope objective to be imaged onto a camera, allowing real-time video diagnostics with temporal resolution of approximately 50 fs and spatial resolution of 1 μm. This video is the primary diagnostic that is used to align the laser–target system and maximize energy into the sheet.

Figure 2 Three-dimensional model of the interaction region. As in Figure 1, the main beam is shown in red and the probe beam is shown in blue. The jets are in white near the top of the image, and the size of the liquid sheet is exaggerated for illustration.

The 107-cm diameter stainless steel chamber is brought to a final vacuum of approximately 1 Torr, which is limited by the vapor pressure of heavy water and the tendency of the liquid target to freeze. This is below the approximately 7-Torr threshold for ion acceleration noted by Snyder et al. ^[ Reference Snyder, Morrison, Feister, Frische, George, Le, Orban, Ngirmang, Chowdhury and Roquemore^30]. This pressure is measured at the edge of the chamber; it is expected that the pressure in the vicinity of the target may be significantly higher than the stated 1 Torr. The distance between the target and pressure transducer is approximately 53 cm.

2.2. Laser-induced deuterium–deuterium fusion

When ultra-intense laser pulses interact with a deuterium-rich target, two processes can give rise to neutron production^[ Reference Willingale, Petrov, Maksimchuk, Davis, Freeman, Joglekar, Matsuoka, Murphy, Ovchinnikov, Thomas, Van Woerkom and Krushelnick^13], from the bulk and from a deuterium-rich catcher placed at the back of the target (the so-called pitcher–catcher scheme). Neutron production in the bulk may be categorized by two general processes. Firstly, as the absorbed laser energy is transferred from hot electrons to ions in the bulk, the local temperature at and near the focal region may become very high, which may cause deuterium–deuterium (D-D) fusion. Secondly, a significant portion of the hot target then explodes after some time, where the exploding high-energy deuterons collide with each other, causing D-D fusion events (a few nanoseconds after the pulse leaves the target; see, for example, the supplementary movie target explosion dynamics captured for a liquid target in Ref. [Reference Ngirmang, Morrison, George, Smith, Frische, Orban, Chowdhury and Roquemore32] for similar intensities). At our laser intensities of high ${10}^{18}\;\mathrm{W}/{\mathrm{cm}}^2$ , energetic deuterons are accelerated outward from the surface primarily via target normal sheath acceleration (TNSA)^[ Reference Hatchett, Brown, Cowan, Henry, Johnson, Key, Koch, Langdon, Lasinski, Lee, Mackinnon, Pennington, Perry, Phillips, Roth, Sangster, Singh, Snavely, Stoyer, Wilks and Yasuike^33, Reference Morrison, Storm, Chowdhury, Akli, Feldman, Willis, Daskalova, Growden, Berger, Ditmire, Van Woerkom and Freeman^34], which may then impinge on a nearby secondary deuteron-rich target. At energies of approximately keV and above, colliding deuterons can fuse together (D-D fusion), with half of the fusion reactions producing ${}^3\mathrm{He}$ and a neutron and the other half producing tritium and a proton:

(1) $$\begin{align}^2\mathrm{H}+{}^2\mathrm{H}\to {}^3\mathrm{H}\mathrm{e}+\mathrm{n}+3.27\ \mathrm{MeV},\nonumber\\^2\mathrm{H}+{}^2\mathrm{H}\to {}^3\mathrm{H}+{}^1\mathrm{H}+4.03\ \mathrm{MeV}.\end{align}$$

In the neutron producing branch, the 3.27 MeV Q-value is distributed as kinetic energy of the two products. With roughly a quarter of the total mass, the free neutron will take roughly three quarters of the energy, or 2.45 MeV. Our experiment is designed to detect neutrons from these D-D fusion events. We also provide evidence in Section 3.2.1 that the neutron energies are 2.45 MeV, as expected.

In the center-of-mass frame, neutrons are emitted without a directional bias. The approximately $5\times {10}^{18}\;\mathrm{W}/{\mathrm{cm}}^2$ average laser intensity of this experiment implies a maximum energy of accelerated deuterons in the 0.1–1 MeV range, which is not expected to produce significant anisotropy in the angular distribution of neutrons.

2.3. Suite of three neutron detection systems

Three independent detection systems are used to verify the generation of neutrons: an EJ-309 organic liquid scintillator (Eljen Technology) coupled to a photomultiplier tube (Hamamatsu, R7724), a ${}^3\mathrm{He}$ proportional counter (Reuter-Stokes) and a set of 36 bubble detectors (Bubble Tech Industries, BDS). The bubble detectors are the fastest to setup but have certain uncertainties in measurement, the EJ-309 has a fast response that could help in measuring the neutron energy and the ${}^3\mathrm{He}$ detector is highly efficient for thermal neutron detection and is inherently gamma-blind^[ Reference Kouzes, Ely, Lintereur and Stephens^35, Reference Mauri, Messi, Kanaki, Hall-Wilton and Piscitelli^36]. Our suite of the three detectors is utilized for high measurement confidence in neutron counts and the results are compared in Section 3.1.

EJ-309 organic scintillator: the EJ-309 scintillator consisted of a 5.08 cm right-hand circular cylindrical liquid cell in a thin aluminum housing. This cell was coupled to a 5.08-cm diameter PMT via a borosilicate glass window and EJ-550 silicone grease. The scintillator and PMT detector system, housed in 1.3-cm thick bismuth box with an additional 1.3-cm thick lead sheet at the front of the bismuth box, was placed 150 cm from the target interaction area. Both photons and neutrons were measured during laser operation even with the shielding in place, indicating an extremely active photon source. Photon and neutron events can be separated in the analysis via PSD due to the different scintillation decay profiles created by recoil electrons (corresponding to photon interactions) and recoil protons (corresponding to neutron interactions). The EJ-309 has good discrimination between photons and neutrons, even in a high gamma-ray environment^[ Reference Stevanato, Cester, Nebbia and Viesti^37]. The detector was biased to –1300 V and events were analyzed with a waveform digitizer (CAEN Technologies, DT5730) and CoMPASS software.

In post-processing, typically a scintillation light-yield threshold is set, below which the neutrons and photons are indistinguishable. Such below-threshold events are discarded from the analysis. A GEometry ANd Tracking (Geant4)^[ Reference Allison, Amako, Apostolakis, Araujo, Dubois, Asai, Barrand, Capra, Chauvie, Chytracek, Cirrone, Cooperman, Cosmo, Cuttone, Daquino, Donszelmann, Dressel, Folger, Foppiano, Generowicz, Grichine, Guatelli, Gumplinger, Heikkinen, Hrivnacova, Howard, Incerti, Ivanchenko, Johnson, Jones, Koi, Kokoulin, Kossov, Kurashige, Lara, Larsson, Lei, Link, Longo, Maire, Mantero, Mascialino, McLaren, Lorenzo, Minamimoto, Murakami, Nieminen, Pandola, Parlati, Peralta, Perl, Pfeiffer, Pia, Ribon, Rodrigues, Russo, Sadilov, Santin, Sasaki, Smith, Starkov, Tanaka, Tcherniaev, Tome, Trindade, Truscott, Urban, Verderi, Walkden, Wellisch, Williams, Wright and Yoshida^38] Monte Carlo radiation transport simulation of the experiment is used to determine the absolute neutron detection efficiency of the EJ-309 scintillator as a function of the light-yield threshold, accounting for geometric effects. In Geant4, the neutrons are modeled as coming from a point source with the scintillator at a distance of 150 cm. Neutrons are emitted isotropically in the Geant4 simulation. Figure 3 shows the simulated estimated absolute efficiency both with and without the environment modeled, highlighting negligible environmental scattering effects above a 0.4 MeVee (MeVee, MeV electron equivalent) light-yield threshold.

Figure 3 Absolute detection efficiency of the EJ-309 scintillator calculated via Geant4 simulation. The vertical axis indicates that roughly ${10}^5$ neutrons are produced for every neutron detected, and the shaded regions indicate error. Light-yield thresholds are equal to or more than 0.4 MeVee, and the difference in the two efficiencies is less than 0.5%. The units are defined such that 1 MeVee (MeV electron equivalent) equals the number of scintillation photons produced by a 1 MeV electron.

A ${}^{137}\mathrm{Cs}$ source is used for the energy calibration of the EJ-309, with a known Compton edge at 478 keV^[ Reference Swiderski, Moszyński, Czarnacki, Iwanowska, Syntfeld-Każuch, Szczęśniak, Pausch, Plettner and Roemer^39]. This value is chosen as the light threshold for the detection efficiency: matching the calibration value minimizes error from an imperfect calibration, and environmental effects at this threshold are negligible. However, multiple-scatter events cloud the Compton edge, and without a Monte Carlo simulation as detailed by Dietze and Klein^[ Reference Dietze and Klein^40], its location cannot be precisely known. In addition, these calculations are based on the MeVee unit, which is flawed in the context of fast neutron detection in organic scintillators due to an inherent assumption of proportionality^[ Reference Laplace, Goldblum, Brown and Manfredi^41]. In calibrating an organic scintillator based on the Compton edge of a known source, care must be taken to minimize these issues.

³He counter: this detector relies on thermal neutron capture on ${}^3\mathrm{He}$ , which creates a proton and a triton that carry the 764 keV Q-value of the reaction as kinetic energy in opposite directions, as given by the following:

(2) $$\begin{align}{}^3\mathrm{H}\mathrm{e}+{\mathrm{n}}_{\mathrm{th}}\to {}^1\mathrm{H}+{}^3\mathrm{H}+764\ \mathrm{keV}.\end{align}$$

A ${}^3\mathrm{He}$ proportional counter leverages the high cross-section of this thermal neutron capture reaction as well as the favorable properties of ${}^3\mathrm{He}$ as a fill gas. When neutron capture occurs, the proton and triton ionize the ${}^3\mathrm{He}$ gas and the created charge is collected via an electric field from an applied bias voltage. This results in a signal with energy that corresponds to the reaction Q-value. As the ${}^3\mathrm{He}$ counter is most sensitive to thermal neutrons (thermal neutron cross-section, 5330 barns^[ Reference Mauri, Messi, Kanaki, Hall-Wilton and Piscitelli^36]), moderation is needed to detect fast neutrons with high efficiency. This moderation must then be carefully modeled in Geant4 to attain an accurate efficiency estimation. A 2.5 cm by 81.3 cm cylindrical ${}^3\mathrm{He}$ proportional counter was positioned approximately 140 cm from the target position. Blocks of paraffin wax were placed in front of the detector for neutron moderation. The detector was biased to 1100 V and signals were recorded using Maestro software.

The integral of the full-energy peak (FEP) at $Q=764$ keV gives the number of detected neutrons. Not all fast neutrons are moderated by the paraffin, and some may elastically scatter off of the ${}^3\mathrm{He}$ nuclei and deposit a wide range of energies. To eliminate these spurious counts, which underlie the FEP, the FEP region is fit to a Gaussian and exponential decay function. The integral of the Gaussian, corresponding to full-energy capture, yields the detected events. The counts in the exponential, corresponding to neutron scatter, are discarded. As the ${}^3\mathrm{He}$ tube is sensitive to neutrons generated by cosmic rays, each experiment is run alongside a natural background measurement, which is then subtracted from the main data, hereafter referred to as background-subtracted data. A Geant4 simulation is used to estimate the ${}^3\mathrm{He}$ absolute neutron detection efficiency, accounting for environmental scattering from the vacuum chamber, EJ-309 detection system, lead and bismuth shielding and low-Z moderating blocks. All other objects in the room were neglected; the model does not account for the room-return impact or additional moderation of the neglected objects. The Geant4 simulated absolute efficiency is $1.53\times {10}^{-4}\pm 0.3\%$ : for every neutron detected, $65,450\pm 170$ neutrons are emitted from the target into 4 $\pi$ .

Bubble detectors: the bubble detector spectrometer (BDS) is a set of 36 detectors rated to measure neutrons above six different energy thresholds: 0.01, 0.1, 0.6, 1, 2.5 and 10 MeV, with six detectors at each threshold. The bubble detectors boast no photon detection and minimal ion/electron sensitivity. In each detector, a polymer gel suspends millimeter-sized super-heated liquid droplets. As a neutron passes through the gel, it deposits its energy into recoil ions; these ions then may pass through a super-heated liquid drop, which quickly vaporizes and expands into a visible bubble^[ Reference Lewis, Smith, Ing, Andrews, Machrafi, Tomi, Matthews, Veloce, Shurshakov, Tchernykh and Khoshooniy^42].

One bubble detector of each energy threshold is placed in a group, with six groups attached at various positions directly to the outside of the chamber. After the experiment, bubbles were counted by eye as a measure of the neutron count, using the bubbles/neutron sensitivity measured by Lewis et al. ^[ Reference Lewis, Smith, Ing, Andrews, Machrafi, Tomi, Matthews, Veloce, Shurshakov, Tchernykh and Khoshooniy^42]. The bubbles can then be compressed, allowing the detectors to be reused. The bubble detectors in these experiments were used 16 times over several months, although they had first been activated two years prior.

The measured data of the neutrons from the laser-driven fusion source and their comparisons are described in the following section.

3.1. Observation with neutron detection suite

EJ-309 organic scintillator: Figure 4 shows the 2D PSD histogram from the EJ-309 liquid scintillator. The PSD metric on the y-axis is the ratio between the integral of the tail of the scintillation event’s pulse to the total pulse integral from scintillation. The x-axis is given by the total pulse integral. Recoil protons from neutron interactions result in more delayed scintillation light compared to recoil electrons from photon interactions, and thus neutron counts have a higher PSD value than photon counts. As such, two separate features form in Figure 4, corresponding to neutrons at the higher-PSD cluster and photons in the lower-PSD cluster. The data in Figure 4 are mapped to a 1D histogram of PSD values in Figure 5. To mitigate environmental scattering effects and avoid misclassifying neutron and gamma-ray signals, a 478 keVee (keVee, keV electron equivalent) light-yield threshold is used for all experimental measurements.

As a control, natural (undeuterated) water is tested under the same conditions, as laser pulses on ${\mathrm{H}}_2\mathrm{O}$ should not produce neutrons. Figure 5 shows the PSD histogram of the ${\mathrm{H}}_2\mathrm{O}$ target in red. No significant neutron feature is present. The ${\mathrm{H}}_2\mathrm{O}$ experiment was half the duration of the ${\mathrm{D}}_2\mathrm{O}$ experiment, causing the size difference of the photon peaks.

Figure 4 Two-dimensional PSD histogram of neutrons and photons in the organic scintillator. The color scale denotes the number of events.

Figure 5 One-dimensional PSD histograms for the EJ-309 scintillator. The blue data are the same as in Figure 4; the red curve shows a shorter-duration experiment with ${\mathrm{H}}_2\mathrm{O}$ instead of ${\mathrm{D}}_2\mathrm{O}$ as a control.

³He counter: Figure 6 shows the background-subtracted data from the ${}^3\mathrm{He}$ proportional counter, comparing the ${\mathrm{D}}_2\mathrm{O}$ (blue) and ${\mathrm{H}}_2\mathrm{O}$ (red) experiments. For ${\mathrm{D}}_2\mathrm{O}$ , a sharp peak corresponding to the thermal neutron capture energy (764 keV) is seen along with a large number of counts at lower energies. For ${\mathrm{H}}_2\mathrm{O}$ , few total counts are recorded with no significant events near the thermal capture energy. Several hundred counts are seen at low energies, most likely due to the X-ray environment from the laser–plasma interactions. This large discrepancy in low-energy counts indicates that the signal at such energy in the ${\mathrm{D}}_2\mathrm{O}$ data is a result of neutron scatter, which is accounted for using the fitting method described in Section 2.3.

Figure 6 Background-subtracted data from the ${}^3\mathrm{He}$ tube. ${\mathrm{D}}_2\mathrm{O}$ (blue) and ${\mathrm{H}}_2\mathrm{O}$ (red) correspond to the same experiments shown in Figure 5.

Bubble detectors: in the analysis of the BDS, the 2.5 and 10 MeV bubble detectors were neglected, as their response to 2.45 MeV neutrons is not well characterized. With each individual detector typically exhibiting 15 or fewer bubbles over an hour of run-time, the set of 24 remaining bubble detectors showed hundreds of bubbles after an experiment with ${\mathrm{D}}_2\mathrm{O}$ . In contrast, over two separate ${\mathrm{H}}_2\mathrm{O}$ experiments totaling 102 minutes of run-time, only five bubbles in total were created. Normalizing by time, this corresponds to 0.05 bubbles/minute for ${\mathrm{H}}_2\mathrm{O}$ as opposed to 5.5 bubbles/minute from a representative ${\mathrm{D}}_2\mathrm{O}$ experiment.

Both the proportional counter and organic scintillator require analysis to eliminate spurious counts: the EJ-309 requires PSD to separate particle types, and the proportional counter requires fitting techniques to remove scattering events. With greater shielding on the EJ-309 and more moderation on the proportional counter, these unwanted events should be reduced, making analysis easier and more accurate. However, increased moderation may not be possible in all applications because of geometry constraints; for example, large amounts of wax are needed surrounding the 85-cm tube. By comparison, the EJ-309 (with PMT and Bi+Pb shielding) is small and additional lead shielding can easily be placed in front of the scintillator, blocking the line-of-sight from the source. Also, the ${}^3\mathrm{He}$ counter is more sensitive to unavoidable room-return effects from extraneous equipment, walls and the floor due to its dependence on thermal neutrons.

The BDS is convenient and easy to use because of its insensitivity to photons; hence, no analysis is needed to distinguish neutrons from other events. To obtain a neutron count from these bubbles even a non-expert can count the bubbles using calibration data provided by the manufacturer. Their small size allows them to be placed almost anywhere; also, they are inexpensive compared to the other two detectors. However, temperature has a significant effect: operating even a few degrees above the recommended temperature of 20 ${}^{\circ }$ C causes an increase in detection efficiency and a decrease in the energy threshold, an effect that is not well characterized in the literature^[ Reference Buckner, Noulty and Cousins^25]. Several hours of compression were needed to reset the bubble detectors for reuse, which may limit how many experiments can be performed. In addition, the bubble detectors have a limited lifespan and their sensitivity shifts after repeated use^[ Reference Vanhavere, Loos and Thierens^26]. Ultimately, the bubble detectors need special care to obtain quantitative measurement of neutrons, especially, maintaining the temperature and keeping track of the measurement history.

The BDS typically measured neutrons in the high ${10}^5$ to low ${10}^6$ neutrons/second range. Due to the close proximity of bubble detectors within a group, it is possible that a neutron may be counted more than once as it scatters multiple times. This may be one reason for this high flux estimation by the BDS, which was significantly higher than the flux recorded by the other two detection systems: the ${}^3\mathrm{He}$ proportional counter typically saw neutrons/second in the low ${10}^5$ regime, with the EJ-309 in the high ${10}^4$ regime. Table 1 displays the neutron counts and uncertainties for a selection of two experiments when all detectors were used simultaneously. The reported error in the BDS is only the statistical (counting) uncertainty, whereas the error in the other two detectors includes both statistical uncertainty and the uncertainty in the Geant4 efficiency calculations.

Table 1 Comparison of neutrons/second and associated statistical errors from the three detection systems across two separate experiments. Note that other non-statistical errors contribute to the uncertainties and are not represented in this table, in particular for the ${}^3$ He detector and bubble detectors.

The unknown shifts in the efficiency of the BDS, not represented in the reported error of Table 1, make it unsurprising that the results do not closely match the results from the other two detectors. However, the scintillator and proportional counter should agree on neutron flux; instead, the ${}^3\mathrm{He}$ counter records a neutron flux two to four times higher than the EJ-309, depending on the experiment. The Geant4 modeling of the ${}^3\mathrm{He}$ counter is likely a large source of this discrepancy due to the detector’s sensitivity to thermal neutrons. Only modeled are the experiment chamber, detectors and wax moderation: all other equipment in the room that may contribute to neutron thermalization is neglected. The uncertainty introduced by not including these features is difficult to quantify and is not included in the error analysis above. Comprehensive modeling of the neutron thermalization process and subsequent thermal neutron scatters is known to be a challenge, in part due to the many possible contributions to thermalization and in part due to deficiencies in knowledge of relevant cross-sections^[ Reference Tran, Marchix, Letourneau, Darpentigny, Menelle, Ott, Schwindling and Chauvin^43]. With a greater amount of thermalization in the experiment than what is modeled, the Geant4 analysis underestimates the absolute efficiency and thus overestimates the total neutrons emitted. Upon adding a thin floor and walls to the model, the estimated neutron emission by the ${}^3\mathrm{He}$ detector was reduced by roughly half, bringing it closer to the predicted flux of the EJ-309. The EJ-309 is insensitive to thermal neutrons and therefore does not depend so strongly on details of the Geant4 model.

Before detection efficiency and solid angle considerations, the proportional counter saw the highest number of raw counts, allowing for lower relative statistical error. However, the significant difference in relative error ( $\sim$ 0.5% for the proportional counter as opposed to $\sim$ 15% for the scintillator) is a result of different considerations being included. The uncertainty in the ${}^3\mathrm{He}$ efficiency only includes the statistical Monte Carlo uncertainty that is standard in Geant4, whereas the uncertainty in the EJ-309 efficiency also includes light-yield uncertainty, as detailed in the literature^[ Reference Bai, Wang, Zhang, Lu, Jiang, Chen and Zhang^44– Reference Egner, Febbraro, Holland and Bevins^46]. In addition, the uncertainty in the EJ-309 efficiency estimation overcomes thermalization uncertainty or environmental neutron scattering effects by using a light-yield threshold of 478 keVee, while the ${}^3\mathrm{He}$ detector response is unable to overcome this issue.

To summarize, the EJ-309 detector likely provides the most reliable measurement of neutron production. Although it is easiest to infer the neutron numbers from the BDS, the results are found to be inconsistent as the detectors have unknown dependencies on many parameters while suffering from poor counting statistics. The ${}^3\mathrm{He}$ proportional counter is too reliant on a detailed environmental model, which can be challenging to develop for a laboratory setup. In comparison, the EJ-309 detectors are well characterized in the literature and they have minimal environmental effects at energy thresholds above 0.4 MeVee.

3.2. Neutron characterization

Next, the neutron energy was measured using the time-of-flight (TOF) between the particle detection pulse at the EJ-309 scintillator and the laser pulse incident on the target. The expected energy is 2.45 MeV. Then, several experimental parameters are varied: the chamber pressure and pre-pulse effects on the neutron yield are studied, and the spatial distribution of the source is measured in a pitcher–catcher scheme.

3.2.1. 2.45 MeV neutrons

With a known distance from the neutron source to the scintillator, the time delay between the impingement of the pulse on the target and the arrival of neutrons can be used to determine the neutrons’ energies. The time at which the neutrons are produced is first estimated by the laser trigger signal, and then corrected with the arrival of the photons in the scintillator, as they travel at a known speed. This TOF analysis is shown in Figure 7. Relativistic forms of all equations are used. The detector electronics bin all events into 4-nanosecond windows, causing the discrete energy data in Figure 7. Still, a sharp peak in energies is seen around the expected 2.45 MeV, confirming D-D fusion.

Figure 7 Energy histogram of emitted neutrons, as measured by the organic scintillator via time-of-flight. Data were collected for 54 minutes with the scintillator subtending 0.0035 steradians, and counts are not scaled with detection efficiency. That is, only the raw counts are shown. The peak corresponds to the expected 2.45 MeV of D-D fusion neutrons.

Detected neutron counts with lower energy are likely a result of neutron scattering from laboratory surrounding features, such as the room’s floor and walls, as well as the nearby paraffin wax surrounding the ${}^3\mathrm{He}$ detector. The low-energy counts may also, in part, be due to deuteron breakup, as high-energy ( $\gtrsim$ 2.2 MeV) deuterons strike nearby materials. In addition, neutrons are emitted with 2.45 MeV in the center-of-momentum (COM) frame: if the initial deuterons have a bias in initial momentum, the scintillator in the lab frame will measure the neutron as faster or slower. Based on Figure 7 (for the configuration in Figure 1), our measurements strongly indicate that fusion neutrons are generated as the ${\mathrm{D}}_2\mathrm{O}$ liquid sheet target is expanding both in the forward and backward directions. The distribution of these COM frames gives a Gaussian profile to the neutron energy peak.

3.2.2. Anisotropy and spatial resolution

To further characterize our neutron source, we examine its angular distribution. Although the neutrons are expected to be produced isotropically, a catcher could introduce a directional bias; the deuterons in a solid ${\mathrm{D}}_2\mathrm{O}$ catcher should function as additional targets for the laser-accelerated deuterons. Placed immediately behind the liquid sheet, deuterons striking the catcher now have a bias in initial momentum, causing a potential bias in the direction of the overall neutron flux. However, as found by Willingale et al. ^[ Reference Willingale, Petrov, Maksimchuk, Davis, Freeman, Joglekar, Matsuoka, Murphy, Ovchinnikov, Thomas, Van Woerkom and Krushelnick^13] with a plastic target and a catcher, catchers are inefficient neutron sources: their low temperature increases the stopping power and makes fusion less likely. Our solid ${\mathrm{D}}_2\mathrm{O}$ catcher (heavy ice) did not significantly affect the isotropy.

Three identical EJ-309 detectors and PMTs were placed at different viewing angles around the chamber, all in the horizontal plane. Each scintillator was shielded from gamma-rays either with lead bricks or a bismuth container. The results are seen in Table 2: no detector measured significantly higher neutron generation than any other, indicating the catcher’s failure to contribute to neutron yield anisotropy. In addition, when compared to an identical experiment without a catcher, no significant difference in total neutron yield was observed. This is consistent with the observation of Ref. [Reference Willingale, Petrov, Maksimchuk, Davis, Freeman, Joglekar, Matsuoka, Murphy, Ovchinnikov, Thomas, Van Woerkom and Krushelnick13], where bulk neutron production was shown to exceed neutron production by the pitcher–catcher method for intensities near the range of intensities described here.

Table 2 Detector angles, estimated source counts and calculated uncertainties for the three-scintillator array. Detectors are named by their viewing angle, with ${0}^{\circ }$ defined as the laser propagation direction.