With b belonging to a new $BMO_{\theta}(\rho )$ space, $L=-\triangle +V$ is a Schrödinger operator on ${\mathbb{R}^{n}}$ with nonnegative potential V belonging to the reverse Hölder class $RH_{n/2}$ . The fractional integral operator associated with L is denoted by ${\mathcal{I}}_{\beta}^{L}$ . We investigate the boundedness of ${\mathcal{I}}_{\beta}^{L}$ and $[b,{\mathcal{I}}_{\beta}^{L}]$ , which are its commutators with $b_{\theta}(\rho )$ on vanishing generalized mixed Morrey spaces $VM_{\vec{p},\varphi}^{\alpha ,V}$ related to Schrödinger operation and generalized mixed Morrey spaces $M_{\vec{p},\varphi}^{\alpha ,V}$ . The boundedness of the operator ${\mathcal{I}}_{\beta}^{L}$ is ensured by finding sufficient conditions on the pair $(\varphi _{1},\varphi _{2})$ , which goes from $M_{\vec{p},\varphi _{1}}^{\alpha ,V}$ to $M_{\vec{q},\varphi _{2}}^{\alpha ,V}$ , and from $VM_{\vec{p},\varphi _{1}}^{\alpha ,V}$ to $VM_{\vec{q},\varphi _{2}}^{\alpha ,V}$ , $\sum \limits _{i=1}^{n}\frac{1}{p_{i}}-\sum \limits _{i=1}^{n}\frac{1}{q_{i}}=\beta $ . When b belongs to $BMO_{\theta}(\rho )$ and $(\varphi _{1},\varphi _{2})$ satisfies some conditions, we also show that the commutator operator $[b,{\mathcal{I}}_{\beta}^{L}]$ is bounded from $M_{\vec{p},\varphi _{1}}^{\alpha ,V}$ to $M_{\vec{q},\varphi _{2}}^{\alpha ,V}$ and from $VM_{\vec{p},\varphi _{1}}^{\alpha ,V}$ to $VM_{\vec{q},\varphi _{2}}^{\alpha ,V}$ .