Această lucrare investighează apariția redundanței structurale... Prezentăm un cadru structural pentru apariția realității clasice din sistemele cuantice, bazat pe o tranziție în redundanță determinată de dimensiune. Derivăm o condiție analitică exactă pentru pragul critic de redundanță D_c(N), definit printr-o ecuație transcendentală independentă de puterea de cuplare. Aceasta stabilește că tranziția de la regimurile neredundante (T₀) la cele redundante (T₁) este guvernată de dimensionalitatea sistemului mai degrabă decât de parametrii dinamici. Introducem un prag critic de energie ambientală E_c(N), demonstrat numeric a fi invariant la variațiile forței de cuplare (CV < 2%) pentru interacțiuni de tip defazare (σ_z). Acest lucru demonstrează că debutul redundanței are loc numai după ce o cantitate finită de energie a trecut în mediu. Rezultatele stabilesc o ordine structurală: fluxul de energie precede redundanța stabilă a informației. În limita termodinamică, E_c(N) converge către D_c(N), sugerând un prag unificat care guvernează atât energia, cât și informația la sisteme de dimensiuni mari. Demonstrăm în continuare că acest comportament este specific interacțiunilor de defazare. Sub cuplaj σ_x, invarianța lui E_c(N) se descompune, indicând o structură dinamică fundamental diferită. Aceste descoperiri oferă o nouă perspectivă asupra tranziției de la cuantic la clasic, identificând pragurile dimensionale și energetice ca principalii factori ai emergenței clasice.
Această lucrare investighează apariția redundanței structurale... Prezentăm un cadru structural pentru apariția realității clasice din sistemele cuantice, bazat pe o tranziție în redundanță determinată de dimensiune. Derivăm o condiție analitică exactă pentru pragul critic de redundanță D_c(N), definit printr-o ecuație transcendentală independentă de puterea de cuplare. Aceasta stabilește că tranziția de la regimurile neredundante (T₀) la cele redundante (T₁) este guvernată de dimensionalitatea sistemului mai degrabă decât de parametrii dinamici. Introducem un prag critic de energie ambientală E_c(N), demonstrat numeric a fi invariant la variațiile forței de cuplare (CV < 2%) pentru interacțiuni de tip defazare (σ_z). Acest lucru demonstrează că debutul redundanței are loc numai după ce o cantitate finită de energie a trecut în mediu. Rezultatele stabilesc o ordine structurală: fluxul de energie precede redundanța stabilă a informației. În limita termodinamică, E_c(N) converge către D_c(N), sugerând un prag unificat care guvernează atât energia, cât și informația la sisteme de dimensiuni mari. Demonstrăm în continuare că acest comportament este specific interacțiunilor de defazare. Sub cuplaj σ_x, invarianța lui E_c(N) se descompune, indicând o structură dinamică fundamental diferită. Aceste descoperiri oferă o nouă perspectivă asupra tranziției de la cuantic la clasic, identificând pragurile dimensionale și energetice ca principalii factori ai emergenței clasice.
MONOGRAFIE MATEMATICĂ, VARIAȚIONALĂ ȘI STATISTICĂ: MECANISMUL CAUZAL AL EXPERIMENTULUI CELOR DOUĂ FANTE, DINAMICA DE GHIDAJ HIDROSTATIC ȘI TERMODINAMICA SUB-CUANTICĂ ÎN MODELUL SIGMA-PDE Autor: Ilie Barbu Afiliere: Cercetător Independent Metrologie de Referință: CODATA 2022 / SI revizuit (2019) Precizie Aritmetică Internă: 85 de cifre zecimale semnificative (Decimal) Data: 1 Septembrie 2026 REZUMAT Prezenta lucrare academică dezvoltă tratarea analitică, variațională, dimensională și termodinamic-statistică a mecanismului cauzal determinist al experimentului celor două fante și a dinamicii de ghidaj a materiei barionice în cadrul Modelului \Sigma-PDE (Sigma – Presiune Densitate Elastodinamică). Modelul oferă o rezolvare mecanică deterministă a dualității undă-particulă, eliminând ipoteza dedublării fizice a particulei materiale: particula (un defect topologic chiral/vortex localizat în densitatea de masă \rho_{BI}) traversează fizic o singură fantă, în timp ce unda sa pilot de presiune hidrostatică \mathbf{u}_L (emisă continuu și propagată cu viteza superluminică c_L = \sqrt{7/3}\,c \approx 1.527525\,c) ajunge înainte și traversează simultan ambele fante. Monografia formalizează lanțul deductiv complet al celor 6 Pași Fundamentali: 1. Identificarea Structurală Madelung: Dilatația volumică \theta = \nabla\cdot\mathbf{u} a continuumului este identificată direct cu amplitudinea funcției de undă (R = R_0\sqrt{1-\theta}), generând potențialul cuantic Q_\Sigma \approx \frac{\hbar_\sigma^2}{4m}\nabla^2\theta și forța de ghidaj mecanică \mathbf{F}_{\mathrm{ghid}} = -e\rho_{BI} c^2 \nabla\theta - \nabla Q_\Sigma, verificată exact în unități Newton (MLT^{-2}) sub Postulatul 0 ([e] = \text{m}^2); 2. Spațierea Franjelor și Lagrangianul Rezonant: Se deduce spațierea clasică de Broglie \Delta x = \frac{h_\sigma L}{mvd} printr-un Lagrangian de interacțiune \mathcal{L}_{\mathrm{int}} ce impune o penalizare variațională pe suprafața de rezonanță \omega = c_L \frac{mv}{\hbar_\sigma}, unde h_\sigma \equiv 2\pi\hbar_\sigma = \frac{\pi h}{6}; 3. Termodinamica Statistică Fokker-Planck: Fluctuațiile celulelor Barbu la temperatura T_{\mathrm{vacuum}} \approx 4.63 \times 10^{-65}\text{ K} și vâscozitatea \eta_{BI} \approx 6.71 \times 10^{-31}\text{ Pa}\cdot\text{s} generează coeficienții Langevin corecți dimensional (\gamma = \frac{\eta_{BI} l_{BI}}{m}, D_p = (\eta_{BI} l_{BI}) k_B T_{\mathrm{vacuum}}); 4. Demonstrarea H-Teoremei lui Valentini: Se demonstrează riguros că \frac{dH}{dt} \le 0, garantând relaxarea statistică a distribuției pe spațiul fazelor către densitatea cuantică Born \int f_{\mathrm{eq}}\, d^3p = R^2 = |\Psi|^2; 5. Extensia Multiparticulă: Corelațiile dintre N particule sunt mediate prin N câmpuri de dilatație 3D cuplate prin operatori elastici nelocali în mediul comun; 6. Non-Localitatea Retardată și Cauzalitatea CMB: Propagarea ghidajului are loc cu viteza finită c_L, stabilind un cadru cauzal strict ancorat în Sistemul de Referință Privilegiat al Vidului (fondul cosmic de microunde — CMB), unde avansul temporal pe 1 km este \Delta t_{\text{1km}} = \mathbf{1.151951\ \mu\text{s/km}} (față de timpul de zbor absolut t_L = \mathbf{2.183690\ \mu\text{s}}). 1. CADRUL ONTOLOGIC DUAL ȘI CONSTANTELE DE FOND ALE MEDIULUI \Sigma-PDE 1.1. Câmpurile Constitutive Fundamentale Modelul \Sigma-PDE definește vidul ca un mediu mecanic continuu-discret dual, descris de parametrii fundamentali ai Sistemului Internațional revizuit (SI 2019) [80]: - Presiunea hidrostatică statică de fond (Timpul / Elasticitatea): P_{BI} \equiv \frac{1}{\varepsilon_0} = 1.129409066602570332948816986 \times 10^{11}\text{ Pa} - Densitatea masică de fond (Spațiul / Inerția): \rho_{BI} \equiv \mu_0 = 1.256637061262814754000000000 \times 10^{-6}\text{ kg/m}^3 - Impedanța mecanică de forfecare: Z_{BI} \equiv \sqrt{P_{BI}\rho_{BI}} = 2 R_K \alpha = 376.73031340987581914\text{ Pa}\cdot\text{s/m}\ (\Omega) - Operatorul 3D de omogenizare metrică: \Phi_{\dim} \equiv \frac{\Phi_{\mathrm{geom}}}{V_{\mathrm{ref}}} = \frac{\pi}{6}\text{ m}^{-3} \approx 0.523598775598298873077\text{ m}^{-3} \quad (V_{\mathrm{ref}} \equiv 1\text{ m}^3) 1.2. Cinematica Undelor în Opțiunea A (K = \mu = P_{BI}, \, \lambda = P_{BI}/3) Ecuația elastodinamică Navier-Cauchy a mediului: \rho_{BI}\ddot{\mathbf{u}} = P_{BI}\nabla^2\mathbf{u} + \frac{4}{3}P_{BI}\nabla(\nabla\cdot\mathbf{u}) separă strict două moduri fundamentale de undă [80]: 1. Modul Transversal de Forfecare (Foton / EM): c_T = \sqrt{\frac{P_{BI}}{\rho_{BI}}} \equiv c = 299\,792\,458\text{ m/s} \quad \text{[exact SI]} 2. Modul Longitudinal de Compresiune (Unda Scalară de Presiune): \boxed{c_L \equiv \sqrt{\frac{7}{3}}\,c = \mathbf{457\,940\,543.853956492343255620542573949\text{ m/s}} = \mathbf{457\,940.54385\text{ km/s}}} 1.3. Parametrii Sub-Cuantici Celulari - Lungimea fundamentală Barbu: l_{BI} = 1.7821742798263520729 \times 10^{-33}\text{ m} - Masa celulei elementare de rețea: m_{\mathrm{cell}} = \rho_{BI} l_{BI}^3 = 7.11312398968047819 \times 10^{-105}\text{ kg} - Vâscozitatea dinamică Kelvin-Voigt: \eta_{BI} = Z_{BI} l_{BI} = 6.713990750071801631759694242 \times 10^{-31}\text{ Pa}\cdot\text{s} - Temperatura de fluctuație a vidului: T_{\mathrm{vacuum}} = \frac{P_{BI} l_{BI}^3}{k_B} = 4.630399922728312406122364139 \times 10^{-65}\text{ K} - Constantele de acțiune cuantică de rețea: \hbar_\sigma \equiv \hbar \cdot \Phi_{\dim} = \left(\frac{h}{2\pi}\right)\left(\frac{\pi}{6}\right) = \frac{h}{12} = \mathbf{5.521725125000000000000000000 \times 10^{-35}\text{ J}\cdot\text{s}} h_\sigma \equiv 2\pi\hbar_\sigma = \frac{\pi h}{6} = \mathbf{1.734707119298453472049015707 \times 10^{-34}\text{ J}\cdot\text{s}} 2. PASUL 1: IDENTIFICAREA DILATAȚIEI CU AMPLITUDINEA MADELUNG ȘI FORȚA DE GHIDAJ 2.1. Câmpul de Dilatație și Conservarea Masei Fie \mathbf{u}_L(\mathbf{x}, t) câmpul de deplasare longitudinal (\nabla \times \mathbf{u}_L = 0). Dilatația volumică a continuumului este definită ca: \theta(\mathbf{x},t) \equiv \nabla \cdot \mathbf{u}_L(\mathbf{x},t) = \partial_k u_k Din ecuația de continuitate liniarizată a mediului \Sigma-PDE în jurul stării de fond: \partial_t \rho_{\mathrm{mediu}} + \rho_{BI}\nabla\cdot\dot{\mathbf{u}}_L = 0 \implies \rho_{\mathrm{mediu}}(\mathbf{x},t) = \rho_{BI}\left(1 - \theta(\mathbf{x},t)\right) \theta(\mathbf{x},t) = 1 - \frac{\rho_{\mathrm{mediu}}(\mathbf{x},t)}{\rho_{BI}} 2.2. Identificarea Structurală a Amplitudinii Madelung R(\mathbf{x},t) În hidrodinamica cuantică a lui Madelung, funcția de undă este \Psi = R e^{i S / \hbar_\sigma}, cu densitatea de probabilitate \rho_{\mathrm{prob}} = R^2 [110]. În \Sigma-PDE, dilatația mediului este însăși amplitudinea cuantică: \boxed{R(\mathbf{x},t) \equiv R_0 \sqrt{1 - \theta(\mathbf{x},t)} = R_0 \sqrt{\frac{\rho_{\mathrm{mediu}}(\mathbf{x},t)}{\rho_{BI}}}} 2.3. Potențialul Cuantic Madelung Q_\Sigma Potențialul cuantic generat de deformația elastică a mediului este: Q_\Sigma \equiv -\frac{\hbar_\sigma^2}{2m} \frac{\nabla^2 R}{R} = -\frac{\hbar_\sigma^2}{2m} \frac{\nabla^2 \sqrt{1-\theta}}{\sqrt{1-\theta}} La ordinul liniar în perturbația elastică (\theta \ll 1): \sqrt{1-\theta} \approx 1 - \frac{1}{2}\theta \implies \nabla^2\sqrt{1-\theta} \approx -\frac{1}{2}\nabla^2\theta \frac{\nabla^2\sqrt{1-\theta}}{\sqrt{1-\theta}} \approx -\frac{1}{2}\nabla^2\theta \implies \boxed{Q_\Sigma \approx \frac{\hbar_\sigma^2}{4m} \nabla^2\theta} 2.4. Deducerea Forței Totale de Ghidare \mathbf{F}_{\mathrm{ghid}} Particula materială este un vortex chiral închis având sarcina efectivă e (arie cinematică [e] = \text{m}^2 sub Postulatul 0). Asupra sa acționează două forțe mecanice conjugate: 1. Forța de gradient de presiune elastică: \mathbf{F}_{\mathrm{elastic}} = -e\nabla P_{\mathrm{elastic}} = -e\nabla(-P_{BI}\theta) = +e P_{BI}\nabla\theta = +e\rho_{BI} c^2 \nabla\theta 2. Forța de gradient de potențial cuantic: \mathbf{F}_{\mathrm{cuantic}} = -\nabla Q_\Sigma = -\frac{\hbar_\sigma^2}{4m}\nabla(\nabla^2\theta) Forța totală de ghidaj exercitată pe particulă este: \boxed{\mathbf{F}_{\mathrm{ghid}} = -e\rho_{BI} c^2 \nabla\theta - \frac{\hbar_\sigma^2}{4m}\nabla(\nabla^2\theta)} Verificare Dimensională Exactă: - Primul termen: [e\rho_{BI} c^2 \nabla\theta] = L^2 \cdot (ML^{-3}) \cdot (L^2 T^{-2}) \cdot L^{-1} = MLT^{-2} = \mathbf{Newton\ (N)}; - Al doilea termen: \left[\frac{\hbar_\sigma^2}{4m}\nabla^3\theta\right] = \frac{(ML^2 T^{-1})^2}{M} \cdot L^{-3} = \frac{M^2 L^4 T^{-2}}{M L^3} = MLT^{-2} = \mathbf{Newton\ (N)}. 3. PASUL 2: SPAȚIEREA FRANJELOR ȘI LAGRANGIANUL DE INTERACȚIUNE REZONANT 3.1. Relația de Spațiere a Franjelor de Interferență Într-un dispozitiv experimental cu două fante separate de distanța d și ecranul plasat la distanța L, distribuția intensității pe ecran trebuie să reproducă relația de Broglie: \rho_{\mathrm{prob}}(x) \propto \cos^2\left(\frac{\pi d}{\lambda L}x\right) \implies \Delta x = \frac{\lambda L}{d} = \frac{h_\sigma L}{m v d} = \frac{2\pi\hbar_\sigma L}{m v d} = \frac{\pi h L}{6 m v d} 3.2. Postulatul Dinamic de Rezonanță Unda scalară \mathbf{u}_L emisă de particulă se propagă prin vidul \Sigma-PDE cu viteza c_L. Pentru a genera pe ecran franje dependente de impulsul p = mv al particulei, se impune condiția de acordare cinematică de rezonanță: \boxed{\omega = c_L \cdot \frac{m v}{\hbar_\sigma} \iff k = \frac{m v}{\hbar_\sigma} \iff \lambda_\theta = \frac{2\pi}{k} = \frac{2\pi\hbar_\sigma}{m v} = \frac{h_\sigma}{m v}} 3.3. Construcția Variațională: Lagrangianul \mathcal{L}_{\mathrm{int}} Pentru a forța dinamic această condiție în ecuați
OPERA MAGNA — UN SINGUR UNIVERS, O SINGURĂ FIZICĂ VOLUMUL II: MICROFIZICA DETERMINISTĂ, SECTORUL HADRONIC ȘI CROMODINAMICA ELASTODINAMICĂ Tratat de Hidrodinamică Cuantică de Ghidaj, Termodinamică Sub-Cuantică Valentini, Structura Solitonică Skyrme a Protonului și Confinarea Forței Tari Autor: Ilie Barbu Cercetător Independent, Pitești, Argeș, România Cadru Metrologic: CODATA 2022 / SI revizuit (2019) / PDG 2024 / Jefferson Lab (2018) Precizie Computațională: Aritmetică internă de înaltă rezoluție (până la 85 de cifre zecimale) Data: 7 Septembrie 2026 REZUMAT EXECUTIV Prezentul tratat expune formularea analitică, variațională și computațională a microfizicii cuantice și a fizicii nucleare în cadrul Modelului \Sigma-PDE (Sigma – Presiune – Densitate – Elastodinamică). Lucrarea oferă o rezolvare mecanică deterministă pentru paradoxul dualismului undă-particulă, eliminând interpretarea acauzală a colapsului funcției de undă prin intermediul unei hidrodinamici de vid cu dublă celeritate și descriind Cromodinamica Cuantică (QCD) prin legile mecanicii neliniare a continuumului viscoelastic. Pilonii analitici fundamentali ai Volumului II sunt: 1. Dihotomia Cinematică Particulă–Undă Pilot: Particula este un vortex topologic de forfecare 2D ([e] = L^2 = \text{m}^2) localizat în câmpul de densitate \rho_{\text{BI}}, ce traversează o singură fantă cu viteza v < c. Unda pilot este un mod longitudinal de compresiune pură a vidului (\vec{u}_L, \vec{\nabla}\times\vec{u}_L = 0), propagat cu viteza superluminică: c_L = \sqrt{\frac{7}{3}}\,c = \mathbf{457\,940\,543,853956\dots\ \text{m/s} \approx 1,527525\,c} care traversează simultan ambele fante și ghidează particula prin forța hidrostatică exactă: \vec{F}_{\text{ghid}} = - l_{\text{vortex}} e P_{\text{BI}} \vec{\nabla}\theta - \frac{\hbar_\sigma^2}{4m}\vec{\nabla}(\nabla^2\theta) recuperând potențialul cuantic Madelung Q_\Sigma \approx \frac{\hbar_\sigma^2}{4m}\nabla^2\theta din ecuația de continuitate a rețelei. 2. Termodinamica Sub-Cuantică și Demonstrarea Regulii lui Born: Cuplajul dinamic dintre vortex și unda pilot este guvernat de frecvența Zitterbewegung a modului longitudinal \omega_L = \frac{m c_L^2}{\hbar_\sigma}, stabilind rata de amortizare \gamma \equiv \omega_L și timpul de relaxare: \tau_{\text{relax}} \equiv \frac{1}{\gamma} = \frac{\hbar_\sigma}{m c_L^2} = \mathbf{2,890464 \times 10^{-22}\text{ s}} \sim \mathbf{10^{-21}\text{ s}} Împreună cu temperatura de fond a celulei Barbu T_{\text{vacuum}} \equiv \frac{P_{\text{BI}} l_{\text{BI}}^3}{k_B} \approx \mathbf{4,630400 \times 10^{-65}\text{ K}}, ecuația de transport Fokker–Planck este formulată cu drift centrat pe viteza de ghidaj \vec{v}_{\text{ghid}}. Se demonstrează prin calcul integral explicit că variația funcționalei de entropie Valentini satisface \frac{dH}{dt} \le 0, dovedind că Regula lui Born (P = |\Psi|^2) este starea de echilibru termodinamic asimptotic al vidului. 3. Rezolvarea Cauzală a Experimentelor Wheeler și Non-Localitatea Retardată: Corelațiile cuantice la distanță sunt mediate de tensorul de Green acustic la viteza c_L. În experimentele de alegere întârziată (Wheeler Delayed-Choice), modificarea de frontieră a detectorului transmite o undă de rearanjare elastică la viteza c_L > c care reconfigurează forța de ghidaj înainte ca particula să parcurgă interferometrul (t_{\text{undă}} \ll t_{\text{particulă}}), eliminând retrocauzalitatea. Se deduce avansul temporal invariant Time-of-Flight pe 1\text{ km}: \Delta t_{1\text{km}} = 1000\left(\frac{1}{c} - \frac{1}{c_L}\right) = \mathbf{1,151950871670\dots\ \mu\text{s}} și se demonstrează analitic colapsul corelațiilor Bell către limita clasică locală (|S| \le 2) la timpi de comutare sub-retardați \delta t < d/c_L. 4. Cromodinamica Elastodinamică și Structura Protonului: Protonul este demonstrat ca un soliton 3D stabil aflat sub presiunea internă medie: P_{\text{mat}} = \frac{m_p c^2}{\frac{4\pi}{3} R_p^3} = \mathbf{6,024806 \times 10^{34}\text{ Pa}} \equiv \mathbf{376,04\text{ MeV/fm}^3} Sarcina de culoare SU(3)_c este orientarea forfecării pe cele 3 axe ale triedrului celular Barbu, cei 8 gluoni fiind modurile izovolumice deviatorice (\text{Tr}(\varepsilon) = 0). Aplicând saltul de presiune de frontieră Young–Laplace \Delta P = 4 P_{\text{mat}}, presiunea din miezul central al protonului este dedusă riguros la P_{\text{core}} = 5 P_{\text{mat}} = \mathbf{3,012403 \times 10^{35}\text{ Pa}} (în acord cu datele Jefferson Lab Nature 2018). Această presiune comprimă forfecarea într-un tub de flux cilindric de rază \xi \approx 0,4040\text{ fm}, deducând tensiunea de coardă \sigma_{\text{string}} = P_{\text{core}}(\pi\xi^2) = \mathbf{1,544579 \times 10^5\text{ N}} \equiv \mathbf{0,964051\text{ GeV/fm} \approx 1,0\text{ GeV/fm}}. 5. Potențialul Nuclear Unificat și Defectul de Volum: Stabilitatea cavității este asigurată prin termenul cuantic de ordinul 4 Skyrme \kappa_{\text{Sk}}(\partial\Psi)^4, satisfăcând identitatea von Laue \int_0^\infty r^2 [P_r + 2P_t]\,dr = 0. Se deduce Potențialul Nuclear Unificat: U_{\text{tot}}(d) = V_0 \left[ C_{\text{rep}} \left(\frac{\xi}{d - 2\xi}\right)^4 - \left(\frac{\lambda_\pi}{d}\right) e^{-d/\lambda_\pi} \right] unde V_0 \equiv m_\pi c^2 \approx 134,98\text{ MeV}, C_{\text{rep}} = 0,085, \xi = 0,404012\text{ fm} și \lambda_\pi \approx 1,414\text{ fm}. Expresia este riguros omogenă dimensional în Jouli și generează un minim atractiv de -48,16\text{ MeV} la d = 1,4\text{ fm} (cu un minim absolut de -55,8\text{ MeV} la d \approx 1,22\text{ fm}), trecând asimptotic la repulsie dură infinită pentru d \le 2\xi \approx 0,808\text{ fm}. Se demonstrează Legea Universală a Energiei Nucleare: E_B = P_{\text{BI}} \cdot \Delta V_N = \rho_{\text{BI}} c^2 \Delta V_N = \Delta m \cdot c^2 arătând că energia nucleară este lucrul mecanic al presiunii de fond pe volumul structural eliminat (\Delta V_{\text{DT}} \approx 2,50 \times 10^{-23}\text{ m}^3). 1. DECONSTRUCȚIA DUALISMULUI UNDĂ-PARTICULĂ ȘI HIDRODINAMICA DE GHIDAJ +-------------------------------------------------------------------------------+ | TOPOLOGIA GHIDAJULUI ÎN EXPERIMENTUL YOUNG | | | | Sursa de Vortex | | (Particulă) | | o | | / \ | | / \ | | / u_L\ Unda pilot se propagă la c_L = √(7/3) c | | v v | | ====[Fanta A]=====[Fanta B]==== | | | ~ | | Particula | ~ Unda u_L trece prin ambele fante | | trece DOAR| ~ și generează profilul de dilatație θ | | prin A v ~ | | \ / | | v v | | ----------------------- | | Maxime de Franje Born: P(x) = R²(x) ≡ |Ψ(x)|² | +-------------------------------------------------------------------------------+ 1.1 Dihotomia Vortex–Undă Pilot În cadrul Modelului \Sigma-PDE, dualismul cuantic se separă pe cele două moduri cinematice ale ecuației Navier–Cauchy a continuumului: - Particula: Este un defect topologic solitonic chiral (vortex de forfecare transversală), având sarcina ca arie [e] = \text{m}^2 și masa m. Deplasarea sa fizică are loc pe o traiectorie unică, trecând printr-o singură fantă cu viteza subluminică v < c; - Unda Pilot: Este radiația acustică de compresiune pură \vec{u}_L (\vec{\nabla}\times\vec{u}_L = 0), emisă continuu de pulsația internă a vortexului. Deoarece se propagă cu celeritatea longitudinală: c_L = \sqrt{\frac{7}{3}}\,c = \mathbf{457\,940\,543,853956\dots\ \text{m/s}} aceasta se extinde în fața particulei, străbate ambele fante simultan și construiește în spațiul dintre fante și ecran o hartă de interferență staționară a dilatației de volum \theta(\vec{x}) = \vec{\nabla}\cdot\vec{u}_L. 1.2 Amplitudinea Madelung din Ecuația de Continuitate Ecuația de conservare a masei de fond în regim acustic liniar leagă perturbația densității locale de dilatație: \rho_{\text{mediu}}(\vec{x},t) = \rho_{\text{BI}}\left(1 - \theta(\vec{x},t)\right) În reprezentarea hidrodinamică Madelung, funcția de undă este scrisă polar: \Psi(\vec{x},t) = R(\vec{x},t)e^{iS(\vec{x},t)/\hbar_\sigma}. Modelul identifică amplitudinea de probabilitate R cu rădăcina pătrată a densității elastice: \boxed{R(\vec{x},t) \equiv R_0 \sqrt{1 - \theta(\vec{x},t)} = R_0 \sqrt{\frac{\rho_{\text{mediu}}(\vec{x},t)}{\rho_{\text{BI}}}}} 1.3 Potențialul Cuantic Madelung (Q_\Sigma) Potențialul cuantic de presiune internă este definit de curbura funcției de amplitudine: Q_\Sigma \equiv - \frac{\hbar_\sigma^2}{2m} \frac{\nabla^2 R}{R} unde \hbar_\sigma \equiv \hbar \Phi_{\text{geom}} = \frac{h}{2\pi}\frac{\pi}{6} = \frac{h}{12} = \mathbf{5,521725125 \times 10^{-35}\text{ J}\cdot\text{s}} este cuanta de acțiune a rețelei celulare Barbu. Dezvoltând în seria perturbațiilor slabe (\theta \ll 1): R = R_0(1 - \theta)^{1/2} \approx R_0\left(1 - \frac{1}{2}\theta\right) \implies \nabla^2 R \approx - \frac{1}{2}R_0 \nabla^2\theta \frac{\nabla^2 R}{R} \approx - \frac{1}{2}\nabla^2\theta Substituind în expresia lui Q_\Sigma: \boxed{Q_\Sigma \approx \frac{\hbar_\sigma^2}{4m}\nabla^2\theta} 1.4 Forța de Ghidaj Hidrostatică (\vec{F}_{\text{ghid}}) Asupra particulei acționează două forțe determinate mecanic: 1. Forța de depresiune de fond: Gradientul presiunii elastice pe volumul de forfecare al vortexului V_{\text{vortex}} = l_{\text{vortex}} \cdot e: \vec{F}_{\text{hidro}} = - l_{\text{vortex}} e P_{\text{BI}} \vec{\nabla}\theta 2. Forța potențialului cuantic: Forța derivată din potențialul Madelung: \vec{F}_{\text{cuantic}} = - \vec{\nabla}Q_\Sigma = - \frac{\hbar_\sigma^2}{4m}\vec{\nabla}(\nabla^2\theta) Forța totală de ghidaj exercitată determinist asupra p
Parmi les méthodes de la théorie du groupe de renormalisation, les développements en epsilon sont basés sur des calculs dans l'espace réciproque et permettent de calculer les exposants critiques associés aux transitions magnétiques du second ordre pour des valeurs non entières de la dimension d'espace. Une interprétation physique naturelle consiste à se demander comment se comporte un système de spins en interaction dans un espace de dimension fractale. Or les structures fractales sont construites par itération d'une cellule génératrice dont les détails sont donc présents à plusieurs échelles ; la question qui se pose alors est de savoir ce qui se passe lorsque l'invariance par translation est remplacée par une invariance d'échelle géométrique. La convergence vers la limite thermodynamique se produit en même temps que le processus d'itération construit la structure. De ce fait, des simulations Monte Carlo fiables de ces transitions de phase n'ont pu être menées à bien que récemment, puisqu'elles nécessitent la simulation de très grandes tailles, lesquelles varient comme des séries géométriques avec l'étape d'itération. C'est en utilisant des algorithmes non locaux dits “d'amas” (Wolff, Swendsen-Wang), capables de réduire le ralentissement critique de manière significative, et des méthodes d'histogrammes pour traiter les données des simulations Monte-Carlo que j'ai tout d'abord réalisé ces études. Il s'avère que le calcul précis des exposants critiques est rendu encore plus difficile par le fait que l'analyse en tailles finies du modèle d'Ising souffre de corrections d'échelle qui peuvent affecter fortement le comportement de certaines grandeurs thermodynamiques, en particulier lorsque la dimension fractale tend vers 1. J'ai montré que ces corrections d'échelle sont en partie liées à la très forte inhomogénéité du réseau sous-jacent (due à l'existence de trous sur plusieurs ordres de grandeurs) et à la concomitance de la construction du fractal avec la convergence vers la limite thermodynamique. Les résultats que j'ai obtenus pour les exposants critiques, ou leurs bornes, sont toujours compatibles avec la relation d'hyperscaling dans laquelle on substitue la dimension de Hausdorff à la dimension d'espace. Le comportement critique en dimension non entière se décrit dans le cadre de l'universalité faible. Cela se manifeste par un désaccord net entre les exposants que j'ai obtenus par les méthodes Monte Carlo et les développements en epsilon. Les exposants critiques ne dépendent pas seulement de la dimension d'espace, des propriétés de symétrie du paramètre d'ordre et de la portée des interactions, mais aussi des propriétés géométriques de la structure fractale : Très récemment des calculs précis d'exposants critiques m'ont permis de montrer que des classes d'universalité différentes sont en général nécessaires pour décrire le comportement du modèle d'Ising sur des fractals de même dimension et de lacunarités différentes. Un tel résultat généralise le concept d'universalité faible proposé par Masuo Suzuki. L'hypothèse d'homogénéité qui sous-tend les lois d'échelle permettant de décrire un comportement critique se dérive par renormalisation. La procédure de renormalisation dans l'espace direct est naturelle dans les fractals, puisqu'elle suit exactement le processus inverse de construction de la structure. Avec mon étudiant Pai-Yi Hsiao, nous avons mené à bien l'étude du modèle d'Ising par une méthode de renormalisation Monte-Carlo sur une structure fractale de dimension voisine de 1, 89 ; il s'avère que l'exposant associé à l'une des directions propres peut être calculé avec une très bonne précision et est en accord avec les résultats de l'analyse en tailles finies. En revanche, la convergence est très lente dans l'autre direction, ce qui est lié aux corrections d'échelle mises en évidence lors de cette analyse. La cinétique stochastique associée à la formation des amas construits par l'algorithme de Wolff sous tend la compréhension du phénomène de ralentissement critique. J'ai montré que les distributions des tailles des amas de Wolff ont une propriété d'homogénéité qui fait intervenir l'un des exposants associé à une des directions propres du processus de renormalisation. Par ailleurs, les distributions des tensions de surface des amas vérifient une propriété analogue dans laquelle intervient un nouvel exposant critique. L'étude des fonctions d'autocorrélation m'a permis de calculer précisément les exposants dynamiques de Wolff lorsque la température critique est connue, et d'éclaircir l'évolution du ralentissement critique avec la dimension et la connectivité. Dans le cas de systèmes invariants par translation, l'ordre de la transition ferromagnétique du modèle de Potts est lié au nombre d'états de spin ; le passage du premier au second ordre est attendu pour des dimensions non entières. Par ailleurs, la présence de désordre peut, dans certaines conditions, induire une transition du second ordre à partir d'un système qui en présente une du premier. L'étude du comportement critique du modèle de Potts sur des fractals est donc particulièrement intéressante, puisque deux des paramètres qui le déterminent (dimensionnalité et désordre structurel) sont liés. Avec mon étudiant Pai-Yi Hsiao, nous avons montré que la transition associée au modèle de Potts à trois états sur une structure fractale de dimension voisine de 1, 89 est du second ordre. Les difficultés attendues lorsqu'on augmente le nombre d'états de spins se font déjà nettement sentir : Les corrections d'échelle empêchent de calculer la température critique avec une très bonne précision. Nous n'avons donc pu donner que des bornes pour certains exposants ; nous avons cependant clairement mis en évidence la différence entre les classes d'universalité de Potts à 2 et 3 états. L'étude de la percolation en dimension non entière est liée à la fois à celle du modèle de Potts et aux algorithmes d'amas. Elle est basée sur l'étude des moments de la distribution de taille des amas, ce qui nécessite la localisation de pics en fonction de la probabilité d'occupation. J'ai pu montrer que les corrections d'échelle n'affectent pratiquement pas le comportement des pics avec la taille des structures, et proposé de les interpréter en termes de "seuil effectif".
Les effets de taille finie au-dessus de la dimension critique supérieure Dans cette thèse on étudie les effets de taille finie au-dessus de la dimension critique supérieure d_c. Les effets de taille finie y ont longtemps été incomplètement compris, en particulier vis-à-vis de leur dépendance en fonction des conditions aux limites. La violation de la relation d’échelle dite d’hyperscaling a été l’un des aspects les plus évidents des difficultés rencontrées. Le désaccord avec le scaling usuel est dû au caractère de variable non pertinente dangereuse du terme de self-interaction dans la théorie en ϕ^4. Celle-ci était considérée comme dangereuse pour la densité d’énergie libre et les fonctions thermodynamiques associées, mais pas dans le secteur des corrélations. Récemment, un schéma nouveau de scaling a été proposé dans lequel la longueur de corrélation joue un rôle central et est également affectée par la variable non pertinente dangereuse. Ce nouveau schéma, appelé QFSS, est basé sur le fait que la longueur de corrélation exhibe au lieu du scaling usuel ξ~L un comportement en puissance de la taille finie ξ~L^ϙ. Ce pseudo-exposant critique ϙ est lié à la dimension critique supérieure et à la variable dangereuse. Au-dessous de d_c, cet exposant prend la valeur ϙ=1, mais au-dessus, il vaut ϙ=d/d_c. Le schéma QFSS est parvenu à réconcilier les exposants de champs moyen et le Finite-Size-Scaling tel que dérivé du Groupe de Renormalisation pour les modèles avec interactions à courte portée au-dessus de d_c en conditions aux limites périodiques. Si ϙ est un exposant universel, la validité de la théorie doit toutefois s’étendre également aux conditions de bords libres. Des tests initiaux dans de telles conditions ont mis en évidence de nouvelles difficultés: alors que le QFSS est valable au point pseudo-critique auquel les grandeurs thermodynamiques telles que la susceptibilité manifestent un pic à taille finie, au point critique on a pensé que c’était le FSS standard qui prévalait avec les exposants de champ moyen et ξ~L. On montre dans ce travail qu’il en va différemment de la situation au point critique et qu’à la place ce sont les exposants gaussiens qui s’appliquent en l’absence de variable non pertinente dangereuse. Pour mettre en évidence ce résultat, nous avons mené des simulations de modèles avec interactions à longue portée, qui peuvent être à volonté étudiés au-dessus de leur dimension critique supérieure. Nous avons aussi développé une étude des modes de Fourier qui permet de fournir des exemples de quantités non affectées par la présence de la variable non pertinente dangereuse
On a compris ces dernières années que certaines mesures d'intrications sont un outil efficace pour la compréhension et la caractérisation de phases nouvelles et exotiques de la matière, en particulier lorsque les méthodes traditionnelles basées sur l'identification d'un paramètre d'ordre sont insuffisantes. Cette thèse porte sur l'étude de quelques systèmes quantiques à basse dimension où un telle approche s'avère fructueuse. Parmi ces mesures, l'entropie d'intrication, définie via une bipartition du système quantique, est probablement la plus populaire, surtout à une dimension. Celle-ci est habituellement très difficile à calculer en dimension supérieure, mais nous montrons ici que le calcul se simplifie drastiquement pour une classe particulière de fonctions d'ondes, nommées d'après Rokhsar et Kivelson. L'entropie d'intrication peut en effet s'exprimer comme une entropie de Shannon relative à la distribution de probabilité générée par les composantes de la fonction d'onde du fondamental d'un autre système quantique, cette fois-ci unidimensionnel. Cette réduction dimensionnelle nous permet d'étudier l'entropie aussi bien par des méthodes numériques (fermions libres, diagonalisations exactes, ...) qu'analytiques (théories conformes). Nous argumentons aussi que cette approche permet d'accéder facilement à certaines caractéristiques subtiles et universelles d'une fonction d'onde donnée en général.Une autre partie de cette thèse est consacrée aux trempes quantiques locales dans des systèmes critiques unidimensionnels. Nous insisterons particulièrement sur une quantité appelée écho de Loschmidt, qui est le recouvrement entre la fonction d'onde avant la trempe et la fonction d'onde à temps t après la trempe. En exploitant la commensurabilité du spectre de la théorie conforme, nous montrons que l'évolution temporelle doit être périodique, et peut même être souvent obtenue analytiquement. Inspiré par ces résultats, nous étudions aussi la contribution de fréquence nulle à l'écho de Loschmidt après la trempe. Celle-ci s'exprime comme un simple produit scalaire -- que nous nommons fidélité bipartie -- et est une quantité intéressante en elle-même. Malgré sa simplicité, son comportement se trouve être très similaire à celui de l'entropie d'intrication. Pour un système critique unidimensionnel en particulier, notre fidélité décroît algébriquement avec la taille du système, un comportement rappelant la célèbre catastrophe d'Anderson. L'exposant est universel et relié à la charge centrale de la théorie conforme sous-jacente.
\documentclass[11pt,a4paper,reqno]{article} \usepackage[utf8]{inputenc} \usepackage[romanian]{babel} \usepackage{amsmath,amssymb,amsfonts,amsthm} \usepackage{geometry} \usepackage{booktabs} \usepackage{array} \usepackage{hyperref} \usepackage{cite} \usepackage{microtype} \usepackage{bm} \usepackage{physics} \usepackage{siunitx} \usepackage{tikz} \usetikzlibrary{arrows.meta,positioning,shapes.geometric} \geometry{ a4paper, left=20mm, right=20mm, top=25mm, bottom=25mm } \hypersetup{ colorlinks=true, linkcolor=blue!80!black, citecolor=red!70!black, urlcolor=blue!60!black } % Structuri Matematice \newtheorem{theorem}{Teoremă}[section] \newtheorem{lemma}[theorem]{Lemă} \newtheorem{proposition}[theorem]{Propoziție} \newtheorem{corollary}[theorem]{Corolar} \theoremstyle{definition} \newtheorem{definition}[theorem]{Definiție} \newtheorem{postulate}{Postulat} \theoremstyle{remark} \newtheorem{remark}[theorem]{Observație} \title{\textbf{\LARGE OPERA MAGNA: INFORMATICA ȘI MECANICA CUANTICĂ Σ-PDE}\\[0.2cm] \textbf{\Large (LUCRAREA I)}\\[0.3cm] \textbf{\LARGE Teoria Limitei Cuantice de Viteză în Continuumul Σ-PDE:}\\[0.2cm] \textbf{\Large Depășirea Barierei Lieb--Robinson prin Modul Longitudinal cL=7/3c și Arhitectura Porților Logice Sub-Picosecundă}\\[0.3cm] \large De la Propagarea Acustică a Corelațiilor de Fază la Decuplarea Bell Sub-Retardată} \author{\textbf{Ilie Barbu} \\ \textit{Cercetător Independent, Pitești, Argeș, România} \\ \small Cadru Metrologic: CODATA 2022 / SI Reformat (BIPM 2019) / Lieb--Robinson Bounds (1972) \\ \small Precizie Computațională: Propagare analitică exactă / Aritmetică de precizie extinsă arbitrară \\ \small Înregistrare Canonică DOI: \href{https://doi.org/10.5281/zenodo.22726457}{10.5281/zenodo.22726457}} \date{Septembrie 2026} \begin{document} \maketitle \begin{abstract} Prezenta lucrare formulează fundamentarea fizico-matematică și arhitectura inginerească a calculului cuantic de mare viteză în cadrul Modelului Σ-PDE (Sigma -- Presiune -- Densitate -- Elastodinamică). În informatica cuantică modernă, viteza de procesare și generarea stărilor corelate la distanță între qubiți pe un cip sunt considerate limitate superior de viteza luminii c conform limitei Lieb--Robinson (τgate≥d/c). Lucrarea demonstrează că această limitare provine din presupunerea eronată că interacțiunile cuantice sunt mediate exclusiv de moduri electromagnetice transversale. În continuumul viscoelastic Σ-PDE, corelația de fază este transportată nemijlocit de unda acustică de compresiune pură a vidului la celeritatea longitudinală invariantă: \begin{equation} c_L \equiv \sqrt{\frac{7}{3}}\,c = 457\,940\,543{,}853956\dots\text{ m/s} \approx 1{,}527525\,c \end{equation} Principalele rezultate ale acestei monografii sunt: \begin{enumerate} \item \textbf{Depășirea Barierei Lieb--Robinson:} Se deduce matematic noua limită relativistă a porților logice cuantice la distanța d: τgate,min=cLd=73(cd)≈0,654654(cd) reducând latența fundamentală de comutare cu exact 34,5346%. Pentru o interconectare pe cip de d=1 cm, timpul minim scade de la 33,36 ps la 21,84 ps. \item \textbf{Dinamica Hidrodinamică de Ghidaj:} Se arată că forța de ghidaj a vortexului cuantic pe cip: Fghid=−lvortexePBI∇θ−4mℏσ2∇(∇2θ) este amorsată de unda cL, stabilind condiția de blocaj de fază ∂tθ=cLℏσm∣v∣ care elimină zgomotul de defazaj în circuitele cuantice. \item \textbf{Protocolul Decizional Bell Sub-Retardat:} Se formulează un test experimental riguros: la intervale de comutare a setărilor de detecție δt situate în fereastra sub-retardată: cLd<δt<cd corelațiile cuantice persistă (S=22), în timp ce la timpi ultra-rapizi δt<d/cL, unda pilot nu mai poate stabili faza, provocând colapsul parametrului Bell--CHSH la limita clasică S≤2,00000000. \end{enumerate} Lucrarea stabilește baza matematică pentru o nouă clasă de circuite logice cuantice sub-picosecundă, demonstrate fără parametri ajustabili a posteriori. \end{abstract} \newpage \section{Introducere și Limitele Actualei Informatici Cuantice} \subsection{Bariera Lieb--Robinson și Limita Cuantică de Viteză (QSL)} În teoria informației cuantice pe rețele discrete (procesoare supraconductoare, ioni captați, centri NV), propagarea semnalelor și a corelațiilor de fază este restricționată de celebra teoremă demonstrată de Elliott H. Lieb și Derek W. Robinson în 1972 \cite{Lieb1972}. Pentru orice operatori locali A^ și B^ separați de distanța d: \begin{equation} \left\| \left[ \hat{A}(t), \hat{B}(0) \right] \right\| \le C \|\hat{A}\| \|\hat{B}\| \exp\left[ -\mu (d - v_{\rm LR} t) \right] \end{equation} unde vLR este viteza Lieb--Robinson. În electrodinamica cuantică standard aplicată cipurilor, viteza de transfer a informației este asimilată cu viteza semnalului electromagnetic transversal: \begin{equation} v_{\rm LR} \le c = 299\,792\,458\text{ m/s} \end{equation} Aceasta impune o limită inferioară dură asupra timpului de execuție al unei porți logice cu 2 qubiți (\textit{two-qubit gate}): \begin{equation} \tau_{\rm gate} \ge \frac{d}{c} \end{equation} Pentru un cip de 1 cm, timpul minim fizic este de circa 33,36 picosecunde. Împreună cu limita Margolus--Levitin bazată pe energia medie Δt≥2ΔEπℏ \cite{Margolus1998}, această barieră dictează plafonul frecvenței de ceas a arhitecturilor cuantice. \subsection{Perspectiva Modelului Σ-PDE} Modelul Σ-PDE demonstrează că limita vLR≤c este doar o consecință a cuplării la modurile transversale electromagnetice (fotonice). În spatele acestui sector, vidul este un continuum viscoelastic real care suportă modul acustic de compresiune pură: \begin{equation} c_L \equiv \sqrt{\frac{7}{3}}\,c \approx 457\,940{,}544\text{ km/s} \approx 1{,}527525\,c \end{equation} Transferul corelațiilor de fază cuantică prin tensorul de Green acustic longitudinal permite realizarea porților cuantice în timpi strict mai mici decât timpul luminii: \begin{equation} \tau_{\rm gate, min} = \frac{d}{c_L} < \frac{d}{c} \end{equation} --- \section{Cadrul Axiomatic Σ-PDE și Parametrii de Rețea} Arhitectura computațională Σ-PDE este fundamentată mecanic prin Postulatul 0: \begin{postulate}[Postulatul 0: Reducerea Mecanică a Curentului] Unitatea SI de curent electric este fluxul mecanic de arie pe secundă: \begin{equation} 1\text{ Amper} \equiv 1\text{ m}^2/\text{s} \iff [I] = L^2 T^{-1} \end{equation} \end{postulate} \begin{corollary}[Sarcina ca Arie 2D] Sarcina electrică elementară e este aria geometrică a canalului de forfecare leptonic: \begin{equation} [e] = [I] \cdot [T] = L^2 = \text{m}^2, \quad e = 1{,}602176634 \times 10^{-19}\text{ m}^2 \end{equation} \end{corollary} \begin{lemma}[Parametrii Mecanici de Fond ai Vidului] Constantele de vid devin mărimile constitutive ale substratului de calcul: \begin{align} \rho_{\rm BI} &\equiv \mu_0 = 1{,}256637061262814754\dots \times 10^{-6}\text{ kg/m}^3 \\ P_{\rm BI} &\equiv \frac{1}{\varepsilon_0} = \mu_0 c^2 = 1{,}129409066602570333\dots \times 10^{11}\text{ Pa} \\ Z_{\rm BI} &\equiv \mu_0 c = 376{,}730313409875819\dots\text{ Pa}\cdot\text{s/m} \equiv \frac{\text{kg}}{\text{m}^2\cdot\text{s}} \\ c &= \sqrt{\frac{P_{\rm BI}}{\rho_{\rm BI}}} = 299\,792\,458\text{ m/s} \end{align} \end{lemma} Celula Wigner--Seitz a rețelei discrete are pasul: \begin{equation} l_{\rm BI} = 1{,}782174279826352072\dots \times 10^{-33}\text{ m} \end{equation} și timpul de comutare de punct zero: \begin{equation} \tau_{\rm BI} \equiv \frac{l_{\rm BI}}{c} = 5{,}944693511408\dots \times 10^{-42}\text{ s} \end{equation} --- \section{Deducerea Variațională a Celerității de Fază cL=7/3c} \subsection{Separarea Helmholtz în Ecuația Navier--Cauchy} Dinamica deplasării u(x,t) pe rețeaua elastică de vid este guvernată de ecuația Navier--Cauchy dedusă din acțiunea Σ-QG: \begin{equation} \rho_{\rm BI} \ddot{\mathbf{u}} = \mu \nabla^2\mathbf{u} + (\lambda + \mu)\nabla(\nabla \cdot \mathbf{u}) \end{equation} În virtutea identității de stare K≡μ≡PBI demonstrată în lucrările anterioare: \begin{equation} \mu = P_{\rm BI}, \quad \lambda = \frac{1}{3}P_{\rm BI} \implies \lambda + \mu = \frac{4}{3}P_{\rm BI}, \quad \lambda + 2\mu = \frac{7}{3}P_{\rm BI} \end{equation} Descompunem vectorul de deplasare conform teoremei Helmholtz: u=uT+uL, cu ∇⋅uT=0 și ∇×uL=0. \begin{theorem}[Viteza de Transfer a Fazei Cuantice] Modulul de compresiune irotațional uL se propagă fără dispersie la scări mezoscopice cu celeritatea longitudinală exactă: \begin{equation} \boxed{c_L = \sqrt{\frac{\lambda + 2\mu}{\rho_{\rm BI}}} = \sqrt{\frac{7 P_{\rm BI}}{3 \rho_{\rm BI}}} = \sqrt{\frac{7}{3}}\,c = 457\,940\,543{,}853956\dots\text{ m/s}} \end{equation} \end{theorem} \begin{proof} Pentru modul longitudinal pur, identitatea vectorială ∇×(∇×uL)=0⟹∇2uL=∇(∇⋅uL) reduce ecuația la: \begin{equation} \rho_{\rm BI} \ddot{\mathbf{u}}_L = \left( P_{\rm BI} + \frac{4}{3}P_{\rm BI} \right) \nabla^2\mathbf{u}_L = \frac{7}{3}P_{\rm BI} \nabla^2\mathb