Search

We introduce a methodology to study the possible matter flows of an ecosystem defined by observational biomass data and realistic biological constraints. The flows belong to a polyhedron in a multi dimensional space that may make statistical exploration difficult in practice; instead, we propose to solve a convex optimization problem. Seven criteria based on ecological network indices have been selected to be used as convex goal functions. Numerical results show that the method is fast and can be used for large systems. Minimum flow solutions are analyzed using flow decomposition in paths and circuits. Their consistency is also tested by introducing a system of differential equations for the biomasses and examining the stability of the biomass fixed point. The method is illustrated and explained throughout the text on an ecosystem toy model. It is also applied to realistic food models., (© 2021. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

We revisit the problem of pinning a reaction-diffusion front by a defect, in particular by a reaction-free region. Using collective variables for the front and numerical simulations, we compare the behaviors of a bistable and monostable front. A bistable front can be pinned as confirmed by a pinning criterion, the analysis of the time independent problem, and simulations. Conversely, a monostable front can never be pinned, it gives rise to a secondary pulse past the defect and we calculate the time this pulse takes to appear. These radically different behaviors of bistable and monostable fronts raise issues for modelers in particular areas of biology, as for example, the study of tumor growth in the presence of different tissues.

We analyze an epidemic model on a network consisting of susceptible-infected-recovered equations at the nodes coupled by diffusion using a graph Laplacian. We introduce an epidemic criterion and examine different isolation strategies: we prove that it is most effective to isolate a node of highest degree. The model is also useful to evaluate deconfinement scenarios and prevent a so-called second wave. The model has few parameters enabling fitting to the data and the essential ingredient of importation of infected; these features are particularly important for the current COVID-19 epidemic., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2020 Elsevier B.V. All rights reserved.)

Background: On January 23, 2020, China imposed a quarantine on the city of Wuhan to contain the SARS-CoV-2 outbreak. Regardless of this measure, the new infection has spread to several countries around the world., Objective: We developed a method to study the dissemination of this infection by airline routes and provide estimations of the time of arrival of the outbreak to different cities., Methods: Using the Kermack and McKendrick model complemented with diffusion on a graph composed of nodes and edges, we made an analysis of COVID-19 dispersion to other cities by air travel., Results: The estimation was accurate in that it was possible to predict in the middle of February 2020 the arrival of the first outbreak in Mexico, which eventually occurred between March 20 and 30. This estimation was robust with respect to small changes in epidemiological parameters at the other nodes., Conclusions: The estimation of the time of arrival of the outbreak from its epicenter, allows for a time period to implement and strengthen preventive measures aimed at the general population as well as to strengthen hospital infrastructure and training of human resources. In the present study, this estimation was accurate, as observed from the real data of the beginning of the outbreak in Mexico City up to April 6, 2020., (Copyright: © 2020 Permanyer.)

To study how nonlinear waves propagate across Y- and T-type junctions, we consider the two-dimensional (2D) sine-Gordon equation as a model and examine the crossing of kinks and breathers. Comparing energies for different geometries reveals that, for small widths, the angle of the fork plays no role. Motivated by this, we introduce a one-dimensional effective model whose solutions agree well with the 2D simulations for kink and breather solutions. These exhibit two different behaviors: a kink crosses if it has sufficient energy; conversely a breather crosses when v>1-ω, where v and ω are, respectively, its velocity and frequency. This methodology can be generalized to more complex nonlinear wave models.

We consider radial sine-Gordon kinks in two, three, and higher dimensions. A full two-dimensional simulation showing that azimuthal perturbations remain small allows us to reduce the problem to the one-dimensional radial sine-Gordon equation. We solve this equation on an interval [r(0),r(1)] and absorb all outgoing radiation. As the kink shrinks toward r(0), before the collision, its motion is well described by a simple law derived from the conservation of energy. In two dimensions for r(0)≤2, the collision disintegrates the kink into a fast breather, while for r(0)≥4 we obtain a kink-breather metastable state where breathers are shed at each kink "return." In three and higher dimensions d, an additional kink-oscillon state appears for small r(0). On the application side, the kink disintegration opens the way for new types of terahertz microwave generators.

The interaction of a Zeldovich-Frank-Kamenetsky reaction-diffusion front with a localized defect is studied numerically and analytically. For the analysis, we start from conservation laws and develop simple, collective variable, ordinary differential equations for the front position and width. Their solutions are in good agreement with the solutions of the full problem. Finally, using this reduced model, we explain the pinning of the front on a large defect and obtain a quantitative criterion.

We analyze the 1D focusing nonlinear Schrödinger equation in a finite interval with homogeneous Dirichlet or Neumann boundary conditions. There are two main dynamics, the collapse which is very fast and a slow cascade of Fourier modes. For the cubic nonlinearity the calculations show no long-term energy exchange between Fourier modes as opposed to higher nonlinearities. This slow dynamics is explained by fairly simple amplitude equations for the resonant Fourier modes. Their solutions are well behaved so filtering high frequencies prevents collapse. Finally, these equations elucidate the unique role of the zero mode for the Neumann boundary conditions.

The spin-transfer effect is investigated for the vortex state of a magnetic nanodot. A spin current is shown to act similarly to an effective magnetic field perpendicular to the nanodot. Then a vortex with magnetization (polarity) parallel to the current polarization is energetically favorable. Following a simple energy analysis and using direct spin-lattice simulations, we predict the polarity switching of a vortex. For magnetic storage devices, an electric current is more effective to switch the polarity of a vortex in a nanodot than the magnetic field.

The propagation of extremely short pulses of an electromagnetic field (electromagnetic spikes) is considered in the framework of a model wherein the material medium is represented by anharmonic oscillators with cubic nonlinearities (Duffing model) and waves can propagate only in the right direction. The system of reduced Maxwell-Duffing equations admits two families of exact analytical solutions in the form of solitary waves. These are bright spikes propagating on a zero background, and bright and dark spikes propagating on a nonzero background. We find that these steady-state pulses are stable in terms of boundedness of the Hamiltonian. Direct simulations demonstrate that these pulses are very robust against perturbations. We find that a high-frequency modulated electromagnetic pulse evolves into a breather-like one. Conversely, a low frequency pulse transforms into a quasiharmonic wave.

The classical stimulated Raman scattering system describing the resonant interaction between two electromagnetic waves and a fast relaxing medium wave is studied by introducing a systematic perturbation approach in powers of the relaxation time. We separate amplitude and phase effects for these complex fields. The analysis of the former shows the existence of a stagnation distance after which monotonic energy transfer begins from one electromagnetic wave to the other, and this quantity is calculated. Concerning phase effects we give the conditions for the formation of a Raman spike from an initial fast and large phase jump in one of the waves. The spike evolution and width estimated from the reduced model agree with the results from numerical simulations of the original system.

We investigate the motion of a nonplanar vortex in a circular easy-plane magnet with a rotating in-plane magnetic field. Our numerical simulations of the Landau-Lifshitz equations show that the vortex tends to a circular limit trajectory, with an orbit frequency which is lower than the driving field frequency. To describe this we develop a new collective variable theory by introducing additional variables which account for the internal degrees of freedom of the vortex core, strongly coupled to the translational motion. We derive the evolution equations for these collective variables and find limit-cycle solutions whose characteristics are in qualitative agreement with the simulations of the many-spin system.