در این پژوهش، از یک مدل مونت کارلو دو دره ای برای شبیه سازی اتصال شاتکی استفاده شده است. در مدل ارائه شده، پراکندگی های ناخالصی و فونونی درنظر گرفته شده اند و یونیزاسیون برخوردی در ماتریس پراکندگی قرار گرفته است. باندهای انرژی غیرسهموی فرض شده اند و سازوکارهای تونل زنی و گسیل گرمایونی مولفه های جریانی هستند. با اضافه کردن یک لایه نازک غیر هم نوع با نیم رسانا نشان داده شده است که شکل گیری یک میدان الکتریکی مخالف با جهت حرکت الکترون در مرز اتصال باعث افزایش ارتفاع موثر سد شاتکی می شود. با تغییر چگالی غلظت ناخالصی این لایه نازک تغییر ارتفاع موثر سد شاتکی و در نتیجه جریان گذرنده شبیه سازی شده است. مقایسه نتایج بدست از شبیه سازی با داده-های معتبر علمی، درستی کارکرد مدل ارائه شده را تایید می کند.