I am fitting the following ARX(1,1)-GARCHX(1,1,1):
\begin{align*} y_t&=c+a_1y_{t-1}+\gamma_1x_t+\varepsilon_t\\ h_t&=\delta+\omega_1h_{t-1}+\theta_1\varepsilon_{t-1}^2+\pi_1x_{1,t} \end{align*}
Delta is negative, suggesting the conditional variance can theoretically be negative. Does anyone know how to get around this problem? My variables are already in logarithmic form, and I have tried standardising them to no avail. Any help would be much appreciated.