Calculate radius gravitational Mass electric Mass and Dachus density for NegativiumL - the first particle (marked as 1) and PositiviumL - the second particle (marked as 2) + constraints - all equations
result=Reduce[{c==299792458,ce==8.987551792*10^9,mu0==1.256637061*10^-6,e==1.602176634*10^-19,g==6.6743*10^-11,qE1==-2/3,qE2==1/3,a1==(3/4)r1,a2==(3/4)r2,7.34456*10^-39==(2m2+2me2+m1+me1)/(1-be),(3(m1+me1))/(4Pir1^3)==(3(m2+me2))/(4Pir2^3),Abs[(8g^2m1^2Pi)/(c^4r1^4)*{((r1^2-(2gm1r1)/c^2+a1^2)/((r1^2+a1^2Cos[θ]^2)))-(1/2)*(((r1^2-(2gm1r1)/c^2+a1^2Cos[θ]^2)*((r1^2+a1^2)^2-a1^2(r1^2-(2gm1r1)/c^2+a1^2)Sin[θ]^2))/(r1^2+a1^2Cos[θ]^2)^3)}]>=Abs[-((36ce^5e^8PiqE1^4)/(c^10me1^4mu0r1^6))+(12ce^3e^4PiqE1^2)/(c^6me1^2mu0r1^4)],me1==(4Pir1^3)/(3c^2)*((ce^2e^2qE1*I)/(2mu0c^2r1^4)+(3ce^4e^6qE1^3*I^3)/(2mu0me1^2c^6r1^6)),r1==((2m1g)/c^2-Sqrt[(2m1g/c^2)^2-4qE1^2(25/16)*(3ce^2e^4)/(me1^2c^4)])/(50/16),me1==(4Pir1^3)/(3c^2)*((ce^2e^2qE1*I)/(2mu0c^2r1^4)+(3ce^4e^6qE1^3*I^3)/(2mu0me1^2c^6r1^6)),Abs[(8g^2m2^2Pi)/(c^4r2^4)*{((r2^2-(2gm2r2)/c^2+a2^2)/((r2^2+a2^2Cos[θ]^2)))-(1/2)*(((r2^2-(2gm2r2)/c^2+a2^2Cos[θ]^2)*((r2^2+a2^2)^2-a2^2(r2^2-(2gm2r2)/c^2+a2^2)Sin[θ]^2))/(r2^2+a2^2Cos[θ]^2)^3)}]>=Abs[-((36ce^5e^8PiqE2^4)/(c^10me2^4mu0r2^6))+(12ce^3e^4PiqE2^2)/(c^6me2^2mu0r2^4)],me2==(4Pir2^3)/(3c^2)*((ce^2e^2qE2*I)/(2mu0c^2r2^4)+(3ce^4e^6qE2^3*I^3)/(2mu0me2^2c^6r2^6)),r2==((2m2g)/c^2-Sqrt[(2m2g/c^2)^2-4qE2^2(25/16)*(3ce^2e^4)/(me2^2c^4)])/(50/16),(*Constraints*)0<be<1,0<θ<Pi,10^-30>m1>0,10^-30>m2>0,10^-30>Abs[me1]>0,10^-30>Abs[me2]>0,10^-30>Abs[r1]>0,10^-30>Abs[r2]>0},{m1,me1,r1,m2,me2,r2,be},Complexes]