In this paper, the direct similarity reduction method is used to reduce the modified variable coefficient Zakharove–Kusnetsov equation (vcmZK), which represents the propagation of quantum dust acoustic waves (QDAWs) in an inhomogeneous magnetized dusty plasma field, into a nonlinear ordinary differential equation (NODE). Moreover, the Jacobi elliptic method is used to solve the NODE, where many novel double periodic, kink, and solitary wave solutions are obtained. Additionally, we have chosen to plot periodic, kink, and soliton QDAWs solutions according to some choices for the variable coefficients. According to the a chieved results, we can conclude that the choices of the functional forms of the variable coefficients can control the QDAWs propagation. Finally, the new discovered results may have many important applications in astrophysics, atmospheric science, nanotechnology, and industrial plasma physics.