First-order Evolution Research Articles

Efficient higher-order matrix product operators for time evolution

We introduce a systematic construction of higher-order matrix product operator (MPO) approximations of the time evolution operator for generic (short and long range) one-dimensional Hamiltonians. We demonstrate the utility of our construction, by showing an order of magnitude speedup in simulation cost compared to conventional first-order MPO time evolution schemes.

An Analytical Model for Tropical Cyclone Outer-Size Expansion on the f Plane

Abstract Tropical cyclones are known to expand to an equilibrium size on the f plane, but the expansion process is not understood. In this study, an analytical model for tropical cyclone outer-size expansion on the f plane is proposed. Conceptually, the storm expands because the imbalance between latent heating and radiative cooling drives a lateral inflow that imports absolute vorticity. Volume-integrated latent heating increases more slowly with size than radiative cooling, and hence, the storm expands toward an equilibrium size. The predicted expansion rate is given by the ratio of the difference in size from its equilibrium value rt,eq to an environmentally determined time scale τrt of 10–15 days. The model is fully predictive if given a constant rt,eq, which can also be estimated environmentally. The model successfully captures the first-order size evolution across a range of numerical simulation experiments in which the potential intensity and f are varied. The model predictions of the dependencies of lateral inflow velocity and expansion rate on latent heating rate are also compared well with numerical simulations. This model provides a useful foundation for understanding storm size dynamics in nature.

The Presence of Chaos in a Viscoelastic Harmonically Forced Von Mises Truss

Abstract This work investigates how viscoelasticity affects the dynamic behavior of a lumped-parameter model of a bistable von Mises truss. The system is controlled by a linear first-order equation and a second-order nonlinear Duffing equation with a quadratic nonlinearity that governs mechanical behavior. The second-order equation controls mechanical oscillations, while the linear first-order equation controls viscoelastic force evolution. Combined, the two equations form a third-order jerk equation that controls system dynamics. Viscoelasticity adds time scales and degrees-of-freedom to material behavior, distinguishing it from viscosity-only systems. Due to harmonic excitation, the system exhibits varied dynamic responses, from periodic to quasi-periodic to chaotic. We explore the dynamics of a harmonically forced von Mises truss with viscous damping to address this purpose. We demonstrate this system's rich dynamic behavior due to driving amplitude changes. This helps explain viscoelastic system behavior. A viscoelastic unit replaces the viscous damper, and we show that, although viscous damping merely changes how fast the trajectory decays to an attractor, viscoelasticity modifies both the energy landscape and the rate of decay. In a conventional linear solid model, three viscoelastic parameters control the system's behavior instead of one, as in pure viscous damping. This adds degrees-of-freedom that affect system dynamics. We present the parameter space for chaotic behavior and the shift from regular to irregular motion. Finally, Melnikov's criteria identify the regular-chaotic threshold. The system's viscous and elastic components affect the chaotic threshold amplitude

Micro- and Macroevolution: A Continuum or Two Distinct Types of Change?

How microevolution and macroevolution are related is one of the major unanswered questions in evolutionary biology. The most prevalent view is that microevolution and macroevolution are part of a continuum of one type of change and that macroevolution is the cumulative result of microevolution. Mathematics, however, distinguishes two fundamentally different, singular types of change: change of a vector in its parameters versus its dimensions. This mathematical distinction may help to articulate the concept of evolution by distinction of two fundamentally different types of evolution: the change of the state vector of an organism in 1) its parameters (= ‘first-order evolution’) and 2) its dimensions (= ‘second-order evolution’). This distinction can be operationalized by identifying genes and regulatory elements in the nucleotide code of an organism as dimensions and the level of expression as parameters of its state vector. This operationalization allows us to substitute the phenotype-based analysis of evolution with a genotype-based analysis and draws attention to the molecular mechanisms that change the parameters or the dimensions of the state vector, respectively. We illustrate the distinction between first- and second-order evolution with a simulation of the adaptive dynamics of a population of digital amoebae. Our genotype-based systems approach reveals that micro- and macroevolution are largely similar to first- and second-order evolution respectively, and are not a continuum of change.

Raznostnye skhemy dekompozitsii na osnove rasshchepleniya resheniya i operatora zadachi

Domain decomposition methods are used for the approximate solution of boundary value problems for partial differential equations on parallel computing systems. The specifics of nonstationary problems is most completely taken into account when using noniterative domain decomposition schemes. Regionally additive schemes are constructed on the basis of various classes of splitting schemes. A new class of domain decomposition schemes with an additive representation of the solution on a new time level is distinguished that is based on splitting the domain into subdomains based on a partition of unity. An example of the Cauchy problem for first-order evolution equations with a positive self-adjoint operator in a finite-dimensional Hilbert space is considered. Unconditionally stable two- and three-level splitting schemes are constructed for the corresponding system of equations.

Shear imposed falling liquid films on a slippery substrate with Marangoni effects: Effect of odd viscosity

We investigate the behavior of a thin fluid with disrupted time-reversal symmetry on a uniformly heated inclined surface under external shear stress using modified Navier–Stokes equations, an energy conservation equation, and incorporating a Navier slip condition. Critical conditions for instability onset are determined by a linear stability analysis within the Orr–Sommerfeld framework. We derive a first-order Benney-type evolution equation to study long-wave instabilities. We find slippery substrate, imposed shear stress along the flow direction, and Marangoni number consistently destabilize the flow, while odd viscosity and imposed counter-flow shear stabilize it. A weakly nonlinear analysis using multiple scales reveal distinct zones of instability. Marangoni number, slip length, odd viscosity, and imposed shear direction significantly impact stability and instability regions. Numerical simulations of the free surface evolution equation of a flow system under consideration provide clear evidence of the contributions of thermocapillary, slip length, odd viscosity, and imposed shear direction. Furthermore, our analysis of linear and weakly nonlinear stability, as well as our numerical simulations, exhibit remarkable consistency.

Accelerated dynamics with dry friction via time scaling and averaging of doubly nonlinear evolution equations

In a Hilbert framework, for convex differentiable optimization, we analyze the long-time behavior of inertial dynamics with dry friction. In classical approaches based on asymptotically vanishing viscous damping (in accordance with Nesterov’s method), the results are expressed in terms of rapid convergence of the values of the function to its minimum value. On the other hand, dry friction induces convergence towards an approximate minimizer, typically the system stops at x when a given threshold ‖∇f(x)‖≤r is satisfied. We obtain rapid convergence results in this direction. In our approach, we start from a doubly nonlinear first-order evolution equation involving two potentials: one is the differentiable function f to be minimized, which acts on the state of the system via its gradient, and the other is the nonsmooth potential dry friction φ(x)=r‖x‖ which acts on the velocity vector via its sub-differential. To highlight the central role played by ∇f(x), we also deal with the dual formulation of these dynamics, which have a Riemannian gradient structure. We then rely on the general acceleration method recently developed by Attouch, Boţ and Nguyen, which consists in applying the method of time scaling and then averaging to a continuous differential equation of the first order in time. We thus obtain fast convergence results for second order time-evolution systems involving dry friction, asymptotically vanishing viscous damping, and Hessian-driven damping in the implicit form.

Line operators in Chern-Simons-Matter theories and Bosonization in Three Dimensions II: Perturbative analysis and all-loop resummation

We study mesonic line operators in Chern-Simons theories with bosonic or fermionic matter in the fundamental representation. In this paper, we elaborate on the classification and properties of these operators using all loop resummation of large N perturbation theory. We show that these theories possess two conformal line operators in the fundamental representation. One is a stable renormalization group fixed point, while the other is unstable. They satisfy first-order chiral evolution equations, in which a smooth variation of the path is given by a factorized product of two mesonic line operators. The boundary operators on which the lines can end are classified by their conformal dimension and transverse spin, which we compute explicitly at finite ’t Hooft coupling. We match the operators in the bosonic and fermionic theories. Finally, we extend our findings to the mass deformed theories and discover that the duality still holds true.

Tikhonov regularization of a nonhomogeneous first-order evolution equation

We consider a nonhomogeneous first-order evolution equation governed by a maximal monotone operator $ A $ in a Hilbert space in the presence of a Tikhonov regularization term. We study the existence and strong convergence of weak solutions to such systems. With boundedness conditions on the central path or on the solution trajectory, without assuming the zero set of $ A $ to be nonempty and with suitable assumptions on the Tikhonov regularization coefficient, we prove that the weak solutions to the system converge strongly to the element of least norm in the zero set of $ A $. As a consequence, we provide sufficient conditions where the boundedness of the weak solutions to the nonhomogeneous Tikhonov system is equivalent to the zero set of $ A $ to be nonempty. Our work is motivated by [4,17,22] and extends some results by those authors.

A coupled problem described by time-dependent subdifferential operator and non-convex perturbed sweeping process

The aim of the present paper is to state the well-posedness of a new dynamical system driven by a differential inclusion involving time-dependent subdifferential operators with integral perturbation and a non-convex perturbed sweeping process. The current study is motivated by recent works on integro-differential sweeping processes and dynamical systems coupled by differential inclusions (involving sweeping processes and maximal monotone operators). In our development, we proceed using a discretization method, combining tools of first-order evolution problems of subdifferential type and those governed by the normal cone of $ r $-prox regular moving sets. We establish a new existence result to the class of dynamical systems under consideration, in the context of infinite dimensional Hilbert setting. This allows us to deal with a Bolza-type problem in optimal control theory.

Shallow water waves generated by a floating object: A control theoretical perspective

<p style='text-indent:20px;'>We consider a control system describing the interaction of water waves with a partially immersed rigid body constraint to move only in the vertical direction. The fluid is modeled by the shallow water equations. The control signal is a vertical force acting on the floating body. We first derive the full governing equations of the fluid-body system in a water tank and reformulate them as an initial boundary value problem of a first-order evolution system. We then linearize the equations around the equilibrium state and we study its well-posedness. Finally we focus on the reachability and stabilizability of the linear system. Our main result asserts that, provided that the floating body is situated in the middle of the tank, any symmetric waves with appropriate regularity can be obtained from the equilibrium state by an appropriate control force. This implies, in particular, that we can project this system on the subspace of states with appropriate symmetry properties to obtain a reduced system which is approximately controllable and strongly stabilizable. Note that, in general, this system is not controllable (even approximately).</p>

Predictive comparison and assessment of ANFIS and ANN, as efficient tools in modeling degradation of total petroleum hydrocarbon (TPH)

Composting is an effective organic waste management practice with a sustainable environmental impact on soil. In most cases, composting materials are agricultural residues constituting a huge environmental challenge. This study investigates the composting biostimulation of total petroleum hydrocarbon (TPH) contaminated soil using organic compost prepared from palm oil mill effluent (CPOME). The performance of CPOME in TPH degradation was evaluated using the first-order kinetic model and CO₂ evolution. The artificial neural network (ANN) and adaptive neuro-fuzzy inference system (ANFIS) modelling approaches were further employed to optimize the experimental response. The experimental designs for ANFIS and ANN modelling were done using the central composite design (CCD). Statistical parameters such as RMSE, MSE, and MAE were used to evaluate model performance and comparisons. The first-order kinetic model and CO₂ released showed that TPH degradation favours contaminated soil with 40 g/L initial TPH concentration. The first-order kinetic model observed a decrease in the constant degradation rate with increasing TPH concentration. The CO₂ evolution, which measures microbial respiration/activity, also showed a decline in released CO₂ with increasing TPH concentration. Lower values of RMSE, MSE, and MAE were obtained for the ANFIS model compared with the ANN. The R² values between the experimental and model-predicted responses showed that the ANFIS (0.9349) model predictions fitted adequately compared to the ANN (0.8862). Finally, comparison and model performances showed that the ANFIS model prediction capability was considerably better than the ANN. • The compost from palm oil mill effluent proved effective in TPH degradation • ANFIS predicted TPH degradation was better than the ANN predicted • Notable organic materials in the compost includes polysaccharides, nitrogen and lignin

Approximate solution of the Cauchy problem for a first-order integrodifferential equation with solution derivative memory

We consider the Cauchy problem for a first-order evolution equation with memory in a finite-dimensional Hilbert space when the integral term is related to the time derivative of the solution. The main problems of the approximate solution of such nonlocal problems are due to the necessity to work with the approximate solution for all previous time moments. We propose a transformation of the first-order integrodifferential equation to a system of local evolutionary equations. We use the approach known in the theory of Volterra integral equations with an approximation of the difference kernel by the sum of exponents. We formulate a local problem for a weakly coupled system of equations with additional ordinary differential equations. We have given estimates of the stability of the solution by initial data and the right-hand side for the solution of the corresponding Cauchy problem. The primary attention is paid to constructing and investigating the stability of two-level difference schemes, which are convenient for computational implementation. The numerical solution of a two-dimensional model problem for the evolution equation of the first order, when the Laplace operator conditions the dependence on spatial variables, is presented.

On constraint preservation and strong hyperbolicity

We use partial differential equations (PDEs) to describe physical systems. In general, these equations include evolution and constraint equations. One method used to find solutions to these equations is the free-evolution approach, which consists in obtaining the solutions of the entire system by solving only the evolution equations. Certainly, this is valid only when the chosen initial data satisfies the constraints and the constraints are preserved in the evolution. In this paper, we establish the sufficient conditions required for the PDEs of the system to guarantee the constraint preservation. This is achieved by considering quasi-linear first-order PDEs, assuming the sufficient condition and deriving strongly hyperbolic first-order partial differential evolution equations for the constraints. We show that, in general, these constraint evolution equations correspond to a family of equations parametrized by a set of free parameters. We also explain how these parameters fix the propagation velocities of the constraints. As application examples of this framework, we study the constraint conservation of the Maxwell electrodynamics and the wave equations in arbitrary space–times. We conclude that the constraint evolution equations are unique in the Maxwell case and a family in the wave equation case.

Line Operators in Chern-Simons–Matter Theories and Bosonization in Three Dimensions

We study Chern-Simons theories at large N with either bosonic or fermionic matter in the fundamental representation. The most fundamental operators in these theories are mesonic line operators, the simplest example being Wilson lines ending on fundamentals. We classify the conformal line operators along an arbitrary smooth path as well as the spectrum of conformal dimensions and transverse spins of their boundary operators at finite 't Hooft coupling. These line operators are shown to satisfy first-order chiral evolution equations, in which a smooth variation of the path is given by a factorized product of two line operators. We argue that this equation together with the spectrum of boundary operators are sufficient to uniquely determine the expectation values of these operators. We demonstrate this by bootstrapping the two-point function of the displacement operator on a straight line. We show that the line operators in the theory of bosons and the theory of fermions satisfy the same evolution equation and have the same spectrum of boundary operators.

On Average Losses of Low-Frequency Sound in a Two-Dimensional Shallow-Water Random Waveguide

For a low-frequency sound signal propagating in a two-dimensionally inhomogeneous shallow-water waveguide, the influence of random bathymetry (fluctuating bottom boundary) was considered based on the local-mode approach and statistical modeling using first-order evolution equations. The study was carried out in shallow sea conditions corresponding to the coastal waveguides of the Russian Arctic seas. Here, a feature was the presence of an almost homogeneous water layer with various characteristics of seabed sediments. To describe the latter, a random model of the impedance was adopted. For the conditions of a strongly penetrable bottom boundary, on average, the calculations predicted adequate weak effects of bathymetry fluctuations on the average sound intensity compared to the effect of fluctuations in the sediment parameters and volumetric random inhomogeneities of the water column. In addition, it was shown that, in terms of statistics, the roughness of the bottom boundary perturbed the average sound intensity in a shallow-water waveguide differently than volumetric fluctuations in the speed of sound. The dependence of the statistical effects (the first and second moments of the signal intensity) on the parameters of the waveguide and the frequency range was studied. As a result of numerical modeling, comparative quantitative estimates of the influence of both the random roughness of the bottom interface and fluctuations of bottom sediment parameters on the average losses of the propagating signal, not presented in the literature, were obtained.

Numerical solution of the heat conduction problem with memory

It is necessary to use more general models than the classical Fourier heat conduction law to describe small-scale thermal conduction processes. The effects of heat flow memory and heat capacity memory (internal energy) in solids are considered in first-order integrodifferential evolutionary equations with difference-type kernels. The main difficulties in applying such nonlocal in-time mathematical models are associated with the need to work with a solution throughout the entire history of the process. The paper develops an approach to transforming a nonlocal problem into a computationally simpler local problem for a system of first-order evolution equations. Such a transition is applicable for heat conduction problems with memory if the relaxation functions of the heat flux and heat capacity are represented as a sum of exponentials. The correctness of the auxiliary linear problem is ensured by the obtained estimates of the stability of the solution concerning the initial data and the right-hand side in the corresponding Hilbert spaces. The study's main result is to prove the unconditional stability of the proposed two-level scheme with weights for the evolutionary system of equations for modeling heat conduction in solid media with memory. In this case, finding an approximate solution on a new level in time is not more complicated than the classical heat equation. The numerical solution of a model one-dimensional in space heat conduction problem with memory effects is presented.

Weak Gravitation in the 4+1 Formalism

The 4+1 formalism in general relativity (GR) prescribes field equations for the spacetime metric γμνx,τ which is local in the spacetime coordinates x and evolves according to an external “worldtime” τ. This formalism extends to GR the Stueckelberg Horwitz Piron (SHP) framework, developed to address the various issues known as the problem of time as they appear in electrodynamics. SHP field theories exhibit a formal 5D symmetry on (x,τ) that is strategically broken to 4+1 representations of the Lorentz group, resulting in a manifestly covariant canonical formalism describing the τ-evolution of spacetime structures as an initial value problem. Einstein equations for γμνx,τ are found by constructing a 5D pseudo-manifold (combining 4D geometry and τ-dynamics) and performing the natural foliation under broken 5D symmetry. This paper discusses weak gravitation in the 4+1 formalism, demonstrating the natural decomposition of the field equations into first-order evolution equations for the unconstrained 4D metric, and the propagation of constraints associated with the Bianchi identity.

Quantum algorithms for quantum dynamics: A performance study on the spin-boson model

Quantum algorithms for quantum dynamics simulations are traditionally based on implementing a Trotter approximation of the time-evolution operator. This approach typically relies on deep circuits and is therefore hampered by the substantial limitations of available noisy and near-term quantum hardware. On the other hand, variational quantum algorithms (VQAs) have become an indispensable alternative, enabling small-scale simulations on present-day hardware. However, despite the recent development of VQAs for quantum dynamics, a detailed assessment of their efficiency and scalability is yet to be presented. To fill this gap, we applied a VQA based on McLachlan's principle to simulate the dynamics of a spin-boson model subject to varying levels of realistic hardware noise as well as in different physical regimes, and discuss the algorithm's accuracy and scaling behavior as a function of system size. We observe a good performance of the variational approach used in combination with a general, physically motivated wave function ansatz, and compare it to the conventional first-order Trotter evolution. Finally, based on this, we make scaling predictions for the simulation of a classically intractable system. We show that, despite providing a clear reduction of quantum gate cost, the variational method in its current implementation is unlikely to lead to a quantum advantage for the solution of time-dependent problems.

Existence Theorem for Noninstantaneous Impulsive Evolution Equations

In this note, the variational form of the classical Lax–Milgram theorem is used for the divulgence of variational structure of the first-order noninstantaneous impulsive linear evolution equation. The existence and uniqueness of the weak solution of the problem is obtained. In future, this constructive theory can be used for the corresponding semilinear problems.