Finite-time blow-up in dynamical systems

A new method to detect finite-time blow-up in systems of ordinary differential equations is presented. This simple algorithmic procedure is based on the analysis of singularities in complex time and amounts to checking the real-valuedness of the leading order term in the asymptotic series describing the behavior of the general solution around movable singularities. Illustrative examples and an application to a magneto-hydrodynamic problem are given.     

Interacting electron gas in a strong magnetic field

The density response function of an electron gas in a strong magnetic field shows logarithmic singularities due to scattering across the Fermi surface. We analyze the parquet equations for the vertex function in leading logarithmic order for a general interaction potential. The parquet equations are solved for a special interaction potential (Schulz and Keiter model). The divergence of the density response function at extremely high fields is discussed in connection with a possible transition to a Wigner lattice.     

Relativistic three-quark equations and low-lying baryon spectroscopy

The relativistic three-quark equations are found in the framework of the dispersion relations technique. The approximate solutions of these equations using the method based on the extraction of leading singularities of the amplitude are obtained. The calculated mass values of two baryonic multipletsJ ^P =1/2⁺, 3/2⁺ are in good agreement with the experimental ones.     

Classification of dark pion multiplets as dark matter candidates and collider phenomenology

In this work we propose to use leading singularities to obtain the classical pieces of amplitudes of two massive particles whose only interaction is gravitational. Leading singularities are generalizations of unitarity cuts. At one-loop we find that leading singularities obtained by multiple discontinuities in the t-channel contain all the classical information. As the main example, we show how to obtain a compact formula for the fully relativistic classical one-loop contribution to the scattering of two particles with different masses. The non-relativistic limit of the leading singularity agrees with known results in the post-Newtonian expansion. We also compute a variety of higher loop leading singularities including some all-loop families and study some of their properties.     

Relativistic three-quark equations and mass spectrum of the (70, 1<Superscript>−</Superscript>) baryon multiplet

A relativistic generalization of three-quark Faddeev equations for the particles of the (70, 1^?) multiplet is constructed within the dispersion-relation approach. An approximate solution to the relativistic three-quark equations is obtained on the basis of the method of isolating leading singularities of scattering amplitudes. The mass spectrum of the (70, 1^?) multiplet is calculated, and the result is in rather good agreement with corresponding experimental data.     

Relating hard QCD processes through universality of mass singularities

Hard QCD processes involving final jets are studied and compared by means of a simple approach to mass singularities. This is based on the Lee-Nauenberg-Kinoshita theorem and on a rather subtle use of gauge invariance in hard collinear gluon bremsstrahlung. One-loop results are easily derived for processes involving any number of initial quarks and/or currents. The method greatly simplifies the computation of higher-order loops at the leading log level and our preliminary results allow us to conclude that the crucial features encountered at the one-loop level will persist. We are thus able to relate different hard processes and to show that suitable ratios of cross sections, being free from mass singularities, can be computed perturbatively, as usually assumed in QCD-inspired parton models. We are also able to relate our universal leading mass singularities to leading scaling violations and to extend therefore the results of the operator product expansion method to processes outside the range of the light-cone analysis. Some delicate points caused by confinement-related singularities (e.g., narrow resonance poles) are also discussed.     

Three-body scattering at positive energies: Solving the problem of moving logarithmic singularities

For three-body scattering at positive total energies, integral equations are obtained whose kernels have no logarithmic singularities on the contour of integration. The corresponding singularities that are present in original integral equations can be circumvented by shifting a part of the contour of integration from the real axis to the complex plane. This is done only for a special auxiliary solution appearing to be an analytic function in this region. The physical amplitude proper is found as one of the solutions to the resulting set of equations. In contrast to conventional techniques, an additional analysis is therefore not required here, so that numerical solutions can be obtained within standard computational schemes.     

Singularities motion equations in 2-dimensional ideal hydrodynamics of incompressible fluid

In this Letter, we have obtained motion equations for a wide class of one-dimensional singularities in 2D ideal hydrodynamics. The simplest of them, are well known as point vortices. More complicated singularities correspond to vorticity point dipoles. It has been proved that point multipoles of a higher order (quadrupoles and more) are not the exact solutions of two-dimensional ideal hydrodynamics. The motion equations for a system of interacting point vortices and point dipoles have been obtained. It is shown that these equations are Hamiltonian ones and have three motion integrals in involution. It means the complete integrability of two-particle system, which has a point vortex and a point dipole.     

Propagation of a linear wave created by a spatially localized perturbation in a regular lattice and punctured Lagrangian manifolds

The following results are obtained for the Cauchy problem with localized initial data for the crystal lattice vibration equations with continuous and discrete time: (i) the asymptotics of the solution is determined by Lagrangian manifolds with singularities (“punctured” Lagrangian manifolds); (ii) Maslov’s canonical operator is defined on such manifolds as a modification of a new representation recently obtained for the canonical operator by the present authors together with A. I. Shafarevich (Dokl. Ross. Akad. Nauk 46 (6), 641–644 (2016)); (iii) the projection of the Lagrangian manifold onto the configuration plane specifies a bounded oscillation region, whose boundary (which is naturally referred to as the leading edge front) is determined by the Hamiltonians corresponding to the limit wave equations; (iv) the leading edge front is a special caustic, which possibly contains stronger focal points. These observations, together with earlier results, lead to efficient formulas for the wave field in a neighborhood of the leading edge front.     

Transport theory for quantum crystals

A correlation function approach is developed to treat non-equilibrium phenomena of quantum crystals at low frequency and long wavelength within the renormalized harmonic approximation (RHA). The derivation of the transport equations is carried out by studying the hierarchy of equations of motion for the retarded Green's functions of a pure, nonprimitive, nonionic, anharmonic lattice. Using a factorization technique to take into account the most important terms due to the particle fluctuations and the leading contributions to the hydrodynamic singularities of the phonon self-energy, we find a differential equation for the displacement field and a generalized transport equation for the phonon gas. The microscopic RHA expressions for the local temperature, the local heat density and the energy current are derived; the quasiparticle parameters (elastic constants, generalized Grüneisen parameters, quasiparticle interaction) entering the equations of motion are shown to be consistent with the RHA. In the hydrodynamic regime the general equations are reduced to two coupled differential equations for the lattice deformations and for the local temperature. Then only the displacement-displacement, the displacement-energy density and the energy density-energy density correlation functions show macroscopic fluctuations; for these functions thermodynamical sum-rules are derived.     

Singularities and spectral changes of Gaussian beams focused by a lens with spherical aberration

We address the effect of the truncation parameter and spherical aberration (SA) on the singularity transformation and spectral behavior of the polychromatic Gaussian beams focused by an aperture lens with SA in detail. The numerical simulation results, based on the derived equations of the intensity and the spectral density, are given. It is found that the axial singularities vanished with the change of the truncated parameter. The intensity and drastic spectral change fade away with an annihilation process of the phase singularities, and the drastic spectral change does not disappear immediately at the moment the phase singularity annihilates. The singularities in the focal region will redistribute with the increment of SA coefficient, some singularities will vanish, some will spilt into two new singularities, and other off-axial singularities will appear and split into two new singularities as well. When SA coefficient changed, we can find that the axial singularities disappear as well with the decreasing value of truncation parameter. These new splitted singularities due to the change of SA coefficient will converge into one singularity again and disappear gradually.     

Recovering singularities of a potential from singularities of scattering data

In this paper we show that the leading singularities of certain potentials can be determined from the leading singularities of the backscattering (as well as other determined sets of scattering data). The potentials in question are conormal with respect to smooth surfaces of arbitrary dimension; the restrictions on their orders allow for unbounded potentials in all dimension greater than or equal to three.Partially supported by NSF Grant DMS-9101298 and an Alfred P. Sloan Research FellowshipPartially supported by NSF Grant DMS-9100178     

Vortex structures with complex points singularities in two-dimensional Euler equations. New exact solutions

In this work we present a new class of exact stationary solutions for two-dimensional (2D) Euler equations. Unlike already known solutions, the new ones contain complex singularities. We consider point singularities which have a vector field index greater than 1 as complex. For example, the dipole singularity is complex because its index is equal to 2. We present in explicit form a large class of exact localized stationary solutions for 2D Euler equations with a singularity whose index is equal to 3. The solutions obtained are expressed in terms of elementary functions. These solutions represent a complex singularity point surrounded by a vortex satellite structure. We also discuss the motion equation of singularities and conditions for singularity point stationarity which provide the stationarity of the complex vortex configuration.     

A Well-Posed Numerical Method to Track Isolated Conformal Map Singularities in Hele-Shaw Flow

We present a new numerical method for calculating an evolving 2D Hele-Shaw interface when surface tension effects are neglected. In the case where the flow is directed from the less viscous fluid into the more viscous fluid, the motion of the interface is ill-posed; small deviations in the initial condition will produce significant changes in the ensuing motion. This situation is disastrous for numerical computation, as small roundoff errors can quickly lead to large inaccuracies in the computed solution. Our method of computation is most easily formulated using a conformal map from the fluid domain into a unit disk. The method relies on analytically continuing the initial data and equations of motion into the region exterior to the disk, where the evolution problem becomes well-posed. The equations are then numerically solved in the extended domain. The presence of singularities in the conformal map outside of the disk introduces specific structures along the fluid interface. Our method can explicitly track the location of isolated pole and branch point singularities, allowing us to draw connections between the development of interfacial patterns and the motion of singularities as they approach the unit disk. In particular, we are able to relate physical features such as finger shape, side-branch formation, and competition between fingers to the nature and location of the singularities. The usefulness of this method in studying the formation of topological singularities (self intersections of the interface) is also pointed out.     

Polynomial integrals of mechanical systems on a torus with a singular potential

The problem on integrability of the equations of motion of a material point on an n-dimensional Euclidean torus under the action of a force field with the potential energy having singularities at a finite number of points is considered. It is assumed that these singularities contain logarithmic coefficients and, consequently, have a more general form in comparison with power features. The potentials having power-type singularities were considered previously by V.V. Kozlov and D.V. Treshchev. In this work, it is proved that the equations of motion in the problem under consideration admit no nontrivial momentum-polynomial first integral with integrable coefficients on this torus.     

Elimination of body-frame singularities in the separation of three collective angles in quantum N-body problems

The appearance of body-frame singularities in gauge potentials when three collective angles are separated by means of Wigner D-functions is a fundamental difficulty in the quantum N-body problem. We show that the use of the overcomplete basis of minimal multipolar harmonics allows these singularities to be avoided at the expense of increasing the dimension of the resulting system of coupled equations describing the internal motion of the system.     

Accurate numerical resolution of transients in initial-boundary value problems for the heat equation

If the initial and boundary data for a PDE do not obey an infinite set of compatibility conditions, singularities will arise in the solution at the corners of the initial time–space domain. For dissipative equations, such as the 1-D heat equation or 1-D convection–diffusion equations, the impacts of these singularities are short lived. However, they can cause a very severe loss of numerical accuracy if we are interested in transient solutions. The phenomenon has been described earlier from a theoretical standpoint. Here, we illustrate it graphically and present a simple remedy which, with only little extra cost and effort, restores full numerical accuracy.     

Are the Singularities Stable?

The spacetime singularities play a useful role in gravitational theories by distinguishing physical solutions from non-physical ones. The problem, we are studying in this paper is whether these singularities are stable. To answer this question, we have analyzed the general problem of stability of the family of the static spherically symmetric solutions of the standard Einstein-Maxwell model coupled to an extra free massless scalar field. We have obtained the equations for the axial and polar perturbations. The stability against axial perturbations has been proven.     

On the exponentiation of leading infrared divergences in massless Yang-Mills theories

We derive, in the axial gauge, the effective U matrix which governs the behaviour of leading infrared singularities in the self-energy functions of Yang-Mills particles. We then show in a very simple manner, that these divergences, which determine the leading singularities in massless Yang-Mills theories, exponentiate.