In FESTIM verification, surface kinetic model gives the code like that,
F.SurfaceKinetics(
k_sb=k_sb,
k_bs=k_bs,
lambda_IS=lambda_IS,
n_surf=n_surf,
n_IS=n_IS,
J_vs=J_vs,
surfaces=1,
initial_condition=exact_solution_cs(t=0),
t=F.t,
),
but this setting is only available in festim1.4, but not found in the festim2 api.
why kinetic surface model is absent for festim new version?
Can alll boundary conditions be set as the plasma dirven permeation case?
Hi @williamhayes123 and welcome to the FESTIM discourse 
Indeed, this surface kinetics model is not implemented in FESTIM2. The reason for that is that FESTIM1 used RealElements in legacy-fenics, which until very recently were not available in FEniCSx.
However, with the (very) recent developments coming up to v2.2 (currently available as a pre-release v2.2rc.2), we introduced manifolds which will allow to simulate surface kinetics.
You can find the relevant tutorials there: Manifold subdomains — FESTIM tutorial
Depending on your actual problem, you may not need the full kinetic surface model if the rates are fast enough and you could get away with dissociation/recombination fluxes instead.
Happy to discuss more