DESTEK OL

LAMMPS

LAMMPS: Atom ve Moleküler Dinamik Simülasyonlarının Gücü — Adım Adım Kullanım Kılavuzu

Bu rehber, LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator) yazılımının atom ve molekül düzeyinde malzemelerin, polimerlerin, nano yapıların ve kimyasal sistemlerin dinamik davranışlarını nasıl simüle edeceğinizi, analiz edeceğinizi ve sonuçları nasıl yorumlayacağınızı adım adım açıklar. Özellikle malzeme bilimi, biyofizik ve kimya mühendisliği alanlarında, faz geçişleri, difüzyon, mekanik özellikler ve reaksiyon kinetiğini incelemek için vazgeçilmez bir araçtır.

💡 Önemli Uyarı: LAMMPS, grafik arayüzü olmayan, komut satırı ve komut dosyaları (input scripts) ile çalışan bir simülasyon motorudur. Kullanımı için temel Linux/Unix komut satırı bilgisi ve metin editörü kullanımı gereklidir. Görselleştirme için VMD veya Ovito gibi programlarla birlikte kullanılmalıdır.

1. LAMMPS Nedir?

LAMMPS, Sandia Ulusal Laboratuvarları tarafından geliştirilen, atom ve moleküler düzeyde sistemlerin dinamiklerini simüle etmek için kullanılan açık kaynaklı bir yazılımdır. Paralel hesaplama mimarilerinde (CPU ve GPU) yüksek performansla çalışacak şekilde optimize edilmiştir. “Potansiyel alan” (force field) tabanlı hesaplamalar yapar.

Temel Özellikler:

  • Atom, molekül, metal, polimer, biyomolekül ve granüler sistemlerin simülasyonu.
  • Çeşitli entegre potansiyel alanlar (Lennard-Jones, Morse, REAXFF, CHARMM, AMBER, AIREBO vb.).
  • NPT, NVT, mikrokanonik (NVE) gibi farklı termodinamik topluluklar.
  • Esnek, sert ve sınırlı moleküler yapılar.
  • Yüksek performanslı paralel hesaplama (MPI ve GPU desteği).
  • Çıktı olarak atom konumları, enerjiler, basınç, sıcaklık ve özel hesaplamalar.

2. Adım Adım Kullanım Kılavuzu

Adım 1: Yazılımı İndirme ve Kurma

LAMMPS tamamen ücretsiz ve açık kaynaklıdır. Resmi GitHub deposundan veya web sitesinden indirebilirsiniz:

Kurulum, işletim sisteminize göre değişir. Linux için genellikle kaynak koddan derlenir. Windows kullanıcıları için önceden derlenmiş .exe dosyaları veya WSL (Windows Subsystem for Linux) üzerinden kurulum önerilir.

Temel Kurulum (Ubuntu):

sudo apt update sudo apt install lammps

Adım 2: İlk Simülasyon — Argon Sıvısının MD Simülasyonu

1. Bir metin editörü (Notepad++, VS Code, nano, vim) ile in.argon adında bir dosya oluşturun.
2. Aşağıdaki komutları dosyaya yapıştırın:

# LAMMPS Input Script: Argon MD Simulation units lj atom_style atomic lattice fcc 0.8442 region box block 0 10 0 10 0 10 create_box 1 box create_atoms 1 box mass 1 1.0 velocity all create 1.0 87287 loop geom pair_style lj/cut 2.5 pair_coeff 1 1 1.0 1.0 2.5 neighbor 0.3 bin neigh_modify delay 0 every 1 check yes fix 1 all nvt temp 1.0 1.0 0.1 thermo 100 thermo_style custom step temp press etotal run 1000

3. Terminalde aşağıdaki komutla simülasyonu çalıştırın:

lmp -in in.argon

Açıklamalar:

  • units lj: Lennard-Jones birim sistemi.
  • pair_style lj/cut: Lennard-Jones potansiyeli.
  • fix nvt: Sıcaklık sabit tutulur (NVT topluluğu).
  • run 1000: 1000 zaman adımı boyunca simülasyon yapar.

Adım 3: Polimer Zincirinin Simülasyonu

1. in.polymer adlı yeni bir dosya oluşturun.
2. Aşağıdaki komutlarla polimer zinciri oluşturun ve simüle edin:

units real atom_style full bond_style harmonic angle_style harmonic special_bonds lj 0.0 1.0 1.0 read_data polymer.data velocity all create 300.0 12345 fix 1 all nvt temp 300.0 300.0 100.0 dump 1 all atom 100 traj.lammpstrj run 10000

Not: Bu simülasyon için önceden hazırlanmış bir polymer.data dosyası gerekir. Bu dosya, atom koordinatlarını, bağları ve açıları içerir. Hazır örnekler LAMMPS web sitesinde mevcuttur.

Adım 4: Sonuçların Görselleştirilmesi ve Analizi

1. Simülasyon çıktısı olarak traj.lammpstrj dosyası üretilir.
2. Bu dosyayı görselleştirmek için VMD veya Ovito programlarını kullanın.
3. VMD’de: File → New Molecule → Browse → traj.lammpstrj → Determine file type: LAMMPS Trajectory
4. Animasyonu oynatın ve atom hareketlerini inceleyin.
5. Enerji, sıcaklık ve basınç verilerini log.lammps dosyasından çıkarıp Excel veya Python ile grafik çizebilirsiniz.

3. Pratik Uygulamalar ve Örnek Senaryolar

Senaryo 1: Su Moleküllerinin Difüzyon Katsayısının Hesaplanması

1. TIP3P su modelini kullanarak bir su kutusu oluşturun.
2. NVT topluluğunda 300K’de simülasyon yapın.
3. MSD (Mean Squared Displacement) hesaplamak için compute msd komutunu kullanın.
4. MSD grafiğinin eğiminden difüzyon katsayısını hesaplayın.

Senaryo 2: Nanotüp Üzerinde Adsorpsiyon Simülasyonu

1. Karbon nanotüp geometrisini oluşturun veya içe aktarın.
2. Nanotüp üzerine gaz molekülleri (örneğin, CO₂) yerleştirin.
3. NPT topluluğunda simülasyon yaparak adsorpsiyon dinamiklerini inceleyin.
4. Potansiyel enerji değişimlerini ve molekül yoğunluk haritalarını analiz edin.

4. Yaygın Sorunlar ve Çözümleri

  • “Unknown command” hatası: Komut yanlış yazılmış veya ilgili paket kurulu değil. make yes-package ile gerekli paketleri derleyin.
  • “Out of range atoms - cannot compute PPPM”: Atomlar simülasyon kutusunun dışına çıkmış. fix shake veya fix recenter kullanın.
  • Simülasyon kararsız: Zaman adımı (timestep) çok büyük olabilir. timestep 0.001 gibi küçük bir değer deneyin.

5. Bilgiyi Test Et

Soru 1: LAMMPS hangi tür simülasyonlar için kullanılır?

Cevap: Atomistik ve Moleküler Dinamik (MD) simülasyonları.

Soru 2: LAMMPS’ta bir simülasyonu çalıştırmak için hangi komut kullanılır?

Cevap: Terminalde lmp -in input_dosyasi.in

Soru 3: LAMMPS çıktılarını görselleştirmek için hangi programlar kullanılır?

Cevap: VMD veya Ovito.