 known process that was found in 1D kinetic simulations [@araneda2008proton; @MatteiniEA2010]. Here we investigate the origin and its possible relation to the strong perpendicular heating that we observe in our simulations.

To explore the plasma heating and the proton beam formation, we have performed the same numerical experiment that consists of a non-monochromatic initial perturbation with different plasma beta. Although the monochromatic case provides good insights into the physics of parametric instabilities, it remains an idealized problem. Previous numerical simulations using MHD and kinetic models [@malara2; @anna1; @MatteiniEA2010; @primavera_2019] have explored the parametric instability using non-monochromatic Alfvén wave packets in 1D and 2D configurations and it is known that it continues operating with a slight reduction of the growth rate, essentially because a broad-band spectrum tends to quench coherent resonances [@malara_Phys_fluids_1996].

We have run 1D, 2D and 3D simulations for the non-monochromatic wave using the same numerical parameters and we have obtained similar results for the three configurations (not shown).

![Comparison of different macroscopic quantities for simulations labeled with run C (**Top**). The normalized standard deviation of density fluctuations. (**Middle**). Parallel magnetic fluctuations and (**Bottom**) total temperature difference for different plasma beta condition.[]{data-label="Fig8"}](Fig8-eps-converted-to.pdf){width="50.00000%"}

The results from the 2D simulations with $\delta b_0=1$ and different plasma beta are summarized in Fig. \[Fig8\]. Density fluctuations are larger for low beta plasma and become progressively smaller as the beta increases. The amplitude of parallel magnetic field fluctuations seems not to be affected by the plasma beta, but the time for them to exponentially grow is affected by the plasma beta.

The proton heating shows different stages also in these non-monochromatic cases. There is an initial strong temperature increase and a beam formation, that correlates with the enhancement of the parallel magnetic field fluctuations. A secular stage is followed where the heating rate is significantly reduced. The final temperature gained by the particles is similar to the monochromatic cases and so are the trends with the plasma beta and pump wave amplitude of $\Delta T$ (and $T_\bot$, $T_\parallel$). Such trends are not affected by the plasma beta, rather, they are
known process that was found in 1D kinetic simulations   [ @araneda2008proton; @MatteiniEA2010 ]. Here we investigate the beginning and its potential relation to the strong vertical heating that we respect in our simulations. 

 To research the plasma heating and the proton beam formation, we have perform the same numerical experiment that consist of a non - monochromatic initial perturbation with different plasma beta. Although the monochromatic sheath provide good insights into the physics of parametric instabilities, it remains an idealized trouble. Previous numerical simulations use MHD and kinetic models   [ @malara2; @anna1; @MatteiniEA2010; @primavera_2019 ] have explored the parametric imbalance using non - monochromatic Alfvén wave packet in 1D and 2D configurations and it is know that it continues operating with a slight decrease of the growth rate, essentially because a broad - band spectrum tends to snuff out coherent resonances   [ @malara_Phys_fluids_1996 ]. 

 We have run 1D, 2D and 3D simulations for the non - monochromatic wave using the same numerical parameters and we have obtained similar results for the three configuration (not shown). 

! [ Comparison of different macroscopic quantity for simulations labeled with run C (* * Top * *). The anneal standard deviation of density fluctuations. (* * Middle * *). Parallel charismatic fluctuations and (* * Bottom * *) total temperature difference for different plasma beta condition.[]{data - label="Fig8"}](Fig8 - eps - converted - to.pdf){width="50.00000% " } 

 The results from the 2D simulations with $ \delta b_0=1 $ and different plasma beta are summarized in Fig.   \[Fig8\ ]. Density fluctuations are large for low beta plasma and become increasingly little as the beta increases. The amplitude of parallel magnetic field fluctuation seems not to be affected by the plasma beta, but the time for them to exponentially originate is affected by the plasma beta. 

 The proton heating shows different stage also in these non - monochromatic case. There is an initial potent temperature increase and a radio beam formation, that correlate with the enhancement of the parallel magnetic sphere fluctuations. A secular stagecoach is stick to where the heating rate is significantly reduced. The final temperature gained by the particles is similar to the monochromatic casing and so are the trends with the plasma beta and pump wave amplitude of $ \Delta T$ (and $ T_\bot$, $ T_\parallel$). Such vogue are not affected by the plasma beta, rather, they are
 knlwn process that was foukd in 1D kinetic simulatmons [@araheda2008protun; @MatteiniEA2010]. Here we investmgatw the origin and its possibue relatiln to thw stcong perpendiculed heatiky thaf we muserve in our slmulations.

Tm explore the [lxsla heating and the proton beam formwtion, wr jave performed the fame numerical experiment that consista of a ion-monochromativ initial perturbation witj divferent plasma betw. Although jge iinochromatic case provpbes good inaights into the physics of paraoetrie instabilijnws, lj remains an ideajized problem. Previogs numetical simulatipns usung MHD and kinetic mmdels [@malara2; @anna1; @MwtteiniEA2010; @pximavera_2019] have exploree rhe pdramatriz inrtagikify usijg ion-monochrojatic Alfvén wave packets in 1D snq 2D configuratjons agd it is known that it continues operatinc wjth a slight reduction if the growth rate, esdentially because a broad-band spectrum tends to quench coharent fesinwvxed [@malara_Phys_fluids_1996].

We have run 1D, 2D and 3D simulaejoms for the non-mokochromatic wave uxijg jhe same numerkcal pcdajeters and we have obtaingd simular resujts gor the three configurationw (not shown).

![Cjnparison of differznt macrosco'ic qusntitoes for simulations labzled wjth run C (**Tlp**). The nodoalized standard deniathon of dtvsity fluctuationf. (**Middle**). 'aralkel magvetiv fluceuations ajd (**Bobdom**) total temperatkre dnfferance for dlfferent plasma beta condition.[]{deva-label="Fig8"}](Fig8-gps-wonnerted-to.pbf){widtm="50.00000%"}

The results frjm the 2D simulctions wnth $\deuta b_0=1$ and siffereit plasma beea are summarheed in Fig. \[Fig8\]. Density fluxtuarions afd larger for lpw beta plasma and vecome progressiveky rjaller as the bzua uncreases. The akplktuqe oh parwnlel magnetiw fidld gluctjations seenr noy to be affected by dhe llasma beta, but thr bime for jhem to evponentially brow is affected bj the plesma brta.

Jhe proton heating shows differsnt staged ajso in these non-nonochromatie cases. There is an initial strong tempecature increase and a bwam formation, that eovrelates witi the qnhancemett of the parallel mqgnetic field flugtuations. A secular stzge is folllwed where the heating rate is significantly reduced. The final temperaturw gainxd by the parficlrs is snmikar to the oonochromatic caxes and so are the trends with vhe plasma bzta and pump wave amplitude og $\Aelta T$ (and $T_\cot$, $T_\parallel$). Such trenda are noy affected by the plasma beta, rsther, they are
known process that was found in 1D [@araneda2008proton; Here we the origin and strong heating that we in our simulations. explore the plasma heating and the beam formation, we have performed the same numerical experiment that consists of a initial perturbation with different plasma beta. Although the monochromatic case provides good insights the of instabilities, remains an idealized problem. Previous numerical simulations using MHD and kinetic models [@malara2; @anna1; @MatteiniEA2010; @primavera_2019] explored the parametric instability using non-monochromatic Alfvén wave in 1D and 2D and it is known that continues with a reduction the rate, essentially because broad-band spectrum tends to quench coherent resonances [@malara_Phys_fluids_1996]. We have run 1D, 2D and 3D simulations for non-monochromatic wave same numerical and have similar results for configurations (not shown). ![Comparison of different simulations labeled with run C (**Top**). The normalized deviation of fluctuations. (**Middle**). Parallel magnetic fluctuations and total temperature difference for different plasma beta condition.[]{data-label="Fig8"}](Fig8-eps-converted-to.pdf){width="50.00000%"} results from the 2D simulations with $\delta b_0=1$ and different plasma beta are summarized in Density fluctuations are larger low beta plasma become smaller the increases. The of parallel magnetic field fluctuations seems not to be affected by plasma beta, but the time for them to exponentially grow by plasma beta. The heating shows different stages in non-monochromatic cases. There is strong increase formation, correlates the enhancement of the magnetic field fluctuations. A secular is followed where the The final temperature gained by the particles is to the monochromatic cases and so are trends with the plasma beta and pump wave amplitude of $\Delta T$ $T_\bot$, $T_\parallel$). are not affected by the plasma beta, rather, are
 known process that was found iN 1D kinetic sImulaTioNs [@aRaNeda2008ProtOn; @MatteiniEA2010]. HeRE we iNvestigate the origin and Its poSsIBle rELaTion tO the strONg PERpeNdIcUlaR hEAtIng thAt wE observE in our simuLatIoNs.

To explore tHE pLasma heatiNg aNd the proton bEam FormatIoN, we HAve peRfoRmed tHe same NUmericAl experimEnT That coNSists of A NOn-MonoChromatic initial pERtURbation with difFerent PlASmA BEta. altHough the moNoChromATic case PRoVIDEs gOOd insights intO the physics OF paRametrIc InsTAbilitIes, it ReMAinS an idealizeD proBlem. PreviOus numERical siMUlationS using mHD And KineTIc MoDelS [@mALarA2; @AnNa1; @MATteIniEA2010; @priMaVeRa_2019] havE expLORED the ParAmetRic inStability usinG noN-monOChrOmatiC AlfvÉn waVe PackeTs in 1D aNd 2D coNfIgurations and it Is knOwn that it ConTiNueS oPeratINg with A slIghT reductIon of thE GroWtH RATe, Essentially because A bROAd-Band specTrum teNDs To QUench cohErEnt ResoNANces [@mAlarA_phYs_fluids_1996].

we have RUn 1d, 2D And 3D simUlAtions FoR thE noN-monoCHromAtic waVe using tHe samE Numerical paramETers and we have OBtAINeD SimiLar Results for tHe thREe coNfigURaTioNS (not sHown).

![COmPArISon of different macroScOpic quAntitIes for simulatIons labeleD WITh run C (**ToP**). The NOrMAlized standard DeviaTion of densITy fluctuAtionS. (**Middle**). PArallel maGNEtic flucTuaTioNs aNd (**BOTToM**) total temperaTURe diFfErence fOr dIfferenT plAsmA beTa cOnDition.[]{datA-label="FiG8"}](FIg8-EpS-cOnvErted-TO.pdf){widtH="50.00000%"}

THe rEsUltS from THe 2D simUlatiOns wItH $\dELta B_0=1$ and difFErENT plaSmA bEta aRe sUmMarizEd in fIg. \[FIg8\]. DensiTy fluctuaTioNS are LaRgEr for loW beta plasma anD bEcome progrEsSivEly smaLLEr as the bEta increases. The amplitudE Of paralLel MagneTic fIeld fluctUatIons seEms NOt to be AffectEd by tHe PlaSMA beta, BUT tHe tImE for them to EXPonEntiaLlY groW is affeCted by the plasma betA.

the Proton heating ShoWs diFFErEnt STaGEs aLsO In tHESe non-monochromaTic cases. ThErE Is An initial sTRonG tEmperatUre incrEase aND a beam fOrmation, tHat correlAtEs wiTH The EnhancemenT of the paRallel magNEtic fIElD flucTuaTions. A SeCulAr staGe is foLLowEd wheRe the hEaTing raTe is sIgNificantLy reduced. The final temperAture gAined By tHe particlEs iS SimIlar to the MonoChromatic cAseS anD so arE thE TrendS witH ThE plASma beTa anD Pump wave aMPlItuDE Of $\delta T$ (and $T_\bOT$, $t_\ParAllel$). sucH Trends Are nOt affected by the plASma beta, rather, tHey aRE
 known process that was fo und in 1Dkinet icsim ul atio ns [ @araneda2008pr o ton;  @MatteiniEA2010]. Her e wein v esti g at e the  origin an d its  p os sib le re latio n t o the s trong perp end ic ular heating th at we obse rve  in our simu lat ions.
T o e x plore  th e pla sma he a ting a nd the pr ot o n beam formati o n ,we h ave performed the sa m e numerical ex perime nt th a t  co nsi sts of a n on -mono c hromati c  i n i t ial perturbationwith differ e ntplasma  b eta .  Altho ugh t he mon ochromaticcase  provides  goodi nsights into th e phys ics  of  par a me tr icin s tab i li tie s , i t remain san  idea lize d p r oble m.Prev iousnumerical sim ula tion s  us ing M HD an d ki ne tic m odels[@mal ar a2; @anna1; @Ma ttei niEA2010;  @p ri mav er a_201 9 ] have  ex plo red the  parame t ric  i n s t ab ility using non-mo no c h ro matic Al fvén w a ve  p a ckets in  1 D a nd 2 D confi gura t io ns and i t is k n ow nthat it  c ontinu es  op era tingw itha slig ht reduc tiono f the growth r a te, essential l yb e ca u se a  br oad-band sp ectr u m te ndst oque n ch co heren tr es o nances [@malara_Phy s_ fluids _1996 ].

We have r un 1D, 2Da n d  3D simu lati o ns for the non-mo nochr omatic wav e  using t he sa me numer ical para m e ters and  we  ha veobt a i ne d similar res u l ts f or  the th ree  config ura tio ns(no tshown).
![Compar is on  o fdif feren t  macrosc op icqu ant ities for si mulat ions  l ab e led  with r u nC (**T op ** ). T heno rmali zeds tan dard de viation o f d e nsit yfl uctuati ons. (**Middl e* *). Parall el  ma gnetic f luctuati ons and (**Bottom**) to t al temp era turediff erence fo r d iffere ntp lasmabeta c ondit io n.[ ] { data- l a be l=" Fi g8"}](Fig8 - e ps- conve rt ed-t o.pdf){ width="50.00000%"} 
Th e results fro m t he 2 D si mul a ti o nswi t h $ \ d elta b_0=1$ and  different  p l as ma beta ar e  su mm arizedin Fig.  \[Fi g 8\]. De nsity flu ctuations  a re l a r ger  for low b eta plas ma and be c ome p r og ressi vel y smal le r a s the  betai ncr eases . Theam plitud e ofpa rallel m agnetic field fluctuati ons se ems n otto be aff ect e d b y the pla smabeta, butthe  ti me fo r t h em to  exp o ne nti a lly g rowi s affecte d  b y t h e  p lasma beta. T heproto n h e atingshow s different stage s  also in these  non - m ono chr o mati ccases. There i s a ni n itial st ro ng temperat ure incr ea s e and  a bea m form ation,t h at correl ates  wi th the en han ce m ent ofth ep aralle l ma gn etic f ield f l uctu a t ions. A secularstage i s fol l owe d whe re  the he a ting  rate is s ignificantl y redu ced.  Thefinal t em peratu rega ined by th e  particle s issimilar  t o th e m onochr omat i c  case s an dsoare the t r e nd s  w it h  th e pl asmabe ta a nd pump w a ve ampli tud e  of $\D el taT $  (and$ T_ \ b ot$, $T_\p ara llel$ ) .  Such tren d s ar e  n o t aff ectedby the  plasma bet a,  rather , t h e y are
 known_process that_was found in 1D_kinetic simulations [@araneda2008proton;_@MatteiniEA2010]._Here we_investigate_the origin and_its possible relation_to the strong perpendicular_heating that we_observe_in our simulations.

To explore the plasma heating and the proton beam formation, we have_performed_the same_numerical_experiment_that consists of a non-monochromatic_initial perturbation with different plasma_beta. Although_the monochromatic case provides good insights into the_physics_of parametric instabilities,_it remains an idealized problem. Previous numerical simulations using_MHD and kinetic models [@malara2; @anna1; @MatteiniEA2010;_@primavera_2019] have explored_the_parametric_instability using non-monochromatic Alfvén_wave packets in 1D and 2D_configurations and it is known that_it continues operating with a slight reduction_of the growth rate, essentially because_a broad-band spectrum tends to_quench coherent_resonances [@malara_Phys_fluids_1996].

We have run 1D, 2D_and 3D simulations_for the_non-monochromatic wave using_the same numerical parameters and we_have obtained similar_results for the three configurations (not_shown).

![Comparison_of different macroscopic_quantities_for_simulations labeled_with run C_(**Top**)._The normalized_standard_deviation of density fluctuations. (**Middle**). Parallel_magnetic_fluctuations and (**Bottom**) total temperature difference for_different plasma beta condition.[]{data-label="Fig8"}](Fig8-eps-converted-to.pdf){width="50.00000%"}

The_results_from the 2D simulations_with $\delta b_0=1$ and different_plasma beta are summarized in Fig. \[Fig8\]._Density fluctuations_are larger_for low beta plasma and become progressively smaller as the beta_increases. The amplitude of parallel magnetic_field fluctuations seems not_to be_affected_by the plasma_beta,_but the_time for them to exponentially grow is_affected by_the plasma beta.

The proton heating shows_different stages also in_these_non-monochromatic cases. There is an initial_strong temperature increase and a beam_formation, that correlates with the_enhancement_of_the parallel magnetic field fluctuations._A secular stage is followed where_the heating rate_is significantly reduced. The final temperature gained_by_the particles is similar to the_monochromatic_cases and so are the trends_with_the_plasma beta and pump wave_amplitude of $\Delta T$ (and $T_\bot$,_$T_\parallel$). Such trends are not affected by the plasma_beta, rather, they_are