q=κ0(zσ′′-θσ+2κ0κz)/(σ)2,