Ve fyzice, zejména v oblasti teoretické astrofyziky a obecné relativity, existují složité matematické struktury, které nám umožňují zjednodušit popis a analýzu interakcí mezi gravitací a elektromagnetismem. Jedním z takových nástrojů je použití komplexního potenciálu, který se ukazuje jako užitečný pro vyjádření rovnic gravitačního a elektromagnetického pole, zejména v prázdném prostoru. Tento přístup vychází z analýzy Maxwellových rovnic a Einsteinových rovnic v kontextu teorie relativity.

Komplexní potenciál je definován jako kombinace elektromagnetického potenciálu A0A_0 a magnetického potenciálu φ\varphi, přičemž pro tento účel zavedeme komplexní skalární pole, které se vyjadřuje jako:

ψ=A0+iφ,\psi = A_0 + i \varphi,

kde ii je imaginární jednotka. Tento potenciál nám umožňuje spojit elektrostatické a magnetické složky elektromagnetického pole do jednoho komplexního objektu. Důležitým krokem v tomto procesu je použití operace „curl“ na vektorový potenciál ω\omega, který je definován až do gradientu, což znamená, že závisí na volbě počátku časové souřadnice x0x_0. Tato volnost může být eliminována právě použitím operace "curl" v metrice γαβ\gamma_{\alpha \beta}, což vede k vytvoření invariantního vektoru.

Rovnice pro tento vektor, označený jako τ\tau, jsou dány vztahem:

τ=curlω.\tau = -\text{curl} \, \omega.

Tato formulace umožňuje vyjádřit Maxwellovy rovnice v jednodušší formě, kde místo složitých komponent elektromagnetického pole získáme vztahy v závislosti na komplexním potenciálu. Maxwellovy rovnice pak mohou být zjednodušeny na jedinou rovnici s parciálními derivacemi pro komplexní skalární pole ψ\psi, což dává:

div(f1)ψ=f2τψ,\text{div} \left( f^{ -1} \nabla \right) \psi = f^{ -2} \tau \cdot \nabla \psi,

kde ff je funkce závislá na metrice prostoru, a τ\tau je vektor definovaný výše. Tento přístup vede k tomu, že celé soustavy rovnic pro elektromagnetické pole a gravitační pole mohou být vyjádřeny prostřednictvím komplexních potenciálů a metriky, což činí analytické řešení problémů efektivnějšími a přehlednějšími.

Dalším krokem je využití Ricciho tenzoru RμνR_{\mu \nu} ve vztahu k vektoru GG, který je definován jako:

G=f+iτ,G = \nabla f + i \tau,

a jeho komponenty mohou být vyjádřeny pomocí diferenciálních operací a operace „curl“. Tento přístup vede k vytvoření rovnic pro energii elektromagnetického pole, které jsou ve formě:

Tμν=14π(gμλFαβFνγ),T_{\mu \nu} = \frac{1}{4\pi} \left( g_{\mu \lambda} F_{\alpha \beta} F_{\nu \gamma} - \ldots \right),

kde TμνT_{\mu \nu} je tenzor energie a hybnosti elektromagnetického pole. Výraz pro tento tenzor lze dále upravit podle konkrétní metriky a použít pro získání analytických řešení v konkrétních geometrických a fyzikálních podmínkách.

Významným přínosem tohoto přístupu je, že umožňuje zjednodušení výpočtů v situacích, kde by jinak bylo obtížné pracovat s klasickými rovnicemi Maxwellovy nebo Einsteinovy teorie. Komplexní potenciály poskytují možnost sjednocení gravitačního a elektromagnetického pole do jedné formy, což vede k celkově elegantnějším a efektivnějším metodám řešení.

Další důležitý bod, který čtenář musí chápat, je to, že volba metriky γαβ\gamma_{\alpha \beta} není trivialní. Metrika musí být kompatibilní s rovnicemi pro čistě gravitační pole, jinak nelze použít tento formalismus. Při aplikaci na konkrétní problémy je nutné pečlivě volit vhodné metriky, které odpovídají konkrétním fyzikálním podmínkám, jak ukazují různé metody řešení v teorii relativity. Metody jako je použití komplexních potenciálů v těchto rovnicích umožňují získat kvalitní analytické výsledky i v případech, kde by běžné metody byly složité nebo neefektivní.

Pokud jde o samotné řešení Einsteinových rovnic pro gravitační pole, tento přístup ukazuje, že rovnice mohou být zjednodušeny a řešení může být získáno analýzou a manipulací s těmito komplexními potenciály. Důležité je, že řešení rovnic ve vakuu pro gravitační pole může být generováno na základě známých výsledků, což poskytuje efektivní způsob, jak rozšířit teorie a aplikace v astrofyzice a kosmologii.

Jak můžeme chápat napětí v materiálech a jeho vzorcích?

Napětí v materiálech lze vidět jako koncept, který přesahuje aplikovanou sílu, což je prostředek k eliminaci síly jako takové z koncepčního arzenálu mechaniky. Napětí aplikované zvenčí je úzce spojeno s ideou síly. Jinak, pokud se vzdálíme od konceptu síly, lze napětí považovat také za hustotu energie, která charakterizuje hmotu, občas i nezávisle na jakékoli síle. Hlavním problémem je zde najít funkci, která implicitně obsahuje fyzikální povahu kontinua, k němuž se vztahují napětí a deformace.

Při hlubší reflexi se objeví specifická vlastnost: zde se setkáváme s nekontrolovatelnými projevy. To je snad hlavní důvod, proč mechanika stále udržuje myšlenkový rámec, který považuje sílu za vektor. Je pravda, že jakákoli naše akce je vykonána silou. Jinými slovy, síly jsou jediným prostředkem, jak můžeme něco ovládat, pokud to má být něco jiného než síla. S výjimkou několika vzácných případů nejsou v mechanice nekontrolovatelné veličiny explicitně brány v úvahu, protože udržování sil je teoreticky nezbytným nástrojem.

Avšak tato myšlenka nekontrolovatelnosti se výborně spojuje s konstitutivním zákonem, který bychom nazvali přirozeným. Tento zákon, který vztahuje napětí k deformacím, má formu:

σ=p0e+p1ε+p2ε2\sigma = p_0 e + p_1 \varepsilon + p_2 \varepsilon^2
kde ee je jednotková matice 3×3. Tento zákon nazýváme přirozeným nejen proto, že by se realizoval v přírodě, ale protože přirozeně vyplývá z našich reprezentací napětí a deformací. V našich modelech napětí a deformací jsou tyto veličiny vyjádřeny jako matice 3×3, a pokud je konstitutivní zákon analytický, rovnice (6.85) musí automaticky platit.

Pak lze vztah mezi těmito maticemi zapsat jako celou řadu v matici deformace, která je vždy zjednodušitelná na druhý stupeň polynomu podle Hamilton-Cayleyovy věty. Vztah může být také dobře napsán s prohozenými místy napětí a deformací, takže deformace zůstávají kvadratickými funkcemi napětí, pouze s jinými koeficienty. Takže pro materiál existují nejen tři charakteristické čísla, ale i šest: tři pro vyjádření napětí v závislosti na deformacích a tři pro vyjádření deformací v závislosti na napětích.

Tento rámec tedy umožňuje materiálu mít alespoň zdánlivou přesnou identitu, protože ho můžeme identifikovat pomocí koeficientů p0p_0, p1p_1, p2p_2, které jsou přístupné experimentálně. To je vlastně to, co se obvykle rozumí pod "charakterizací materiálu". Příliš často se však v experimentální praxi tyto koeficienty považují za vlastnosti čistoty materiálu, jako je hustota nebo magnetismus, což vede k záměně pojmů, někdy s vážnými důsledky, zejména v inženýrských problémech.

Tato skutečnost ukazuje, že každý z těchto materiálových parametrů lze získat z experimentů, při nichž je materiál zatěžován, buď v tahu, nebo v tlaku. Bez ohledu na povahu tohoto zatížení se hlavní směry napětí shodují s hlavními směry deformací. Pokud jsou σ1,σ2,σ3\sigma_1, \sigma_2, \sigma_3 hlavní hodnoty matice napětí a ε1,ε2,ε3\varepsilon_1, \varepsilon_2, \varepsilon_3 hlavní hodnoty matice deformace, pak je konstitutivní zákon (6.85) ekvivalentní systému:

σ1=p0+p1ε1+p2ε12,\sigma_1 = p_0 + p_1 \varepsilon_1 + p_2 \varepsilon_1^2,
σ2=p0+p1ε2+p2ε22,\sigma_2 = p_0 + p_1 \varepsilon_2 + p_2 \varepsilon_2^2,
σ3=p0+p1ε3+p2ε32.\sigma_3 = p_0 + p_1 \varepsilon_3 + p_2 \varepsilon_3^2.

Pro argumentaci bychom mohli provést experimenty, které by nám umožnily současně měřit všechny tři hlavní hodnoty deformací a napětí. Takový experiment není prakticky proveditelný, ale teoretický argument naznačuje, že je to vždy možné. Výsledky těchto experimentů nám umožní spočítat vlastnosti materiálu, které jsou obsaženy v koeficientech p0,p1,p2p_0, p_1, p_2 z uvedeného systému. Protože materiál je jedinečný, čekáme jedinečné řešení systému, které je získáno pouze tehdy, pokud je jeho determinant nenulový.

I když jsou tyto koeficienty jedinečné a velmi vhodné pro charakterizaci materiálů, získané hodnoty nejsou čistě materiálovými vlastnostmi, protože závisí na stavu napětí, které bylo na materiál aplikováno. Jsme tedy nuceni specifikovat, co vlastně znamenají čisté materiálové vlastnosti, a to je, jak tomu vždy bylo, problém. Tento problém lze vyřešit tím, že si uvědomíme, že formalismus ukazuje, že deformace se vyskytují i v případě, kdy na materiál nejsou aplikována žádná napětí. Protože neznáme jejich původ, můžeme je oprávněně považovat za vnitřní vlastnosti, tedy čisté, materiálu. Nazvěme je však pouze vnitřními, protože by mohly být generovány silami, jejichž přítomnost v současnosti neznáme.

V tomto případě můžeme tuto skutečnost popsat následujícím systémem rovnic:

0=p0+p1ε1+p2ε12,0 = p_0 + p_1 \varepsilon_1 + p_2 \varepsilon_1^2,
0=p0+p1ε2+p2ε22,0 = p_0 + p_1 \varepsilon_2 + p_2 \varepsilon_2^2,
0=p0+p1ε3+p2ε32.0 = p_0 + p_1 \varepsilon_3 + p_2 \varepsilon_3^2.

Vnitřní charakterizace materiálu experimentem je tedy nyní delegována na hledání řešení tohoto lineárního a homogenního systému, pokud existují. Ve skutečnosti vždy existují, musíme pouze rozhodnout, kolik jich je, a to závisí na tom, co jsme schopni měřit v praxi. Pokud vždy měříme tři různé deformace ve třech ortogonálních směrech prostoru, materiál na aplikovaná napětí nereaguje. To už víme z naší historie: je to primární vlastnost éteru, který šíří světlo okamžitě. Pokud však současné měření tří hlavních hodnot matice je možné pouze v naší představivosti, kvalita materiálu musí být stejně produktem naší představivosti. V reálném světě můžeme současně měřit nanejvýš dvě vlastní hodnoty, které nám ukazují, že materiál stále může reagovat na napětí, jinými slovy, jeho deformace skutečně provází napětí.

Pokud tedy měříme jednu a tutéž hodnotu deformace v jakémkoli směru prostoru, máme dvojitou nekonečnost stavů napětí materiálu, závislých na dvou materiálových parametrech. Pokud měříme dvě hodnoty deformace, například v jednom směru a v jeho kolmé rovině, pak máme stavy napětí materiálu závislé na jediném materiálovém parametru. Předpokládejme, že můžeme zahrnout jeden z materiálových parametrů do měřené hodnoty, nejobecnější konstitutivní zákon, který je splněn materiálem vykazujícím napětí při deformacích, bude mít tvar:

σ=K(εε1e)(εε2e),\sigma = K(\varepsilon - \varepsilon_1 e)(\varepsilon - \varepsilon_2 e),

kde KK je libovolná konstanta, která má rozměry napětí. Takový materiál má tři nekontrolovatelné veličiny, z nichž dvě jsou měřitelné.

Pokud se soustředíme pouze na měřitelné veličiny, pohodlný způsob

Jak popsat materiály bez deformací při napětí

Pokud se rozhodneme charakterizovat vnitřní vlastnosti materiálu jako deformace, musíme vzít v úvahu, že tato konvence je nezávislá na materiálovém popisu a musí být zajištěna naší měřicí schopností. Proto můžeme dojít k následujícímu obecnému závěru: kdykoli materiál volně deformuje, tj. bez působení zjevné síly, musí být jeho deformační matice ve formě (6.90), přičemž všechny speciální případy jsou zahrnuty. Deformace, stejně jako doprovázející napětí, se pak projeví jako ortogonální tensory.

Dále můžeme diskutovat o kategorii materiálů, které jsou schopny udržovat napětí bez odpovědi v podobě deformací. Tento jev je v podstatě klíčovým vjemem pro naše smysly. Pro vyjádření tohoto jevu musíme zvážit opačný zákon, totiž vztah:

ε=q0e+q1σ+q2σ2.(6.91)\varepsilon = q_0 e + q_1 \sigma + q_2 \sigma^2. \tag{6.91}

V tomto případě může být σ\sigma nazváno napětím pouze zneužitě; ve skutečnosti představuje hustotu vnitřní energie materiálu. Stav tohoto éteru bude charakterizován systémem rovnic:

0=q0+q1σ1+q2σ120 = q_0 + q_1 \sigma_1 + q_2 \sigma_1^2
0=q0+q1σ2+q2σ220 = q_0 + q_1 \sigma_2 + q_2 \sigma_2^2
0=q0+q1σ3+q2σ320 = q_0 + q_1 \sigma_3 + q_2 \sigma_3^2

Tento systém odpovídá situaci, ve které není pozorována žádná deformace éteru. Opět je charakterizace tohoto éteru závislá na počtu netriviálních řešení systému, přičemž nejběžnější forma deformační matice je zde:

ε=K11(σσ1e)(σσ2e).(6.92)\varepsilon = K^{ -1}_1 (\sigma - \sigma_1 e)(\sigma - \sigma_2 e). \tag{6.92}

Kde konstanta K1K_1 má rozměry napětí. Bell [16] uvedl, že vztah (6.92) s σ1+σ2=0\sigma_1 + \sigma_2 = 0 je charakteristický pro kovy, a to jak při malých, tak i velkých deformacích: kovy vždy ovlivňují svou tvrdost. Pokud se tedy zaměříme pouze na měřitelné veličiny pro charakterizaci tohoto materiálu, jeho vnitřní napětí budou mít následující pohodlné tensoriální zobrazení, analogické rovnici (6.91):

σij=σ2δij+(σ1σ2)mimj,i,j=1,2,3.(6.93)\sigma_{ij} = \sigma_2 \delta_{ij} + (\sigma_1 - \sigma_2) \cdot m_i m_j, \quad i,j = 1, 2, 3. \tag{6.93}

Kde mm je jednotkový vektor odpovídající hlavní hodnotě σ1\sigma_1. Lze tedy oprávněně říci, že obecná vlastnost materiálu, který neprojevuje deformace při napětí, je zakotvena v podobě (6.92), přičemž všechny konkrétní případy jsou zahrnuty.

Případ rovnic (6.90) a (6.93) je specifický pro tensory, které bychom zde rádi nazvali ekvivalentními vektorovému poli. Tuto ekvivalenci chápeme následujícím způsobem: nechť vv je vektorové pole, pomocí něhož sestavíme matici:

vij=αδij+βvivj.(6.94)v_{ij} = \alpha \delta_{ij} + \beta v_i v_j. \tag{6.94}

Je zřejmé, že pokud vkv_k jsou složky vektoru a α\alpha a β\beta jsou skaláry, pak vijv_{ij} automaticky představují složky tensora. Jedna z jeho hlavních hodnot, totiž α\alpha, je dvojnásobná. Druhá hlavní hodnota, která se liší od α\alpha, je dána vzorcem:

α=α+βv2.(6.95)\alpha' = \alpha + \beta v^2. \tag{6.95}

Je zajímavé si povšimnout několika vlastností tohoto tensora. Prvně, pokud je buď β\beta, nebo vkv_k nulový, vijv_{ij} je čistě kulový tensor. Dále, pokud vypočteme vlastní vektor vv odpovídající vlastní hodnotě (6.95), zjistíme, že je to právě vektor vv, až na normalizační faktor. Tato vlastnost je nezávislá na parametru α\alpha, a právě tato skutečnost nám umožňuje definovat zmíněnou ekvivalenci: pro daný vektor vv můžeme přímo konstrukčně vytvořit vektor jako rodinu tensorů závislých na dvou libovolných parametrech, které mají tento vektor jako vlastní vektor.

Je nutné si uvědomit, že hlavní výhodou existence takového konstitutivního zákona je, že nám umožňuje vyjádřit, co jsme právě řekli o fraktálním médiu ve prostoru a fraktálním médiu v materiálu, v algebrické podobě. Poté je zřejmé, že forma napětí, která není doprovázena deformacemi, je charakteristická pro případ klasického elektrického pole, a tvar deformací, které nejsou doprovázeny napětím, je charakteristický pro klasické magnetické pole, nebo naopak.

Nicméně, jelikož neexistuje pouze fraktální médium v materiálu nebo pouze fraktální médium v prostoru, pravé fraktální médium musí být kombinací těchto dvou druhů fraktálních médií. Nejjednodušší kombinace odhaluje matematickou strukturu, která popisuje elektromagnetické záření a tedy světlo.

Fraktální médium můžeme charakterizovat tak, že není popsáno jediným tenzorem ve tvaru (5.38), ale dvěma, odpovídajícími dvěma charakteristickým vektorům, řekněme uu a ν\nu. Tenzor, který by potom popisoval fraktální médium, by mohl mít tvar:

wij=αδij+βuiuj+γvivj.(6.96)w_{ij} = \alpha \delta_{ij} + \beta u_i u_j + \gamma v_i v_j. \tag{6.96}

Tento vztah je možno přepsat do symetričtější formy:

wij=λuij+μvij,(6.97)w_{ij} = \lambda u_{ij} + \mu v_{ij}, \tag{6.97}

kde λ\lambda a μ\mu jsou reálné parametry, které popisují "prostorové" a "materiální" složky fraktálního média, přičemž matice uu a vv jsou definovány následovně:

uij=uiuju23δij,vij=vivjv23δij.(6.98)u_{ij} = u_i u_j - \frac{u^2}{3} \delta_{ij}, \quad v_{ij} = v_i v_j - \frac{v^2}{3} \delta_{ij}. \tag{6.98}

Pokud bychom tuto matici rozšířili, získáme následující výsledek:

wij=λuiuj+μvivj13(λu2+μv2)δij.(6.99)w_{ij} = \lambda u_i u_j + \mu v_i v_j - \frac{1}{3} (\lambda u^2 + \mu v^2) \delta_{ij}. \tag{6.99}

Tento tenzor má tři reálné a různé hlavní hodnoty, jak ukazuje charakteristická rovnice matice:

w1=e,w2,3=±e2g2.(6.102)w_1 = e, \quad w_2,3 = \pm \sqrt{e^2 - g^2}. \tag{6.102}

Pochopení těchto vztahů nám umožňuje lépe porozumět interakcím v materiálech a jejich chování v různých podmínkách.

Jak Lyapunovův exponent závisí na H a ΩB?

Lyapunovův exponent je klíčovým nástrojem pro studium stability dynamických systémů. Je to hodnota, která vyjadřuje rychlost konvergence nebo divergence trajektorií v fázovém prostoru, a jejíž velikost závisí na počátečních podmínkách. V zásadě určuje, jak rychle se dva velmi blízké body v počátečním stavu systému oddělují. Počet Lyapunovových exponentů odpovídá dimenzi fázového prostoru. Nejvyšší exponent, pokud je kladný, obvykle určuje nejdůležitější dynamiku systému, jelikož indikuje vysokou citlivost na počáteční podmínky, což je charakteristické pro chaotické systémy. Tento charakteristický rys je klíčový pro porozumění chaosu, jelikož znamená, že i malé změny v počátečních podmínkách mohou vést k dramatickým změnám v chování systému.

Při analýze Lyapunovových exponentů je jedním z přístupů definice exponentu prostřednictvím malé změny v počáteční pozici. Tato změna, označená jako δx0, se vypočítává jako rozdíl mezi dvěma velmi blízkými počátečními stavy, x1(0) a x2(0). V čase t se tato změna vyjadřuje jako δx0(t) = x2(t) − x1(t). V našem výzkumu jsme použili tuto definici exponentu a prováděli numerické simulace pro případy, kdy t/T = 10³, což představuje dostatečně dlouhý časový interval pro stabilní měření, a δx0 byla dostatečně malá, aby bylo možné analyzovat chování systému. Po analýze s vyšší přesností jsme nezjistili žádné zásadní změny v výsledcích, což ukazuje na stabilitu numerických metod při dostatečné přesnosti.

Bifurkační diagram je dalším důležitým nástrojem pro studium přechodu k chaosu. V tomto případě se systematicky pozoruje evoluce trajektorií směrem k chaotickému chování prostřednictvím překrývání rezonancí. Tento proces je doprovázen vznikem stochastických vrstev, ve kterých částice vykonávají náhodné pohyby, a které jsou odděleny invariántními křivkami. Některé z těchto křivek mohou zaniknout, což vede k vytváření složité fraktální struktury v fázovém portrétu systému. Tento jev je příkladem bifurkace, kdy se stabilní oblasti, reprezentující periodické chování systému, transformují do oblastí s chaotickým pohybem. Tento přechod k chaosu se může opakovat v podobě kaskády bifurkací, což je typické pro Hamiltonovské systémy.

Při analýze bifurkačních diagramů je nutné určit frekvence oscilací, které jsou funkcemi parametru H. Pro numerickou analýzu jsme použili metody, které mapují hodnoty X, Px a Py vůči času. Tento proces zahrnuje hledání maximálních hodnot v určitém časovém intervalu a identifikaci frekvenčních změn, což umožňuje rozpoznat fenomén dvojnásobení frekvencí a s ním související chaotické chování. Fourierova transformace nám poté pomáhá identifikovat oblasti, kde dochází k překrývání rezonancí, což značí přechod do chaotického režimu.

Dalším zajímavým aspektem chování systému je multifraktální analýza. Z této analýzy vyplývá, že pohyb částice není pouze chaotický, ale vykazuje vlastnosti fraktální geometrie. Trajektorie pohybu částice jsou popsány křivkami, které jsou kontinuální, ale nejsou diferenciovatelné – tedy fraktálními křivkami. Tato vlastnost křivek je definována jejich samo-podobností na různých měřítkách. Tato samo-podobnost je charakteristickým rysem fraktálních struktur a je jedním z hlavních znaků, které umožňují modelovat trajektorie pohybu částic v systému.

V předchozích studiích jsme ukázali, že pohyb částic se odehrává na křivkách s konstantní fraktální dimenzí, která je rovna hodnotě D = 1,5. Tato dimenze je zajímavá i v širším kosmologickém kontextu, protože novější studie ukazují, že hmota ve vesmíru je rozložena fraktálně, s dimenzí, která se asymptoticky blíží právě hodnotě 1,5. Tato zjištění ukazují na vysokou koherenci mezi teoretickými modely a experimentálními pozorováními. Fraktální analýza tedy poskytuje důležité nástroje pro pochopení dynamických vlastností systémů, které vykazují složité nelineární interakce, jako jsou ty, které zkoumáme v rámci pohybu částic v gravitačně elektromagnetických polích.

Je třeba si uvědomit, že přechod od lineárního k nelineárnímu chování, který je pozorován v tomto typu dynamických systémů, není vždy snadno předvídatelný. Zatímco ve fázi lineárního chování je možné použít běžné metody analýzy, jako je Fourierova transformace, jakmile se systém dostane do chaotického režimu, tyto metody již nejsou dostatečné. Zde je nezbytné aplikovat složitější přístupy, jako je multifraktální analýza, která dokáže zachytit složitost a samo-podobnost trajektorií, což jsou klíčové vlastnosti pro studium chaotických systémů.