<p><span style="padding: 3px 10px; border-radius: 5px; color: #ffffff; font-weight: bold; display: inline-block; background-color: #ff0000;">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;External Email - Use Caution&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</span></p><p></p><div dir="ltr"><div class="gmail_default" style="font-family:verdana,sans-serif"><div id="gmail-:2bj" class="gmail-ii gmail-gt" style="font-size:0.875rem;direction:ltr;margin:8px 0px 0px;padding:0px;font-family:Roboto,RobotoDraft,Helvetica,Arial,sans-serif"><div id="gmail-:2dq" class="gmail-a3s gmail-aXjCH" style="overflow:hidden;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:small;line-height:1.5;font-family:Arial,Helvetica,sans-serif"><div dir="ltr"><div class="gmail_default" style="font-family:verdana,sans-serif"><div class="gmail_default">I&#39;ve found something unusual in the lead fields computed using MNE-Python that I can&#39;t explain. When arrays containing a large number of magnetometer sensors are used with the &#39;sample&#39; subject BEM (sample-5120-5120-5120-bem.fif) and the ico4 source space to create a forward model, there are some unusual strong spikes in the lead field. These spikes are not present when I use a spherical BEM model instead. </div><div class="gmail_default"><br></div><div class="gmail_default">To demonstrate, I created three different theoretical arrays of point magnetometers (coil_id: 2000) measuring roughly the normal component of the magnetic field with respect to the scalp surface. Each model contains an increasing number of total sensors (102 sensors, 1632 sensors, and 3672 sensors). </div><div class="gmail_default"><br></div><div class="gmail_default">Figure 1 shows the lead field for a particular dipole (#364 using ico4 source space) for each array using a spherical BEM model. A smooth dipolar response is shown for each that is roughly equivalent over the arrays. The peak response for each array doesn&#39;t change appreciably, but there is an improvement in detail of the lead field. All evoked topomap plots are fixed to the same color bar limits (400 fT).</div><div class="gmail_default"><br></div><div class="gmail_default">Figure 2 shows the same set of lead fields but with the &#39;sample&#39; subject BEM instead. It shows curiously strong spikes present in the two high density arrays. These spikes vary in their location and magnitude. The main dipolar response is still present, and appears to be fairly consistent. However, due to the spikes, the peak response for each array is vastly different. Again, evoked topomap plots use a fixed color bar limit at 400 fT.</div><div class="gmail_default"><br></div><div class="gmail_default">Are these lead field patterns in Figure 2 explainable? They certainly seem like erroneous values, but perhaps something about the BEM model creates a focusing effect for the magnetic field response when we sample at different spatial frequencies? I&#39;m happy to provide some sample forward models for inspection. </div><div class="gmail_default"><br></div><div class="gmail_default">Thanks,</div><div class="gmail_default">Taylor Williams</div></div></div></div></div></div></div>