MCPlas, düşük sıcaklıklı ve termal olmayan plazmalar için akışkan-Poisson modellerinin MATLAB üzerinden otomatik olarak oluşturulup COMSOL Multiphysics içinde çalıştırılmasını sağlayan açık kaynaklı bir araç kutusudur. Modeli tanımlayan parçacık türleri, reaksiyon kinetiği, taşınım ve hız katsayıları, yüzey koşulları ve sayısal ayarlar yapılandırılmış JSON belgelerinde tutulur; MATLAB kodu bu girdileri okuyarak LiveLink for MATLAB aracılığıyla denklem tabanlı COMSOL modelini otomatik oluşturur. Çalışmadaki DC ve RF argon parıltı boşalması testlerinde, MCPlas ile COMSOL Plasma Module aynı fiziksel denklemler ve sınır koşullarıyla yapılandırıldığında sonuçlar son derece yakın biçimde çakışmıştır. Buna karşılık MCPlas’a özgü gelişmiş elektron taşınım yaklaşımı ve sınır koşulları etkinleştirildiğinde özellikle elektrot yakınındaki katot-kılıf bölgelerinde belirgin farklılıklar ortaya çıkmıştır. Sonuçlar bir yazılım/model doğrulama ve yeniden üretilebilirlik çalışmasına dayanmaktadır; deneysel bir plazma ölçümü veya bütün plazma koşullarına genellenebilecek evrensel bir doğrulama sunmamaktadır.

Çalışmanın ikinci önemli sonucu, reaksiyon kinetiğinin yazılıma elle tekrar girilmesi yerine şema doğrulamalı JSON biçiminde taşınabilmesidir. MCPlas, dört parçacık türünden oluşan basit argon modeli ile 23 parçacık türü ve 409 çarpışma/radyatif süreç içeren ayrıntılı reaksiyon kinetiği modeli arasında giriş dosyaları üzerinden geçiş yapabilmektedir. Aynı JSON tabanlı reaksiyon kinetiği verileri ayrıca PLASIMO ve FEDM araçlarında kullanılmış ve üç modelleme sistemi arasında yalnızca birkaç yüzde düzeyindeki en büyük farklılıklarla yüksek uyum elde edilmiştir. Böylece çalışma, MATLAB’ın yalnızca hesaplama dili olarak değil; model tanımı, veri standardizasyonu, COMSOL otomasyonu ve farklı simülasyon ortamları arasında yeniden üretilebilirlik sağlayan bir orkestrasyon katmanı olarak kullanılabileceğini göstermektedir.

Araştırmanın temel problemi nedir?

Düşük sıcaklıklı plazmaların akışkan modelleri; elektronların, iyonların ve nötr parçacıkların yoğunluklarını, elektron enerjisini, elektrik alanı ve çok sayıda kimyasal/radyatif reaksiyonu birlikte çözmeyi gerektirir. Model karmaşıklaştıkça yalnız diferansiyel denklemleri çözmek değil, onlarca parçacık türüne ait fiziksel özellikleri ve yüzlerce reaksiyonu doğru biçimde yazılıma aktarmak da önemli bir sorun hâline gelir.

Çalışmanın vurguladığı temel sorun, farklı modelleme araçlarında reaksiyon kinetik modellerinin çoğunlukla ayrı ayrı ve elle tanımlanmasıdır. Bu yaklaşım veri giriş hatalarına, modeller arasında küçük fakat bilimsel olarak önemli farklılıklara ve yayımlanmış bir simülasyonun başka bir araştırma grubu tarafından tam olarak yeniden kurulmasının zorlaşmasına yol açabilir.

MCPlas bu problemi iki katmanda ele almaktadır. Birinci katman, fiziksel modeli MATLAB fonksiyonlarıyla açık ve düzenlenebilir biçimde tanımlamaktır. İkinci katman ise modele ait bütün temel giriş bilgisini makine tarafından okunabilen ve şemayla doğrulanabilen JSON belgelerine taşımaktır. Böylece fiziksel modelin tanımı COMSOL proje dosyasının içine gizlenmiş tekil ayarlardan mümkün olduğunca ayrılmaktadır.

MATLAB bu sistemde tam olarak ne yapıyor?

MATLAB, MCPlas mimarisinde COMSOL için yalnızca yardımcı bir hesap makinesi değildir. Araç kutusunun ana kontrol ve model üretim katmanıdır. Kullanıcı genel model girdilerini ve reaksiyon kinetiği verilerini JSON biçiminde sağlar. Ana MCPlas.m betiği bu dosyaları okur, ilgili değişkenleri ve fiziksel katsayıları tanımlar, denklemleri oluşturur, sınır koşullarını ekler, geometri ve ağı kurar, sayısal çözüm ayarlarını belirler ve LiveLink for MATLAB üzerinden COMSOL modelini derler.

Çalışmanın Şekil 2’sindeki akış dört temel aşamadan oluşmaktadır: giriş verilerinin hazırlanması, MATLAB ile modelin kurulması, LiveLink üzerinden COMSOL modelinin derlenmesi ve kullanıma hazır model dosyasının oluşturulması. Bu iki orta aşama MCPlas.m çalıştırıldıktan sonra otomatik yürütülür.

MCPlas hangi MATLAB fonksiyonlarından oluşuyor?

Çalışmanın Tablo 1’i, model oluşturma zincirinin parçalarını açık biçimde göstermektedir. Fonksiyon isimlerinin fiziksel modelin hangi aşamasını oluşturduğu doğrudan okunabildiği için bu yapı yeniden üretilebilirlik açısından önemlidir.

Bu liste, modelin hangi fiziksel bileşeninin hangi kod parçası tarafından oluşturulduğunu görünür hâle getirir. Bu nedenle MCPlas’ın bilimsel önemi yalnız otomasyon hızından değil, denklem ve sınır koşullarının doğrudan MATLAB kaynak kodunda izlenebilmesinden gelir.

Desteklenen model geometrileri nelerdir?

Şekil 1, araç kutusunun dört geometrik modelleme seçeneğini gösterir. 1D seçeneği Kartezyen, 1p5D seçeneği polar koordinatlı tek boyutlu modeli; 2D Kartezyen ve 2p5D ise silindirik koordinatlı iki boyutlu modeli temsil eder. Geometriler düzlemsel-paralel ve koaksiyel plazma kaynaklarına uygun olacak şekilde düzenlenmiştir. Elektrot yüzeylerinde isteğe bağlı dielektrik katmanlar da modele dahil edilebilir.

JSON neden bilimsel olarak önemli?

MCPlas’ın giriş verileri iki ana gruptadır: reaksiyon kinetiği modeli ve genel model kurulumu. Reaksiyon kinetiği tarafında üst düzey JSON yapısı üç ana bileşen içerir: references, states ve processes. References alanı kullanılan veri kaynaklarını; states plazmadaki parçacık ve enerji durumlarını; processes ise çarpışma ve radyatif reaksiyonları tanımlar.

Çalışmada bu yapı yeni nesil LXCat veri modeliyle ilişkilendirilmiş ve JSON şemalarıyla doğrulanabilir hâle getirilmiştir. Genel model girdileri ise plazma kaynağı, plazma ortamı ve tanısal/modelleme özelliklerini içerir. Bu bilgiler için Plasma-MDS tabanlı şema ve şema tabanlı metadata oluşturma yaklaşımı kullanılmaktadır.

Bilimsel açıdan önemli nokta şudur: reaksiyon kinetiği, COMSOL arayüzündeki elle doldurulmuş bir dizi ayar olmaktan çıkarak bağımsız, makine tarafından okunabilir ve başka kodlara aktarılabilir bir veri nesnesine dönüşmektedir.

Akışkan-Poisson modeli nasıl kuruluyor?

MCPlas’ın temel fiziksel modeli, her parçacık türü için sayı yoğunluğu dengesi, elektron enerji dengesi ve elektrostatik Poisson denkleminden oluşur:

\[ \frac{\partial n_j}{\partial t}+\nabla\cdot\Gamma_j=S_j \]

\[ \frac{\partial w_e}{\partial t}+\nabla\cdot Q_e=-e_0\Gamma_e\cdot E+\widetilde{S}_e \]

\[ -\nabla\cdot(\varepsilon_r\varepsilon_0\nabla\phi)=\sum_j q_j n_j \]

Burada \(n_j\), \(j\) türünün parçacık sayı yoğunluğunu; \(\Gamma_j\) parçacık akısını; \(S_j\) reaksiyonlardan kaynaklanan net parçacık üretim veya kayıp terimini; \(w_e=n_eu_e\) elektron enerji yoğunluğunu; \(Q_e\) elektron enerji akısını; \(\widetilde{S}_e\) elektron enerji kaynak/kayıp terimini; \(\phi\) elektrik potansiyelini ve \(E=-\nabla\phi\) elektrik alanını temsil eder.

Ağır parçacıkların akısı sürüklenme-difüzyon yaklaşımıyla

\[ \Gamma_h=\operatorname{sgn}(q_h)n_hb_hE-D_h\nabla n_h \]

şeklinde tanımlanır. \(b_h\) hareketlilik, \(D_h\) ise ağır parçacığın difüzyon katsayısıdır.

MCPlas’ın üç elektron taşınım seçeneği

Araç kutusu elektron akısı ve elektron enerji akısı için üç farklı yaklaşım sunmaktadır: DDAc, DDA53 ve DDAn.

Geleneksel DDAc yaklaşımında elektron akısı

\[ \Gamma_e=-n_eb_eE-\nabla(D_en_e) \]

şeklindedir. DDA53, elektron parçacık akısında aynı formu kullanırken elektron enerji akısını yaygın \(5/3\) katsayılı yaklaşım ile tanımlar:

\[ Q_e=-\frac{5}{3}w_eb_eE-\frac{5}{3}\nabla(D_ew_e) \]

MCPlas’ın özellikle öne çıkardığı DDAn seçeneği ise elektron hız dağılım fonksiyonunun Legendre polinomlarıyla açılımından ve elektron Boltzmann denkleminin ilk dört momentinden türetilmiş daha ayrıntılı bir taşınım tanımıdır. Elektron parçacık akısının temel biçimi:

\[ \Gamma_e= -\frac{e_0}{m_e\nu_e}\nabla\left[(\xi_0+\xi_2)n_e\right] -\frac{e_0}{m_e\nu_e}En_e \]

ve elektron enerji akısı:

\[ Q_e= -\frac{e_0}{m_e\widetilde{\nu}_e} \nabla\left[(\widetilde{\xi}_0+\widetilde{\xi}_2)w_e\right] -\frac{e_0}{m_e\widetilde{\nu}_e} \left(\frac{5}{3}+\frac{2}{3}\frac{\xi_2}{\xi_0}\right)Ew_e \]

olarak verilmektedir. Burada \(\nu_e\) ve \(\widetilde{\nu}_e\) momentum ve enerji akısı sönüm frekanslarını, \(\xi_0\), \(\xi_2\), \(\widetilde{\xi}_0\) ve \(\widetilde{\xi}_2\) ise elektron hız dağılım fonksiyonundan hesaplanan taşınım katsayılarını temsil eder.

Bu yaklaşımın sonucu, elektron taşınımının yalnız standart hareketlilik-difüzyon katsayılarıyla değil, elektron hız dağılımının daha ayrıntılı moment bilgisiyle ilişkilendirilmesidir. Çalışmanın testleri, bu seçimin özellikle elektrot sınırlarına yakın bölgelerde hesaplanan plazma özelliklerini belirgin biçimde değiştirebildiğini göstermektedir.

Reaksiyon kaynak terimleri nasıl otomatik oluşturuluyor?

Bir parçacık türünün reaksiyonlardan doğan kaynak terimi çalışmada

\[ S_j=\sum_{l=1}^{N_r}(G_{jl}-L_{jl})R_l \]

ve bir reaksiyonun hızı

\[ R_l=k_l\prod_{i=1}^{N_s}n_i^{\beta_{il}} \]

ile tanımlanmıştır. \(G_{jl}\) ve \(L_{jl}\), ilgili reaksiyonda tür \(j\) için üretim ve kayıp stokiyometrisini; \(k_l\) reaksiyon hız katsayısını; \(\beta_{il}\) kısmi reaksiyon derecesini ifade eder.

MCPlas bu üretim ve kayıp matrislerini JSON tabanlı reaksiyon kinetiği girdilerinden otomatik oluşturur. Bunun pratik sonucu, daha karmaşık bir kimya modeline geçildiğinde yüzlerce reaksiyonun COMSOL içine tek tek yeniden girilmesine gerek kalmamasıdır.

Sayısal kararlılık nasıl destekleniyor?

Araç kutusu iki tür kararlılık tekniği sunmaktadır. Bunlardan biri parçacık yoğunlukları ve ortalama elektron enerjisine logaritmik dönüşüm uygulanmasıdır. Genel denge denklemi

\[ \frac{\partial P}{\partial t}+\nabla\cdot F=S \]

için \(p=\ln(P)\) dönüşümü kullanıldığında zamansal terim \(e^p\partial p/\partial t\) biçimine gelir. Çalışmaya göre bu yaklaşım çözümün pozitif kalmasına yardım eder ve çok düşük yoğunluklu veya keskin gradyanlı bölgelerdeki salınımları bastırır.

İkinci yöntem, çok düşük çözüm değerlerinde sayısal kararsızlığı dengelemek üzere kaynak terimine eklenen bir stabilizasyon terimidir:

\[ N_A\exp[-\xi\ln(p)] \]

Burada \(N_A\) Avogadro sabiti, \(\xi\) ise 0,25 ile 1 arasında seçilebilen stabilizasyon parametresidir; kaynakta varsayılan değer 1 olarak verilmektedir.

Türkiye açısından ne ifade ediyor?

Çalışma Türkiye’ye özgü deney, kurum, endüstriyel uygulama veya yerel plazma verisi içermemektedir. Bu nedenle sonuçlardan Türkiye için doğrudan performans veya başarı oranı çıkarılamaz. Bununla birlikte çalışma yöntemsel düzeyde MATLAB ve COMSOL kullanan araştırma grupları açısından önemlidir: reaksiyon kinetiğinin şema doğrulamalı ve yeniden kullanılabilir bir veri yapısına dönüştürülmesi, aynı modeli farklı araştırmacıların veya farklı sayısal araçların daha kontrollü biçimde yeniden kurabilmesini hedeflemektedir. Bunun Türkiye’deki bir araştırma ortamındaki gerçek yararı ancak yerel iş akışlarında uygulanıp doğrulandığında değerlendirilebilir.

Çalışma ne söylüyor, ne söylemiyor?

Çalışmanın desteklediği sonuçlar

• MATLAB tabanlı MCPlas, JSON girdilerinden denklem tabanlı COMSOL akışkan-Poisson modelleri otomatik oluşturabilmektedir.
• MCPlas ile COMSOL Plasma Module aynı model denklemleri ve aynı sınır koşulları kullanıldığında DC ve RF argon testlerinde çok yakın sonuçlar üretmektedir.
• Elektron taşınım yaklaşımı ve sınır koşullarının seçimi, özellikle elektrot kılıf bölgelerinde model sonuçlarını önemli ölçüde değiştirebilmektedir.
• 23 tür ve 409 süreç içeren ayrıntılı argon reaksiyon kinetiği modeli MCPlas ile otomatik olarak ele alınabilmiştir.
• Aynı JSON reaksiyon kinetiği girdilerinin MCPlas, PLASIMO ve FEDM arasında yeniden kullanılabilmesi test edilmiştir.
• Bu üç yazılımın sonuçları, incelenen DC ve RF argon referans durumlarında yalnızca birkaç yüzde düzeyinde en büyük farklılıklar göstermiştir.

Çalışmanın kanıtlamadığı sonuçlar

• MCPlas’ın bütün gazlarda, bütün basınçlarda veya bütün plazma geometrilerinde otomatik olarak en doğru sonucu verdiği gösterilmemiştir.
• DDAn yaklaşımının her fiziksel durumda diğer bütün elektron taşınım modellerinden üstün olduğu bu çalışmayla evrensel olarak kanıtlanmamıştır.
• Çalışma yeni bir deneysel plazma ölçüm kampanyası değildir.
• COMSOL ile uyum, gerçek fiziksel sistemin her ayrıntısının deneysel olarak doğrulandığı anlamına gelmez.
• Dört tür ve 23 tür modeller arasındaki farklar, daha ayrıntılı modelin bütün kullanım koşullarında zorunlu olduğu anlamına gelmez; model karmaşıklığının gerekliliği araştırılan fiziksel büyüklüğe bağlıdır.
• İncelenen 1 cm, 1 Torr argon DC/RF referans durumlarının sonuçları doğrudan başka gazlara, basınçlara, boyutlara veya endüstriyel plazma sistemlerine aktarılamaz.

Çalışmanın Yöntemi ve Bulguları

Referans simülasyon koşulları

MCPlas iki düşük basınçlı argon parıltı boşalmasıyla sınanmıştır. Her iki testte de uzaysal olarak tek boyutlu geometri kullanılmış, sol elektrot topraklanmış ve sağ elektrot beslenmiştir.

İlk doğrulama: MCPlas ve COMSOL Plasma Module aynı modeli çözerse ne oluyor?

İlk doğrulama aşamasında MCPlas, COMSOL Plasma Module ile mümkün olduğunca aynı fiziksel model tanımını kullanacak biçimde yapılandırılmıştır. Elektron parçacık akısında DDAc biçimi ve enerji akısında DDA53 yaklaşımı kullanılmış; sınır koşulları da COMSOL Plasma Module’ün varsayılan modeline uyarlanmıştır.

DC boşalmasında Şekil 9’daki çizgiler ve semboller parçacık sayı yoğunlukları, ortalama elektron enerjisi ve elektrik alan için neredeyse üst üste gelmektedir. Aynı durum RF boşalması için Şekil 11’de görülmektedir. Yazarlar bu uyumu, MCPlas’ın amaçlanan denklemleri COMSOL ortamında doğru biçimde oluşturduğuna yönelik kod doğrulaması olarak yorumlamaktadır.

Bu sonuç özellikle önemlidir çünkü karşılaştırma “hangi fizik modeli daha doğrudur?” sorusuna değil, “aynı fizik modeli iki farklı kurulum yolu ile verdiğimizde aynı sayısal davranışı elde ediyor muyuz?” sorusuna yanıt vermektedir.

DDAn ve gelişmiş sınır koşulları etkinleştirildiğinde ne değişiyor?

İkinci aşamada MCPlas’ın daha gelişmiş DDAn elektron taşınım modeli ile çalışmada verilen elektron, elektron enerjisi ve ağır parçacık sınır koşulları kullanılmıştır. Bu noktada MCPlas ve standart COMSOL Plasma Module sonuçları artık aynı fiziksel formülasyonu temsil etmemektedir.

DC durumunda Şekil 10, parçacık sayı yoğunluklarının boşluk boyunca farklılaştığını; en belirgin ayrışmanın katot kılıfına yakın bölgede oluştuğunu göstermektedir. MCPlas modeli bu bölgede daha düşük sayı yoğunlukları ve daha geniş bir katot-kılıf bölgesi öngörmektedir. Ortalama elektron enerjisi ve elektrik alan profilleri de büyüklük ve uzaysal dağılım açısından farklılaşmaktadır.

RF durumundaki Şekil 12’de de benzer bir fiziksel duyarlılık görülmektedir. Yüklü parçacık yoğunlukları elektrot kılıflarında farklılaşırken, uyarılmış argon durumu için hesaplanan yoğunluklar boşluğun büyük bölümünde bir mertebeden fazla ayrılabilmektedir. Ortalama elektron enerjisi ile elektrik alan dağılımları da iki formülasyon arasında belirgin fark göstermektedir.

Bu bulgu, modelin “küçük” görünen elektron taşınımı ve sınır koşulu tercihlerinin yalnız sayısal ayrıntı olmadığını; hesaplanan plazma bileşimini ve alan dağılımını değiştirebilen fiziksel model tercihleri olduğunu göstermektedir.

Dört türden 23 türe geçildiğinde ne oluyor?

Basit reaksiyon kinetiği modeli dört tür içermektedir: elektronlar, temel durumdaki argon atomu, Ar⁺ iyonu ve Ar[1s5] ile Ar[1s3] metastabil durumlarını tek bir toplulaştırılmış uyarılmış durum olarak temsil eden Ar*.

Ayrıntılı model ise toplam 23 parçacık türü kullanır. Modelde atomik ve moleküler iyonlar, farklı uyarılmış moleküler durumlar, bireysel 1s ve 2p argon seviyeleri ve yüksek enerjili durumları temsil eden toplulaştırılmış bir seviye bulunur. Bu 23 tür, toplam 409 çarpışma süreci ve radyatif geçiş içinde yer almaktadır.

DC sonuçlarında elektron ve Ar⁺ yoğunluklarının genel uzaysal profilleri dört ve 23 tür modellerinde benzer kalırken, uyarılmış durumların yoğunluklarında çok daha büyük farklar ortaya çıkmaktadır. Kaynakta, dört tür modelindeki Ar* yoğunluğu ile 23 tür modelindeki Ar[1s5] ve Ar[1s3] seviyelerinin toplamı arasında bir mertebeden fazla fark bulunduğu belirtilmektedir.

RF durumunda yüklü türlerin sonuçları yine birbirine daha yakın olmakla birlikte, ortalama elektron enerjisinde ve uyarılmış durum yoğunluklarında belirgin sapmalar vardır. Ayrıca 23 tür modeli, dört tür modelinde hiç bulunmayan moleküler argon türlerinin uzaysal dağılımını da sağlayabilmektedir.

Aynı JSON verisi başka plazma kodlarında kullanılabilir mi?

Çalışmanın yeniden üretilebilirlik açısından en güçlü testlerinden biri, aynı dört tür ve 23 tür argon reaksiyon kinetiği JSON dosyalarının iki farklı modelleme aracı olan PLASIMO ve FEDM’de de kullanılmasıdır.

Karşılaştırma yalnız farklı yazılım paketleri arasında değildir; kullanılan sayısal yöntemler de farklıdır. PLASIMO sonlu hacim yaklaşımına, FEDM ise sonlu eleman yaklaşımına dayanmaktadır. Buna rağmen aynı fiziksel model ve aynı JSON reaksiyon kinetiği girdileri kullanıldığında MCPlas, PLASIMO ve FEDM sonuçları hem DC hem RF testlerinde yüksek derecede uyumludur.

Çalışma, üç araç arasındaki en büyük farkın yalnızca birkaç yüzde düzeyinde olduğunu belirtmektedir. Bu sonuç, JSON dosyasının yalnız MCPlas’a özel bir konfigürasyon dosyası olmadığını; reaksiyon kinetiğini farklı modelleme ortamlarına taşıyabilen yeniden kullanılabilir bir bilimsel veri yapısı olarak işlev görebildiğini desteklemektedir.

Çalışmanın güçlü yönleri

• Fiziksel denklemler ve sınır koşulları MATLAB kaynak kodunda açık biçimde görülebilmektedir.
• Reaksiyon kinetiği ve model tanımı, yazılım projesinden ayrılmış şema doğrulamalı JSON veri yapılarında tutulmaktadır.
• Basit ve çok kapsamlı reaksiyon kinetiği modelleri aynı otomasyon mimarisiyle ele alınmıştır.
• Model doğrulaması yalnız tek bir örneğe değil, DC ve RF olmak üzere iki referans boşalmaya uygulanmıştır.
• JSON veri yeniden kullanılabilirliği iki bağımsız plazma modelleme kodu ile sınanmıştır.
• Model karşılaştırması farklı sayısal çözüm yöntemlerini kullanan yazılımlar arasında gerçekleştirilmiştir.

Temel sınırlılıklar

• Sunulan referans testleri düşük basınçlı argon ve belirli 1D DC/RF boşalma koşullarıyla sınırlıdır.
• Çalışmanın ana doğrulaması model-kod karşılaştırmasına dayanır; bu makale yeni bir deneysel ölçüm kampanyası sunmamaktadır.
• Şekil 9–17’de çok sayıda uzaysal profil gösterilmesine rağmen eğrilerin bütün ham sayısal noktaları makale tablolarında verilmemektedir.
• Farklı elektron taşınım modelleri arasında görülen farklar, belirli bir yaklaşımın bütün plazma rejimlerinde üstün olduğunu tek başına kanıtlamaz.
• Diğer MATLAB ve COMSOL sürümleri için yazarlar kısıtlama beklemediklerini belirtse de çalışma farklı yazılım sürümlerinin sistematik uyumluluk testini raporlamamaktadır.

Kaynak ve Yöntem Notu

Tam özgün çalışma adı: MCPlas, a MATLAB toolbox for reproducible plasma modelling with COMSOL

Yazarlar ve sıralama: Marjan N. Stankov; Daan Boer; Wouter Graef; Kevin van ‘t Veer; Aleksandar P. Jovanović; Florian Sigeneger; Detlef Loffhagen; Jan van Dijk; Markus M. Becker.

Sorumlu yazar: Marjan N. Stankov. Yüklenen sürümde eş birinci yazar/eş katkı beyanı bulunmamaktadır.

Kurumlar: Leibniz Institute for Plasma Science and Technology (INP), Greifswald, Almanya; Eindhoven University of Technology, Eindhoven, Hollanda; Plasma Matters B.V., Eindhoven, Hollanda.

Kaynak türü: MATLAB tabanlı bilimsel yazılım geliştirme, plazma modelleme ve sayısal doğrulama çalışması.

Yüklenen sürümün yayın durumu: Yüklenen PDF, 19 Aralık 2025 tarihinde Computer Physics Communications için gönderilmiş ve üzerinde açıkça hakem değerlendirmesinden geçmemiş preprint olduğu belirtilmiş sürümdür. Bu çalışma hakem değerlendirmesinden geçmemiş bir preprint sürümü olarak incelenmiştir; bilimsel içerik bu yüklenen sürüme dayanmaktadır.

Preprint DOI: 10.2139/ssrn.6056661.

Daha sonra yayımlanan sürüm: Bibliyografik doğrulamada aynı başlıklı çalışma 2026 yılında Computer Physics Communications dergisinde, Cilt 327, makale 110248 olarak yayımlanmış görünmektedir.

Hakemli yayın DOI: 10.1016/j.cpc.2026.110248.

Dergi: Computer Physics Communications.

Yayınevi: Elsevier B.V.

Lisans: Yayımlanmış sürüm açık erişim ve Creative Commons lisansı altında listelenmektedir; bibliyografik kayıtlarda CC BY 4.0 olarak belirtilmektedir. MCPlas yazılımının kendi lisansı kaynakta MIT olarak verilmektedir.

Yazılım: Araç kutusunun programlama dili MATLAB’dır. Preprint program özetinde geliştirme ve test için MATLAB R2023a, COMSOL Multiphysics 6.2 ve LiveLink for MATLAB kullanıldığı belirtilmektedir.

Kaynak içi sürüm tutarsızlığı: Program özetinde MATLAB R2023a yazarken kaynakçadaki MATLAB girdisi R2023b olarak verilmiştir. Kaynak bu iki ifadeyi uzlaştıran ayrı bir açıklama sunmadığı için burada sessizce tek bir sürüme dönüştürülmemiştir.

Model: Elektron, iyon ve nötr parçacık sayı yoğunlukları; ortalama elektron enerjisi; elektrik potansiyeli ve elektrik alan için zaman ve uzaya bağlı akışkan-Poisson modeli.

Referans koşulları: 1 cm boşluk, 1 Torr argon ve 300 K gaz sıcaklığı. DC testinde 350 V sabit gerilim; RF testinde 350 V genlik ve 13,56 MHz frekans kullanılmıştır.

Reaksiyon kinetiği: Dört türden oluşan basit argon modeli ile 23 tür ve 409 çarpışma/radyatif süreç içeren genişletilmiş argon modeli değerlendirilmiştir.

Veri erişilebilirliği: Yüklenen preprint, çalışma verilerinin makale kabulünden sonra kamuya açılacağını belirtmektedir. Güncel bibliyografik doğrulamada ilgili veri seti 10.34711/inptdat.1033 DOI’siyle yayımlanmış durumdadır.

Kod erişilebilirliği: Preprint, MCPlas kaynak kodu ve giriş verilerinin INP-SDT/MCPlas geliştirici deposunda yayımlandığını ve incelenen sürüm için 0b67eea commit kimliğini vermektedir.

Finansman: Çalışma kısmen Deutsche Forschungsgemeinschaft tarafından 213099267, 368502453, 407462159, 504701852, 509169873 ve 535827833 proje numaraları kapsamında desteklenmiştir.

Çıkar çatışması: Yüklenen 40 sayfalık preprintte ayrı bir çıkar çatışması beyanı tespit edilmemiştir; bundan çıkar çatışması olmadığı sonucu çıkarılmamıştır.

CRediT/yazar katkıları: Yüklenen sürümde ayrı bir CRediT veya ayrıntılı yazar katkı beyanı bulunmamaktadır.

Bilimsel içerik sınırı: Bu Verianla makalesindeki yöntem, denklemler, model koşulları ve bilimsel bulgular yüklenen çalışmaya dayanmaktadır. Dış doğrulama yalnız başlık, DOI, yayımlanmış sürüm, dergi, veri erişimi ve yayın kimliği gibi bibliyografik alanlar için kullanılmış; dış kaynaklardan yeni deneysel veya sayısal bulgu eklenmemiştir.

Temel yorum sınırı: Çalışma bir yazılım ve sayısal modelleme doğrulamasıdır. MCPlas’ın bütün plazma koşullarında fiziksel olarak en doğru model olduğu sonucunu desteklemez. Sunulan temel karşılaştırmalar düşük basınçlı argon için 1D DC ve RF referans durumlarıyla sınırlıdır.