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

Metaclass IV:

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

By Theorem 3.5 we may assume that a2 = a1 and b2 = b1

sub={a2->a1,b2->b1};
kappa=kappa//.sub;

Let us also assume that a1 != a3

Medium characteristics

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

Define coefficients in Theorem 5.1:

rho = (1/(8 (a3-a1) a4));
mu=-a1;
lambda=b1^2;
sigma=((((a4^2)-(((a3-a1))^2)))^2)+(b1^2) (2 (a4^2)+(b1^2)+2 (((a3-a1))^2));
Abivector=(({{0, 0, 0,  ((a3-a1)^2+(a4^2)+(b1^2)+(Sqrt[sigma])) }, {0, 0, 2 (a3-a1) a4 , 0}, {0, -2 (a3-a1) a4 , 0, 0}, {- ((a1-a3)^2+(a4^2)+(b1^2)+(Sqrt[sigma])) , 0, 0, 0}}));
Bbivector=(({{0, 0, 0, ((a1-a3)^2+(a4^2)+(b1^2)-(Sqrt[sigma]))}, {0, 0, 2 (-a1+a3) a4, 0}, {0, 2 (a1-a3) a4, 0, 0}, {-((a1-a3)^2+(a4^2)+(b1^2)-(Sqrt[sigma])), 0, 0, 0}}));
Simplify[Abivector+Transpose[Abivector]]
Simplify[Bbivector+Transpose[Bbivector]]

Verify claim:

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

Note: 
* kappa + mu Id is not 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/(-a1+a3)) Bbivector;
(* 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]]]