////////////////////////////////////////////////////////////////////////////////////////// // // input for lpic // ////////////////////////////////////////////////////////////////////////////////////////// &pulse_front ------------------------------------------------------------------------------------------ Q = 1 # switch ON (Q=1) OFF (Q=0) amplitude = 1 # dimensionless laser field amplitude amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic phase2 = 0 # 2nd harmonic's phase with respect to fundamental [degree] phase3 = 0 # 3rd harmonic's phase [degree] angle = 45 # in degree polarization = 2 # s=1, p=2, c=3 shape = 3 # linear=1, sin=2, sin^2=3, Gauss=4, sech=5 raise = 15 # pulse raise/fall time in periods duration = 30 # pulse duration in periods delay = 0 # propagation delay in cycles pulse_save = 1 # save pulse shape? yes=1, no=0 pulse_save_step = 0.004 # time step in periods &pulse_rear ------------------------------------------------------------------------------------------ Q = 0 # switch ON (Q=1) OFF (Q=0) amplitude = 0.0 # dimensionless laser field amplitude amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic phase2 = 0 # 2nd harmonic's phase with respect to fundamental phase3 = 0 # 3rd harmonic's phase angle = 0 # in degree polarization = 1 # s=1, p=2, c=3 shape = 3 # linear=1, sin=2, sin^2=3 raise = 15 # pulse raise/fall time in periods duration = 30 # pulse duration in periods delay = 0 # propagation delay in cycles pulse_save = 1 # save pulse shape? yes=1, no=0 pulse_save_step = 0.02 # time step in periods &propagate ------------------------------------------------------------------------------------------ prop_start = 0 # start time in periods prop_stop = 40 # stop time in periods &box ------------------------------------------------------------------------------------------ cells_per_wl = 1000 # cells per wavelength, lab frame cells = 8000 # total number of cells cells_left = 3000 # cells vacuum left cells_plasma = 2000 # occupied cells for cells_ramp=0 cells_ramp = 200 # cells in the linear ramp region cells_ramp_form = 1 # 0 = linear; 1 = exponential n_ion_over_nc = 400 # maximum density/critical density box_save = 1 # save configuration? yes=1, no=0 &electrons ------------------------------------------------------------------------------------------ fix = 0 # 0->not fixed 1->fixed ppc = 200 # max. number of MacroParticles per cell vtherm = 0.0001 # thermal velocity in units of C &ions ------------------------------------------------------------------------------------------ fix = 1 # 0->not fixed 1->fixed z = 1 # q/e m = 50000 # m/m_e ppc = 200 # max. number of MacroParticles per cell vtherm = 0.00 # thermal velocity in units of C &output ------------------------------------------------------------------------------------------ path = ../data_eff5_0b # output path &energy Q = 1 # energy plot? t_start = 0 # start time in periods t_stop = 40 # stop time in periods t_step = 1 # time step in periods &flux Q = 1 # flux plot? t_start = 0 # start time in periods t_stop = 40 # stop time in periods t_step = 1 # time step in periods &reflex Q = 1 # reflectivity plot? t_start = 0 # start time in periods t_stop = 40 # stop time in periods t_step = 1 # time step in periods &snapshot Q = 1 # snapshots? t_start = 0 # start time in periods t_stop = 40 # stop time in periods t_step = 1 # time step in periods &el_phasespace Q = 0 # phasespace plots? t_start = 0 # start time in periods t_stop = 40 # stop time in periods t_step = 0.05 # time step in periods &ion_phasespace Q = 0 # phasespace plots? t_start = 0 # start time in periods t_stop = 40 # stop time in periods t_step = 0.05 # time step in periods &el_velocity Q = 1 # electron velocity distributions? t_start = 0 # start time in periods t_stop = 40 # stop time in periods t_step = 1 # time step in periods &ion_velocity Q = 1 # ion velocity distributions? t_start = 0 # start time in periods t_stop = 40 # stop time in periods t_step = 1 # time step in periods &de Q = 1 # electron density plots? t_start = 0 # start time in periods t_stop = 40 # stop time in periods x_start = 0 # left boundary in cells x_stop = 100000 # right boundary in cells &di Q = 1 # ion density plots? t_start = 0 # start time in periods t_stop = 40 # stop time in periods x_start = 0 # left boundary in cells x_stop = 100000 # right boundary in cells &jx Q = 0 # jx plots? t_start = 0 # start time in periods t_stop = 40 # stop time in periods x_start = 0 # left boundary in cells x_stop = 100000 # right boundary in cells &jy Q = 0 # jy plots? t_start = 0 # start time in periods t_stop = 40 # stop time in periods x_start = 0 # left boundary in cells x_stop = 100000 # right boundary in cells &jz Q = 0 # jz plots? t_start = 0 # start time in periods t_stop = 40 # stop time in periods x_start = 0 # left boundary in cells x_stop = 100000 # right boundary in cells &ex Q = 0 # ex plots? t_start = 0 # start time in periods t_stop = 40 # stop time in periods x_start = 0 # left boundary in cells x_stop = 100000 # right boundary in cells &ey Q = 0 # ey plots? t_start = 0 # start time in periods t_stop = 40 # stop time in periods x_start = 0 # left boundary in cells x_stop = 100000 # right boundary in cells &ez Q = 0 # ez plots? t_start = 0 # start time in periods t_stop = 40 # stop time in periods x_start = 0 # left boundary in cells x_stop = 100000 # right boundary in cells &bx Q = 0 # bx plots? t_start = 0 # start time in periods t_stop = 40 # stop time in periods x_start = 0 # left boundary in cells x_stop = 100000 # right boundary in cells &by Q = 0 # by plots? t_start = 0 # start time in periods t_stop = 40 # stop time in periods x_start = 0 # left boundary in cells x_stop = 100000 # right boundary in cells &bz Q = 0 # bz plots? t_start = 0 # start time in periods t_stop = 40 # stop time in periods x_start = 0 # left boundary in cells x_stop = 100000 # right boundary in cells &edens Q = 1 # field energy density plots? t_start = 0 # start time in periods t_stop = 40 # stop time in periods x_start = 0 # left boundary in cells x_stop = 100000 # right boundary in cells &traces Q = 1 # traces? t_start = 0 # start time in periods t_stop = 40 # stop time in periods traces = 7 # of traces at fixed positions x: t0=2, t1=2500, t2=3000, t3=3500, t4=4000, t5=4500, t6=7998 ////////////////////////////////////////////////////////////////////////////////////////// &restart ------------------------------------------------------------------------------------------ Q = 0 # restart from intermediate stage? file = restart # start file Q_save = 1 # save intermediate stages periodically? file_save = restart # save file ¶llel ------------------------------------------------------------------------------------------ N_domains = 1 # number of parallel processes Q_reo = 1 # periodic reorganizations? delta_reo = 1 # laser cycles between reo's ////////////////////////////////////////////////////////////////////////////////////////// thermal velocity v/c vs. thermal energy [eV] M ( vx^2 + vy^2 + vz^2 ) = 3 k T ------------------------------------------------------------------------------------------ 1e-4 0.005 3e-4 0.046 5e-4 0.128 1e-3 0.511 4e-3 8.2 1e-2 51 2e-2 204 3e-2 460 4e-2 818 4.42e-2 1000 1e-1 5110 choose: v_th/c = (lambda_d/dx) * (omega_p/omega) * (2pi dx/lambda) :=1 = sqrt(z_ion * n_ion_over_nc) = 2pi/cells_per_wl //////////////////////////////////////////////////////////////////////////////////////////