Ana Sayfa
Fizik

Klasik Mekanik Serisi Bölüm 1

Bu bölümde klasik mekaniğin geometrik temelleri; afin uzaylar, konfigürasyon manifoldları ve en az eylem ilkesini anlattık

Eren Sivri20 Mayıs 202622 dk okuma süresi
Klasik Mekanik Serisi Bölüm 1

Bölüm 1: Newton'un Kafesinden Kaçış: Konfigürasyon Uzayı ve Manifoldlar

Teorik fiziğin derinliklerine inildikçe, birinci sınıf amfilerindeki standart eğik düzlem problemlerinin yerini, doğanın temel simetrileri ve bu simetrilerin dokuduğu geometrik yapılar alır. Bu noktada klasik mekaniğin asıl zarafeti, Newton'un oklarla çizilmiş o kısıtlayıcı kafesinden (kartezyen vektörlerden) çıkıp, evrenin geometrik dilini keşfetmektir.


1 Koordinatlardan Kurtulmak: Afin Uzay Olarak Evren

Fiziğin en temel ilkesi şudur: Evrende ayrıcalıklı bir nokta yoktur. Yani evrenin bir "orijini" (0,0,0)(0,0,0) noktası yoktur. Eğer uzayı doğrudan bir vektör uzayı olan R3\mathbb{R}^3 olarak tanımlarsak, bu uzayın doğasında olan sıfır vektörüne (orijine) haksız bir kutsallık atfetmiş oluruz. Arnold bu felsefi problemi, uzayı bir Afin Uzay (Affine Space) olarak tanımlayarak çözer.

Afin uzayı, "başlangıç noktası unutulmuş bir vektör uzayı" olarak düşünebilirsiniz. Noktalar ve vektörler iki farklı hayvandır. Uzayımız olan A3\mathbb{A}^3'te sadece "noktalar" (örneğin A,B,CA, B, C) vardır. İki noktayı birbirinden çıkardığımızda bir vektör elde edersiniz, ancak iki noktayı toplayamazsınız (iki konumu toplamanın fiziksel bir anlamı yoktur!).

Matematiksel olarak uzayımız, bir R3\mathbb{R}^3 vektör uzayının etki ettiği bir A3\mathbb{A}^3 afin uzayıdır. Yani bir AA noktasına bir v\mathbf{v} vektörü eklerseniz, yeni bir BB noktası elde edersiniz:

B=A+v    v=BA(1)B = A + \mathbf{v} \implies \mathbf{v} = B - A \tag{1}

İşte bu basit tanım, uzayın homojenliğini (her yerin aynı kurallara tabi olmasını) doğrudan geometriye kodlar. Orijin yoktur, sadece noktalar arası ilişkiler vardır. Uzay, yönsüz ve hiyerarşisiz boş bir sayfadır.


İki Uzay Arasındaki Köprü: Orijini Seçmek ve Unutmak

Peki bu yönsüz, hiyerarşisiz A3\mathbb{A}^3 afin uzayından, bildiğimiz o hesaplanabilir R3\mathbb{R}^3 vektör uzayına (veya tam tersine) nasıl geçeriz? Cevap çok basittir ama fizikte derin bir anlam taşır: Bir orijin (referans noktası) seçerek.

Afin uzayı "orijinini unutmuş bir vektör uzayı" olarak tanımlamıştık. Eğer bu boş sayfa olan A3\mathbb{A}^3 içinde tamamen keyfi bir OO noktası seçer ve buna "orijin" dersek, afin uzaydaki herhangi bir PP noktası, OO'dan PP'ye çizilen bir konum vektörü ile (r=PO)(\mathbf{r} = P - O) birebir eşleşir. İşte o an, noktalar vektörlere dönüşür ve afin uzayımız bir R3\mathbb{R}^3 vektör uzayına çöker. Artık elimizde standart kartezyen eksenleri oturtabileceğimiz bir temel vardır.

Tersine, R3\mathbb{R}^3 vektör uzayından A3\mathbb{A}^3 afin uzayına geçmek istiyorsak, orijini "unuturuz". Vektörlerin mutlak konumlarını siler, uzayı sadece nesnelerin birbirine göre olan uzaklıkları ve yönleriyle (fark vektörleriyle) tanımlarız.
Arnold'un mekanik vizyonundaki felsefe şudur: Evrenin gerçek dokusu A3\mathbb{A}^3'tür; evrenin kendisi bir orijine sahip değildir. Ancak biz ölümlü fizikçiler denklemleri yazabilmek ve ölçüm yapabilmek için uzaya bir gözlemci (referans sistemi) yerleştiririz ve uzayı geçici olarak R3\mathbb{R}^3'e indirgeriz. Klasik mekaniğin asıl gücü, bulduğumuz fiziksel yasaların, bizim sonradan uzaya çaktığımız bu keyfi OO noktasına asla bağlı olmamasından gelir. Peki bu bağımsızlık (veya değişmezlik) uzay-zamanda kendini nasıl gösterir?


2 Mutlak Zaman ve Galilei Uzay-Zamanı

Sadece uzay yetmez, olayların ne zaman gerçekleştiğini de bilmeliyiz. Newtonyen (veya Galilei) mekaniğinde zaman, uzaydan tamamen bağımsız, tek boyutlu bir afin uzaydır: A1\mathbb{A}^1. Zamanın da bir başlangıcı (orijini) yoktur; "M.Ö. 10.000" veya "Big Bang'den sonra 3 saniye" demek, zamanı bir vektör uzayı yapmaz, sadece referans noktaları seçmektir.

Peki, uzay-zaman (Spacetime) dediğimiz şey nedir? Arnold uzay-zamanı, noktaları "olaylar" (events) olan 4 boyutlu bir A4\mathbb{A}^4 afin uzayı olarak tanımlar. Ancak bu herhangi bir 4 boyutlu uzay değildir. Bu uzayın içinde "zaman" fonksiyonu adını verdiğimiz çok özel bir izdüşüm (projection) vardır:

t:A4A1(2)t : \mathbb{A}^4 \to \mathbb{A}^1 \tag{2}

Uzay-zamandaki iki olay aa ve bb olsun. Bu olaylar arasındaki zaman farkı t(b)t(a)t(b) - t(a) ile verilir.

Burada Arnold'un klasik mekanik inşasındaki en kritik kavrama geliyoruz: Eşzamanlılık (Simultaneity). İki olayın aynı anda gerçekleşmesi demek, t(a)=t(b)t(a) = t(b) olması demektir. Galilei evreninde uzaklık, sadece ve sadece eşzamanlı olaylar arasında tanımlıdır. Yani "Şu an benim bulunduğum yer ile, 10 saniye sonra Mars'ta olacak olan patlama arasındaki mesafe nedir?" sorusu klasik mekanikte anlamsızdır. Çünkü hareket eden bir gözlemciye göre uzaysal mesafe değişir. Ancak iki olay aynı anda gerçekleşiyorsa, evrendeki tüm gözlemciler bu iki olay arasındaki uzaysal mesafenin (Öklidyen metrik) aynı olduğu konusunda hemfikirdir.

İşte Arnold tarzı uzay-zaman: Eşzamanlılık dilimlerinden (her biri birer A3\mathbb{A}^3) oluşan bir krep yığını gibi düşünebilirsiniz. İleride Görelilik Teorisi'ne veya Kuantum Alan Teorisi'ne geçtiğinizde, Einstein'ın bu "krep yığınını" nasıl büküp erittiğini (mutlak zamanı yıkıp, mesafeyi ve zamanı tek bir Minkowski metriğinde nasıl birleştirdiğini) çok daha iyi takdir edeceksiniz.


3 Galilei Grubu: Fiziğin Kurallarını Belirleyen Simetri

Sahneyi kurduk. Peki bu sahnede izin verilen "kamera hareketleri" nelerdir? Fizik yasalarının değişmediği referans sistemi dönüşümlerine Galilei Grubu denir. Bu grup şunları içerir:

  • Uzaysal Ötelemeler: rr+a\mathbf{r} \to \mathbf{r} + \mathbf{a} (Uzayın homojenliği)
  • Zamansal Ötelemeler: tt+τt \to t + \tau (Zamanın homojenliği)
  • Uzaysal Döndürmeler: rRr\mathbf{r} \to R\mathbf{r} (Uzayın izotropisi, yönden bağımsızlığı)
  • Eylemsiz (Inertial) Dönüşümler (Galilei Boost): rr+vt\mathbf{r} \to \mathbf{r} + \mathbf{v}t (Sabit hızla hareket eden sistemlerin denkliği)

4 Eylemsizlik Referans Sistemleri ve Vektörlerin Yetersizliği

Mekanik olayları incelemek için öncelikle bir referans sistemi seçilmelidir; zira keyfi bir sistem seçildiğinde uzay homojen ve izotropik olmaz. Ancak Galileo'nun görelilik ilkesine (Galileo's relativity principle) göre, uzayın homojen (her noktasının eşdeğer) ve izotropik (her yönünün eşdeğer), zamanın ise homojen olduğu "eylemsizlik sistemleri" (inertial frames) bulmak daima mümkündür. Eylemsiz bir sistemde serbest bir cisim, büyüklüğü ve yönü sabit bir hızla hareket eder.

Ancak elimizdeki sistem serbest parçacıklardan değil de bağlara (constraints) sahip cisimlerden oluştuğunda ne olur? Newtonyen mekanik bu noktada bizi zorlu bir sürece sokar. Bir küre yüzeyine hapsedilmiş bir sarkacın hareketini R3\mathbb{R}^3 uzayındaki (x,y,z)(x, y, z) vektörleriyle tanımlamaya çalışırsak, cismin yüzeyden ayrılmasını engelleyen ve yönü büyüklüğü sürekli değişen bağ kuvvetlerini (constraint forces) hesaba katmak zorunda kalırız. Oysa bu tepki kuvvetleri fiziksel bir etkileşimden ziyade, sistemin geometrik sınırlarının bir sonucudur. İdeal holonomik bağlara (holonomic constraints) sahip sistemlerin analizi için, sistemi tanımlayan denklemleri bu görünmez ve işe yaramaz kuvvetlerin yükünden kurtarmamız gerekir.


5 Neden R3\mathbb{R}^3 Uzayını Terk Edip Toruslara ve Kürelere Sığınıyoruz?

Newton bize bir başlangıç noktası verdi ama bizi düz, sıkıcı ve kartezyen bir R3\mathbb{R}^3 uzayına hapsetti. Oysa evrensel kısıtlamalar altında (örneğin küre yüzeyinde hareket eden bir sarkaç veya bir ray üzerindeki trende) kartezyen koordinatlar çaresiz kalır. Denklemler, bağ kuvvetleri (constraint forces) ile dolar taşar. Bağ kuvvetleri denmesinin nedeni de kavramsal olarak bu kuvvetler, sistemi konfigürasyon manifolduna bağlar ve sistemin o manifoldun (yüzeyin veya eğrinin) dışına çıkmasını geometrik olarak kısıtlar.

Arnold, "Mathematical Methods of Classical Mechanics" kitabında bu kısıtlar karşısında olaya getirilen yeniliği şöyle belirtir: Fiziksel sistemin yaşayacağı alanı, sistemin kendi serbestlik dereceleriyle (degrees of freedom) yeniden tanımlamak.

Eğer sistemimiz NN parçacığa sahipse ve kk tane holonomik yani integral edilebilir kısıtlaması varsa, sistemin özgürce hareket edeceği alan n=3Nkn = 3N - k olur. İşte bu nn adet bağımsız koordinatla (q1,q2,,qn)(q_1, q_2, \ldots, q_n) tanımladığımız soyut uzaya Konfigürasyon Manifoldu (MM) diyoruz. Bir çiftli sarkaç düşünün; uzayda hareket eden iki kütle olmasına rağmen, sistemi sadece iki açıyla (θ1,θ2)(\theta_1, \theta_2) tamamen tanımlayabilirsiniz. Yani bu çiftli sarkacın konfigürasyon uzayı aslında bir Torus'tur (T2=S1×S1T^2 = S^1 \times S^1). Arnold'un geometrik diliyle konuşursak, sistemin her bir anlık konumu bu manifold üzerinde bir noktadır.

Bu manifold üzerindeki bir nokta qMq \in M sistemin o anki tam durumunu (tüm parçacıkların konumunu) tek başına ifade eder. Zaman aktıkça (t)(t), bu nokta manifold üzerinde pürüzsüz bir eğri çizer:

γ:RM,tq(t)(3)\gamma : \mathbb{R} \to M, \quad t \mapsto q(t) \tag{3}

Peki ya hız? Hız dediğimiz şey, manifold üzerindeki bu eğriye teğet olan vektördür. Ancak bu vektör manifoldun içinde yaşamaz; manifoldun o noktasındaki Teğet Uzay'da (Tangent Space, TqMT_qM) yaşar. Bir noktanın teğet uzayı, o noktadan geçebilecek tüm olası hız vektörlerinin kümesidir. Konumları ve hızları aynı anda ele almak istersek, manifoldun tüm noktalarındaki teğet uzayları birleştirip yepyeni, 2n2n boyutlu devasa bir uzay inşa ederiz: Teğet Demeti (Tangent Bundle, TMTM).

TM=qMTqM(4)TM = \bigcup_{q \in M} T_q M \tag{4}

İşte Lagrange fonksiyonu dediğimiz o efsanevi LL, ne gökten inmiş bir formüldür ne de rastgele bir skalerdir. LL, bu Teğet Demeti üzerinden ve zamandan gerçel sayılara giden pürüzsüz (diferansiyellenebilir) bir fonksiyondur:

L:TM×RRL : TM \times \mathbb{R} \to \mathbb{R} (q,q˙,t)L(q,q˙,t)(5)(q, \dot{q}, t) \mapsto L(q, \dot{q}, t) \tag{5}

Yani sisteme konumunu (q)(q), hızını (q˙)(\dot{q}) ve zamanı (t)(t) verirsin, o sana skaler bir sayı verir. Bütün fizik bu skaler fonksiyonun içindedir.


6 Örnek: Helisel Rollercoaster Sarkacı

Bu soyut geometrik yapıyı tamamen içselleştirmek için Newtonyen kuvvetlerin içinden çıkılamaz bir kaos yaratacağı, ancak analitik mekaniğin bir sanat eserine dönüştüreceği bir sistem tasarlayalım: Helisel Rollercoaster Sarkacı.

Sistemimiz şu şekilde olsun: Uzayda sabitlenmiş sarmal (helis) şeklinde bir telimiz var. m1m_1 kütleli bir boncuk bu tel üzerinde sürtünmesizce kayabiliyor. Bu m1m_1 kütlesine ise ll uzunluğunda kütlesiz bir çubukla m2m_2 kütleli ikinci bir parçacık bağlanmış (sarkaç). Sarkacın sadece radial-dikey (kütle merkezinden dışarı doğru) düzlemde sallanmasına izin verelim.

Koordinatların (Manifoldun) Tanımlanması:

Sistemimizin iki serbestlik derecesi vardır. Manifoldumuz, helisin tek boyutlu yapısı (S1S^1 benzer bir topoloji) ve sarkacın açısal hareketi (S1S^1) ile oluşur. Koordinatları tanımlayalım:

  • q1q_1: Helis üzerindeki m1m_1 kütlesinin konumunu belirleyen azimut açısı.
  • q2q_2: Sarkacın düşeyle yaptığı açı.

Helisin parametreleri yarıçap RR ve adım aralığı sabiti c=h/2πc = h/2\pi olsun. m1m_1 kütlesinin uzaydaki konumu:

r1=(Rcosq1, Rsinq1, cq1)(6)\mathbf{r}_1 = (R\cos q_1,\ R\sin q_1,\ cq_1) \tag{6}

m2m_2 kütleli sarkacın konumu, m1m_1 bulunduğu noktadan dışarıya doğru radyal olarak lsinq2l\sin q_2 kadar, aşağı doğru (Z ekseninde) ise lcosq2l\cos q_2 kadar yer değiştirmiştir:

r2=((R+lsinq2)cosq1, (R+lsinq2)sinq1, cq1lcosq2)(7)\mathbf{r}_2 = \big((R + l\sin q_2)\cos q_1,\ (R + l\sin q_2)\sin q_1,\ cq_1 - l\cos q_2\big) \tag{7}

Hızlar ve Kinetik Enerji:

Kinetik enerji, genelleştirilmiş hızların her zaman karesel (quadratic) bir fonksiyonudur ve T=12aik(q)q˙iq˙kT = \frac{1}{2}\sum a_{ik}(q)\dot{q}_i\dot{q}_k şeklinde ifade edilir. Koordinatların zamana göre türevini alarak hızların karelerini buluruz.

r˙12=(R2+c2)q˙12(8)\dot{\mathbf{r}}_1^2 = (R^2 + c^2)\dot{q}_1^2 \tag{8}

r˙22\dot{\mathbf{r}}_2^2 terimi için türev alıp kareleri topladığımızda, manifoldun eğriliğinden kaynaklanan çapraz terimler (cross terms) ortaya çıkar:

r˙22=q˙12[(R+lsinq2)2+c2]+l2q˙22+2clsinq2q˙1q˙2(9)\dot{\mathbf{r}}_2^2 = \dot{q}_1^2\left[(R + l\sin q_2)^2 + c^2\right] + l^2\dot{q}_2^2 + 2cl\sin q_2\,\dot{q}_1\dot{q}_2 \tag{9}

Toplam kinetik enerji T=12m1r˙12+12m2r˙22T = \frac{1}{2}m_1\dot{\mathbf{r}}_1^2 + \frac{1}{2}m_2\dot{\mathbf{r}}_2^2 düzenlendiğinde şu forma ulaşır:

T=12[m1(R2+c2)+m2((R+lsinq2)2+c2)]a11q˙12+12[m2l2]a22q˙22+m2clsinq2a12=a21q˙1q˙2(10)T = \frac{1}{2}\underbrace{\left[m_1(R^2+c^2) + m_2\left((R+l\sin q_2)^2+c^2\right)\right]}_{a_{11}}\dot{q}_1^2 + \frac{1}{2}\underbrace{\left[m_2 l^2\right]}_{a_{22}}\dot{q}_2^2 + \underbrace{m_2 cl\sin q_2}_{a_{12}=a_{21}}\dot{q}_1\dot{q}_2\tag{10}

Metrik Tensör (Riemannian Metric):

Arnold'a göre, eylemsiz hareket eden bir sistemin konfigürasyon uzayındaki her teğet uzayda (tangent space) tanımlı pozitif tanımlı karesel form, bir Riemann metriği (ds2)(ds^2) yaratır. Bu metrik doğrudan kinetik enerjinin katsayı matrisidir:

gij=(m1(R2+c2)+m2((R+lsinq2)2+c2)m2clsinq2m2clsinq2m2l2)(11)g_{ij} = \begin{pmatrix} m_1(R^2+c^2) + m_2\left((R+l\sin q_2)^2+c^2\right) & m_2 cl\sin q_2 \\ m_2 cl\sin q_2 & m_2 l^2 \end{pmatrix} \tag{11}

İşte R3\mathbb{R}^3'ten kaçışın zaferi budur. O karmaşık tepki kuvvetleriyle boğuşmak yerine, sistemin yaşadığı iki boyutlu manifoldun "metrik tensörünü" elde ettik. Sarkaç sallandıkça (q2q_2 değiştikçe), metrik tensörün g11g_{11} ve g12g_{12} bileşenleri değişir; yani manifoldun eğriliği bizzat sistemin konfigürasyonu tarafından şekillendirilir. Sistemin serbest hareketleri, artık bu eğri uzaydaki jeodeziklerden (geodesics) başka bir şey değildir.

Ancak burada çok zarif bir detayı gözden kaçırmamalıyız: "Serbest hareket" tanımı, sistemin üzerinde hiçbir fiziksel potansiyel enerjinin (V(q)=0V(q) = 0) olmadığı durumlar için geçerlidir. Kurguladığımız sistemde, manifoldun z-eksenindeki helisel eğiminden ve sarkacın düşey düzlemdeki salınımından dolayı işin içine kaçınılmaz olarak kütleçekim potansiyel enerjisi girer. Sistemin mutlak z-koordinatlarını (z1=cq1z_1 = cq_1 ve z2=cq1lcosq2z_2 = cq_1 - l\cos q_2) dikkate aldığımızda sistemin gerçek potansiyel enerjisi şu şekildedir:

V(q)=m1g(cq1)+m2g(cq1lcosq2)=(m1+m2)gcq1m2glcosq2(12)V(q) = m_1 g(cq_1) + m_2 g(cq_1 - l\cos q_2) = (m_1 + m_2)gcq_1 - m_2 gl\cos q_2 \tag{12}

Potansiyel enerjinin varlığında, parçacığın yörüngesi artık sadece kütle dağılımından doğan o saf gijg_{ij} metriğinin jeodeziklerine itaat etmez; kütleçekimi, konfigürasyon manifoldunun geometrisini kendi etrafında bükmeye başlar. Bu kuvvet altındaki hareketi yeniden pürüzsüz bir "jeodezik problemi" olarak geometrik bir çerçeveye oturtmak istersek, kinetik enerjiden elde ettiğimiz metrik tensörü, sistemin toplam enerjisini (EE) ve potansiyelini barındıran Maupertuis-Jacobi Metriği ile yeniden tanımlamamız gerekir:

ds2=2(EV(q))aik(q)dqidqk(13)ds^2 = 2(E - V(q))\sum a_{ik}(q)\,dq_i\,dq_k \tag{13}

Böylece klasik mekanik, potansiyel enerjiyi bile kuvvet vektörleri yerine "uzayın eğriliğine" dönüştürerek geometrikleştirir. Dinamiğin gerçek kalbine inmek için artık tek yapmamız gereken, teğet demetindeki bu hızları (q˙)(\dot{q}) momentumlara (p)(p) çevirip Faz Uzayı'na (TMT^*M) adım atmaktır.


7 Doğanın Tembelliği: En Az Eylem İlkesi (Principle of Least Action)

Manifoldumuzu kurduk, kinetik enerjinin yarattığı Riemann metriğini tanımladık ve potansiyelin uzayı nasıl büktüğünü gördük. Peki ama bir parçacık AA noktasından BB noktasına giderken, bu karmaşık geometri içinde hangi yolu seçeceğine nasıl karar verir? Parçacıklar ilerideki engelleri "görüp" ona göre mi rota çizerler, yoksa anlık kuvvetlerin körü körüne itmesiyle mi hareket ederler?

İşte tüm klasik mekaniği, kuantum mekaniğini ve genel göreliliği tek bir çatı altında birleştiren o keskin cevap: En Az Eylem İlkesi'dir.

Teğet demeti (TMTM) üzerinde tanımladığımız L(q,q˙,t)L(q, \dot{q}, t) Lagrange fonksiyonunun, sistemin belirli bir t1t_1 anından t2t_2 anına kadar uzanan bir yörünge (γ)(\gamma) boyunca integrali alındığında ortaya "Eylem" (Action, SS) adını verdiğimiz bir skaler çıkar:

S(γ)=t1t2L(q(t),q˙(t),t)dt(14)S(\gamma) = \int_{t_1}^{t_2} L(q(t), \dot{q}(t), t)\,dt \tag{14}

Bu bir fonksiyon değil, bir fonksiyoneldir (fonksiyonların fonksiyonu). Sistemin konfigürasyon uzayında çizebileceği sonsuz sayıda olası yörünge vardır. Ancak doğa son derece tutumlu, hatta "tembeldir". Evren, parçacığın izleyeceği gerçek yörüngeyi öyle bir seçer ki, bu SS eylem integrali bir ekstremum (genellikle bir minimum) değerini alır. Yani sistemin hareketi, varyasyonu sıfır olan (δS=0\delta S = 0) rotadır. Parçacık sanki tüm olası yolları aynı anda "koklar" ve eylemi minimize eden o eşsiz jeodeziği seçip süzülür.

Varyasyon Nedir? Yörüngelerin Sanal Dansı

Burada "varyasyonu sıfır olan (δS=0\delta S = 0) rotadır" dedik, ancak bu δ\delta (varyasyon) operatörünün fiziksel ve matematiksel anlamını netleştirmek zorundayız.

Normal diferansiyel hesapta (d/dt)(d/dt), zaman akarken bir parçacığın konumundaki değişimi inceleriz. Varyasyonlar hesabında (δ)(\delta) ise durum bambaşkadır: Zamanı dondurur ve sistemin konfigürasyon uzayında çizdiği o pürüzsüz yörüngeyi bütünüyle ele alıp, onu hayali olarak bükeriz.

Diyelim ki doğanın seçtiği o kusursuz, gerçek yörünge q(t)q(t) olsun. Bu yörüngeye çok küçük, sanal bir sapma (virtual displacement) ekleyerek komşu bir yörünge hayal edelim:

q~(t)=q(t)+δq(t)(15)\tilde{q}(t) = q(t) + \delta q(t) \tag{15}

Buradaki δq(t)\delta q(t) fonksiyonuna, yörüngenin varyasyonu denir. Ancak doğanın bu sanal denemelere izin verirken koyduğu çok katı bir kural vardır: Sistemin nereden yola çıktığı (t1)(t_1) ve nereye varacağı (t2)(t_2) kesindir. Uç noktalar duvara çivilenmiştir, yani varyasyonlar başlangıç ve bitiş anlarında sıfır olmak zorundadır:

δq(t1)=0,δq(t2)=0(16)\delta q(t_1) = 0, \quad \delta q(t_2) = 0 \tag{16}

İşte δS=0\delta S = 0 (Eylemin durağanlığı) ilkesinin bize söylediği derin gerçek budur: Uç noktaları sabit tutarak gerçek yörünge üzerinde yapacağınız her türlü sonsuz küçük hayali bükülme (varyasyon δq\delta q), eylem integralinde (SS) birinci mertebeden hiçbir değişime yol açmaz. Nasıl ki standart kalkülüste bir fonksiyonun minimum noktasında teğetin eğimi sıfırsa (df=0df = 0), konfigürasyon manifoldunun üzerindeki o devasa fonksiyonel uzayda da doğanın seçtiği yörüngenin varyasyonu sıfırdır. Parçacık, eylemin topografyasında adeta düzleştiği o ekstremum noktasını bulup oradan süzülür.


8 Neden L=TVL = T - V? Simetrilerin Katı Dayatması

Analitik mekanik ders kitaplarında veya teorik fizik amfilerindeki ilk derslerde genellikle tepeden inme bir kabulle başlarız: Lagrange fonksiyonu kinetik enerji eksi potansiyel enerjidir (L=TVL = T - V). Peki ama neden? Evren neden L=T+VL = T + V veya L=T2/VL = T^2/V gibi bir formülü değil de bu spesifik yapıyı seçmiştir?

İşte Lev Landau'nun mekanik inşasındaki en büyük dehası burada yatar. Landau, L=TVL = T - V denkleminin gökten zembille inmediğini, yazının başında tanımladığımız Galilei Grubu simetrilerinin matematiksel bir zorunluluğu olduğunu ispatlar.

Bunu matematiksel olarak adım adım görelim. Eylemsiz bir referans sisteminde, tamamen serbest (üzerinde hiçbir kuvvet olmayan) bir parçacık hayal edelim. Bu parçacığı tanımlayacak olan Lagrange fonksiyonu başlangıçta konumun (r)(\mathbf{r}), hızın (v=r˙)(\mathbf{v} = \dot{\mathbf{r}}) ve zamanın (t)(t) rastgele bir fonksiyonu olabilir: L(r,v,t)L(\mathbf{r}, \mathbf{v}, t).

Fakat doğanın simetrileri bu fonksiyona anında kısıtlamalar getirmeye başlar:

  • Uzayın Homojenliği: Evrende ayrıcalıklı bir nokta yoktur. Parçacığı uzayda Δr\Delta\mathbf{r} kadar kaydırmak fiziksel hiçbir yasayı değiştirmez. Bu yüzden LL, konum vektörüne (r)(\mathbf{r}) bağlı olamaz: Lr=0\frac{\partial L}{\partial \mathbf{r}} = 0.

  • Zamanın Homojenliği: Fizik kuralları dün neyse bugün de odur. Zamanın mutlak bir başlangıcı olmadığı için LL, zamana açıkça (explicitly) bağlı olamaz: Lt=0\frac{\partial L}{\partial t} = 0.

  • Uzayın İzotropisi: Uzayda ayrıcalıklı bir yön yoktur. Dolayısıyla LL, hız vektörünün (v)(\mathbf{v}) yönüne değil, sadece onun skaler büyüklüğüne, daha doğrusu büyüklüğünün karesine (v2=vv)(v^2 = \mathbf{v} \cdot \mathbf{v}) bağlı olabilir.

Bütün bu geometrik kısıtlamalar sonucunda, serbest parçacığın Lagrangian'ı oldukça sadeleşir: L=L(v2)L = L(v^2).

Peki fonksiyonun tam formu nedir? Bunu bulmak için Galilei Görelilik İlkesi'ni kullanmalıyız. Kurallar, birbirine göre sabit bir ϵ\boldsymbol{\epsilon} hızıyla hareket eden iki eylemsiz sistemde (KK ve KK') tamamen aynı olmalıdır. KK' sistemindeki hız v=v+ϵ\mathbf{v}' = \mathbf{v} + \boldsymbol{\epsilon} olur. Bu çok küçük (infinitesimal) bir ϵ\boldsymbol{\epsilon} hızı için yeni Lagrangian'ı Taylor serisine açarsak:

L(v2)=L ⁣((v+ϵ)2)=L(v2+2vϵ+ϵ2)L(v2)+Lv2(2vϵ)(17)L(v'^2) = L\!\left((\mathbf{v}+\boldsymbol{\epsilon})^2\right) = L(v^2 + 2\mathbf{v}\cdot\boldsymbol{\epsilon} + \epsilon^2) \approx L(v^2) + \frac{\partial L}{\partial v^2}(2\mathbf{v}\cdot\boldsymbol{\epsilon}) \tag{17}

İki referans sisteminin fiziksel olarak tamamen aynı hareket denklemlerini vermesi için, bu iki Lagrangian arasındaki farkın, konum ve zamanın tam bir zamansal türevi (dFdt)\left(\frac{dF}{dt}\right) olması şarttır. Fark terimine yakından bakalım:

δL=2Lv2vϵ=2Lv2drdtϵ(18)\delta L = 2\frac{\partial L}{\partial v^2}\mathbf{v}\cdot\boldsymbol{\epsilon} = 2\frac{\partial L}{\partial v^2}\frac{d\mathbf{r}}{dt}\cdot\boldsymbol{\epsilon} \tag{18}

Bu ifadenin saf bir ddt()\frac{d}{dt}(\ldots) formunda yazılabilmesi ancak ve ancak Lv2\frac{\partial L}{\partial v^2} teriminin hıza bağlı olmayan bir sabit olmasıyla mümkündür. Bu sabiti 12m\frac{1}{2}m olarak adlandırırsak (buradaki mm parçacığın kütlesidir), türevi geri integre ettiğimizde muazzam bir kesinlikle şu sonuca varırız:

L(v2)=12mv2=T(19)L(v^2) = \frac{1}{2}mv^2 = T \tag{19}

Gördüğümüz gibi kinetik enerjinin formu (v2v^2'ye orantılı olması), tamamen uzayın izotropisi ve Galilei dönüşümlerinin bir diktesidir.

Eğer uzayda parçacığın serbestliğini bozan, uzayın homojenliğini kıran bir etkileşim alanı varsa (örneğin kütleçekimi), bu alan koordinatlara bağlı skaler bir fonksiyon olarak eyleme katılmalıdır. Doğanın basitlik ilkesi gereği, bu etkileşim terimi sadece konuma bağlı bir V(r)-V(\mathbf{r}) eklentisi olarak Lagrangian'a girer. İşte Galilei görelilik ilkesi ve evrenin homojenliği, bizi o kaçınılmaz formüle böyle ulaştırır:

L=12mv2V(r)=TV(20)L = \frac{1}{2}mv^2 - V(\mathbf{r}) = T - V \tag{20}

9 Geometriden Dinamiğe: Euler-Lagrange Denklemlerinin Türetilişi

Eylemi minimize etme felsefesini (δS=0\delta S = 0) ve sanal yer değiştirmeleri (δq)(\delta q) tanımladık. Şimdi bu geometrik fikri katı bir matematiğe dökelim. Eylem integralinin varyasyonunu alırsak:

δS=δt1t2L(q,q˙,t)dt=t1t2(Lqδq+Lq˙δq˙)dt=0(21)\delta S = \delta \int_{t_1}^{t_2} L(q, \dot{q}, t)\,dt = \int_{t_1}^{t_2} \left(\frac{\partial L}{\partial q}\delta q + \frac{\partial L}{\partial \dot{q}}\delta\dot{q}\right)dt = 0 \tag{21}

Burada kritik bir adım atmamız gerekiyor. δq˙\delta\dot{q} terimi, aslında varyasyonun zamana göre türevidir, yani δq˙=ddt(δq)\delta\dot{q} = \frac{d}{dt}(\delta q). Bunu integralin içine yerleştirip, ikinci terime kısmi integrasyon (integration by parts) uygularsak:

t1t2Lq˙ddt(δq)dt=[Lq˙δq]t1t2t1t2ddt ⁣(Lq˙)δqdt(22)\int_{t_1}^{t_2} \frac{\partial L}{\partial \dot{q}}\frac{d}{dt}(\delta q)\,dt = \left[\frac{\partial L}{\partial \dot{q}}\delta q\right]_{t_1}^{t_2} - \int_{t_1}^{t_2} \frac{d}{dt}\!\left(\frac{\partial L}{\partial \dot{q}}\right)\delta q\,dt \tag{22}

Hatırlayalım, doğa sanal varyasyonlara izin verirken uç noktaları duvara çivilemişti! Yani t1t_1 ve t2t_2 anlarında varyasyon sıfırdır (δq(t1)=0\delta q(t_1) = 0 ve δq(t2)=0\delta q(t_2) = 0). Bu katı kural sayesinde, sağ taraftaki sınır değer terimi (köşeli parantez içindeki kısım) tamamen sıfırlanarak yok olur. Geriye kalan integrali orijinal denklemde yerine koyarsak:

δS=t1t2[Lqddt ⁣(Lq˙)]δqdt=0(23)\delta S = \int_{t_1}^{t_2} \left[\frac{\partial L}{\partial q} - \frac{d}{dt}\!\left(\frac{\partial L}{\partial \dot{q}}\right)\right]\delta q\,dt = 0 \tag{23}

İşte varyasyonlar hesabının temel leması (Fundamental Lemma of Calculus of Variations) burada devreye girer. t1t_1 ve t2t_2 arasında seçtiğimiz sanal sapma (δq)(\delta q) tamamen keyfidir. Eğer keyfi bir δq\delta q fonksiyonu ile çarpılan bir integral her zaman sıfır sonucunu veriyorsa, parantez içindeki ifadenin kendisi zorunlu olarak her yerde sıfır olmalıdır! Böylece konfigürasyon manifoldunun koordinatlarından tamamen bağımsız olan o meşhur diferansiyel denklemler ortaya çıkar:

ddt ⁣(Lq˙i)Lqi=0(i=1,2,,n)(24)\frac{d}{dt}\!\left(\frac{\partial L}{\partial \dot{q}_i}\right) - \frac{\partial L}{\partial q_i} = 0 \quad (i = 1, 2, \ldots, n) \tag{24}

Newton mekaniğinde (F=maF = ma) kartezyen koordinatlardan küresel veya silindirik koordinatlara geçmek bir kabustur; ivme terimleri akıl almaz çapraz çarpımlarla dolar. Oysa bu türettiğimiz Euler-Lagrange denklemleri kovaryanttır. Yani uzayımızı tanımlamak için kullandığımız qiq_i koordinatları ister kartezyen olsun, ister açı, ister bizim rollercoaster örneğimizdeki gibi bir helis uzunluğu olsun, bu diferansiyel denklemin formu asla değişmez. Konfigürasyon uzayındaki metrik ne kadar eğri büğrü olursa olsun, bu denklem o uzaydaki doğru yörüngeyi kusursuzca bulur.


10 Sistemin Uyanışı: Hareket Denklemlerinin Çözümü

Artık her şeye sahibiz. 6. bölümde inşa ettiğimiz o karmaşık Helisel Rollercoaster Sarkacı sistemini "uyandırma" zamanı geldi. Kinetik enerjimizi (T)(T) ve potansiyel enerjimizi (V)(V) bulmuştuk. Sistemin devasa Lagrangian'ı şu şekildedir:

L=12[m1(R2+c2)+m2((R+lsinq2)2+c2)]q˙12+12m2l2q˙22+m2clsinq2q˙1q˙2V(q1,q2)(25)L = \frac{1}{2}\left[m_1(R^2+c^2) + m_2((R+l\sin q_2)^2+c^2)\right]\dot{q}_1^2 + \frac{1}{2}m_2 l^2\dot{q}_2^2 + m_2 cl\sin q_2\,\dot{q}_1\dot{q}_2 - V(q_1, q_2) \tag{25}

Bu sistem için hareket denklemlerini bulmak, karmaşık serbest cisim diyagramları veya 3 boyutlu vektör izdüşümleri çizmeyi gerektirmez. Tek yapmamız gereken bu LL fonksiyonunu Euler-Lagrange makinesine atmaktır. Örneğin, sarkacın düşey salınımını veren q2q_2 koordinatını ele alalım:

ddt ⁣(Lq˙2)Lq2=0(26)\frac{d}{dt}\!\left(\frac{\partial L}{\partial \dot{q}_2}\right) - \frac{\partial L}{\partial q_2} = 0 \tag{26}

Önce hızlara göre kısmi türevi (genelleştirilmiş momentumu) bulalım ve zamana göre türevini alalım:

Lq˙2=m2l2q˙2+m2clsinq2q˙1    ddt()=m2l2q¨2+m2clcosq2q˙1q˙2+m2clsinq2q¨1(27)\frac{\partial L}{\partial \dot{q}_2} = m_2 l^2\dot{q}_2 + m_2 cl\sin q_2\,\dot{q}_1 \implies \frac{d}{dt}(\ldots) = m_2 l^2\ddot{q}_2 + m_2 cl\cos q_2\,\dot{q}_1\dot{q}_2 + m_2 cl\sin q_2\,\ddot{q}_1 \tag{27}

Şimdi de konuma (q2)(q_2) göre kısmi türevi (genelleştirilmiş kuvveti) alalım. Potansiyel enerjiden gelen türevin Vq2=m2glsinq2-\frac{\partial V}{\partial q_2} = -m_2 gl\sin q_2 olduğuna dikkat edin:

Lq2=m2l(R+lsinq2)cosq2q˙12+m2clcosq2q˙1q˙2m2glsinq2(28)\frac{\partial L}{\partial q_2} = m_2 l(R + l\sin q_2)\cos q_2\,\dot{q}_1^2 + m_2 cl\cos q_2\,\dot{q}_1\dot{q}_2 - m_2 gl\sin q_2 \tag{28}

Bu iki ifadeyi birbirinden çıkardığımızda o harika matematiksel "sadeleşme" gerçekleşir; q˙1q˙2\dot{q}_1\dot{q}_2'li çapraz terimler birbirini yok eder ve geriye saf hareket denklemi kalır:

m2l2q¨2+m2clsinq2q¨1m2l(R+lsinq2)cosq2q˙12+m2glsinq2=0(29)m_2 l^2\ddot{q}_2 + m_2 cl\sin q_2\,\ddot{q}_1 - m_2 l(R + l\sin q_2)\cos q_2\,\dot{q}_1^2 + m_2 gl\sin q_2 = 0 \tag{29}

Bakın, Newton'un kuvvet diyagramlarında asla bu kadar kolay göremeyeceğimiz her şey bir çırpıda karşımızda: Sarkacın eylemsizlik momenti (m2l2q¨2)(m_2 l^2\ddot{q}_2), helisteki ivmelenmeden kaynaklı etkileşim kuplajı (m2clsinq2q¨1)(m_2 cl\sin q_2\,\ddot{q}_1), sistem dönerken sarkacı dışarı savuran merkezkaç terimi (m2l(R+lsinq2)cosq2q˙12)(m_2 l(R + l\sin q_2)\cos q_2\,\dot{q}_1^2) ve kütleçekiminin o ince izdüşümü (m2glsinq2)(m_2 gl\sin q_2). Hepsi salt geometriden doğdu!

Geometrik Bir Çıkarım: Vektör Alanları ve Faz Akışı

İşte Arnold, elde ettiğimiz bu ivme (q¨)(\ddot{q}) denklemlerine sadece "çözülecek bir cebir problemi" olarak bakmaz. Bu diferansiyel denklemler dizisi, aslında 2n2n-boyutlu Teğet Demeti (TM)(TM) üzerinde tanımlı devasa bir vektör alanıdır (vector field).

Sistemin t1t_1 anından t2t_2 anına evrimi, anlık durum olan (q,q˙)(q, \dot{q}) noktasının, bu vektör alanının akıntılarına kapılarak bir sıvı damlası gibi sürüklenmesidir. Arnold bu duruma Faz Akışı (Phase Flow) adını verir. Doğanın temel yasaları, uzayı bu akışın integral eğrileri ile baştan aşağı dokur.


11 Enerjinin Doğuşu ve Hamilton'a İlk Bakış: Jacobi İntegrali

Peki bu vektör alanında sürüklenen sistemde değişmeden kalan hiçbir şey yok mudur? İşte bu noktada Arnold, zamanın homojenliğini kullanarak enerjinin korunumunu o gökten inme formüllerle değil, bizzat Euler-Lagrange üzerinden ispatlar.

Eğer sistemde dış koşullar zamanla değişmiyorsa (Lt=0)\left(\frac{\partial L}{\partial t} = 0\right), Lagrangian'ın toplam zaman türevini zincir kuralıyla açalım:

dLdt=iLqiq˙i+iLq˙iq¨i(30)\frac{dL}{dt} = \sum_i \frac{\partial L}{\partial q_i}\dot{q}_i + \sum_i \frac{\partial L}{\partial \dot{q}_i}\ddot{q}_i \tag{30}

Euler-Lagrange denklemi gereği Lqi\frac{\partial L}{\partial q_i} yerine ddt ⁣(Lq˙i)\frac{d}{dt}\!\left(\frac{\partial L}{\partial \dot{q}_i}\right) yazarsak:

dLdt=i[ddt ⁣(Lq˙i)q˙i+Lq˙iq¨i]=ddt ⁣(iLq˙iq˙i)(31)\frac{dL}{dt} = \sum_i \left[\frac{d}{dt}\!\left(\frac{\partial L}{\partial \dot{q}_i}\right)\dot{q}_i + \frac{\partial L}{\partial \dot{q}_i}\ddot{q}_i\right] = \frac{d}{dt}\!\left(\sum_i \frac{\partial L}{\partial \dot{q}_i}\dot{q}_i\right) \tag{31}

Bütün türevleri tek bir tarafa topladığımızda karşımıza muazzam bir matematiksel sabitlik, yani Jacobi İntegrali çıkar:

ddt ⁣(i=1nq˙iLq˙iL)=0    E=i=1nq˙iLq˙iL=Sabit(32)\frac{d}{dt}\!\left(\sum_{i=1}^{n} \dot{q}_i\frac{\partial L}{\partial \dot{q}_i} - L\right) = 0 \implies E = \sum_{i=1}^{n} \dot{q}_i\frac{\partial L}{\partial \dot{q}_i} - L = \text{Sabit} \tag{32}

Doğa, zamanın homojenliğini kullanarak işte bu devasa EE fonksiyonunu korur. Ancak bu korunum yasasının içine gizlenmiş matematiksel form (q˙pL)(\sum \dot{q}p - L), basit bir enerjiden çok daha fazlasıdır. Bu form, bizi hızların teğet dünyasından koparıp mekaniğin asıl kalbinin attığı yer olan Kotanjant Demeti'ne (TM)(T^*M), diğer bir deyişle Faz Uzayı'na taşıyacak olan gizli bir köprüdür.

Bu köprünün adını koymak ve bizi o yepyeni evrene nasıl sıçrattığını görmek için, bir sonraki bölümde Hamilton'un kusursuz geometrik devrimine yelken açacağız.

Kaynakça

1)Arnol'd, V. I. (1989). Mathematical Methods of Classical Mechanics (2nd ed.). Springer-Verlag.
2)Landau, L. D., & Lifshitz, E. M. (1976). Mechanics (Vol. 1, 3rd ed.). Butterworth-Heinemann

E

Eren Sivri

Yazar