We assign a Riemannian metric to a system of nonlinear equations that describe the one-dimensional propagation of long magnetoacoustic waves (also called magnetosonic waves) in a cold plasma under the inference of a transverse magnetic field. The metric, which in general is expressed in terms of the density of the plasma and its speed across the magnetic field, when specialized to a particular solution of the nonlinear system (the Gurevich-Krylov (G-K) solution) is mapped explicitly to a Jackiw-Teitelboim (J-T) black hole metric, which is the main result. Dilaton fields, constructed from data involved in the G-K solution, are presented - which with the plasma metric provide for elliptic function solutions of the J-T equations of motion in 2d dilaton gravity. A correspondence between solutions of the nonlinear plasma system (whose Galilean invariance is also established) and certain solutions of a resonant nonlinear Schrödinger equation is set up, along with some other general background material to render an expository tone in the presentation.
Citation: Floyd L. Williams. From a magnetoacoustic system to a J-T black hole: A little trip down memory lane[J]. Communications in Analysis and Mechanics, 2023, 15(3): 342-361. doi: 10.3934/cam.2023017
[1] | Yonghang Chang, Menglan Liao . Nonexistence of asymptotically free solutions for nonlinear Schrödinger system. Communications in Analysis and Mechanics, 2024, 16(2): 293-306. doi: 10.3934/cam.2024014 |
[2] | Vladimir Rovenski . Generalized Ricci solitons and Einstein metrics on weak $ K $-contact manifolds. Communications in Analysis and Mechanics, 2023, 15(2): 177-188. doi: 10.3934/cam.2023010 |
[3] | Efstratios Stratoglou, Alexandre Anahory Simoes, Leonardo J. Colombo . Reduction in optimal control with broken symmetry for collision and obstacle avoidance of multi-agent system on Lie groups. Communications in Analysis and Mechanics, 2023, 15(2): 1-23. doi: 10.3934/cam.2023001 |
[4] | Sergey A. Rashkovskiy . Nonlinear Schrödinger equation with a short-range compensating field. Communications in Analysis and Mechanics, 2025, 17(2): 520-549. doi: 10.3934/cam.2025021 |
[5] | Chunyou Sun, Junyan Tan . Attractors for a Navier–Stokes–Allen–Cahn system with unmatched densities. Communications in Analysis and Mechanics, 2025, 17(1): 237-262. doi: 10.3934/cam.2025010 |
[6] | Erlend Grong, Irina Markina . Harmonic maps into sub-Riemannian Lie groups. Communications in Analysis and Mechanics, 2023, 15(3): 515-532. doi: 10.3934/cam.2023025 |
[7] | Chen Yang, Chun-Lei Tang . Sign-changing solutions for the Schrödinger-Poisson system with concave-convex nonlinearities in $ \mathbb{R}^{3} $. Communications in Analysis and Mechanics, 2023, 15(4): 638-657. doi: 10.3934/cam.2023032 |
[8] | Tianfang Wang, Wen Zhang . Existence and concentration of homoclinic orbits for first order Hamiltonian systems. Communications in Analysis and Mechanics, 2024, 16(1): 121-146. doi: 10.3934/cam.2024006 |
[9] | Senhao Duan, Yue MA, Weidong Zhang . Conformal-type energy estimates on hyperboloids and the wave-Klein-Gordon model of self-gravitating massive fields. Communications in Analysis and Mechanics, 2023, 15(2): 111-131. doi: 10.3934/cam.2023007 |
[10] | Zhigang Wang . Serrin-type blowup Criterion for the degenerate compressible Navier-Stokes equations. Communications in Analysis and Mechanics, 2025, 17(1): 145-158. doi: 10.3934/cam.2025007 |
We assign a Riemannian metric to a system of nonlinear equations that describe the one-dimensional propagation of long magnetoacoustic waves (also called magnetosonic waves) in a cold plasma under the inference of a transverse magnetic field. The metric, which in general is expressed in terms of the density of the plasma and its speed across the magnetic field, when specialized to a particular solution of the nonlinear system (the Gurevich-Krylov (G-K) solution) is mapped explicitly to a Jackiw-Teitelboim (J-T) black hole metric, which is the main result. Dilaton fields, constructed from data involved in the G-K solution, are presented - which with the plasma metric provide for elliptic function solutions of the J-T equations of motion in 2d dilaton gravity. A correspondence between solutions of the nonlinear plasma system (whose Galilean invariance is also established) and certain solutions of a resonant nonlinear Schrödinger equation is set up, along with some other general background material to render an expository tone in the presentation.
We consider a direct connection of black holes in the Jackiw-Teitelboim model of two-dimensional dilaton gravity to the dynamics of two-component, cold collisionless plasma in the presence of an external transverse magnetic field. The propagation of long magnetoacoustic waves in the cold plasma (under a uni-axial propagation assumption) is described by the key nonlinear system of equations (2.5), where ρ is the plasma density, u is its speed across the magnetic field →B, and β>0 is a parameter that arises in a particular power series expansion of the magnetic field strength B(x,t) when a shallow water approximation is imposed. Such dynamics are discussed more generally in [1,2,3], for example, where a general set of equations is involved. These include the Maxwell equations ∇⋅→B=0 and ∂→B∂t=−curl →E (Gauss and Faraday laws), where the electric field →E is eliminated from the discussion by way of a few simplifying assumptions, including the above uni-axial propagation assumption - namely that →B=B(x,t)→ez and →u=u(x,t)→ex, for the velocity vector →u of the plasma. Later, u is also referred to as the "velocity" of the plasma. Getting back to the parameter β>0 in the second equation of (2.5) below, and the statement about the power series expansion of B, one can say that the latter expansion has the form
B=ρ+β2∂∂x(1ρ∂ρ∂x)+O(β4). | (1.1) |
The first equation in (2.5) is a continuity equation. The system of two equations in (2.5) is the reduction of an initial system of seven equations that are next reduced to a system of three equations, and then to the two equations by way of a shallow water approximation [1,4]. Actually, we show that the magnetoacoustic system (2.5), which is of fundamental interest here, can be derived directly and quickly from the equations of motion of a suitable Lagrangian density L, which is given in definition (2.1).
In [1], the system (2.5) is reduced to a resonant nonlinear Schrödinger equation (RNLSE) and particular solutions are generated that are descriptive of the interaction of solitonic magnetoacoustive waves. One can also obtain, conversely, solutions of the system (2.5) from solutions of the RNLSE, a matter that is discussed in section 3. The RNLSE is obtained from the nonlinear Schrödinger equation (NLSE) by the addition of a diffraction type term |ψ|xx/|ψ| called the de Broglie potential (see equation (3.1)), and it (like the NLSE) also admits both bright and dark soliton solutions. Equation (3.30) is an example of a dark soliton solution. The NLSE is known as an indispensable tool that facilitates the description of a multiplicity of nonlinear phenomena. The latter ranges, for example, from plasma physics and hydrodynamics to molecular biology (the nonlinear dynamics of DNA) and the propagation of light in optical fibers. There are applications to the modeling of extreme deep water rogue waves. The beautiful connection of the RNLSE to cold plasma physics is discussed, for example, in the papers [1,4,5,6,7], and in the new book [8]. The work in [6] was, and continues to be, the direct inspiration afforded by the author in [1,4,5,9].
Some other topics or ideas associated with the RNLSE - cold plasma physics connection are, for example, a reaction-diffusion system (RDS), a Madelung fluid system, and the continuous Heisenberg model. Of particular importance is the gauge equivalence (by way of the construction of suitable Lax pairs) of a particular RDS and a Heisenberg model. In the present paper, a direct connection of the cold plasma system (2.5) to a Jackiw-Teitelboim (J-T) black hole is emphasized. As is discussed in section 4, we use the continuous Heisenberg model to assign to the cold plasma system (2.5) a Riemannian metric gplasma. Then, in section 5, an explicit transformation of variables is presented by which the metric gplasma is mapped exactly to a J-T black hole metric gbh. gplasma and gbh are given in equations (4.19) and (5.1) respectively, where the black hole mass M in (5.1) is given by equation (5.7). m there is also given by equation (5.7), so that, in fact, the cosmological constant Λ in the J-T theory is given by Λ=−m2, which is negative. Section 5 contains the main result, where also some dilaton fields Φ(j)plasma, j=1,2,3, are presented for which the pairs (gplasma, Φ(j)plasma) provide for elliptic function solutions of the J-T gravitational field equations.
Of central interest for the discussions that follow is a system of two nonlinear partial differential equations that describe the propagation of one-dimensional long magneto-acoustic waves in a cold plasma of density ρ(x,t)>0 with a velocity u(x,t) across a magnetic field. See equations (2.5) below. The references [1,2,3], for example, provide background material and further details. Here, we are contented to start with two functions S(x,t) and ρ(x,t)>0, and to consider the Lagrangian density given by
L=L(x,t,S,Sx,St,ρ,ρx,ρt)def.=ρ(St−(Sx)2)−ρ24+β24(ρx)2ρ | (2.1) |
where β>0 is a fixed real number, and where the subscripts x,t,xx denote partial differentiation as usual. The corresponding equations of motion
∂∂t(∂L∂ft)+∂∂x(∂L∂fx)=∂L∂f, f=S,ρ, | (2.2) |
respectively, are directly computed:
ρt−2(ρSx)x=0,−St+(Sx)2+ρ2+β22[ρxxρ−12(ρxρ)2]=0, | (2.3) |
the first equation here being a continuity equation and the second one being a Hamilton-Jacobi equation. Next, we differentiate the second equation in (2.3) with respect to x, assuming that Sx=−u/2 for a third function u(x,t). Using the equality of mixed partial derivatives, we get Stx=Sxt=−ut/2, which with the equations
∂∂x(Sx)2=2SxSxx=(−u)(−ux/2) | (2.4) |
and (2.3) lead to the system
ρt+(ρu)x=0,ut+uux+ρx+β2[ρxxρ−12(ρxρ)2]x=0. | (2.5) |
This is the magnetoacoustic system (MAS) of interest for what follows.
We note that this system is Galilean invariant. That is, for c0∈R the field of real numbers, c0 fixed, define the Galilean transforms ˆu, ˆρ of u, ρ by
ˆu(x,t)=c0+u(x−c0t,t),ˆρ(x,t)=ρ(x−c0t,t). | (2.6) |
Then, by the chain rule,
(ˆρt+(ˆρˆu)x)(x,t)=(ρt+(ρu)x)(x−c0t,t),(ˆut+ˆuˆux+ˆρx)(x,t)=(ut+uux+ρx)(x−c0t,t),(ˆρxxˆρ−12(ˆρxˆρ)2)(x,t)=(ρxxρ−12(ρxρ)2)(x−c0t,t),[ˆρxxˆρ−12(ˆρxˆρ)2]x(x,t)=[ρxxρ−12(ρxρ)2]x(x,t). | (2.7) |
Thus if (u,ρ) is a solution of the MAS (2.5) then the pair (ˆu,ˆρ) is also a solution.
The standard Jacobi elliptic functions sn(x,κ), cn(x,κ), dn(x,κ) with elliptic modulus κ will be needed. Their basic properties are the following [8,10]:
sn2(x,κ)+cn2(x,κ)=1, dn2(x,κ)+κ2sn2(x,κ)=1,sn(x,1)=tanhx, cn(x,1)=dn(x,1)=sechx,ddxsn(x,κ)=cn(x,κ)dn(x,κ), ddxcn(x,κ)=−sn(x,κ)dn(x,κ)ddxdn(x,κ)=−κ2sn(x,κ)cn(x,κ). | (2.8) |
In particular, some attention will be given to the following traveling wave solution of the system (2.5), due to A.Gurevich and A.Krylov [11], expressed in terms of the elliptic function dn(x,κ). For u0∈R, u0>0 fixed, and for choices α1, α2, α3∈R with α3>α2⩾α1⩾0
ρ(x,t)=α1+(α3−α1)dn2((α3−α1)122β(x−u0t),κ)>0,u(x,t)=u0+Cρ(x,t), κdef.=+(α3−α2α3−α1)12, Cdef.=+(α1α2α3)12. | (2.9) |
A convenient way to express ρ(x,t) and C is
ρ(x,t)=α1+4a2β2dn2(a(x−βvt),κ),C=α121[4a2β2(1−κ2)+α1]12[4a2β2+α1]12 | (2.10) |
for
adef.=+(α3−α1)122β>0, vdef.=u0β>0, | (2.11) |
since
1−κ2def.=1−(α3−α2)(α3−α1)=α2−α1α3−α1,4a2β2def.=α3−α1, (4a2β2+α1)12=α123,4a2β2(1−κ2)=(α3−α1)(α2−α1α3−α1)=α2−α1⇒[4a2β2(1−κ2)+α1]12=α122. | (2.12) |
For plasma physics, convenient choices for α1 and α2 are the values α2=α1=1 since then κ=1, and by (2.8), (2.10) the solution
ρ(x,t)=1+4a2β2sech2a(x−βvt) | (2.13) |
for the plasma density does achieve the convenient value of 1 as |x|→∞. Also by (2.8), (2.10) note that
ρx(x,t)=−8a3β2κ2(sncndn)(a(x−βvt),κ),ρt(x,t)=8a3β3κ2v(sncndn)(a(x−βvt),κ)⇒ρt+u0ρxdef.=ρt+vβρx=0. | (2.14) |
For the choice α1=0, C=0 and u=u0 by (2.9). Also one can check that
[ρxxρ−12(ρxρ)2](x,t)=2√ρ(x,t)∂2∂x2(√ρ(x,t))for α1=0=22aβdn(a(x−βvt),κ)∂2∂x2(2aβdn(a(x−βvt),κ)=−2a2κ2[−sn2+cn2](a(x−βvt),κ)=a2[−2dn2(a(x−βvt),κ)+2−κ2]2 | (2.15) |
by (2.8), which with (2.10) gives (again for α1=0)
(ρ2+β22[ρxxρ−12(ρxρ)2])(x,t)=β2a2(2−κ2). | (2.16) |
As Sx=−u2=−u02=−vβ2
ρt−2(ρSx)x=ρt+u0ρx=0 | (2.17) |
(by (2.14)), which is the first equation in (2.3). By (2.16), the second equation in (2.3) is
−St+v2β24+β2a2(2−κ2)=0⇒S(x,t)=[v2β24+β2a2(2−κ2)]t+f(x), | (2.18) |
for some function of integration f(x).
−vβ2=Sx(x,t)∴=f′(x)⇒f(x)=−vβ2x+b,b∈R,⇒S(x,t)=[v2β24+β2a2(2−κ2)]t−vβx2+b. | (2.19) |
Thus, in summary, for the choice α1=0 in (2.9), so that C=0 and u(x,t)=u0=vβ (by (2.11)), S(x,t) given by (2.19), for b∈R, which satisfies Sx=−u/2, is a solution of the equations of motion given by the system (2.3) for the Gurevich-Krylov (G-K) solution ρ(x,t) given in (2.9), or in (2.10), of the MAS (2.5).
Given a triple of functions (S, ρ, u), we set up the Lagrangian density L in (2.1) for the pair (S,ρ) and we differentiated with respect to x the second corresponding equation of motion in (2.3) to derive the MAS of (2.5) for the pair (ρ, u) assuming that S was a velocity potential for u - namely that Sx=−u/2. Conversely, given a pair (ρ, u) that solves the system (2.5), there always exists a velocity potential S for u that solves the system (2.3). To see this, start with any velocity potential S0 of u whatsoever: S0x=−u/2. Then, again, S0tx=S0xt=−ut/2 and ((S0x)2)x=uux/2 (as in (2.4)) so that for
Fdef.=−2S0t+2(S0x)2+ρ+β2[ρxxρ−12(ρxρ)2], | (2.20) |
we get by (2.5) that
Fx=−2(−ut/2)+ρx+β2[ρxxρ−12(ρxρ)2]x=0, | (2.21) |
or that (F/2)x=0, which says that
F(x,t)/2=ϕ(t) | (2.22) |
for some function ϕ(t) of integration. Now choose any function h(t) such that h′(t)=ϕ(t). Then
S(x,t)def.=S0(x,t)+h(t)⇒Sx=S0xdef.=−u/2, St=S0t+h′def.=S0t+ϕdef.=S0t+F/2⇒−St+(Sx)2+ρ2+β22[ρxxρ−12(ρxρ)2]=−S0t−F2+(S0x)2+ρ2+β22[ρxxρ−12(ρxρ)2]=0 | (2.23) |
by the definition of F in (2.20). That is, S defined in (2.23) solves the second equation in (2.3), S is a velocity potential for u (Sx=−u/2 by (2.23)), and S solves the first equation in (2.3) since by (2.5), 0=ρt+(ρu)x=ρt+ρ(−2Sx)x. Thus the converse assertion regarding the existence of S is established. For the G-K solution in (2.10), we know that for α1=0, S is given by (2.19). Note that the first equation in (2.3) can be written as
0=ρt−2(ρSxx+ρxSx)2ρ=ρt2ρ−Sxx−ρxSxρ. | (2.24) |
The purpose of this section is to set up an explicit correspondence between solutions (u, ρ>0) of the magnetoacoustic system (MAS) in (2.5) and certain solutions ψ of the resonant nonlinear Schrödinger (RNLS) equation
iψt+ψxx+γ|ψ|2ψ=δ|ψ|xxψ|ψ|,γ,δ∈R, | (3.1) |
with de Broglie quantum potential |ψ|xx/|ψ|. The choice for δ will be δdef.=1+β2, for β in (2.5). Besides its occurrence in plasma physics, the RNLS equation occurs in quite many other studies as well. It occurs in the study of nonlinear fiber optics, for example. It's stability and dynamic properties, which are not discussed here, are considered in [12], for example. The Galilean invariance established for the MAS (2.5) is considered for equation (3.1).
Given a pair of real-valued functions (S, ρ>0) of (x, t) and c>0, c∈R, set
ψ=√cρ e−iS | (3.2) |
for i2=−1. Then
ψt=ψ(−iSt+ρt2ρ), ψx=ψ(−iSx+ρx2ρ),ψxx=ψ(12[ρxxρ−12(ρxρ)2]−(Sx)2+i[−Sxx−Sxρxρ])⇒iψt+ψxx=ψ(12[ρxxρ−12(ρxρ)2]−(Sx)2+St+i[−Sxx−Sxρxρ+ρt2ρ]). | (3.3) |
Also (as noted in (2.15))
2√ρ∂2∂x2(√ρ)=ρxxρ−12(ρxρ)2, |ψ|2=cρ⇒ρ2=|ψ|22c; 1√ρ∂2∂x2(√ρ)=√c|ψ|∂2∂x2(|ψ|√c)=|ψ|xx|ψ|⇒|ψ|xx|ψ|=12[ρxxρ−12(ρxρ)2]⇒iψt+ψxx+|ψ|2ψ−2c=ψ(|ψ|xx|ψ|+St−(Sx)2−ρ2+i[−Sxx−Sxρxρ+ρt2ρ]) | (3.4) |
by (3.3).
Now suppose (u, ρ>0) solves the MAS (2.5). Then from section 2 we know that S with Sx=−u/2 can be chosen to solve the system (2.3); see (2.23). The first equation in (2.3) is the equation
ρt2ρ−Sxx−ρxSxρ=0 | (3.5) |
by (2.24). Then by (3.4) and the second equation in (2.3)
iψt+ψxx+|ψ|2ψ−2c=ψ(|ψ|xx|ψ|+St−(Sx)2−ρ2)=ψ(|ψ|xx|ψ|+β22[ρxxρ−12(ρxρ)2])=ψ(|ψ|xx|ψ|+β2|ψ|xx|ψ|)=ψ(1+β2)|ψ|xx|ψ|, | (3.6) |
which is the RNLS equation (3.1) for δ=1+β2>1 and γ=−1/2c<0,
Conversely, suppose we are given a solution ψ of the RNLS equation
iψt+ψxx+|ψ|2ψ−2c=(1+β2)|ψ|xxψ|ψ|,c∈R, c>0, | (3.7) |
where ψ is of the form
ψ=eR−iS | (3.8) |
for real-valued functions R, S of (x, t). Define
u=−2Sx, ρ=e2Rc>0. | (3.9) |
Then we claim that the pair (u, ρ) solves the MAS (2.5). √cρ=eR⇒ψ=√cρe−iS, which means that the formulas in (3.4) apply. By (3.4) and (3.7)
ψ(|ψ|xx|ψ|+St−(Sx)2−ρ2+i[−Sxx−Sxρxρ+ρt2ρ])=iψt+ψxx+|ψ|2ψ−2c∴=ψ(1+β2)|ψ|xx|ψ|⇒|ψ|xx|ψ|+St−(Sx)2−ρ2+i[−Sxx−Sxρxρ+ρt2ρ]=(1+β2)|ψ|xx|ψ|=|ψ|xx|ψ|+β2|ψ|xx|ψ|⇒St−(Sx)2−ρ2=β2|ψ|xx|ψ|, −Sxx−Sxρxρ+ρt2ρ=0, | (3.10) |
where we equate real and imaginary parts. Again by (2.24), the last equation in (3.10) is the first equation in (2.3). Also by (3.4) again, the next to the last equation in (3.10) is the equation
St−(Sx)2−ρ2=β22[ρxxρ−12(ρx2)2]. | (3.11) |
In other words, the last two equations in (3.10) are precisely the two equations in the system (2.3) - which are the equations of motion for the Lagrangian density L in (2.1). Moreover, since Sx=−u/2 (by definition (3.9)), we have already shown that differentiation of equation (3.11) with respect to x leads exactly to the second equation in (2.5). Thus, conversely, (u, ρ) solves (2.5).
In summary, the following has been established. Given a solution(u, ρ>0) of the MAS (2.5), we can always choose a velocity potential S for u (ie. Sx=−u/2) such that the pair (S, ρ) solves the system (2.3). Given γ∈R, γ<0, set
ψdef.=√cρe−iS, c=1−2γ>0. | (3.12) |
Then ψ is a solution of the RNLS equation (3.1) for δdef.=1+β2>1. Conversely, suppose for real-valued functions R, S that
ψdef.=eR−iS | (3.13) |
is a solution of equation (3.1), with δ=1+β2, γ<0. Define
udef.=−2Sx,ρdef.=e2Rc>0,c=1−2γ>0. | (3.14) |
ψ=√cρ e−iS, the pair (S, ρ) solves the system (2.3) and the pair (u, ρ) solves the MAS (2.5).
For the Gurevich-Krylov (G-K) solution (2.10) with the choice α1=0, we know that S is given by (2.19). Therefore by (3.12) one obtains the solution
ψ(x,t)=2aβ√−2γdn(a(x−βvt),k)e−i(β2[v24+a2(2−κ2)]t−vβx2+b) | (3.15) |
of (3.1), again for δ=1+β2>1, b∈R, α1=0, b∈R. In particular for the choice κ=1, this solution reduces to the 1-soliton solution
ψ(x,t)=2aβ√−2γ sech(a(x−βvt))e−i(β2[v24+a2]t−vβx2+b) | (3.16) |
of (3.1), by (2.8).
The resonant nonlinear Schrödinger equation (3.1) is the equation of motion
∂∂t(∂L∂¯ψt)+∂∂x(∂L∂¯ψx)=∂L∂¯ψ | (3.17) |
for the Lagrangian density L given by
Ldef.=i2(¯ψψt−¯ψtψ)−¯ψxψx+δ(|ψ|x)2+γ2|ψ|4. | (3.18) |
The main point here is the key formula
∂∂x[∂∂ψx(|ψ|x)2]=|ψ|xx¯ψ|ψ|+∂∂ψ(|ψ|x)2. | (3.19) |
Namely, one has that
∂∂t(∂L∂¯ψt)=∂∂t(−iψ2)=−i2ψt, ∂L∂¯ψx=−ψx+δ∂∂¯ψx(|ψ|x)2⇒∂∂x(∂L∂¯ψx)=−ψxx+δ∂∂x[∂∂¯ψx(|ψ|x)2]=−ψxx+δ|ψ|xxψ|ψ|+δ∂∂¯ψ(|ψ|x)2, | (3.20) |
by conjugation of the key formula (3.19), Also
∂L∂¯ψ=i2ψt+δ∂∂¯ψ(|ψ|x)2+γ¯ψψ2 | (3.21) |
since γ2|ψ|4=γ2ˉψ2ψ2. Thus equation (3.17) is the equation
−i2ψt−ψxx+δ|ψ|xxψ|ψ|+δ∂∂¯ψ(|ψ|x)2=i2ψt+δ∂∂¯ψ(|ψ|x)2+γ¯ψψψ. | (3.22) |
That is,
δ|ψ|xxψ|ψ|=iψt+ψxx+γ|ψ|2ψ | (3.23) |
which is exactly the RNLS equation (3.1).
For c0∈R fixed, the Galilean transforms ˆu, ˆρ of u, ρ were defined in (2.6). Now Sx=−u/2, and we can define the Galilean transform ˆS of S by
ˆS(x,t)def.=−c02x+c204t+S(x−c0t,t), | (3.24) |
since then
ˆS(x,t)=−c02+Sx(x−c0t,t)=−c02−u2(x−c0t,t)=−ˆu(x,t)2, | (3.25) |
by (2.6). For ψ=eR−iS, with R, S= real functions, we saw (following (3.9)) that ψ=√ρe−iS for ρ=e−2R. This suggests that we set
ˆψ(x,t)=√ˆρ(x,t)e−iˆS(x,t); | (3.26) |
that is, by (2.6), (3.24),
ˆψ(x,t)=√ρ(x−c0t,t)e−i[−c0x/2+c20t/4+S(x−c0t,t)]=ei[c0x/2−c20t/4]ψ(x−c0t,t). | (3.27) |
The point is that indeed for any solution ψ of the RNLS equation (3.1), with ψ of the form eR−iS for real functions R, S, its Galilean transform ˆψ defined by the last equation in (3.27) will also be a solution of equation (3.1). Throughout, one may regard c0 as a velocity parameter.
As a simple example, for ρ0, γ∈R, ρ0>0
ψ0(x,t)def.=√ρ0eiγρ0t | (3.28) |
is a solution of (3.1) - a ground state on condensate solution. It's Galilean transform
ˆψ0(x,t)=√ρ0ei[c02x+(γρ0−c20/4)t], c0∈R | (3.29) |
therefore is also a solution.
We have connected equations (2.5) and (3.1) by the choice δ=1+β2>1. For δ<1, there are also important solutions of (3.1) of independent interest. For example, in [12] for δ<1 and μ<0 the dark (or topological) soliton solution
ψ(x,t)def.=eiμt√−μ tanh√−μ2(1−δ)x | (3.30) |
is considered, in case γ=−1 in (3.1)
The ultimate goal of our discussion is to establish a connection between the MAS (2.5) and a J-T black hole. For this, an initial key point is the assignment of a suitable Riemannian metric gplasma to the system (2.5). Fortunately, this can be done by way of the classical continuous (hyperbolic) Heisenberg model, to which we turn our attention. It is the genius and beauty of mathematics that often enough there exist startling, connective threads between seemingly disparate and unrelated topics or ideas. This happens in the present case here.
<, > will denote the Minkowski inner product on R3 given by
<X,Y>def.=−x1y1+x2y2−x3y3 | (4.1) |
for X=(x1,x2,x3), Y=(y1,y2,y3)∈R3. The Heisenberg model of interest is given by real-valued functions S1(x,t), S2(x,t), S3(x,t) for which
H:(x,t)→H(x,t)def.=(S1(x,t), S2(x,t), S3(x,t))∈R3 | (4.2) |
satisfies certain equations of motion(see (4.6)), and
<H(x,t), H(x,t)>=−1:−S21(x,t)+S22(x,t)−S23(x,t)=−1, | (4.3) |
by (4, 1). Thus the points H(x,t) lie on a single-sheeted hyperboloid. The function H:R2→R3 provides for a natural induced metric (or fundamental form) gH on the model:
gH=[g11g12g21g22]def.=[⟨Hx,Hx⟩⟨Hx,Ht⟩⟨Hx,Ht⟩⟨Ht,Ht⟩]. | (4.4) |
In other words, by (4.1), (4.2)
g11=−(∂S1∂x)2+(∂S2∂x)2−(∂S3∂x)2g21=g12=−∂S1∂x∂S1∂t+∂S2∂x∂S2∂t−∂S3∂x∂S3∂tg22=−(∂S1∂t)2+(∂S2∂t)2−(∂S3∂t)2. | (4.5) |
It was mentioned that certain equations of motion are satisfied by the model. They are
St=12i[S,Sxx]def.=12i(SSxx−SxxS)for Sdef.=i[S3S1−S2S1+S2−S3]. | (4.6) |
One can write
S1=S++S−2, S2=S+−S−2, for S+def.=S1+S2, S−def.=S1−S2. | (4.7) |
Then the metric gH and the equations of motion assume the form
g11=−∂S+∂x∂S−∂x−(∂S3∂x)2g21=g12=−∂S+∂x∂S−∂t−∂S+∂t∂S−∂x2−∂S3∂x∂S3∂tg22=−∂S+∂t∂S−∂t−(∂S3∂t)2,∂S3∂t=12(S−∂2S+∂x2−S+∂2S−∂x2), ∂S−∂t=S3∂2S−∂x2−S−∂2S3∂x2∂S+∂t=S+∂2S3∂x2−S3∂2S+∂x2. | (4.8) |
As an example [5], for a,v∈R with a≠0 and 4a2−v2≠0, define
ϕ(x,t)=(a2+v2/4)t−vx/2S+(x,t)def.=8a24a2−v2eϕ(x,t)(secha(x−vt))[tanha(x−vt)+v2a]S−(x,t)def.=8a24a2−v2e−ϕ(x,t)(secha(x−vt))[tanha(x−vt)−v2a]S3(x,t)def.=8a24a2−v2sech2a(x−vt)−1. | (4.9) |
Then
S3(x,t)2+S+(x,t)S−(x,t)=1. | (4.10) |
That is, (4.3) holds since S1, S2 satisfy (4.7). Also the equations of motion for S3, S+, S− in (4.8) hold, which are equivalent to the equation of motion for S1, S2, S3 in (4.6), as has been noted.
Of utmost practical importance is the fact that the metric gH in (4.4) has constant Ricci scalar curvature. It is given by
R(gH)=2, | (4.11) |
which can be verified by a Maple program(tensor), for example. In particular, for a positive multiple bgH of gH, b>0, b∈R, one has that
R(bgH)=2b. | (4.12) |
We move now to the main point of this section that connects the MAS (2.5) and the continuous Heisenberg model - a connection explicated by way of the metric gH in (4.4). Fix γ<0, γ∈R, γ as in (3.1) for example. Let (u, ρ>0) be a solution of the system (2.5). Then it is possible to construct a solution (r,s) of the reaction-diffusion system
rt−rxx+Br2s=0, st+sxx−Brs2=0, Bdef.=−γβ2. | (4.13) |
Namely, for S0(x,t) with S0x=−u/2 (as usual)
r(x,t)def.=[ρ(x,t/β)/(−2γ)]12eϕ0(x,t)s(x,t)def.=−[ρ(x,t/β)/(−2γ)]12e−ϕ0(x,t)ϕ0(x,t)def.=S0(x,t/β)/β. | (4.14) |
From the solution (r,s) in (4.14), one can construct, moreover, a solution H(x,t) in (4.3) of the equations of motion (4.6) (or (4.8)) that define the continuous Heisenberg model. The construction of H(x,t) is a bit involved. It amounts to the construction of a suitable Lax pair to establish in fact that the reaction-diffusion system and Heisenberg model are gauge equivalent. Details are presented in [5,8], for example. In the end, from a solution (u,ρ>0) of the MAS (2.5) one obtains a metric gH (see (4.4)) such that the multiple
gplasmadef.=gH(−2γ/β2) | (4.15) |
of gH is given explicitly by the following formulas:
(gplasma)11(x,t)=ρ(x,t/β)/2γ(gplasma)21(x,t)=(gplasma)12(x,t)=ρ(x,t/β)u(x,t/β)−4γβ(gplasma)22(x,t)=ρ(x,t/β)−8γ[(1ρ∂ρ∂x)2−u2β2](x,t/β) | (4.16) |
Also by (4.12), gplasma has constant Ricci scalar curvature given by
R(gplasma)=−4γβ2>0. | (4.17) |
In particular, the formulas in (4.16) apply to the Gurevich-Krylov solution in (2.10), where (gplasma)22 is a little messy to compute. In the end, one obtains the following formulas, where details appear in section 15 of [8]:
(gplasma)11(x,t)=α1+4a2β2dn2(a(x−vt),κ)2γ(gplasma)21(x,t)=(gplasma)12(x,t)=−vγa2β2dn2(a(x−vt),κ)−vα14γ−C4γβ(gplasma)22(x,t)=4a2β2[a2κ4−2γ(sncn)2(a(x−vt),κ)+v28γdn2(a(x−vt),κ)]+16α1a4β2κ4(sncn)2(a(x−vt),κ)+C2β28γ[α1+4a2β2dn2(a(x−vt),κ)]+v2α18γ+vC4γβ. | (4.18) |
Again C is given by (2.9), or by (2.10).
Finally, although (apriori) gplasma is non-diagonal (ie.(gplasma)12≠0) there is a suitable change of variables (x,t)→(τ,δ) [6,13] by which gplasma assumes the diagonal form
gplasma:ds2=A(δ)dτ2−4a4β4κ4(sncndn)2(δ,κ)A(δ)γ2dδ2, | (4.19) |
where
A(δ)=4a2β2[a2κ4−2γ(sncn)2(δ,k)+v28γdn2(δ,κ)]+16α1a4β2κ4(sncn)2(δ,k)+C2β28γ[α1+4a2β2dn2(δ,κ)]+v2α18γ−vC4γβ, | (4.20) |
with conditions of course that A(δ) never vanishes. This is the case if v2 is sufficiently large. More precisely, as is shown in [6], if α1≠0, then A(δ) never vanishes if
v2>4a2κ4 and v2⩾4C2α21β2. | (4.21) |
If α1=0, then the single condition
v2>4a2κ4 | (4.22) |
suffices for the non-vanishing of A(δ). We point out that in [6], the notation A(ρ) is used for the A(δ) here. The g11 in [6,7,13] corresponds to the notation g22 used here and vice versa. In [13], α1=0, γ=−1/2, and b2 there is 4β2 here.
As in section 2, the various results assume a much simpler form for the choice α1=0, mainly as then C=0 by definition (2.9). For example, (4.20) reduces to
A(δ)=4a2β2[a2κ4−2γ(sncn)2(δ,κ)+v28γdn2(δ,κ)]. | (4.23) |
Moreover, for the choice κ=1 and for γ of the form γ=Λ/4 for a suitable cosmological constant Λ, A(δ) in (4.23) assumes the form
A(δ)=−8a2β2Λ(sech2δ)[a2tanh2δ−v24], | (4.24) |
(by (2.8)), and consequently (4.19) reduces to
gplasma:ds2=−8a2β2Λ(sech2δ)[a2tanh2δ−v24]dτ2+8a2β2tanh2δsech2δΛ[a2tanh2δ−v24]dδ2=−8a2β2Λ(sech2δ)[(a2tanh2δ−v24)dτ2−tanh2δ(a2tanh2δ−v24)dδ2], | (4.25) |
which, up to the factor β2, is exactly the black hole metric ds2 in equation (3.14) of [9], with a horizon singularity at
a2tanh2δ−v24=0: tanhδ=±v2a. | (4.26) |
Note that by formula (4.24), A(δ)≠0 indeed for v2>4a2, as asserted in (4.22).
Given a solution ψ of equation (3.1) of the form (3.13) (with δ=1+β2, γ<0 in (3.1)), we saw that the prescription
u=−2Sx, ρ=e2Rc, c=1−2γ>0 | (4.27) |
in (3.14) provided for solutions (S,ρ>0), (u,ρ>0) of the systems (2.3), (2.5), respectively. By the equations in (4.16) we can express the components gij of the plasma metric g=gplasma in terms of ψ. We proceed as follows. Since |ψ|2=cρ, the 1st formula in (4.16) gives
g11(x;t)=−cρ(x,t|β)=−|ψ(x,t/β)|2. | (4.28) |
Next by (3.3) and (4.27)
ψx=ψ(−iSx+ρx2ρ)⇒(ˉψ)x=ˉψ(iSx+ρx2ρ)⇒ψx(ˉψ)x=ψˉψ[(Sx)2+14(ρxρ)2]=cρ[u24+14(ρxρ)2]=cρ4[u2+(ρxρ)2], ψ(ˉψ)x−ˉψψx=|ψ|2(2iSx)=−icρu. | (4.29) |
By the last equation here and the 2nd equation in (4.16)
g21(x,t)=g12(x,t)=(ρu)(x,t|β)−4γβ=c2β(ρu)(x,t/β)=i2β[ψ(ˉψ)x−ˉψψx](x,t/β). | (4.30) |
Finally, again as |ψ|2=cρ, we see that cρx=2|ψ||ψ|x⇒c2(ρx)2=4cρ(|ψ|x)2⇒
4c(|ψ|x)2=(ρx)2ρ, | (4.31) |
and by the next to last equation in (4.29)
(1+1β2)4c(∣ψ|x)2−4cβ2ψx(ˉψ)x=(1+1β2)(ρx)2ρ−4cβ2cρ4[u2+(ρxρ)2]=(ρx)2ρ+1β2(ρx)2ρ−ρu2β2−(ρx)2β2ρ=(ρx)2ρ−ρu2β2=ρ[(ρxρ)2−u2β2]⇒c4ρ[(ρxρ)2−u2β2]=(1+1β2)(|ψ|x)2−1β2ψx(ˉψ)x, | (4.32) |
where c/4=1/−8γ. In other words, by the last equation in (4.16)
g22(x,t)=[(1+1β2)(|ψ|x)2−1β2ψx(ˉψ)x](x,t|β). | (4.33) |
Formulas (4.28), (4.30) and (4.33) provide for the expression of the cold plasma metric components gij=(gplasma)ij in (4.16) in terms of the solution ψ in (3.8) of the RNLS equation (3.1), for δ=1+β2, γ<0. Here for (u,ρ>0) in (4.16), we have the MAS ↔ RNLS equation correspondence (u,ρ>0)↔ψ discussed in section 3.
Formula (4.25) provides for a realization of the cold plasma metric gplasma as a black hole metric in case of the special choices α1=0 and κ=1 in (2.9). For α1⩾0 arbitrary and for an arbitrary elliptic modulus κ, where gplasma assumes the general form given in (4.19), in the variables τ, δ with A(δ) given by (4.20), we can in fact map gplasma, more specifically, to a Jackiw-Teitelboim (J-T) black hole metric gbh given by
gbh=−(m2r2−M)dτ2+dr2(m2r2−M) | (5.1) |
in the variables τ, r. Here an explicit transformation of variables (τ,δ)→(τ,r) is presented by which gplasma in (4.19) is indeed mapped to gbh in (5.1), which is the main result of this section, where the black hole parameters m,M are expressed in terms of the magnetoacoustic parameters in the G-K solution (2.9) (or (2.10)).
We begin first with some brief, contextual remarks regarding the J-T model of 2d gravity. Any metric g whatsoever on a two-dimensional space-time M2 automatically solves the Einstein gravitational vacuum field equations with a zero matter tensor [14]. Given this well known fact, R. Jackiw and C. Teitlboim set out to construct a non-trivial theory of gravity for M2 that involved in addition to g a scalar field Φ on M2 [15,16]. Φ is called a dilaton field. The action integral for the theory is given, up to some constant, by
S(g,Φ)=∫M2(R(g)−2l2)Φ√|detg|d2x, | (5.2) |
where (as in section 4) R(g) denotes the Ricci scalar curvature of g, and where Λ=−1/l2 is a (negative) cosmological constant in the theory. The corresponding equations of motion for the pair (g,Φ) are given by
R(g)=2l2=−2Λ,∇i∇jΦ=gijΦl2,1≤i,j≤2, | (5.3) |
where for local coordinates (x1,x2) on M2, and for the Christoffel symbols (of the second kind) Γkij for g, the Hessian in (5.3) is given by
∇i∇jΦ=∂2Φ∂xi∂xj−2∑k=1Γkij∂Φ∂xk. | (5.4) |
Thus g has constant Ricci scalar curative, as does gplasma in (4.17). We note that here, and throughout, our sign convention for scalar curvature is the negative of that employed by J-T in [15,16], and possibly by others in the literature. A key solution of the theory, for example, is of course the J-T black hole solution gbh given in the Lorentzian form (5.1) with coordinates (x1,x2)=(τ,r), where m=1/l:
R(g)=2m2,Φ(τ,r)def.=mr,Λ=−m2, | (5.5) |
with M being a black hole mass parameter.
Returning to the remarks that followed equation (5.1), we now state the main result. Namely, the transformation of variables (τ,δ)→(τ,r) by which the cold plasma metric gplasma in (4.19) is mapped to the J-T black hole metric in (5.1) is given by
r=Ψ(δ)def.=a2β2dn2(δ,κ)−γ+a2β2(2−κ2)2γ−(v2β2−α1)8γ, | (5.6) |
where m and M in (5.1) are given by
m=+(−2γ)1/2βM=β2−2γ[vC2β3+A−α1(3v28β2+3α116β4+a2(2−κ2)2β2)],for Adef.=−a2v22(2−κ2)+a4κ4+v416. | (5.7) |
We keep the assumptions in (4.21) in order to have A(δ)≠0∀δ. Again see (2.9), (2.10), (2.11) for the notation in (5.7). Also see (4.15) where we chose γ∈R, γ<0 arbitrary. The result (5.6), (5.7) are proved in [6], but the more compact version of the black hole mass M in (5.7) is noted in [7]. In [6] it is shown that indeed M>0 for
v2>8a2(2−κ2)+6α1β2,vC4β−3α2132β2−α1a22⩾0. | (5.8) |
For α1=0, C=0 (again by definition (2.9)) and so the second inequality in (5.8) is automatic. For α1>0 it is the statement
v⩾α1C[3α18β+2a2β],α1≠0. | (5.9) |
For α1=0, (5.7) and (5.8) reduce to
M=β2−2γA,v2>8a2(2−κ2). | (5.10) |
Now 2κ2+κ4⩽2+1<4(as κ⩽1)⇒κ4<4−2κ2=2(2−κ2)⇒4a2κ4<8a2(2−κ2)⇒v2−4a2κ4>v2−8a2(2−κ2). That is, for α1=0 the inequality in (5.10) that implies that M>0 also implies that v2>4a2κ4, which is the condition in (4.22) for the non-vanishing of A(δ) in (4.23).
Since the choice α1=0 is so practical and important as well, we provide a direct argument that A>0 if v2>8a2(2−κ2) (as in (5.8) or (5.10)), and hence also M>0 by (5.10), again as γ<0. Note that by definition (5.7) we can write
A=−a2v2+a2κ2v22+a4κ4+v416=2γβ2[−a2β2v22γ+a2κ2v2β24γ+a4κ4β22γ+v4β232γ]. | (5.11) |
Now if v2>8a2(2−κ2) then for γ<0
v4β232γ=(v2β232γ)v2<(v2β232γ)8a2(2−κ2)=(v2β232γ)(16a2−8a2κ2)=v2β2a22γ−v2β2a2κ24γ⇒v4β232γ−v2β2a22γ+v2β2a2κ24γ<0. | (5.12) |
That is, for the bracket in (5.11), the sum of the 1st, 2nd and 4th terms is <0. The 3rd term a4κ4β2/2γ in the bracket is also <0 (since γ<0)⇒ the bracket <0, which multiplied by the negative number 2γ/β2 gives that A>0 if v2>8a2(2−κ2), as claimed, in which case M>0 also in (5.10).
By (5.5) and (5.7) the cosmological constant Λ is given by Λ=2γ/β2<0. Also by (5.5) and (5.6), for the dilaton field Φ(1)plasma given by
Φ(1)plasma(τ,δ)def.=mΨ(δ)=m[a2β2dn2(δ1κ)−γ+a2β2(2−κ2)2γ−(v2β2−α1)8γ], | (5.13) |
we obtain a solution (gplasma, Φ(1)plasma) of the J-T field equations (5.3):
R(gplasma)+2Λ=0,∇i∇jΦ(1)plasma+ΛgijΦ(1)plasma=0, | (5.14) |
where the first equation here is equation (4.17), and the second set of equations are argued for in [6]. There are two more plasma dilaton fields Φ(j)plasma, j=2,3, for gplasma in (4, 19) that are computed in [7] (not in [6]), and also in [13] by a different method. For the choice α1=0, equation (5.13) can be written as
Φ(1)plasma(τ,δ)=mβ2−2γ[2a2dn2(δ,κ)−a2(2−κ2)+v24], | (5.15) |
and the other two dilaton fields are given by
Φ(2)plasma(τ,δ)=√2−γaβdn(δ,κ)[v24−a2κ4(sncndn)2(δ,k)]12sinh(√Aτ)Φ(3)plasma(τ,δ)=√2−γaβdn(δ,κ)[v24−a2κ4(sncndn)2(δ,κ)]12cosh(√Aτ), | (5.16) |
for A in (5.7), where we checked directly that A>0, given the standing assumption that v2>8a2(2−κ2) in (5.8). Also, one has the inequality
(sncndn)2(x,κ)⩽1 | (5.17) |
so that in (5.16) (where already γ<0) we always have that
[v24−a2κ4(sncndn)2(δ,κ)]⩾v24−a2κ4>0. | (5.18) |
Here as we have seen, in the argument that followed (5.10), the assumption v2>8a2(2−κ2)⇒v2/4−a2κ4>0. Thus the scalar fields in (5.16) are real, as they should be, and they with the field in (5.15) and the cold plasma metric provide for new solutions in Jackiw-Teitelboim dilaton gravity. An additional reference for section 4 and for this section is [17].
We have maintained a special interest in the magnetoacoustic system (2.5) that describes the propagation of one-dimensional magnetoacoustic waves in a cold plasma of density ρ with a speed v across a transverse magnetic field. Using a positive multiple of the Riemannian metric gH in (4.4), (4.5) (also see (4.8)) attached to the classical, continuous (hyperbolic) Heisenberg model, we assigned to the nonlinear system (2.5) the metric gplasma defined in (4.15), where β>0 is the parameter in the second equation of (2.5), which arises from the magnetic field strength expansion (1.1), and where γ is an arbitrary negative real number. Given the correspondence that was set up in section 3 between solutions of the system (2.5) and certain solutions of the resonant nonlinear Schrödinger equation in (3.1), a choice sometimes for γ is the coefficient of |ψ|2ψ in (3.1). gplasma being a real number multiple of gH also has constant Ricci scalar curvature given by R(gplasma)=−4γ/β2 in (4.17), and its components (gplasma)ij, 1⩽i, j⩽2, are given in terms of the cold plasma density and speed (ρ and u) by the concrete formulas in (4.18). If (ρ,u) is the traveling wave solution of the system (2.5) given by Gurevich and Krylov in (2.9) (or in (2.10)), then the plasma metric components are given by the concrete formulas in (4.18), As was shown in [6,13], a change of variables (from (x,t) in (4.18) to new variables (τ,δ)) exists by which the plasma metric assumes the diagonal form (4.19), with A(δ) defined in (4.20). The main result is that by another change of variables (τ,δ)→(τ,r), where r is given explicitly in terms of δ by (5.6), the cold plasma metric gplasma diagonalized in (4.19) is mapped exactly to the Jackiw-Teitelboim (J-T) black hole metric gbh given in (5.1), where m and M in (5.1) are given by the formulas in (5.7). M is the black hole mass, and −m2 is the negative J-T cosmological constant Λ.
By way of the Gurevich-Krylov solution and the continuous Heisenberg model, we have established a direct connection of the magnetoacoustic system (2.5) for a cold plasma to a black hole solution gbh in the J-T theory of 2d dilaton gravity. Dilaton fields Φ(j)plasma, 1⩽j⩽3, in terms of elliptic functions are presented in (5.15), (5.16) for which the pairs (gplasma, Φ(j)plasma) (again for gplasma in (4.19)) are solutions of the J-T gravitational field equations (5.3).
The author is most grateful to Yaping Yuan for the assistance she has provided in the excellent preparation of this manuscript.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
The author declares there is no conflict of interest.
[1] |
J. H. Lee, O. K. Pashaev, C. Rogers, W. K. Schief, The resonant nonlinear Schrödinger equation in cold plasma physics. Applications of Bäcklund-Darboux transformations and superposition principles, J. Plasma Physics, 73 (2007), 257–272. https://doi.org/10.1017/S0022377806004648 doi: 10.1017/S0022377806004648
![]() |
[2] | V. I. Karpman, Nonlinear waves in dispersive media, Pergamon, Oxford, 1975. https://doi.org/10.1016/B978-0-08-017720-5.50008-7 |
[3] | A. I. Akhiezer, Plasma electrodynamics, Pergamon, Oxford, 1975. |
[4] |
J. H. Lee, O. K. Pashaev, Solitons of the resonant nonlinear Schrödinger equation with nontrivial boundary conditions and Hirota bilinear method, Theoretical and Mathematical Physics, 152 (2007), 991–1003. https://doi.org/10.1007/s11232-007-0083-3 doi: 10.1007/s11232-007-0083-3
![]() |
[5] |
O. K. Pashaev, J. H. Lee, Resonance solitons as black holes in Madelung fluid, Modern Physics Letters A, 17(2002), 1601–1619. https://doi.org/10.1142/S0217732302007995 doi: 10.1142/S0217732302007995
![]() |
[6] | F. Williams, Exploring a cold plasma - 2d black hole connection, Advances in Mathematical Physics, Hindawi. |
[7] |
J. D'Ambroise, F. Williams, Relating some nonlinear systems to a cold plasma magnetoacoustic system, J. Modern physics, 11 (2020), 886–906. https://doi.org/10.4236/jmp.2020.116054 doi: 10.4236/jmp.2020.116054
![]() |
[8] | F. Williams, Some musings on theta, eta, and zeta: From E8 to cold plasma to an inhomogeneous universe, accepted for publication in Springer Book series: Mathematical Physics Studies, Springer, 2023. |
[9] | L. Martina, O. K. Pashaev, G. Soliani, Bright solitons as black holes, Physical Rev. D, 58 (1998). https://doi.org/10.1103/PhysRevD.58.084025 |
[10] | K. Chandrasekharan, Elliptic functions, Grundlehren der mathematischen Wissenschaften, Springer-Verlag, Berlin, Heidelberg, New York, Tokyo, 1985. https://doi.org/10.1007/978-3-642-52244-4 |
[11] | A. V. Gurevich, A. L. Krylov, A shock wave in dispersive hydrodynamics, Soviet Physics Doklady, 32 (1988), 73–74. |
[12] | T. P. Horikis, P. G. Kevrekidis, F. Tsitoura, F. Williams, Solitary waves in the resonant nonlinear Schrödinger equation: Stability and dynamical properties, Physics Letters A, 2020. http://doi.org/10.1016/j-physleta.2020.126441 |
[13] | J. D'Ambroise, F. Williams, Elliptic function solutions in Jackiw-Teitelboim dilaton gravity, Hindawi. https://doi.org/10.1155/2017/2154784 |
[14] | F. Williams, Some selected thoughts old and new on soliton-black hole connections in 2d dilaton gravity, in: The sine-Gordon Model and its Applications: From Pendula and Josephson Junctions to Gravity and High Energy Physics, J.C.-Maraver, P.G.Kevrekidis, F. Williams, Editors, Springer (2014), 171–205. https://doi.org/10.1007/978-3-319-06722-3-8 |
[15] | R. Jackiw, A two dimensional model for gravity, in: Quantum Theory of Gravity, S.Christensen, Editor, Adam Hilger Ltd. (1984), 403–420. |
[16] | C. Teitelboim, The Hamiltonian structure of two-dimensional spacetime and its relation with the conformal anomaly, in: Quantum Theory of Gravity, S.Christensen, Editor, Adam Hilger Ltd. (1984), 327–344. |
[17] | L. Martina, O. K. Pashaev, G. Soliani, Integrable dissipative structures in the gauge theory of gravity, Classical and Quantum Gravity, 14 (1977). https://doi.org/10.1088/0264-9381/14/12/005 |