Since we are interested in finding the unit cost let the output needed to be produced be equal to 1 i.e., $y(i)\overset{set}=1$
We need to solve: $$\begin{align}
\min_{l(i),x(i),h(i)\geq0} \quad & w_Ll(i)+x(i)+w_{H}h(i)\\
\textrm{s.t.} \quad & [l(i)^{\frac{\epsilon - 1}{\epsilon}} +\alpha(i)(\tilde{\gamma}x(i))^{\frac{\epsilon - 1}{\epsilon}}]^{\frac{\epsilon \beta}{\epsilon-1}}h(i)^{1-\beta}=1\\
\textrm{and} \quad & \alpha(i)=\begin{cases}1 & \text{if } i\text{ is automated}\\ 0 & \text{otherwise}\end{cases}
\end{align}$$
I am going to use the indicator function and solve for two different cases and then try to combine them using the indicator.
Case 1: If $i$ is not automated i.e., $ \alpha(i)=0$
If $i$ is not automated, then the cost minimization problem becomes a standard one with a Cobb-Douglas production function
$$\begin{align}
\min_{l,x,h\geq 0} \quad & w_Ll+x+w_{H}h\\
\textrm{s.t.} \quad & l^{\beta}h^{1-\beta}=1\\
\end{align}$$
This is just a standard cost minimization problem with a cobb-douglas production function, so I think you can verify yourself that the above problem gives: $l(i)=\left(\frac{w_H\beta}{w_L(1-\beta)}\right)^{1-\beta}, \quad x(i)=0, \quad h(i)= \left(\frac{w_L(1-\beta)}{w_{H} \beta}\right)^{\beta}$
therefore, $$\begin{eqnarray}
& c=w_L^\beta\left(\frac{w_H\beta}{1-\beta}\right)^{1-\beta}+w_H^{1-\beta}\left(\frac{w_L(1-\beta)}{\beta}\right)^{\beta}\\\\
\implies & {\boxed{c=w_L^\beta w_H^{1-\beta} \beta^{-\beta}(1-\beta)^{\beta-1}}}\tag{a}
\end{eqnarray}$$
Case 2: If $i$ is automated
$$\begin{align}
\min_{l,x,h\geq 0} \quad & w_Ll+x+w_{H}h\\
\textrm{s.t.} \quad & [l^{\frac{\epsilon - 1}{\epsilon}} +(\tilde{\gamma}x)^{\frac{\epsilon - 1}{\epsilon}}]^{\frac{\epsilon \beta}{\epsilon-1}}h^{1-\beta}=1\\
\end{align}$$
let us write the above problem as:
$$\begin{align}
{\min_{0\leq h \leq 1}
\underbrace{\begin{pmatrix} \underset{l,x}{\min} \quad & w_Ll+w_Hh+x \\
\textrm{s.t.} \quad & [l^{\frac{\epsilon - 1}{\epsilon}} +(\tilde{\gamma}x)^{\frac{\epsilon - 1}{\epsilon}}]^{\frac{\epsilon}{\epsilon-1}}=h^\frac{\beta -1}{\beta} \end{pmatrix}}_{\text{auxiliary problem}}} \tag{1}
\end{align}$$
We can solve the above problem in two stages:
Stage 1: first solve the auxiliary problem for $l$ and $x$ holding $h$ fixed
Stage 2: Substitute the solutions of the auxiliary problem$-$ $l(h)$ and $x(h)$ $-$ into $(1)$ and solve the problem to find optimal $h$
Stage 1:
we need to solve: $$\begin{align}
\underset{l,x}{\min} \quad & w_Ll+w_Hh+x \\
\textrm{s.t.} \quad & [l^{\frac{\epsilon - 1}{\epsilon}} +(\tilde{\gamma}x)^{\frac{\epsilon - 1}{\epsilon}}]^{\frac{\epsilon}{\epsilon-1}}=h^\frac{\beta -1}{\beta}\\
\text{gives:} \quad & {\left.\begin{matrix} l(h)=\frac{h^\frac{\beta -1}{\beta}}{[1+(w_L\tilde \gamma)^{\epsilon -1}]^\frac{\epsilon}{\epsilon-1}} \\ x(h)=\frac{w_L^\epsilon \tilde{\gamma}^{\epsilon -1} h^\frac{\beta -1}{\beta}}{[1+(w_L\tilde \gamma)^{\epsilon -1}]^\frac{\epsilon}{\epsilon-1}}
\end{matrix}\right\}}\tag{2}
\end{align}$$
Because we are given that $\epsilon>1$, You can verify that isoquants will be convex. Thus, I was able to solve the above by equating ratios of MP to MC for $l$ and $x$ and substituting them into the constraint.
Stage 2: substituting $(2)$ in $(1)$ we get the final problem solving which will give optimal $h$
$$\begin{align}
\min_{0\leq h\leq 1} \quad & w_Hh+w_L[1+(w_L\tilde \gamma)^{\epsilon -1}]^\frac{1}{1-\epsilon} h^\frac{\beta -1}{\beta}\tag{3}
\end{align}$$
the above objective is convex in $h \; \because \beta<1$, so assuming the parametric condition gives us a stationary point that lies in the interval $(0,1)$ we can solve for the optimal h in the above problem and substitute it in the problem itself to get $c$ for $\alpha=1$:
$$h^*=\left(\frac{(1-\beta)w_L[1+(w_L\tilde \gamma)^{\epsilon -1}]^\frac{1}{1-\epsilon}}{w_H\beta}\right)^\beta$$
substituting $h$ in $(3)$ give us:
$$\boxed{c=w_H^{1-\beta}(1-\beta)^{\beta-1}\beta^{-\beta}\left[w_L^{1-\epsilon}+\tilde \gamma^{\epsilon-1}\right]^\frac{\beta}{1-\epsilon}}\tag{b}$$
we can combine the solutions to case 1 and case 2 , i.e., $(a)$ and $(b)$
Therefore, we have
$c=\begin{cases} w_H^{1-\beta} \beta^{-\beta}(1-\beta)^{\beta-1}w_L^\beta & \text{if } \alpha(i)=0 \\
w_H^{1-\beta}(1-\beta)^{\beta-1}\beta^{-\beta}\left[w_L^{1-\epsilon}+\tilde \gamma^{\epsilon-1}\right]^\frac{\beta}{1-\epsilon} & \text{if } \alpha(i)=1
\end{cases}$
using $\alpha(i)$, the above can also be written as:
$$\boxed{c=w_H^{1-\beta}\beta^{-\beta}(1-\beta)^{\beta-1}\left(w_L^{1-\epsilon}+\tilde \gamma^{\epsilon-1} \alpha(i) \right)^\frac{\beta}{1-\epsilon}}$$