Fig. 1. Physical system of nanodiamond cantilever coupled to NV center

Fig. 2. Time evolution of negativity between mechanical model a and NV center and state populations

Fig. 3. Time evolution of entanglement between mechanical mode a and NV center under ideal conditions (g=1 and κ1=κ2=γ=0). (a) θ2=π/4, θ3=π/2; (b) θ1=π/2, θ3=π/2; (c) θ1=π/2, θ2=π/4

Fig. 4. Time evolution of entanglement between mechanical mode a and NV center under condition that dissipation is considered (g=1, κ1=0.1g, κ2=0.1g, γ=0.1g). (a) θ2=π/4, θ3=π/2; (b) θ1=π/2, θ3=π/2; (c) θ1=π/2, θ2=π/4

Fig. 5. Time evolution of entanglement between mechanical mode a and NV center. (a) Influence of decay rate of mechanical mode on entanglement; (b) influence of spontaneous decay rate of NV center on entanglement

Fig. 6. Time evolution of entanglement between mechanical mode b and NV center (g=1, κ1=κ2=γ=0). (a) θ2=π/2, θ3=0; (b) θ1=π/4, θ3=0; (c) θ1=π/4, θ2=π/2

Fig. 7. Time evolution of entanglement between mechanical mode b and NV center (g=1, κ1=0.1g, κ2=0.1g, γ=0.1g). (a) θ2=π/2, θ3=0; (b) θ1=π/4, θ3=0; (c) θ1=π/4, θ2=π/2

Fig. 8. Time evolution of entanglement between mechanical mode b and NV center. (a) Influence of decay rate of mechanical mode on entanglement; (b) influence of spontaneous decay rate of NV center on entanglement

Fig. 9. Time evolution of negativity between mechanical modes a and b and state populations

Fig. 10. Time evolution of entanglement between mechanical modes a and b (g=1, κ1=κ2=γ=0). (a) θ2=0, θ3=π/4; (b) θ1=π/2, θ2=0

Fig. 11. Time evolution of entanglement between mechanical modes a and b (g=1, κ1=0.1g, κ2=0.1g, γ=0.1g). (a) θ2=0, θ3=π/4; (b) θ1=π/2, θ2=0

Fig. 12. Time evolution of entanglement between mechanical modes a and b. (a) Influence of decay rate of mechanical mode on entanglement; (b) influence of spontaneous decay rate of NV center on entanglement