Dipoles switching

Hi Everyone,
I wrote this code to get dipoles switching after applying an external electric field, but I’m getting errors on lost atoms after specific time steps. Sometimes at 11 ps and sometimes at 17 ps.
Can you please take a look at my code and tell me what am I doing wrong?
I have changed electric field magnitude and directions, but I’m getting same errors at the same time steps.
My other question is that how can I define specific temperature and keep it fixed during the simulation for sphere atome style?

units metal
dimension 3
boundary p p p
atom_style hybrid dipole sphere
atom_modify map array

neighbor 2.0 bin
neigh_modify delay 10 check yes

read_data BiFeO3.dat

pair_style hybrid/overlay coul/streitz 12.0 wolf 0.20 meam/c lj/cut/dipole/cut 10.0
pair_coeff * * meam/c BiFeO3zigzag.library.meam Fe Bi O BiFeO3zigzag.parameter.meam Fe Bi O
pair_coeff * * coul/streitz FeBiO3.streitz Fe Bi O
pair_coeff * * lj/cut/dipole/cut 1.0 1.0

fix 1 all spring/self 1000.0
fix_modify 1 energy yes
fix 3 all qeq/slater 1 12.0 1.0e-6 100 coul/streitz
dump 1 all custom 10000 dump.atom x y z mux muy muz
dump 2 all xyz 10000 dump.xyz
timestep 0.0001
run 10000
fix 5 all efield -0.0003 0.0 -0.0003
fix 6 all nve/sphere update dipole/dlm
compute 1 all property/atom mux muy muz
run 200000
unfix 5
run 200000

Data File:
1 3 1.12238351 1.13820203 1.67962192 -0.81466667 1 1 1 1.2 1.141
2 3 1.12238351 1.13820203 15.7995499 -0.81466667 1 1 1 1.2 1.141
3 3 1.12238351 1.13820203 29.9194779 -0.81466667 1 1 1 1.2 1.141
4 3 1.12238351 1.13820203 44.0394059 -0.81466667 1 1 1 1.2 1.141
5 3 1.12238351 1.13820203 58.1593339 -0.81466667 1 1 1 1.2 1.141