arXiv:2011.03805v1 [physics.plasm-ph] 7 Nov 2020 oee,i sdffiutt nraeteitniyo CP thresh- J/cm of fluence few the a intensity typically to the with due increase old optics 11]. traditional to [10, using difficult ef- lasers whistler is and standing it [9] by However, shocks heating by collisionless acceleration ion or ion fective laser-driven [8] for pressure means thought radiation ideal are an amplitude be relativistic to with (CP) lasers polarized CP circularly to . laser (LP) polarized linearly 6] [7]. [5, production amplification harmonic pulse higher as and for such Fruit- possibilities devices 4]. the [3, up rou- various energy open laser is interactions of grating laser-plasma flux ful the Plasma deflect to [2]. example, that used shutter contrast For tinely optical pulse ultrafast the [1]. an enhances as optics works traditional mirror the plasma of beyond magnitude point of manipulation orders breaking many device for intensities photonic high medium at Plasma attractive an purpose. as each wavesemerges electromagnetic for the plasmas of and to lasers situations intense sci- desired control to produce beam is The applications quantum astrophysics. these laser of as insight essence laboratory and such new fusion, laser fields brings ence, research which various conditions, for plasma extreme oe eei nteueo netra antcfield. magnetic external an of use pro- the scheme a in of is difference Compared here key posed the [16]. methods, rela- these that laser with known the is alters distributions drastically It anisotropic by transparency polarization. that tivistic circular laser numerically LP into transform demonstrate light can Lehmann crystal 15] photonic mixing. plasma two- [14, transient optical Spatschek wave- of plasma use and as the device by photonic Michel novel intensity. plate a laser propose the 13] on [12, limit severe although a is interactions, there laser-plasma utilizing generation ∗ [email protected] nti ae,w osdrapam eiet ovr a convert to device plasma a consider we paper, this In investigate to tool unique a are lasers Ultra-intense hr r eea da fpam eie o CP-wave for devices plasma of ideas several are There lsacnetfrgnrtn iclryplrzdelectr polarized circularly generating for concept Plasma mltd.Frvrosrsac upss nes Plight laboratory. CP the intense in purposes, plasmas inte density research generating energy various for advantage For an waveamplitude. has CP two device into exter photonic laser an polarized plasma linearly under a foil short- separate intense plasma to thin plate for a device that st photonic simulations field plasma numerical magnetic unique the a and design density to plasma the by determined are rpgto etrso iclryplrzd(P electro (CP) polarized circularly of features Propagation .INTRODUCTION I. aaoh Sano, Takayoshi 1 nttt fLsrEgneig sk nvriy ut,O Suita, University, Osaka Engineering, Laser of Institute Dtd o ,22;acpe o ulcto nPyia Re Physical in publication for accepted 2020; 6, Nov (Dated: 2 . 1, ∗ uueTatsumi, Yusuke eaiitcamplitude relativistic tal. et 1 aaauHata, Masayasu npamsa hslrwvswe h tegho an critical of The strength large. the sufficiently when strength is waves field magnetic whistler external as propagate waves plasmas Electromagnetic in magne- 18]. and impor- [17, lights plasmas the laser tized between on interactions light the shed of tance experiments laser in generated can am- it relativistic the but exceeding design, plitude. lasers simple intense quite to applied a be has device plasma Our nL ae ycnrliglsrpam neatosun- interactions laser-plasma controlling by laser LP an inta h aefrequency wave the that tion rcini lodsusd ial,orcnlsosare ex- conclusions wave our L Finally, VI. the Sec. of in discussed. summarized possibility also field The is magnetic V. and traction laser Sec. the the in examine between We the directions angle and the detail. of strength in effects elucidated field are magnetic intensity transmittance the laser the on of per- reflectivity as dependences is device and The plasma simulations filter. In the the CP-wave of III. a of evaluation functional Sec. survey for systematic in formed a described IV, are Numeri- Particle- Sec. simulations (2D) field. two-dimensional (PIC) of magnetic and in-Cell method external foil and an trans- setup plasma along waves cal a CP laser between in LP interaction is an the interest for Our after strategy II. mitted is design amplitude Sec. The relativistic in out- with follows. waves explained The CP as of is generation wave. paper the whistler this of the line of characteristics unique [23]. waves effect CP Faraday generate and the [22] by waves CP manage of to transmission devices the plasma magnetized stimulates feature ylto rqec,i.e., , cyclotron lcrnms and mass electron sn uo est o h hslrwv,i a propa- density can it any wave, at whistler plasmas the in for gate density cutoff no is h civmn fqaisai togmgei fields magnetic strong quasi-static of achievement The ecnie osbewyt xrc Pwvsfrom waves CP extract to way possible a consider We ecnie lsapooi eiebsdo the on based device photonic plasma a consider We rvln nteopst ieto.This direction. opposite the in traveling s antcwvsi antzdplasmas magnetized in waves magnetic a antcfil ok sapolarizing a as works field magnetic nal us aes ehv eosrtdby demonstrated have We lasers. pulse r togyrqie ocet high create to required strongly are s s Pwvsee iharelativistic a with even waves CP nse B egh hspoet a eapplied be can property This rength. I NELIGCONCEPT UNDERLYING II. c 1 o h hslrwv sgvnb condi- a by given is wave whistler the for n auioSentoku Yasuhiko and aa5507,Japan 565-0871, saka e steeeetr hre ic there Since charge. elementary the is mgei ae with waves omagnetic iwE) view B c = ω m 0 [19–21] e qast h electron the to equals ω 0 /e 1 where , hsessential This . m e sthe is 2 der a uniform magnetic field. When an external magnetic the field strength as nL = Bext +1, the parameter range field is applied to a plasma target along the direction of of the region (I) is given by laser propagation, electromagnetic waves can travel in e e plasma as either the R wave or the L wave. The R (L) wave is a right-hand (left-hand) CP light with respect to 1 < Bext < ne − 1 . (5) the magnetic field direction. The propagation features of e e these CP waves depend on the electron density ne of the The transmittance of the R and L wave at the vacuum- plasma and the magnetic field strength Bext. Assuming plasma boundary is evaluated by the refractive indexes the laser frequency is ω0, the refractive indexes for the R according to the Fresnel formula [19]. The energy frac- and L waves are given by tion of the transmitted wave relative to the injected en- 2 ergy of each CP wave is derived as T∗ =4N∗/(1 + N∗) , 2 ne NR =1 − , (1) where the subscript ∗ denotes R or L. The reflectivity is 1 − B 2 2 e ext then R∗ = (1 − N∗) /(1 + N∗) . 2 ne NL =1 − e , (2) Figure 1(b) shows the transmittance TR and TL for a 1+ Bext e plasma density with ne = 10. The region (I) corresponds 2 2 that 1 < Bext < 9 for this case, which is highlighted where ne = ωpe/ω0 = ne/nc is thee electron density nor- by light-gray in the figure.e The most ideal situation for malized by the critical density n , and B = ω /ω = c ext ce 0 the R-wavee generation is T = 1 and T = 0. In the B /Be is the ratio of the external magnetic field to the R L ext c region (I), the L wave is perfectly reflected, whereas the critical value B . Here, ω and ω aree the electron c pe ce transmittance for the R wave increases with the magnetic plasma frequency and cyclotron frequency, respectively. field strength. Then the highest transmittance for the We focus on the transmittance of the R and L waves whistler wave will be realized at just below the L-cutoff, at the interface between the vacuum and a plasma tar- that is B ∼ n − 1. This magnetic field strength is get. Since the electromagnetic waves can propagate in ext e the optimized condition by which the plasma foil target plasmas only when N 2 is positive in the linear amplitude will worke as an R-wavee filter most efficiently. The higher limit. Figure 1 indicates the propagation features along transmittance could be realized for much denser targets a magnetic field summarized in the density and magnetic as well, although the optimal strength of the magnetic field diagram (n vs. B ), which is the so-called CMA e ext field becomes larger B ∼ n ≫ 1. (Clemmow-Mullaly-Allis) diagram [24]. Here, the cold ext e e e By contrast, only the L wave is allowed to propagate plasma approximation is adopted, and the ion mass is as- e e sumed to be much higher than the electron mass. The pa- in the region (II) of Fig. 1(a). The parameter range of rameter space is divided into several regions by the cutoff the region (II) can be expressed as conditions for the R and L waves (ωR = ω0 and ωL = ω0) and the cyclotron resonance condition (ωce = ω0), where the cutoff frequency for each wave is written as |ne − 1| < Bext < 1 . (6)

1/2 e e ωR 1 2 = Bext + Bext +4ne , (3) The required strength of the magnetic field is much ω0 2     smaller compared to that of the region (I), which could e e 1/2 ωL 1 2 e be an advantage in terms of practical application. The = −Bext + Bext +4ne . (4) ω0 2     transmittance for the case of ne = 1 is depicted in e e e Fig. 1(c). The region (II) for this case is the range of This diversity in the propagation features is quite differ- e Bext < 1. The transmittance of the L wave increases ent from the normal cutoff condition (ωpe = ω0 or ne = 1) with Bext, whereas the R wave starts to enter the plasma without the external magnetic field. e when B > 1. Then the best performance as an L-wave The forbidden regions for the R and L wave aree dif- e ext filter (T = 0 and T ∼ 1) will be achieved at just be- ferent. An LP light is regarded as a combined wave of a e R L right-hand CP (RCP) and a left-hand CP (LCP) wave. low the resonance condition, that is Bext ∼ 1. Because of the lower density of the region (II), the relativistic effects If one of the CP waves is allowed to enter the plasma e target, but the other is forbidden, the plasma target will may not be ignored when the incident laser amplitude is work as a filter to generate a CP wave from an LP light. relativistic (see more Sec. V). Two regions satisfy such conditions, which are labeled by By controlling such propagation properties, we can de- (I) and (II) in Fig. 1(a). sign a plasma device that converts an LP light into CP When the magnetic field strength is larger than the waves. In this paper, we concentrate our discussion on critical Bc, the whistler-wave mode of the R wave prop- the region (I) where only the whistler wave gets into the agates in plasmas at any density. On the other hand, plasma. In the following, we will search for the opti- there is a cutoff density for the L-wave component. Only mal conditions of the plasma photonic device to generate the R wave is allowed to pass through the plasma target CP lights with the help of multi-dimensional PIC simu- in the region (I). Since the L-cutoff density depends on lations. 3

102 1

ωL R (I) 0.5 L c

n 10 /

e (b) Transmittance n 0 1 1 (II) Density 0.5 L R ωR ωce (a) (c) 0.1 Transmittance 0 2 0.1 1 10 102 0.1 1 10 10 Magnetic Field Strength Bext/Bc Magnetic Field Strength Bext/Bc

FIG. 1. (a) Propagation properties of the electromagnetic waves in plasmas along an external magnetic field. The allowance of the propagation for the R and L waves depends on the plasma density ne/nc and the magnetic field strength Bext/Bc. The region (I) stands for where the R wave propagates but the L wave is forbidden, and the region (II) is the opposite. Both waves are allowed to propagate in the white regions, whereas both are forbidden in the dark-gray region. The boundary curves of the separated regions are the cutoff conditions (ωR = ω0 and ωL = ω0) and the resonance condition (ωce = ω0). (b,c) Transmittance of the R and L waves, TR (red) and TL (black), as a function of the magnetic field strength for the cases of (b) ne/nc = 10 and (c) ne/nc = 1.

III. NUMERICAL SETUP the x boundaries, whereas the periodic boundary con- ditions are adopted in the y direction. The spatial and −3 The interaction between a short-pulse laser and a dense time resolutions are ∆x = ∆t = 5 × 10 and the par- ticle number is 100 per each grid at the beginning. The plasma foil is examined in our numerical experiment. We e e consider a simple situation in which a thin foil target of third order interpolation algorithm is used to suppress solid hydrogen is irradiated by a Gaussian-shape laser. the numerical heating. The time resolution ∆t should The plasma target is assumed to be fully ionized ini- be sufficiently smaller than the electron gyration time as tially. The number densities of the and ions, well as the time to avoid the unphysi- cal numerical heating in strongly magnetized plasmas. ne = ni, are spatially constant inside the target. For simplicity, the foil is in contact directly with the vacuum The transmittance and reflectivity of the electromag- without any preplasma. The thickness of the target L netic waves are assessed by integrating the energy of the waves passing through each side of the x bound- is set to be the same as the laser wavelength λ0 in the vacuum. Hereafter the length scale is normalized by the ary. The calculations are continued long enough until the transmittance and reflectivity become constant with laser wavelength such as L = L/λ0 = 1. time, which is typically tend = 25 after the injected laser The injected laser is an LP light traveling in the x di- e hits the target. The key parameters of this system are rection. The pulse duration at half maximum amplitude e the target density ne, the magnetic field strength Bext, is τ0 = τ0/t0 = 3, where the unit of time scale is given and the laser amplitude a0. The fiducial values repre- by the laser period t0 = λ0/c and c is the speed of light. e e Thee incident angle of the laser is normal to the target senting the region (I) are selected as ne = 10, Bext = 7, and a = 1. The dimensional quantities of these fidu- surface, and the focal spot is W = 3 in diameter. The 0 cial parameters are estimated as a functione ofe the laser peak laser intensity I0 is parameterized by the normal- f wavelength λ0. For instance, the laser intensity is given ized vector potential a = eE /(m cω ) via the relation − 0 0 e 0 by I =1.4 × 1018a2 λ 2 W/cm2, where λ is the laser of I = ǫ cE2/2, where E is the laser electric field and 0 0 µm µm 0 0 0 0 wavelength in microns. The critical density is calculated ǫ0 is the vacuum permittivity. A uniform external mag- 27 −2 −3 as nc = 1.1 × 10 λµm m , and the critical magnetic netic field Bext is applied parallel to the laser propagation −1 direction. field strength is Bc = 11 λµm kT. For this case, the unit time corresponds to t =3.3 λ fs. Numerical simulations on the laser-plasma interaction 0 µm in Cartesian coordinates (x, y) are calculated by a PIC scheme, PICLS [25]. The collision term is ignored in this analysis, unless otherwise mentioned. The compu- IV. NUMERICAL RESULTS tational domain is set to be (Dx, Dy) = (13, 9), and the plasma target is located at the middle −0.5 ≤ x ≤ 0.5. The transmittance and reflectivity of the injected LP The escape boundary conditionse fore the electromagnetic laser light are evaluated by 2D PIC simulations. In this waves and plasma particles are applied at bothe side of section, the density of the plasma target is fixed to be 4

....6 5 4 3 ⇐ 2NRL/λ0 1 1 0.8

0.8 0

0 0.6 0.6 E/E E/E 0.4 0 4 . Energy

Energy 0.2 0.2 0 0 0.1 1 10 2 0.1 1 10 10 Foil Thickness L/λ0 Magnetic Field Strength Bext/Bc FIG. 3. Dependence of the transmittance (black circles) FIG. 2. Transmittance (black) and reflectivity (red) after and refectivity (red triangles) on the thickness of target foil L the interaction between laser and magnetized plasma target obtained by 1D PIC simulations. The kinetic energy of ions evaluated numerically by a series of PIC simulations. The ki- (green squares) and electrons (blue crosses) are also shown. netic energies of ions (green) and electrons (blue) at the end The parameters except for the foil thickness are the density of each run are also indicated. The energy fractions are nor- ne/nc = 10, the magnetic field strength Bext/Bc = 7, and the malized by the incident laser energy E0. The plasma density is laser amplitude a0 = 1. The transmittance is independent of 0 > assumed to be ne/nc = 10 and the laser amplitude is a0 = 1. the foil thickness when L/λ ∼ 0.5. The circles are obtained from 2D runs (120 runs). The same calculations performed in the 1D and 3D geometries are also shown by the thin solid curves (600 runs) and cross marks (1 run), respectively. The gray dashed curve is the theoretical surface is taken into account in the theoretical curve. The prediction of the transmittance. The interference condition numerical results of the transmittance seem to oscillate < < 2NRL/λ0 is indicated at the top of the figure, where NR has around the predicted curve in the range of 1 ∼ Bext ∼ a dependence of Bext. 10. The oscillation is due to the interference by multiple reflections inside the target. The oscillation peakse in the transmittance appear when an integer multiple of the whistler wavelength (λ0/N ) coincides with the optical ne = 10, and then we will seek conditions that only the R whistler wave is selectively transmitted. The foil thick- path difference of 2L (see the labels at the top of Fig. nesse is L = 1, and the amplitude and pulse duration of 2). When the target thickness becomes much thicker > the laser are a0 = 1 and τ0 = 3. After enough time has than the whistler wavelength (e.g., L ∼ 3), the oscillatory e behavior disappears, and the numerical results exhibit passed since the laser-plasma interaction, the incident e laser energy should be transferrede to four different forms: an excellent agreement with the theory. If the target the transmitted and reflected electromagnetic-wave ener- is sufficiently thicker than the electron skin depth c/ωpe, gies and the absorbed ion and electron kinetic energies. the transmittance feature is independent of L (see Fig. 3). The outcome of the energy redistribution after the In this analysis, we look for the physical parameters of laser-plasma interaction depends on the strength of the a plasma device that is transparent for the R wave and external magnetic field. The transmittance and reflec- entirely reflective for the L wave. It is confirmed that tivity in terms of the incident laser energy are depicted the ideal situation as for a CP wave generator is real- by the black and red circles in Fig. 2. The final state of ized when the magnetic field strength is just below the the ion and electron kinetic energies are also shown by L-cutoff condition, that is Bext ∼ 9. Figure 4 shows snap- the green and blue circles, respectively. The energy con- shot images of the electromagnetic waves in one of the version to plasmas is not significant except for the cases best parameters as the plasmae device. The field strength around Bext ∼ 1, where the electrons absorb the wave of this run is Bext = 7. energy through the cyclotron resonance [26]. Then, most e The spatiale distribution of the perpendicular compo- of the laser energy remains in the form of electromagnetic nent of the electric field E⊥ is indicated in Fig. 4(a). The waves. target foil is located at the center, and the incident laser The numerically obtained transmittance is consistent comes from the left. The energy conversion efficiency with the theoretical estimation shown by the gray dashed from the injected LP light to the transmitted wave is curve in Fig. 2. Because an LP light is a sum of an RCP 0.47. At the time of the snapshot, the transmitted wave light and an LCP light, each having a half component, and the reflected wave exist almost symmetrically to the the total transmittance at a given Bext is basically the target. Figure 4(b) is a 1D distribution of the electric average of T and T shown in Fig. 1(b). Note that the field along the x-axis. The phase delay of π/2 between R L e effect of multiple reflections on both sides of the target Ey and Ez indicates that the transmitted wave is the R 5

4 E

0 ⊥

λ 2

/ (c) y 0 0 5 . 0.6

−2 0.4 0 E Position /

4 (a) 0 2 ⊥ E 0 − 0 . E E / / 0

z 0 0.5 z E Ey Ez E 0

E 0 /

E −0.5 0.5 0 E /E0 −15 y −0.5 (b) −10 −5 0 −0.5 5 10 15 −15 −10 −5 0 5 10 15 Position x/λ0 Position x/λ0

2 2 1/2 FIG. 4. (a) Spatial distribution of the perpendicular electric field E⊥ =(Ey + Ez ) of the transmitted (x> 0) and reflected (x < 0) waves, where the plasma foil is located at −0.5 ≤ x/λ0 ≤ 0.5 (gray dashed line). In this 2D fiducial run, the plasma density is ne/nc = 10, the magnetic field strength is Bext/Bc = 7, and the laser amplitude is a0 = 1. The pulse duration of the incident laser is τ0/t0 = 3 and the focal spot size at the foil surface is W/λ0 = 3 in the y-direction. This snapshot is taken at a time t/t0 = 11.2 after the incident laser reaches the target surface. Only for this run, a wider domain size in the x direction is used (Dx/λ0 = 31). (b) 1D distributions of the electric fields, Ey (black) and Ez (red), along the x-axis extracted from the same snapshot shown in (a). (c) 3D plot of the direction of E⊥ using the same data as in (b). The color denotes the amplitude of the perpendicular electric field. wave and the reflected wave is the L wave. The ampli- The laser-plasma interaction considered here is inde- tude of the transmitted and reflected wave is about half pendent of the spatial dimensions of the simulations. In of the incident laser amplitude, aT ∼ aR ∼ a0/2. The Fig. 2, the 1D and 3D PIC results are plotted by the polarization is mostly circular, but slightly elliptical. For solid curves and cross marks, respectively. As seen from example, the deviation of the amplitude between Ey and the figure, all the data exhibit the same trend, and which Ez is around 20% in the transmitted wave, which might indicates that the effect of the finite laser spot must be be due to the contamination of components other than negligible. The numerical resolution of the 1D runs is the whistler-wave mode. Figure 4(c) is a 3D plot of the identical to that of the 2D simulations. The conclusions direction of the E⊥ vector shown in Fig. 4(b). It is evi- illustrated by Fig. 2 are unaffected by the numerical res- dent that two CP waves are generated successfully after olutions, which has been tested by the higher resolution the laser-plasma interaction. runs in 1D with ∆x = 10−3. It is also verified by 1D runs Note that both of the transmitted and reflected CP that the inclusion of the collisional effects has little in- e waves are RCP waves with respect to the propagation fluence on the transmittance in our system. The fiducial run is also performed in 3D, where the grid resolution is direction. The magnetic field direction is parallel to the −2 laser injection. The transmitted wave is the R wave so ∆x = 10 and fewer particle numbers (20 per grid) are that the circular polarization is right-hand to the propa- solved. These numerical results demonstrate that the dis- e gation direction. The reflected wave, on the other hand, covered feature of this plasma device is robust and quite is the L wave, and then it has left-hand polarity to the practical. magnetic field direction. Because the propagation direc- The most significant advantage of this plasma device tion is opposite to the magnetic field, the reflected wave is is that it is applicable even for the relativistic amplitude also right-hand to the propagation direction. Then, this of the wave. If the amplitude of the injected wave is rel- plasma device generates two RCP waves propagating in ativistic, the transmitted wave amplitude could also be the opposite direction. It might be possible to establish nearly relativistic. Figure 5 denotes the dependence of an optics system for combining these two CP lasers in the transmittance on the laser amplitude a0. The fidu- practical experiments. cial parameters, ne = 10 and Bext = 7, are adopted for > Interestingly, if the magnetic field direction is anti- these 2D simulations. When a0 ∼ 10, some particles have escaped from thee computationale domain, but those are parallel to the laser propagation (i.e., Bext = −7), two LCP waves will be formed by the same process. There- not counted in the energies of this figure. The escaped fore, the polarity of the generated CP wavese can be se- energy is at most 10% of the injected laser energy. lected by the direction of the external magnetic field. The transmittance is independent of the laser inten- 6

1 V. DISCUSSION

0 8 . For practical use of this plasma photonic device, it

0 might be difficult to align the magnetic field direction 0.6 precisely the same as the laser direction. Then, it would E/E be important to examine the angle dependence of the 0.4 transmittance. The laser injection is kept to be along the x-axis in the following analysis. We define ϕ as an Energy 0.2 angle between the laser and magnetic field directions in the unit of degree. Here the angle is assumed to vary 0 within the x-y plane. The angle ϕ = 90 (0) denotes that 0.1 1 10 102 the magnetic field direction is along the y-axis (x-axis). 2D PIC simulations reveal the angle dependence of the Wave Amplitude a0 efficiency of our plasma device. The angle dependence of the phase velocity in the region (I) is characterized by a FIG. 5. Dependence of the transmittance (black) and re- figure-eight shape of the wave normal surface [24]. Thus, flectivity (red) on the laser intensity. The kinetic energies absorbed by ions (green) and electrons (blue) are also shown. the transmittance of the R-wave branch should have a The energy fractions obtained from 2D PIC simulations are dependence on the angle ϕ. The polarization angle of given by the ratio to the incident laser energy E0. The fiducial the incident LP wave is another independent parameter, E e E e parameters are used for the plasma density ne/nc = 10 and so that we examine two cases, 0 k y and 0 k z. In the magnetic field strength Bext/Bc = 7. the former cases, the electric field of the incident laser has a parallel component to the external magnetic field. However, the laser electric field is always perpendicular to Bext in the latter cases. < Figures 6(a) and 6(b) show the angle dependence of sity until a0 ∼ 3. The polarization property of the transmitted wave is also unchanged for this range of the energy redistribution for the cases that the initial electric field is parallel to the y-axis and z-axis, respec- a0. These facts confirm the applicability of this mech- anism in the relativistic intensity regime. It should be tively. Except for the angle ϕ, the model parameters are noticed that even when the amplitude of the transmit- identical to the fiducial run shown by Fig. 4, in which the plasma density is ne = 10, the magnetic field strength is ted whistler wave is relativistic, aT =2a0/(NR + 1) > 1, the electron quiver velocity remains still non-relativistic, Bext = 7, and the incident laser amplitude is a0 = 1. The fiducial run correspondse to the angle ϕ = 0 case. vq/c = aT /(Bext − 1). It is because the magnetic field e considered here is sufficiently large enough. It is found that the transmittance feature is roughly e independent of the angle ϕ when the angle is less than If the laser amplitude becomes further higher, there are about 15 degrees. These results suggest the angle tol- two relativistic effects that prevent the CP wave genera- erance for the installation of this plasma device. The tion. One is the electron heating through the cyclotron transmittance of the R-wave branch decreases with the resonance. The external magnetic field in these runs is increase of the angle. When the angle is ϕ = 90, the far beyond the resonance condition. However, due to the injected wave is totally reflected. As a reference, the the- relativistic correction, the resonance condition becomes oretical curve is indicated by the gray dashed curve in ω0 ∼ ωce/γ, where γ is the Lorentz factor of electrons. Fig. 6, which corresponds to half of the transmittance for The electrons in the target gain the energy efficiently the R-wave branch. There is a critical angle ϕc = 69.9 > 2 1/2 where the transmittance, or the phase velocity of the through the resonance when a0 ∼ (Bext − 1) ∼ 7 [26]. In fact, the electromagnetic energy entered into the tar- electromagnetic wave, becomes zero [24]. When the elec- get is largely absorbed by the electrons,e and the transmit- tric field has the y component [Fig. 6(a)], a fairly large < < amount of the plasma heating occurs near the critical tance decreases drastically in the range of 7 ∼ a0 ∼ 10. angle ϕc. The parallel component of the electric field in The sudden rise in the transmittance at a0 ∼ 10 is terms of the external magnetic field is essential for the caused by the other effect, that is the relativistic trans- resonant absorption by electrons for this case. The 1D parency [27, 28]. The cutoff density for the relativistic- and 3D PIC results in Fig. 6 make sure that the angle amplitude laser is approximately given by γnc. The min- dependence is unaffected by the spatial dimensions of the imum intensity for the relativistic transparency is ob- simulations. > 2 1/2 tained by a0 ∼ [2(ne − 1)] ∼ 14, where the quiver en- In Sec. IV, the transmittance of the R wave is exam- 2 1/2 ergy γ = (1+a0/2) is assumed. This is also consistent ined intensively. However, as discussed in Sec. II, the with the numericale results shown in Fig. 5. Therefore, extraction of only the L wave may also be possible in the maximum intensity valid for the CP wave generator the parameter region (II). A benefit of the L-wave fil- could be managed by the field strength and the target ter is that the required strength of the magnetic field is density. relatively smaller, so that it could be more feasible for 7

1 1 (a) (b) 0.8 0.8 0 0 0.6 0.6 E/E E/E

0.4 0.4 Energy Energy 0.2 0.2

0 0 0 15 30 45 60 75 90 0 15 30 45 60 75 90 Angle ϕ [deg] Angle ϕ [deg]

FIG. 6. Angle dependence of the transmittance (black), reflectivity (red), and the energy fraction of ions (green) and electrons (blue). The definition of ϕ is the angle between the laser direction and the magnetic field orientation. The fiducial parameters are used for the plasma density ne/nc = 10, the magnetic field strength Bext/Bc = 7, and the laser amplitude a0 = 1. The laser-plasma interaction is affected by the electric field direction of the incident LP laser. The results when the laser electric field is in the y-direction and z-direction are shown in (a) and (b), respectively. The indicated marks are the same as in Fig. 2.

1 1 (a) (b) 0.8 0.8 0 0 0.6 0.6 E/E E/E

0.4 0.4 Energy Energy 0.2 0.2

0 0 0 15 30 45 60 75 90 0 15 30 45 60 75 90 Angle ϕ [deg] Angle ϕ [deg]

FIG. 7. Angle dependence of the transmittance (black), reflectivity (red), and the energy fraction of ions (green) and electrons (blue). The parameters representing the region (II) in Fig. 1(a) are chosen to be the plasma density ne/nc = 1, the magnetic field strength Bext/Bc = 0.7, and the laser amplitude a0 = 0.1. The results when the laser electric field is in the y-direction and z-direction are shown in (a) and (b), respectively. The indicated marks are the same as in Fig. 4. practical use. On the other hand, a defect of the L wave When the laser injection is parallel to the magnetic field extraction is a weaker limit in the laser intensity. Because (ϕ = 0), both of the transmittance and reflectivity are > of the relativistic transparency, the cases of a0 ∼ 1 are in- nearly half of the injected energy. In this case, only the appropriate for the LCP wave generation by our method. L-wave component is transparent, and thus it works as Furthermore, it turned out that there is a severe angle a converter to LCP waves from an LP laser. However, dependence in the transmittance of the L wave. The if the angle is not zero, the energy absorption by the ideal transmittance of the L wave is realized only when electrons cannot be ignored. If underdense preplasmas the laser injection is strictly parallel to the external mag- form at the ablation side of the target, the electrons are netic field. Otherwise, a large fraction of the laser energy accelerated effectively by the laser field through the v×B goes to the thermal energy of electrons. force. The wave propagation is allowed at any angle for Figure 7 shows the angle dependence of the energy this case, which is the L wave at ϕ = 0 and the X wave fraction after the laser-plasma interaction for a typical at ϕ = 90. The angle dependence of the transmittance case in the region (II). The model parameters calculated and reflectivity exhibits complicated behavior, as seen in Figs. 7(a) and 7(b). here is the plasma density ne = 1, the magnetic field strength Bext = 0.7, and the incident laser amplitude In this paper, the incident angle of the laser is assumed a0 =0.1. The other parameterse are the same as in Fig. 6. to be always normal to the target surface. Even when e 8

1 der a strong magnetic field. When a thin plasma foil is irradiated by an LP laser under some appropriate condi- 0.8 tions, the transmitted and reflected lights become cir-

0 cularly polarized and have nearly the same intensity. 0.6 The whistler-mode of the R wave passes through a thin E/E < < plasma target when 1 ∼ Bext ∼ ne − 1, while the L- 0.4 wave component is totally reflected. To function as a conversion filter, the densitye of thee plasma foil should Energy be higher than n > 2. The optimal condition requires a 0.2 e strong external magnetic field with a size of Bext ∼ ne −1 along the incidente light direction. The polarity of the CP 0 e e 1 10 102 waves can be switched by the direction of the magnetic field. The transmittance is unaffected by a small angle Pulse Duration τ0/t0 difference between the wave injection and the magnetic field less than about 15 degrees. FIG. 8. Dependence of the transmittance and refectivity on the pulse duration τ0 of the incident laser obtained by 1D PIC A significant result of this work is that linear anal- simulations. The indicated marks are the same as in Fig. 3. ysis of this problem is quite adequate for assessing the The fiducial parameters are used in these runs, but the foil behavior of this concept even for the relativistic laser in- thickness is assumed to be L/λ0 = 10. The transmittance decreases dramatically if the pulse duration becomes longer tensity. The advantage of this plasma photonic device about τ0/t0 ∼> 10. is a large amplitude of CP wave, which is valid even for the relativistic amplitude. Therefore, the high-intensity CP waves generated by this plasma device could be use- the laser has an incident angle to the target, if the mag- ful for high-energy-density plasma experiments, such as netic field is parallel to the laser direction, the transmit- laser-driven ion acceleration [8, 9] and ion heating in over- ted light should have circular polarization. It might be dense plasmas [10, 11]. It would be possible to propose interesting to investigate the dependence of the trans- plasma photonic devices for various different uses by ap- mittance on the laser incident angle together with the ϕ propriately selecting the plasma environment for the den- dependence. We ignore the effects of preplasmas in this sity and magnetic field strength. In the future, it will analysis, because the preplasma has little influence on the be meaningful to pursue further device developments for propagation of the whistler wave. Throughout this anal- ultra-intense lasers while utilizing numerical simulations. ysis, we assume a relatively shorter pulse duration for the injected laser, τ0 = 3. However, the stimulated Brillouin scattering could prevent whistler-wave propagation if the pulse duratione becomes longer (see Fig. 8). The critical pulse length is determined by the growth time of the Bril- louin instability, and thus it would depend largely on the laser intensity [11]. The dependence of the transmittance ACKNOWLEDGMENTS on the pulse duration would also be an exciting task to be studied. We thank Ryosuke Kodama, Kunioki Mima, and Yoshitaka Mori for useful discussions. This research VI. CONCLUSIONS was partially supported by JSPS KAKENHI Grant No. JP19KK0072, JP20H00140, and JSPS Core-to- We have investigated a unique mechanism of CP wave Core Program, B. Asia-Africa Science Platforms No. generation by controlling the laser-plasma interaction un- JPJSCCB20190003.

[1] R. Kodama, Y. Sentoku, Z. L. Chen, G. R. Kumar, [2] H. C. Kapteyn, A. Szoke, R. W. Falcone, and M. M. S. P. Hatchett, Y. Toyama, T. E. Cowan, R. R. Free- Murnane, Opt. Lett. 16, 490 (1991). man, J. Fuchs, Y. Izawa, M. H. Key, Y. Kitagawa, [3] P. Michel, L. Divol, E. A. Williams, S. Weber, C. A. K. Kondo, T. Matsuoka, H. Nakamura, M. Nakatsut- Thomas, D. A. Callahan, S. W. Haan, J. D. Salmonson, sumi, P. A. Norreys, T. Norimatsu, R. A. Snavely, R. B. S. Dixit, D. E. Hinkel, M. J. Edwards, B. J. Mac- Stephens, M. Tampo, K. A. Tanaka, and T. Yabuuchi, Gowan, J. D. Lindl, S. H. Glenzer, and L. J. Suter, Nature 432, 1005 (2004). Phys. Rev. Lett. 102, 025004 (2009). 9

[4] J. D. Moody, P. Michel, L. Divol, R. L. Berger, E. Bond, [17] S. Fujioka, Z. Zhang, K. Ishihara, K. Shigemori, D. K. Bradley, D. A. Callahan, E. L. Dewald, S. Dixit, Y. Hironaka, T. Johzaki, A. Sunahara, N. Yamamoto, M. J. Edwards, S. Glenn, A. Hamza, C. Haynam, D. E. H. Nakashima, T. Watanabe, H. Shiraga, H. Nishimura, Hinkel, N. Izumi, O. Jones, J. D. Kilkenny, R. K. Kirk- and H. Azechi, Sci. Rep. 3, 1170 (2013). wood, J. L. Kline, W. L. Kruer, G. A. Kyrala, O. L. [18] J. J. Santos, M. Bailly-Grandvaux, M. Ehret, A. V. Landen, S. LePape, J. D. Lindl, B. J. MacGowan, N. B. Arefiev, D. Batani, F. N. Beg, A. Calisti, S. Ferri, Meezan, A. Nikroo, M. D. Rosen, M. B. Schneider, D. J. R. Florido, P. Forestier-Colleoni, S. Fujioka, M. A. Strozzi, L. J. Suter, C. A. Thomas, R. P. J. Town, Gigosos, L. Giuffrida, L. Gremillet, J. J. Honrubia, K. Widmann, E. A. Williams, L. J. Atherton, S. H. Glen- S. Kojima, P. Korneev, K. F. F. Law, J.-R. Marqu`es, zer, and E. I. Moses, Nat. Phys. 8, 344 (2012). A. Morace, C. Moss´e, O. Peyrusse, S. Rose, M. Roth, [5] V. M. Malkin, G. Shvets, and N. J. Fisch, S. Sakata, G. Schaumann, F. Suzuki-Vidal, V. T. Phys.l Rev. Lett. 82, 4448 (1999). Tikhonchuk, T. Toncian, N. Woolsey, and Z. Zhang, [6] Y. Shi, H. Qin, and N. J. Fisch, Phys. Plasmas 25, 056705 (2018). Phys. Rev. E 95, 023211 (2017). [19] K. G. Budden, Radio waves in the ionosphere (Cam- [7] R. L. Carman, D. W. Forslund, and J. M. Kindel, bridge University Press, New York, 1961). Phys. Rev. Lett. 46, 29 (1981). [20] X. H. Yang, W. Yu, H. Xu, M. Y. Yu, Z. Y. Ge, B. B. [8] T. Esirkepov, M. Borghesi, S. V. Bulanov, G. Mourou, Xu, H. B. Zhuo, Y. Y. Ma, F. Q. Shao, and M. Borghesi, and T. Tajima, Phys. Rev. Lett. 92, 175003 (2004). Appl. Phys. Lett. 106, 224103 (2015). [9] F. Fiuza, A. Stockem, E. Boella, R. A. Fon- [21] S. X. Luan, W. Yu, F. Y. Li, D. Wu, Z. M. Sheng, M. Y. seca, L. O. Silva, D. Haberberger, S. Tochit- Yu, and J. Zhang, Phys. Rev. E 94, 053207 (2016). sky, C. Gong, W. B. Mori, and C. Joshi, [22] G. Ma, W. Yu, M. Y. Yu, S. Luan, and D. Wu, Phys. Rev. Lett. 109, 215001 (2012). Phys. Rev. E 93, 053209 (2016). [10] T. Sano, M. Hata, D. Kawahito, K. Mima, and Y. Sen- [23] S. Weng, Q. Zhao, Z. Sheng, W. Yu, S. Luan, M. Chen, toku, Phys. Rev. E 100, 053205 (2019). L. Yu, M. Murakami, W. B. Mori, and J. Zhang, [11] T. Sano, S. Fujioka, Y. Mori, K. Mima, and Y. Sentoku, Optica 4, 1086 (2017). Phys. Rev. E 101, 013206 (2020). [24] T. H. Stix, Waves in plasmas (Springer-Verlag, New [12] P. Michel, L. Divol, D. Turnbull, and J. D. Moody, York, 1992). Phys. Rev. Lett. 113, 205001 (2014). [25] Y. Sentoku and A. J. Kemp, [13] P. Michel, E. Kur, M. Lazarow, T. Chapman, L. Divol, J. Comp. Phys. 227, 6846 (2008). and J. S. Wurtele, Phys. Rev. X 10, 021039 (2020). [26] T. Sano, Y. Tanaka, N. Iwata, M. Hata, [14] G. Lehmann and K. H. Spatschek, K. Mima, M. Murakami, and Y. Sentoku, Phys. Rev. Lett. 116, 225002 (2016). Phys. Rev. E 96, 043209 (2017). [15] G. Lehmann and K. H. Spatschek, [27] P. Kaw and J. Dawson, Phys. Fluids 13, 472 (1970). Phys. Rev. E 97, 063201 (2018). [28] H. Sakagami and K. Mima, [16] D. J. Stark, C. Bhattacharjee, A. V. Arefiev, Phys. Rev. E 54, 1870 (1996). T. Toncian, R. D. Hazeltine, and S. M. Mahajan, Phys. Rev. Lett. 115, 025002 (2015), 1412.1865.