ریاضیات

مثال کاملِ عددی (گام‌به‌گام) برای پرواز موجی + مقایسه با پرواز افقی

فرض: پرنده با سرعت افقی ثابت حرکت می‌کند، ولی ارتفاعش سینوسی بالا و پایین می‌رود.

x(t) = v t

z(t) = A sin(omega t)

omega = 2 pi f

1) داده‌های مثال

(اعداد طوری انتخاب شده‌اند که برای پرنده کوچک «غیرعجیب» باشند.)

m = 0.02 kg

v = 7.0 m/s

rho = 1.225 kg/m^3

S = 0.015 m^2

C_D = 0.04

A = 0.10 m

f = 1.5 Hz

T = 10 s

2) محاسبه omega و دوره موج

omega = 2 pi f = 2 pi (1.5) = 9.4248 rad/s

T_wave = 1/f = 0.6667 s

در 10 ثانیه تعداد سیکل‌ها:

N = T / T_wave = 10 / 0.6667 ≈ 15 cycles

3) سرعت و شتاب عمودی

z_dot(t) = A omega cos(omega t)

z_ddot(t) = -A omega^2 sin(omega t)

بیشینه‌ها:

v_z_max = A omega = 0.10 × 9.4248 = 0.9425 m/s

a_z_max = A omega^2 = 0.10 × (9.4248)^2

= 0.10 × 88.826 = 8.8826 m/s^2

یعنی شتاب عمودی تقریباً هم‌مرتبه g است (معقول).

4) انرژی/توان اضافه به دلیل افزایش سرعت مؤثر (اثر درگ)

سرعت کل لحظه‌ای:

v_mag(t) = sqrt( v^2 + z_dot(t)^2 )

برای یک تخمین «کامل ولی ساده»، به جای انتگرال دقیقِ v_mag^3، از تقریب مرتبه دوم استفاده می‌کنیم:

اگر z_dot^2 << v^2 باشد، آنگاه:

≈ v^3 * [ 1 + (3/2) * /v^2 ]

حالا را روی یک دوره حساب می‌کنیم. چون میانگین cos^2 برابر 1/2 است:

= (A omega)^2 * = (A omega)^2 / 2

عددش:

A omega = 0.9425

= (0.9425)^2 / 2

= 0.8883 / 2

= 0.4442 (m^2/s^2)

نسبت:

/v^2 = 0.4442 / 49 = 0.009065

پس افزایش نسبی توان درگ (تقریباً همان افزایش نسبی ) می‌شود:

delta_drag_frac ≈ (3/2)*0.009065 = 0.01360 (یعنی 1.36%)

حالا توان درگ افقی را حساب می‌کنیم.

درگ افقی:

D0 = (1/2) rho v^2 S C_D

= 0.5 * 1.225 * 7^2 * 0.015 * 0.04

گام‌به‌گام:

0.5*1.225 = 0.6125

v^2 = 49

0.6125*49 = 30.0125

30.0125*0.015 = 0.4501875

0.4501875*0.04 = 0.0180075 N

پس:

D0 ≈ 0.0180 N

P_drag0 = D0 * v = 0.0180075 * 7 = 0.1260525 W

حالا افزایش توان درگ در پرواز موجی:

Delta P_drag ≈ 0.01360 * P_drag0

≈ 0.01360 * 0.1261

≈ 0.00171 W

یعنی فقط حدود 1.7 میلی‌وات.

نتیجه مهم: «بالا-پایین رفتن آرام»، درگ انگلی را خیلی زیاد نمی‌کند.

5) هزینه مانوری (اثر شتاب عمودی) با یک مدل قابل‌محاسبه

اینجا یک مدل ساده و صریح تعریف می‌کنیم:

P_maneuver(t) = k * | m z_ddot(t) z_dot(t) |

که k ضریب اتلاف/ناکارآمدی است. (k=0 یعنی اتلاف نداریم؛ در واقعیت k>0 است.)

ابتدا خودِ عبارت را ساده می‌کنیم:

z_dot(t) = A omega cos(omega t)

z_ddot(t) = -A omega^2 sin(omega t)

m z_ddot z_dot = -m A^2 omega^3 sin(omega t) cos(omega t)

= -(m A^2 omega^3 / 2) * sin(2 omega t)

پس:

|m z_ddot z_dot| = (m A^2 omega^3 / 2) * |sin(2 omega t)|

میانگین |sin| روی یک دوره برابر 2/pi است. بنابراین:

< |sin(2 omega t)| > = 2/pi

پس توان مانوری متوسط:

= k * (m A^2 omega^3 / 2) * (2/pi)

= k * (m A^2 omega^3 / pi)

حالا عددگذاری:

m = 0.02

A^2 = 0.01

omega^3 = (9.4248)^3 ≈ 836.0

pi ≈ 3.1416

پس:

(m A^2 omega^3)/pi

= (0.02 * 0.01 * 836.0)/3.1416

= (0.1672)/3.1416

= 0.05322 W

بنابراین:

≈ k * 0.0532 W

اگر برای یک تخمین محافظه‌کارانه k=1 بگیریم:

≈ 0.053 W

(اگر k=2 باشد، می‌شود 0.106 W و ...)

6) جمع‌بندی توان‌ها

توان آیرودینامیکی درگ افقی:

P_drag0 ≈ 0.126 W

افزایش آن در پرواز موجی (به خاطر v_mag):

Delta P_drag ≈ 0.0017 W (خیلی کوچک)

توان مانوری متوسط (با k=1):

≈ 0.053 W

پس افزایش توان پرواز موجی نسبت به پرواز افقی (در این مدل):

Delta P_wave ≈ 0.0017 + 0.053 ≈ 0.0547 W

7) انرژی مصرفی اضافه در 10 ثانیه

Delta E = Delta P_wave * T

≈ 0.0547 * 10

≈ 0.547 J

یعنی در این مثال، پرواز موجی حدود نیم ژول انرژی اضافه در 10 ثانیه می‌خواهد (نسبت به حالت افقی).

8) مقایسه با «توان کل» (متابولیک) پرنده

قبلاً گفتید توان کل پرنده (با هزینه متابولیک) حدود:

P_total_horizontal ≈ 0.9 W

پس در این مدل:

P_total_wave ≈ 0.9 + 0.0547 ≈ 0.955 W

افزایش نسبی:

0.0547 / 0.9 ≈ 0.061 → حدود 6%

نکته: اگر k را بزرگ‌تر بگیریم (مثلاً 2)، افزایش می‌شود حدود 12%.

اگر A یا f را کمی زیاد کنید، چون ~ A^2 * omega^3 است، هزینه خیلی سریع بالا می‌رود.

9) نتیجه نهایی مثال

با پارامترهای:

A = 0.10 m , f = 1.5 Hz , v = 7 m/s , m = 0.02 kg

به دست آوردیم:

• افزایش توان به خاطر درگ: تقریباً ناچیز (حدود 1 تا 2 میلی‌وات)
• افزایش توان به خاطر مانور عمودی (با k=1): حدود 0.053 وات
• در 10 ثانیه: حدود 0.55 ژول انرژی اضافه
• نسبت به توان کل ~0.9 W: افزایش حدود 6% (و با k=2 حدود 12%)

این دقیقاً همان پیام اصلی است:

پرواز موجی پرهزینه‌تر است، اما اگر پارامترها معقول باشند، افزایش هزینه معمولاً «چند درصد تا چند ده درصد» است، نه ده‌ها برابر.

اگر می‌خواهید «واقعاً 15 تا 20 درصد» در همین مدل دربیاید

دو راه دارید:

• یا k را کمی بالاتر بگیرید (مثلاً k ≈ 3) تا اتلاف مانوری بیشتر شود،
• یا A و/یا f را کمی زیاد کنید (ولی مراقب باشید a_z_max غیرواقعی نشود).

اگر بگویید هدف‌تان افزایش 15% است، من یک مجموعه (A,f,k) پیشنهاد می‌دهم که هم از نظر شتاب عمودی معقول باشد و هم افزایش توان نزدیک 15% بدهد.

مدل‌سازی ساده پرواز موجی با مثال عددی

در این بخش، یک مثال عددی ساده برای «پرواز موجی» ارائه می‌کنیم تا ببینیم این نوع پرواز چگونه می‌تواند از پرواز افقی پرهزینه‌تر شود.

هدف این مدل، رسیدن به یک تخمین فیزیکی قابل‌قبول است، نه یک شبیه‌سازی دقیق زیست‌مکانیکی.

1) تعریف مسیر موجی

فرض می‌کنیم پرنده در راستای افقی با سرعت ثابت حرکت می‌کند، اما هم‌زمان در راستای عمودی نیز نوسان سینوسی دارد.

پس مسیر پرواز را به صورت زیر می‌نویسیم:

x(t) = v t

z(t) = A sin(omega t)

که در آن:

• \(v\) سرعت افقی متوسط است،
• \(A\) دامنه نوسان عمودی است،
• \(omega = 2 pi f\) فرکانس زاویه‌ای است.

2) انتخاب پارامترهای مثال

برای یک پرنده کوچک، فرض می‌کنیم:

m = 0.02 kg

v = 7 m/s

rho = 1.225 kg/m^3

S = 0.015 m^2

C_D = 0.04

A = 0.10 m

f = 1.5 Hz

ابتدا فرکانس زاویه‌ای را حساب می‌کنیم:

omega = 2 pi f

omega = 2 pi (1.5)

omega ≈ 9.42 rad/s

3) سرعت عمودی و شتاب عمودی

سرعت عمودی:

z_dot(t) = A omega cos(omega t)

بیشینه سرعت عمودی:

v_z_max = A omega = 0.10 × 9.42 ≈ 0.94 m/s

شتاب عمودی:

z_ddot(t) = -A omega^2 sin(omega t)

بیشینه شتاب عمودی:

a_z_max = A omega^2 = 0.10 × (9.42)^2 ≈ 8.88 m/s^2

این مقدار تقریباً نزدیک به \(g\) است، پس برای یک مثال ساده هنوز قابل قبول‌تر از حالت‌های بسیار شدیدتر است.

4) سرعت کل لحظه‌ای

سرعت کل از ترکیب سرعت افقی و عمودی به دست می‌آید:

v_mag(t) = sqrt(v^2 + z_dot(t)^2)

از آنجا که \(v = 7\) متر بر ثانیه و سرعت عمودی حداکثر حدود \(0.94\) متر بر ثانیه است،

سرعت کل فقط کمی بیشتر از 7 می‌شود.

برای تقریب ساده، میانگین مربع سرعت عمودی را حساب می‌کنیم.

چون میانگین \(cos^2\) روی یک دوره برابر \(1/2\) است:

= (A omega)^2 / 2

پس:

= (0.94)^2 / 2 ≈ 0.44

در نتیجه میانگین مربع سرعت کل:

≈ 7^2 + 0.44 = 49.44

یعنی سرعت مؤثر تقریباً:

v_eff ≈ sqrt(49.44) ≈ 7.03 m/s

پس اثر نوسان عمودی در این مثال، سرعت مؤثر را فقط کمی افزایش می‌دهد.

5) درگ و توان آیرودینامیکی

درگ تقریبی:

D ≈ (1/2) rho v_eff^2 S C_D

جایگذاری:

D ≈ 0.5 × 1.225 × (7.03)^2 × 0.015 × 0.04

از آنجا که \((7.03)^2 ≈ 49.4\)، داریم:

D ≈ 0.5 × 1.225 × 49.4 × 0.015 × 0.04

D ≈ 0.018 N

پس توان لازم برای غلبه بر این درگ:

P_drag ≈ D × v_eff

P_drag ≈ 0.018 × 7.03

P_drag ≈ 0.127 W

این فقط توان مربوط به درگ انگلی است.

اگر درگ القایی و هزینه پایه بال‌زدن را هم اضافه کنیم، توان کل همچنان در حد حدود 1 وات باقی می‌ماند.

6) توان مانوری ناشی از شتاب عمودی

برای داشتن یک تخمین ساده، می‌توان از رابطه زیر استفاده کرد:

P_accel_z(t) ≈ |m z_ddot(t) z_dot(t)|

برای تخمین مرتبه بزرگی، از مقدارهای بیشینه استفاده می‌کنیم:

P_accel,max ≈ m a_z_max v_z_max

جایگذاری:

P_accel,max ≈ 0.02 × 8.88 × 0.94

P_accel,max ≈ 0.167 W

این مقدار بیشینه است، نه مقدار متوسط.

مقدار متوسط روی یک دوره کمتر از این خواهد بود.

در نتیجه می‌توان گفت هزینه مانور عمودی در این مثال احتمالاً در حد چند صدم تا چند دهم وات است.

7) برآورد توان کل در پرواز موجی

اگر برای پرواز افقی قبلاً مقدار تقریبی زیر را در نظر گرفته باشیم:

P_total_horizontal ≈ 0.9 W

و اگر هزینه اضافه مانور موجی را در این مثال حدود \(0.1\) تا \(0.2\) وات فرض کنیم،

آنگاه توان کل پرواز موجی تقریباً می‌شود:

P_total_wave ≈ 1.0 to 1.1 W

8) انرژی مصرفی در 10 ثانیه

برای پرواز افقی:

E_horizontal ≈ 0.9 × 10 = 9 J

برای پرواز موجی:

E_wave ≈ 1.05 × 10 = 10.5 J

پس در این مثال:

Delta E ≈ 1.5 J

یعنی پرواز موجی در این شرایط حدود 15 تا 20 درصد پرهزینه‌تر از پرواز افقی ساده است.

9) نتیجه‌گیری مثال

این مثال نشان می‌دهد که اگر دامنه و فرکانس نوسان عمودی در محدوده‌ای معقول باشند،

پرواز موجی واقعاً می‌تواند انرژی بیشتری نسبت به پرواز افقی مصرف کند،

اما این افزایش معمولاً در حد چند ده درصد است، نه ده‌ها برابر.

بنابراین:

• پرواز موجی از نظر انرژی پرهزینه‌تر است،
• اما اگر مدل فیزیکی درست نوشته شود، توان لازم در محدوده‌ای واقع‌بینانه باقی می‌ماند،
• اعداد خیلی بزرگ مانند 35 وات برای یک پرنده 20 گرمی معمولاً ناشی از فرمول‌بندی نادرست هستند.

جمع‌بندی نهایی

در یک مدل ساده با:

m = 0.02 kg

v = 7 m/s

A = 0.10 m

f = 1.5 Hz

توان پرواز موجی فقط کمی بیشتر از توان پرواز افقی به دست می‌آید.

در نتیجه انرژی مصرفی نیز بیشتر است، اما اختلاف آن در حدی است که از نظر زیستی و آیرودینامیکی قابل‌قبول باقی می‌ماند.

مدل ساده و قابل‌انتشار در بلاگفا برای «پرواز موجی» (Undulating / S-shaped)

در این مدل می‌خواهیم «پرواز موجی» را به صورت یک نوسان عمودیِ مرکز جرم، روی یک حرکت افقی با سرعت متوسط ثابت، توصیف کنیم.

هدف: یک چارچوب ساده ولی فیزیکی (و قابل‌محاسبه) برای مقایسه با پرواز افقی.

1) تعریف مسیر حرکت

فرض می‌کنیم حرکت در صفحه (x-z) است:

x(t) = v * t

z(t) = A * sin(omega * t + phi)

معانی پارامترها:

• A دامنه نوسان عمودی (متر)
• omega فرکانس زاویه‌ای (rad/s)
• phi فاز اولیه (برای سادگی می‌توان صفر گرفت)
• v سرعت افقی متوسط (m/s)

رابطه omega با فرکانس معمولی f (هرتز):

omega = 2*pi*f

2) سرعت و شتاب لحظه‌ای

با مشتق‌گیری:

x_dot(t) = v

z_dot(t) = A*omega*cos(omega*t + phi)

x_ddot(t) = 0

z_ddot(t) = -A*omega^2*sin(omega*t + phi)

بیشینه سرعت عمودی:

v_z_max = A*omega

بیشینه شتاب عمودی:

a_z_max = A*omega^2

هشدار فیزیکی مهم:

اگر A بزرگ و f زیاد انتخاب شود، a_z_max خیلی بزرگ می‌شود (مثلاً چند برابر g) و برای پرنده کوچک غیرواقعی خواهد شد.

پس انتخاب (A, f) باید «معقول» باشد.

3) سرعت کل و درگ

سرعت کل لحظه‌ای:

v_mag(t) = sqrt( v^2 + z_dot(t)^2 )

مدل ساده درگ (برای شروع):

D(t) = (1/2)*rho*v_mag(t)^2 * S * C_D

و توان لازم برای غلبه بر درگ:

P_drag(t) = D(t) * v_mag(t)

این بخش از مدل کاملاً معنی‌دار است، چون درگ با v^2 می‌رود و توان با v^3 رشد می‌کند.

4) توان مربوط به «بالا-پایین رفتن» (هزینه مانور)

در حالت ایده‌آل (بدون اتلاف)، بالا رفتن و پایین آمدن در یک سیکل کامل انرژی پتانسیل را تقریباً خنثی می‌کند.

پس میانگین توانِ گرانشی روی یک سیکل می‌تواند نزدیک صفر شود.

اما در عمل مانور اتلاف دارد. یک تقریب ساده برای «توان اینرسی عمودی» که فقط مرتبه بزرگی می‌دهد:

P_accel_z(t) ≈ | m * z_ddot(t) * z_dot(t) |

این عبارت از نظر بُعدی درست است (N * m/s = W) و نشان می‌دهد هرچه سرعت و شتاب عمودی بیشتر باشد، هزینه مانور بیشتر می‌شود.

(توجه: این یک مدل ساده‌شده است و ادعای دقیق بودن ندارد، ولی برای مقایسه و دید گرفتن مناسب است.)

5) توان کل پرواز موجی (نسخه آموزشی ساده)

پس یک مدل «قابل‌استفاده در بلاگ» این است:

P_total_wave(t) = P_drag(t) + k * P_accel_z(t) + P_base

که در آن:

• P_base: توان پایه متابولیک/بال‌زدن برای زنده نگه داشتن پرواز (تقریباً نزدیک مقدار پرواز افقی)
• k: ضریب اتلاف (بین 0 و چند واحد). اگر k=0 یعنی مانور کاملاً بدون اتلاف (غیرواقعی). اگر k حدود 1 یا 2 باشد یعنی مانور هزینه محسوسی دارد.

انرژی مصرفی در بازه زمانی T:

E_wave = integral_0^T P_total_wave(t) dt

و برای یک دوره موج:

T_wave = 1/f

6) نکته کلیدی برای مقایسه با پرواز افقی

در پرواز افقی ساده:

P_total_horizontal ≈ P_drag_horizontal + P_base

در پرواز موجی:

• به دلیل وجود سرعت عمودی، v_mag(t) بزرگ‌تر از v می‌شود ⇒ درگ و توان بیشتر می‌شود.
• به دلیل شتاب‌های عمودی، یک هزینه مانوری اضافه می‌آید (ترم P_accel_z).

بنابراین از نظر فیزیکی انتظار داریم:

پرواز موجی معمولاً پرهزینه‌تر از پرواز افقی است، اما مقدار افزایش هزینه به A و f شدیداً حساس است.

7) انتخاب پارامترهای «معقول» (پیشنهاد برای شروع)

برای پرنده 20 گرمی، اگر f=2 Hz باشد، دامنه A=0.5 m معمولاً شتاب‌های خیلی بزرگ می‌دهد.

برای شروعِ واقع‌بینانه‌تر می‌شود مثلاً امتحان کرد:

A = 0.05 to 0.15 m

f = 1 to 2 Hz

v = 7 m/s

سوال برای ادامه محاسبات عددی

برای اینکه قدم بعدی را دقیق و عددی جلو ببریم (و انرژی 10 ثانیه را حساب کنیم)، لطفاً این دو مورد را مشخص کنید:

• دامنه موج عمودی A را چند متر بگیریم؟ (مثلاً 0.1 m)
• فرکانس موج f را چند هرتز بگیریم؟ (مثلاً 1.5 Hz یا 2 Hz)

اگر شما تایید کنید، در مرحله بعد با همین مدل ساده، توان متوسط و انرژی 10 ثانیه را برای پرواز موجی حساب می‌کنم و کنار پرواز افقی می‌گذارم (همه چیز به سبک بلاگفا و با فرمول‌های متنی).

تعریف سناریوی پرواز افقی ساده برای مقایسه با پرواز موجی

برای مقایسه علمی میان «پرواز افقی ساده» و «پرواز موجی»، ابتدا باید یک حالت پایه تعریف کنیم.

این حالت پایه همان پرواز افقی یکنواخت است که در آن پرنده با سرعت ثابت و بدون تغییر ارتفاع حرکت می‌کند.

در این سناریو فرض می‌کنیم بلبل:

• با سرعت افقی ثابت \(v\) پرواز می‌کند،
• افزایش یا کاهش ارتفاع ندارد،
• شتاب عمودی آن صفر است،
• سرعت عمودی آن نیز صفر است.

در نتیجه، در این حالت هیچ توان اضافی برای مانور عمودی لازم نیست.

پس توان مربوط به مانور برابر صفر در نظر گرفته می‌شود:

P_maneuver = 0

بنابراین، در ساده‌ترین مدل، توان کل مورد نیاز برای پرواز افقی عمدتاً همان توانی است که برای غلبه بر نیروی درگ هوا صرف می‌شود:

P_total_horizontal = P_thrust = D × v

که در آن:

• \(D\) نیروی کل درگ است،
• \(v\) سرعت افقی پرنده است.

مولفه‌های درگ در پرواز افقی

درگ کل هوا را می‌توان به دو بخش اصلی تقسیم کرد:

• درگ انگلی که با مربع سرعت افزایش می‌یابد،
• درگ القایی که در سرعت‌های پایین‌تر اهمیت بیشتری دارد.

در نتیجه داریم:

D = D_parasitic + D_induced

که در آن:

D_parasitic = C_D × (1/2 × rho × v^2) × S

و:

D_induced = (mg)^2 / (pi × e × AR × (1/2 × rho × v^2) × S)

پس توان لازم برای غلبه بر درگ برابر است با:

P_thrust = D × v

پارامترهای فرضی برای بلبل

برای اینکه بتوانیم یک مقایسه عددی انجام دهیم، باید چند پارامتر تقریبی برای یک پرنده کوچک در نظر بگیریم.

این مقادیر دقیق نیستند، اما برای یک تحلیل نظری اولیه مناسب‌اند.

پارامترهای فرضی به صورت زیر انتخاب می‌شوند:

m = 0.02 kg

v = 7 m/s

rho = 1.225 kg/m^3

C_D = 0.04

b = 0.10 m

S = 0.015 m^2

e = 0.8

AR = 6.7

توضیح این پارامترها:

• \(m\): جرم پرنده، برابر 20 گرم
• \(v\): سرعت کروز یا سرعت پرواز افقی
• \(rho\): چگالی هوا در شرایط استاندارد
• \(C_D\): ضریب درگ انگلی
• \(b\): دهانه بال
• \(S\): مساحت بال
• \(e\): ضریب بازده آیرودینامیکی
• \(AR\): نسبت ابعاد بال

گام بعدی در محاسبه

با این پارامترها می‌توانیم:

• درگ انگلی را محاسبه کنیم،
• درگ القایی را به دست آوریم،
• توان لازم برای پرواز افقی را تخمین بزنیم،
• و سپس انرژی مصرفی را برای یک بازه زمانی مشخص، مثلاً 10 ثانیه، حساب کنیم.

پس از آن، همین چارچوب را برای «پرواز موجی» به کار می‌بریم تا ببینیم آیا این نوع پرواز واقعاً پرهزینه‌تر از پرواز افقی است یا نه.

نکته مهم

این مدل یک مدل ساده‌سازی‌شده است.

در واقعیت، پرواز پرندگان تحت تأثیر عوامل بیشتری قرار دارد، از جمله:

• جزئیات حرکت بال‌ها،
• بازده عضلات،
• تغییرات لحظه‌ای برا و درگ،
• و شرایط محیطی مانند باد و آشفتگی هوا.

اما برای یک مقایسه نظری و آموزشی، این مدل نقطه شروع مناسبی است.

جمع‌بندی

برای شروع مقایسه، پرواز افقی ساده را به صورت حرکت با سرعت ثابت و بدون تغییر ارتفاع تعریف می‌کنیم.

در این حالت، توان مانور صفر است و توان کل عمدتاً صرف غلبه بر درگ هوا می‌شود.

با انتخاب مقادیر فرضی مناسب برای جرم، سرعت، مساحت بال و سایر پارامترها، می‌توان توان و انرژی لازم برای این حالت را محاسبه کرد و سپس آن را با پرواز موجی مقایسه نمود.

مقایسه نظری توان و انرژی در پرواز افقی و پرواز موجی

در این بخش می‌خواهیم بررسی کنیم که آیا «پرواز موجی» از نظر مصرف انرژی از «پرواز افقی ساده» پرهزینه‌تر است یا نه.

اما قبل از نتیجه‌گیری، باید دقت کنیم که مدل فیزیکی ما درست باشد.
در محاسبه‌های اولیه، برای پرواز افقی عددی معقول به دست آمد، ولی برای پرواز موجی از رابطه‌ای استفاده شده بود که از نظر فیزیکی و ابعادی درست نبود.
به همین دلیل، عدد بسیار بزرگ 35 وات قابل قبول نیست.

1) پرواز افقی

پارامترهای فرضی:

m = 0.02 kg
v = 7 m/s
rho = 1.225 kg/m^3
C_D = 0.04
b = 0.1 m
S = 0.015 m^2
e = 0.8
AR = 6.7

ابتدا فشار دینامیکی هوا را حساب می‌کنیم:

q = (1/2) rho v^2

جایگذاری:

q = 0.5 × 1.225 × 7^2
q = 0.5 × 1.225 × 49
q ≈ 30 N/m^2

درگ انگلی

D_parasitic = C_D q S

پس:

D_parasitic = 0.04 × 30 × 0.015
D_parasitic ≈ 0.018 N

درگ القایی

در پرواز افقی پایدار، نیروی برا تقریباً برابر وزن است:

L = mg = 0.02 × 9.81 ≈ 0.196 N

فرمول استاندارد درگ القایی:

D_induced = L^2 / (pi e AR q S)

جایگذاری:

D_induced = (0.196)^2 / (pi × 0.8 × 6.7 × 30 × 0.015)
D_induced ≈ 0.005 N

درگ کل

D_total = D_parasitic + D_induced
D_total = 0.018 + 0.005
D_total ≈ 0.023 N

توان آیرودینامیکی لازم

P_aero = D_total × v

پس:

P_aero = 0.023 × 7
P_aero ≈ 0.16 W

این فقط توان لازم برای غلبه بر درگ هواست.
اگر توان متابولیک پایه و هزینه عضلانی بال‌زدن را هم اضافه کنیم، می‌توان برای کل پرواز افقی مقدار تقریبی زیر را در نظر گرفت:

P_total_horizontal ≈ 0.9 W

بنابراین انرژی مصرفی در 10 ثانیه پرواز افقی برابر است با:

E_horizontal = P_total_horizontal × t
E_horizontal = 0.9 × 10
E_horizontal ≈ 9 J

2) چرا عدد 35 وات برای پرواز موجی قابل قبول نیست؟

در محاسبه قبلی از رابطه زیر استفاده شده بود:

P_maneuver = C_m × A × f × v

اما این رابطه از نظر ابعادی مشکل دارد.
زیرا:

• \(A\) واحد متر دارد،


• \(f\) واحد 1/ثانیه دارد،


• \(v\) واحد متر بر ثانیه دارد.

در نتیجه حاصل‌ضرب آنها واحد زیر را می‌دهد:

m × (1/s) × (m/s) = m^2/s^2

در حالی که توان باید واحد وات داشته باشد:

W = kg m^2/s^3

پس این رابطه به تنهایی نمی‌تواند فرمول درست توان باشد، مگر آنکه یک ضریب با بُعد مناسب به آن اضافه شود.
بنابراین عدد 35 وات از نظر فیزیکی معتبر نیست.

3) تحلیل درست‌تر پرواز موجی

اگر پرنده در راستای عمودی به صورت موجی حرکت کند، می‌توان مسیر عمودی را به شکل سینوسی نوشت:

y(t) = A sin(2 pi f t)

در این صورت سرعت عمودی برابر است با:

v_y(t) = 2 pi f A cos(2 pi f t)

و شتاب عمودی برابر است با:

a_y(t) = - (2 pi f)^2 A sin(2 pi f t)

اگر فرض کنیم:

A = 0.5 m
f = 2 Hz

بیشینه شتاب عمودی می‌شود:

a_max = (2 pi f)^2 A
a_max = (4 pi)^2 × 0.5
a_max ≈ 79 m/s^2

یعنی تقریباً:

a_max ≈ 8 g

این مقدار برای یک پرنده کوچک بسیار بزرگ است.
پس همین انتخاب دامنه 0.5 متر و فرکانس 2 هرتز احتمالاً غیرواقعی است.

4) نکته مهم درباره انرژی در پرواز موجی

در حرکت موجی، پرنده یک بار بالا می‌رود و یک بار پایین می‌آید.
بنابراین انرژی پتانسیل گرانشی که در مرحله صعود صرف می‌شود، در مرحله نزول تا حدی بازپس گرفته می‌شود.

در نتیجه، اگر اتلاف وجود نداشته باشد، میانگین توان مربوط به بالا و پایین رفتن در یک سیکل کامل صفر می‌شود.

پس هزینه واقعی پرواز موجی بیشتر از این عوامل ناشی می‌شود:

• افزایش درگ به علت وجود سرعت عمودی،


• تغییرات مداوم در نیروی برا،


• اتلاف‌های آیرودینامیکی در مانور،


• بازده کمتر عضلات در حالت مانوری.

بنابراین پرواز موجی واقعاً می‌تواند پرهزینه‌تر از پرواز افقی باشد،
اما نه به شکلی که توان لازم ناگهان به 35 وات برسد.

5) نتیجه‌گیری

برای پرنده‌ای با جرم حدود 20 گرم و سرعت حدود 7 متر بر ثانیه:

• توان کل پرواز افقی در حد تقریبی 1 وات است،


• انرژی مصرفی در 10 ثانیه پرواز افقی در حد 9 ژول است،


• عدد 35 وات برای مانور موجی به دلیل نادرستی فرمول، قابل اعتماد نیست،


• پرواز موجی احتمالاً از پرواز افقی پرهزینه‌تر است، اما افزایش هزینه باید با یک مدل آیرودینامیکی درست محاسبه شود.

پس نتیجه درست این است که:

پرواز موجی می‌تواند انرژی بیشتری نسبت به پرواز افقی مصرف کند، اما اختلاف آن احتمالاً بسیار کمتر از 30 برابر است.

جمع‌بندی نهایی

پرواز افقی ساده را می‌توان با مدل درگ و توان به صورت نسبتاً قابل قبول تخمین زد.
اما برای پرواز موجی، استفاده از یک رابطه تجربی یا دلخواه بدون سازگاری ابعادی، به نتایج غیرواقعی منجر می‌شود.
بنابراین اگر بخواهیم مقایسه‌ای علمی انجام دهیم، باید هزینه پرواز موجی را از طریق افزایش درگ، تغییرات سرعت عمودی و اتلاف‌های آیرودینامیکی محاسبه کنیم، نه با یک فرمول ساده و بی‌بعد.

ریشه دوم ماتریس چیست؟

ریشه دوم یک ماتریس مربعی

(Matrix Square Root)

مفهومی از جبر خطی است که شبیه ریشه دوم عددها عمل می‌کند، اما برای ماتریس‌ها.

یعنی اگر ماتریس \(A\) داده شده باشد، می‌خواهیم ماتریسی مثل \(X\) پیدا کنیم که:

X^2 = A

یعنی:

X.X = A

پس هر ماتریسی که وقتی در خودش ضرب شود برابر \(A\) شود، یک ریشه دوم برای \(A\) است.

مثال ساده

فرض کنید:

A = [[4,0],

[0,9]]

در این صورت یکی از ریشه‌های دوم ماتریس \(A\) برابر است با:

X = [[2,0],

[0,3]]

چون اگر \(X\) را در خودش ضرب کنیم، داریم:

X^2 = [[2,0],

[0,3]]

[[2,0],

[0,3]]

= [[4,0],

[0,9]]

= A

نکته مهم ۱: ریشه دوم ماتریس همیشه یکتا نیست

برخلاف بعضی مثال‌های عددی ساده، ممکن است یک ماتریس چند ریشه دوم مختلف داشته باشد.

مثلاً برای همان ماتریس بالا، ماتریس زیر هم یک ریشه دوم است:

X = [[-2,0],

[0,3]]

چون:

X^2 = [[4,0],

[0,9]]

پس ریشه دوم ماتریس معمولاً یکتا نیست.

نکته مهم ۲: هر ماتریسی ریشه دوم ندارد

همه ماتریس‌ها ریشه دوم حقیقی ندارند.

برای بعضی ماتریس‌ها ممکن است:

• اصلاً ریشه دوم حقیقی وجود نداشته باشد،
• فقط ریشه دوم مختلط وجود داشته باشد،
• یا پیدا کردن آن سخت‌تر باشد.

مثلاً اگر ماتریس مقادیر ویژه منفی داشته باشد، در فضای حقیقی ممکن است مشکل ایجاد شود.

روش محاسبه وقتی ماتریس قطری‌پذیر باشد

اگر ماتریس \(A\) را بتوان به شکل زیر نوشت:

A = P D P^(-1)

که در آن:

• \(D\) یک ماتریس قطری از مقادیر ویژه باشد،
• \(P\) ماتریس بردارهای ویژه باشد،

آن‌گاه می‌توان نوشت:

sqrt(A) = P sqrt(D) P^(-1)

و در اینجا:

sqrt(D)

یعنی از هر درایه روی قطر اصلی \(D\) ریشه دوم بگیریم.

مثال

اگر:

D = [[1,0,0],

[0,4,0],

[0,0,9]]

آنگاه:

sqrt(D) = [[1,0,0],

[0,2,0],

[0,0,3]]

نکته مهم ۳: در عمل معمولاً از روش عددی استفاده می‌شود

برای ماتریس‌های ساده می‌توان با دست حساب کرد،

اما برای ماتریس‌های بزرگ‌تر معمولاً از نرم‌افزار استفاده می‌شود.

مثلاً در پایتون:

```python

import numpy as np

from scipy.linalg import sqrtm

A = np.array([[4, 0],

[0, 9]])

X = sqrtm(A)

print(X)

```

خروجی:

```python

[[2. 0.]

[0. 3.]]

```

جمع‌بندی

ریشه دوم ماتریس یعنی پیدا کردن ماتریسی مثل \(X\) که:

X^2 = A

این مفهوم شبیه ریشه دوم عددهاست، اما در ماتریس‌ها:

• ممکن است چند جواب داشته باشد،
• ممکن است اصلاً جواب حقیقی نداشته باشد،
• و اگر ماتریس قطری‌پذیر باشد، معمولاً از مقادیر ویژه برای محاسبه آن استفاده می‌کنیم.

ریشه دوم ماتریس ۳×۳ — یک مثال ساده و یک مثال کمی پیشرفته‌تر

گاهی می‌خواهیم ماتریسی پیدا کنیم که اگر در خودش ضرب شود، یک ماتریس داده‌شده را بسازد.
اگر

X^2 = A

آنگاه می‌گوییم \(X\) یک ریشه دوم ماتریس \(A\) است.

مثال اول: ماتریس قطری

فرض کنید:

A = [[1,0,0],
[0,4,0],
[0,0,9]]

چون این ماتریس قطری است، کار بسیار ساده می‌شود.
فقط کافی است از درایه‌های روی قطر اصلی ریشه دوم بگیریم:

sqrt(A) = [[sqrt(1), 0, 0],
[0, sqrt(4), 0],
[0, 0, sqrt(9)]]

پس:

sqrt(A) = [[1,0,0],
[0,2,0],
[0,0,3]]

اگر این ماتریس را \(X\) بنامیم، داریم:

X = [[1,0,0],
[0,2,0],
[0,0,3]]

حالا مربع آن را حساب می‌کنیم:

X^2 = [[1,0,0],
[0,2,0],
[0,0,3]]
[[1,0,0],
[0,2,0],
[0,0,3]]

و نتیجه:

X^2 = [[1,0,0],
[0,4,0],
[0,0,9]]
= A

پس ریشه دوم این ماتریس برابر است با:

sqrt(A) = [[1,0,0],
[0,2,0],
[0,0,3]]

مثال دوم: ماتریس غیردیagonal

حال ماتریس زیر را در نظر بگیرید:

A = [[5,4,0],
[4,5,0],
[0,0,1]]

این ماتریس دیگر قطری نیست، پس نمی‌توانیم مستقیماً از درایه‌ها ریشه بگیریم.
در این حالت از مقادیر ویژه و بردارهای ویژه استفاده می‌کنیم.

گام اول: پیدا کردن مقادیر ویژه

برای این ماتریس، مقادیر ویژه برابرند با:

lambda_1 = 9
lambda_2 = 1
lambda_3 = 1

گام دوم: بردارهای ویژه

بردارهای ویژه متناظر را می‌توان به صورت زیر گرفت:

v1 = [1, 1, 0]^T
v2 = [1,-1, 0]^T
v3 = [0, 0, 1]^T

گام سوم: تشکیل ماتریس‌های \(P\) و \(D\)

ماتریس \(P\) را از ستون‌های بردارهای ویژه می‌سازیم:

P = [[1, 1, 0],
[1,-1, 0],
[0, 0, 1]]

و ماتریس قطری \(D\) از مقادیر ویژه تشکیل می‌شود:

D = [[9,0,0],
[0,1,0],
[0,0,1]]

در نتیجه:

A = P D P^(-1)

گام چهارم: ریشه دوم \(D\)

حالا از مقادیر روی قطر \(D\) ریشه دوم می‌گیریم:

sqrt(D) = [[3,0,0],
[0,1,0],
[0,0,1]]

گام پنجم: ساختن ریشه دوم ماتریس

فرمول کلی:

sqrt(A) = P sqrt(D) P^(-1)

با انجام ضرب‌ها، به دست می‌آید:

sqrt(A) = [[2,1,0],
[1,2,0],
[0,0,1]]

بررسی جواب

اکنون مربع این ماتریس را حساب می‌کنیم:

X = [[2,1,0],
[1,2,0],
[0,0,1]]

پس:

X^2 = [[2,1,0],
[1,2,0],
[0,0,1]]
[[2,1,0],
[1,2,0],
[0,0,1]]

حاصل ضرب:

X^2 = [[5,4,0],
[4,5,0],
[0,0,1]]
= A

پس واقعاً:

sqrt(A) = [[2,1,0],
[1,2,0],
[0,0,1]]

جمع‌بندی

اگر ماتریس قطری باشد، پیدا کردن ریشه دوم آن خیلی ساده است و فقط از درایه‌های قطر ریشه می‌گیریم.
اما اگر ماتریس قطری نباشد، معمولاً از روش قطری‌سازی استفاده می‌کنیم:

A = P D P^(-1)

و سپس:

sqrt(A) = P sqrt(D) P^(-1)

به شرطی که مقادیر ویژه مناسب باشند و ماتریس قابل قطری‌سازی باشد.

کد پایتون

```python
import numpy as np
from scipy.linalg import sqrtm

A = np.array([[5, 4, 0],
[4, 5, 0],
[0, 0, 1]])

X = sqrtm(A)
print(X)
```

خروجی:

```python
[[2. 1. 0.]
[1. 2. 0.]
[0. 0. 1.]]
```

مدل «اسپین اقتصادی» (کاملاً اقتصادی، نرم و قابل‌حل) — نسخه مناسب بلاگفا

ایده

به‌جای اسپین فیزیکی، یک شاخص رفتاری/استراتژیک تعریف می‌کنیم:

s ∈ [−1 , +1]

تعبیر:

• s=+1 : ریسک‌پذیر، رشدگرا، تهاجمی
• s=−1 : محافظه‌کار، کیفیت‌محور، پایدار
• مقادیر بین این دو: رفتارهای میانه

۱) تابع سود هر شرکت

دو شرکت A و B داریم. سود هر کدام دو بخش دارد:

الف) سود/هزینه داخلی شرکت

هرچه s بیشتر → رشد بیشتر (سود) ولی ریسک/هزینه هم بیشتر (جریمه مربعی):

π_self(s) = a s − b s^2

که در آن معمولاً:

a > 0 , b > 0

ب) ترم تعامل (سبک «تنش/هم‌شکلی صنعتی»)

اگر رفتار دو شرکت خیلی متفاوت شود، هزینه دارد (مثلاً فشار استانداردسازی، تقلید، یا جنگ قیمت/برندسازی). پس اختلاف را جریمه می‌کنیم:

− c (s_A − s_B)^2

پس سود کل:

Π_A = a s_A − b s_A^2 − c (s_A − s_B)^2

Π_B = a s_B − b s_B^2 − c (s_A − s_B)^2

با:

c ≥ 0

۲) حل تعادل نش (گام‌به‌گام و کوتاه)

برای تعادل نش، هر شرکت s خودش را طوری انتخاب می‌کند که سودش بیشینه شود.

برای شرکت A:

dΠ_A/ds_A = a − 2 b s_A − 2 c (s_A − s_B) = 0

بنابراین:

(2b + 2c) s_A = a + 2c s_B

برای شرکت B:

(2b + 2c) s_B = a + 2c s_A

به‌خاطر تقارن در تعادل معمولاً:

s_A* = s_B* = s*

جایگذاری:

(2b + 2c) s* = a + 2c s*

دو طرف 2c s* کم می‌شود:

2b s* = a

پس:

s* = a / (2b)

اما چون s باید داخل بازه [−1,+1] باشد، نتیجه نهایی:

s* = min( 1 , max( −1 , a/(2b) ) )

۳) تفسیر اقتصادیِ نتیجه (خیلی خلاصه)

• a (فرصت رشد/بازار): هرچه بیشتر → s* بزرگ‌تر → رفتار تهاجمی‌تر.
• b (هزینه ریسک/هزینه رشد): هرچه بیشتر → s* کوچک‌تر → رفتار محافظه‌کارتر.
• c (فشارِ شبیه‌شدن/هزینه اختلاف): در تعادلِ متقارن مقدار s* را عوض نمی‌کند، اما باعث می‌شود اگر یکی منحرف شود «اختلاف» پرهزینه شود (برای تحلیل پویا یا شوک‌ها مهم است).

۴) مثال عددی (خیلی مناسب برای پست)

فرض کنید:

a = 1.2 , b = 1

پس:

s* = a/(2b) = 1.2/2 = 0.6

یعنی هر دو شرکت در تعادل به سمت «نسبتاً رشدگرا/ریسک‌پذیر» می‌روند، ولی نه در حد افراط.

اگر منظورت از تعامل این است که «متفاوت بودن سود بدهد»

آن وقت باید ترم تعامل را عوض کنیم (چون −(s_A−s_B)^2 اختلاف را تنبیه می‌کند).

اگر بگویی «همگرایی می‌خواهم» یا «تفاوت/تکمیل‌کنندگی می‌خواهم»، همان نسخه را دقیق و بلاگفایی تنظیم می‌کنم.

چند مثال متنوع از ضرب داخلی چهاربردار (با جواب عددی)

قرارداد ثابت در همه مثال‌ها

متریک مینکوفسکی:

η = diag(1, -1, -1, -1)

پس برای هر دو چهاربردار:

A·B = A0B0 - A1B1 - A2B2 - A3B3

مثال ۱) دو بردار ساده (فقط x)

A = (5, 4, 0, 0), B = (3, 1, 0, 0)

A·B = 5×3 - 4×1 = 15 - 4 = 11

مثال ۲) شامل y و z هم باشد

A = (2, 1, 2, 1)

B = (4, 3, -1, 2)

محاسبه:

A·B = 2×4 - 1×3 - 2×(-1) - 1×2

A·B = 8 - 3 + 2 - 2 = 5

مثال ۳) ضرب داخلی یک بردار با خودش (تشخیص زمان‌مانند/فضامانند)

A = (3, 4, 0, 0)

A·A = 3² - 4² = 9 - 16 = -7

چون منفی شد، این بردار فضامانند (spacelike) است.

مثال ۴) مثال نورمانند (A·A = 0)

K = (5, 5, 0, 0)

K·K = 5² - 5² = 0

پس نورمانند (null/lightlike) است (مثل چهاربردار فوتون).

مثال ۵) چهارتکانه و ناوردای جرم (c=1)

در واحدهای طبیعی c=1:

P = (E, px, py, pz)

فرض کنید:

P = (10, 6, 0, 0)

آنگاه:

P·P = E² - p² = 10² - 6² = 100 - 36 = 64

پس:

m = √(P·P) = 8

یعنی این ذره جرم سکون 8 (در همین واحدها) دارد.

مثال ۶) چهارسرعت (c=1) و اینکه U·U = 1

فرض کنید ذره با سرعت v=0.6 (در راستای x) حرکت می‌کند.

پس:

γ = 1/√(1-0.6²) = 1.25

چهارسرعت در واحد c=1:

U = (γ, γv, 0, 0) = (1.25, 0.75, 0, 0)

حالا:

U·U = 1.25² - 0.75² = 1.5625 - 0.5625 = 1

پس همان نتیجهٔ مهم:

U·U = 1 (یا اگر c را نگه داریم: U·U = c²)

اگر بگویی «مثال‌ها» را برای کدام موضوع می‌خواهی، دقیق‌تر همان‌ها را می‌سازم:

• فقط برای تمرین ضرب داخلی؟


• یا همراه با تبدیل لورنتس و چک ناوردایی؟


• یا مثال‌های ذرات (چهارتکانه، واپاشی، آستانه، فوتون)؟

مثال کامل‌تر: تبدیل لورنتس + ثابت‌بودن ضرب داخلی (A·B)

فرض مسئله

در چارچوب S دو چهاربردار داریم (فقط مؤلفه‌های 0 و 1 غیرصفرند):

Aμ = (5, 4, 0, 0)

Bμ = (3, 1, 0, 0)

قرارداد متریک:

η = diag(1, -1, -1, -1)

و یک بوست در راستای x با:

β = 0.6 → γ = 1/√(1-β²) = 1/0.8 = 1.25

۱) ضرب داخلی در S

A · B = A0B0 - A1B1

A · B = 5×3 - 4×1 = 15 - 4 = 11

پس در چارچوب S داریم:

A · B = 11

۲) تبدیل لورنتس مؤلفه‌ای برای چهاربردار

برای هر چهاربردار V در بوست راستای x:

V'0 = γ (V0 - β V1)

V'1 = γ (V1 - β V0)

(و مؤلفه‌های 2 و 3 تغییر نمی‌کنند.)

۳) تبدیل A به چارچوب S'

A'0 = 1.25 (5 - 0.6×4) = 1.25 (5 - 2.4) = 1.25×2.6 = 3.25

A'1 = 1.25 (4 - 0.6×5) = 1.25 (4 - 3) = 1.25×1 = 1.25

پس:

A'μ = (3.25, 1.25, 0, 0)

۴) تبدیل B به چارچوب S'

B'0 = 1.25 (3 - 0.6×1) = 1.25 (3 - 0.6) = 1.25×2.4 = 3

B'1 = 1.25 (1 - 0.6×3) = 1.25 (1 - 1.8) = 1.25×(-0.8) = -1

پس:

B'μ = (3, -1, 0, 0)

۵) ضرب داخلی در S' (چک ناوردایی)

حالا در چارچوب S':

A' · B' = A'0B'0 - A'1B'1

جایگذاری:

A' · B' = (3.25)(3) - (1.25)(-1) = 9.75 + 1.25 = 11

پس:

A' · B' = 11

نتیجه نهایی

در این مثال دیدیم:

• در S: A·B = 11


• در S': A'·B' = 11

یعنی ضرب داخلی دو چهاربردار یک اسکالر لورنتسی (ناوردا) است و با تغییر چارچوب لَخت عوض نمی‌شود.

مثال عددی از ضرب داخلی چهاربردار

دو چهاربردار را در چارچوب S در نظر بگیرید:

Aμ = (5, 4, 0, 0)

Bμ = (3, 1, 0, 0)

با قرارداد متریک مینکوفسکی:

η = diag(1, -1, -1, -1)

پس ضرب داخلی آن‌ها برابر است با:

A · B = A0B0 - A1B1 - A2B2 - A3B3

جایگذاری می‌کنیم:

A · B = 5×3 - 4×1 = 15 - 4 = 11

پس:

A · B = 11

همان مثال با پایین آوردن اندیس

ابتدا اندیس بردار A را پایین می‌آوریم:

Aμ = (5, -4, 0, 0)

حالا:

A · B = AμBμ

یعنی:

A · B = 5×3 + (-4)×1 = 15 - 4 = 11

باز هم همان نتیجه به‌دست آمد.

نتیجه

در این مثال، ضرب داخلی دو چهاربردار برابر 11 شد، و این مقدار یک اسکالر ناوردا است؛ یعنی اگر به چارچوب لخت دیگری برویم، مقدارش عوض نمی‌شود.

نسخه خیلی کوتاه برای خودِ پست:

Aμ = (5,4,0,0), Bμ = (3,1,0,0)

A·B = 5×3 - 4×1 = 11

پس:

A·B = 11

اگر بخواهی، در پیام بعدی یک مثال کامل‌تر می‌زنم که:

• هم تبدیل لورنتس انجام شود،
• هم نشان دهیم قبل و بعد از تبدیل، A·B ثابت می‌ماند.

ضرب داخلی (Inner Product) چهاربردار در نسبیت خاص — مناسب بلاگفا

۱) تعریف ضرب داخلی در فضای مینکوفسکی

برخلاف فضای اقلیدسی، در نسبیت خاص ضرب داخلی با متریک مینکوفسکی تعریف می‌شود. با قرارداد رایج:

ημν = diag(1, -1, -1, -1)

اگر دو چهاربردار داشته باشیم:

Aμ = (A0, A1, A2, A3)

Bμ = (B0, B1, B2, B3)

ضرب داخلی (اسکالر لورنتسی) آن‌ها این است:

A · B = ημν Aμ Bν

که در این قرارداد به صورت مؤلفه‌ای می‌شود:

A · B = A0B0 - A1B1 - A2B2 - A3B3

این کمیت یک ناوردا است (در همه چارچوب‌های لَخت یک مقدار می‌دهد).

۲) راه خیلی ساده‌تر: پایین آوردن اندیس

اول اندیس یکی از بردارها را پایین می‌آوریم:

Aμ = ημν Aν = (A0, -A1, -A2, -A3)

بعد ضرب داخلی فقط می‌شود:

A · B = Aμ Bμ = A0B0 + A1B1 + A2B2 + A3B3

دقت کنید علامت‌های منفی «داخلِ Aμ» نشسته‌اند.

۳) مثال عددی سریع

فرض کنید:

• Aμ = (2, 1, 0, 0)
• Bμ = (5, 3, 0, 0)

پس:

A · B = 2×5 - 1×3 = 10 - 3 = 7

۴) مربع طول یک چهاربردار (Self-inner product)

اگر A=B، داریم:

A · A = AμAμ = (A0)² - (A1)² - (A2)² - (A3)²

این کمیت (با همین قرارداد علامت‌ها) می‌تواند:

• مثبت باشد → بردار زمان‌مانند (timelike)
• صفر باشد → بردار نورمانند (lightlike/null)
• منفی باشد → بردار فضامانند (spacelike)

۵) دو کاربرد خیلی مهم در فیزیک

الف) چهارسرعت U

برای ذرهٔ مادی:

Uμ = (γc, γvx, γvy, γvz)

و همیشه:

U · U = UμUμ = c²

ب) چهارتکانه P

چهارتکانه:

Pμ = (E/c, px, py, pz)

و ناوردای مشهور انرژی-تکانه:

P · P = PμPμ = (E/c)² - p² = m²c²

که در آن:

p² = px² + py² + pz²

۶) نکتهٔ مهم (قرارداد علامت‌ها)

بعضی کتاب‌ها برعکسِ این قرارداد را می‌گیرند:

η = diag(-1, +1, +1, +1)

در آن صورت علامت‌های ضرب داخلی عوض می‌شود. ولی فیزیک یکی است؛ فقط باید در کل متن یک‌دست بمانید.

اگر دوست داری، برای بلاگفا یک مثال «کامل‌تر» هم می‌نویسم که در آن:

• یک چهاربردار را به چارچوب دیگر تبدیل کنیم،
• بعد نشان بدهیم مقدار A·B قبل و بعد دقیقاً یکی است.

تبدیل لورنتس بدون ماتریس (متن‌محور و مناسب بلاگفا) + مثال عددی

چون بلاگفا معمولاً نمایش ماتریس‌ها را به‌هم می‌ریزد، همان内容 را کاملاً بدون ماتریس و با فرمول‌های خطی می‌نویسم.

۱) تبدیل لورنتس (Boost در راستای x)

دو چارچوب لَخت داریم: S و S'.

فرض کنید S' با سرعت v در راستای +x نسبت به S حرکت می‌کند.

تعریف می‌کنیم:

β = v/c , γ = 1 / √(1-β²)

تبدیل لورنتس برای مختصات (فقط راستای x):

x' = γ (x - v t)

t' = γ (t - v x / c²)

و در راستاهای دیگر:

y' = y , z' = z

اگر c=1 بگیریم (واحد نسبیتی)، این دو رابطه ساده‌تر می‌شوند:

x' = γ (x - β t)

t' = γ (t - β x)

۲) مثال عددی گام‌به‌گام (چهارمکان)

فرض کنید در چارچوب S یک رویداد این مختصات را دارد:

• t = 10
• x = 6
• y = 0
• z = 0

و سرعت نسبی دو چارچوب:

β = 0.6 → γ = 1/√(1-0.36) = 1/0.8 = 1.25

حالا تبدیل می‌دهیم:

محاسبه x':

x' = γ (x - β t) = 1.25 (6 - 0.6×10) = 1.25 (6 - 6) = 0

محاسبه t':

t' = γ (t - β x) = 1.25 (10 - 0.6×6) = 1.25 (10 - 3.6) = 1.25×6.4 = 8

پس در چارچوب S':

• t' = 8
• x' = 0
• y' = 0
• z' = 0

۳) چک ناوردایی فاصله-زمان (بدون ماتریس)

با قرارداد متریک مینکوفسکی (+,-,-,-)، فاصله-زمان:

s² = c² t² - x² - y² - z²

اگر c=1:

s² = t² - x² - y² - z²

در چارچوب S:

s² = 10² - 6² = 100 - 36 = 64

در چارچوب S':

s'² = 8² - 0² = 64

پس:

s'² = s²

یعنی فاصله-زمان ناوردا است.

۴) تبدیل لورنتس برای چهاربردارها (بدون ماتریس)

اگر یک چهاربردار هموردا داشته باشیم:

Vμ = (V0, V1, V2, V3)

در Boost راستای x:

V'0 = γ (V0 - β V1)

V'1 = γ (V1 - β V0)

V'2 = V2

V'3 = V3

(همان محتوای ماتریسی است، فقط به شکل مؤلفه‌ای نوشته شده.)

۵) مثال عددی برای چهارتکانه (بدون ماتریس)

در واحد c=1، چهارتکانه را می‌نویسیم:

Pμ = (E, px, py, pz)

فرض کنید در S:

• E = 10
• px = 8
• py = 0
• pz = 0

با همان β=0.6 و γ=1.25:

انرژی در چارچوب جدید:

E' = P'^0 = γ (E - β px) = 1.25 (10 - 0.6×8) = 6.5

تکانه x در چارچوب جدید:

p'x = P'^1 = γ (px - β E) = 1.25 (8 - 0.6×10) = 2.5

پس:

• E' = 6.5
• p'x = 2.5

جرم ناوردا (بررسی):

m² = E² - p²

اینجا:

m² = 10² - 8² = 36

و در چارچوب جدید:

m² = (6.5)² - (2.5)² = 42.25 - 6.25 = 36

پس جرم ناوردا ثابت می‌ماند.

اگر می‌خواهی دقیقاً با قالب بلاگفا هماهنگ‌تر شود:

بگو «فرمول‌ها را با چه مدلی بهتر می‌بینی؟» یکی از این‌ها را انتخاب کن تا همان را یک‌دست کنم:

• ۱) فقط متن ساده (بدون
  و بدون خط‌های لاتک)
• ۲) همین مدل (HTML + خط‌های لاتین در جهت چپ)
• ۳) فرمول‌ها به صورت عکس (اگر خودت در بلاگفا عکس فرمول می‌گذاری)

تبدیل لورنتس با ماتریس + مثال عددی گام‌به‌گام — مناسب برای بلاگفا

در این پست می‌خواهیم تبدیل لورنتس را به زبان ماتریسی بنویسیم و بعد با یک مثال عددی ببینیم چطور یک چهاربردار از یک چارچوب به چارچوب دیگر می‌رود.

۱) ایده اصلی

در نسبیت خاص، وقتی از یک چارچوب لَخت به چارچوب لَخت دیگری می‌رویم، مختصات و چهاربردارها با یک تبدیل خطی عوض می‌شوند که به آن تبدیل لورنتس می‌گوییم.

اگر یک چهاربردار هموردا داشته باشیم:

Vμ =
\[
\begin{pmatrix}
V^0\\
V^1\\
V^2\\
V^3
\end{pmatrix}
\]

آنگاه در چارچوب جدید:

V'μ = \Lambdaμν Vν

یا به‌صورت ماتریسی:

\[
\begin{pmatrix}
V'^0\\
V'^1\\
V'^2\\
V'^3
\end{pmatrix}
=
\begin{pmatrix}
\gamma & -\beta\gamma & 0 & 0\\
-\beta\gamma & \gamma & 0 & 0\\
0 & 0 & 1 & 0\\
0 & 0 & 0 & 1
\end{pmatrix}
\begin{pmatrix}
V^0\\
V^1\\
V^2\\
V^3
\end{pmatrix}
\]

که این ماتریس مربوط به Boost در راستای x است.

۲) تعریف β و γ

اگر چارچوب S' با سرعت v در راستای x نسبت به S حرکت کند، داریم:

\[
\beta = \frac{v}{c}, \qquad \gamma = \frac{1}{\sqrt{1-\beta^2}}
\]

۳) شکل صریح مؤلفه‌ها

اگر ماتریس بالا را در بردار ضرب کنیم، به روابط آشنا می‌رسیم:

\[
V'^0 = \gamma(V^0 - \beta V^1)
\]
\[
V'^1 = \gamma(V^1 - \beta V^0)
\]
\[
V'^2 = V^2
\]
\[
V'^3 = V^3
\]

یعنی در Boost راستای x فقط مؤلفه‌های زمانی و x با هم مخلوط می‌شوند.

۴) مثال عددی اول: چهارمکان

چهارمکان یک رویداد را به‌صورت زیر می‌نویسیم:

\[
x^\mu = (ct,\; x,\; y,\; z)
\]

برای راحتی واحدهایی می‌گیریم که c=1. پس:

\[
x^\mu = (t,\; x,\; y,\; z)
\]

فرض کنید در چارچوب S داریم:

\[
x^\mu =
\begin{pmatrix}
10\\
6\\
0\\
0
\end{pmatrix}
\]

و سرعت نسبی بین دو چارچوب:

\[
\beta = 0.6
\]

پس:

\[
\gamma = \frac{1}{\sqrt{1-0.6^2}}=\frac{1}{\sqrt{0.64}}=1.25
\]

۵) ساخت ماتریس لورنتس

برای این مقدار از \beta و \gamma، ماتریس لورنتس می‌شود:

\[
\Lambda =
\begin{pmatrix}
1.25 & -0.75 & 0 & 0\\
-0.75 & 1.25 & 0 & 0\\
0 & 0 & 1 & 0\\
0 & 0 & 0 & 1
\end{pmatrix}
\]

چون:

\[
\beta\gamma = 0.6 \times 1.25 = 0.75
\]

۶) ضرب ماتریسی گام‌به‌گام

حالا می‌خواهیم حساب کنیم:

\[
x'^\mu = \Lambda x^\mu
\]

یعنی:

\[
\begin{pmatrix}
x'^0\\
x'^1\\
x'^2\\
x'^3
\end{pmatrix}
=
\begin{pmatrix}
1.25 & -0.75 & 0 & 0\\
-0.75 & 1.25 & 0 & 0\\
0 & 0 & 1 & 0\\
0 & 0 & 0 & 1
\end{pmatrix}
\begin{pmatrix}
10\\
6\\
0\\
0
\end{pmatrix}
\]

حالا سطر به سطر ضرب می‌کنیم:

سطر اول:

\[
x'^0 = (1.25)(10) + (-0.75)(6) = 12.5 - 4.5 = 8
\]

سطر دوم:

\[
x'^1 = (-0.75)(10) + (1.25)(6) = -7.5 + 7.5 = 0
\]

سطر سوم:

\[
x'^2 = 0
\]

سطر چهارم:

\[
x'^3 = 0
\]

پس جواب نهایی:

\[
x'^\mu =
\begin{pmatrix}
8\\
0\\
0\\
0
\end{pmatrix}
\]

۷) تفسیر فیزیکی این نتیجه

در چارچوب اول، رویداد در مختصات (t,x)=(10,6) بود، اما در چارچوب دوم به صورت (8,0) دیده شد.

یعنی ناظر دوم این رویداد را دقیقاً روی مبدأ مکانی خودش دیده است، ولی در زمان متفاوت.

این همان چیزی است که از اختلاط زمان و مکان در نسبیت انتظار داریم.

۸) بررسی ناوردایی با ماتریس متریک

در نسبیت خاص، کمیت زیر باید ثابت بماند:

\[
x^\mu x_\mu = x^T \eta x
\]

که در آن:

\[
\eta =
\begin{pmatrix}
1 & 0 & 0 & 0\\
0 & -1 & 0 & 0\\
0 & 0 & -1 & 0\\
0 & 0 & 0 & -1
\end{pmatrix}
\]

برای بردار اولیه:

\[
x^\mu x_\mu = 10^2 - 6^2 = 100 - 36 = 64
\]

برای بردار نهایی:

\[
x'^\mu x'_\mu = 8^2 - 0^2 = 64
\]

پس کمیت ناوردا حفظ شده است.

۹) مثال دوم: چهارتکانه با ماتریس

حالا دقیقاً همین روش را برای چهارتکانه انجام می‌دهیم.

فرض کنید:

\[
P^\mu =
\begin{pmatrix}
10\\
8\\
0\\
0
\end{pmatrix}
\]

آنگاه:

\[
P'^\mu = \Lambda P^\mu
\]

یعنی:

\[
\begin{pmatrix}
P'^0\\
P'^1\\
P'^2\\
P'^3
\end{pmatrix}
=
\begin{pmatrix}
1.25 & -0.75 & 0 & 0\\
-0.75 & 1.25 & 0 & 0\\
0 & 0 & 1 & 0\\
0 & 0 & 0 & 1
\end{pmatrix}
\begin{pmatrix}
10\\
8\\
0\\
0
\end{pmatrix}
\]

حساب می‌کنیم:

\[
P'^0 = 1.25(10) - 0.75(8) = 12.5 - 6 = 6.5
\]
\[
P'^1 = -0.75(10) + 1.25(8) = -7.5 + 10 = 2.5
\]

پس:

\[
P'^\mu =
\begin{pmatrix}
6.5\\
2.5\\
0\\
0
\end{pmatrix}
\]

و جرم ناوردا:

\[
P^\mu P_\mu = 10^2 - 8^2 = 36
\]
\[
P'^\mu P'_\mu = 6.5^2 - 2.5^2 = 42.25 - 6.25 = 36
\]

باز هم ثابت ماند.

۱۰) یک رابطه خیلی مهم ماتریسی

تمام دلیل اینکه کمیت‌های ناوردا ثابت می‌مانند، این است که ماتریس لورنتس شرط زیر را ارضا می‌کند:

\[
\Lambda^T \eta \Lambda = \eta
\]

این معادله یعنی تبدیل لورنتس متریک مینکوفسکی را حفظ می‌کند.

به زبان ساده:
لورنتس تبدیل‌هایی هستند که فاصله-زمان را تغییر نمی‌دهند.

جمع‌بندی نهایی

• هر چهاربردار هموردا با رابطه V'^\mu=\Lambda^\mu_{\ \nu}V^\nu تبدیل می‌شود.


• برای Boost در راستای x فقط مؤلفه‌های زمانی و x با هم ترکیب می‌شوند.


• ضرب ماتریسی، همان فرمول‌های معمول تبدیل لورنتس را تولید می‌کند.


• کمیت‌هایی مثل x^\mu x_\mu و P^\mu P_\mu ناوردا هستند.


• علت اصلی ناوردایی این است که \Lambda^T \eta \Lambda = \eta.

اگر بخواهی، در پیام بعدی برای بلاگفا یکی از این‌ها را هم آماده می‌کنم:

• چهارسرعت و ناوردایی آن


• رپیدیتی و سرعت ویژه


• بردار کوواریانت و پایین‌آوردن اندیس با مثال عددی

تبدیل لورنتس با ماتریس + مثال عددی گام‌به‌گام — مناسب برای بلاگفا

در این پست می‌خواهیم تبدیل لورنتس را به زبان ماتریسی بنویسیم و بعد با یک مثال عددی ببینیم چطور یک چهاربردار از یک چارچوب به چارچوب دیگر می‌رود.

۱) ایده اصلی

در نسبیت خاص، وقتی از یک چارچوب لَخت به چارچوب لَخت دیگری می‌رویم، مختصات و چهاربردارها با یک تبدیل خطی عوض می‌شوند که به آن تبدیل لورنتس می‌گوییم.

اگر یک چهاربردار هموردا داشته باشیم:

Vμ =
\[
\begin{pmatrix}
V^0\\
V^1\\
V^2\\
V^3
\end{pmatrix}
\]

آنگاه در چارچوب جدید:

V'μ = \Lambdaμν Vν

یا به‌صورت ماتریسی:

\[
\begin{pmatrix}
V'^0\\
V'^1\\
V'^2\\
V'^3
\end{pmatrix}
=
\begin{pmatrix}
\gamma & -\beta\gamma & 0 & 0\\
-\beta\gamma & \gamma & 0 & 0\\
0 & 0 & 1 & 0\\
0 & 0 & 0 & 1
\end{pmatrix}
\begin{pmatrix}
V^0\\
V^1\\
V^2\\
V^3
\end{pmatrix}
\]

که این ماتریس مربوط به Boost در راستای x است.

۲) تعریف β و γ

اگر چارچوب S' با سرعت v در راستای x نسبت به S حرکت کند، داریم:

\[
\beta = \frac{v}{c}, \qquad \gamma = \frac{1}{\sqrt{1-\beta^2}}
\]

۳) شکل صریح مؤلفه‌ها

اگر ماتریس بالا را در بردار ضرب کنیم، به روابط آشنا می‌رسیم:

\[
V'^0 = \gamma(V^0 - \beta V^1)
\]
\[
V'^1 = \gamma(V^1 - \beta V^0)
\]
\[
V'^2 = V^2
\]
\[
V'^3 = V^3
\]

یعنی در Boost راستای x فقط مؤلفه‌های زمانی و x با هم مخلوط می‌شوند.

۴) مثال عددی اول: چهارمکان

چهارمکان یک رویداد را به‌صورت زیر می‌نویسیم:

\[
x^\mu = (ct,\; x,\; y,\; z)
\]

برای راحتی واحدهایی می‌گیریم که c=1. پس:

\[
x^\mu = (t,\; x,\; y,\; z)
\]

فرض کنید در چارچوب S داریم:

\[
x^\mu =
\begin{pmatrix}
10\\
6\\
0\\
0
\end{pmatrix}
\]

و سرعت نسبی بین دو چارچوب:

\[
\beta = 0.6
\]

پس:

\[
\gamma = \frac{1}{\sqrt{1-0.6^2}}=\frac{1}{\sqrt{0.64}}=1.25
\]

۵) ساخت ماتریس لورنتس

برای این مقدار از \beta و \gamma، ماتریس لورنتس می‌شود:

\[
\Lambda =
\begin{pmatrix}
1.25 & -0.75 & 0 & 0\\
-0.75 & 1.25 & 0 & 0\\
0 & 0 & 1 & 0\\
0 & 0 & 0 & 1
\end{pmatrix}
\]

چون:

\[
\beta\gamma = 0.6 \times 1.25 = 0.75
\]

۶) ضرب ماتریسی گام‌به‌گام

حالا می‌خواهیم حساب کنیم:

\[
x'^\mu = \Lambda x^\mu
\]

یعنی:

\[
\begin{pmatrix}
x'^0\\
x'^1\\
x'^2\\
x'^3
\end{pmatrix}
=
\begin{pmatrix}
1.25 & -0.75 & 0 & 0\\
-0.75 & 1.25 & 0 & 0\\
0 & 0 & 1 & 0\\
0 & 0 & 0 & 1
\end{pmatrix}
\begin{pmatrix}
10\\
6\\
0\\
0
\end{pmatrix}
\]

حالا سطر به سطر ضرب می‌کنیم:

سطر اول:

\[
x'^0 = (1.25)(10) + (-0.75)(6) = 12.5 - 4.5 = 8
\]

سطر دوم:

\[
x'^1 = (-0.75)(10) + (1.25)(6) = -7.5 + 7.5 = 0
\]

سطر سوم:

\[
x'^2 = 0
\]

سطر چهارم:

\[
x'^3 = 0
\]

پس جواب نهایی:

\[
x'^\mu =
\begin{pmatrix}
8\\
0\\
0\\
0
\end{pmatrix}
\]

۷) تفسیر فیزیکی این نتیجه

در چارچوب اول، رویداد در مختصات (t,x)=(10,6) بود، اما در چارچوب دوم به صورت (8,0) دیده شد.

یعنی ناظر دوم این رویداد را دقیقاً روی مبدأ مکانی خودش دیده است، ولی در زمان متفاوت.

این همان چیزی است که از اختلاط زمان و مکان در نسبیت انتظار داریم.

۸) بررسی ناوردایی با ماتریس متریک

در نسبیت خاص، کمیت زیر باید ثابت بماند:

\[
x^\mu x_\mu = x^T \eta x
\]

که در آن:

\[
\eta =
\begin{pmatrix}
1 & 0 & 0 & 0\\
0 & -1 & 0 & 0\\
0 & 0 & -1 & 0\\
0 & 0 & 0 & -1
\end{pmatrix}
\]

برای بردار اولیه:

\[
x^\mu x_\mu = 10^2 - 6^2 = 100 - 36 = 64
\]

برای بردار نهایی:

\[
x'^\mu x'_\mu = 8^2 - 0^2 = 64
\]

پس کمیت ناوردا حفظ شده است.

۹) مثال دوم: چهارتکانه با ماتریس

حالا دقیقاً همین روش را برای چهارتکانه انجام می‌دهیم.

فرض کنید:

\[
P^\mu =
\begin{pmatrix}
10\\
8\\
0\\
0
\end{pmatrix}
\]

آنگاه:

\[
P'^\mu = \Lambda P^\mu
\]

یعنی:

\[
\begin{pmatrix}
P'^0\\
P'^1\\
P'^2\\
P'^3
\end{pmatrix}
=
\begin{pmatrix}
1.25 & -0.75 & 0 & 0\\
-0.75 & 1.25 & 0 & 0\\
0 & 0 & 1 & 0\\
0 & 0 & 0 & 1
\end{pmatrix}
\begin{pmatrix}
10\\
8\\
0\\
0
\end{pmatrix}
\]

حساب می‌کنیم:

\[
P'^0 = 1.25(10) - 0.75(8) = 12.5 - 6 = 6.5
\]
\[
P'^1 = -0.75(10) + 1.25(8) = -7.5 + 10 = 2.5
\]

پس:

\[
P'^\mu =
\begin{pmatrix}
6.5\\
2.5\\
0\\
0
\end{pmatrix}
\]

و جرم ناوردا:

\[
P^\mu P_\mu = 10^2 - 8^2 = 36
\]
\[
P'^\mu P'_\mu = 6.5^2 - 2.5^2 = 42.25 - 6.25 = 36
\]

باز هم ثابت ماند.

۱۰) یک رابطه خیلی مهم ماتریسی

تمام دلیل اینکه کمیت‌های ناوردا ثابت می‌مانند، این است که ماتریس لورنتس شرط زیر را ارضا می‌کند:

\[
\Lambda^T \eta \Lambda = \eta
\]

این معادله یعنی تبدیل لورنتس متریک مینکوفسکی را حفظ می‌کند.

به زبان ساده:
لورنتس تبدیل‌هایی هستند که فاصله-زمان را تغییر نمی‌دهند.

جمع‌بندی نهایی

• هر چهاربردار هموردا با رابطه V'^\mu=\Lambda^\mu_{\ \nu}V^\nu تبدیل می‌شود.


• برای Boost در راستای x فقط مؤلفه‌های زمانی و x با هم ترکیب می‌شوند.


• ضرب ماتریسی، همان فرمول‌های معمول تبدیل لورنتس را تولید می‌کند.


• کمیت‌هایی مثل x^\mu x_\mu و P^\mu P_\mu ناوردا هستند.


• علت اصلی ناوردایی این است که \Lambda^T \eta \Lambda = \eta.

اگر بخواهی، در پیام بعدی برای بلاگفا یکی از این‌ها را هم آماده می‌کنم:

• چهارسرعت و ناوردایی آن


• رپیدیتی و سرعت ویژه


• بردار کوواریانت و پایین‌آوردن اندیس با مثال عددی

بردار هموردا و کوواریانت (بالا/پایین اندیس) + مثال عددی — مناسب برای بلاگفا

۱) اول یک نکته مهم درباره واژه‌ها

در نسبیت (و کلاً فیزیک ریاضی) سه اصطلاح را قاطی نکنیم:

• هموردا (Contravariant): مؤلفه‌ها با اندیس بالا نوشته می‌شوند: Vμ


• کوواریانت (Covariant): مؤلفه‌ها با اندیس پایین نوشته می‌شوند: Vμ


• ناوردا (Invariant): کمیتی که تحت تبدیل لورنتس ثابت می‌ماند، مثل VμVμ

۲) متریک مینکوفسکی (برای بالا/پایین بردن اندیس)

برای اینکه از هموردا به کوواریانت برویم (یا برعکس)، از متریک استفاده می‌کنیم. من این قرارداد رایج را می‌گیرم:

ημν = diag(1, -1, -1, -1)

پس اگر:

Vμ = (V0, V1, V2, V3)

آنگاه:

Vμ = ημν Vν = (V0, -V1, -V2, -V3)

یعنی با این قرارداد: مؤلفه زمانی ثابت می‌ماند و مؤلفه‌های مکانی علامتشان عوض می‌شود.

۳) تبدیل لورنتس برای حرکت در راستای x (Boost)

اگر چارچوب S' با سرعت v در راستای x نسبت به S حرکت کند، تعریف می‌کنیم:

β = v/c , γ = 1 / √(1-β²)

برای هر چهاربردار هموردا:

V'μ = Λμν Vν

که در این Boost (در راستای x) یعنی:

V'0 = γ (V0 - β V1)
V'1 = γ (V1 - β V0)
V'2 = V2
V'3 = V3

برای بردار کوواریانت تبدیل با معکوس انجام می‌شود، ولی راه ساده این است:

اول V'μ را حساب کن، بعد اندیس را پایین بیاور و V'μ را بساز.

۴) مثال کامل و عددی (چهارمکان)

چهارمکان یک رویداد:

xμ = (ct, x, y, z)

برای سادگی واحدهایی می‌گیریم که c=1، پس x0=t.

فرض کنید در چارچوب S رویداد:

xμ = (10, 6, 0, 0)

و سرعت نسبی:

β = 0.6 → γ = 1/√(1-0.36) = 1.25

تبدیل هموردا:

x'0 = 1.25(10 - 0.6×6) = 1.25(6.4) = 8
x'1 = 1.25(6 - 0.6×10) = 1.25(0) = 0

پس:

x'μ = (8, 0, 0, 0)

کوواریانتِ همان بردار:

xμ = (10, -6, 0, 0)
x'μ = (8, 0, 0, 0)

۵) چک ناوردایی (Invariant) فاصله-زمان

کمیت زیر ناوردا است:

s² = xμ xμ

در چارچوب S:

xμxμ = 10×10 + 6×(-6) = 100 - 36 = 64

در چارچوب S':

x'μx'μ = 8×8 + 0 = 64

نتیجه: مقدار xμxμ در هر دو چارچوب یکی است؛ یعنی ناورداست.

۶) مثال دوم (چهارتکانه) + جرم ناوردا

چهارتکانه:

Pμ = (E, px, py, pz) (با c=1)

فرض کنید:

Pμ = (10, 8, 0, 0)

پس:

Pμ = (10, -8, 0, 0)

کمیت ناوردای جرم:

m² = PμPμ = 10² - 8² = 36 → m = 6

حالا به همان چارچوب با β=0.6 برویم:

P'0 = 1.25(10 - 0.6×8) = 6.5
P'1 = 1.25(8 - 0.6×10)= 2.5

پس:

P'μ = (6.5, 2.5, 0, 0)
P'μ = (6.5, -2.5, 0, 0)

و دوباره:

P'μP'μ = 6.5² - 2.5² = 42.25 - 6.25 = 36

پس m ثابت ماند.

جمع‌بندی خیلی کوتاه

• هموردا: Vμ و تبدیلش با Λ


• کوواریانت: Vμ که با متریک از روی هموردا ساخته می‌شود


• ناوردا: چیزهایی مثل VμVμ (عدد ثابت در همه چارچوب‌ها)

اگر دوست داری، در قالب بلاگفا یک پست جدا هم می‌نویسم که فقط «ماتریس لورنتس» را با مثال عددی (گام‌به‌گام) توضیح بدهد.

بردارهای ناوردا (Covariant) و هم‌وردا (Contravariant) در نسبیت خاص

در فضازمان مینکوفسکی (که اساس نسبیت خاص است)، ما با «چارچوب‌های مرجع» مختلفی سروکار داریم که نسبت به هم با سرعت ثابت حرکت می‌کنند. اشیاء فیزیکی (مثل چهاربردارها) باید طوری توصیف شوند که قوانین فیزیک در همهٔ این چارچوب‌ها یکسان باقی بمانند. اینجا نقش «بردارهای ناوردا» و «هم‌وردا» پررنگ می‌شود.

۱) چرا به این مفاهیم نیاز داریم؟

یک بردار معمولی (مثلاً در فضای سه‌بعدی اقلیدسی) را در نظر بگیرید. اگر شما دستگاه مختصات را بچرخانید، مؤلفه‌های بردار عوض می‌شوند، ولی خودِ بردار (طول و جهت آن) ثابت می‌ماند. در نسبیت، قضیه از این هم پیچیده‌تر است؛ چون علاوه بر چرخش، «تبدیل لورنتس» (که شامل افزایش سرعت است) را هم داریم.

تبدیل لورنتس، مختصات فضا-زمان را از یک چارچوب S به چارچوب S' تغییر می‌دهد. مؤلفه‌های یک «چهاربردار» تحت این تبدیل، به شیوه‌ای خاص تغییر می‌کنند که این شیوه، بین بردارهای «هم‌وردا» و «ناوردا» فرق دارد.

۲) تبدیل لورنتس و ماتریس آن

فرض کنید S و S' دو چارچوب باشند که S' با سرعت v در راستای محور x نسبت به S حرکت می‌کند.

یک رویداد (event) در فضا-زمان با مختصات (x⁰, x¹, x², x³) = (ct, x, y, z) نمایش داده می‌شود.

تبدیل لورنتس (برای حرکت در راستای x) به شکل ماتریسی:

x'μ = Σν=03 Λμν xν

که ماتریس Λ (لامبدا) این است (با β = v/c و γ = 1/√(1-β²)):

Λ =

⎛ γ -βγ 0 0 ⎞

⎜ -βγ γ 0 0 ⎟

⎜ 0 0 1 0 ⎟

⎝ 0 0 0 1 ⎠

این یعنی:

ct' = γ(ct - βx)

x' = γ(x - βct)

y' = y

z' = z

۳) بردارهای هم‌وردا (Contravariant Vectors)

این‌ها همان «چهاربردار»هایی هستند که معمولاً با آن‌ها کار می‌کنیم (مانند چهارمکان، چهارسرعت، چهارتکانه). مؤلفه‌های یک چهاربردار هم‌وردا، مثلاً Vμ، تحت تبدیل لورنتس به صورت زیر تبدیل می‌شوند:

V'μ = Σν=03 Λμν Vν

این یعنی مؤلفه‌های V' با استفاده از همان ماتریس Λ از مؤلفه‌های V به دست می‌آیند.

مثال: چهارمکان xμ = (ct, x, y, z) یک بردار هم‌وردا است.

۴) بردارهای ناوردا (Covariant Vectors)

این بردارها «دوگان» (dual) بردارهای هم‌وردا هستند. مؤلفه‌های یک بردار ناوردا، مثلاً Wμ، تحت تبدیل لورنتس به صورت زیر تبدیل می‌شوند:

W'μ = Σν=03 (Λ-1)νμ Wν

که Λ-1 ماتریس معکوس Λ است.

برای ماتریس تبدیل لورنتس (در راستای x)، ماتریس معکوس آن (یعنی تبدیل از S' به S) با تغییر علامت v (یا β) به دست می‌آید:

Λ-1 =

⎛ γ βγ 0 0 ⎞

⎜ βγ γ 0 0 ⎟

⎜ 0 0 1 0 ⎟

⎝ 0 0 0 1 ⎠

پس تبدیل برای مؤلفه‌های ناوردا این است:

W'⁰ = γ(W₀ + βW₁)

W'¹ = γ(βW₀ + W₁)

W'² = W₂

W'³ = W₃

۵) متریک فضازمان (Minkowski Metric)

برای رفتن از مختصات هم‌وردا به ناوردا (و بالعکس)، از «تنسور متریک» ημν استفاده می‌کنیم. در مختصات معمول کارتزین:

ημν =

⎛ -1 0 0 0 ⎞

⎜ 0 1 0 0 ⎟

⎜ 0 0 1 0 ⎟

⎝ 0 0 0 1 ⎠

(گاهی تعریف η00 = +1 هم رایج است؛ مهم این است که قرارداد ثابت باشد).

با استفاده از ημν داریم:

- تبدیل از مؤلفهٔ بالانویس (هم‌وردا) به زیرنویس (ناوردا):

Wμ = Σν=03 ημν Wν

مثال: W₀ = -W⁰, W₁ = W¹, W₂ = W², W₃ = W³

- تبدیل از مؤلفهٔ زیرنویس (ناوردا) به بالانویس (هم‌وردا):

Wμ = Σν=03 ημν Wν

که ημν ماتریس معکوس ημν است، که در این مورد با خودش برابر است (ماتریس قطری با ۱، ۱، ۱، ۱-).

۶) حاصلضرب داخلی (Inner Product)

حاصلضرب داخلی یک بردار هم‌وردا با یک بردار ناوردا، یک «اسکالر» (عدد تنها) است که تحت تبدیل لورنتس تغییر نمی‌کند (ناوردا است).

S = V ⋅ W = Σμ=03 Vμ Wμ

= Σμ,ν Vμ ημν Wν

اگر دو بردار هم‌وردا داشته باشیم Vμ و Uμ، حاصلضرب داخلی آن‌ها (که یک اسکالر ناوردا است) به صورت زیر تعریف می‌شود:

V ⋅ U = Σμ,ν ημν Vμ Uν

این همان چیزی است که به عنوان «جرم ناوردای» یک ذره یا سیستم محاسبه می‌کردیم:

(m c)² = P ⋅ P / c² = (Pμ Pμ) / c²

یا با c=1:

m² = P ⋅ P = Pμ Pμ = Σμ Pμ Pμ

خلاصه:

• هم‌وردا (Contravariant): مؤلفه‌ها تحت تبدیل لورنتس Λ تغییر می‌کنند. (مثال: xμ)
• ناوردا (Covariant): مؤلفه‌ها تحت تبدیل معکوس Λ-1 تغییر می‌کنند. (مثال: ∂μ، گرادیان)
• متریک ημν: ابزاری برای تبدیل بین این دو نوع بردار و محاسبهٔ اسکالرهای ناوردا.
• اسکالرها/ناورداها: کمیت‌هایی که تحت تبدیل لورنتس تغییر نمی‌کنند (مانند جرم ناوردا، فاصله‌زمان ناوردا).

این تمایز برای اطمینان از اینکه فیزیک در تمام چارچوب‌های لخت یکسان بیان می‌شود، حیاتی است.

مثال‌های واپاشی دو/سه‌ذره‌ای با محاسبهٔ انرژی‌ها در «چارچوب مرکز جرم / سکون ذرهٔ مادر»

فرمول‌های طلایی (برای استفاده در همه مثال‌ها)

فرض کنید ذرهٔ مادر A در چارچوب سکون خودش (همان مرکز جرمِ واپاشی) واپاشی کند.

A → 1 + 2

در چارچوب سکون A داریم:

E1 = (m_A^2 + m_1^2 - m_2^2) c^2 / (2 m_A)
E2 = (m_A^2 + m_2^2 - m_1^2) c^2 / (2 m_A)

و اندازهٔ تکانه‌ها برابر است:

|p1| = |p2| = √(E1^2/c^2 - m_1^2 c^2) = √(E2^2/c^2 - m_2^2 c^2)

در واحدهای رایج ذره‌ای (c=1):

E1 = (m_A^2 + m_1^2 - m_2^2)/(2m_A)
E2 = (m_A^2 + m_2^2 - m_1^2)/(2m_A)
p = √(E1^2 - m_1^2) = √(E2^2 - m_2^2)

برای واپاشی سه‌ذره‌ای:

A → 1 + 2 + 3

در مرکز جرم، انرژی‌ها «تک‌مقداری ثابت» نیستند و یک بازه دارند؛ اما در حالت‌های خاص (مثلاً هر سه ذره بی‌جرم، یا دو ذره را به صورت یک «سیستم مؤثر» جمع کنیم) می‌توان مقدار/بازه را دقیق حساب کرد.

مثال ۱) واپاشی دوذره‌ای: π⁰ → γ + γ (دو فوتون)

می‌دانیم:

m_{π0} c^2 ≈ 135 MeV
m_γ = 0

در چارچوب سکون π⁰، دو فوتون باید پشت‌به‌پشت باشند و چون جرمشان صفر است، انرژی‌ها برابر می‌شوند:

E_γ = (m_{π0} c^2)/2 ≈ 67.5 MeV

تکانهٔ هر فوتون:

p_γ = E_γ/c → (با c=1: p_γ = E_γ)

نتیجه: در مرکز جرم π⁰ هر فوتون انرژی 67.5 MeV دارد.

مثال ۲) واپاشی دوذره‌ای با جرم‌های متفاوت: K⁺ → μ⁺ + ν_μ

(برای ساده‌سازی، جرم نوترینو را صفر می‌گیریم.)

جرم‌ها:

m_K c^2 ≈ 493.7 MeV
m_μ c^2 ≈ 105.7 MeV
m_ν ≈ 0

در چارچوب سکون K:

E_μ = (m_K^2 + m_μ^2 - m_ν^2)/(2m_K)
= (m_K^2 + m_μ^2)/(2m_K)

عددگذاری:

m_K^2 ≈ 493.7^2 ≈ 243,740
m_μ^2 ≈ 105.7^2 ≈ 11,170
جمع ≈ 254,910
2m_K ≈ 987.4

پس:

E_μ ≈ 254,910 / 987.4 ≈ 258.2 MeV

انرژی نوترینو:

E_ν = (m_K^2 - m_μ^2)/(2m_K)
≈ (243,740 - 11,170)/987.4
≈ 232,570 / 987.4
≈ 235.6 MeV

تکانه‌ها (در c=1):

p = √(E_μ^2 - m_μ^2)
≈ √(258.2^2 - 105.7^2)
≈ √(66,670 - 11,170)
≈ √55,500
≈ 235.6 MeV

نتیجه:

در مرکز جرم K، میون انرژی حدود 258 MeV دارد و نوترینو حدود 236 MeV، و اندازهٔ تکانهٔ هر دو یکسان است.

مثال ۳) واپاشی دوذره‌ای معروف: Z⁰ → e⁺ + e⁻

جرم Z:

m_Z c^2 ≈ 91.19 GeV

جرم الکترون در مقیاس GeV ناچیز است:

m_e c^2 ≈ 0.000511 GeV (قابل صرف‌نظر)

پس تقریباً:

E_{e+} ≈ E_{e-} ≈ m_Z/2 ≈ 45.6 GeV
p ≈ E (چون m خیلی کوچک است)

نتیجه: در مرکز جرم Z، هرکدام از e⁺ و e⁻ تقریباً 45.6 GeV انرژی دارند.

مثال ۴) واپاشی سه‌ذره‌ای: μ⁻ → e⁻ + ν̄_e + ν_μ (طیف انرژی الکترون)

این واپاشی سه‌جسمی است، پس انرژی e ثابت نیست و طیف دارد.

برای اینکه «صرفاً سینماتیکی» یک بازهٔ ممکن به دست آوریم (بدون دینامیک برهم‌کنش ضعیف)، جرم نوترینوها را صفر و جرم الکترون را هم (تقریباً) صفر می‌گیریم.

در چارچوب سکون میون:

m_μ c^2 ≈ 105.7 MeV

حداکثر انرژی الکترون زمانی است که دو نوترینو در یک جهت حرکت کنند و تکانهٔ مؤثرشان «تا حد ممکن» مخالف الکترون شود؛ نتیجهٔ سینماتیکی:

E_e,max ≈ (m_μ c^2)/2 ≈ 52.85 MeV

حداقل انرژی الکترون (در تقریب m_e≈0) می‌تواند به نزدیک صفر برسد:

E_e,min ≈ 0

پس:

0 ≤ E_e ≤ 52.85 MeV

نکته: اگر جرم الکترون را دقیق لحاظ کنیم، حداقل انرژی دقیقاً m_e c^2 است و سقف کمی تغییر می‌کند، ولی مرتبهٔ همان ~53 MeV می‌ماند.

مثال ۵) واپاشی سه‌ذره‌ای با «ذرهٔ مؤثر»: A → 1 + 2 + 3 و محاسبهٔ بازهٔ انرژی 1

یک تکنیک خیلی کاربردی: دو ذره را با هم به صورت یک سیستم مؤثر در نظر می‌گیریم:

X ≡ (2+3)
پس: A → 1 + X

در چارچوب سکون A، این مثل واپاشی دوذره‌ای است، با این تفاوت که m_X ثابت نیست و بین یک حداقل و حداکثر تغییر می‌کند.

اگر جرم‌ها را m_A, m_1, m_2, m_3 بگیریم، داریم:

E1 = (m_A^2 + m_1^2 - m_X^2)/(2m_A)

و بازهٔ ممکن برای جرم مؤثر X:

m_X,min = m_2 + m_3
m_X,max = m_A - m_1

پس انرژی 1 بین دو حد قرار می‌گیرد:

- وقتی m_X کمینه است، E1 بیشینه می‌شود:

E1,max = (m_A^2 + m_1^2 - (m_2+m_3)^2)/(2m_A)

- وقتی m_X بیشینه است، E1 کمینه می‌شود:

E1,min = (m_A^2 + m_1^2 - (m_A - m_1)^2)/(2m_A)

که دومی ساده می‌شود:

E1,min = m_1

(در واحد c=1) یعنی در چارچوب سکون A، کمترین انرژی ذرهٔ 1 همان انرژی سکون خودش است، وقتی تقریباً بی‌حرکت تولید شود.

یک عددگذاری کوچک برای جا افتادن:

فرض کنید:

m_A = 10
m_1 = 2
m_2 = 3
m_3 = 1 (همه در واحد GeV/c² و با c=1)

پس:

m_X,min = 3+1 = 4
→ E1,max = (10^2 + 2^2 - 4^2)/(2×10) = (100+4-16)/20 = 88/20 = 4.4 GeV

E1,min = m_1 = 2 GeV

نتیجه:

2 ≤ E1 ≤ 4.4 (GeV)

اگر بگویی دقیقاً کدام سبک را بیشتر می‌خواهی، ادامه می‌دهم:

• کاملاً عددی و واقعی (با جرم‌های π, K, μ, τ, Λ, …)


• مرحله‌به‌مرحله با چهارتکانه و استخراج فرمول‌ها


• تستی/چهارگزینه‌ای از همین واپاشی‌ها

۵ مثالِ دیگر (متنوع‌تر و کمی عمیق‌تر) از ناوردایی چهارتکانه

مثال ۷) محاسبهٔ انرژی آستانه برای تولید ذرهٔ جدید (کاربرد مستقیمِ ناوردای s)

فرآیند:

p + p (at rest) → p + p + π⁰

یک پروتون در آزمایشگاه به پروتونِ ساکن برخورد می‌کند. می‌خواهیم «کمترین انرژی لازمِ پرتابه» را پیدا کنیم. ایده: کمیت ناوردا

s = (P₁ + P₂)·(P₁ + P₂)

در آستانه، در چارچوب مرکز جرم، محصولات «تا حد ممکن ساکن» هستند؛ بنابراین:

√s_threshold = (2m_p + m_π) c²

در چارچوب آزمایشگاه:
- پروتون هدف ساکن است: P₂ = (m_pc, 0)
- پروتون پرتابه: P₁ = (E/c, p)

پس:

s = (E + m_pc²)²/c²? (بهتر: با انرژی‌ها)
s c⁴ = (E + m_pc²)² - (pc)²

ولی چون برای پرتابه:

E² - (pc)² = m_p² c⁴

می‌گیریم:

s c⁴ = m_p²c⁴ + m_p²c⁴ + 2m_pc² E
= 2m_p²c⁴ + 2m_pc² E

در آستانه:

s_threshold c⁴ = (2m_p + m_π)² c⁴

پس:

2m_p²c⁴ + 2m_pc² E_th = (2m_p + m_π)² c⁴
→ E_th = [ (2m_p + m_π)² - 2m_p² ] c² / (2m_p)

عددگذاری (واحد MeV، با c=1):

m_p = 938 MeV
m_π0 = 135 MeV

(2m_p + m_π) = 2×938 + 135 = 2011 MeV
(2m_p + m_π)² = 2011² ≈ 4,044,121
2m_p² = 2×938² = 2×879,844 ≈ 1,759,688

E_th ≈ (4,044,121 - 1,759,688) / (2×938) MeV
≈ 2,284,433 / 1,876
≈ 1,217 MeV

نتیجه: انرژی کل پرتابه در آستانه حدود

E_th ≈ 1.22 GeV

است. (اگر «انرژی جنبشی» بخواهیم: K_th = E_th - m_p ≈ 279 MeV.)

مثال ۸) از ناوردای چهارتکانه، زاویهٔ برخورد هم مهم می‌شود (دو ذرهٔ پرانرژی با زاویه)

دو ذرهٔ بی‌جرم (مثلاً دو فوتون) با انرژی‌های زیر داریم:

E₁ = 3 MeV
E₂ = 5 MeV

و زاویه بین بردارهای تکانه‌شان θ است.

برای ذرات بی‌جرم:

|p| = E/c

جرم ناوردای سیستم:

M²c^4 = (P₁+P₂)² c⁴ = 2E₁E₂(1 - cosθ)

اگر روبه‌رو باشند θ=π:

1 - cosπ = 2
→ M²c^4 = 4E₁E₂
→ Mc² = 2√(E₁E₂) = 2√(15) ≈ 7.746 MeV

اگر هم‌جهت باشند θ=0:

1 - cos0 = 0
→ M = 0

نتیجه: حتی با انرژی‌های ثابت، فقط با تغییر زاویه، جرم ناوردای سیستم از ۰ تا یک مقدار بزرگ تغییر می‌کند.

مثال ۹) حداقل انرژی فوتون برای تولید زوج e⁺e⁻ کنار یک هسته (تقریب آستانه)

به طور ساده، شرط لازم (از دید انرژی) برای تولید زوج:

γ → e⁺ + e⁻

در خلأ ممکن نیست (به‌خاطر پایستگی چهارتکانه). ولی اگر یک هستهٔ سنگین حضور داشته باشد:

γ + N → e⁺ + e⁻ + N

هسته می‌تواند تکانهٔ لازم را بگیرد.

آستانهٔ انرژی فوتون (وقتی جرم هسته خیلی بزرگ است و پس‌زنی کم است) تقریباً:

E_γ,th ≈ 2 m_e c²

با:

m_e c² = 0.511 MeV

پس:

E_γ,th ≈ 1.022 MeV

نکتهٔ ارتباط با ناوردایی: در خلأ چون فوتون P·P=0 دارد، نمی‌تواند به دو ذرهٔ پرجرم با P_tot·P_tot > 0 تبدیل شود مگر اینکه جسم سومی در تبادل چهارتکانه نقش داشته باشد.

مثال ۱۰) «جرم ناوردای سیستم» همیشه جمع جرم‌ها نیست

دو ذرهٔ هم‌جرم با

m = 1 GeV/c²

را در نظر بگیر که در چارچوب آزمایشگاه هر کدام انرژی

E = 3 GeV

دارند و روبه‌روی هم حرکت می‌کنند (تکانه‌ها مخالف هم).

برای هر ذره:

p c = √(E² - m²) = √(9 - 1) = √8 ≈ 2.828 GeV

برای کل سیستم:

E_tot = 6 GeV
p_tot = 0
→ Mc² = E_tot = 6 GeV

ولی جمع جرم‌های سکون فقط:

m_1 + m_2 = 2 GeV/c²

نتیجه: جرم ناوردای سیستم (۶) خیلی بزرگ‌تر از جمع جرم‌های سکون (۲) است، چون انرژی جنبشی هم در «جرم ناوردای سیستم» سهم دارد.

مثال ۱۱) چک ناوردایی با «تغییر چارچوب» (یک عددگیری واقعی)

ذره‌ای با جرم سکون

m = 4 GeV/c²

در چارچوب S انرژی و تکانه‌اش:

E = 10 GeV
p = 8 GeV/c

ناوردای جرم:

m = √(E² - (pc)²) = √(100 - 64) = 6 GeV/c²

می‌بینید که این با m=4 سازگار نیست؛ پس این داده‌ها برای «یک ذره با جرم ۴» نمی‌تواند درست باشد!

حالا یک دادهٔ سازگار می‌سازیم: اگر m=4 و p=8 باشد، انرژی باید:

E = √( (pc)² + (mc²)² ) = √(8² + 4²) = √80 ≈ 8.944 GeV

حالا فرض کن به چارچوبی برویم که با سرعت v=0.6c حرکت می‌کند (γ=1.25). تبدیل انرژی-تکانه:

E' = γ(E - v p c?) (با c=1: E' = γ(E - β p))
p' = γ(p - β E)

با β=0.6:

E' = 1.25(8.944 - 0.6×8) = 1.25(8.944 - 4.8) = 5.180 GeV
p' = 1.25(8 - 0.6×8.944) = 1.25(8 - 5.366) = 3.292 GeV/c

چک ناوردا:

√(E'² - p'²) = √(5.180² - 3.292²)
≈ √(26.83 - 10.84)
≈ √15.99
≈ 3.999 ≈ 4 GeV

نتیجه: با عوض شدن چارچوب، E و p عوض شدند، اما m (یعنی ناوردای چهارتکانه) ثابت ماند.

اگر بگویی «مثال‌ها را بیشتر ذره‌ای می‌خواهی یا کلاسیکی» می‌توانم دقیق‌تر هدف‌گیری کنم؛ مثلاً:

• مثال‌های کاملاً تستی (چهارگزینه‌ای)


• مثال‌های واپاشی دو/سه‌ذره‌ای با محاسبهٔ انرژی‌ها در چارچوب مرکز جرم


• مثال‌های پرتاب از دید دو ناظر برای ارتباط مستقیم با تبدیل لورنتس

۳ مثالِ دیگر از ناوردایی چهارتکانه

مثال ۴) چارچوبِ سکونِ سیستم را با ناوردا پیدا کن

فرض کنید دو ذره در یک بُعد حرکت می‌کنند:

Particle 1:

E₁ = 8 GeV

p₁ = +6 GeV/c

Particle 2:

E₂ = 5 GeV

p₂ = -3 GeV/c

چهارتکانهٔ کل سیستم:

E_tot = E₁ + E₂ = 13 GeV

p_tot = p₁ + p₂ = 3 GeV/c

جرم ناوردای سیستم:

M²c^4 = E_tot² - (p_tot c)²

پس:

M²c^4 = 13² - 3² = 169 - 9 = 160

بنابراین:

Mc² = √160 = 4√10 ≈ 12.65 GeV

نتیجه:

جرم ناوردای این سیستم برابر است با:

M ≈ 12.65 GeV/c²

حالا سرعت چارچوب مرکز جرم از رابطهٔ

V = p_tot c² / E_tot

به دست می‌آید. در واحدهایی که c=1:

V/c = p_tot c / E_tot = 3/13 ≈ 0.231

یعنی اگر با سرعت حدود 0.231c در جهت مثبت x حرکت کنیم، به چارچوب سکون سیستم می‌رسیم.

مثال ۵) واپاشی یک ذرهٔ ساکن به دو فوتون

فرض کنید ذره‌ای با جرم سکون

M = 2 MeV/c²

در حال سکون واپاشی می‌کند به دو فوتون:

X → γ + γ

چون ذرهٔ اولیه ساکن است:

P_initial = (Mc, 0)

پس:

E_initial = Mc² = 2 MeV

p_initial = 0

برای اینکه تکانه پایسته بماند، دو فوتون باید در دو جهت مخالف حرکت کنند و چون فوتون‌ها:

E = pc

دارند، انرژی دو فوتون باید برابر باشد:

E₁ = E₂ = 1 MeV

و تکانه‌هایشان:

p₁ = +1 MeV/c

p₂ = -1 MeV/c

حالا جمع چهارتکانهٔ نهایی:

E_tot = 1 + 1 = 2 MeV

p_tot = 1 - 1 = 0

پس:

M²c^4 = E_tot² - (p_tot c)² = 2² - 0 = 4

یعنی:

Mc² = 2 MeV

نتیجه: جرم ناوردای سیستم نهایی همان جرم ذرهٔ اولیه است؛ پس ناوردایی چهارتکانه با پایستگی انرژی-تکانه کاملاً سازگار است.

مثال ۶) دو فوتونِ روبه‌رو می‌توانند یک جرم ناوردا بسازند

این مثال خیلی قشنگ است، چون هر فوتون به‌تنهایی جرم سکون صفر دارد، ولی مجموعِ دو فوتون می‌تواند جرم ناوردا داشته باشد.

فرض کنید دو فوتون روبه‌روی هم حرکت می‌کنند:

Photon 1: E₁ = 4 MeV, p₁ = +4 MeV/c

Photon 2: E₂ = 6 MeV, p₂ = -6 MeV/c

جمع انرژی و تکانه:

E_tot = 4 + 6 = 10 MeV

p_tot = 4/c - 6/c = -2 MeV/c

جرم ناوردای سیستم:

M²c^4 = E_tot² - (p_tot c)²

= 10² - 2²

= 100 - 4

= 96

پس:

Mc² = √96 = 4√6 ≈ 9.80 MeV

نتیجه:

با اینکه هر فوتون جداگانه

m = 0

دارد، ولی سیستم دوتایی آن‌ها جرم ناوردا دارد:

M ≈ 9.80 MeV/c²

نکتهٔ خیلی مهم

برای یک ذره:

(E/c)² - p² = m²c²

برای یک سیستم چندذره‌ای:

M²c² = P_tot·P_tot

یعنی اول چهارتکانه‌ها را جمع می‌کنیم، بعد ناوردا را می‌گیریم.

اگر بخواهی، مرحلهٔ بعدی می‌تواند یکی از این‌ها باشد:

1) مثال تستی چهارگزینه‌ای

2) مثال سخت‌تر از انرژی آستانه

3) واپاشی و برخورد در چارچوب مرکز جرم

4) تفاوت جرم سکون، جرم ناوردا، و جرم نسبیتی

۳ مثالِ دیگر از ناوردایی چهارتکانه (هم عددی، هم مفهومی)

مثال ۱) از روی انرژی و تکانه، جرم سکون را پیدا کن (ناورداییِ m)

فرض کن در یک چارچوب آزمایشگاه داریم:

E = 10 GeV

p = 8 GeV/c

ناوردای چهارتکانه می‌گوید:

(E/c)² − p² = m²c²

پس (با واحدهای طبیعی GeV/c²):

m²c^4 = E² − (pc)²

m = √(E² − (pc)²) / c²

اینجا:

pc = 8 GeV

m = √(10² − 8²) GeV/c²

= √(100 − 64)

= √36

= 6 GeV/c²

نتیجه: جرم سکونِ ذره ۶ GeV/c² است و این مقدار ناورداست (هر ناظری همین را به دست می‌آورد).

مثال ۲) برخورد فوتون + الکترونِ ساکن → جرم ناوردا/انرژی مرکز جرم

این مثال خیلی مهم است چون «جمع چهارتکانه‌ها» هم یک چهاربردار است و ناوردای مجموع معنی فیزیکی دارد (جرم ناوردای سیستم).

فرض کن یک فوتون با انرژی

E_γ = 1 MeV

به یک الکترونِ ساکن می‌خورد. جرم سکون الکترون:

m_e c² = 0.511 MeV

چهارتکانه‌ها:

- فوتون:

P_γ = (E_γ/c , p_γ)

و چون فوتون بی‌جرم است: p_γ = E_γ/c

- الکترونِ ساکن:

P_e = (m_e c , 0)

پس چهارتکانه کل:

P_tot = P_e + P_γ = ( (m_e c² + E_γ)/c , E_γ/c )

جرم ناوردای سیستم از:

M²c² = P_tot·P_tot = (E_tot/c)² − p_tot²

→ M²c^4 = E_tot² − (p_tot c)²

اینجا:

E_tot = m_e c² + E_γ = 0.511 + 1 = 1.511 MeV

p_tot c = E_γ = 1 MeV

پس:

M c² = √(E_tot² − (p_tot c)²)

= √(1.511² − 1²) MeV

= √(2.283 − 1) MeV

= √1.283 MeV

≈ 1.133 MeV

نتیجهٔ فیزیکی:

سیستم «فوتون+الکترون» یک جرم ناوردا دارد:

M c² ≈ 1.133 MeV

این مقدار همان «انرژی مرکز جرم» (در واحد انرژی) است و برای همهٔ چارچوب‌ها یکسان است.

مثال ۳) فوتون‌ها هم چهارتکانه دارند ولی ناوردای‌شان صفر است

فرض کن یک فوتون با انرژی:

E = 3 eV

حرکت می‌کند. چون فوتون:

p = E/c

ناوردای چهارتکانه:

(E/c)² − p² = (E/c)² − (E/c)² = 0

نتیجه: برای فوتون‌ها

P·P = 0

و این «صفر بودن» هم ناوردا است (هر ناظر هر انرژی‌ای ببیند، باز هم این ناوردا صفر می‌ماند).

اگر بگویی مثال‌ها را «کنکوری/تستی» می‌خواهی یا «سخت‌تر» (مثلاً واپاشی دوذره‌ای، یا محاسبهٔ انرژی آستانه در تولید ذرات)، همان سبک را ادامه می‌دهم.

اثبات ناوردایی «چهارتکانه» (یعنی چرا (E/c)² − p² ثابت می‌ماند)

هدف

می‌خواهیم نشان دهیم کمیت زیر (در ۱+۱ بعد برای سادگی) تحت تبدیل لورنتس ناورداست:

P·P = (E/c)² − p_x²

و برابر است با:

(E/c)² − p_x² = m²c²

۱) تعریف چهارتکانه

تعریف استاندارد:

P^μ = (E/c, p_x)

و از طرف دیگر:

P^μ = m U^μ

که U^μ چهارسرعت است.

۲) اثبات کوتاهِ مفهومی (با چهارسرعت)

می‌دانیم:

U·U = c²

پس:

P·P = (mU)·(mU) = m² (U·U) = m² c²

این نشان می‌دهد اندازهٔ چهاربردار چهارتکانه ناورداست، چون «ضرب داخلی مینکوفسکی» تحت تبدیل لورنتس ثابت می‌ماند.

اما اگر بخواهی «کاملاً جبری و قدم‌به‌قدم» مثل فاصله‌زمان ببینی، بخش بعدی دقیقاً همان است.

۳) اثبات جبریِ قدم‌به‌قدم (با تبدیل‌های انرژی و تکانه)

اگر چارچوب S' نسبت به S با سرعت v در راستای +x حرکت کند، تبدیل لورنتس برای انرژی-تکانه (۱+۱ بعد) چنین است:

E' = γ (E - v p_x)
p'_x = γ (p_x - v E/c²)

می‌خواهیم نشان دهیم:

(E'/c)² − p'_x² = (E/c)² − p_x²

طرف چپ:

(E'/c)² − p'_x²
= [γ(E - v p_x)/c]² − [γ(p_x - vE/c²)]²
= γ²{ (E - vp_x)²/c² − (p_x - vE/c²)² }

حالا دو مربع را باز می‌کنیم:

الف) جملهٔ اول:

(E - vp_x)²/c²
= (E² - 2Evp_x + v²p_x²)/c²
= E²/c² - 2vEp_x/c² + (v²/c²)p_x²

ب) جملهٔ دوم:

(p_x - vE/c²)²
= p_x² - 2p_x(vE/c²) + (v²E²/c^4)
= p_x² - 2vEp_x/c² + (v²/c^4)E²

حالا تفریق می‌کنیم (دقت کنید جمله‌های میانی دقیقاً حذف می‌شوند):

(E - vp_x)²/c² − (p_x - vE/c²)²
= [E²/c² + (v²/c²)p_x²] - [p_x² + (v²/c^4)E²]

جدا می‌کنیم:

= E²/c² (1 - v²/c²) − p_x² (1 - v²/c²)
= (1 - v²/c²)(E²/c² − p_x²)

پس:

(E'/c)² − p'_x²
= γ² (1 - v²/c²)(E²/c² − p_x²)

اما:

γ² = 1/(1 - v²/c²)

پس:

γ²(1 - v²/c²) = 1

نتیجه:

(E'/c)² − p'_x² = (E/c)² − p_x²

یعنی کمیت (E/c)² − p² تحت لورنتس ناورداست.

۴) نتیجهٔ فیزیکی: چرا برابر m²c² است؟

این ناوردا را بهترین جا می‌توان در چارچوبِ سکون ذره حساب کرد.

در چارچوب سکون:

p = 0
E = mc²

پس:

(E/c)² − p² = (mc²/c)² − 0 = m²c²

چون این کمیت ناورداست، در هر چارچوب دیگری هم همین مقدار می‌شود.

۵) چک عددی خیلی کوتاه (برای اطمینان)

فرض کنید ذره‌ای با m = 2 kg و v_particle = 0.6c داریم.

γ = 1/√(1-0.6²) = 1.25
E = γmc² = 1.25×2×c² = 2.5 c²
p = γmv = 1.25×2×0.6c = 1.5 c

پس:

(E/c)² − p² = (2.5c)² − (1.5c)²
= (6.25 - 2.25)c²
= 4c²
= (2²)c²
= m²c²

جمع‌بندی یک‌خطی

چهارتکانه یک چهاربردار است و «طول مینکوفسکی» آن ناورداست؛ بنابراین:

(E/c)² − p² = m²c² (Lorentz invariant)

مثال از «ناوردایی» با اثبات (فاصله‌زمان مینکوفسکی)

صورت مسئله

نشان دهید کمیت زیر تحت تبدیل لورنتس (بوست در راستای x) ناورداست:

s² = c²t² − x²

(برای سادگی، فقط ۱+۱ بعد: زمان و یک بعد مکانی)

تبدیل لورنتس

x' = γ (x - vt)
t' = γ (t - vx/c²)
γ = 1/sqrt(1 - v²/c²)

می‌خواهیم ثابت کنیم:

c²t'² − x'² = c²t² − x²

اثبات (قدم‌به‌قدم)

ابتدا طرف چپ را می‌نویسیم:

c²t'² − x'²
= c²[γ(t - vx/c²)]² − [γ(x - vt)]²

گاما را فاکتور می‌گیریم:

= γ²{ c²(t - vx/c²)² − (x - vt)² }

دو مربع را باز می‌کنیم:

1) جملهٔ زمانی:

c²(t - vx/c²)²
= c²(t² - 2t(vx/c²) + (v²x²/c^4))
= c²t² - 2vtx + (v²x²/c²)

2) جملهٔ مکانی:

(x - vt)²
= x² - 2vtx + v²t²

حالا تفریق می‌کنیم:

c²(t - vx/c²)² − (x - vt)²
= (c²t² - 2vtx + v²x²/c²) - (x² - 2vtx + v²t²)

ترم‌های -2vtx و +2vtx حذف می‌شوند:

= c²t² + (v²x²/c²) - x² - v²t²

جدا می‌کنیم:

= (c²t² - v²t²) - (x² - v²x²/c²)
= t²(c² - v²) - x²(1 - v²/c²)

حالا توجه کنید:

c² - v² = c²(1 - v²/c²)

پس:

= c²(1 - v²/c²)t² - (1 - v²/c²)x²
= (1 - v²/c²)(c²t² - x²)

برمی‌گردیم به عبارت اصلی که جلویش γ² داشت:

c²t'² − x'²
= γ² (1 - v²/c²)(c²t² - x²)

اما:

γ² = 1/(1 - v²/c²)

پس:

γ²(1 - v²/c²) = 1

نتیجه:

c²t'² − x'² = c²t² − x²

یعنی فاصله‌زمان مینکوفسکی ناوردا است.

تعبیر فیزیکی خیلی کوتاه

ممکن است t و x جداگانه با عوض شدن ناظر تغییر کنند، اما ترکیب خاص
c²t² − x²
مثل «طول» در هندسهٔ چهاربعدی ثابت می‌ماند؛ همان چیزی که هندسهٔ فضا-زمان را تعریف می‌کند.

یک مثال عددی سریع (چک کردن ناوردایی)

فرض کنید در چارچوب S:

t = 5 μs → ct = 1500 m
x = 1200 m

پس:

s² = (1500)² - (1200)² = 810000 (m²)

حالا به چارچوب S' با v=0.6c برویم:

γ = 1/√(1-0.6²) = 1.25

x' = γ(x - vt)
= 1.25(1200 - 0.6c×5μs)
= 1.25(1200 - 900)
= 375 m

t' = γ(t - vx/c²)
= 1.25(5μs - 0.6×1200/c)

چون:

1200/c = 4 μs

پس:

t' = 1.25(5μs - 2.4μs) = 3.25 μs
→ ct' = 975 m

اکنون:

s'² = (975)² - (375)²
= 950625 - 140625
= 810000 (m²)

نتیجه: s'² = s² و ناوردایی عددی هم تأیید شد.

اگر دوست داری، یک «اثبات ناوردایی چهارتکانه» هم می‌توانم بدهم (اینکه چرا (E/c)² - p² ثابت می‌ماند) که خیلی درک تقارن لورنتسی را کامل می‌کند.

چند مثال برای فرمول‌نامهٔ نسبیت خاص

مثال ۱) اتساع زمان

سفینه‌ای با سرعت v = 0.8c حرکت می‌کند. در چارچوب زمین، مدت سفر
Δt = 10 s
است. زمان ویژهٔ سفینه را پیدا کنید.

γ = 1 / sqrt(1 - 0.8²)
= 1 / sqrt(0.36)
= 1.667

پس:

Δτ = Δt / γ = 10 / 1.667 ≈ 6 s

نتیجه:

زمانی که سرنشین سفینه احساس می‌کند، حدود 6 s است.

مثال ۲) انقباض طول

میله‌ای در چارچوب خودش طول
L0 = 5 m
دارد. اگر با سرعت
v = 0.6c
از کنار ناظر عبور کند، طول دیده‌شده چقدر است؟

γ = 1 / sqrt(1 - 0.6²)
= 1 / sqrt(0.64)
= 1.25

پس:

L = L0 / γ = 5 / 1.25 = 4 m

نتیجه:

ناظر طول میله را 4 m می‌بیند.

مثال ۳) نسبیت هم‌زمانی

دو رویداد در چارچوب
S
هم‌زمان‌اند:

Δt = 0

و فاصلهٔ مکانی آن‌ها:

Δx = 300 m

چارچوب
S'
با سرعت
v = 0.6c
حرکت می‌کند. اختلاف زمان در
S'
را بیابید.

γ = 1.25

فرمول:

Δt' = -γ(vΔx/c²)

جایگذاری:

Δt' = -1.25 × (0.6c × 300 / c²)
= -1.25 × (180 / c)

با
c = 3×10^8 m/s:

Δt' = -1.25 × (180 / 3×10^8)
= -7.5×10^-7 s
= -0.75 μs

نتیجه:

در چارچوب
S'
این دو رویداد هم‌زمان نیستند.

مثال ۴) جمع سرعت‌ها

در چارچوب
S
ذره‌ای با سرعت
u = 0.9c
حرکت می‌کند. چارچوب
S'
با سرعت
v = 0.6c
در همان راستا حرکت می‌کند. سرعت ذره در
S'
چقدر است؟

u' = (u - v)/(1 - uv/c²)

جایگذاری:

u' = (0.9c - 0.6c)/(1 - 0.9×0.6)
= 0.3c / 0.46
≈ 0.652c

نتیجه:

سرعت جدید حدود 0.652c است.

مثال ۵) چهارتکانه

ذره‌ای با جرم سکون
m = 2 kg
و سرعت
v = 0.6c
حرکت می‌کند. انرژی و تکانهٔ نسبیتی آن را حساب کنید.

ابتدا:

γ = 1.25

انرژی:

E = γmc²
= 1.25 × 2 × (3×10^8)²
= 2.25×10^17 J

تکانه:

p = γmv
= 1.25 × 2 × 0.6 × 3×10^8
= 4.5×10^8 kg·m/s

پس چهارتکانه:

P^μ = (E/c , p, 0, 0)
= (7.5×10^8 , 4.5×10^8 , 0 , 0)

مثال ۶) ناوردایی انرژی-تکانه

برای همان ذرهٔ مثال قبل، بررسی کنید که:

(E/c)² - p² = m²c²

داریم:

(E/c)² = (7.5×10^8)² = 5.625×10^17
p² = (4.5×10^8)² = 2.025×10^17

پس:

(E/c)² - p² = 5.625×10^17 - 2.025×10^17
= 3.6×10^17

و:

m²c² = 2² × (3×10^8)² = 4 × 9×10^16 = 3.6×10^17

نتیجه:

رابطهٔ ناوردای انرژی-تکانه درست است.

جمع‌بندی

این مثال‌ها تقریباً مهم‌ترین بخش‌های فرمول‌نامه را پوشش می‌دهند:


- اتساع زمان

- انقباض طول

- نسبیت هم‌زمانی

- جمع سرعت‌ها

- انرژی و تکانهٔ نسبیتی

- ناوردایی چهارتکانه

اگر بخواهی، در پیام بعدی می‌توانم همین‌ها را به صورت «مثال‌های تستی کوتاه با جواب نهایی» هم برای وبلاگ آماده کنم.

فرمول‌نامهٔ خلاصهٔ نسبیت خاص (مناسب برای وبلاگ/بلاگفا)

۰) نمادگذاری سریع

c = speed of light β = v/c γ = 1/sqrt(1-β²)

۱) تبدیل‌های لورنتس (حرکت S' با سرعت v در راستای +x نسبت به S)

x' = γ (x - v t) t' = γ (t - v x / c²) y' = y z' = z

فرم معکوس:

x = γ (x' + v t') t = γ (t' + v x' / c²)

۲) فاصله‌زمان (ناوردای لورنتسی)
برای دو رویداد:

Δs² = c²Δt² − Δx² − Δy² − Δz²

شرط نوع جدایی:

timelike: c²Δt² > Δr² lightlike: c²Δt² = Δr² spacelike: c²Δt² < Δr²

که:

Δr² = Δx² + Δy² + Δz²

۳) زمان ویژه (Proper Time)

dτ = dt * sqrt(1 - v²/c²) = dt/γ

برای سرعت ثابت:

Δτ = Δt/γ

ارتباط با فاصله‌زمان (برای مسیرهای timelike):

c² dτ² = ds²

۴) اتساع زمان و انقباض طول
اتساع زمان (ساعت متحرک کندتر):

Δt = γ Δτ

انقباض طول (طول در راستای حرکت):

L = L0/γ

۵) نسبیت هم‌زمانی
اگر در S هم‌زمان باشند (Δt=0):

Δt' = - γ (v Δx / c²)

۶) جمع سرعت‌ها (یک‌بعدی)
اگر در S، جسم با سرعت u حرکت کند و S' با v:

u' = (u - v) / (1 - uv/c²)

(فرم معکوس)

u = (u' + v) / (1 + u'v/c²)

۷) راپیدیتی (Rapidity) — جمع‌پذیر و تمیز

β = tanh(η) γ = cosh(η) γβ = sinh(η)

مزیت:

η_total = η1 + η2

۸) سرعت ویژه (Proper Velocity)

w = dx/dτ = γ v

برعکس:

v = w / sqrt(1 + (w/c)²) γ = sqrt(1 + (w/c)²)

نکته: v همیشه کمتر از c است، ولی w می‌تواند از c بزرگ‌تر شود.

۹) چهاربردار مکان

X^μ = (ct, x, y, z)

متریک مینکوفسکی (قرارداد +---):

A·B = A^0 B^0 - A^1 B^1 - A^2 B^2 - A^3 B^3

۱۰) چهارسرعت (Four-velocity)

U^μ = dX^μ/dτ = (γc, γv_x, γv_y, γv_z)

ناوردا:

U·U = c²

۱۱) چهارتکانه (Four-momentum)

P^μ = m U^μ = (E/c, p_x, p_y, p_z)

روابط مهم:

E = γ m c² p = γ m v

ناوردای انرژی-تکانه:

P·P = (E/c)² - p² = m² c²

فرم رایج:

E² = (pc)² + (mc²)²

۱۲) ایدهٔ کلیدی: «ناوردایی» و «تقارن» در یک خط
- ناوردایی لورنتسی = ثابت ماندن کمیت تحت تبدیل‌های لورنتس.
- تقارن لورنتس ⇒ قوانین فیزیک در همهٔ چارچوب‌های لخت یکسان‌اند.
- (نکتهٔ مفهومی) تقارن‌های انتقالی/چرخشی ⇒ پایستگی‌ها (قضیه نودر).

جمع‌بندی خیلی کوتاه

γ = 1/sqrt(1-v²/c²) x' = γ(x-vt) t' = γ(t - vx/c²) s² = c²t² - x² - y² - z² (invariant) U·U = c² (E/c)² - p² = m²c²

ناوردایی از دیدگاه «تقارن» (چرا ناوردایی مهم است؟)

در فیزیک مدرن، یک اصل بسیار طلایی وجود دارد:

«هر ناوردایی، ریشه در یک تقارن دارد و هر تقارن، منجر به یک قانون بقا (پایستگی) می‌شود.» (قضیه نودر - Emmy Noether).

وقتی می‌گوییم فضا-زمان دارای تقارن است، یعنی اگر ما سیستم را در فضا جابه‌جا کنیم، بچرخانیم یا ناظرمان را عوض کنیم، قوانین بنیادی فیزیک نباید عوض شوند.

۱. تقارن لورنتس (Lorentz Symmetry)

این تقارن سنگ‌بنای نسبیت خاص است. یعنی فضا و زمان طوری به هم بافته شده‌اند که اگر شما با سرعت ثابت حرکت کنید (بوست لورنتس)، ساختار کلی فیزیک تغییر نمی‌کند.

ناوردایِ این تقارن، همان فاصله‌زمان ($s^2$) است.

- **شهود:** همان‌طور که در هندسه اقلیدسی، چرخاندن یک خط‌کش باعث تغییر طول آن نمی‌شود (ناوردایی تحت چرخش)، در نسبیت هم حرکت با سرعت بالا باعث تغییر "طول چهاربعدی" رویدادها نمی‌شود.

۲. انواع تقارن‌های فضا-زمانی و ناورداهای آن‌ها

در نسبیت خاص، ما با «گروه پوانکاره» (Poincaré group) سر و کار داریم که شامل تقارن‌های زیر است:

- **تقارن در انتقال مکانی (Translation in Space):**

اگر آزمایشی را امروز در تهران انجام دهید یا در مریخ، قوانین فیزیک یکی است.

نتیجه: پایستگی تکانه (Momentum).

- **تقارن در انتقال زمانی (Translation in Time):**

اگر آزمایش را الان انجام دهید یا ۲ ساعت بعد، قوانین فیزیک فرقی نمی‌کند.

نتیجه: پایستگی انرژی (Energy).

- **تقارن چرخشی (Rotation):**

اگر دستگاه آزمایش را ۹۰ درجه بچرخانید، فیزیک عوض نمی‌شود.

نتیجه: پایستگی تکانه زاویه‌ای (Angular Momentum).

- **تقارن بوست (Boost):**

اگر با سرعت ثابت حرکت کنید، قوانین فیزیک یکسان باقی می‌مانند.

نتیجه: ناوردایی لورنتسی.

۳. چرا ناوردایی برای فیزیک‌دان‌ها حیاتی است؟

۱. **کشف قوانین جهانی:** فیزیک‌دان‌ها به دنبال فرمول‌هایی هستند که "ناوردا" باشند. اگر فرمولی بنویسیم که فقط برای ناظر ساکن درست باشد، آن فرمول یک قانون بنیادی نیست. قانونی بنیادی است که در همه چارچوب‌ها (با تبدیل لورنتس) شکلش حفظ شود.

۲. **ساده‌سازی محاسبات:** وقتی می‌دانیم کمیتی مثل $P \cdot P = m^2c^2$ ناوردا است، می‌توانیم محاسبات سخت را در چارچوبی انجام دهیم که ذره در آن ساکن است (جایی که $p=0$) و بعد مطمئن باشیم که نتیجه در تمام چارچوب‌های دیگر هم معتبر است.

۴. هندسه مینکوفسکی: فراتر از اقلیدس

در هندسه معمولی (اقلیدسی)، فاصله ناوردا این است:

d² = x² + y² + z²

در این هندسه، زمان یک پارامتر بیرونی است.

اما در نسبیت، فضا و زمان با هم یک "موجود" به نام **متریک مینکوفسکی** می‌سازند که ناوردای آن این است:

s² = (ct)² − (x² + y² + z²)

علامت منفی پشتِ بخش مکانی، نشان‌دهنده تقارن خاص فضا-زمان است. این منفی است که اجازه می‌دهد زمان منبسط و طول منقبض شود تا $s^2$ ثابت بماند.

خلاصه کلام:

تقارن یعنی «تغییر بدون تغییر». یعنی شما چیزی را عوض می‌کنید (سرعت ناظر)، اما یک چیز دیگر (قوانین فیزیک و ناورداها) ثابت می‌ماند. بدون این تقارن‌ها، جهان غیرقابل پیش‌بینی و آشفته بود.

اگر موافق باشی، می‌توانیم این بحث را با یک «اینفوگرافیک متنی» یا یک «خلاصه فرمول‌نامه برای وبلاگ» جمع‌بندی کنیم که تمام این مفاهیم (ناوردایی، چهاربردارها و تقارن) را در یک جا داشته باشی. چطور است؟

تفاوت «ناوردایی» و «پایا بودن»

این دو واژه خیلی به هم نزدیک‌اند، ولی دقیقاً یکی نیستند. در فیزیک اگر این تفاوت را خوب بفهمیم، خیلی از ابهام‌ها برطرف می‌شود.

۱) ناوردایی چیست؟

ناوردایی یعنی یک کمیت یا یک رابطه، وقتی چارچوب مرجع یا مختصات را عوض می‌کنیم، تغییر نکند.

یعنی:

Q' = Q

مثلاً در نسبیت خاص:

s² = c²Δt² − Δx² − Δy² − Δz²

یک کمیت ناوردا است، چون تحت تبدیل لورنتس ثابت می‌ماند.

پس ناوردایی بیشتر به این سؤال مربوط است:

«اگر ناظر را عوض کنیم، چه چیزی ثابت می‌ماند؟»

۲) پایا بودن چیست؟

پایا بودن معمولاً یعنی یک کمیت در طول زمان تغییر نکند.

مثلاً اگر انرژی یک سیستم در طول حرکت ثابت بماند، می‌گوییم انرژی پایا یا باقی‌مانده است.

یعنی:

dQ/dt = 0

پس پایا بودن بیشتر به این سؤال مربوط است:

«آیا این کمیت با گذشت زمان تغییر می‌کند یا نه؟»

فرق اصلی در یک جمله

- ناوردایی → مربوط به تغییر ناظر یا دستگاه مختصات است.

- پایا بودن → مربوط به تحول زمانی یک سیستم است.

مثال برای روشن شدن تفاوت

مثال ۱) فاصله‌زمان

s² = c²Δt² − Δx² − Δy² − Δz²

این کمیت ناوردا است، چون همهٔ ناظرهای لخت مقدار یکسانی برایش به دست می‌آورند.

اما اینکه بگوییم «پایا» است یا نه، اصلاً بحث دیگری است؛ چون اینجا داریم یک فاصله بین دو رویداد را حساب می‌کنیم، نه تحول زمانی یک کمیت در یک سیستم.

مثال ۲) انرژی یک ذرهٔ آزاد

انرژی ممکن است در یک چارچوب خاص در طول زمان ثابت باشد، پس پایا باشد.

اما خود انرژی معمولاً ناوردا نیست، چون ناظرهای مختلف انرژی‌های متفاوتی برای همان ذره اندازه می‌گیرند.

یعنی:

- انرژی می‌تواند پایا باشد،

- ولی ناوردا نباشد.

مثال ۳) جرم سکون

جرم سکون ناوردا است، چون همهٔ ناظرها همان مقدار را به دست می‌آورند.

در بسیاری از مسائل هم برای یک ذرهٔ معین پایا هم هست، چون با زمان عوض نمی‌شود.

پس یک کمیت می‌تواند هم‌زمان هم ناوردا باشد و هم پایا.

جدول مقایسه

ناورداهای لورنتسی یعنی چه؟

ناورداهای لورنتسی کمیت‌هایی هستند که تحت تبدیل لورنتس تغییر نمی‌کنند.

یعنی اگر از چارچوب S به چارچوب S' برویم، مقدار آن‌ها ثابت می‌ماند.

مهم‌ترین نمونه‌ها:

۱) فاصله‌زمان

s² = c²Δt² − Δx² − Δy² − Δz²

۲) اندازهٔ چهارسرعت

U·U = c²

۳) اندازهٔ چهارتکانه

P·P = (E/c)² − p² = m²c²

۴) زمان ویژه

dτ² = dt² − (dx²+dy²+dz²)/c²

یا به شکل فشرده‌تر:

c²dτ² = ds²

یک نکتهٔ مهم

هر کمیت پایا، لزوماً ناوردا نیست؛

و هر کمیت ناوردا هم لزوماً موضوع «پایا بودن» نیست.

پس این دو مفهوم را نباید با هم یکی گرفت.

خلاصهٔ خیلی کوتاه

- ناوردایی = ثابت ماندن نسبت به تغییر ناظر

- پایا بودن = ثابت ماندن در طول زمان

نتیجه

در نسبیت خاص، وقتی می‌گوییم یک کمیت «ناوردا» است، منظورمان این است که همهٔ ناظرهای لخت بر مقدار آن توافق دارند. اما وقتی می‌گوییم یک کمیت «پایا» است، منظورمان این است که آن کمیت در فرایند زمانی تغییر نمی‌کند.

اگر بخواهی، در پیام بعدی می‌توانم «ناوردایی را از دید تقارن» هم توضیح بدهم؛ یعنی اینکه چرا ناوردایی‌ها به گروه لورنتس و ساختار فضا-زمان ربط دارند.

ناوردایی یعنی چه؟

در فیزیک، ناوردایی یعنی یک کمیت یا رابطه، با عوض شدن دستگاه مختصات یا چارچوب مرجع، تغییر نکند.

یعنی اگر دو ناظر مختلف یک پدیده را توصیف کنند، ممکن است بعضی کمیت‌ها برایشان فرق کند، اما بعضی چیزها برای هر دو یکسان بماند. به این ویژگی می‌گوییم ناوردایی.

مثال ساده در نسبیت خاص

مختصات زمان و مکانِ یک رویداد در دو چارچوب مختلف عوض می‌شود:

t, x, y, z → t', x', y', z'

اما کمیت

s² = c²Δt² − Δx² − Δy² − Δz²

برای همهٔ ناظرهای لخت یکسان می‌ماند.

یعنی:

s² = s'²

پس s² یک کمیت ناوردا است.

معنی شهودی

مثلاً ممکن است:

- یک ناظر بگوید فاصلهٔ مکانی بزرگ‌تر است،

- ناظر دیگر بگوید فاصلهٔ زمانی بزرگ‌تر است،

اما ترکیب خاص آن‌ها در فرمول فاصله‌زمان ثابت می‌ماند.

یعنی خودِ زمان و مکان جداگانه ناوردا نیستند، ولی فاصله‌زمان ناوردا است.

مثال‌های مهم از ناوردایی در نسبیت

۱) ناوردایی فاصله‌زمان

s² = c²Δt² − Δx² − Δy² − Δz²

۲) ناوردایی اندازهٔ چهارسرعت

U·U = c²

۳) ناوردایی اندازهٔ چهارتکانه

P·P = m²c²

تعریف خیلی خلاصه

ناوردایی یعنی:

«کمیتی که با تبدیل لورنتس تغییر نمی‌کند.»

یک بیان ریاضی

اگر کمیتی در چارچوب S و S' مقدار یکسان داشته باشد، می‌گوییم ناوردا است:

Q' = Q

تفاوت با کمیت‌های معمولی

مثلاً این‌ها معمولاً ناوردا نیستند:

- زمان t

- مکان x

- سرعت معمولی v

اما این‌ها ناوردا هستند:

- فاصله‌زمان

- جرم سکون

- اندازهٔ چهاربردارها

نتیجه

در نسبیت خاص، ناوردایی یکی از مهم‌ترین ایده‌هاست، چون نشان می‌دهد با اینکه ناظرها ممکن است زمان و مکان را متفاوت ببینند، بعضی کمیت‌های بنیادی برای همه یکسان‌اند.

اگر بخواهی، در پیام بعدی می‌توانم تفاوت «ناوردایی» و «پایا بودن» یا «ناورداهای لورنتسی» را هم خیلی دقیق توضیح بدهم.

نسخهٔ تمیزتر و یکدست‌ترِ مثال‌های ۴بُعدی سخت‌تر

مثال ۴) ناوردا بودن فاصله‌زمان

فرض کنید چارچوب S' با سرعت

v = 0.6c

در راستای محور x

نسبت به چارچوب S حرکت می‌کند.

پس داریم:

β = 0.6

γ = 1 / √(1 − 0.6²) = 1.25

دو رویداد در چارچوب S:

A : (t=0, x=0, y=0, z=0)

B : (t=5 μs, x=1200 m, y=0, z=0)

گام ۱) محاسبهٔ فاصله‌زمان در S

Δt = 5 μs

Δx = 1200 m

Δy = 0

Δz = 0

ابتدا:

cΔt = (3×10^8)(5×10^-6) = 1500 m

پس:

s² = c²Δt² − Δx² − Δy² − Δz²

s² = (1500)² − (1200)² = 2,250,000 − 1,440,000 = 810,000 m²

بنابراین:

s = 900 m

و زمانِ ویژه:

Δτ = s/c = 900 / (3×10^8) = 3×10^-6 s = 3 μs

گام ۲) تبدیل رویداد B به چارچوب S'

x' = γ(x − vt)

t' = γ(t − vx/c²)

y' = y

z' = z

ابتدا:

vt = (0.6c)(5 μs) = 0.6 × 1500 = 900 m

پس:

x'_B = 1.25(1200 − 900) = 1.25 × 300 = 375 m

اکنون:

x/c = 1200 / (3×10^8) s = 4 μs

در نتیجه:

vx/c² = 0.6 × 4 μs = 2.4 μs

پس:

t'_B = 1.25(5 − 2.4) μs = 1.25 × 2.6 μs = 3.25 μs

گام ۳) فاصله‌زمان در S'

Δt' = 3.25 μs

Δx' = 375 m

Δy' = 0

Δz' = 0

ابتدا:

cΔt' = (3×10^8)(3.25×10^-6) = 975 m

پس:

s'² = c²Δt'² − Δx'²

s'² = (975)² − (375)² = 950625 − 140625 = 810000 m²

یعنی:

s'² = s²

نتیجه: فاصله‌زمان در همهٔ چارچوب‌های لخت ناوردا است.

مثال ۵) چهارسرعت

ذره‌ای در چارچوب S با سرعت

v = 0.8c

در راستای x حرکت می‌کند.

چهارسرعت آن را بیابید.

تعریف چهارسرعت:

U^μ = dX^μ / dτ = (c dt/dτ, dx/dτ, dy/dτ, dz/dτ)

و چون:

dt/dτ = γ

داریم:

U^μ = (γc, γv, 0, 0)

برای

v = 0.8c:

γ = 1 / √(1 − 0.8²) = 1 / √0.36 = 1.6667

پس:

U^μ = (1.6667c, 1.6667×0.8c, 0, 0)

U^μ = (1.6667c, 1.3333c, 0, 0)

بررسی ناوردای چهارسرعت:

U·U = (γc)² − (γv)²

= γ²(c² − v²)

= γ²c²(1 − v²/c²)

= c²

نتیجه: اندازهٔ مینکوفسکی چهارسرعت همیشه برابر c² است.

مثال ۶) چهارتکانه و انرژی نسبیتی

جرم سکون یک ذره

m = 2 kg

است و این ذره با سرعت

v = 0.6c

در راستای محور x حرکت می‌کند.

چهارتکانهٔ آن را بیابید.

تعریف چهارتکانه:

P^μ = mU^μ = (E/c, p_x, p_y, p_z)

پس:

E = γmc²

p_x = γmv

ابتدا:

γ = 1 / √(1 − 0.6²) = 1.25

انرژی:

E = 1.25 × 2 × c² = 2.5c²

اگر

c = 3×10^8 m/s

قرار دهیم:

E = 2.5 × 9×10^16

= 2.25×10^17 J

تکانه:

p_x = 1.25 × 2 × 0.6c

= 1.5c

عددی:

p_x = 1.5 × 3×10^8

= 4.5×10^8 kg.m/s

پس:

P^μ = (E/c, p_x, 0, 0)

یعنی:

E/c = (2.25×10^17) / (3×10^8) = 7.5×10^8

بنابراین:

P^μ = (7.5×10^8, 4.5×10^8, 0, 0)

بررسی ناوردای چهارتکانه:

P·P = (E/c)² − p²

از طرف دیگر باید داشته باشیم:

P·P = m²c²

برای این مثال:

m²c² = (2²)(3×10^8)² = 4 × 9×10^16 = 3.6×10^17

و اگر مستقیم هم حساب کنیم:

(E/c)² − p²

= (7.5×10^8)² − (4.5×10^8)²

= 56.25×10^16 − 20.25×10^16

= 36×10^16

= 3.6×10^17

پس رابطه درست است.

جمع‌بندی

در نسبیت خاص:

- s² = c²t² − x² − y² − z² یک کمیت ناوردا است.

- چهارسرعت برابر U^μ = (γc, γv_x, γv_y, γv_z) است.

- چهارتکانه برابر P^μ = (E/c, p_x, p_y, p_z) است.

- ناوردای چهارتکانه: (E/c)² − p² = m²c²

اگر بخواهی، در پیام بعدی می‌توانم همین‌ها را به صورت تستی یا به صورت فرمول‌نامهٔ جمع‌وجور برای بلاگفا هم بنویسم.