Exact Conservation Laws Research Articles

A Legendre–Fourier spectral method with exact conservation laws for the Vlasov–Poisson system

We present the design and implementation of an L2-stable spectral method for the discretization of the Vlasov- Poisson model of a collisionless plasma in one space and velocity dimension. The velocity and space dependence of the Vlasov equation are resolved through a truncated spectral expansion based on Legendre and Fourier basis functions, respectively. The Poisson equation, which is coupled to the Vlasov equation, is also resolved through a Fourier expansion. The resulting system of ordinary differential equation is discretized by the implicit second-order accurate Crank-Nicolson time discretization. The non-linear dependence between the Vlasov and Poisson equations is iteratively solved at any time cycle by a Jacobian-Free Newton-Krylov method. In this work we analyze the structure of the main conservation laws of the resulting Legendre-Fourier model, e.g., mass, momentum, and energy, and prove that they are exactly satisfied in the semi-discrete and discrete setting. The L2-stability of the method is ensured by discretizing the boundary conditions of the distribution function at the boundaries of the velocity domain by a suitable penalty term. The impact of the penalty term on the conservation properties is investigated theoretically and numerically. An implementation of the penalty term that does not affect the conservation of mass, momentum and energy, is also proposed and studied. A collisional term is introduced in the discrete model to control the filamentation effect, but does not affect the conservation properties of the system. Numerical results on a set of standard test problems illustrate the performance of the method.

Exact solutions and conservation laws of the system of two-dimensional viscous Burgers equations

Fluid turbulence is one of the phenomena that has been studied extensively for many decades. Due to its huge practical importance in fluid dynamics, various models have been developed to capture both the indispensable physical quality and the mathematical structure of turbulent fluid flow. Among the prominent equations used for gaining in-depth insight of fluid turbulence is the two-dimensional Burgers equations. Its solutions have been studied by researchers through various methods, most of which are numerical. Being a simplified form of the two-dimensional Navier–Stokes equations and its wide range of applicability in various fields of science and engineering, development of computationally efficient methods for the solution of the two-dimensional Burgers equations is still an active field of research. In this study, Lie symmetry method is used to perform detailed analysis on the system of two-dimensional Burgers equations. Optimal system of one-dimensional subalgebras up to conjugacy is derived and used to obtain distinct exact solutions. These solutions not only help in understanding the physical effects of the model problem but also, can serve as benchmarks for constructing algorithms and validation of numerical solutions of the system of Burgers equations under consideration at finite Reynolds numbers. Independent and nontrivial conserved vectors are also constructed.

Darboux transformation of the general Hirota equation: multisoliton solutions, breather solutions, and rogue wave solutions

In this paper, we investigate the exact solutions and conservation laws of a general Hirota equation. Firstly, the N-fold Darboux transformation of this equation is proposed. Then by choosing three kinds of seed solutions, the multisoliton solutions, breather solutions, and rogue wave solutions of the general Hirota equation are obtained based on the Darboux transformation. Finally, the conservation laws of this equation are derived by using its linear spectral problem. The results in this paper may be useful in the study of ultrashort optical solitons in optical fibers.

Another look at synchronized neutrino oscillations

In dense neutrino backgrounds present in supernovae and in the early Universe neutrino oscillations may exhibit complex collective phenomena, such as synchronized oscillations, bipolar oscillations and spectral splits and swaps. We consider in detail possible decoherence effects on the simplest of these phenomena – synchronized neutrino oscillations that can occur in a uniform and isotropic neutrino gas. We develop an exact formalism of spectral moments of the flavour spin vectors describing such a system and then apply it to find analytical approaches that allow one to study decoherence effects on its late-time evolution. This turns out to be possible in part due to the existence of the (previously unknown) exact conservation law satisfied by the quantities describing the considered neutrino system. Interpretation of the decoherence effects in terms of neutrino wave packet separation is also given, both in the adiabatic and non-adiabatic regimes of neutrino flavour evolution.

Exact Solutions and Conservation Laws for a Generalized Double Combined sinh–cosh–Gordon Equation

In this paper we study a generalized double combined sinh–cosh–Gordon equation, which appears in a wide range of physical applications. It admits geometric interpretation as the differential equation which determines time-like surfaces of constant positive curvature in the same spaces. We first compute the optimal system of one-dimensional subalgebras and then use it to obtain the optimal system of group-invariant solutions for the equation. We employ Lie group method along with the simplest equation method to investigate travelling wave solutions for this equation. In addition conservation laws are constructed using two different techniques, namely, the new conservation theorem and the multiplier method.

Lagrangian and Hamiltonian constraints for guiding-center Hamiltonian theories

A consistent guiding-center Hamiltonian theory is derived by Lie-transform perturbation method, with terms up to second order in magnetic-field nonuniformity. Consistency is demonstrated by showing that the guiding-center transformation presented here satisfies separate Jacobian and Lagrangian constraints that have not been explored before. A new first-order term appearing in the guiding-center phase-space Lagrangian is identified through a calculation of the guiding-center polarization. It is shown that this new polarization term also yields a simpler expression of the guiding-center toroidal canonical momentum, which satisfies an exact conservation law in axisymmetric magnetic geometries. Finally, an application of the guiding-center Lagrangian constraint on the guiding-center Hamiltonian yields a natural interpretation for its higher-order corrections.

Energetically consistent collisional gyrokinetics

We present a formulation of collisional gyrokinetic theory with exact conservation laws for energy and canonical toroidal momentum. Collisions are accounted for by a nonlinear gyrokinetic Landau operator. Gyroaveraging and linearization do not destroy the operator's conservation properties. Just as in ordinary kinetic theory, the conservation laws for collisional gyrokinetic theory are selected by the limiting collisionless gyrokinetic theory.

Multi-dimensional, fully-implicit, spectral method for the Vlasov–Maxwell equations with exact conservation laws in discrete form

A spectral method for the numerical solution of the multi-dimensional Vlasov–Maxwell equations is presented. The plasma distribution function is expanded in Fourier (for the spatial part) and Hermite (for the velocity part) basis functions, leading to a truncated system of ordinary differential equations for the expansion coefficients (moments) that is discretized with an implicit, second order accurate Crank–Nicolson time discretization. The discrete non-linear system is solved with a preconditioned Jacobian–Free Newton–Krylov method.It is shown analytically that the Fourier–Hermite method features exact conservation laws for total mass, momentum and energy in discrete form. Standard tests involving plasma waves and the whistler instability confirm the validity of the conservation laws numerically. The whistler instability test also shows that we can step over the fastest time scale in the system without incurring in numerical instabilities. Some preconditioning strategies are presented, showing that the number of linear iterations of the Krylov solver can be drastically reduced and a significant gain in performance can be obtained.

Irreversible energy flow in forced Vlasov dynamics

A recent paper [Phys. Plasmas 20, 032304 (2013)] considered the forced linear Vlasov equation as a model for the quasi-steady state of a single stable plasma wavenumber interacting with a bath of turbulent fluctuations. This approach gives some insight into possible energy flows without solving for nonlinear dynamics. The central result of the present work is that the forced linear Vlasov equation exhibits asymptotically zero (irreversible) dissipation to all orders under a detuning of the forcing frequency and the characteristic frequency associated with particle streaming. We first prove this by direct calculation, tracking energy flow in terms of certain exact conservation laws of the linear (collisionless) Vlasov equation. Then we analyze the steady-state solutions in detail using a weakly collisional Hermite-moment formulation, and compare with numerical solution. This leads to a detailed description of the Hermite energy spectrum, and a proof of no dissipation at all orders, complementing the collisionless Vlasov result.

Quasilocal conservation laws in XXZ spin-1/2 chains: Open, periodic and twisted boundary conditions

A continuous family of quasilocal exact conservation laws is constructed in the anisotropic Heisenberg (XXZ) spin-1/2 chain for periodic (or twisted) boundary conditions and for a set of commensurate anisotropies densely covering the entire easy plane interaction regime. All local conserved operators follow from the standard (Hermitian) transfer operator in fundamental representation (with auxiliary spin s=1/2), and are all even with respect to a spin flip operation. However, the quasilocal family is generated by differentiation of a non-Hermitian highest weight transfer operator with respect to a complex auxiliary spin representation parameter s and includes also operators of odd parity. For a finite chain with open boundaries the time derivatives of quasilocal operators are not strictly vanishing but result in operators localized near the boundaries of the chain. We show that a simple modification of the non-Hermitian transfer operator results in exactly conserved, but still quasilocal operators for periodic or generally twisted boundary conditions. As an application, we demonstrate that implementing the new exactly conserved operator family for estimating the high-temperature spin Drude weight results, in the thermodynamic limit, in exactly the same lower bound as for almost conserved family and open boundaries. Under the assumption that the bound is saturating (suggested by agreement with previous thermodynamic Bethe ansatz calculations) we propose a simple explicit construction of infinite time averages of local operators such as the spin current.

Exact Solutions and Conservation Laws of the Drinfel’d-Sokolov-Wilson System

We study the Drinfel'd-Sokolov-Wilson system, which was introduced as a model of water waves. Firstly we obtain exact solutions of this system using the(G′/G)-expansion method. In addition to exact solutions we also construct conservation laws for the underlying system using Noether's approach.

Exact solutions and conservation laws of a ( 3 + 1 ) -dimensional B-type Kadomtsev-Petviashvili equation

In this paper we study a -dimensional generalized B-type Kadomtsev-Petviashvili (BKP) equation. This equation is an extension of the well-known Kadomtsev-Petviashvili equation, which describes weakly dispersive and small amplitude waves propagating in quasi-two-dimensional media. We first obtain exact solutions of the BKP equation using the multiple-exp function and simplest equation methods. Furthermore, the conservation laws for the BKP equation are constructed by using the multiplier method.

New exact solutions and conservation laws of a coupled Kadomtsev–Petviashvili system

This paper obtains exact solutions of a new coupled Kadomtsev–Petviashvili system, which arises in the analysis of various problems in fluid mechanics, theoretical physics and many scientific applications. Lie symmetry method along with the (G′/G)-expansion method is employed to find the travelling wave solutions of the underlying system. In addition, we derive the conservation laws of the coupled Kadomtsev–Petviashvili system using the multiplier method.

Many-Body Localization in One Dimension as a Dynamical Renormalization Group Fixed Point

We formulate a dynamical real space renormalization group (RG) approach to describe the time evolution of a random spin-1/2 chain, or interacting fermions, initialized in a state with fixed particle positions. Within this approach we identify a many-body localized state of the chain as a dynamical infinite randomness fixed point. Near this fixed point our method becomes asymptotically exact, allowing analytic calculation of time dependent quantities. In particular, we explain the striking universal features in the growth of the entanglement seen in recent numerical simulations: unbounded logarithmic growth delayed by a time inversely proportional to the interaction strength. This is in striking contrast to the much slower entropy growth as loglogt found for noninteracting fermions with bond disorder. Nonetheless, even the interacting system does not thermalize in the long time limit. We attribute this to an infinite set of approximate integrals of motion revealed in the course of the RG flow, which become asymptotically exact conservation laws at the fixed point. Hence we identify the many-body localized state with an emergent generalized Gibbs ensemble.

Exact Solutions and Conservation Laws of a (2+1)-Dimensional Nonlinear KP-BBM Equation

In this paper, we study the two-dimensional nonlinear Kadomtsov-Petviashivilli-Benjamin-Bona-Mahony (KP-BBM) equation. This equation is the Benjamin-Bona-Mahony equation formulated in the KP sense. We first obtain exact solutions of this equation using the Lie group analysis and the simplest equation method. The solutions obtained are solitary waves. In addition, the conservation laws for the KP-BBM equation are constructed by using the multiplier method.

Exact Solutions and Conservation Laws of a Two-Dimensional Integrable Generalization of the Kaup-Kupershmidt Equation

We study a two-dimensional integrable generalization of the Kaup-Kupershmidt equation, which arises in various problems in mathematical physics. Exact solutions are obtained using the Lie symmetry method in conjunction with the extended tanh method and the extended Jacobi elliptic function method. In addition to exact solutions we also present conservation laws which are derived using the multiplier approach.

Symmetry reduction, exact group-invariant solutions and conservation laws of the Benjamin–Bona–Mahoney equation

We show that the Benjamin–Bona–Mahoney (BBM) equation with power law nonlinearity can be transformed by a point transformation to the combined KdV–mKdV equation, that is also known as the Gardner equation. We then study the combined KdV–mKdV equation from the Lie group-theoretic point of view. The Lie point symmetry generators of the combined KdV–mKdV equation are derived. We obtain symmetry reduction and a number of exact group-invariant solutions for the underlying equation using the Lie point symmetries of the equation. The conserved densities are also calculated for the BBM equation with dual nonlinearity by using the multiplier approach. Finally, the conserved quantities are computed using the one-soliton solution.

Thermodyamic Bounds on Drude Weights in Terms of Almost-conserved Quantities

We consider one-dimensional translationally invariant quantum spin (or fermionic) lattices and prove a Mazur-type inequality bounding the time-averaged thermodynamic limit of a finite-temperature expectation of a spatio-temporal autocorrelation function of a local observable in terms of quasi-local conservation laws with open boundary conditions. Namely, the commutator between the Hamiltonian and the conservation law of a finite chain may result in boundary terms only. No reference to techniques used in Suzuki’s proof of Mazur bound is made (which strictly applies only to finite-size systems with exact conservation laws), but Lieb-Robinson bounds and exponential clustering theorems of quasi-local C* quantum spin algebras are invoked instead. Our result has an important application in the transport theory of quantum spin chains, in particular it provides rigorous non-trivial examples of positive finite-temperature spin Drude weight in the anisotropic Heisenberg XXZ spin 1/2 chain (Prosen, in Phys Rev Lett 106:217206, 2011).

Exact solutions and conservation laws of a coupled integrable dispersionless system

In this paper we study the coupled integrable dispersionless system (CIDS), which arises in the analysis of several problems in applied mathematics and physics. Lie symmetry analysis is performed on CIDS and symmetry reductions and exact solutions with the aid of simplest equation method are obtained. In addition, the conservation laws of the CIDS are also derived using the multiplier (and homotopy) approach.

Exact solutions and conservation laws of Zakharov–Kuznetsov modified equal width equation with power law nonlinearity

This paper obtains solutions to the Zakharov–Kuznetsov modified equal width equation with power law nonlinearity. The Lie symmetry approach and the simplest equation method are used to obtain these solutions. Moreover, conservation laws are derived for the underlying equation by employing two approaches: the new conservation theorem and the multiplier method.