Complex Singularities Research Articles

Tracking complex singularities of fluids on log-lattices

Abstract In 1981, Frisch and Morf (1981 Phys. Rev. A 23 2673–705) postulated the existence of complex singularities in solutions of Navier–Stokes equations. Present progress on this conjecture is hindered by the computational burden involved in simulations of the Euler equations or the Navier–Stokes equations at high Reynolds numbers. We investigate this conjecture in the case of fluid dynamics on log-lattices, where the computational burden is logarithmic concerning ordinary fluid simulations. We analyze properties of potential complex singularities in both 1D and 3D models for lattices of different spacings. Dominant complex singularities are tracked using the singularity strip method to obtain new scalings regarding the approach to the real axis and the influence of normal, hypo and hyper dissipation.

Proceedings of the 17th International Workshop on Real and Complex Singularities

Proceedings of the 17th International Workshop on Real and Complex Singularities

Disconnected gauge groups in the infrared

Gauging a discrete 0-form symmetry of a QFT is a procedure that changes the global form of the gauge group but not its perturbative dynamics. In this work, we study the Seiberg-Witten solution of theories resulting from the gauging of charge conjugation in 4d N\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{N} $$\\end{document} = 2 theories with SU(N) gauge group and fundamental hypermultiplets. The basic idea of our procedure is to identify the ℤ2 action at the level of the SW curve and perform the quotient, and it should also be applicable to non-lagrangian theories. We study dynamical aspects of these theories such as their moduli space singularities and the corresponding physics; in particular, we explore the complex structure singularity arising from the quotient procedure. We also discuss some implications of our work in regards to three problems: the geometric classification of 4d SCFTs, the study of non-invertible symmetries from the SW geometry, and the String Theory engineering of theories with disconnected gauge groups.

Complex ODEs, singularity theory and dynamics

Complex ODEs, singularity theory and dynamics

Analyticity theorems for parameter-dependent plurisubharmonic functions

In this paper, we first show that a union of upper-level sets associated to fibrewise Lelong numbers of plurisubharmonic functions is in general a pluripolar subset. Then we obtain analyticity theorems for a union of sub-level sets associated to fibrewise complex singularity exponents of some special (quasi-)plurisubharmonic functions. As a corollary, we confirm that, under certain conditions, the logarithmic poles of relative Bergman kernels form an analytic subset when the (quasi-)plurisubharmonic weight function has analytic singularities. In the end, we give counterexamples to show that the aforementioned sets are in general non-analytic even if the plurisubharmonic function is supposed to be continuous.

A note on complex plane curve singularities up to diffeomorphism and their rigidity

We prove that if two germs of plane curves (C, 0) and (C′,0)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$(C',0)$$\\end{document} with at least one singular branch are equivalent by a (real) smooth diffeomorphism, then C is complex isomorphic to C′\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$C'$$\\end{document} or to C′¯\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\overline{C'}$$\\end{document}. A similar result was shown by Ephraim for irreducible hypersurfaces before, but his proof is not constructive. Indeed, we show that the complex isomorphism is given by the Taylor series of the diffeomorphism. We also prove an analogous result for the case of non-irreducible hypersurfaces containing an irreducible component that is non-factorable. Moreover, we provide a general overview of the different classifications of plane curve singularities.

Generalized Drude–Lorentz Model Complying with the Singularity Expansion Method

AbstractDeriving analytical expressions of dielectric permittivities is required for numerical and physical modeling of optical systems and the soar of non‐Hermitian photonics motivates their prolongation in the complex plane. Analytical models are based on the association of microscopic models to describe macroscopic effects. However, the question is to know whether the resulting Debye–Drude–Lorentz models are not too restrictive. Here, it is shown that the permittivity must be treated as a meromorphic transfer function that complies with the requirements of complex analysis. This function can be naturally expanded on a set of complex singularities. This singularity expansion of the dielectric permittivity allows to derive a generalized expression of the Debye–Drude–Lorentz model that complies with the requirements of complex analysis and the constraints of physical systems. It is shown that the complex singularities and other parameters of this generalized expression can be retrieved from experimental data acquired along the real frequency axis. The accuracy of this expression is assessed for a wide range of materials including metals, 2D materials and dielectrics, and it is shown how the distribution of the retrieved poles helps in characterizing the materials.

Topological atom optics and beyond with knotted quantum wavefunctions

Atom optics demonstrates optical phenomena with coherent matter waves, providing a foundational connection between light and matter. Significant advances in optics have followed the realization of structured light fields hosting complex singularities and topologically non-trivial characteristics. However, analogous studies are still in their infancy in the field of atom optics. Here, we investigate and experimentally create knotted quantum wavefunctions in spinor Bose–Einstein condensates which display non-trivial topologies. In our work we construct coordinated orbital and spin rotations of the atomic wavefunction, engineering a variety of discrete symmetries in the combined spin and orbital degrees of freedom. The structured wavefunctions that we create map to the surface of a torus to form torus knots, Möbius strips, and a twice-linked Solomon’s knot. In this paper we demonstrate close connections between the symmetries and underlying topologies of multicomponent atomic systems and of vector optical fields—a realization of topological atom-optics.

Singularity structure simplification for hex mesh via integer linear program

Topology optimization of hexahedral (hex) meshes has been a widely studied topic, with the primary goal of optimizing the singularity structure. Previous works have focused on simplifying complex singularity structures by collapsing sheets/chords. However, these works require a large number of checks during the process to prevent illegal operations. Moreover, the employed simplification strategies are not based on the topological characteristics of the structure, but rather on the rank of the components that can be simplified. To overcome these problems, we analyze how topology operations affect the degree of edges in hex meshes, and introduce a fast and automatic algorithm to simplify the singularity structure of hex meshes. The algorithm relies on sheet operations, using mesh volume as a metric to assess the degree of simplification. Moreover, it designs constraints to prevent illegal operations and employs integer linear program to plan the overall optimization strategy for a mesh. After that, we relax the singularity constraints to further simplify the structure, and handle unreasonable singularities via sheet inflation operation. Our algorithm can also improve singularity structure without merging singularities by adjusting the singularity constraint conditions. Numerous experiments demonstrate the effectiveness and efficiency of our algorithm.

Angle-dependent phononic dynamics for data-driven source localization

The source angle localization problem is studied based on scattering of elastic waves in two dimensions by a phononic array and the exceptional points of its band structure. Exceptional points are complex singularities of a parameterized eigen-spectrum, where two modes coalesce with identical mode shapes. These special points mark the qualitative transitions in the system behavior and have been proposed for sensing applications. The equi-frequency band structures are analyzed with focus on the angle-dependent modal behaviors. At the exceptional points and critical angles, the eigen-modes switch their energy characteristics and symmetry, leading to enhanced sensitivity as the scattering response of the medium is inherently angle-dependent. An artificial neural network is trained with randomly weighted and superposed eigen-modes to achieve deep learning of the angle-dependent dynamics. The trained algorithm can accurately classify the incident angle of an unknown scattering signal, with minimal sidelobe levels and suppressed main lobewidth. The neural network approach shows superior localization performance compared with standard delay-and-sum technique. The proposed application of the phononic array highlights the physical relevance of band topology and eigen-modes to a technological application, adds extra strength to the existing localization methods, and can be easily enhanced with the fast-growing data-driven techniques.

Locating complex singularities of Burgers’ equation using exponential asymptotics and transseries

Burgers’ equation is an important mathematical model used to study gas dynamics and traffic flow, among many other applications. Previous analysis of solutions to Burgers’ equation shows an infinite stream of simple poles born at t = 0 + , emerging rapidly from the singularities of the initial condition, that drive the evolution of the solution for t > 0 . We build on this work by applying exponential asymptotics and transseries methodology to an ordinary differential equation that governs the small-time behaviour in order to derive asymptotic descriptions of these poles and associated zeros. Our analysis reveals that subdominant exponentials appear in the solution across Stokes curves; these exponentials become the same size as the leading order terms in the asymptotic expansion along anti-Stokes curves, which is where the poles and zeros are located. In this region of the complex plane, we write a transseries approximation consisting of nested series expansions. By reversing the summation order in a process known as transasymptotic summation, we study the solution as the exponentials grow, and approximate the pole and zero location to any required asymptotic accuracy. We present the asymptotic methods in a systematic fashion that should be applicable to other nonlinear differential equations.

Asymmetric phase modulation of light with parity-symmetry broken metasurfaces

The design of wavefront-shaping devices is conventionally approached using real-frequency modeling. However, since these devices interact with light through radiative channels, they are by default non-Hermitian objects having complex eigenvalues (poles and zeros) that are marked by phase singularities in a complex frequency plane. Here, by using temporal coupled mode theory, we derive analytical expressions allowing to predict the location of these phase singularities in a complex plane and as a result, allowing to control the induced phase modulation of light. In particular, we show that spatial inversion symmetry breaking—implemented herein by controlling the coupling efficiency between input and output radiative channels of two-port components called metasurfaces—lifts the degeneracy of reflection zeros in forward and backward directions, and introduces a complex singularity with a positive imaginary part necessary for a full 2π-phase gradient. Our work establishes a general framework to predict and study the response of resonant systems in photonics and metaoptics.

Breakdown of homoclinic orbits to L3 in the RPC3BP (II). An asymptotic formula

The Restricted 3-Body Problem models the motion of a body of negligible mass under the gravitational influence of two massive bodies called the primaries. If one assumes that the primaries perform circular motions and that all three bodies are coplanar, one has the Restricted Planar Circular 3-Body Problem (RPC3BP). In rotating coordinates, it can be modeled by a two degrees of freedom Hamiltonian, which has five critical points called the Lagrange points L1,…,L5.The Lagrange point L3 is a saddle-center critical point which is collinear with the primaries and beyond the largest of the two. In this paper, we obtain an asymptotic formula for the distance between the stable and unstable manifolds of L3 for small values of the mass ratio 0<μ≪1. In particular we show that L3 cannot have (one round) homoclinic orbits.If the ratio between the masses of the primaries μ is small, the hyperbolic eigenvalues of L3 are weaker, by a factor of order μ, than the elliptic ones. This rapidly rotating dynamics makes the distance between manifolds exponentially small with respect to μ. Thus, classical perturbative methods (i.e. the Melnikov-Poincaré method) can not be applied.The obtention of this asymptotic formula relies on the results obtained in the prequel paper [10] on the complex singularities of the homoclinic of a certain averaged equation and on the associated inner equation.In this second paper, we relate the solutions of the inner equation to the analytic continuation of the parameterizations of the invariant manifolds of L3 via complex matching techniques. We complete the proof of the asymptotic formula for their distance showing that its dominant term is the one given by the analysis of the inner equation.

Numerical analytic continuation

Let f be an analytic function on a simply-connected compact continuum E of the complex z-plane. This might be an interval of the real line, where f might be real analytic. How can we calculate good estimates of the analytic continuation of f to other points zin {mathbb C}? How can we estimate the locations of real or complex singularities of f? We review both the theory and the practice of some existing methods for these problems and propose that excellent results can be obtained from the computation of rational approximations of f by the AAA algorithm. In the case of analytic functions of two or more variables, the rational approximations are applied along line segments or other analytic arcs.

The role of exponential asymptotics and complex singularities in self-similarity, transitions, and branch merging of nonlinear dynamics

We study a prototypical example in nonlinear dynamics where transition to self-similarity in a singular limit is fundamentally changed as a parameter is varied. Here, we focus on the complicated dynamics that occur in a generalised unstable thin-film equation that yields finite-time rupture. A parameter, n, is introduced to model more general disjoining pressures. For the standard case of van der Waals intermolecular forces, n=3, it was previously established that a countably infinite number of self-similar solutions exist leading to rupture. Each solution can be indexed by a parameter, ϵ=ϵ1>ϵ2>⋯>0, and the prediction of the discrete set of solutions requires examination of terms beyond-all-orders in ϵ. However, recent numerical results have demonstrated the surprising complexity that exists for general values of n. In particular, the bifurcation structure of self-similar solutions now exhibits branch merging as n is varied. In this work, we shall present key ideas of how branch merging can be interpreted via exponential asymptotics.

Las Leñas International Tunnel – Geomechanical challenges in large‐scale tunneling projects across The Andes

AbstractDeveloping large‐scale tunnelling projects in orogenic regions implies enormous technical challenges. Las Leñas International Tunnel Project through the Andes Main Cordillera of Chile and Argentina has been a long‐awaited infrastructure project, whose geological and topographic settings suggest several complex singularities to be considered, analyzed, and assessed during its development. Throughout the course of this feasibility study, the complex geological subsurface interpretation offers a broad spectrum of technical demands to be fulfilled. The overburden conditions, the presence of certain lithologies such as expansive clays or soluble evaporitic sequences, along with regional stress assessment, are major technical challenges. Through extensive geological, geotechnical, and hydrogeological mapping, subsurface 3D modelling interpretation, along with numerical analysis and its relationship with similar projects in the region, brittle, and plastic behaviour are analyzed. A wide range of modern techniques are used to assess those mentioned tasks, which enable an extensive comprehension of the geological and geotechnical conditions of the area. The purpose of this contribution is to enrich the knowledge on geological and geomechanical challenges faced in large tunnelling infrastructure projects in high‐mountain environments and what specific lessons may arise from them for similar future projects.

Design, Analysis, and Optimization of a Kinematically Redundant Parallel Robot

It has been nearly 60 years since the introduction of the Gough–Stewart manipulator (GSM). With the advantages of superior load capacity and high precision, the GSM still plays an important role in many fields. However, the GSM has limitations such as a small workspace and complex singularities. To overcome these problems, a novel kinematically redundant parallel robot is designed with three redundant actuators added on the basis of the GSM. First, the structure of the proposed robot is introduced, and the kinematics of the proposed robot is established. Second, the workspaces of the proposed robot are analyzed, and the results show that the position workspace volume and the maximum torsion and orientation angles of the proposed robot can be improved effectively. Third, the singularities of the proposed robot and the GSM are compared based on the Jacobian matrices. Finally, the multi-parameter, multi-objective optimization is used to optimize the geometric parameters of the proposed robot, and an experimental prototype of the proposed robot is designed based on the optimization results.

Burgers’ equation in the complex plane

Burgers’ equation is a well-studied model in applied mathematics with connections to the Navier–Stokes equations in one spatial direction and traffic flow, for example. Following on from previous work, we analyse solutions to Burgers’ equation in the complex plane, concentrating on the dynamics of the complex singularities and their relationship to the solution on the real line. For an initial condition with a simple pole in each of the upper- and lower-half planes, we apply formal asymptotics in the small- and large-time limits in order to characterise the initial and later motion of the singularities. The small-time limit highlights how infinitely many singularities are born at t=0 and how they orientate themselves to lie increasingly close to anti-Stokes lines in the far field of the inner problem. This inner problem also reveals whether or not the closest singularity to the real axis moves toward the axis or away. For intermediate times, we use the exact solution, apply method of steepest descents, and implement the AAA approximation to track the complex singularities. Connections are made between the motion of the closest singularity to the real axis and the steepness of the solution on the real line. While Burgers’ equation is integrable (and has an exact solution), we deliberately apply a mix of techniques in our analysis in an attempt to develop methodology that can be applied to other nonlinear partial differential equations that do not.

On the use of asymptotically motivated gauge functions to obtain convergent series solutions to nonlinear ODEs

AbstractWe examine the power series solutions of two classical nonlinear ordinary differential equations of fluid mechanics that are mathematically related by their large-distance asymptotic behaviours in semi-infinite domains. The first problem is that of the ‘Sakiadis’ boundary layer over a moving flat wall, for which no exact analytic solution has been put forward. The second problem is that of a static air–liquid meniscus with surface tension that intersects a flat wall at a given contact angle and limits to a flat pool away from the wall. For the latter problem, the exact analytic solution—given as distance from the wall as a function of meniscus height—has long been known (Batchelor, G. K. (1967) An Introduction to Fluid Dynamics, chapter 1: The physical properties of fluids. Cambridge). Here, we provide an explicit solution as meniscus height versus distance from the wall to elucidate structural similarities to the Sakiadis boundary layer. Although power series solutions are readily obtainable to the governing nonlinear ordinary differential equations, we show that—in both problems—the series diverge due to non-physical complex or negative real-valued singularities. In both cases, these singularities can be moved by expanding in exponential gauge functions motivated by their respective large distance asymptotic behaviours to enable series convergence over their full semi-infinite domains. For the Sakiadis problem, this not only provides a convergent Taylor series (and conjectured exact) solution to the ODE, but also a means to evaluate the wall shear parameter (and other properties) to within any desired precision. Although the nature of nonlinear ODEs precludes general conclusions, our results indicate that asymptotic behaviours can be useful when proposing variable transformations to overcome power series divergence.

A multiscale scaled boundary finite element method solving steady-state heat conduction problem with heterogeneous materials

A new numerical model is proposed by bridging small and large scales via the numerical base functions constructed by Scaled Boundary Finite Element Method (SBFEM) for the multiscale heterogeneous analysis of steady-state heat conduction problems. Utilizing the advantages of polygonal Scaled Boundary Finite Element (SBFE) and image-based quadtree gridding, the construction of base functions can be conveniently and efficiently conducted particularly for the heterogeneous media with complex geometries or singularity of heat flux at the small-scale, and the solution accuracy at the large-scale can be improved by increasing nodes of coarse element only without increasing new nodes inside. Consequently, the SBFEM based solutions at the large-scale can be solved with a reduced order and sufficient accuracy, and the solutions at the small-scale can be achieved using SBFEM based base functions. Numerical examples are provided and the effectiveness of the proposed approach is stressed at both the large and small scales.