Alice-Bob Physics: Coherent Solutions of Nonlocal KdV Systems

In natural and social science, many events happened at different space-times may be closely correlated. Two events, A (Alice) and B (Bob) are defined correlated if one event is determined by another, say, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{B}}=\hat{{\boldsymbol{f}}}A$$\end{document}B=fˆA for suitable \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hat{{\boldsymbol{f}}}$$\end{document}fˆ operators. Taking KdV and coupled KdV systems as examples, we can find some types of models (AB-KdV systems) to exhibit the existence on the correlated solutions linked with two events. The idea of this report is valid not only for physical problems related to KdV systems but also for problems described by arbitrary continuous or discrete models. The parity and time reversal symmetries are extended to shifted parity and delayed time reversal. The new symmetries are found to be useful not only to establish AB-systems but also to find group invariant solutions of numerous AB-systems. A new elegant form of the N-soliton solutions of the KdV equation and then the AB-KdV systems is obtained. A concrete AB-KdV system derived from the nonlinear inviscid dissipative and barotropic vorticity equation in a β-plane channel is applied to the two correlated monople blocking events which is responsible for the snow disaster in the winter of 2007/2008 happened in Southern China.

It is known that various events happened in different places and times are closely correlated and/or entangled 1 . For example, the decrease of ice at Arctic sea in the summer of 2007 (event A, A = A(x, t)) is considered to be responsible for the heavy freezing rain in winter 2008 in Southern China (event B, = ′ ′ =B B x t fA ( , ) ) 2, 3 . The strong El Nino event occurred in 1997 (event A) induced the heavy Yangtze River flooding in 1998 (event ′ ′ =B x t f A ( , ) ) 4 . In the communication field, in order to safely get information from Alice (event A), Bob should have a private Key f such that he can find the thing he wants ( =B f A). It is reported that the event B, the detection of a signal of gravitational wave, is due to the event A, the merging of two black holes from 1.3 billion light years away 5 ; In quantum physics, many faraway particles (atoms) may construct a completely entangled state, a measurement of one particle (event A) will affect the state of the other (event =B f A) 6 .
Traditionally, physical models are locally established around a single space-time point, say, {x, t}. To describe the correlated events happened in two or more places, one has to establish some kinds of new physically significant models. In this paper, we restrict ourselves to study only two-place physical models, Alice-Bob (AB) models for simplicity. We say Alice-Bob physics if the physics is related to two correlated/entangled events occurred in two places {x, t} and {x′, t′} respectively. The event at {x, t} is called Alice event (AE) (denoted by A(x, t)) and the event at {x′, t′} is called Bob event (BE) (denoted by B(x′, t′)). The events AE and BE are called correlated/entangled if AE happens, the BE can be determined at once by the correlation condition f for suitable operator(s) f . Usually, {x′, t′} is not neighbour to {x, t}. Thus, the intrinsic two-place models, the Alice-Bob systems (ABSs), are nonlocal. Some special types of two-place nonlocal models have been proposed, say, the nonlocal nonlinear Schrödinger (NLS) equation, , t x x 2 was firstly proposed by Ablowitz and Musslimani 7 . This type of NLS systems are also related to the significant PT symmetric Schrödinger equations 8 . The operators P and Ĉ are the usual parity and charge conjugation. Some types of other nonlocal nonlinear systems such as the coupled nonlocal NLS systems 9 , the nonlocal modified KdV systems 10,11 , the discrete nonlocal NLS system 12 and the nonlocal Davey-Stewartson systems [13][14][15] were also proposed. Though the nonlocal NLS and the other known models mentioned above have been studied by many authors, they are all proposed in the mathematical way from the integrability requirement without any possible physical applications.
In this paper, we attempt to provide a possible physical application: They are special AB systems and can be possibly used to describe some special types of AB physics. In the second section, we discuss the general aspect of the AB systems via shifted parity and the delayed time reversal. In the third section, the detailed group invariant solutions for general AB-KdV systems and symmetry breaking solutions for a special AB-KdV system are studied. The multiple soliton solutions of the usual integrable systems have been applied in almost all the physical fields and there are some types of formally different complicated expressions such as the Hirota forms, determinant forms, Pffafion forms and so on 16 . Here, a much more elegant form of the solutions for the KdV type systems is presented. In the application section, a special AB-KdV system is derived from a rotating fluid system via multiple scale method. A symmetry breaking soliton solution of the special KdV system is applied to a two-atmospheric-blocking-events which is related to the heavy snow disaster in the winter of 2007/2008 in the south of China. The last section is a summary and discussion.

General aspect of AB systems
There exist two possible ways to establish ABSs which may be used to describe AB physics in the future.
(i) The first way is to solve out two events AE and BE from traditional basic principle equation(s), for suitable U, where u may be a scalar, a vector, a matrix or a tensor. Equation (2)  To find out some concrete results, we introduce a principle 17 , the AB-BA equivalence principle, which requires U(A, B) = U(B, A) and Substituting the two event assumption (3) into the original principle (2) and splitting the resulting equation to two, we may have Two equations of (5) should be reduced back to one equation because of the correlated condition (4). To reduce two equations of (5) to one we can apply the operator f on one of them. The results should be , the superscription T of a matrix means the transposition of the matrix, we have the following theorem and conjecture.

Theorem. The linear transformation
0 0 possesses the solution x, x 0 , M}→{Ux′, Ux, Ux 0 , UMU −1 } is also a solution of (7), the conjecture implies that the solution (8) is similarity transformation invariant up to redefinition of the arbitrary matrices Γ, Ω and x 0n .

AB-KdV systems and their coherent structures
To make it more clear, we take the well known KdV equation 18 as a simple example. The KdV equation was first introduced by Boussinesq 19 and rediscovered by Diederik Korteweg and Gustav de Vries 20 . The KdV equation has various connections to physical problems. It approximately describes the evolution of long, one-dimensional waves in many physical settings, including shallow-water waves with weakly non-linear restoring forces, long internal waves in a density-stratified ocean, ion acoustic waves in a plasma, acoustic waves on a crystal lattice 18 and the 2-dimensional quantum gravity 21 . The KdV equation can be solved using the inverse scattering transform and other methods such as those applied to other integrable systems 22 .
Furthermore, for simplicity, we take the symmetric U as a simplest form = + U A B ( ) 1 2 . From the equation (9), it is clear that the KdV equation is parity-time reversal invariant, space-time translation invariant and complex conjugate invariant. Thus, we can take f as for the real KdV system and =ˆˆf PT C s d , into the KdV equation (9), we have which can be split to two equations t xxx x where G(A, B) may be an arbitrary functional of A and B. Applying f defined in Eq. (10) on Eq. (12) will lead to a compatibility condition In other words, G(A, B) is not an arbitrary functional but an arbitrary f invariant functional. Finally, we obtain a quite general nonlocal AB-KdV system t xxx . It is known that the original KdV equation is a famous integrable model. However, the integrability of the AB-KdV system (14) is not very clear because the introduction of the arbitrary f invariant functional G(A, B) though there are some different formal Lax pairs. The first formal Lax pair possesses the form, Scientific RepoRts | 7: 869 | DOI:10.1038/s41598-017-00844-y and Φ is just the recursive operator of the classical KdV equation with the potential A − B.
If the arbitrary f invariant functional G is fixed, one may obtain more concrete Lax pairs, say, if we select x the related AB-KdV systems are just ABKdV i , i = 2, 3, 4 listed in ref. 17. The AB-KdV systems (14) with the conditions (18) can be considered as special reductions of the coupled KdV system derived from two-layer fluid systems 23 .
In addition to the AB-KdV system (14), there are other types of AB-KdV systems which can not be directly solved out from the KdV equation. For instance, the following AB-KdV system can be obtained from the reduction of coupled KdV systems, say, the Hirota-Satsuma system 24 . The Lax pair of the AB-KdV system (19) can be directly read out from that of the Hirota- 17,25 .
The most general real AB-KdV systems may have the more general form where KdV is defined in Eq. (9). Now a natural question is how to solve the AB-KdV systems listed in this paper. To solve complicated nonlinear systems, the symmetry methods play essential roles. In fact, the powerful Darboux and Bäcklund transformations are essentially special symmetry approaches 26,27 . Usually, one uses continuous symmetries to find group invariant solutions of nonlinear systems. Here we have derived many AB-systems by using discrete symmetries, the shifted parity and the delayed time reversal. Thus it is interesting that if we can find some nontrivial shifted parity and delayed time reversal invariant solutions for all the AB-KdV systems (20). The group invariant condi- (20) is equivalent to find the ˆP T s d invariant solutions of the usual KdV equation (9). For nonlinear systems, soliton excitations are most important exact solutions. For the KdV equation (9), its well known multiple soliton solutions possesses the form 28 where the summation of μ should be done for all permutations of µ = 0, 1 The reason is that the KdV equation is space-time translation invariant and then every soliton of the KdV equation can be located at anywhere ξ 0j . However, the AB-KdV is space-time translation symmetry breaking. Thus we have to fix the arbitrary constants ξ 0j such that the multiple soliton solution is ˆP T s d invariant but not the space-time translation symmetry invariant. After finishing some detailed calculations, we find that if we rewrite ξ j as 0 oj solves all the AB-KdV system (20) including Eqs (14) and (19). Similarly, all the ˆP T s d invariant solutions of the KdV equation are also solutions of all AB-KdV systems. Here we list two more examples, the Painlevé II reduction and the soliton-cnoidal periodic wave interaction solutions.
The ˆP T s d invariant Painlevé II reduction possesses the form, where U = U(ξ) satisfies which is equivalent to the Painlevé II equation.
A simple ˆP T s d invariant soliton-cnoidal periodic wave interaction solution [29][30][31] can be written as 0 0 with arbitrary constants v 1 , v 2 and m. It should be also emphasized that though we have obtained many f invariant solutions for many models, there are various other solutions which are f symmetry breaking. Usually, f symmetry breaking solutions may not be same for different models. Here we just write down an f symmetry breaking soliton solution for the AB-KdV equation (19), where k and a are arbitrary constants. More about f symmetry breaking solutions will be reported in our near future studies.

Derivation and application of a special AB-KdV system in atmospheric dynamics
In ref.

32, Jia et al. established a multiple vortex interaction model
iy jx , ≠ j i denotes the i-jth stream-vorticity interactions (SVI) and ω ω [ , ] i j denotes the i-jth vorticity-vorticity interactions (VVI), C is related to the strength of VVI and the β term comes from the coriolis force. In (30), the small parameter  is introduced by considering the fact that the stream-vorticity interaction between two faraway (both in space and time) events should be small.
For N = 1, the model (30) is just the well known (2 + 1)-dimensional rotating fluid model (Euler equation with rotating effect), such as the atmospheric and oceanic systems, which can be used to describe the nonlinear inviscid dissipative and equivalent barotropic vorticity equation (NIDEBE) in a β-plane channel 33 . The model (30) is also derived from NIDEBE by neglecting higher order smaller interactions among different vortices.
In fact, the model can also be used to describe N-event problems. For the N = 2 case, Eq. (30) becomes [ , ] [ , ] 0 (33)  To get some approximate analytic solutions of the AB-EBVE, we utilize the multiple scale method (MSM) to derive an AB-KdV system from Eqs (31)-(33) with conditions (34).
As in the standard MSM, by introducing some slow variables, while the model parameters, C and β, can also be expanded in some series of . In this paper, we just take 2   which means for two far away events, the vorticity-vorticity interaction is in the same order as for the stream-vorticity interaction while the effect of the coriolis force is smaller in the next order. Substituting the expansion (36) with Eqs (34) and (37) into the original model (32) and (33), we have yyy y yy yyy y y y y y y yy yy y y y yyy yyy yy    To solve A 2 from (45), we can take the following variable separation form     possesses the following periodic wave solution The atmospheric and oceanic phenomena are rich and multifarious. The richness and multifariousness of our approximate solution comes from the introduction of arbitrary functions … u u u , , , 0 1 7 . By selecting the arbitrary functions appropriately, we can explain various two-place correlated events in atmospheric and oceanic dynamics. For example, Fig. 1a displays the theoretic result of A-event (atmospheric blocking) described by ψ 1 of (51) happened near 40 degree east longitude while the Fig. 1b is the corresponding real weather chart of geopotential height at 500 hPa on 19th November, 2007, from the National Centers for Environmental Prediction/the National Center for Atmospheric Research (NCEP/NCAR) reanalysis data 34 . Figure 2a shows the theoretic density and contour plot of the B-event described by ψ 2 of (52) appeared near 130 degree west longitude while the Fig. 2b is the corresponding plot of the real weather chart happened after 65 days later after the A-event. The function and parameter selections related Figs 1a and 2a are as follows,

Summary and Discussions
In summary, a new physical fields, AB-Physics or even AB-Science, is opened. This is a first attempt to study two-place physics and multi-place physics, which may be widely and deeply affect all other scientific fields. The reason to introduce Alice-Bob physics is because there exist various two-place (and multi-place) correlated physical events in almost all natural scientific fields. For instance, all the teleconnection phenomena in atmospheric and oceanic science belong to the category of the AB-physics/science.
A simple method is established to find huge amount of new physically meaningful nonlinear models, Alice-Bob systems. The method can be applied to any known physics to find novel phenomena though here we take only the celebrating KdV model as an illustration to find infinitely many Alice-Bob KdV systems. A special type of AB-KdV systems, the functional KdV equations, with an arbitrary ˆP T s d invariant functional G (A, B), can be derived directly from the original KdV system. Some other types of AB-KdV equations can be directly obtained from coupled KdV systems by using correlated relations between two components. The method can be applied to any original principle physical models such as the Maxwell equations, Navier-Stokes equations, Schrödinger equations, Einstein field equation and their integrable and nonintegrable derivative systems. In the preprint paper 17 , various other types of integrable AB systems have been listed by directly applying suitable correlated relations on coupled two component models.
A special AB-KdV system is derived from the multiple vorticity interaction model which is related to a standard atmospheric and oceanic dynamic system, the nonlinear inviscid dissipative and barotropic vorticity equation in a β-plane channel. The ˆP T s d symmetry breaking soliton solution of the derived AB-KdV system is used to qualitatively describe the two real events, the atmospheric blocking happened in November 2007 and January 2008 respectively while the atmospheric blockings are responsible for the heavy snow disaster in Southern China in the winter 2007/2008. Some types of physically important new symmetries such as shifted parity and delayed time reversal, are found. It is well known that parity and time reversal are two very important symmetries in physics. Now these two symmetries are extended to more general symmetries. These new symmetries exist in various physical fields  Fig. 1, but for B-event (a) described by ψ 2 of (52), and real case (b) of geopotential height at 500 hPa on 23rd January, 2008 from NCEP/NCAR reanalysis data. The contour interval is 4 gpdm in (b), which is generated by using the charge-free software GrADS 2.0 download at http://grads.iges.org/grads. which include not only the fields where the KdV equation is valid but also those where the Maxwell equations, (linear and nonlinear) Schrödinger equations, Navier-Stokes equation and Einstein equation are valid. In addition to the shifted parity and delayed time reversal symmetries, we have found more newly general symmetries given in theorem 1. These symmetries will be very useful in two place physics and should be studied further.
Some types of group (shifted parity and delayed time reversal symmetry group) invariant solutions including N-soliton solutions (arbitrary N), Painlevé II reductions and soliton cnoidal wave interactions are obtained for all AB-KdV systems. The physical meaning of the group invariant solutions is that the event happened at {x, t} will happen also at {x′, t′}.
One special shifted parity (P s ) and delayed time reversal (T d ) symmetry breaking soliton solution is also given for special AB-KdV systems (14) and (48). The physical meaning of the group symmetry breaking solution is that for Alice-Bob systems there are real physical phenomena where the event B at {x′, t′} is different from the event A at {x, t}.
Infinitely many nonlinear excitations are found to be solutions of infinitely many models. This fact indicates and emphasizes that to conclude a theory should be very careful even if you have observed various (may be infinitely many) facts are valid for a theory without any parameters. For instance, even if you have observed N-soliton solutions for an arbitrary N in a real system, there are still infinitely many candidate theories to describe this system.
It is indicated that additional nonlinearities can be introduced by shifted parity and delayed time reversal correlated constrained conditions. That means shifted parity and delayed time reversal correlation conditions are nonlinear constraints. This implies that we can introduce useful nonlinearities from useful nonlocal nonlinear symmetry invariant constraints. In our future studies, we will prove that various physically important nonlinear systems (such as KdV, KP, NLS, sine-Gordon etc.) can be derived from linear systems by means of this types of nonlinear and nonlocal constraints.
It is also indicated that there are infinitely many nonlocal physical functional models which have not yet been studied by present physics and lack of effective methods to solve nonlocal and nonlinear functional physical models. This fact arouses many challenges to the present physical society.
Though the multiple soliton solutions for many integrable systems have been obtained in some types of formally different complicated expressions, in this paper, a much more elegant form for the KdV type systems is proposed. The similar elegant forms for many other types of soliton systems can also be obtained 17 by using the same approach.
In addition to the beauty of the new soliton expressions, our forms display two new physical phenomena. The first one is that the resonant soliton solution is analytical when the wave numbers k i , and k j are very closed to each other, however, it is not analytical for the traditional multiple soliton expressions when resonance happens. The second one is that if we take the positions of the solitons are wave numbers dependent, then we may obtain various types of quite different resonant solutions which need to be further explored in future.