Search

The well-conditioned multi-product formula (MPF), proposed by [Low, Kliuchnikov, and Wiebe, 2019], is a simple high-order time-independent Hamiltonian simulation algorithm that implements a linear combination of standard product formulas of low order. While the MPF aims to simultaneously exploit commutator scaling among Hamiltonians and achieve near-optimal time and precision dependence, its lack of a rigorous error bound on the nested commutators renders its practical advantage ambiguous. In this work, we conduct a rigorous complexity analysis of the well-conditioned MPF, demonstrating explicit commutator scaling and near-optimal time and precision dependence at the same time. Using our improved complexity analysis, we present several applications of practical interest where the MPF based on a second-order product formula can achieve a polynomial speedup in both system size and evolution time, as well as an exponential speedup in precision, compared to second-order and even higher-order product formulas. Compared to post-Trotter methods, the MPF based on a second-order product formula can achieve polynomially better scaling in system size, with only poly-logarithmic overhead in evolution time and precision., Comment: 28 pages

For the problem of the non-isentropic compressible Euler Equations coupled with a nonlinear Poisson equation with the electric potential satisfying the Dirichlet boundary condition in three spatial dimensions with a general free boundary not restricting to a graph, we identify suitable stability conditions on the electric potential and the pressure under which we obtain a priori estimates on the Sobolev norms of the fluid and electric variables and bounds for geometric quantities of free surface. The stability conditions in this case for a general variable entropy are that the outer normal derivative of the electric potential is positive on the free surface, whereas that on the pressure is negative. In the isentropic case, the stability condition reduces to a single one, the outer normal derivative of the difference of the enthalpy and the electric potential is negative on the free surface. For the free boundary problem of the non-isentropic compressible Euler equations with variable entropy without coupling with the nonlinear Poisson equation, the corresponding higher-order estimates are also obtained under the Taylor sign condition. It is also found that one less derivative is needed to close the energy estimates for the problem for the non-isentropic compressible Euler Equations coupled with a nonlinear Poisson equation when the electric potential satisfies the Dirichlet boundary condition under the stability conditions on the electric potential and the pressure, compared with the problem of the non-isentropic compressible Euler equations.

High-dimensional fractional reaction-diffusion equations have numerous applications in the fields of biology, chemistry, and physics, and exhibit a range of rich phenomena. While classical algorithms have an exponential complexity in the spatial dimension, a quantum computer can produce a quantum state that encodes the solution with only polynomial complexity, provided that suitable input access is available. In this work, we investigate efficient quantum algorithms for linear and nonlinear fractional reaction-diffusion equations with periodic boundary conditions. For linear equations, we analyze and compare the complexity of various methods, including the second-order Trotter formula, time-marching method, and truncated Dyson series method. We also present a novel algorithm that combines the linear combination of Hamiltonian simulation technique with the interaction picture formalism, resulting in optimal scaling in the spatial dimension. For nonlinear equations, we employ the Carleman linearization method and propose a block-encoding version that is appropriate for the dense matrices that arise from the spatial discretization of fractional reaction-diffusion equations.

In a previous paper [Jayanti, P.C., Trivisa, K. Local Existence of Solutions to a Navier-Stokes-Nonlinear-Schr\"odinger Model of Superfluidity. J. Math. Fluid Mech. 24, 46 (2022)], the authors proved the existence of local-in-time weak solutions to a model of superfluidity. The system of governing equations was derived by Pitaevskii in 1959 and couples the nonlinear Schr\"odinger equation (NLS) and the Navier-Stokes equations (NSE). In this article, we prove a weak-strong type uniqueness theorem for these weak solutions. Only some of their regularity properties are used, allowing room for improved existence theorems in the future, with compatible uniqueness results., Comment: 20 pages. Major revisions from previous versions including improved estimates and the removal of assumptions (small data/low energy/short time) for uniqueness

In a 1959 paper by Pitaevskii, a macroscopic model of superfluidity was derived from first principles, to describe the interacting dynamics between the superfluid and normal fluid phases of Helium-4. The model couples two of the most fundamental PDEs in mathematics: the nonlinear Schr\"odinger equation (NLS) and the Navier-Stokes equations (NSE). In this article, we show the local existence of weak solutions to this system (in a smooth bounded domain in 3D), by deriving the required a priori estimates. (We will also establish an energy inequality obeyed by the weak solutions constructed in Kim's 1987 paper for the incompressible, inhomogeneous NSE.) To the best of our knowledge, this is the first rigorous mathematical analysis of a bidirectionally coupled system of the NLS and NSE., Comment: 48 pages; added some references in the Introduction section

Nonlinear differential equations model diverse phenomena but are notoriously difficult to solve. While there has been extensive previous work on efficient quantum algorithms for linear differential equations, the linearity of quantum mechanics has limited analogous progress for the nonlinear case. Despite this obstacle, we develop a quantum algorithm for dissipative quadratic $n$-dimensional ordinary differential equations. Assuming $R < 1$, where $R$ is a parameter characterizing the ratio of the nonlinearity and forcing to the linear dissipation, this algorithm has complexity $T^2 q~\mathrm{poly}(\log T, \log n, \log 1/\epsilon)/\epsilon$, where $T$ is the evolution time, $\epsilon$ is the allowed error, and $q$ measures decay of the solution. This is an exponential improvement over the best previous quantum algorithms, whose complexity is exponential in $T$. While exponential decay precludes efficiency, driven equations can avoid this issue despite the presence of dissipation. Our algorithm uses the method of Carleman linearization, for which we give a novel convergence theorem. This method maps a system of nonlinear differential equations to an infinite-dimensional system of linear differential equations, which we discretize, truncate, and solve using the forward Euler method and the quantum linear system algorithm. We also provide a lower bound on the worst-case complexity of quantum algorithms for general quadratic differential equations, showing that the problem is intractable for $R \ge \sqrt{2}$. Finally, we discuss potential applications, showing that the $R < 1$ condition can be satisfied in realistic epidemiological models and giving numerical evidence that the method may describe a model of fluid dynamics even for larger values of $R$., Comment: 36 pages, 1 figure. Published in PNAS

We investigate the global existence of weak solutions to a free boundary problem governing the evolution of finitely extensible bead-spring chains in dilute polymers. The free boundary in the present context is defined with regard to a density threshold of \r{ho} = 1, below which the fluid is modeled as compressible and above which the fluid is modeled as incompressible. The present article focuses on the physically relevant case in which the viscosity coefficients present in the system depend on the polymer number density, extending the earlier work [8]. We construct the weak solutions of the free boundary problem by perform ing the asymptotic limit as the adiabatic exponent \gamma goes to \infty for the macroscopic model introduced by Feireisl, Lu and Suli in [10] (see also [6]). The weak sequential stability of the family of dissipative (finite energy) weak solutions to the free boundary problem is also established., Comment: 24 pages

The Hall-Vinen-Bekharevich-Khalatnikov (HVBK) equations are a macroscopic model of superfluidity at non-zero temperatures. For smooth, compactly supported data, we prove the global well-posedness of strong solutions to these equations in $\mathbb{R}^2$, in the incompressible and isothermal case. The proof utilises a contraction mapping argument to establish local well-posedness for high-regularity data, following which we demonstrate global regularity using an analogue of the Beale-Kato-Majda criterion in this context. In the appendix, we address the sufficient conditions on a 2D vorticity field, in order to have a finite kinetic energy., Comment: Typo fixed in equation (13): missing absolute value added

We investigate the global existence of weak solutions to a free boundary problem governing the evolution of finitely extensible bead-spring chains in dilute polymers. We construct weak solutions of the two-phase model by performing the asymptotic limit as the adiabatic exponent $\gamma$ goes to $\infty$ for a macroscopic model which arises from the kinetic theory of dilute solutions of nonhomogeneous polymeric liquids. In this context the polymeric molecules are idealized as bead-spring chains with finitely extensible nonlinear elastic (FENE) type spring potentials. This class of models involves the unsteady, compressible, isentropic, isothermal Navier-Stokes system in a bounded domain $\Omega$ in $\mathbb{R}^d,$ $d=2, 3$ coupled with a Fokker-Planck-Smoluchowski-type diffusion equation (cf. Barrett and S\"{u}li [4], [5], [9]). The convergence of these solutions, up to a subsequence, to the free-boundary problem is established using weak convergence methods, compactness arguments which rely on the monotonicity properties of certain quantities in the spirit of [19].

This article establishes the global existence of weak solutions to a model proposed by Rosensweig (Rosensweig, Ferrohydrodynamics (1985)) for the dynamics of ferrofluids. The system is expressed by the conservation of linear momentum, the incompressibility condition, the conservation of angular momentum, and the evolution of the magnetization. The existence proof is inspired by the DiPerna-Lions theory of renormalized solutions. In addition, the rigorous relaxation limit of the equations of ferrohydrodynamics towards the quasi-equilibrium is investigated. The proof relies on the relative entropy method, which involves constructing a suitable functional, analyzing its time evolution and obtaining convergence results for the sequence of approximating solutions.

If two probability density functions (PDFs) have values for their first $n$ moments which are quite close to each other (upper bounds of their differences are known), can it be expected that the PDFs themselves are very similar? Shown below is an algorithm to quantitatively estimate this "similarity" between the given PDFs, depending on how many moments one has information about. This method involves the concept of functions behaving "similarly" at certain "length scales", which is also precisely defined. This technique could find use in data analysis, to compare a data set with a PDF or another data set, without having to fit a functional form to the data.

We investigate the long-time behavior of solutions to a stochastically forced one-dimensional Navier-Stokes system, describing the motion of a compressible viscous fluid, in the case of linear pressure law. We prove existence of an invariant measure for the Markov process generated by strong solutions. We overcome the difficulties of working with non-Feller Markov semigroups on non-complete metric spaces by generalizing the classical Krylov-Bogoliubov method, and by providing suitable polynomial and exponential moment bounds on the solution, together with pathwise estimates., Comment: 23 pages

We investigate the dynamics of a nonlinear system modeling tumor growth with drug application. The tumor is viewed as a mixture consisting of proliferating, quiescent and dead cells as well as a nutrient in the presence of a drug. The system is given by a multi-phase flow model: the densities of the different cells are governed by a set of transport equations, the density of the nutrient and the density of the drug are governed by rather general diffusion equations, while the velocity of the tumor is given by Darcy's equation. The domain occupied by the tumor in this setting is a growing continuum $\Omega$ with boundary $\partial \Omega$ both of which evolve in time. Global-in-time weak solutions are obtained using an approach based on the vanishing viscosity of the Brinkman's regularization. Both the solutions and the domain are rather general, no symmetry assumption is required and the result holds for large initial data., Comment: arXiv admin note: substantial text overlap with arXiv:1408.4794

We investigate the stability and global existence of weak solutions to a free boundary problem governing the evolution of polymeric fluids. We construct weak solutions of the two-phase model by performing the asymptotic limit of a macroscopic model governing the suspensions of rod-like molecules (known as Doi-Model) in compressible fluids as the adiabatic exponent $\gamma$ goes to $\infty.$ The convergence of these solutions, up to a subsequence, to the free-boundary problem is established using techniques in the spirit of Lions and Masmoudi \cite{LionsMasmoudi-1999}.

We investigate the evolution of tumor growth relying on a nonlinear model of partial differential equations which incorporates mechanical laws for tissue compression combined with rules for nutrients availability and drug application. Rigorous analysis and simulations are presented which show the role of nutrient and drug application in the progression of tumors. We construct an explicit convergent numerical scheme to approximate solutions of the nonlinear system of partial differential equations. Extensive numerical tests show that solutions exhibit a necrotic core when the nutrient level falls below a critical level in accordance with medical observations. The same numerical experiment is performed in the case of drug application for the purpose of comparison. Depending on the balance between nutrient and drug both shrinkage and growth of tumors can occur. The role of inhomogeneous boundary conditions, vascularization and anisotropies in the development of tumor shape irregularities are discussed., Comment: arXiv admin note: substantial text overlap with arXiv:1504.05982

Mechanical models for tumor growth have been used extensively in recent years for the analysis of medical observations and for the prediction of cancer evolution based on imaging analysis. This work deals with the numerical approximation of a mechanical model for tumor growth and the analysis of its dynamics. The system under investigation is given by a multi-phase flow model: The densities of the different cells are governed by a transport equation for the evolution of tumor cells, whereas the velocity field is given by a Brinkman regularization of the classical Darcy's law. An efficient finite difference scheme is proposed and shown to converge to a weak solution of the system. Our approach relies on convergence and compactness arguments in the spirit of Lions (Mathematical Topics in Fluid Dynamics, 1998).

We investigate the dynamics of a nonlinear system modeling tumor growth with drug application. The tumor is viewed as a mixture consisting of proliferating, quiescent and dead cells as well as a nutrient in the presence of a drug. The system is given by a multi-phase flow model: the densities of the different cells are governed by a set of transport equations, the density of the nutrient and the density of the drug are governed by rather general diffusion equations, while the velocity of the tumor is given by Brinkman's equation. The domain occupied by the tumor in this setting is a growing continuum $\Omega$ with boundary $\partial \Omega$ both of which evolve in time. Global-in-time weak solutions are obtained using an approach based on penalization of the boundary behavior, diffusion and viscosity in the weak formulation. Both the solutions and the domain are rather general, no symmetry assumption is required and the result holds for large initial data. This article is part of a research program whose aim is the investigation of the effect of drug application in tumor growth., Comment: 22 pages, 1 figure. arXiv admin note: substantial text overlap with arXiv:1408.4606; and text overlap with arXiv:1203.1215 by other authors

We investigate the dynamics of a nonlinear model for tumor growth within a cellular medium. In this setting the "tumor" is viewed as a multiphase flow consisting of cancerous cells in either proliferating phase or quiescent phase and a collection of cells accounting for the "waste" and/or dead cells in the presence of a nutrient. Here, the tumor is thought of as a growing continuum $\Omega$ with boundary $\partial \Omega$ both of which evolve in time. The key characteristic of the present model is that the total density of cancerous cells is allowed to vary, which is often the case within cellular media. We refer the reader to the articles \cite{Enault-2010}, \cite{LiLowengrub-2013} where compressible type tumor growth models are investigated. Global-in-time weak solutions are obtained using an approach based on penalization of the boundary behavior, diffusion, viscosity and pressure in the weak formulation, as well as convergence and compactness arguments in the spirit of Lions \cite{Lions-1998} (see also \cite{Feireisl-book, DT-MixedModel-2013})., Comment: 30 pages, 1 figure. arXiv admin note: text overlap with arXiv:1203.1215 by other authors

This article deals with relaxation approximations of nonlinear systems of hyperbolic balance laws. We introduce a class of relaxation schemes and establish their stability and convergence to the solution of hyperbolic balance laws before the formation of shocks, provided that we are within the framework of the compensated compactness method. Our analysis treats systems of hyperbolic balance laws with source terms satisfying a special mechanism which induces weak dissipation in the spirit of Dafermos [C.M. Dafermos J. Hyp. Diff. Equations, 3, 505-527, 2006], as well as hyperbolic balance laws with more general source terms. The rate of convergence of the relaxation system to a solution of the balance laws in the smooth regime is established. Our work follows in spirit the analysis presented in [S. Jin, X. Xin, Comm. Pure. Appl. Math. (1995), 48] and [Ch. Arvanitis, Ch. Makridakis, and A.E. Tzavaras, SIAM J. on Num. Anal. (2005), 42-4] for systems of hyperbolic conservation laws without source terms., Comment: 31 pages

We present a general framework for the approximation of systems of hyperbolic balance laws. The novelty of the analysis lies in the construction of suitable relaxation systems and the derivation of a delicate estimate on the relative entropy. We provide a direct proof of convergence in the smooth regime for a wide class of physical systems. We present results for systems arising in materials science, where the presence of source terms presents a number of additional challenges and requires delicate treatment. Our analysis is in the spirit of the framework introduced by Tzavaras [A. Tzavaras, Commun. Math. Sci., 3-2, 2005] for systems of hyperbolic conservation laws.

The hydrodynamic limit of a kinetic Cucker-Smale model is investigated. In addition to the free-transport of individuals and the Cucker-Smale alignment operator, the model under consideration includes a strong local alignment term. This term was recently derived as the singular limit of an alignment operator due to Motsch and Tadmor. The model is enhanced with the addition of noise and a confinement potential. The objective of this work is the rigorous investigation of the singular limit corresponding to strong noise and strong local alignment. The proof relies on the relative entropy method and entropy inequalities which yield the appropriate convergence results. The resulting limiting system is an Euler-type flocking system., Comment: 23 pages

In two recent papers the authors study the existence of weak solutions and the hydrodynamic limit of kinetic flocking equations with strong local alignment. The introduction of a strong local alignment term to model flocking behavior was formally motivated in these papers as a limiting case of an alignment term proposed by Motsch and Tadmor. In this paper, we rigorously justify this limit, and show that the considered equation is indeed a limit of the Motsch-Tadmor model when the radius of interaction goes to zero. The analysis involves velocity averaging lemmas and several $L^p$ estimates., Comment: 10 pages

We establish the global existence of weak solutions to a class of kinetic flocking equations. The models under consideration include the kinetic Cucker-Smale equation with possibly non-symmetric flocking potential, the Cucker-Smale equation with additional strong local alignment, and a newly proposed model by Motsch and Tadmor. The main tools employed in the analysis are the velocity averaging lemma and the Schauder fixed point theorem along with various integral bounds., Comment: 27 pages

This article deals with the issues of global-in-time existence and asymptotic analysis of a fluid-particle interaction model in the so-called bubbling regime. The mixture occupies the physical space $\Omega \subset \mathbb{R}^3$ which may be unbounded. The system under investigation describes the evolution of particles dispersed in a viscous compressible fluid and is expressed by the conservation of fluid mass, the balance of momentum and the balance of particle density often referred as the Smoluchowski equation. The coupling between the dispersed and dense phases is obtained through the drag forces that the fluid and the particles exert mutually by action-reaction principle. We show that solutions exist globally in time under reasonable physical assumptions on the initial data, the physical domain, and the external potential. Furthermore, we prove the large-time stabilization of the system towards a unique stationary state fully determined by the masses of the initial density of particles and fluid and the external potential., Comment: 35 pages