Restringir una expresión en LMFIT a 0.

Estoy empezando a usar LMFIT en un problema magnetostático y quiero asegurar que el producto punto de dos vectores sea cero. He definido los parámetros (componentes del vector) n0_x, n0_y, n0_z, a0_x, a0_y y a0_z y les he dado valores iniciales. Luego defino el parámetro:

param.add(‘a0dotn0’, expr=’n0_xa0_x+n0_ya0_y+n0_z*a0_z’)

pero a0dotn0 toma el valor del producto punto de los vectores iniciales a0 y n0. ¿Cómo puedo asegurarme de que la restricción de que a0dotn0 siempre sea cero se cumpla?