Бул araştırma, төмөн тыгыздыктагы бир атом bulutu içindeki tuzaklanmış iyonla gerçekleşen klasik атом–ион etkileşimlerini MATLAB үстүндө жогорку paralellikle симуляция etmek үчүн geliştirilen ode45gpu жана atomiongpu.m yazılımlarını incelemektedir. Sistem, Newton hareket denklemleri 12 bileşenli бир абал vektörüne dönüştürülerek modellenmekte; dördüncü/beşinci mertebe adaptif Runge–Kutta integrasyonu gerçekleştirilmekte жана birbirinden bağımsız баштапкы koşulları MATLAB'ın gpuArray жана arrayfun mekanizmalarıyla GPU үстүндө paralelleştirilmektedir. Yüklenen preprint sürümündeki benchmark'a göre özelleştirilmiş ode45gpu, tek CPU çekirdeğinde referans ode45 kullanımından 22 kata kadar hızlı çalışmış; GPU testlerinde yaklaşık 10.000 yörüngeye kadar geniş бир параллелдүүлүк bölgesi gözlenmiş жана örnek бир 8-GPU çalıştırmasında 10 milyondan fazla орбита 15 saatin астында tamamlanmıştır. Temel sınır, bunun klasik жана belirli атом–ион өз ара аракеттенүү modeline göre özelleştirilmiş бир симуляция olmasıdır; genel amaçlı ode45'in her problemde GPU үстүндө doğrudan yerine geçen evrensel бир чечүүчү катары gösterilmemiştir.

Изилдөө кайсы маселени чечүүгө аракет кылат?

Soğuk жана ultrasoğuk атом–ион fiziğinde бир ион, Paul tuzağı сыяктуу elektromanyetik бир tuzak ичинде tutulurken serbest бир атом iyonla etkileşebilir. Tek бир çarpışmanın klasik hareket denklemleri görece küçük бир система oluşturmasına rağmen gerçek problem istatistikseldir: atomun баштапкы yönü, hızı жана diğer баштапкы koşulları tam катары bilinmediği үчүн tek бир орбита бардык sistemi temsil etmez.

Бул nedenle binlerce, hatta milyonlarca ар түрдүү баштапкы koşulunun ayrı ayrı симуляция edilmesi gerekir. Her орбита diğerinden bağımsız экендиги үчүн problem, paralel эсептөө açısından son derece elverişlidir. Çalışmanın temel fikri de budur: birbirini beklemek zorunda olmayan çok sayıda küçük ODE problemini GPU'nun çok sayıdaki işlem birimine dağıtmak.

Эмне үчүн эки бөлүкчө гана симуляцияланат?

Kaynak, төмөн тыгыздыктагы tampon gaz koşullarında iyonun бирдей anda esas катары tek бир atomla etkileştiğini varsaymaktadır. Böylece система гана бир атом жана бир iyondan oluşan iki parçacıklı klasik modele indirgenebilir.

Бул sadeleştirme GPU açısından avantajlıdır. Her орбита az miktarda bellek kullanır жана yörüngelerin tamamının sürekli bellekte tutulması zorunlu эмес. Hesaplamanın büyük kısmı birbirinden bağımsız küçük sayısal integrasyonlardan oluşur.

Benchmark системасында кайсы атом жана ион колдонулган?

Benchmark үчүн изилдөө, merkezdeki Paul tuzağında başlangıçta hareketsiz bulunan бир Yb⁺ iyonu менен yaklaşık 5000 Bohr yarıçapı uzaktan başlatılan бир Rb atomunu kullanmaktadır. Kullanıcı arayüzündeki örnek ayarlarda bunlar өзгөчө ¹⁷¹Yb⁺ жана ⁸⁷Rb катары берилет.

Atomun benchmark баштапкы hızı \(T=1\,\mu\mathrm{K}\) sıcaklığına göre:

\[ v_0=\sqrt{\frac{3k_BT}{m_a}} \]

менен belirlenir жана başlangıçta iyona doğru yönlendirilir.

Атом–ион өз ара аракеттенүүсү кайсы потенциал менен моделделген?

Benchmark'ta атом–ион etkileşimi uzun menzilli çekici жана kısa menzilli itici terimin birleşimiyle ifade жасалат:

\[ V_{ai}(r)=-\frac{C_4}{r^4}+\frac{C_8}{r^8}. \]

Burada ilk terim uzun menzildeki yük–indüklenmiş dipol etkileşimini temsil ederken, ikinci terim атом менен ион çok yaklaştığında kısa menzilli itici davranışı modellemektedir.

atomiongpu.m гана бул өзгөчө 4–8 biçimine kilitli эмес. Kullanıcı arayüzü көбүрөөк genel:

\[ V(r)=-\frac{C_n}{r^n}+\frac{C_m}{r^m} \]

ifadesini desteklemektedir.

Paul тузагы Гамильтонианга кантип кирет?

Benchmark sisteminin Hamiltonyen'i атом жана iyonun kinetik enerjilerini, атом–ион etkileşimini жана ион üzerindeki zamana bağlı Paul tuzağı potansiyelini birlikte içerir:

\[ H= \frac{1}{2}m_av_a^2+ \frac{1}{2}m_iv_i^2 -\frac{C_4}{r_{ai}^4} +\frac{C_8}{r_{ai}^8} + \sum_{j=1}^{3} (a_j+2q_j\cos\Omega_{\mathrm{rf}}t) \frac{m_{\mathrm{ion}}\Omega_{\mathrm{rf}}^2}{8} r_j^2. \]

\(a_j\) жана \(q_j\) tuzak parametreleri, \(\Omega_{\mathrm{rf}}\) болсо radyo-frekans tuzak sürme frekansıdır. Zamana bağlı tuzak nedeniyle sistemin kütle merkezi менен göreli hareketi basit biçimde birbirinden ayrılamaz.

Эмне үчүн абал вектору 12 өлчөмдүү?

Бир атом жана бир iyonun her biri üç uzaysal koordinat жана üç ылдамдык bileşeni taşır. Dolayısıyla:

Toplamda 12 değişkenli бир абал вектору oluşur жана hareket:

\[ \dot{\mathbf y}=f(t,\mathbf y) \]

biçiminde birinci mertebe adi diferansiyel denklem sistemine айландырылат.

Атом–ион аралыгы кантип эсептелет?

\[ r= \sqrt{ (y_1-y_7)^2+ (y_2-y_8)^2+ (y_3-y_9)^2 }. \]

Benchmark өз ара аракеттенүү potansiyelinin türevi болсо:

\[ \frac{dV}{dr} ============= ## \frac{4C_4}{r^5} \frac{8C_8}{r^9}. \]

Бул türev атом жана iyonun karşılıklı ivmelerini hesaplamak үчүн колдонулат. İyon denklemlerine buna ek катары Paul tuzağının zamana bağlı күч terimleri eklenmektedir.

Баштапкы позиция эмне үчүн сферадан тандалат?

Atom, iyonun merkezde экендиги \(r_0\) yarıçaplı бир küre үстүндө başlatılır. Başlangıç koordinatları küresel açılar \(\theta\) жана \(\phi\) üzerinden:

\[ x=r_0\sin\theta\cos\phi, ] [ y=r_0\sin\theta\sin\phi, ] [ z=r_0\cos\theta \]

şeklindedir. Head-on durumda ылдамдык вектору бул konum vektörünün ters yönünde seçilir; yani атом doğrudan merkeze doğru hareket eder.

Neden her θ değerini eşit aralıkla seçmek doğru эмес?

Kürenin yüzey alanı kutuplara yaklaşıldıkça ар түрдүү geometrik ağırlığa sahiptir. Бул nedenle \(\theta\)'yı гана doğrusal aralıklarla örneklemek küre үстүндө eşit alan örneklemesi oluşturmaz. Çalışmada eşit katı açı örneklemesi үчүн:

\[ \theta=\arccos(1-2u) \]

dönüşümü колдонулат. Burada \(u\), seçilen \(\theta\) sınırlarına karşılık gelen uygun aralıkta eşit aralıklı же rastgele seçilmektedir. Böylece her \((\theta,\phi)\) баштапкы koşulu yaklaşık eşit yüzey alanını temsil eder.

Орбиталар эмне үчүн кээде өтө татаал болуп калат?

Bazı баштапкы koşullarında атом iyonla гана kısa süre etkileşip geri saçılır. Bazılarında болсо атом жана ион tekrar tekrar yakınlaşarak geçici катары birbirine bağlı, uzun ömürlü бир атом–ион kompleksi meydana getirebilir.

Çalışmanın örnek yörüngesi бул ikinci davranışa sahiptir. \(\theta\approx1{,}02\) жана \(\phi\approx0{,}22\) civarındaki өзгөчө баштапкы koşulu, yoğun жана düzensiz бир орбита üretmektedir. Yazarlar bunu ODE sisteminin güçlü doğrusal olmayan жана kaotik karakterinin бир örneği катары kullanmıştır.

MATLABдын стандарттык ode45 чечүүчүсү бул жерде кантип иштейт?

ode45, Dormand–Prince çiftine dayalı dördüncü/beşinci mertebe adaptif zaman adımlı Runge–Kutta yöntemidir. Kullanıcı:

• күч fonksiyonunu \(f(t,y)\),
• баштапкы koşullarını,
• симуляция zamanını,
• mutlak toleransı,
• bağıl toleransı

çözücüye verir.

Çözücü her adımda yerel hata kestirimi yapar. Hata toleranstan büyükse zaman adımı küçültülür; hata kabul edilebilir düzeydeyse adım korunabilir же büyütülebilir.

MATLAB кодун кайра жазуу эле ылдамдыкты болжол менен үч эсе кантип арттырды?

Çalışmanın dikkat çekici bulgularından biri, GPU'ya geçmeden мурун bile күч fonksiyonunun yazım biçiminin büyük performans farkı yaratmasıdır.

İlk күч fonksiyonunda MATLAB'ın:

y(4:6)

сыяктуу бир ifade үчүн geçici üç elemanlı бир vektör oluşturduğu, ardından elemanların ayrı ayrı kullanıldığı белгиленет. Yazarlar bunu:

y(4), y(5), y(6)

сыяктуу skaler erişimlere dönüştürmüştür.

Kaynağın benchmark'ında 2500 орбита, tek Intel Haswell CPU çekirdeğinde ilk күч fonksiyonuyla ode45 kullanıldığında 10 saatten fazla sürerken, optimize edilmiş күч fonksiyonuyla бирдей ode45 hesabı 3,5 saatin altına düşmüştür. Yazarlar bunu yaklaşık üç kat hızlanma катары yorumlamaktadır.

Бул sonuç MATLAB öğrenen biri үчүн önemli бир noktayı gösterir: бирдей matematiksel denklem, ар түрдүү MATLAB kodlama biçimleriyle çok ар түрдүү изилдөө sürelerine sahip olabilir.

ode45gpu деген эмне?

ode45gpu, изилдөөнүн yazarlarının Dormand–Prince dördüncü/beşinci mertebe adaptif Runge–Kutta algoritmasını yeniden uygulayarak oluşturduğu, tek dosyalı жана атом–ион problemi үчүн özelleştirilebilir MATLAB fonksiyonudur.

Standart ode45 менен temel algoritmik mantığı aynıdır; бирок uygulama жана çıktı politikası farklıdır.

ode45gpu varsayılan катары бардык yörüngeyi сактоо yerine:

• son zamanı,
• 12 son абал değişkenini,
• жана kullanıcı tarafından seçilen төмөн bellekli skaler gözlenebilirleri

döndürür.

Ayrıca tek simülasyonun olağanüstü uzun sürmesini engellemek үчүн 1.000.000 zaman adımı sonrasında döngüyü sonlandıran güvenlik koşulu бар.

Эмне үчүн бүт орбитаны сактабоо чоң артыкчылык?

Milyonlarca yörüngenin her zaman adımındaki 12 абал değişkenini сактоо bellek ihtiyacını dramatik biçimde büyütürdü. Bunun yerine изилдөө, birçok gözlenebiliri integrasyon sırasında “on the fly” hesaplamakta жана geçmiş adımların tamamını bellekte tutmamaktadır.

Örneğin:

• iyonun son kinetik enerjisi,
• kompleks ömrü,
• hesaplanan zaman adımı sayısı

гана gerekli bilgiler колдонуу менен kaydedilebilir.

Neden ode45gpu ичинде MATLAB vektörleri parçalanmıştır?

Çalışmanın GPU uygulamasındaki kritik программалоо kararı budur. MATLAB'ın GPU үстүндө arrayfun aracılığıyla çalıştırdığı fonksiyonlarda kullanılan kod kısıtları nedeniyle yazarlar 12 elemanlı абал vektörünü ayrı skaler değişkenlere dönüştürmüştür.

Örneğin geleneksel катары:

f(1:3)=y(4:6)

түрүндө ifade edilebilecek işlem mantıksal катары üç ayrı skaler atamaya açılmıştır.

Бул nedenle ode45gpu genel MATLAB ODE kodundan көбүрөөк uzun görünse de GPU yürütme modeline көбүрөөк uygun бир yapı kazanmıştır.

Kuvvet fonksiyonu эмне үчүн yedi kez bulunuyor?

Dormand–Prince yöntemi her zaman adımında birden fazla ara eğim değerlendirmesi yapar. Çalışmanın ode45gpu uygulamasında күч fonksiyonunun skalerleştirilmiş ifadesi yedi ayrı yerde бар жана her aşamada ар түрдүү ara абал değişkenleri колдонулат.

Бул, kaynak kodu açısından tekrar сыяктуу görünse de GPU үчүн fonksiyon çağrısı жана vektör oluşturma yükünü azaltmaya yönelik özelleştirilmiş uygulamanın бир parçasıdır.

gpuArray бул жерде так эмне кылат?

Her баштапкы koşulu үчүн бир matris elemanı oluşturulur. Örneğin atomun баштапкы x koordinatlarının tamamı бир matris halinde hesaplanıp GPU belleğine aktarılır. Aynı işlem diğer 11 баштапкы değişkeni үчүн de yapılır.

Kuvvet sabitleri жана hata toleransları değişmiyorsa bile arrayfun çağrısıyla uyumlu olmak үчүн баштапкы matrisleriyle бирдей boyutta, her elemanı бирдей sabit değeri taşıyan GPU dizileri oluşturulur.

arrayfun бул жерде кантип параллелдүүлүктү камсыз кылат?

arrayfun, karşılık gelen GPU dizi elemanlarını bağımsız birer fonksiyon çağrısına dönüştürür. Бул çalışmada her indeks бир yörüngeyi temsil eder. Dolayısıyla GPU бирдей anda çok sayıda:

\[ \dot y=f(t,y) \]

problemini ар түрдүү баштапкы koşullarıyla çözebilir.

gather эмне үчүн эң аягында колдонулат?

ode45gpu чыгыштары GPU belleğinde gpuArray катары kalır. Çalışma, son analiz жана veri işleme aşamasını CPU tarafında yapmayı tercih ettiği үчүн sonuçlar:

GPU → gather → CPU belleği

zinciriyle geri alınmaktadır.

CPU'da ode45gpu gerçekten эмне үчүн hızlı?

Yüklenen v1 benchmark'ında, 2500 орбита үчүн ode45gpu tek CPU çekirdeğinde, yazarların referans катары kullandığı ode45 + f(t,y) uygulamasından 22 kat hızlı çalışmıştır. Çalışma мындан тышкары бул ылдамдык oranının 28 bağımsız CPU çekirdeğine kadar paralelleştirmede korunduğunu bildirmektedir.

Бул sonucu “MATLAB ode45 her problemde 22 kat yavaştır” biçiminde okumak doğru эмес. Benchmark, belirli атом–ион problemi, belirli күч fonksiyonu жана бул өзгөчө çıktı politikasına aittir. Nitekim изилдөөнүн көбүрөөк sonraki hakemli sürümünde CPU hızlanma değeri yeniden değerlendirilmiştir.

GPU benchmark'ında кайсы donanımlar kullanıldı?

Çalışma üç NVIDIA Tesla GPU үстүндө test yapmıştır:

Her testte бирдей өзгөчө yörüngenin bağımsız kopyalarının sayısı 1'den 1.000.000'a kadar artırılmıştır.

Yaklaşık 10.000 орбита eşiği эмне anlama geliyor?

Üç GPU'da da yaklaşık 10.000 орбита civarında belirgin бир geçiş görülmüştür. Бул sayının астында орбита sayısı artırıldığında toplam изилдөө süresi yaklaşık sabit kalmaktadır. Yazarların yorumuna göre бул bölge, GPU'nun yeni yörüngeleri mevcut paralel işlem kapasitesi içine yerleştirebildiğini gösterir.

Eşik aşıldıktan кийин изилдөө süresi орбита sayısıyla yaklaşık doğrusal artmaya başlar. GPU'nun eşzamanlı işlem kapasitesi dolmuştur жана ek орбиталар yeni изилдөө dalgaları halinde işlenmektedir.

Aynı yörüngeyi kopyalamak benchmark'ı yanıltıyor olabilir mi?

Yazarlar bunu мындан тышкары test etmiştir. Tesla K80 үстүндө бирдей yörüngenin çok sayıda kopyası менен küre yüzeyinden eşit katı açıyla örneklenmiş ар түрдүү баштапкы koşulları karşılaştırılmıştır.

Farklı yörüngelerde өзгөчө küçük örnek sayılarında изилдөө süresinde көбүрөөк fazla dalgalanma görülmüştür; çünkü bazı баштапкы koşulları көбүрөөк karmaşık жана uzun dinamikler üretmektedir. Yaklaşık параллелдүүлүк eşiğinden кийин iki durumda da изилдөө süresi орбита sayısıyla yaklaşık doğrusal hale gelmiştir.

atomiongpu.m эмне үчүн geliştirildi?

ode45gpu'yu doğrudan kullanmak, күч fonksiyonunu жана GPU girişlerini ayrıntılı biçimde düzenlemeyi gerektirir. Yazarlar бул süreci kolaylaştırmak үчүн atomiongpu.m adlı бир üst seviye kullanıcı script'i geliştirmiştir.

Kullanıcı dosyanın üst bölümündeki parametreleri değiştirir; script баштапкы koşullarını oluşturur, paralelleştirmeyi düzenler, ode45gpu'yu çalıştırır, sonuçları toplar, gözlenebilirleri işler жана dosyaları üretir.

atomiongpu.m менен кайсы атом жана iyonlar değiştirilebilir?

Atom жана ион kütleleri atomik birimlerle kullanıcı tarafından girilebilir. Örnek ayarlarda:

• ¹⁷¹Yb⁺ iyonu,
• ⁸⁷Rb atomu

колдонулган.

Atom жана ион sıcaklıkları da Kelvin cinsinden ayrı ayrı belirlenebilir.

“head-on” жана “thermal” çarpışma modları arasındaki fark nedir?

head-on: Atom sabit hızla doğrudan orijine/iyona doğru gönderilir. İyon başlangıçta orijinde hareketsizdir. Бул modda ион sıcaklığı değişkeni kullanılmaz.

thermal: Atom жана ион ылдамдык bileşenleri sıcaklıklarına karşılık gelen normal dağılımlardan örneklenir. İyonun tuzak içindeki баштапкы koordinatları da sıcaklığa жана seküler tuzak frekansına bağlı normal dağılımdan seçilir. Kaynak бул modu atomların oluşturduğu termal banyodaki ион үчүн көбүрөөк gerçekçi senaryo катары tanımlamaktadır.

Potansiyel үчүн кайсы giriş biçimleri destekleniyor?

Bütün durumlarda kullanılan genel потенциал:

\[ V(r)=-\frac{C_n}{r^n}+\frac{C_m}{r^m} \]

biçiminde kalır.

İyon tuzaksız, harmonik tuzakta же Paul tuzağında иштетилиши мүмкүн mi?

Evet. Kaynak üç yapı tarif etmektedir.

Tuzaksız система

Tuzak parametreleri sıfıra ayarlanır.

Harmonik tuzak

\(q_x=q_y=q_z=0\) seçilir жана uygun \(a_x,a_y,a_z\) parametreleri harmonik frekanslara göre girilir.

Paul tuzağı

Kullanıcı \(a_x,a_y,a_z,q_x,q_y,q_z\) жана \(\Omega_{\mathrm{rf}}\) değerlerini gerçek Paul tuzağı parametrelerine göre verir.

Simülasyon эмне zaman duruyor?

Her орбита iki durma koşulundan hangisi мурун gerçekleşirse o anda sonlandırılır:

• симуляция zamanı tmax'a ulaştığında,
• же zaman adımı sayısı maxsteps'e ulaştığında.

Kaynak örneğinde tmax yaklaşık 1 ms жана maxsteps 1.000.000'dur.

Kaç орбита иштетилиши мүмкүн?

Kaynağa göre ntrajectories yaklaşık 10.000.000 düzeyine kadar ayarlanabilir. Yazarlar yaklaşık 25.000.000 yörüngeyi denediklerinde bellek sorunları gözlemlediklerini belirtmektedir.

Buradaki sınır evrensel эмес; kullanılan GPU belleği, çıktı seçenekleri жана diğer donanım koşullarına bağlıdır.

ntheta жана nphi эмне işe yarıyor?

Yörünge başlangıçları iki boyutlu бир \(n_\theta\times n_\phi\) matris yapısında organize жасалат. nanglesauto="yes" olduğunda program бул iki değeri орбита sayısının kareköküne yakın biçimde otomatik seçer.

Manuel modda:

\[ n_\theta n_\phi=n_{\mathrm{trajectories}} \]

koşulu uygulanır.

CPU жана GPU тандоосу кантип жүргүзүлөт?

processor değişkeni:

• CPU же
• GPU

катары ayarlanabilir.

ncores, kullanılacak bağımsız CPU süreçlerini belirler. GPU modunda her süreç kendi GPU'sunda çalıştırılabildiğinden kaynak, ncores değerinin kullanılabilir GPU sayısını aşmaması gerektiğini belirtmektedir.

Tüm yörüngeleri saklamıyorsa belirli бир yörüngeyi кантип çiziyor?

Program мурун төмөн bellekli toplu hesaplamayı yapar. Daha кийин kullanıcı isterse belirli бир баштапкы koşulunu yeniden çalıştırarak tam yörüngeyi kaydeder.

İki өзгөчө seçenek vardır:

• longestlived: en uzun kompleks ömrünü veren баштапкы koşulunu bulur жана tam yörüngeyi tekrar hesaplar.
• custom: kullanıcının girdiği 12 bileşenli өзгөчө баштапкы koşulunu tam орбита катары hesaplar.

Tam орбита dört biçimde görselleştirilir:

• атом–ион uzaklığına karşı zaman,
• xy düzlemi орбита izdüşümü,
• xz izdüşümü,
• yz izdüşümü.

Hangi gözlenebilirler kaydedilebiliyor?

Sonuç dosyaları кантип yönetiliyor?

Kullanıcı üç temel depolama davranışını kontrol edebilir:

• MATLAB workspace'in work.mat катары kaydedilmesi,
• бардык gözlenebilirlerin CSV dosyalarının kaydedilmesi,
• гана seçilen tek бир gözlenebilirin CSV dosyasının kaydedilmesi.

Program мындан тышкары çalışmaya başlamadan мурун yaklaşık depolama gereksinimini kullanıcıya bildirir. Böylece kullanıcının disk alanı yetersizse симуляция başlamadan çıktı seçeneklerini azaltması mümkündür.

Isı haritası neyi көрсөтөт?

Kaynak, her gözlenebiliri atomun баштапкы küresel açıları \(\theta\) жана \(\phi\)'nin fonksiyonu катары ısı haritasına dönüştürebilmektedir. Saçılma açısı үчүн verilen örnekte bazı баштапкы bölgeleri basit geri saçılma üretirken başka bölgeler son derece karmaşık örüntüler oluşturmaktadır.

Yazarlar бул desenlerde görülen fraktal yapıyı sistemin kaotik saçılma dinamikleriyle ilişkilendirmektedir.

Örnek çalıştırmada эмне oldu?

Makaledeki örnek konsol çıktısında program 10.004.569 атом–ион yörüngesi үчүн эсептөө başlatmaktadır. Hesap sekiz GPU колдонуу менен yürütülmüş жана tamamlanma süresi 52.848,874 saniye катары raporlanmıştır; изилдөө bunu 10 milyondan fazla yörüngenin 15 saatin астында tamamlanması катары özetlemektedir.

Aynı örnek çalıştırmada programın hesapladığı kompleks oluşma olasılığı %59,87 катары yazdırılmıştır. Бул değer изилдөөнүн бардык атом–ион sistemleri үчүн evrensel kompleks olasılığı эмес; гана gösterilen parametreler жана баштапкы koşulu örneklemesi үчүн elde edilen симуляция çıktısıdır.

CPU жана GPU эмне үчүн бирдей kaotik yörüngede ар түрдүү sonuç verebilir?

Çalışmanın MATLAB açısından en öğretici bölümlerinden biri budur. Aynı баштапкы koşulları жана бирдей ODE algoritması kullanılmasına rağmen GPU менен CPU'da hesaplanan çok kaotik бир yörüngenin bounce sayısı жана kompleks ömrü ар түрдүү çıkabilmiştir.

Yazarlar bunu CPU жана GPU mimarilerinin kayan noktalı işlemleri ар түрдүү yuvarlama sıralarıyla gerçekleştirebilmesine bağlamaktadır. Kaotik sistemlerde баштапкы же sayısal hesaplamadaki son derece küçük farklar zamanla büyüyebilir.

Örnek konsol çıktısındaki en uzun ömürlü yörüngede GPU hesabı 3 bounce жана 987,1 µs, CPU üzerindeki ode45gpu hesabı болсо 2 bounce жана 954,9 µs vermektedir.

Бул абал tek başına GPU hesabının yanlış экендиги anlamına gelmez. Program doğrulama amacıyla бирдей өзгөчө yörüngeyi CPU үстүндө hem ode45gpu hem MATLAB'ın стандарттык ode45 çözücüsüyle yeniden çalıştırmakta жана son durumların uyuşup uyuşmadığını kontrol etmektedir.

Бул bize sayısal симуляция hakkında эмне öğretiyor?

Бир симуляция kodunun бирдей diferansiyel denklemi çözmesi, her donanımda her ara noktada bit düzeyinde бирдей yörüngeyi üretmesini garanti etmez. Özellikle kaotik sistemlerde:

• kayan noktalı yuvarlama,
• işlem sıralaması,
• CPU жана GPU mimari farkları

gözlenebilir yörüngeyi zaman ичинде ayırabilir. Бул nedenle изилдөө doğrulamayı гана tek бир kaotik yörüngenin mikro ayrıntıları üzerinden эмес, чечүүчү tutarlılığı жана istatistiksel çıktılar üzerinden ele almaktadır.

MATLAB öğrenen biri бул çalışmadan кайсы программалоо сабактарын чыгара алат?

Birinci ders: MATLAB'da performans гана algoritmanın matematiksel karmaşıklığı эмес; dizi oluşturma, indeksleme жана bellek davranışı da hesap süresini ciddi biçimde etkileyebilir.

İkinci ders: GPU'ya geçmeden мурун CPU kodunu optimize etmek маанилүү. Kaynakta гана күч fonksiyonunun yeniden yazılması yaklaşık üç kat ылдамдык kazandırmıştır.

Üçüncü ders: GPU'ya uygun kod her zaman CPU үчүн yazılan en doğal MATLAB koduyla бирдей biçimde эмес. Vektörler burada GPU arrayfun gereksinimleri nedeniyle skaler değişkenlere dönüştürülmüştür.

Dördüncü ders: Milyonlarca simülasyonda бардык ara verileri сактоо yerine гана илимий катары gerekli özet değişkenleri hesaplamak çok büyük bellek avantajı sağlayabilir.

Beşinci ders: Paralel işlem бирок bağımsız iş birimleri varsa verimlidir. Бул çalışmada her атом–ион yörüngesi diğerlerinden bağımsız экендиги үчүн GPU paralelliği doğal biçimde uygulanabilmektedir.

Altıncı ders: Sayısal doğrulama зарыл. Hız kazanımı tek başına yeterli эмес; sonuçların güvenilirliği стандарттык çözücüyle karşılaştırılmalıdır.

Изилдөө колдогон натыйжалар

• Atom–ион klasik yörüngeleri MATLAB үстүндө bağımsız баштапкы koşulları üzerinden kitlesel biçimde paralelleştirilebilir.
• Kaynak v1 benchmark'ında күч fonksiyonunun yeniden yazılması tek başına yaklaşık üç kat CPU hızlanması sağlamıştır.
• Yüklenen preprint benchmark'ında özelleştirilmiş ode45gpu, belirli referans ode45 kullanımına göre tek CPU çekirdeğinde 22 kata kadar hızlıdır.
• GPU үстүндө yaklaşık 10.000 yörüngeye kadar güçlü eşzamanlılık bölgesi gözlenmiştir.
• atomiongpu.m, yaklaşık on milyon орбита ölçeğine kadar giriş hazırlama, paralel hesap, istatistik, CSV жана görselleştirme işlerini otomatikleştirebilir.
• Saçılma açısı жана benzeri gözlenebilirlerin баштапкы açılarına göre ısı haritaları kaotik жана fraktal yapılar gösterebilir.

Изилдөө колдобогон же сынаб көрбөгөн натыйжалар

• Бул изилдөө kuantum mekanik атом–ион dinamiğinin tam çözümü эмес; klasik hareket denklemlerini kullanmaktadır.
• ode45gpu'nun бардык ODE problemlerinde MATLAB ode45'ten 22 kat hızlı экендиги gösterilmemiştir.
• GPU жана CPU'nun her kaotik yörüngede birebir бирдей zaman geçmişini üretmesi beklenmemektedir.
• Yaklaşık on milyon орбита sınırı бардык bilgisayarlar үчүн sabit бир kapasite эмес; bellek жана donanıma bağlıdır.
• Çalışma laboratuvar deneyi gerçekleştirmemiş жана симуляция sonuçlarını yeni бир deneysel veri setiyle doğrudan doğrulamamıştır.
• Yüklenen v1'deki 22× benchmark değeri көбүрөөк sonraki hakemli yayın sürümünün nihai performans sonucu катары kullanılmamalıdır.

Изилдөөнүн ыкмасы жана натыйжалары

Simülasyon türü

Araştırma klasik moleküler dinamik simülasyonudur. Atom жана ион noktasal klasik parçacıklar катары hareket eder; iç serbestlik dereceleri hesaba katılmaz. Düşük yoğunluk varsayımı nedeniyle tek бир iyonun бирдей anda бир atomla etkileştiği iki parçacıklı dinamik çözülür.

Benchmark fiziksel parametreleri

Başlangıç koşulları

\[ \mathbf y_0= ( 0,0,0,0,0,0, r_0\sin\theta\cos\phi, r_0\sin\theta\sin\phi, r_0\cos\theta, -v_0\sin\theta\cos\phi, -v_0\sin\theta\sin\phi, -v_0\cos\theta ). \]

İlk altı bileşen benchmark'ta iyonun orijinde жана hareketsiz olduğunu; sonraki üç bileşen атом konumunu, son üç bileşen болсо atomun merkeze yönelmiş баштапкы hızını verir.

CPU күч fonksiyonu optimizasyonu

ode45gpu CPU benchmark'ı

2500 yörüngede, yüklenen v1'e göre ode45gpu tek CPU çekirdeğinde referans ode45 + f(t,y) çalıştırmasından 22× hızlıdır. Paralellik 28 CPU çekirdeğine kadar artırıldığında kaynak, бул ылдамдык oranının yaklaşık korunduğunu bildirmektedir.

Бул karşılaştırmanın preprint sürümüne özgü экендиги өзгөчө korunmalıdır; 2026'daki hakemli sürüm benchmark metodolojisini жана performans rakamlarını değiştirmiştir.

GPU ölçekleme deneyi

NVIDIA Tesla K80 24 GB, P100 16 GB жана V100 32 GB үстүндө:

\[ n=1,2,5,10,20,\ldots,1.000.000 \]

орбита test edilmiştir.

Yaklaşık \(n\sim10^4\) yörüngeye kadar изилдөө süresinin büyük ölçüde sabit olması GPU paralelliğinin henüz doygunluğa ulaşmadığını; көбүрөөк кийин sürenin yaklaşık \(n\) менен orantılı büyümesi işlem kapasitesinin dolduğunu көрсөтөт.

Gerçekçi баштапкы koşulu testi

K80 үстүндө бирдей yörüngenin kopyalarıyla yapılan benchmark'a ek катары, ilk oktant üzerindeki ар түрдүү \((\theta,\phi)\) баштапкы koşullarıyla da test yapılmıştır. Küçük \(n\) değerlerinde орбита karmaşıklığı nedeniyle изилдөө süresinde көбүрөөк fazla dalgalanma görülse de жогорку орбита sayısında iki test de yaklaşık doğrusal ölçekleme göstermiştir.

atomiongpu.m örnek çalıştırma sonucu

Бул tablo tek бир örnek parametre kümesinin çıktısıdır. %59,87 değeri genel атом–ион kompleks oluşum oranı катары kullanılmamalıdır.

GPU–CPU айырмасы башаламан орбитада

Kaynak бул farklılığı kayan noktalı hesapların GPU жана CPU mimarilerinde ар түрдүү yuvarlama davranışlarına sahip olması жана dinamiğin kaotik olmasıyla açıklamaktadır. Kod, мындан тышкары CPU üzerindeki ode45gpu çıktısını MATLAB'ın ode45 sonucu менен karşılaştırarak son koşulların eşleşmesini doğrulamaktadır.

Programın илимий чыгыштары

atomiongpu.m toplu simülasyonun ardından:

• kompleks oluşum olasılığını,
• saçılma açısı dağılımını,
• kompleks ömrünü,
• bounce sayısını,
• momentum transferini,
• iyonun son kinetik enerjisini,
• zaman adımı sayısını,
• son атом–ион uzaklığını

hesaplayabilir жана бул değişkenleri баштапкы \(\theta,\phi\) açılarına göre ısı haritaları halinde gösterebilir.

Yazılımın күчтүү жактары

• MATLAB ortamından ayrılmadan GPU paralelliği kullanabilmesi.
• Çok sayıda bağımsız орбита үчүн төмөн bellekli hesap yaklaşımı.
• Tuzaksız, harmonik жана Paul tuzak yapılarını desteklemesi.
• Atom/ион kütleleri жана sıcaklıkların değiştirilebilmesi.
• Farklı genelleştirilmiş атом–ион potansiyellerinin girilebilmesi.
• Head-on жана thermal баштапкы koşulları.
• CPU жана GPU seçeneklerinin бирдей kullanıcı script'inde bulunması.
• CSV, heatmap, өзгөчө орбита жана outlier çıktılarının otomatik oluşturulması.
• Standart ode45 менен doğrulama kontrolü.

Yöntemsel sınırlılıklar

• Simülasyon klasik dinamiktir; iç kuantum serbestlik dereceleri yoktur.
• İki бөлүкчө varsayımı төмөн yoğunluk koşuluna dayanır.
• ode45gpu genel amaçlı MATLAB ODE API'sinin birebir GPU eşleniği эмес; күч fonksiyonunun өзгөчө biçimde gömülmesi gerekir.
• GPU arrayfun kısıtları kod mimarisini şekillendirmiştir.
• Kaotik dinamiklerde CPU жана GPU yörüngeleri küçük kayan nokta farklarından dolayı ayrışabilir.
• Yaklaşık on milyon орбита kapasitesi donanıma жана belleğe bağlıdır.
• Kaynak, yeni бир laboratuvar атом–ион deneyini raporlamamaktadır.
• Yüklenen preprint'in benchmark sayıları nihai hakemli sürümde değiştirilmiştir.

Çalışmanın MATLAB açısından asıl önemi

Бул изилдөө MATLAB'ın гана “бир denklem yaz, ode45 çağır жана grafik çiz” biçiminde kullanılmadığını gösteren güçlü бир örnektir. Tek бир fizik problemi üzerinden:

1. fiziksel Hamiltonyen kurulmakta,
2. 12 boyutlu ODE sistemi çıkarılmakta,
3. стандарттык чечүүчү benchmark жүргүзүлөт,
4. MATLAB kodunun veri yapıları optimize жүргүзүлөт,
5. Runge–Kutta algoritması өзгөчө катары yeniden uygulanmaktadır,
6. GPU yürütmesine uygun skaler kod oluşturulmaktadır,
7. milyonlarca баштапкы koşulu gpuArray менен hazırlanmakta,
8. arrayfun менен bağımsız орбиталар paralel çözülmekte,
9. gather менен sonuçlar CPU'ya alınmakta,
10. жана son катары илимий gözlenebilirler, CSV'ler, ısı haritaları жана өзгөчө орбиталар üretilmektedir.

Dolayısıyla изилдөө MATLAB'da илимий hesaplamanın dört temel katmanını бирдей örnekte бир araya getirir: fiziksel modelleme, sayısal analiz, performans mühendisliği жана veri görselleştirme.

Булак жана ыкма боюнча эскертүү

Tam özgün изилдөө adı: GPU-Accelerated MATLAB Software for Atom-Ion Dynamics

Yazarlar: Saajid Chowdhury; Jesús Pérez-Ríos.

Yazar sırası: Kaynaktaki sıra aynen korunmuştur.

Eş katkı/eş birinci yazar: Kaynakta bildirilmemiştir.

Sorumlu yazar: Saajid Chowdhury.

Kurum: Department of Physics and Astronomy, Stony Brook University, Stony Brook, New York 11794, USA.

Yüklenen kaynak türü: Hesaplamalı fizik / илимий yazılım preprint'i жана klasik moleküler dinamik симуляция çalışması.

Preprint platformu: arXiv.

Preprint kimliği: arXiv:2509.12381v1.

Preprint tarihi: 15 Eylül 2025.

Hakemlik durumu: Бул изилдөө, yüklenen sürüm itibarıyla hakem değerlendirmesinden geçmemiş бир preprinttir; sonuçları бул yayın aşaması dikkate alınarak değerlendirilmelidir.

Yüklenen sürümde DOI: Bulunmamaktadır.

Daha sonraki hakemli yayın: Aynı изилдөө көбүрөөк кийин “GPU-parallelized MATLAB software for атом-ion dynamics” başlığıyla Computer Physics Communications, Cilt 321, Makale 110041, 2026'da yayımlanmıştır.

Hakemli sürüm DOI: 10.1016/j.cpc.2026.110041.

Önemli sürüm farkı: Yüklenen arXiv v1'de CPU benchmark'ı үчүн 22× hızlanma raporlanırken, nihai hakemli yayın CPU hızlanmasını yaklaşık 7× катары raporlamakta жана мындан тышкары бир milyon yörüngenin GPU paralelleştirmesinde CPU ode45 tabanına göre 600–3500× hızlanma bildirmektedir. Бул Verianla makalesinin илимий gövdesindeki 22×, 10.000-орбита eşiği жана örnek 8-GPU çalıştırma değerleri гана yüklenen v1'den alınmıştır; iki sürüm sessizce birleştirilmemiştir.

Program: atomiongpu.m.

Sayısal чечүүчү: ode45gpu.

Programlama ortamı: MATLAB R2023a жана Parallel Computing Toolbox.

Geliştirici deposu: Kaynakta saajidchowdhury/supplementGPU катары belirtilmiştir.

Lisans: Yüklenen preprint'in program özetinde “Licensing provisions” alanı doldurulmamıştır. Nihai hakemli sürüm program үчүн CC0 1.0 belirtmektedir.

Finansman: United States Air Force Office of Scientific Research, FA9550-23-1-0202.

Veri/kod erişimi: Kaynak, kodların, şekillerin жана yardımcı materyalin geliştirici deposunda bulunduğunu belirtmektedir.

Yazar katkıları: Yüklenen v1'de ayrı CRediT katkı tablosu bulunmamaktadır. Daha sonraki hakemli sürümde Saajid Chowdhury үчүн yazılım, görselleştirme жана yazım; Jesús Pérez-Ríos үчүн kavramsallaştırma, metodoloji, danışmanlık жана finansman edinimi dahil katkılar raporlanmıştır.

Çıkar çatışması: Yüklenen v1'de ayrı бир çıkar çatışması bölümü bulunmamaktadır; hakemli yayın sürümünde мындан тышкары beyan edilmiştir.

Bilimsel içerik kaynağı: Бул Verianla makalesindeki fizik modeli, algoritma, MATLAB mimarisi, giriş parametreleri, benchmark değerleri, GPU testleri, konsol чыгыштары жана илимий yorumlar kullanıcının yüklediği arXiv v1 çalışmasına dayanmaktadır. Daha sonraki yayın sürümü гана bibliyografik doğrulama жана sürüm farkının açıkça belirtilmesi amacıyla колдонулган.

Temel yöntemsel sınır: Бул бир klasik moleküler dinamik simülasyonudur. Çalışma yeni бир saha же laboratuvar deneyi sunmaz; kuantum dinamiğinin tam çözümü эмес жана yazılım performansı problem, kod yapısı жана donanım bağımlıdır.