SetDirectory["~/KappaLib"];
<<kappaLib-1.2.m

Metaclass I:

kappa=emMatrixToKappa[(({{a1, 0, 0, -b1, 0, 0}, {0, a2, 0, 0, -b2, 0}, {0, 0, a3, 0, 0, -b3}, {b1, 0, 0, a1, 0, 0}, {0, b2, 0, 0, a2, 0}, {0, 0, b3, 0, 0, a3}}))];

By Theorem 3.5 we may assume that a2 = a3 and b2 = b3

sub={a3->a2,b3->b2}

In addition, let us assume that a1 = a2

sub=Append[sub,a2->a1]

kappa=kappa//.sub;

Medium characteristics

emKappaToMatrix[kappa]//MatrixForm
FullSimplify[emDet[kappa]]
Simplify[emTrace[kappa]]

Define coefficients in Theorem 5.1:

Note: Since a1=a2, we may assume that b1 != b2 (by representation theorem 
of medium with double light cone)

rho =((b2^2)-(b1^2))/2;
mu=-a1;
lambda=b2^2;
Abivector=(({{0, 1, 0, 0}, {-1, 0, 0, 0}, {0, 0, 0, 0}, {0, 0, 0, 0}}));
Bbivector=(({{0, 0, 0, 0}, {0, 0, 0, 0}, {0, 0, 0, 1}, {0, 0, -1, 0}}));

Verify claim:

eta=kappa+mu emIdentityKappa[];
LHS =emCompose[eta,eta];
RHS= -lambda emIdentityKappa[]+ emBiProduct[rho, Abivector, Bbivector]+ emBiProduct[rho, Bbivector, Abivector];
Union[Flatten[LHS-RHS]]

Note: 
* kappa + mu Id is trace-free in this case
* lambda > 0
* rho, A, B are all non-zero

emTrace[eta]

Solvability of equation for D

Define constants alpha, beta, gamma that appear in definition of algebraically decomposable medium

alpha = lambda + mu^2;
beta = mu;
gamma = 1;

Explicitly verify that kappa is algebraically decomposable:

LHS =alpha   emIdentityKappa[] +beta (kappa +emPoincare[kappa]) + gamma emCompose[emPoincare[kappa],kappa]; 
RHS=emBiProduct[rho,Abivector,Bbivector]+emBiProduct[rho,Bbivector,Abivector];
eqs=Union[Simplify[Flatten[(LHS-RHS)]]]

Existence of D:

Dbivector=-1/b1 Abivector;

(* contract kappa with a bivector from the left *)
contract[biv_,kappa_]:=Table[
1/2 Sum[
biv[[i]][[j]]emReadNormal[kappa,a,b,i,j]
,
{i,1,4},{j,1,4}
]
,
{a,1,4},{b,1,4}
]
dLHS = contract[Dbivector,kappa+beta emIdentityKappa[]];
dRHS = 1/2 emTrace[emBiProduct[rho,Dbivector,Dbivector]]Abivector+Bbivector;
deqs=Union[Flatten[Simplify[dLHS-dRHS]]]