Verlet entegrasyonu

Verlet entegrasyonu, Newton'un hareket denklemlerini uygulamak için kullanılan nümerik yöntemlerden biridir. Genellikle Moleküler dinamik simülasyonlarında parçacıkların bir sonraki zaman dilimindeki konumlarını belirlemek için kullanılır. Hız hesaplaması yerine sadece o anki konum, önceki konum ve o anki ivmeyi kullanan bu yöntem Euler yönteminden daha isabetlidir ve gerektirdiği işlem sayısı pek farklı değildir. İlk defa 1791 yılında Delambre tarafından kullanılmıştır ve o zamandan beri çok kez yeniden keşfedilmiştir: 1909'da Cowell and Crommelin tarafından Halley kuyruklu yıldızı'nın yörüngesini hesaplamak için veya 1907'de Carl Størmer tarafından manyetik alandaki elektrik yüklü parçacıkların yörüngesini incelemek için kullanılması gibi (ayrıca bu yüzden Störmer yöntemi de denir). Daha sonra 1960'larda Loup Verlet tarafından moleküler dinamikte kullanıldı.

 Temel Størmer-Verlet entegrasyonu 

İkinci dereceden bir diferansiyel denklem olan  denklemi, başlangıç koşulları  ve  ile birlikte,  zaman aralığında yaklaşık olarak, ,  aşağıdaki yöntemle elde edilebilir:
 Birinci eleman 'in değeri için: ,
 Serinin diğer elemanları için, n = 1, 2, ...  olmak şartıyla: 

 Hareket Denklemleri 
Enerjinin korunduğu fiziksel sistemler için Newton'un hareket denklemi:

ya da ayrı ayrı her bir kütle için:

olarak ifade edilebilir. Burada:
  zamanı ifade eder,
 ,  sayıdaki cisimlerin konum vektörlerinin kümesidir,
  skaler potansiyel fonksiyonudur,
 , parçacıklar üzerindeki kuvvetlerin kümesini veren Potansiyel gradyanının negatifidir,
 , her bir parçacık için  kütlelerine sahip bloklarla genellikle diyagonal olan kütle matrisidir.

Bu denklem, potansiyel fonksiyonu 'nin farklı  seçimleri için, etkileşen moleküllerin hareketinden gezegenlerin yörüngelerine kadar çeşitli fiziksel sistemlerin hareketini tanımlamak için kullanılabilir.

Kütleyi denklemin sağ tarafına almak ve birden fazla parçacığın etkileşimli yapısını ihmal etmek için yapılan dönüşümlerden sonra, denklem şu şekilde basitleştirilebilir: 

Burada , konuma bağlı ivmeyi temsil eden sisteme uygun bir vektör değerli fonksiyondur. Genellikle, başlangıç konumu  ve başlangıç hızı  ifadeleri ile verilir.

Verlet integrasyonu (hızlar olmadan)
Bu başlangıç değer problemini sayılabilir aralıklara bölerek  çözmek için, bir zaman adımı  seçilir ve zaman dizisi ,  olarak ele alınır. Amaç, kesin çözümün yörüngesi üzerindeki  noktalarını yakından takip eden bir  noktalar dizisi oluşturmaktır.

Euler yöntemi, birinci dereceden diferansiyel denklemlerde birinci türev için ileri fark yaklaşımını kullanırken, Verlet integrasyonu, ikinci türev için merkezi fark yaklaşımını kullanır:

Verlet integrasyonu, Størmer yöntemi olarak kullanılan formda,web sayfası  Størmer yöntemi açıklamasıyla. hızı kullanmadan, yalnızca önceki iki konum vektörünü kullanarak bir sonraki konum vektörünü elde etmek için şu denklemi kullanır: 

Burada çözüme ulaşmak için  terimi yalnız bırakılmıştır.

Hızların Hesaplanması – Størmer–Verlet yöntemi

Hızlar, temel Størmer denklemlerinde açıkça verilmemiştir, ancak genellikle kinetik enerji gibi bazı fiziksel niceliklerin hesaplanması için gereklidir. Bu, moleküler dinamik simülasyonlarında teknik zorluklar yaratabilir, çünkü bir sistemdeki kinetik enerji ve anlık sıcaklıklar,  zamanındaki konumlar bilinmeden hesaplanamaz. Bu eksiklik, ya velocity Verlet algoritması kullanılarak ya da konum terimleri ve ortalama değer teoremi kullanılarak hızın tahmin edilmesiyle giderilebilir.

Bu hız terimi, konum teriminden bir adım geridedir, çünkü bu  anındaki hıza değil,  anındaki hıza karşılık gelir. Yani,  ifadesi,  için ikinci dereceden bir yaklaşımdır.

 Verlet Hız Algoritması 

Sistem ile ilgili ve daha yaygın olarak kullanılan bir diğer algoritma ise Verlet hız algoritmasıdır. Bu yöntem, hız hesabı için yukarıda verilen yöntem ile benzer bir yaklaşıma sahiptir fakat hızı açıkça dahil ederek, temel Verlet algoritmasındaki adım problemini çözer:

 Kaynakça 

Kategori:Sayısal diferansiyel denklemler
Kategori:Hesaplamalı fizik
Kategori:Moleküler dinamik