Search

The group-theoretic method for constructing symmetric isometric embeddings is used to describe all possible four-dimensional surfaces in flat $(1,9)$-dimensional space, whose induced metric is static and spherically symmetric. For such surfaces, we propose a classification related to the dimension of the elementary blocks forming the embedding function. All suitable 52 classes of embeddings are summarized in one table and analyzed for the unfolding property (wich means that the surface does not belong locally to some subspace of the ambient space), as well as for the presence of smooth embeddings of the Minkowski metric. The obtained results are useful for the analysis of the equations of motion in the Regge-Teitelboim embedding gravity, where the presence of unfolded embeddings of the Minkowski metric is essential., Comment: LaTeX, 16 pages

This paper examines the reconstruction of a family of dynamical systems with neuromorphic behavior using a single scalar time series. A model of a physiological neuron based on the Hodgkin-Huxley formalism is considered. Single time series of one of its variables is shown to be enough to train a neural network that can operate as a discrete time dynamical system with one control parameter. The neural network system is created in two steps. First, the delay-coordinate embedding vectors are constructed form the original time series and their dimension is reduced with by means of a variational autoencoder to obtain the recovered state-space vectors. It is shown that an appropriate reduced dimension can be determined by analyzing the autoencoder training process. Second, pairs of the recovered state-space vectors at consecutive time steps supplied with a constant value playing the role of a control parameter are used to train another neural network to make it operate as a recurrent map. The regimes of thus created neural network system observed when its control parameter is varied are in very good accordance with those of the original system, though they were not explicitly presented during training., Comment: 15 pages, 15 figures, 3 tables

We study a perturbation theory for embedding gravity equations in a background for which corrections to the embedding function are linear with respect to corrections to the flat metric. The arbitrariness remaining after solving the linearized field equations is fixed by an assumption that the solution is static in the second order. A nonlinear differential equation is obtained, which makes it possible to find the gravitational potential for a spherically symmetric case if a background embedding is given. An explicit form of a spherically symmetric background parameterized by one function of radius is proposed. It is shown that this function can be chosen in such a way that the gravitational potential is in a good agreement with the observed distribution of dark matter in a galactic halo., Comment: LaTeX, 14 pages. This version corresponds to the published one

We consider Hodgkin-Huxley-type model that is a stiff ODE system with two fast and one slow variables. For the parameter ranges under consideration the original version of the model has unstable fixed point and the oscillating attractor that demonstrates bifurcation from bursting to spiking dynamics. Also a modified version is considered where the bistability occurs such that an area in the parameter space appears where the fixed point becomes stable and coexists with the bursting attractor. For these two systems we create artificial neural networks that are able to reproduce their dynamics. The created networks operate as recurrent maps and are trained on trajectory cuts sampled at random parameter values within a certain range. Although the networks are trained only on oscillatory trajectory cuts, it also discover the fixed point of the considered systems. The position and even the eigenvalues coincide very well with the fixed point of the initial ODEs. For the bistable model it means that the network being trained only on one brunch of the solutions recovers another brunch without seeing it during the training. These results, as we see it, are able to trigger the development of new approaches to complex dynamics reconstruction and discovering. From the practical point of view reproducing dynamics with the neural network can be considered as a sort of alternative method of numerical modeling intended for use with contemporary parallel hard- and software., Comment: 17 pages, 12 figures, 2 tables

Faddeev variant of embedding theory is an example of using the embedding approach for the description of gravity. In the original form of the embedding approach, the gravity is described by an embedding function of a four-dimensional surface representing our spacetime. In Faddeev variant, the independent variable is a non-square vielbein, which is a derivative of embedding function in embedding theory. We study the possibility of the existence of extra solutions in Faddeev variant, which makes this theory non-equivalent to GR. To separate the degrees of freedom corresponding to extra matter, we propose a formulation of this theory as GR with an additional contribution to the action. We analyze the equations of motion for a specific class of solutions corresponding to a weak gravitational field. We construct a simple exact solution corresponding to arbitrary matter and nontrivial torsion, which is an extra solution in Faddeev variant in the absence of real matter., Comment: LaTeX, 7 pages, based on a talk given at the International Conference PhysicA.SPb/2021, Saint Petersburg, Russia, 18-22 October 2021

A system of three non-identical Josephson junctions connected via an RLC circuit is considered. The method of Lyapunov exponents charts is used, which makes it possible to identify the main types of dynamics of the system and to analyze the dependence of its properties on parameters. The possibility of both two and three-frequency invariant tori is demonstrated. Saddle-node bifurcations of resonant tori are studied with the use of instantaneous Lyapunov exponents. The dependence of the charts on the type of coupling in the system is discussed., Comment: 17 pages, 9 figures, 1 table

We suggest a universal map capable to recover a behavior of a wide range of dynamical systems given by ODEs. The map is built as an artificial neural network whose weights encode a modeled system. We assume that ODEs are known and prepare training datasets using the equations directly without computing numerical time series. Parameter variations are taken into account in the course of training so that the network model captures bifurcation scenarios of the modeled system. Theoretical benefit from this approach is that the universal model admits using common mathematical methods without needing to develop a unique theory for each particular dynamical equations. Form the practical point of view the developed method can be considered as an alternative numerical method for solving dynamical ODEs suitable for running on contemporary neural network specific hardware. We consider the Lorenz system, the Roessler system and also Hindmarch-Rose neuron. For these three examples the network model is created and its dynamics is compared with ordinary numerical solutions. High similarity is observed for visual images of attractors, power spectra, bifurcation diagrams and Lyapunov exponents., Comment: 17 pages, 10 figures

(Abridged) Modern studies and industrial applications related to the design, radiation protection, and reliability of nuclear facilities, medical applications, as well as space research and exploration are relying on extensive simulations and modeling. Computer codes realizing semi-classical and quantum molecular dynamics (QMD) approaches are often employed to make up for the lack of accelerator data on many nuclear reactions at intermediate and high energies (>10s of MeV/n) and are in high demand. This contribution focuses on the methodology of generating reliable proton-induced cross sections in the energy range below 1 GeV. We developed a problem-oriented computer framework based on MCNPX and CASCADE/INPE codes to calculate activation cross section data at intermediate and high energies. Goodness of the fits of nucleon-nucleus interaction models to the existing data is evaluated based on elaborated algorithms. The method is based on the analysis of a large set of data and calculated cross sections for different targets and residual nuclei in a wide range of proton energies using numerous criteria. In practice, this could be done by tuning the model parameters and algorithms to fit required experimental data subset, or through achieving unification and consistency of fundamental parameters for all considered reactions. The presented framework is pursuing the latter approach. We use proton-induced reactions on Si and Fe nuclei to illustrate the application of the developed multi-criteria algorithm, where we use all data below 1 GeV available from the EXFOR data library and the IAEA CRP "Benchmark of Spallation Models." We show that the analysis of the predictive power of various intermediate and high-energy models based on the multi-criteria algorithm allows more sophisticated selection of appropriate models in a given energy range and residual nuclei domain., Comment: A poster E1.16-0085-21 presented at an event E1.16 "Origin of Cosmic Rays," 43th COSPAR Scientific Assembly (hybrid), 28 Jan - 4 Feb 2021, Sydney, Australia. For a full agenda of the event E1.16, see https://www.cospar-assembly.org/admin/session_cospar.php?session=903

We consider a self-oscillator whose excitation parameter is varied. Frequency of the variation is much smaller then the natural frequency of the oscillator so that oscillations in the system are periodically excited and decay. Also a time delay as added such that when the oscillations start to grow at a new excitation stage they are influenced via the delay line by the oscillations at the penultimate excitation stage. Due to a nonlinearity the seeding from the past arrives with a doubled phase so that oscillation phase changes from stage to stage according to chaotic Bernoulli-type map. As a result, the system operates as two coupled hyperbolic chaotic subsystems. Varying the relation between the delay time and the excitation period we affect a coupling strength between these subsystems as well as intensity of the phase doubling mechanism responsible for the hyperbolicity. Due to this, a transition from non-hyperbolic to hyperbolic hyperchaos occurs. The following steps of the transition scenario are revealed and analyzed: (a) an intermittency as alternation of long staying near a fixed point at the origin and short chaotic bursts; (b) chaotic oscillations with often visits to the fixed point; (c) plain hyperchaos without hyperbolicity after termination visiting the fixed point; (d) transformation of hyperchaos to hyperbolic form., Comment: 13 pages, 11 figures

Dimension of an inertial manifold for a chaotic attractor of spatially distributed system is estimated using autoencoder neural network. The inertial manifold is a low dimensional manifold where the chaotic attractor is embedded. The autoencoder maps system state vectors onto themselves letting them pass through an inner state with a reduced dimension. The training processes of the autoencoder is shown to depend dramatically on the reduced dimension: a learning curve saturates when the dimension is too small and decays if it is sufficient for a lossless information transfer. The smallest sufficient value is considered as a dimension of the inertial manifold, and the autoencoder implements a mapping onto the inertial manifold and back. The correctness of the computed dimension is confirmed by its remarkable coincidence with the one obtained as a number of covariant Lyapunov vectors with vanishing pairwise angles. These vectors are called physical modes. Unlike never having zero angles residual ones they are known to span a tangent subspace for the inertial manifold., Comment: 10 pages and 9 figures

In this paper we analyze local structure of several chaotic attractors recently suggested in literature as pseudohyperbolic. The absence of tangencies and thus the presence of the pseudohyperbolicity is verified using the method of angles that includes computation of distributions of the angles between the corresponding tangent subspaces. Also we analyze how volume expansion in the first subspace and the contraction in the second one occurs locally. For this purpose we introduce a family of instant Lyapunov exponents. Unlike the well known finite time ones, the instant Lyapunov exponents show expansion or contraction on infinitesimal time intervals. Two types of instant Lyapunov exponents are defined. One is related to ordinary finite time Lyapunov exponents computed in the course of standard algorithm for Lyapunov exponents. Their sums reveal instant volume expanding properties. The second type of instant Lyapunov exponents shows how covariant Lyapunov vectors grow or decay on infinitesimal time. Using both instant and finite time Lyapunov exponents we demonstrate that specific to the pseudohyperbolicity average expanding and contracting properties are typically violated on infinitesimal time. Instantly volumes from the first subspace can sometimes be contacted, directions in the second subspace can sometimes be expanded, and the instant contraction in the first subspace can sometimes be stronger than the contraction in the second subspace., Comment: 25 pages, 26 figures

We develop an extension of the fast method of angles for hyperbolicity verification in chaotic systems with an arbitrary number of time-delay feedback loops. The adopted method is based on the theory of covariant Lyapunov vectors and provides an efficient algorithm applicable for systems with high-dimensional phase space. Three particular examples of time-delay systems are analyzed and in all cases the expected hyperbolicity is confirmed., Comment: 16 pages, 3 figures

After having announced the statistically significant observation (5.6~$\sigma$) of the new exotic $\pi K$ atom, the DIRAC experiment at the CERN proton synchrotron presents the measurement of the corresponding atom lifetime, based on the full $\pi K$ data sample: $\tau = (5.5^{+5.0}_{-2.8}) \cdot 10^{-15}s$. By means of a precise relation ($<1\%$) between atom lifetime and scattering length, the following value for the S-wave isospin-odd $\pi K$ scattering length $a_0^{-}~=~\frac{1}{3}(a_{1/2}-a_{3/2})$ has been derived: $\left|a_0^-\right| = (0.072^{+0.031}_{-0.020}) M_{\pi}^{-1}$., Comment: 18 pages, 17 figures