Introduction

A large energy, high-power, high-intensity laser can create a large volume of ultra-high-energy-density plasma for particle accelerators1, planetary science2,3, astrophysics4,5, and nuclear physics6,7. Inertial confinement fusion (ICF) is an ultimate application of such extreme plasmas8,9,10.

A few millimeter-scale spherical capsule used in the ICF contains a deuterium–tritium (DT) fusion fuel ice layer. In the laser indirect-drive approach, the capsule surface is irradiated by X-rays in a high-Z metal enclosure (hohlraum). The X-rays drive a sequence of converging shock waves and compressing the fusion fuel  > 1000 times solid density. Adiabatic compression heats up the DT gas initially filling the capsule interior, which becomes the ignition spark at the final stage of the compression.

The scientific breakeven, where the energy released by fusion reaction exceeds the energy contained in the compressed fusion fuel, was achieved in National Ignition Facility (NIF)11; however, the pathway to ICF ignition is still unclear. One of the crucial problems of the current central ignition scheme is the hot spark mixing with the cold dense fuel because of significant growth of hydrodynamic instabilities during the compression.

Fast isochoric heating, also known as fast ignition12, of a pre-compressed core, was proposed as an alternative approach to the ICF ignition that avoids the ignition quench caused by the mixing because the hot spark is generated not by the adiabatic compression but by the external energy injection whose timescale is shorter than the hydrodynamic timescale (<100 ps). Relativistic intensity laser pulses (\(I_{\mathrm{L}}\lambda _{\mathrm{L}}^2\) > 1.37 × 1018 W μm2 cm−2) efficiently produce relativistic electron beams (REB) via laser–plasma interactions13,14,15; here, IL and λL are laser intensity [W cm−2] and wavelength [μm], respectively. The REB travels in a plasma from its generation zone to the core. A part of the REB kinetic energy is deposited into the core, and then the heated region becomes the hot spark to trigger the fusion ignition.

Several breakthroughs were described in past articles on fast heating research: the invention of the cone-in-shell target16, high areal density core formation with the cone-in-shell target17, and efficient laser-to-core coupling after reduction of pre-plasma filling in the cone18. Our achievement, namely enhanced laser-to-core energy coupling with the magnetized fast isochoric heating enabled by an application of external kilo-Tesla-level magnetic field, is also a milestone as a demonstration of the efficient energy coupling, as a means to reduce the large radial spread of the REB. This is considered essential to secure scalability of the fast heating to the ignition19.

Three critical problems have been recognized as obstacles to the efficient fast heating with the REB. The first problem is that the REB generated in a long-scale-length pre-plasma filled in the cone becomes too energetic to heat the core, the second one is that a part of the REB is scattered and absorbed in a high-Z cone tip20. The third one is that the REB has a large divergence of 100° as a typical full-angle, so that only a small fraction of the diverged REB can collide with the core21.

The long-scale-length pre-plasma filling the cone is produced by prepulse and foot pulses of the heating laser and also by the cone breakup due to high pressure of a plasma surrounding the cone. The heating efficiency of 7% was attained at OMEGA laser facility with an ignition-scale large areal density (ρR ~0.3 g cm−2) core17 after reducing the pre-plasma formation in the cone18; however, the aforementioned second and third critical problems still remain to be resolved.

Here, we have introduced two experimental techniques. A solid ball target is used for making an open-tip cone utilizable along with the plasma compression. The solid ball compression does not generate shocks and rarefactions traveling ahead of the dense shell; therefore, a closed tip is not required for preventing a hot plasma from filling inside the cone. In addition, a relatively cold and dense core can be produced stably by using the solid ball target. The cold core enables us to visualize REB transport region18 in a dense core by using monochromatic Cu-Kα imaging technique without a significant energy shift of the Cu-Kα photon energy due to ionization of Cu atoms. The use of a solid ball target has another benefit for producing a moderate guiding field as discussed later.

Another technique is to use a laser-driven capacitor–coil target22 to generate a strong magnetic field. The strength of the magnetic field was measured on GEKKO-XII23,24, LULI200025, Shengguang-II26, and OMEGA-EP laser facilities27,28.  As discussed in the Supplementary Discussion and summarized in Supplementary Table 1, experimental results revealed that a 600–700-T magnetic field was generated by using a tightly focused kilo-joule and nano-second infrared (λL = 1.053 μm) laser beam.

Application of external magnetic fields to the path of a REB is expected to guide the diverging REB to the dense core19. For example, the gyroradius of a 1-MeV electron under the influence of a 1-kT magnetic field is 5 μm, which is smaller than a typical core radius (20 μm); therefore, a kilo-Tesla-level magnetic field is sufficient to guide the REB to the core. The guidance of the REB by the laser-produced external magnetic field has already been demonstrated experimentally in an uncompressed-planar geometry at LULI2000 facility (Ecole Polytechnique, France)29.

In this study, we have demonstrated the enhancement of the laser-to-core energy coupling with the magnetized fast isochoric heating. The maximum coupling of 7.7 ± 1.2% was achieved even with a relatively small radial areal density core (ρR = 0.08 g cm−2), which could be higher than that has been achieved by the fast isochoric heating without applying a magnetic field18. This value is fairly consistent with the simple evaluation, which reveals that a higher areal density core leads to higher laser-to-core coupling.

Results

Laser-to-core coupling measurement

The laser-to-core energy couplings were experimentally measured by systematically varying the experimental conditions: heating laser energy, injection timing of the heating laser, application or not of the external magnetic field, and open- or closed-tip cones.

Table 1 summarizes laser-to-core coupling efficiency obtained in this experiment. The data are separated into three groups by the rules according to the experimental conditions, where the external magnetic field was applied or not, and the cone-tip was open or closed. The data are sorted by the laser-to-core coupling efficiency in each group.

Table 1 Summary of laser-to-core coupling efficiencies

The coupling was calculated from the absolute number of Cu-Kα X-ray (8.05 keV) photons emitted from Cu atoms contained in the laser-compressed hydrocarbon core. Cross-sections of electron-impact K-shell ionization have a similar dependence on electron energy as collisional energy loss. The two are essentially the same process but with a different threshold energy. Collisional deposition of REB energy (J) in a Cu-contained core can be obtained from the number of Cu-Kα photons (photons/sr) emitted from the core18 with a correlation factor. Details of the correlation factor derivation are described in the Methods section.

Figure 1a shows an experimental layout. The experiment was conducted on the GEKKO-LFEX laser facility at the Institute of Laser Engineering, Osaka University. The fusion fuel surrogate was made of a 200-μm-diameter solid Cu(II) oleate solid ball [Cu(C17H33COO)2]30, whose surface was coated with a 25-μm-thick polyvinyl alcohol (PVA) layer to prevent the Cu atoms from being ionized directly by the compression laser beams. The Cu(II) oleate ball contains 9.7% Cu atoms in weight for visualization of the relativistic electron beam (REB) transport in the core and for measurement of the laser-to-core energy coupling (9.7% is a little bit smaller than the expected value (10.1%) because of a contamination inclusion in the Cu(II) oleate). The fuel surrogate was attached to a Au cone, whose open angle, wall thickness, and tip diameter were 45°, 7 μm, and 100 μm, respectively. The outer surface of the Au cone was coated with a 50-μm-thick PVA layer to delay the cone breakup time31. Open-tip or closed-tip Au cones were used in the experiments; the tips of the closed cones were covered with a 7-μm-thick Au layer.

Fig. 1
figure 1

Experimental layout for the laser-to-core energy coupling measurement. a A schematic drawing of the experimental layout for the magnetized fast isochoric heating. Three and six of GEKKO-XII laser beams were used for generation of the magnetic field using the capacitor–coil target and compression of the solid ball, respectively. Four LFEX laser beams were irradiated on the tip of the cone to produce a REB. The X-ray image using a spherically bent crystal and the X-ray spectrometer using a HOPG crystal were used for visualization of the REB transport and measurement of the laser-to-core energy coupling, respectively. b Timing chart of the magnetic-field-generation laser (pink broken line), fuel compression laser (green solid line), and laser-produced magnetic field (blue circular marks) pulses. The error bars show the uncertainty of magnetic field measurement described in ref. 24. The hatching area indicates the injection timing of the heating laser

The solid ball was compressed by six of GEKKO-XII laser beams arranged in a quasi-cylindrical geometry, whose wavelength, pulse shape, pulse duration, and energy were 0.526 μm, Gaussian, 1.3 ns full-width at half-maximum (FWHM), and 240 ± 15 J per beam, respectively. The center of the nickel-made coil having a 500-μm diameter was located 230 μm from the center of the ball to apply a strong magnetic field near the REB generation zone and the solid ball. The first disk of the laser-driven capacitor was irradiated through the hole of the second disk by three tightly focused GEKKO-XII laser beams, whose wavelength, pulse shape, pulse duration, and energy were 1.053 μm, Gaussian, 1.3 ns (FWHM), and 600 ± 20 J per beam, respectively, yielding 7 × 1015 W cm−2 of the peak intensity.

The tip of the cone was irradiated to produce a REB by four LFEX laser beams, whose wavelength, pulse shape, and pulse duration were 1.053 μm, Gaussian, and 1.8 ± 0.3 ps (FWHM), respectively. The total energy of the four LFEX beams on the tip was varied from 630 to 1520 J. The focal spot diameter was 50 μm (FWHM) containing 30% of the total energy, yielding an intensity of 1.3 × 1019 W cm−2 at the maximum energy shot.

The laser energy and injection timing of the heating laser are also summarized in Table 1. The injection timings were measured with an X-ray streak camera with an accuracy of ± 0.02 ns. Figure 1b shows a time chart of the magnetic-field-generation laser, compression laser, and laser-produced magnetic field pulses. The time origin (t = 0 ns) is defined hereafter as the peak of the compression laser pulse. The peak of the magnetic-field-generation laser pulse was set at t = −1.5 ns; therefore, the magnetic field strength reaches its maximum value before the compression beam irradiation. The heating lasers (four LFEX beams) were injected around the maximum compression timing t = 0.38–0.72 ns shown as the hatching area in Fig. 1b.

Copper Kα X-rays (8.05 keV) were imaged using a spherically bent quartz (2131) crystal to visualize the transport of the REB in the pre-compressed core from the direction perpendicular to the LFEX incident axis. The magnification, spatial resolution, and spectral bandwidth were 20, 13 μm (FWHM), and 5 eV (FWHM), respectively.

The X-ray spectrometer using a planar highly oriented pyrolytic graphite (HOPG), was installed at 40° from the LFEX incident axis to measure the absolute Cu-Kα yield. The absolute integral reflectance of the HOPG was measured using an X-ray diffractometer. The HOPG has  ± 16% of spatial nonuniformity of the integral reflectance. The spectral resolution of the spectrometer was 17.9 eV (FWHM).

Figure 2 shows an example of Cu-Kα spectra. Cu-Kα X-ray yield produced during the compression process (green dotted line) was negligibly weak compared with those produced with the heating lasers (red solid and black dashed lines). The green dotted and red solid lines are spectra obtained with application of the external magnetic field. The vertical error bar corresponds to spatial nonuniformity of the integral reflectance, and the horizontal error bar is equal to the spectral resolution. The Cu-Kα photon yields were integrated within the energy range of 8.0–8.1 keV after background subtraction.

Fig. 2
figure 2

An example of Cu-Kα spectra peaked at 8.05 keV emitted from the compressed plasma core. Red solid, black dashed, and green dotted lines are, respectively, spectra obtained by heating with application of an external magnetic field, heating without the application, and only fuel compression with application of an external magnetic field. The vertical error bar corresponds to 16% of spatial nonuniformity of the integral reflectance of the HOPG. The horizontal error bar of 17.9 eV is the spectral resolution of the spectrometer. The Cu-Kα photon yields were integrated within the energy range of 8.0–8.1 keV

Initial magnetic field profile calculation

An externally applied magnetic field penetrates diffusively into the target from the outside. The diffusion timescale is determined by the electrical conductivity and geometry of the targets. Based on a previous study, the externally applied magnetic field is guaranteed to penetrate rapidly into an insulator hydrocarbon29; however, the diffusion dynamics of the magnetic field into the gold cone remained unclear.

The magnetic diffusion time (tdif) is expressed as tdif = μ0σaL/2 for a cylindrical tube that is topologically similar to a cone, where μ0, σ, a, and L are the permeability, electrical conductivity, radius, and diffusion layer thickness of the tube, respectively. The temporal change in magnetic field strength drives an inductive current in the gold cone, and the current ohmically heats the gold. The electron conductivity of the gold depends on its temperature. For a 7-μm-thick gold cone wall at 50-μm radius cone tip, the diffusion times are tdif = 18 ns, 860 ps, and 420 ps at room temperature (σ = 4 × 107 S m−1), 0.1 eV (σ = 2 × 106 S m−1), and 1 eV (σ = 1 × 106 S m−1), respectively32. The small temperature increment helps to rapidly soak the cone in the magnetic field.

For more detailed estimation, we developed an electro-magneto dynamics simulation code with consideration of the inductive heating and the temperature dependence of conductivity to calculate the spatial and temporal distribution of a magnetic field in the cone-in-ball target33. A current of 250 kA was generated with a capacitor–coil target driven by one GEKKO-XII beam as measured using proton radiography24. The current of 250 kA was used in the following analysis, although three GEKKO-XII beams were used in this integrated experiment.

Figure 3 shows the two-dimensional profile of the magnetic field calculated at the maximum magnetic field strength timing generated using a 250-kA Gaussian current pulse with a 1-ns (FWHM) pulse width. The diffused magnetic field profiles were calculated using a constant electrical conductivity σ = 4 × 107 of a cold gold (Fig. 3a) or temperature-dependence conductivity of gold (Fig. 3b). The magnetic field strengths at the tip are Bint = 96 and 335 T, respectively. Temporal changes of the magnetic field profile are discussed in Supplementary Methods and are shown in Supplementary Fig. 2.

Fig. 3
figure 3

Two-dimensional profiles of the magnetic field calculated with different electrical conductivities. The magnetic field diffusion at the peak timing of current flow in the coil was calculated using (a) the constant electrical conductivity [σ = 4 × 107 S m−1] and (b) with consideration of temporal changes in temperature and conductivity of the gold cone due to the inductive heating

In this study, the magnetic field lines are bent due to magnetic field compression associated with the plasma compression. If a magnetic mirror is formed in the transport region, reflection of the REB by the magnetic mirror could degrade the laser-to-core coupling. The magnetic field strength at the ball center is B0 = 225 T before the compression, as shown in Fig. 3b. The ball radius at the maximum compression timing was about rmin = 50 μm, as shown in Fig. 4b. The magnetic field strength at the maximum compression timing (Bmax) can be estimated with Bmax = \(B_0\left( {r_0{\mathrm{/}}r_{{\mathrm{min}}}} \right)^{2\left( {1 - 1/{\mathrm{Re}}_{\mathrm{m}}} \right)}\), where r0 = 125 μm and Rem ~ 231,34 are initial radius of the ball and magnetic Reynolds number, respectively, which gives Bmax = 560 T. Therefore, the mirror ratio along the REB path, the ratio of magnetic field strengths at the peak and the generation zone, is Rm = Bmax/Bint = 1.7.

Fig. 4
figure 4

The core density profiles of the pre-compressed core at three different timings. a Experimental layout of the pre-compressed core density measurement experiment. The LFEX laser was used to generate a Cu-Kα backlight flash. b Density profiles measured at t = 0.38, 0.72, and 0.92 ns after the peak of the compression beam pulse. The numbers on the contour lines represent mass density (g cm−3)

The REB transport was simulated in the mirror geometry with several mirror ratios (Rm) from 0 to 2035, and it was determined that a moderate mirror (Rm < 10) can focus on the REB without a significant loss caused by the mirror effect. This mirror ratio Rm = 1.7 is small enough to guide the REB efficiently to the core in this system with avoiding the mirror effect.

Two-dimensional density profile measurement

Flash X-ray backlighting with a monochromatic imager17,36,37,38 was used to measure two-dimensional density profiles of pre-compressed Cu(II) oleate solid balls under an external magnetic field.

The experimental layout is shown in Fig. 4a. The X-ray shadows of compressed cores were recorded using imaging plates with the same spherically bent quartz crystal used in the laser-to-core coupling experiment. The solid ball specifications and the laser parameters of the compression and magnetic-field-generation beams were also identical to those used in the laser-to-core coupling measurement. The LFEX laser was used for flash Cu-Kα X-ray backlight generation in this experiment. In order to generate a large-format backlight, the LFEX laser was defocused to produce a 350-μm-diameter spot on a 20-μm-thick Cu foil at 3 mm behind the solid ball along the line of sight of the crystal imager.

An X-ray shadow is converted to an X-ray transmittance profile by interpolating the two-dimensional backlight X-ray intensity profile within the core region from the outside of the core region. The areal density of the pre-compressed core was calculated from the X-ray transmittance profile with a calculated opacity of 100 eV Cu(II) oleate for 8.05-keV X-rays39. A two-dimensional density profile of the core was obtained after applying an inverse Abel transformation to the areal density profile, assuming rotational symmetry of the core along the cone axis.

Figure 4b shows the core density profiles at t = +0.38, +0.72, and +0.92 ns. The numbers on the contour lines represent mass density (g cm−3). The converging shock wave was still traveling to the center of the ball at t = +0.38 ns, maximum compression was reached at around t = +0.72 ns, and the core had already begun to disassemble at t = +0.92 ns. The areal mass densities (ρL) and average mass density (ρ) of the core along the REB path length (L) were, respectively, ρL = 0.08 g cm−2 and ρ = 5.7 g cm−3 at t = +0.38 ns, ρL = 0.16 g cm−2, and ρ = 11.3 g cm−3 at t = +0.72 ns. These values were used in the calculation of the correlation factor described in the Methods section.

Discussion

Figure 5 shows the dependence of the measured coupling efficiency on the heating laser intensity (bottom axis) and energy (top axis). Solid and open marks represent the couplings of two injection-timing groups t = 0.61–0.72 ns and t = 0.37–0.42 ns. The solid and dashed lines are fitted, as an eye guide, to the couplings measured with or without application of the external magnetic field neglecting the injection timing difference.

Fig. 5
figure 5

Dependence of laser-to-core energy coupling on heating laser intensity and energy. The blue rectangular, green triangle, and red circle marks represent laser-to-core coupling efficiencies obtained with the following conditions: no application of an external magnetic field with the open-tip cone, application of an external magnetic field with the closed-tip cone, and application of an external magnetic field with the open-tip cone, respectively. The error bars show the uncertainty of the experimental measurement described in the text. Solid and open marks represent the couplings of two injection-timing groups t = 0.61–0.72 ns and t = 0.37–0.42 ns. The error bars are the spectrometer uncertainty due to nonuniformity of the the integral reflectance. The solid and dashed lines are fitted, as an eye guide, to the couplings with neglecting the injection timing difference

It is not possible to identify the exact area of the HOPG that diffracts the detected Cu-Kα X-rays. Therefore, we considered the ±16% of nonuniformity of the integral reflectance over its full surface to impact the uncertainties in evaluations of the deposited energy and the laser-to-core coupling. Nonetheless, it is important to refer here that the HOPG alignment in respect to the targets was not changed during these measurements, and the detected Kα photons were always reflected from the same area in the crystal. This fact means that the true coupling points do not scatter randomly within the error bars but are at certain points keeping relations between other points within the error bars.

The figure reveals clearly that the REB guiding by the external magnetic field (red circular marks) shows 1.8 times higher coupling than that without an external magnetic field (blue rectangular marks), and that electron scattering in the closed cone-tip (green triangular marks) reduces the coupling 0.8 times compared to the open-tip values (red circular marks). Furthermore, the coupling was degraded gradually by increasing the heating laser energy with keeping both pulse duration and spot diameter unchanged (1.8 ps and 50 μm), because a higher-intensity laser produces higher-temperature REB, as shown in Table 2 and this results in less guiding and less energy deposition. We anticipate that the pulse duration should be extended to sustain this efficient coupling for higher laser energy.

Table 2 Summary of correlation factors used in the analysis

As shown in Fig. 2, Cu-Kα X-ray yield produced during the compression process (green dotted line) was negligibly weak compared with those produced during the heating process (red solid and black dashed lines) even with application of the external magnetic field. This result means that Cu-Kα enhancement is not due to an improved confinement of thermal electrons in the compressed plasma but to the guiding of the REB.

Figure 6 shows two-dimensional Cu-Kα emission profiles and density profiles of the pre-compressed core at two different timings (t = 0.40 ± 0.03 ns and 0.69 ± 0.04 ns), and also the comparison of Cu-Kα emission profiles with and without the external magnetic field application at the two different timings. In the Cu-Kα emission profiles, the numbers on the contour lines represent Cu-Kα emissivity normalized with the heating laser energy [×109 photons sr−1 cm−3 J−1], which was calculated from the combination of the absolute Cu-Kα spectra and the two-dimensional Cu-Kα monochromatic images after applying an inverse Abel transformation to the line-integrated emission profiles assuming rotational symmetry of the emissivity along the cone axis.

Fig. 6
figure 6

Two-dimensional emission and density profiles. Two-dimensional profiles of Cu-Kα emission (a, c, e, g) and mass density (b, d). These profiles were measured in the experiments with the application of the external magnetic field at t = 0.40 ± 0.03 ns (a, b, e) and t = 0.69 ± 0.04 ns (c, d, g). The numbers on the contour lines represent heating-laser energy-normalized Cu-Kα emissivity (×109 photons sr−1 cm−3 J−1) and mass density (g cm−3), respectively. Cu-Kα emission profiles are compared between those obtained with (e, g) and without (f, h) application of the external magnetic field at t = 0.40 ± 0.03 ns and t = 0.69 ± 0.04 ns. These images were obtained after applying an inverse Abel transformation to the line-integrated emission profile, assuming rotational symmetry of the core along the cone axis

At t = 0.40 ± 0.03 ns, the converging shock wave was still traveling to the center of the ball, and the shock front is clearly observable in Fig. 6b. Strong Cu-Kα emission region locates behind the shock front in Fig. 6a and e. Externally applied magnetic field lines were accumulated in the downstream of the shock wave; therefore, REB was guided to the shock-compressed region along the field lines. This feature was less prominent without the external magnetic field as shown in Fig. 6f.

At t = 0.69 ± 0.04 ns, the solid ball reached the maximum compression (Fig. 6d). The strong Cu-Kα emission spot appeared at 50-μm longitudinal distance, which was more than 100 μm away from the REB generation zone, namely the cone tip as shown in Fig. 6c and g. This Cu-Kα emission feature is the evidence of the long-distance REB guiding by the externally applied magnetic field as well as significant enhancement of the laser-to-core coupling. This strong emission spot disappeared, when the external magnetic field was not applied, as shown in Fig. 6h. Energy shift of the Cu-Kα X-ray due to ionization of Cu atoms in a hot core could be the reason why the Cu-Kα emission was weak in the core central region40. Note that the strong emission spot that appeared in Fig. 6c and g was found also in the other shot (ID 40556) performed with similar injection timing (t = 0.61 ± 0.02 ns).

The laser-to-core coupling (η) can be simplified as a product of laser-to-REB energy conversion efficiency (ηREB), REB collision probability (ηcoll), and energy deposition rate of REB in the core (ηdep)41.

$$\eta = \eta _{{\mathrm{REB}}} \cdot \eta _{{\mathrm{col}}} \cdot \eta _{{\mathrm{dep}}}$$
(1)

The previous experiments show that ηREB = 0.441 and ηcoll = 0.729. ηdep was calculated by the simplified model41 with the measured areal density and REB temperatures shown in Table 2.

$$\eta _{{\mathrm{dep}}} = {\textstyle{{AT_{{\mathrm{REB1}}}^2} \over {AT_{{\mathrm{REB1}}}^2 + (1 - A)T_{{\mathrm{REB2}}}^2}}} \cdot {\textstyle{{\rho L} \over {0.6T_{{\mathrm{REB1}}}}}} + {\textstyle{{(1 - A)T_{{\mathrm{REB2}}}^2} \over {AT_{{\mathrm{REB1}}}^2 + (1 - A)T_{{\mathrm{REB2}}}^2}}} \cdot {\textstyle{{\rho L} \over {0.6T_{{\mathrm{REB2}}}}}}.$$
(2)

As an approximated relation of Eq. (11) in ref.42, RREB[g cm−2] = 0.6 fR TREB [MeV] is used to calculate the RREB from the experimentally measurable parameter (TREB), here, fR is an adjustable parameter, set to 1 in the standard model. Finally, the simple model yields η = ηREB · ηcoll · ηdep = 6.2%. This simple evaluation seems fairly consistent with the measured coupling (7.7 ± 1.2%), and the simple evaluation reveals that higher areal density core leads to higher laser-to-core coupling. An ultra-high-energy-density state could be efficiently created by the magnetized fast isochoric heating. Integrated REB transport simulation result is described in the Supplementary Discussion and Supplementary Fig. 1. The simulation result is fairly consistent with the experimental results.

We have demonstrated experimentally the magnetized fast isochoric heating scheme for creating an ultra-high-energy-density state as a potential path to the ignition spark formation. The enhancement of the laser-to-core coupling as well as the strong Cu-Kα emission spot located at 100 μm away from the cone tip are the evident features produced by guiding of the diverged REB with the externally applied magnetic field in the long-transport distance. An energy-density increment of the heated core is close to 1 Gbar, which corresponds to 50 J of the energy deposition in a 100-μm-diameter spherical volume. Plasma hydrodynamics, generation and transport of electron/ion beams, thermal conduction, and α-particle transport will be able to be controlled by the externally applied strong magnetic field. There is no doubt that laser–plasma experiments with strong magnetic fields contain a lot of unexplored physics; therefore, this research also stimulates spin-off sciences in the field of atomic physics, nuclear physics, and astrophysics which act to broaden inertial fusion sciences and high-energy-density sciences.

Methods

Derivation of the correlation factor between Cu-K α photons and deposited REB energy

The Solodov model43 was used to calculate the stopping power S(E,ρ) [J s−1] of the REB in a core, where E is electron kinetic energy. Both Davies44 and Hombourger45 models were used to calculate the electron-impact K-shell ionization cross-section \(\sigma _{K_\alpha }(E)\) [m−2]. The ratio between the stopping power and K-shell ionization cross-section is the correlation factor (C [J photons−1 sr]). Note that the Hombourger model, which was used in the previous experiment18, gives a 1.3 times higher correlation factor than the Davies model because \(\sigma _{K_\alpha }(E)_{}^{}\) calculated with the Davies model is larger than that with the Hombourger one. The correlation factor C is defined as follows:

$$C = \frac{{{\int}_0^L {\kern 1pt} {\int}_{10{\kern 1pt} {\mathrm{keV}}}^{1{\kern 1pt} {\mathrm{GeV}}} {\kern 1pt} \varepsilon _{{\mathrm{dep}}}(E,x)\mathrm{d}E\mathrm{d}x}}{{{\int}_0^L {\kern 1pt} {\int}_{10{\kern 1pt} {\mathrm{keV}}}^{1{\kern 1pt} {\mathrm{GeV}}} {\kern 1pt} P_{K_\alpha }(E,x)\mathrm{d}E\mathrm{d}x}}$$
(3)
$$= \frac{{{\int}_0^L {\kern 1pt} {\int}_{10{\kern 1pt} {\mathrm{keV}}}^{1{\kern 1pt} {\mathrm{GeV}}} {\kern 1pt} v(E,x)f(E,x)S(E,\rho )\mathrm{d}E\mathrm{d}x}}{{{\int}_0^L {\kern 1pt} {\int}_{10{\kern 1pt} {\mathrm{keV}}}^{1{\kern 1pt} {\mathrm{GeV}}} {\kern 1pt} v(E,x)f(E,x)n_{{\mathrm{Cu}}}\sigma _{K_\alpha }(E){\mathrm{/}}4\pi \mathrm{d}E\mathrm{d}x}}$$
(4)

where εdep(E, x), PKα(E, x), v(E, x), f(E, x), and nCu are the collisionally deposited energy by the REB to the dense core [J s−1], the probability of Cu-Kα emission [photons s−1 sr−1], the relativistic electron velocity [m s−1], the REB energy distribution, and the number density of Cu atoms in a core [m−3], respectively. The initial REB energy distribution was a Boltzmann function as f(E, 0) = exp(−E/TREB) at the generation zone (x = 0), where TREB is the slope temperature of the energy distribution. The distribution at x (f(E, x)) is calculated by considering the collisional slowing down during the transport.

Figure 7 shows the calculated correlation factor at different REB slope temperatures (TREB) for a ρ = 11.3 g cm−3 and ρL = 0.16 g cm−2 Cu(II) oleate plasma, which are equal to those observed at t = +0.72 ns. The correlation factor also depends weakly on the core density. Dependence of the correlation factor on TREB was considered in the coupling evaluation. The average values of the two correlation factors calculated with the two models were used in the data analysis.

Fig. 7
figure 7

Correlation factor. Dependence of the correlation factor between the deposited energy (J) and Cu-Kα yield (photons/sr) on the REB slope temperature. The data were calculated using the Davies (blue circular marks) and Hombourger (red triangular marks) models of electron-impact K-shell ionization for a 11.3 g cm−3 and 0.16 g cm−2 Cu(II) oleate, which correspond to the average core density and areal density at t = 0.72 ns (maximum compression timing). The average values of the two correlation factors calculated were used in the analysis

Energy distributions of electrons escaped from the plasma into vacuum were measured using an electron energy analyzer positioned along the LFEX incident axis. Although the energy distribution of the so-called vacuum electrons is not exactly identical to that at the generation zone due to scattering, absorption, and reflection by the cone, core, and spontaneous electromagnetic fields, we found that the slope temperatures of the escaped electrons are close to those in the transport region estimated from Bremsstrahlung X-ray spectra41. Therefore, the energy distribution of the measured vacuum electrons was used in the correlation factor calculation. The energy distribution of the escaped electrons was fitted with a two-temperature Boltzmann distribution function as f(E) = A exp(−E/TREB1) + (1 − A) exp(−E/TREB2), where A and E are the intercept and electron energy, and TREB1 < TREB2. The correlation factor calculated for each shot is summarized in Table 2 along with other measured parameters. The error in the deposited energy was evaluated taking into account 16% of spatial nonuniformity of the integral reflectance of the HOPG.

Code availability

The computer codes used in the current study are accessible from the corresponding author upon reasonable request.