1. INTRODUCTION
Research on ion beam heating of matter is a feasible method to investigate and generate warm dense matter (WDM) and high energy density (HED) matter, which can be applied in basic physics (Brambrink et al., Reference Brambrink, Roth, Blazevic and Schlegel2006; Patel et al., Reference Patel, Mackinnon, Key, Cowan, Foord, Allen, Price and Ruhl2003; Tahir et al., Reference Tahir, Schmidt, Shutov, Lomonosov, Piriz, Hoffmann, Deutsch and Fortov2009), astrophysical processes (Nettelmann et al., Reference Nettelmann, Holst, Kietzmann, French, Redmer and Blaschke2008), heavy ion beam research (Hoffmann et al., Reference Hoffmann, Tahir, Udreal, Rosmej, Meister, Varentsov, Roth, Schaumann, Frank, Blažević, Ling, Hug, Menzel, Hessling, Harres, Günther, El-Moussatil, Schumacher and Imran2010; Stokl et al., Reference Stöckl, Frankenheim, Roth, Suß, Wetzler, Seelig, Kulish, Dornik, Laux, Spiller, Stetter, Stöwe, Jacoby and Hoffmann1996; Tahir et al., Reference Tahir, Kim, Matvechev, Ostrik, Lomonosov, Piriz, Cela and Hoffmann2007), and high energy density physics (HEDP) (Hoffmann et al., Reference Hoffmann, Fortov, Kuster, Mintsev, Sharkov, Tahir, Udrea, Varentsov and Weyrich2009; Nellis, Reference Nellis2006; Tahir et al., Reference Tahir, Deutsch, Fortov, Gryaznov, Hoffmann, Kulish, Lomonosov, Mintsev, Ni, Nikolaev, Piriz, Shilkin, Spiller, Shutov, Temporal, Ternovoi, Udrea and Varentsov2005). In particular, the beam-matter heating techniques are of considerable potential interest for driving an inertial confinement fusion (ICF) (Deutsch, Reference Deutsch1986; Reference Deutsch1992; Deutsch et al., Reference Deutsch, Maynard, Bimbot, Gardes, Dellanegra, Dumail, Kubica, Richard, River, Sernagean, Fleurier, Sanba, Hoffmann, Weyrich and Wahl1989; Hora, Reference Hora2007; Hoffmann, Reference Hoffmann2008), where a target matter can be irradiated by intense ion beams and be transformed into a plasma of high density and high temperature, which can initiate fusion reactions (Deutsch, Reference Deutsch1992; Tahir et al., Reference Tahir, Deutsch, Fortov, Gryaznov, Hoffmann, Kulish, Lomonosov, Mintsev, Ni, Nikolaev, Piriz, Shilkin, Spiller, Shutov, Temporal, Ternovoi, Udrea and Varentsov2005). Much work has been made to study the beam-matter heating including heavy ion (Hoffmann et al., Reference Hoffmann, Tahir, Udreal, Rosmej, Meister, Varentsov, Roth, Schaumann, Frank, Blažević, Ling, Hug, Menzel, Hessling, Harres, Günther, El-Moussatil, Schumacher and Imran2010; Zhao et al., Reference Zhao, X, Xu, Zhao, Xia, Jin, Ma, Liu, Yang, Zhang, Wang, Li, Zhao, Zhan, Xu, Zhao, Li and Chen2009), proton (Patel et al., Reference Patel, Mackinnon, Key, Cowan, Foord, Allen, Price and Ruhl2003; Tahir et al., Reference Tahir, Schmidt, Shutov, Lomonosov, Piriz, Hoffmann, Deutsch and Fortov2009), and laser beams (Zhao et al., Reference Zhao and Shin2012). Intense ion beams have emerged as efficient and flexible tool to heat a target to HED matter, where beam ions propagate in a target and volumetrically heat it over macroscopic volumes.
Two general approaches to beam-matter heating have been identified. First, at GSI Darmstadt, long (~3 mm) targets have been heated by very high energetic beams of several GeV Uranium (Tahir et al., Reference Tahir, Schmidt, Shutov, Lomonosov, Piriz, Hoffmann, Deutsch and Fortov2009; Reference Tahir, Shutov, Piriz, Lomonosov, Deutsch, Spiller and Stöhlker2011), where the target temperatures can reach an order of eV. A cylindrical target has been facially irradiated by the large hadron collider (LHC) proton beam along the axis, where the sample was heated to HED matter (Tahir et al., Reference Tahir, Schmidt, Shutov, Lomonosov, Piriz, Hoffmann, Deutsch and Fortov2009). Second, a number of researchers (Barnard et al., Reference Barnard, Armijo, More, Friedman, Kaganovich, Logan, Marinak, Penn, Sefkow, Santhanam, Stoltz, Veitzer and Wurtele2007) at the Heavy Ion Fusion Virtual National Laboratory (HIF VNL) have studied the feasibility of using beam ions at lower energy (~1 MeV/u) and shorter pulse time (≤1 ns) with a focal spot (~1 mm), where the ion range (~1–100 µm) is much smaller than the radius of the focal spot. The U.S. HIF VNL is planning to heat thin foils using low- to medium-mass ions of much lower energy to heat a 1 mm radius volume of a target to electron temperatures of a few eV in 1 ns. A 48 μm thick 10% solid density aluminum foam, has been heated in simulations by a 20 MeV, 330 A, 1 ns, ${\rm N}_e^+$ ion beam with a 1 mm radius, to reach 5 eV central temperature at the U.S. HIF VNL.
In this context, the near term research point is WDM, since the laser and accelerator technology has improved sufficiently to provide intense ion beams to heat samples of target and reach the regime of WDM. Studies of beam-matter heating, which play an essential role in generation of WDM, have been encouraged to be explored thoroughly. Due to this reason, the need to comprehend the complexity of target heating driven by an intense ion beam, has prompted a number of theoretical and experimental studies and numerical simulations in this field. Despite these studies, the underlying physics and effects of WDM have not fully been understood yet.
In particular, quantum hydrodynamic (QHD) model (Chen et al., Reference Chen, Cockburn, Gardner and Jerome1995; Haas et al., Reference Haas, Manfredi and Feix2000; Manfredi et al., Reference Manfredi and Haas2001) has been developed by solving the nonlinear Schördinger-Poisson or the Wigner-Poisson kinetic models, which can be devoted to studying the properties of WDM. The QHD method has its advantages that the QHD equations are quite simple for numerical studies and have a straightforward interpretation in terms of fluid quantities that are employed in classical physics. Moreover, this model can obtain time-evolution of a target temperature and density quantitatively, and can explicitly simulate an ion beam propagating process in the target. In our previous works, we have implemented the QHD theory to simulate a test point charge ((Zhang et al., Reference Zhang, Song and Wang2011b) and a test ion-cluster beam. To our best knowledge, a two-dimensional (2D) QHD model have not been performed to study the interaction process of continuous beam and a solid target carefully.
In this paper, the 2D hydrodynamic approach is employed to study the dynamics of two-componant (ions and electrons) aluminum solid target for an isochoric heating of the target under a proton beam interaction. Here 2D QHD model is used for electron dynamics, and 2D classical hydrodynamic formalism is considered for target ion and beam ion dynamics. The proton beam can heat the solid target to warm dense states on a picosecond time scale in which hydrodynamic expansion of the target is negligible. The electron energy-balance equation is supplied to calculate the target electron temperature. We shall explicitly express the time-evolution of the target temperature and density in order to show how the WDM is generated from a cold target. The nonlinear quantum and classical hydrodynamic equations are solved by employing the flux corrected transport (FCT) method (Boris et al., Reference Boris, Landsberg, Oran and Gardner1993), which has proven to be an accurate and easy to use algorithm to solve nonlinear time-dependent conservative continuity equations of hydrodynamics. The paper is organized as follows. In Section 2, the QHD formalisms for target electrons and the classical hydrodynamic model for target ions and beam ions, coupled with the Poisson's equation, are presented to describe the dynamics of beam ions and the two-component target. In Section 3, numerical results are shown and discussed. Finally, a short summary is given in Section 4. Gauss units will be adopted throughout the paper except in specific definitions.
2. QUANTUM HYDRODYNAMIC MODEL AND SIMULATION METHOD
A 2D hydrodynamic model is performed to investigate the time-dependent dynamic polarization and heating of two-component (ion and electron) aluminum solid target when irradiated by a proton beam pulse. In a cylindrical coordinate system with {r,z,φ}, φ is neglected in this work, since the system is symmetric with respect to z axis (r = 0). The 2D QHD model is employed to describe the target electron dynamics, with electron velocity field ue (r,z,t), electron number density n e (r,z,t), and electron energy density W e (r,z,t). The dynamics of target ions and beam ions are described by a 2D classical hydrodynamic model, with the target ion velocity field ui (r,z,t) and number density n i (r,z,t), the beam ion number density n b (r,z,t), and velocity field ub (r,z,t). The hydrodynamic methods involve calculation of the conservation equations of particle number, momentum and energy, where a self-consistent determination of these conservative quantities is required with the help of Poisson's equation. A proton beam pulse with a particle energy of E b (1.5–10 MeV) and an intensity of N (7.9 × 108) ions that is the total particle number delivered in a single pulse time 10 ps, propagates along the z direction and irradiates perpendicularly on the target. The transverse Gaussian distribution is considered with a focal radius r b (23 μm), which is defined as a multiple of the root mean square radius. The beam parameters here are selected to be relevant with the experiment carried by Patel et al. (Reference Patel, Mackinnon, Key, Cowan, Foord, Allen, Price and Ruhl2003).
In what follows, the simulation starts at t = 0 at a room temperature. At t = 0, the target ions and electrons are at equilibrium states and the velocities of them are all zero (ue (r,z,t = 0) = 0,ui (r,z,t = 0) = 0). The beam is projected from the axis z = 0 at t = 0 with a constant velocity and a initial density n b0. Our simulation results are obtained for a semi-infinite long cylindrical solid aluminum target with initially equilibrium electron density n e (r,z,t = 0) = n0 = 1.18 × 1023 and ion density n i satisfying Z in i (r,z,t = 0) = n0. Z i = 3 is the effective charge state of aluminum target ion, which is equal to the number of valence electrons provided that the bound electrons are not ionized. At the axis z = 0 fixed boundary condition is applied for protons and open boundary condition is used for target ions and electrons. Equilibrium boundary conditions are implemented for z → ∞ or r → ∞ where there is no beam propagation and the target ions and electrons are not disturbed. At the axis r = 0 symmetric boundary condition is presented.
The target electron dynamics are described as:


and the energy-balance equation (Chen et al., Reference Chen, Cockburn, Gardner and Jerome1995)

with $W_e={\textstyle{3 \over 2}}n_e T_e+{\textstyle{1 \over 2}}m_e n_e u_e^2 - {\textstyle{{\hbar ^2 n_e } \over {24m_e }}}\nabla ^2\, log\lpar n_e \rpar +O\lpar \hbar ^4 \rpar $. Thus, the electron temperature is
$T_e={\textstyle{2 \over {3n_e }}}\left[{W_e - {\textstyle{1 \over 2}}m_e n_e u_e^2+ {\textstyle{{\hbar ^2 n_e } \over {24m_e }}}\nabla ^2\, log\lpar n_e \rpar - O\lpar \hbar ^4 \rpar } \right]$. In particular, the terms proportional to
$\hbar ^2 $ in Eqs. (2) and (3) are regarded as the additional quantum terms. In Eq. (2), the quantum term
${\textstyle{{\hbar ^2 } \over {2m_e^2 }}}\nabla \lpar {\textstyle{1 \over {\sqrt {n_e } }}}\nabla ^2 \sqrt {n_e } \rpar $ is induced by a quantum pressure
$P^Q= {\textstyle{{\hbar ^2 } \over {2m_e }}}\lpar \lpar \nabla \sqrt {n_e } \rpar ^2 - \sqrt {n_e } \nabla ^2 \sqrt {n_e } \rpar $, which is related to the derivative of the root density
$\sqrt {n_e } $. The target ion dynamics are described as:


The non-relativistic beam fluid equations are:


The above equations are closed by adding Poisson's equation

where $\nabla={\textstyle{\partial \over {\partial r}}}{\bf e}_r+{\textstyle{\partial \over {\partial z}}}{\bf e}_z $. Z b = 1 is the charge state of beam ion as a proton considered here. Φ is the total electric potential generated by the breaking of local neutrality inside the target by the irradiated beam ions and the collective response of the target ions and electrons to the propagating beam ions, which is usually generated in and around the beam ion position.
Here, m e is the electron mass, m i is the target ion mass, m b is the beam proton mass, e is the elementary charge, qe, T e, P e and P i are the electron heat flux, temperature, pressure and ion pressure, respectively. γe, γi, γb and, τw obtained from our previous work (Zhang et al., Reference Zhang, Song and Wang2011a), are the frictional frequencies with τw = 2γe, where the subscription of e, i, and b express the electron, ion and beam ion, respectively. The quotidian equation of state (QEOS) model (More et al., Reference More, Warren, Young and Zimmerman1988) is supplied to calculate the electron pressure P e by finite temperature Thomas-Fermi theory and the ion pressure by modified Cowan model. In particular, the heat flux is specified by Fourier's law ${\bf q}_e=- {\rm \kappa} \nabla T_e $, where the thermal conductivity κ is defined by (Chen et al., Reference Chen, Cockburn, Gardner and Jerome1995). In the right-hand side, the second term of Eq. (2) is semiclassical statistics pressure term, and the last terms of Eqs. (2), (3), (5), and (7) are the frictional forces.
The above Eqs. (1)–(8) involve a self-consistent determination of the target ion and electron density, electron temperature, beam ion density, and the total electric field ${\bf E}=- \nabla \Phi $ with the components E r and E z. The Poisson's equation is solved by the successive over relaxation (SOR) method at every time step. The 2D FCT method (Boris et al., Reference Boris, Landsberg, Oran and Gardner1993) is adopted to numerically solve the Eqs. (1)–(7) by split-time integration from the initial time t = 0 when the values of all quantities are known, whose general idea is as follows. For solving Eqs. (1)–(7) a minimum amount of numerical diffusion must be added to assure stability, since these equations are subject to a numerical instability. An antidiffusion term that can ensure stability, which should not generate new maxima or minima in the solution, nor accentuate already existing extrema, is considered to remove this additional diffusion and correct the flux in FCT.
For convenience we introduce the dimensionless variables u/u 0 → u, l/l 0 → l, t/t 0 → t, where u, l, and t stand for any quantity of velocity, length and time, respectively. The units u 0 = 109, l 0 = 10−3 and t 0 = l 0/u 0 are chosen to bestly meet numerical stability.
3. NUMERICAL RESULTS
Here, we will take the beam particle energy E b = 1.5 MeV, intensity N = 7.9 × 108, the focal radius r b = 23 µm, and a fixed pulse time 10 ps in all figures except in specific definition. Note that during the pulse time 10 ps, the beam is consecutive.
Figure 1 shows the distribution of the electron density and temperature along the axis z (r = 0) with and without quantum effect at three different times 0.05, 0.5, and 1 ps. Here the beam parameters are E b = 1.5 MeV, n b0 = n 0/103. It is interesting to see that at 0.05 ps and 1 ps, there is no quantum effect to influence the target electron density and temperature, while at 0.5 ps the quantum effect is evident. The essence of quantum effect is to make the electrons back to equilibrium states and when the electrons are at a quasi-equilibrium state the quantum effect is very weak, therefore, as seen from Figure 1b, at 0.05 ps the quantum effect is so small (attributed to its physical nature) that can not affect the electron density and temperature. At 1 ps, the numerical characteristics of the quantum terms ${\textstyle{{\hbar }^2 \over {2m_e^2 }}}\nabla \lpar {\textstyle{1 \over {\sqrt {n_e } }}}\nabla ^2 \sqrt {n_e } \rpar $ and
$ - \textstyle{{{\hbar }^2 n_e } \over {24m_e }}\nabla ^2\, log\lpar n_e \rpar $ lead the quantum effect to be failing to influence the interacting process. Overall, the quantum effect is obvious only at the very beginning the beam-plasma interaction, however, the afterward plasma evolution will definitely be affected by this process. Furthermore, we also found that including quantum effect can stabilize the numerical calculation here. Therefore, the quantum effect is considered in all following results and figures.

Fig. 1. (Color online) Electron density (a) and temperature (b) versus z at r = 0, without quantum effect (solid line) and with quantum effect (dashed line) at three different times 0.05 ps, 0.5 ps and 1 ps. Initial beam proton energy is E b = 1.5 MeV and density is nb0 = n 0/103.
Figure 2 represent the z directional (at r = 0) evolution of beam proton density Figure 2a, electron density Figure 2b, and temperature Figure 2c for three different times 10 ps (solid line), 6 ps (dashed line) and 4 ps (dotted line). E b = 1.5 MeV. At t = 10 ps all beam ions are into the target and when t < 10 ps the only part of the beam ions are into the target. In order to express more clearly, beam ion density and electron density are normalized by initial equilibrium densities. The normalized beam proton densities are similar, while the electron density and temperature are variable at different times with a constant proton energy E b = 1.5 MeV. As shown in Figure 2c that the pulse duration 10 ps is short enough to allow the central (at r = 0) temperature to reach its maximum value of about 5 eV without undergoing hydrodynamic expansion that would lower the peak temperature from 5 eV. Therefore, the target electrons already reach the WDM states. The WDM temperature range is roughly T e = 1.0−20 eV for solid aluminum. One can observe that the electron density depends sensitively on the distribution of beam ions, and the electron temperature depends closely on the electron density. The electron density at the proton distributed position forms a peak, therefore, the electrons at the left and right sides of the beam proton positions move to the protons and two valleys are formed. The fluctuation of electron density increases and electron temperature increases too, as time elapses. It is interesting to find that beam proton for the first time to reach its range (z = 20 μm) at 8.2 ps, however, more protons are stopped very near the boundary z = 0 attributed to the very high electron density. Correspondingly, there is a strong interaction process between high density protons and electrons in several μm region near the boundary z = 0, which produces a maxima in electron temperature.

Fig. 2. (Color online) Beam proton density (a), electron density (b) and temperature (c) versus z at r = 0, with 10 ps (solid line), 6 ps (dashed line) and 4 ps (dotted line). Beam proton energy is E b = 1.5 MeV.
The 1D figures described above can clearly express the spatial variation of the density and temperature in quantities, while 2D spatial density and temperature spectrums are displayed more visually in Figure 3. Here the time is fixed to 10 ps and E b = 10 MeV. Beam protons impinging on the target let the electrons move to them, which leads to changes in the electron density as well as in the temperature. The beam proton density has its maximum value at the point (r = 0, z = 0), that is, the center of the focus of the incident ion beam, and electron density and temperature have their maximum values near the point. The temperature takes values of T e = 0.5−4 eV in all heated region where WDM is generated, which is in favor of important experiments on HEDP.

Fig. 3. (Color online) Spacial distribution of beam proton density (a), electron density (b), and temperature (c) at 10 ps. E b = 10 MeV.
To understand the time and spatial evolutions of the beam-target heating, in Figure 4, we illustrate beam proton density Figure 4a, electron density Figure 4b, and temperature Figure 4c at r = 0. They demonstrate a clear time and space variations and a larger volume of target has been disturbed as time increases. A highest n b produces a maximum n e, as well as a maximum T e, since the beam protons in the target pull the electrons moving to them and also the hydrodynamic expansion is neglected in a picosecond time scale. The number of protons entered in the target is gradually increased during the pulse time 10 ps, as a result, the changes of electron density and temperature increase and the changes can reach maximum values at 10 ps.

Fig. 4. (Color online) Space (z) and time (t) resolved beam proton density (a), electron density (b), and temperature (c) at r = 0. E b = 1.5 MeV.
To analyze the average influence of the beam and the collective response of the target, an average electron temperature is calculated as a space average in 50 μm × 50 μm (r,z) plane. The influences of N 5(a), r b 5(b), and E b 5(c) on the temporal evolution of this average T e are shown in Figures 5. In Figure 5a one observes a significant increase in the average T e as N increases. As r b decreases in Figure 5b, the distribution volume of the total beam particles becomes small leading to a rise in n b, so that the electron temperature as well as the average T e increases distinctly. The effect of increasing E b, seen in Figure 5c, introduces a reduction in T e, since the beam-target interaction time becomes shorter for higher E b. Note that the electron temperature of the 50 μm × 50 μm heated target region can reach 0.2 eV for E b = 1.5 MeV, N = 7.9 × 108, and r b = 23 μm.

Fig. 5. (Color online) Electron temperature is averaged by a 2D volume of 50 μm × 50 μm at (r,z) plane, defined as Average T e. Time evolution of average T e with r b = 23 μm, E b = 1.5 MeV (a), E b = 1.5 MeV, N = 7.9 × 108 (b), and N = 7.9 ×108, r b = 23 μm (c).
4. SUMMARY
We study a cold solid target heated by a proton beam pulse of particle energy E b, intensity N and focal radius r b of transverse Gaussian distribution, with a pulse time 10 ps. The dynamics of target electrons are presented by the 2D QHD theory and the target ions and beam ions dynamics are examined by a 2D classical hydrodynamic model. Numerical simulations for this beam-target heating are carried out by employing the FCT method. As the beam ions penetrate into the target, the target electrons are collectively polarized and the target is isochorically heated in a picosecond time scale. It is seen that the target is heated to several eV, as a result, WDM is generated in the target on a picosecond time scale. Therefore, one can access the very interesting WDM states in this work.
ACKNOWLEDGMENTS
This work was supported by the National Basic Research Program of China (Grants No. 2010CB832901), National Natural Science Foundation of China (Grant No.10975028), and the Program for New Century Excellent Talents in University (NCET-08-0073).