  Motivated by some recent potential theoretic results on subordinate killed
L\'evy processes in open subsets of the Euclidean space, we study processes in
an open set $D\subset {\mathbb R}^d$ defined via Dirichlet forms with jump
kernels of the form $J^D(x,y)=j(|x-y|)\mathcal{B}(x,y)$ and critical killing
functions. Here $j(|x-y|)$ is the L\'evy density of an isotropic stable process
(or more generally, a pure jump isotropic unimodal L\'evy process) in
$\mathbb{R}^d$. The main novelty is that the term $\mathcal{B}(x,y)$ tends to 0
when $x$ or $y$ approach the boundary of $D$. Under some general assumptions on
$\mathcal{B}(x,y)$, we construct the corresponding process and prove that
non-negative harmonic functions of the process satisfy the Harnack inequality
and Carleson's estimate. We give several examples of boundary terms satisfying
those assumptions. The examples depend on four parameters, $\beta_1, \beta_2,
\beta_3$, $\beta_4$, roughly governing the decay of the boundary term near the
boundary of $D$.
  In the second part of this paper, we specialise to the case of the half-space
$D=\mathbb{R}_+^d=\{x=(\widetilde{x},x_d):\, x_d>0\}$, the $\alpha$-stable
kernel $j(|x-y|)=|x-y|^{-d-\alpha}$ and the killing function$\kappa(x)=c
x_d^{-\alpha}$, $\alpha\in (0,2)$, where $c$ is a positive constant. Our main
result in this part is a boundary Harnack principle which says that, for any
$p>(\alpha-1)_+$, there are values of the parameters $\beta_1, \beta_2,
\beta_3$, $\beta_4$, and the constant $c$ such that non-negative harmonic
functions of the process must decay at the rate $x_d^p$ if they vanish near a
portion of the boundary. We further show that there are values of the
parameters $\beta_1, \beta_2, \beta_3$, $\beta_4$, for which the boundary
Harnack principle fails despite the fact that Carleson's estimate is valid.
