We report a lattice QCD calculation of the strange quark contribution to the nucleon's magnetic moment and charge radius. This analysis presents the first direct determination of strange electromagnetic form factors including at the physical pion mass. We perform a model-independent extraction of the strange magnetic moment and the strange charge radius from the electromagnetic form factors in the momentum transfer range of 0.051 GeV 2 < ∼ Q 2 < ∼ 1.31 GeV 2 . The finite lattice spacing and finite volume corrections are included in a global fit with 24 valence quark masses on four lattices with different lattice spacings, different volumes, and four sea quark masses including one at the physical pion mass. We obtain the strange magnetic moment G s M (0) = −0.064(14)(09) µN . The four-sigma precision in statistics is achieved partly due to low-mode averaging of the quark loop and low-mode substitution to improve the statistics of the nucleon propagator. We also obtain the strange charge radius r 2 s E = −0.0043(16)(14) fm 2 .The determination of the strange (s) quark contribution to nucleon electromagnetic (EM) form factors is of immense importance since this is a pure sea quark effect. A nonzero value of the strange Sachs electric form factor (FF) G s E at any Q 2 = 0 would mean that the spatial distributions of s ands quarks are not the same in the nucleon. Since the extraction of the vector strange matrix elements N |sγ µ s|N was proposed in Refs. [1-3] via parity-violating e − N scattering for which the dominant contribution arises from interference between photon (γ) and weak boson (Z) exchanges by the following relation assuming isospin symmetry:a considerable number of experimental efforts by the SAMPLE, HAPPEX, G0, and A4 [4-15] Collaborations have been going on for the past two decades. The world data constrains that G s M (0) contributes less than 6% and r 2 s E contributes less than 5% to the magnetic moment and the mean-square charge radius of the proton respectively [16]. However, all these experimental results are limited by rather sizable error bars. Three different global analyses give G s M (Q 2 = 0.1 (GeV/c) 2 ) consistent with zero within uncertainties and differ in sign in their central values [17][18][19].Despite tremendous theoretical efforts, e.g. [20][21][22][23], a detailed convincing understanding about the sign and magnitude of strange EM FFs is still lacking. A detailed review of these theoretical efforts can be found in [24].Since the direct calculation of the s-quark loop in the (8) µ N and, for the first time, a nonzero signal for G s E (Q 2 ) which gave r 2 s E = −0.0067(25) fm 2 . However, one still has to perform the calculation at the physical pion mass and on several lattices to consider volume and finite cutoff corrections and over all beat down the noise to obtain a convincing result which will substantially sharpen our picture of strange quark contributions to the nucleon's EM structure.Conventionally, we omit the unit nucleon magneton µ N for G s M in the rest of the lette...