Initial value problem of the Whitham equations for the Camassa-Holm equation

We study the Whitham equations for the Camassa-Holm equation. The equations are neither strictly hyperbolic nor genuinely nonlinear. We are interested in the initial value problem of the Whitham equations. When the initial values are given by a step function, the Whitham solution is self-similar. When the initial values are given by a smooth function, the Whitham solution exists within a cusp in the x-t plane. On the boundary of the cusp, the Whitham solution matches the Burgers solution, which exists outside the cusp.     

The initial value problem for the Whitham averaged system

We study the initial value problem for the Whitham averaged system which is important in determining the KdV zero dispersion limit. We use the hodograph method to show that, for a generic non-trivial monotone initial data, the Whitham averaged system has a solution within a region in thex-t plane for all time bigger than a large time. Furthermore, the Whitham solution matches the Burgers solution on the boundaries of the region. For hump-like initial data, the hodograph method is modified to solve the non-monotone (inx) solutions of the Whitham averaged system. In this way, we show that, for a hump-like initial data, the Whitham averaged system has a solution within a cusp for a short time after the increasing and decreasing parts of the initial data beging to interact. On the cusp, the Whitham and Burgers solutions are matched.     

Solution of the Bethe-Goldstone equation for the reaction matrix in finite nuclei

A method for solving the BG equation for the reaction matrixt in finite nuclei is presented. The application of this method is demonstrated for a one-dimensional case, which is similar to the problem where the internucleon potential acts only in the relatives-state. The single particle potential has a harmonic oscillator form and the phenomenological internucleon potentialv(r) contains a hard core and an attractive part of the Yukawa type. By taking the exclusion principle into account exactly an infinite system of integral equations is obtained. It is proved that the solution of the corresponding finite system converges to the exact solution. An iteration method for solving such a finite system with an arbitrary number of equations is developed. Its main feature consists in the exclusion of the dependence on the hard core part ofv(r) (which is treated as the limit case of a rectangular repulsive potential with a variable heightv ₀). This exclusion transforms the original system to a system of integral equations depending only on the attractive part ofv(r) and to a linear algebraic system. Both these systems can be solved by iteration for all values ofv ₀ as well as for v₀= +. The numerical results confirm the rapid convergence of the proposed iteration method and demonstrate that the solution of the finite system with a sufficiently large number of equations approximates the exact solution very precisely.     

Symplectic and hyperkähler structures in a dimensional reduction of the Seiberg-Witten equations with a Higgs field

In this paper we show that the dimensionally reduced Seiberg-Witten equations lead to a Higgs field and we study the resulting moduli spaces. The moduli space arising out of a subset of the equations, shown to be non-empty for a compact Riemann surface of genus g ≥ 1, gives rise to a family of moduli spaces carrying a hyperkähler structure. For the full set of equations the corresponding moduli space does not have the aforementioned hyperkähler structure but has a natural symplectic structure. For the case of the torus, g = 1, we show that the full set of equations has a solution, different from the “vortex solutions”.     

Theory of Gravitational-Inertial Field of Universe. I. Gravitational-Inertial Field Equations

In the series of present articles the original proposition is a generalization of the real world tensor by the introduction of a inertial field tensor. From this generalization it follows, particularly, that ?_ig_lm ? g_lm;i ≠ 0. This allows to use as a Lagrangian density of the field the expression A_g = k₁ g_lm;ig^lm _;kg^ik. On the basis of variational equations a system of more general covariant equations of the gravitational-inertial field is obtained. In the Einstein approximation these equations reduce to the field equations of Einstein. The solution of fundamental problems in the general theory of relativity by means of the new equations gives the same results as the solution by means of Einstein's equations. However, application of these equations to the cosmologic problem gives a result different from that obtained by Friedmann's theory. In particular, the solution gives the Hubble law as the law of motion of a free body in the inertial field - in contrast to Galileo-Newton's law.     

Physics-Based Preconditioning and the Newton–Krylov Method for Non-equilibrium Radiation Diffusion

An algorithm is presented for the solution of the time dependent reaction-diffusion systems which arise in non-equilibrium radiation diffusion applications. This system of nonlinear equations is solved by coupling three numerical methods, Jacobian-free Newton–Krylov, operator splitting, and multigrid linear solvers. An inexact Newton's method is used to solve the system of nonlinear equations. Since building the Jacobian matrix for problems of interest can be challenging, we employ a Jacobian–free implementation of Newton's method, where the action of the Jacobian matrix on a vector is approximated by a first order Taylor series expansion. Preconditioned generalized minimal residual (PGMRES) is the Krylov method used to solve the linear systems that come from the iterations of Newton's method. The preconditioner in this solution method is constructed using a physics-based divide and conquer approach, often referred to as operator splitting. This solution procedure inverts the scalar elliptic systems that make up the preconditioner using simple multigrid methods. The preconditioner also addresses the strong coupling between equations with local 2×2 block solves. The intra-cell coupling is applied after the inter-cell coupling has already been addressed by the elliptic solves. Results are presented using this solution procedure that demonstrate its efficiency while incurring minimal memory requirements.     

Random Matrices, Graphical Enumeration and the Continuum Limit of Toda Lattices

In this paper we derive analytic characterizations for and explicit evaluations of the coefficients of the matrix integral genus expansion. The expansion itself arises from the large N asymptotic expansion of the logarithm of the partition function of N × N Hermitian random matrices. Its g ^th coefficient is a generating function for graphical enumeration on Riemann surfaces of genus g. The case that we particularly consider is for an underlying measure that differs from the Gaussian weight by a single monomial term of degree 2ν. Our results are based on a hierarchy of recursively solvable differential equations, derived through a novel continuum limit, whose solutions are the coefficients we want to characterize. These equations are interesting in their own right in that their form is related to partitions of 2g + 1 and joint probability distributions for conditioned random walks. K. D. T-R McLaughlin was supported in part by NSF grants DMS-0451495 and DMS-0200749, as well as a NATO Collaborative Linkage Grant “Orthogonal Polynomials: Theory, Applications, and Generalizations” Ref no. PST.CLG.979738. N. M. Ercolani and V. U. Pierce were supported in part by NSF grants DMS-0073087 and DMS-0412310.     

Bipartite Matching and Van der Waerden Conjecture

In this paper I show that the free energy F and the cost C associated to a bipartite matching problem can be explicitly estimated in term of the solution of a suitable system of equations (cavity equations in the following). The proof of these results relies on a well known result in combinatorics: the Van der Waerden conjecture (Egorychev–Falikman Theorem). Cavity equations, derived by a mean field argument by Mèzard and Parisi, can be considered as a smoothed form of the dual formulation for the bipartite matching problem. Moreover cavity equation are the Euler–Lagrange equations of a convex functional G parameterized by the temperature T. In term of their unique solution it is possible to define a free-energy-like function of the temperature g(T). g is a strictly decreasing concave function of T and C=g(0). The convexity of G allows to define an explicit algorithm to find the solution of the cavity equations at a given temperature T. Moreover, once the solution of the cavity equations at a given temperature T is known, the properties of g allow to find exact estimates from below and from above of the cost C.     

Yang—Mills instantons over Riemann surfaces

Exact solutions to the self-dual Yang—Mills equations over Riemann surfaces of arbitrary genus are constructed. They are characterized by the conformal class of the Riemann surface. They correspond to U(1) instantonic solutions for an Abelian-Higgs system. A functional action of a genus g Riemann surface is constructed, with minimal points being the two-dimensional self-dual connections. The exact solutions may be interpreted as connecting initial and final nontrivial vacuum states of a conformal theory, in the sense of Segal, with a Feynman functor constructed from the functional integral of the action.     

Coisotropic Deformations of Associative Algebras and Dispersionless Integrable Hierarchies

The paper is an inquiry of the algebraic foundations of the theory of dispersionless integrable hierarchies, like the dispersionless KP and modified KP hierarchies and the universal Whitham hierarchy of genus zero. It stands out for the idea of interpreting these hierarchies as equations of coisotropic deformations for the structure constants of certain associative algebras. It discusses the link between the structure constants and Hirota’s tau function, and shows that the dispersionless Hirota bilinear equations are, within this approach, a way of writing the associativity conditions for the structure constants in terms of the tau function. It also suggests a simple interpretation of the algebro-geometric construction of the universal Whitham equations of genus zero due to Krichever.     

Theory of Gravitational-Inertial Field of Universe

In this paper the basic proposition is a generalization of the metric tensor by introduction of an inertial field tensor satisfying ?_ig_lm ? g_lm;i ≠ 0. On the basis of variational equations a system of more general covariant equations of gravitational-inertial field is obtained. In Einstein's approximation these equations reduce to the field equations of Einstein. The solution of fundamental problems of generl taheory of relativity by means of the new equations give the same results as Einstein's equations. However application of these equations to the cosmologic problem leads to following results: 1. All Galaxies in the Universe (actually all bodies if gravitational attraction is not considered) “disperse” from each other according to Hubble's law. Thus contrary to Friedmann's theory (according to which the “expansion of Universe” began from the singular state with an infinite velocity) the velocity of “dispersion” of bodies begins from the zero value and in the limit tends to the velocity of light. 2. The “dispertion” of bodies represents a free motion in the inertial field and Hubble's law represents a law of motion of free bodies in the inertial field - the law of inertia. All critical systems (with Schwarzschild radius) are specific because they exist in maximal inertial and gravitational potentials. The Universe represents a critical system, it exists under the Schwarzschild radius. In the high-potential inertial and gravitational fields the material mass in a static state or in the process of motion with decelleration is subject to an inertial and gravitational “annihilation”. Under the maximal value of inertial and gravitational potentials (= c²) the material mass is completely “evaporated” transforming into a radiation mass. The latter is concentrated in the “horizon” of the critical system. All critical systems –“black holes”- represent geon systems, i.e., the local formations of gravitational-electromagnetic radiations, held together by their own gravitational and inertial fields. The Universe, being a critical system, is “wrapped” in a geon crown. The Universe is in a state of dynamical equilibrium. Near the external part of its boundary surface a transformation of matter into electromagnetic-gravitational-neutrineal energy (geon mass) takes place. Inside the Universe, in the galaxies takes place the synthesis of matter from geon mass, penetrating from the external part of the world (from geon crown) by means of a tunneling mechanism. The geon system may be considered as a natural entire cybernetic system.     

A unified grand canonical description of the nonextensive thermostatistics of the quantum gases: Fractal and fractional approach

In this paper, the particles of quantum gases, that is, bosons and fermions are regarded as g-ons which obey fractional exclusion statistics. With this point of departure the thermostatistical relations concerning the Bose and Fermi systems are unified under the g-on formulation where a fractal approach is adopted. The fractal inspired entropy, the partition function, distribution function, the thermodynamics potential and the total number of g-ons have been found for a grand canonical g-on system. It is shown that from the g-on formulation; by a suitable choice of the parameters of the nonextensivity q, the parameter of the fractional exclusion statistics g, nonextensive Tsallis as well as extensive (q=1) standard thermostatistical relations of the Bose and Fermi systems are recovered. Received 17 September 1999     

The Local Gromov–Witten Theory of $${\mathbb{C}\mathbb{P}^1}$$ and Integrable Hierarchies

In this paper we begin the study of the relationship between the local Gromov–Witten theory of Calabi–Yau rank two bundles over the projective line and the theory of integrable hierarchies. We first of all construct explicitly, in a large number of cases, the Hamiltonian dispersionless hierarchies that govern the full-descendent genus zero theory. Our main tool is the application of Dubrovin’s formalism, based on associativity equations, to the known results on the genus zero theory from local mirror symmetry and localization. The hierarchies we find are apparently new, with the exception of the resolved conifold in the equivariantly Calabi–Yau case. For this example the relevant dispersionless system turns out to be related to the long-wave limit of the Ablowitz–Ladik lattice. This identification provides us with a complete procedure to reconstruct the dispersive hierarchy which should conjecturally be related to the higher genus theory of the resolved conifold. We give a complete proof of this conjecture for genus g ≤ 1; our methods are based on establishing, analogously to the case of KdV, a “quasi-triviality” property for the Ablowitz–Ladik hierarchy at the leading order of the dispersive expansion. We furthermore provide compelling evidence in favour of the resolved conifold/Ablowitz–Ladik correspondence at higher genus by testing it successfully in the primary sector for g = 2.     

Density of the SO(3) TQFT Representation of Mapping Class Groups

We show that for an odd prime r>3 and an integer g>1, in the projective representation given by the SO(3) Witten-Reshitikhin-Turaev theory at an r^th root of unity, the image of the mapping class group of a surface of genus g is dense. Partially supported by NSF DMS 0100537 and DMS 0354772. Partially supported by NSF EIA 0130388, and DMS 0354772, and ARO.     

Generation of soliton packets in a two-level laser

A variant of perturbation theory is constructed for a system of nearly integrable equations. Perturbations of a special type are considered, which makes it possible to represent the system in the form of compatibility condition for “deformed” linear systems. The corresponding deformation of the Whitham equations is found. The mathematical apparatus is used to theoretically examine the generation of a sequence of solitons in a two-level laser. The generation process is described by a system of Maxwell-Bloch equations with pumping of the upper level and with allowance for some relaxation effects. The dynamics of the transformation of the initial perturbation into a sequence of solitons under pumping is studied. Finally, the various generation regimes are analyzed and compared with the experimental data. Zh. éksp. Teor. Fiz. 112, 2237–2251 (December 1997)     

Multiplication for solutions of the equation

Linear first-order systems of partial differential equations (PDEs) of the form f=Mg, where M is a constant matrix, are studied on vector spaces over the fields of real and complex numbers. The Cauchy–Riemann equations belong to this class. We introduce on the solution space a bilinear *-multiplication, playing the role of a nonlinear superposition principle, that allows for algebraic construction of new solutions from known solutions. The gradient equation f=Mg is a simple special case of a large class of systems of PDEs, admitting a *-multiplication of solutions. We prove that any gradient equation has the exceptional property that the general analytic solution can be expressed as *-power series of certain simple solutions.     

Solution and ellipticity properties of the self-duality equations of Corriganet al. in eight dimensions

We study the two sets of self-dual Yang-Mills equations in eight dimensions proposed in 1983 by E. Corriganet at. and show that one of these sets forms an elliptic system under the Coulomb gauge condition, and the other (overdetermined) set can have solutions that depend at most onN arbitrary constants, whereN is the dimension of the gauge group, hence the global solutions of both systems are finite dimensional. We describe a subvarietyP ₈ of the skew-symmetric 8 x 8 matrices by an eigenvalue criterion and we show that the solutions of the elliptic equations of Corriganet al. are among the maximal linear submanifolds ofP ₈. We propose an eighth-order action for which the elliptic set is a maximum.