The paper presents a decoupled, conservative virtual element method for inductionless magnetohydrodynamic equations. A semi-implicit Euler formulation is employed for time derivative terms, while the mixed virtual element method is used for spatial discretization. The approach utilizes the \(\textbf{H}^1\) -conforming Stokes-like virtual element and discontinuous piecewise polynomials to approximate velocity and pressure, \(\textbf{H}(\text{ div})\) -conforming virtual element and discontinuous piecewise polynomials for the current density and electric potential, respectively. In the fully discretized scheme, the velocity and pressure are decoupled from the current density and electric potential. Additionally, the discrete velocity and current density are pointwise divergence-free. Optimal error estimates are established for the velocity, current density, and electric potential. Finally, numerical experiments are conducted to validate the effectiveness of the theoretical analysis.
