משוואת לנז'בן

בפיזיקה, משוואת לנז'בן, היא משוואה דיפרנציאלית סטוכסטית, המתארת התפתחות לפי זמן של מערכת בעלת דרגות חופש. דרגות חופש אלו, מאופיינות בדרך כלל על ידי משתנים בעלי תכונות קולקטיביות הנעים בפרקי זמן ארוכים מאוד ביחס למשתנים אחרים במערכת, המתנהגים באופן מהיר מאוד או לא סדור.

לעיתים קרובות, ניתן לתאר באמצעות משוואה זו מערכות מקרוסקופיות בעלות חלק מיקרוסקופי, המשתנה בפרקי זמן קצרים מאוד ביחס למערכת הנמדדת. משוואה זו נוסחה לראשונה על ידי פול לנז'בן ב1908, בניסיונו לתאר תנועה בראונית, הכוללת תנועה אקראית של חלקיק בנוזל עקב תנודות תרמיות סביב שיווי משקל. כיום, נעשה שימוש נרחב במשוואת לנז'בן, עבור מערכות בעלות משתנים סטוכסטיים בתחומי הפיזיקה הסטטיסטית, ביופיזיקה וכימיה בפרט.

היסטוריה

תנועה בראונית

התיעוד הראשון של תנועה בראונית היתה על ידי רוברט בראון (בוטניקאי) ב-1827. בראון ראה, כי אבקנים שנעו על פני מים, תחת מיקרוסקופ, מבצעים תנועה לא יציבה ורנדומלית. הראשון לתת הסבר מתמטי ופיזיקלי לכך היה אלברט איינשטיין, אשר הסביר זאת במאמרו ב-1905. באותו מאמר, תוך הנחת שיווי משקל בנוזל, שילב איינשטיין בין התפלגות מקסוול בולצמן לתהליך של הילוך מקרי.

איינשטיין טען, כי כאשר חלקיק נמצא בתוך נוזל, ולא חל עליו כוח גרר, אז בהינתן התנגשות עם מולקולה אחרת, מהירותו תשתנה. לעומת זאת, אם תווך הנוזל הוא בעל רמה מספיקה של צמיגות, השינוי במהירות ידעך מהר מאוד, כך שסך התוצאה של ההתנגשות יצור שינוי בהעתק החלקיק בלבד. למעשה, איינשטיין הסיק שתרומה של כל ההתנגשויות תביא ל"קפיצות" רנדומליות במיקום החלקיק הבראוני- כך שהחלקיק יבצע הילוך מקרי. בלקיחת "קפיצות" אינפיטסימליות, הוא הגיע למשוואה דיפרנציאלית חלקית במימד אחד- שהיא משוואת דיפוזיה. איינשטיין הצליח למצוא את הקשר בין קבוע הדיפוזיה המקרוסקופי, לתכונותיו האטומיות של החומר:

הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle D={k_BT \over 6\pi \eta r}}

כאשר D הינו קבוע הדיפוזיה, הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle k_B} הינו קבוע בולצמן, T מייצגת את הטמפרטורה, הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \eta }  הינו מקדם הצמיגות של הנוזל, ו- r הוא הרדיוס של החלקיק הבראוני. אופי המשוואה הדיפרנציאלית החלקית אליה הגיע איינשטיין (כמו גם מריון סמולוצ'ובסקי, שעל שמו נקראת המשוואה), הינה מקרה פרטי של משוואה מסוג המוכר כמשוואות פוקר-פלאנק - המראות התפתחות בזמן של צפיפות הסתברות של מערכת. בפתרונו, הצליח איינשטיין להראות, שריבוע הערך הממוצע של העתק החלקיק הוא פרופורציונלי לקבוע של הדיפוזיה. בכך הגיע איינשטיין לתוצאה מהפכנית, והראה (תוך שימוש בהתפלגות מקסוול-בולצמן בשיווי משקל), כי המשוואה עבור צפיפות ההסתברות תתפתח באופן ליניארי בזמן. פרט זה חשוב במיוחד, מכיוון שכך התפתחות זו תלויה בטמפרטורה וצמיגות הנוזל בלבד.

הצגת משוואת לנז'בן

התיאוריה שפיתחו סמולוצ'ובסקי ואיינשטיין עבור התנועה הבראונית, נשענה כמעט לחלוטין על התפלגות צפיפות ההסתברות, ועל משוואות פוקר-פלאנק. ב-1908, הציג פול לנז'בן במאמרו "על התיאוריה של תנועה בראונית", דרך שונה לפתרון בעיה זו. לנז'בן, ניסח את הבעיה באמצעות משוואת כוחות ניוטונית, המכילה חלק "דטרמיניסטי" (המכיל את כוח הגרר של חלקיק הנמצא בנוזל, ופוטנציאלים נוספים וידועים) וחלק אקראי- אשר מייצג את ההתנגשויות של חלקיקי הנוזל בחלקיק הבראוני. אותו חלק רנדומלי, בדיעבד אינו ניתן למדידה, עקב פרקי הזמן הקצרים ביותר בהם נעשות ההתנגשויות השונות. גישתו של לנז'בן הייתה מהפכנית בזמנו, וקידמה משמעותית את נושא המחקר של משוואות דיפרנציאליות סטוכסטיות, הן בפיזיקה והן במתמטיקה. בשנות ה-20 של המאה ה-20, נתן נורברט וינר את הניסוח המתמטי והגדיר תנועה בראונית, כמקרה פרטי של מערכת סטוכסטית. כמו כן הוכיח שתנועה בראונית מוגדרת על ידי פונקציה רציפה, שאינה גזירה באף נקודה. ב-1942, הגדיר קיושו איטו את האינטגרל הסטוכסטי- המבוסס גם הוא על תנועה בראונית.

משוואת לנז'בן אינרטית במימד אחד

משוואת לנז'בן במימד אחד, על פי החוק השני של ניוטון עבור תנועה של חלקיק בראוני בנוזל, הינה:

הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle m{dx^2 \over dt^2} = -\xi {dx\over dt} -{dU \over dx} +R(t)}

כאשר החלק הדטרמיניסטי מכיל את האיברים:

• הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle -{dU \over dx}} – הכוח המאפיין פוטנציאלים שונים (כוח אלקטרומגנטי, כוח גרביטציוני וכו'). חלק זה אינו מופיע בפיתוח המשוואה עבור חלקיק בראוני בנוזל של לנז'בן, עקב תרומתן הזניחה של השפעות אלו.
• הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle -\xi {dx\over dt}} – כוח הגרר בנוזל, בדרך כלל מאופיין על ידי חוק סטוקס, הקובע כי כוח הגרר מתכונתי למהירות החלקיק (ראה חוק סטוקס). בצורה זו הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \xi \equiv 6\pi \eta R}  , כאשר הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \eta}  הוא מקדם הצמיגות של התווך הנוזלי, ו- R הינו רדיוס הכדורי של החלקיק.

החלק האקראי של המשוואה הינו R(t), המייצג את ההתנגשות האקראית של חלקיקי הנוזל בחלקיק הבראוני. לחלק זה המאפיינים נובעים מהשיקולים הפיזיקליים הבאים:

• R(t) אינו תלוי במפורש ב-הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle dx(t) \over dt} .
• R(t) משתנה בפרקי זמן קצרים בהרבה מ הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle dx(t) \over dt} .
• מכיוון שהתנגשות חלקיקי הנוזל בחלקיק הבראוני הינה אקראית ואיזוטרופית למרחב, ערך הממוצע של פונקציה זו על אנסמבל של חלקיקים שווה לאפס. הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle <R(t)>=0}
• מכיוון ש- R(t) משתנה בפרקי זמן קצרים בהרבה בהשוואה למיקום, ניתן להניח שכל התנגשות הינה מיידית. שינוי מהיר זה ניתן לבטא כ: הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle <R_i (t)R_j(t')>=2 \xi \delta_{ij} \delta(t-t')k_BT }

כאשר הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \delta(t-t') } הינה הדלתא של דיראק, ו- הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \delta_{ij} }  היא הדלתא של קרוניקר. שיקול זה נובע כמובן מהקירוב עבור הבדלי סדרי הגודל של פרק הזמן של פונקציה R(t) ל- הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle dx(t) \over dt}  ואף להעתק, והופך למדויק רק עבור טווחי זמן גדולים- המתארים תנועה עבור חלקיקים בראונים "מקרוסקופיים" (כאשר מסת החלקיק גדולה משמעותית ממסת חלקיקי הנוזל). יש לציין כי הסתכלות דיפרנציאלית על קשר הקורלציה אינו נכון לחלוטין, שכן גם מהירות החלקיק אינה מוגדרת היטב בגבול זה. בטיפול במשוואות דיפרנציאליות סטוכסטיות כגון משוואת לנז'בן, נהוג לקחת אינטגרל כך ש-

הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle m{dx \over dt} = \int_{}^{t} -{\xi {dx\over dt} -{dU \over dx} +R(t)}}

• ההסתברות למציאת ערך R(t) כלשהו, מתפלגת התפלגות גאוסיאנית, כך ש-

הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle P(R_1,R_2,...,R_n)= {1 \over \sqrt {2\pi <R^2>}} \exp {-\sum_{i}^nR_i^2 \over 2<R^2>}}

כאשר בגבול הרצף

הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \sim {1 \over \sqrt {2\pi <R^2>}} \exp {-\int_{t_1}^{t_n} R(t)^2 dt \over 2\pi <R^2> \Delta t}}

כאשר הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Delta t= t_{j+1} -t_j} .

החלק האקראי של המשוואה נקרא לעיתים כוח לנז'בן, ועבור חלקיק בראוני, חלק זה מראה כיצד שינויים בטמפרטורה מכתיבים את אופייה הסטוכסטי של התנועה הבראונית.

לדוגמה, נסתכל פתרון של משוואת לנז'בן עבור חלקיק בראוני. אם נכפיל את שני צידי המשוואה ב x(t):

הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle m{dx^2 \over dt^2}x(t) = -\xi {dx\over dt}x(t) +R(t)x(t)}

כאשר :

הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle x(t){d^2x(t) \over dt^2} = {1 \over 2} {d \over dt}({dx(t)^2 \over dt}) - {dx(t)^2 \over dt}}
 ו- הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle x(t){dx(t) \over dt} = {1 \over 2} {d \over dt}({dx(t)^2 \over dt})}
 

בהצבה במשוואה, ולקיחת ערך ממוצע עבור אנסמבל חלקיקים, נקבל:

הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle {1 \over 2} m{d \over dt}<{dx(t)^2 \over dt}> -m<{dx(t)^2 \over dt}> = -{\xi \over 2}<{dx(t)^2 \over dt}> +<R(t)x(t)>}

מכיוון שעבור החלק האקראי ערך ממוצע של אנסמבל חלקיקים מתאפס, חלק זה של המשוואה נופל. מחוק החלוקה השווה אנו יודעים הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle m<{dx(t) \over dt}> = k_BT} .

נסמן הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle <{dx(t)^2 \over dt}>=z} ,

ונקבל: הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle {m \over 2}<{dz \over dt}> + {\xi \over 2}z = k_BT}

עבור פרקי זמן גדולים מ הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle m \over \xi} , הפתרון הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle <{dx(t)^2 \over dt}>} מתפתח ליניארית בזמן. כלומר:

הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle <x(t)^2>-<x(0)^2> \approx {2k_BT \over \xi} t \longrightarrow 2Dt}

ובכך, הצטמצמנו למקרה הפרטי של פתרון איינשטיין לתנועה בראונית, ושהעתק החלקיק מתכונתי לשורש הטמפרטורה.

אפיון כללי של גישת לנז'בן

ניתן לתת למשוואות לנז'בן פורמליזם כללי יותר. בהנחה וישנה מערכת שהתנהגותה המקרוסקופית ידועה, וניתן לראות כי חייב להתקיים בה תנודות מסוימות. בשביל תיאור מלא של התנודות הללו יש לבצע את שלושת הצעדים הבאים:

צעד ראשון: יש לרשום את משוואות התנועה המקרוסקופיות של המערכת.

צעד שני: יש להוסיף כוח לנז'בן, המקיים את התכונות הבאות (חלקן אופיינו מוקדם יותר בפיתוח המקרה של תנועה בראונית):

1. במשוואת הכוחות ישנו גורם דטרמיניסטי, ליניארי עם מקדם ריסון, וגורם אקראי סטוכסטי (הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle R(t)} ) . הגורם הסטוכסטי הינו אי-רגולרי, אך תכונותיו הממוצעות כאנסמבל חלקיקים פשוטות, בצורה שאפשר להניח: שהמרחקים בין החלקיקים המרכיבים אותו הינו גדול כך שאין אינטראקציה ביניהם או שהמדידות שנלקחות הינם בפרקי זמן גדולים מספיק ביניהן.
2. התכונות הסטוכסטיות של הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle R(t)} ניתנות ללא קשר למהירות החלקיק, כך שהוא מתפקד ככוח חיצוני. בצורה זו, הממוצע שלו מתאפס הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle <{\displaystyle R(t)}> = 0}  
3. הכוח הסטוכסטי מתבטא עקב שינויים מהירים מאוד בפרקי הזמן הנתונים, ולכן ניתן להתייחס לפונקציית הקורלציה שלו: הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle <R (t)R(t')> \Gamma \delta (t-t')} . בתצורה זו הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Gamma} הינו קבוע. באופן מעשי, פונקציית הדלתא צריכה להיות פונקציה של הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle |t-t'|}  עם רוחב המתאים לזמן של התנגשות/ שינוי יחיד. ניתן להתייחס לפונקציה זו כפונקציה זו כל פעולה זו קצרה משמעותית בכמה סדרי גודל מכל שאר סדרי הגודל הזמניים בבעיה.
4. הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle R(t)}  הינה פונקציה עם התפלגות גאוסיאנית, כך שהמומנטים החיוביים שלה מקיימים: הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle <R (t_1)R(t_2)><R (t_3)R(t_4)>+... = \Gamma ^2 [\delta (t_1 - t_2)\delta (t_3 - t_4)+\delta (t_1 - t_3)\delta (t_2 - t_4) + \delta (t_1 - t_4)\delta (t_2 - t_3)]} באופן זה, R(t) מאופיינת כרעש לבן גאוסיאני.

צעד שלישי: יש לדרוש שהקבוע מתאים כך שהפתרון עבור תהליך סטוכסטי סטציונרי, ונותן בחזרה את ערכי ממוצע של התנודות כפי שהן ידועות ממכניקה סטטיסטית (או משיקולים פיזיקליים נוספים).

שלושת צעדים אלו מרכיבים את גישת לנז'בין, בה משתמשים בתחומים רבים בפיזיקה, כימיה וביולוגיה. התנודות אינן בהכרח צריכות לנבוע מתנודות תרמיות או מטבעם הדיסקרטי של חלקיקים. התנודות הסטוכסטיות יכולות לנבוע ממידול אות אקראי הנכנס לקו תמסורת, גדילה של מין ביולוגי תחת השפעת מזג אוויר, עומס אקראי על גשרים וכו'.

דוגמאות

משוואת לנז'בן עבור פוטנציאל הרמוני

עבור חלקיק בעל מסה m בממד אחד, הנע תחת פונטציאל הרמוני (הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle U(x) = {kx^2 \over 2}} ) משוואות לנז'בן הינם:

הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle {dx(t) \over dt}= v}

הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle {dv \over dt}= -{\xi \over m}v -\omega_0^2x +{1 \over m}R(t) }

כאשר הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \omega_0^2= {k \over m}} . הכוח האקראי הינו תהליך אקראי גאוסיאני, עבורו:

הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle <{\displaystyle R(t)}> = 0}

הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle <R (t_1)R(t_2)>=2 \xi k_BT \delta(t_1-t_2)}

בהסתכלות על צפיפות ספקטרום התדרים בהתמרת פורייה, נקבל:

הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle S_R(\omega)=\textstyle \int\limits_{-\infty}^{\infty} dt e^{i\omega t} <R(t)R(0)> = 2 \xi k_BT }

המאופיין כרעש לבן. אם נבצע את התמרת פורייה על משוואות לנז'בין ולא רק על קשר הקורולציה שלנו:

הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle -{i\omega x(\omega)} = v(\omega)}

הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle -{i\omega v(\omega)}= -{\xi \over m}v(\omega) -\omega_0^2x(\omega) +{1 \over m}R(\omega)}

מכן ניתן לפתור עבור הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle x(\omega) }  :

הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle x(\omega) = {1 \over m}{R(\omega) \over \omega _0 ^2 - \omega ^2 -{\xi \over m}i\omega} }

צפיפות ספקטרום התדרים הינה פרופורציונלית ל- הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle |x(\omega)^2| } ולכן :

הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle S_x(\omega)= {1 \over m^2}{S_R(\omega) \over {|\omega_0 ^2 -\omega^2 -{\xi \over m}i\omega}|^2} = {2\xi k_BT \over m^2 [(\omega _0^2 -\omega ^2)^2 +({\xi \over m} \omega)^2 ]} }

בהתמרה חזרה למרחב הזמן נקבל:

הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle C_x(t) = {1 \over 2\pi}\int\limits_{-\infty}^{\infty} d\omega e^{-i\omega t}S_x(\omega) = {2\xi k_BT \over 2\pi m^2} \int\limits_{-\infty}^{\infty} d\omega e^{-i\omega t} {1 \over [(\omega _0^2 -\omega ^2)^2 +({\xi \over m} \omega)^2 ]} }

אינטגרל זה ניתן לפתור במישור המרוכב, תוך שימוש במשפט השאריות. תוצאתו הסופית תהיה:

הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle C_x(t) = {k_BT \over m\omega_0^2}e^{-{\xi \over 2m}t} [cos (\sqrt{\omega_0^2 -{\xi^2 \over 4m^2} }t) + {\xi \over 2m\sqrt{\omega_0^2 -{\xi^2 \over 4m^2}}}sin(\sqrt{\omega_0^2 -{\xi^2 \over 4m^2}}t)] }

נשים לב שעבור t=0, נקבל חזרה את חוק החלוקה השווה - הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle {1 \over 2}m\omega_0^2 <x_0 ^2> = {1 \over 2} k_BT } .

תנודות במעגל חשמלי RC

אם נסתכל על מעגל חשמלי RC עם התנגדות ליניארית R. המטען Q על קבל (באופן מאקרוסקופי) חווה ריסון עקב ההתנגדות, אליו נוסיף חלק אקראי של רעש L(t).

הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle {dQ \over dt} = -{Q \over RC} +L(t) }

בשיווי משקל תרמי, האנרגיה האלקטרוסטטית הפוטנציאלית ידועה, כמו כן מחוק החלוקה השווה נקבל:

הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle {<Q^2> \over 2C} = {k_BT \over 2} }

הטמפרטורה T הינה של המעגל כולו, אך ניתן להניח שזוהי הטמפרטורה של הנגד בלבד, שכן רק שם מיוצר הרעש התרמי. אם כן, אפשר להסתכל על הנגד כמקור לתנודות בזרם הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Delta I } המתווספות לזרם המקרוסקופי הכולל. ניתן לראות כי הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Delta I = -L(t) }  ושקבוע הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Gamma } צריך להיות שווה ל- הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle {2k_BT \over R} } .ולכן נוכל לראות את הקשר:

הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle <\Delta I(t)\Delta I(t')> = {2k_BT \over R}\delta (t-t') }

באופן דומה ניתן לבצע עבור המתח כך ש הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Delta V = -RL(t) } משם ניתן לראות:

הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle <\Delta V(t)\Delta V(t')> = {2k_BTR}\delta (t-t') }

ואנחנו מקבלים באופן מיידי רעש ג'ונסון-נייקוויסט.

הקשר לאינטגרלי איטו וסטרטונוביץ'

כפי שהוזכר, קשר הקורלציה עבור R(t) הינו פרופורציונלי לפונקציית דלתא של דיראק, אך אינו בהכרח פונקציית דלתא, אלא פונקציה בעלת עליה חדה עם רוחב המתאים לפרק זמן קצר ביותר (הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \tau_c \rightarrow0 } ) . באופן כללי יותר, R(t) הינה פונקציה סטוכסטית מוגדרת היטב, כך שנוכל להגדיר את משוואת לנז'בן גם בצורה לא ליניארית. כל עוד R(t) אינה סינגולרית- ניתן לקחת את הגבול ללא בעיה. כאשר פונקציית הדלתא בקשר הקורלציה, אשר מפורשת עבור קפיצה קטנה אך לא אינסופית, ניתן להשתמש באינטגרלי סטרטונוביץ'. עבור שימוש באינטגרלי איטו, יש לדרוש מפורשות שפרק זמן זה הינו אפס. (למידע נוסף ראה הלמה של איטו).

יתרונות

ישנם יתרונות לשימוש בגישת לנז'בן על פני שיטות דומות, כגון משוואת פוקר-פלאנק:

1. באופן כללי, קל יותר להבין באופן אינטואיטיבי את שיטה זו, שכן היא מבוססת על התפתחות בזמן של משתנה אקראי- במקום התפתחות בזמן של פונקציית צפיפות הסתברות של משתנה כלשהו.
2. לפעמים יש קושי מסוים בהפרדת משתנים במשוואות כגון משוואת פוקר-פלאנק, ויש צורך בהרחבה עם פונקציות אורתוגונליות המהוות פתרון (למשל פולינומי פורייה -הרמיט)

ראו גם

• משוואת מסטר
• משוואת פוקר-פלאנק
• תהליך ווינר
• תהליך סטוכסטי

ביבליוגרפיה

• The Langevin Equation: With Applications to Stochastic Problems in Physics, Chemistry and Electrical Engineering (Third edition), World Scientific Series in Contemporary Chemical Physics - Vol 27.
• van Kampen, Stochastoc Processes in Physics and Chemistry
• Progress of Theoretical Physics Supplement, Volume 130, January 1998, Pages 17–27,https://doi.org/10.1143/PTPS.130.17
• ; reviewed by D. S. Lemons & A. Gythiel: Paul Langevin’s 1908 paper "On the Theory of Brownian Motion" [...], Am. J. Phys. 65, 1079 (1997),
• Reif, F. Fundamentals of Statistical and Thermal Physics, McGraw Hill New York, 1965. See section 15.5 Langevin Equation

הערות שוליים

הערך באדיבות ויקיפדיה העברית, קרדיט,
רשימת התורמים
רישיון cc-by-sa 3.0

31418964משוואת לנז'בן