begin:constant
intens = 1e28 # t-averaged intensity Wm^-2 (I_max for CP, not LP)
n0 = 1e29 # particle density m^-3
lambda = 1.0e-6 # laser wavelength
omega = 2.0 * pi * c / lambda # laser frequency
t_fwhm = 25 * femto # Pulse length FWHM (Intensity)
t_w = t_fwhm / sqrt(2*loge(2)) # 1/e value for pulse
t_Imax = 1.5 * t_w # time from simulation start that intensity reaches it's maximum
# convert to E-field magnitude - varies for LP/CP
# linear E_field
ampl = sqrt(2 * intens / 0.00265442)
# circular E_field
#ampl = sqrt(intens / 0.00265442)
end:constant
begin:control
# number cells
nx = 1024 / 2
# total number of macroparticles to load
npart = nx*256 / 2
# max number of iterations - set to -1 to run until finished
nsteps = -1
# final time of simulation - longer sims will be harder to finish because of pair/photons
t_end = 40e-15
# size of domain
x_min = -6.0e-6
x_max = 0
dt_multiplier = 0.95
smooth_currents = T # to help improve self-heating - just delete if you don't want it
check_stop_frequency = 1 # for stopping - outputs a restart dump (check manual)
simplify_deck = T
# restart_snapshot =
dlb_threshold = 0.01
end:control
begin:boundaries
bc_x_min = simple_laser
bc_x_max = reflect
end:boundaries
begin:collisions
use_collisions = T
coulomb_log = auto
collide = all
end:collisions
begin:qed
use_qed = T
qed_start_time = 0
produce_photons = T
photon_energy_min = 1 * mev
produce_pairs = T
photon_dynamics = T
end:qed
begin:species
name = Electron
charge = -1
mass = 1
frac = 0.5
dump = T
rho = if (abs(x) lt 0.5e-6, n0, 0)
identify:bw_electron
end:species
begin:species
name = Ion
charge = 1
mass = 1836.0
frac = 0.5
dump = T
temp = 0
rho = rho(Electron)
end:species
begin:species
name = Photon
npart = 0
dump = T
identify:photon
end:species
begin:species
name = Positron
charge = 1
npart = 0
dump = T
identify:bw_positron
end:species
begin:output
dt_snapshot = 1e-15
full_dump_every = 5
force_final_to_be_restartable = T
particles = full + single
px = full + single
py = full + single
pz = full + single
grid = always
ex = always
ey = always
ez = always
bx = always
by = always
bz = always
ekbar = always + species
ekflux = always + species
number_density = always + species
poynt_flux = always
absorption = always
total_energy_sum = always
end:output
begin:laser
boundary = x_min
amp = ampl
omega = omega
pol_angle = 0.0
phase = 0.0
t_profile = supergauss(time,t_Imax,t_w,6) # temporal supergaussian, e^(-(x/w)^n), n integer
end:laser