BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Florian Méhats appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
In this work we prove the nonlinear instability of inhomogeneous steady states solutions to the Hamiltonian Mean Field (HMF) model. We first study the linear instability of this model under a simple criterion by adapting the techniques developed in [19]. In a second part, we extend to the inhomogeneous case some techniques developed in [14, 17, 18] and prove a nonlinear instability result under the same criterion.
In this article, we propose a new numerical model for computation of the transport of electrons in a graphene device. The underlying quantum model for graphene is a massless Dirac equation, whose eigenvalues display a conical singularity responsible for non adiabatic transitions between the two modes. We first derive a kinetic model which takes the form of two Boltzmann equations coupled by a collision operator modeling the non-adiabatic transitions. This collision term includes a Landau-Zener transfer term and a jump operator whose presence is essential in order to ensure a good energy conservation during the transitions. We propose an algorithmic realization of the semi-group solving the kinetic model, by a particle method. We give analytic justification of the model and propose a series of numerical experiments studying the influences of the various sources of errors between the quantum and the kinetic models.
We apply the two-scale formulation approach to propose uniformly accurate (UA) schemes for solving the nonlinear Dirac equation in the nonrelativistic limit regime. The nonlinear Dirac equation involves two small scales $\varepsilon$ and $\varepsilon^2$ with $\varepsilon\to0$ in the nonrelativistic limit regime. The small parameter causes high oscillations in time which brings severe numerical burden for classical numerical methods. We transform our original problem as a two-scale formulation and present a general strategy to tackle a class of highly oscillatory problems involving the two small scales $\varepsilon$ and $\varepsilon^2$. Suitable initial data for the two-scale formulation is derived to bound the time derivatives of the augmented solution. Numerical schemes with uniform (with respect to $\varepsilon\in (0,1]$) spectral accuracy in space and uniform first order or second order accuracy in time are proposed. Numerical experiments are done to confirm the UA property.
We study dimension reduction for the three-dimensional Gross-Pitaevskii equation with a long-range and anisotropic dipole-dipole interaction modeling dipolar Bose-Einstein condensation in a strong interaction regime. The cases of disk shaped condensates (confinement from dimension three to dimension two) and cigar shaped condensates (confinement to dimension one) are analyzed. In both cases, the analysis combines averaging tools and semiclassical techniques. Asymptotic models are derived, with rates of convergence in terms of two small dimensionless parameters characterizing the strength of the confinement and the strength of the interaction between atoms.
This work is devoted to the analysis of the quantum Liouville-BGK equation. This equation arises in the work of Degond and Ringhofer on the derivation of quantum hydrodynamical models from first principles. Their theory consists in transposing to the quantum setting the closure strategy by entropy minimization used for kinetic equations. The starting point is the quantum Liouville-BGK equation, where the collision term is defined via a so-called quantum local equilibrium, defined as a minimizer of the quantum free energy under a local density constraint. We then address three related problems: we prove new results about the regularity of these quantum equilibria; we prove that the quantum Liouville-BGK equation admits a classical solution; and we investigate the long-time behavior of the solutions. The core of the proofs is based on a fine analysis of the properties of the minimizers of the free energy.
We consider the possibility that dark matter halos are described by the Fermi-Dirac distribution at finite temperature. This is the case if dark matter is a self-gravitating quantum gas made of massive neutrinos at statistical equilibrium. This is also the case if dark matter can be treated as a self-gravitating collisionless gas experiencing Lynden-Bell's type of violent relaxation. In order to avoid the infinite mass problem and carry out a rigorous stability analysis, we consider the fermionic King model. In this paper, we study the non-degenerate limit leading to the classical King model. This model was initially introduced to describe globular clusters. We propose to apply it also to large dark matter halos where quantum effects are negligible. We determine the caloric curve and study the thermodynamical stability of the different configurations. Equilibrium states exist only above a critical energy $E_c$ in the microcanonical ensemble and only above a critical temperature $T_c$ in the canonical ensemble. For $E<E_c$, the system undergoes a gravothermal catastrophe and, for $T<T_c$, it undergoes an isothermal collapse. We compute the profiles of density, circular velocity, and velocity dispersion. We compare the prediction of the classical King model to the observations of large dark matter halos. Because of collisions and evaporation, the central density increases while the slope of the halo density profile decreases until an instability takes place. We show that large dark matter halos are relatively well-described by the King model at, or close to, the point of marginal microcanonical stability. At that point, the King model generates a density profile that can be approximated by the modified Hubble profile. This profile has a flat core and decreases as $r^{-3}$ at large distances, like the observational Burkert profile. Less steep halos are unstable.
We discuss the nature of phase transitions in the fermionic King model which describes tidally truncated quantum self-gravitating systems. This distribution function takes into account the escape of high energy particles and has a finite mass. On the other hand, the Pauli exclusion principle puts an upper bound on the phase space density of the system and stabilizes it against gravitational collapse. As a result, there exists a statistical equilibrium state for any accessible values of energy and temperature. We plot the caloric curves and investigate the nature of phase transitions as a function of the degeneracy parameter in both microcanonical and canonical ensembles. We consider stable and metastable states and emphasize the importance of the latter for systems with long-range interactions. Phase transitions can take place between a "gaseous" phase unaffected by quantum mechanics and a "condensed" phase dominated by quantum mechanics. The phase diagram exhibits two critical points, one in each ensemble, beyond which the phase transitions disappear. There also exist a region of negative specific heats and a situation of ensemble inequivalence for sufficiently large systems. We apply the fermionic King model to the case of dark matter halos made of massive neutrinos. The gaseous phase describes large halos and the condensed phase describes dwarf halos. Partially degenerate configurations describe intermediate size halos. We argue that large dark matter halos cannot harbor a fermion ball because these nucleus-halo configurations are thermodynamically unstable (saddle points of entropy). Large dark matter halos may rather contain a central black hole resulting from a dynamical instability of relativistic origin occurring during the gravothermal catastrophe.
This paper is devoted to the cubic nonlinear Schrödinger equation in a two dimensional waveguide with shrinking cross section of order $ε$. For a Cauchy data living essentially on the first mode of the transverse Laplacian, we provide a tensorial approximation of the solution $ψ^ε$ in the limit $ε\to 0$, with an estimate of the approximation error, and derive a limiting nonlinear Schrödinger equation in dimension one. If the Cauchy data $ψ^ε_0$ has a uniformly bounded energy, then it is a bounded sequence in $H^1$ and we show that the approximation is of order $\mathcal O(\sqrtε)$. If we assume that $ψ^ε_0$ is bounded in the graph norm of the Hamiltonian, then it is a bounded sequence in $H^2$ and we show that the approximation error is of order $\mathcal O(ε)$.
We study the fermionic King model which may provide a relevant model of dark matter halos. The exclusion constraint can be due to quantum mechanics (for fermions such as massive neutrinos) or to Lynden-Bell's statistics (for collisionless systems undergoing violent relaxation). This model has a finite mass. Furthermore, a statistical equilibrium state exists for all accessible values of energy. Dwarf and intermediate size halos are degenerate quantum objects stabilized against gravitational collapse by the Pauli exclusion principle. Large halos at sufficiently high energies are in a gaseous phase where quantum effects are negligible. They are stabilized by thermal motion. Below a critical energy they undergo gravitational collapse (gravothermal catastrophe). This may lead to the formation of a central black hole that does not affect the structure of the halo. This may also lead to the formation of a compact degenerate object surrounded by a hot massive atmosphere extending at large distances. We argue that large dark matter halos should not contain a degenerate nucleus (fermion ball) because these nucleus-halo structures are thermodynamically unstable. We compare the rotation curves of the classical King model to observations of large dark matter halos (Burkert profile). Because of collisions and evaporation, the central density increases while the slope of the halo density profile decreases until an instability takes place. We find that the observations are compatible with a King profile at, or close to, the point of marginal stability in the microcanonical ensemble. At that point, the King profile can be fitted by the modified Hubble profile that has a flat core and a halo in which the density decreases as $r^{-3}$. This is qualitatively similar to the Burkert profile.
We present and analyze two mathematical models for the self consistent quantum transport of electrons in a graphene layer. We treat two situations. First, when the particles can move in all the plane $\RR^2$, the model takes the form of a system of massless Dirac equations coupled together by a selfconsistent potential, which is the trace in the plane of the graphene of the 3D Poisson potential associated to surface densities. In this case, we prove local in time existence and uniqueness of a solution in $H^s(\RR^2)$, for $s > 3/8$ which includes in particular the energy space $H^{1/2}(\RR^2)$. The main tools that enable to reach $s\in (3/8,1/2)$ are the dispersive Strichartz estimates that we generalized here for mixed quantum states. Second, we consider a situation where the particles are constrained in a regular bounded domain $Ω$. In order to take into account Dirichlet boundary conditions which are not compatible with the Dirac Hamiltonian $H_{0}$, we propose a different model built on a modified Hamiltonian displaying the same energy band diagram as $H_{0}$ near the Dirac points. The well-posedness of the system in this case is proved in $H^s_{A}$, the domain of the fractional order Dirichlet Laplacian operator, for $1/2\leq s<5/2$.
We investigate the controllability of quantum electrons trapped in a two-dimensional device, typically a MOS field-effect transistor. The problem is modeled by the Schrödinger equation in a bounded domain coupled to the Poisson equation for the electrical potential. The controller acts on the system through the boundary condition on the potential, on a part of the boundary modeling the gate. We prove that, generically with respect to the shape of the domain and boundary conditions on the gate, the device is controllable. We also consider control properties of a more realistic nonlinear version of the device, taking into account the self-consistent electrostatic Poisson potential.
In this paper, we present the Stroboscopic Averaging Method (SAM), recently introduced in [7,8,10,12], which aims at numerically solving highly-oscillatory differential equations. More specifically, we first apply SAM to the Schrödinger equation on the 1-dimensional torus and on the real line with harmonic potential, with the aim of assessing its efficiency: as compared to the well-established standard splitting schemes, the stiffer the problem is, the larger the speed-up grows (up to a factor 100 in our tests). The geometric properties of SAM are also explored: on very long time intervals, symmetric implementations of the method show a very good preservation of the mass invariant and of the energy. In a second series of experiments on 2-dimensional equations, we demonstrate the ability of SAM to capture qualitatively the long-time evolution of the solution (without spurring high oscillations).
This work is devoted to the numerical simulation of nonlinear Schrödinger and Klein-Gordon equations. We present a general strategy to construct numerical schemes which are uniformly accurate with respect to the oscillation frequency. This is a stronger feature than the usual so called "Asymptotic preserving" property, the last being also satisfied by our scheme in the highly oscillatory limit. Our strategy enables to simulate the oscillatory problem without using any mesh or time step refinement, and the orders of our schemes are preserved uniformly in all regimes. In other words, since our numerical method is not based on the derivation and the simulation of asymptotic models, it works in the regime where the solution does not oscillate rapidly, in the highly oscillatory limit regime, and in the intermediate regime with the same order of accuracy. In the same spirit as in [5], the method is based on two main ingredients. First, we embed our problem in a suitable "two-scal" reformulation with the introduction of an additional variable. Then a link is made with classical strategies based on Chapman-Enskog expansions in kinetic theory despite the dispersive context of the targeted equations, allowing to separate the fast time scale from the slow one. Uniformly accurate (UA) schemes are eventually derived from this new formulation and their properties and performances are assessed both theoretically and numerically.
We consider the semiclassical limit for the nonlinear Schrodinger equation. We introduce a phase/amplitude representation given by a system similar to the hydrodynamical formulation, whose novelty consists in including some asymptotically vanishing viscosity. We prove that the system is always locally well-posed in a class of Sobolev spaces, and globally well-posed for a fixed positive Planck constant in the one-dimensional case. We propose a second order numerical scheme which is asymptotic preserving. Before singularities appear in the limiting Euler equation, we recover the quadratic physical observables as well as the wave function with mesh size and time step independent of the Planck constant. This approach is also well suited to the linear Schrodinger equation.
This work is devoted to the numerical simulation of a Vlasov-Poisson model describing a charged particle beam under the action of a rapidly oscillating external electric field. We construct an Asymptotic Preserving numerical scheme for this kinetic equation in the highly oscillatory limit. This scheme enables to simulate the problem without using any time step refinement technique. Moreover, since our numerical method is not based on the derivation of the simulation of asymptotic models, it works in the regime where the solution does not oscillate rapidly, and in the highly oscillatory regime as well. Our method is based on a "double-scale" reformulation of the initial equation, with the introduction of an additional periodic variable.
We introduce a new micro-macro decomposition of collisional kinetic equations in the specific case of the diffusion limit, which naturally incorporates the incoming boundary conditions. The idea is to write the distribution function $f$ in all its domain as the sum of an equilibrium adapted to the boundary (which is not the usual equilibrium associated with $f$) and a remaining kinetic part. This equilibrium is defined such that its incoming velocity moments coincide with the incoming velocity moments of the distribution function. A consequence of this strategy is that no artificial boundary condition is needed in the micro-macro models and the exact boundary condition on $f$ is naturally transposed to the macro part of the model. This method provides an 'Asymptotic preserving' numerical scheme which generates a very good approximation of the space boundary values at the diffusive limit, without any mesh refinement in the boundary layers. Our numerical results are in very good agreement with the exact so-called Chandrasekhar value, which is explicitely known in some simple cases.
In this work, we study the orbital stability of steady states and the existence of blow-up self-similar solutions to the so-called Vlasov-Manev (VM) system. This system is a kinetic model which has a similar Vlasov structure as the classical Vlasov-Poisson system, but is coupled to a potential in $-1/r- 1/r^2$ (Manev potential) instead of the usual gravitational potential in $-1/r$, and in particular the potential field does not satisfy a Poisson equation but a fractional-Laplacian equation. We first prove the orbital stability of the ground states type solutions which are constructed as minimizers of the Hamiltonian, following the classical strategy: compactness of the minimizing sequences and the rigidity of the flow. However, in driving this analysis, there are two mathematical obstacles: the first one is related to the possible blow-up of solutions to the VM system, which we overcome by imposing a sub-critical condition on the constraints of the variational problem. The second difficulty (and the most important) is related to the nature of the Euler-Lagrange equations (fractional-Laplacian equations) to which classical results for the Poisson equation do not extend. We overcome this difficulty by proving the uniqueness of the minimizer under equimeasurabilty constraints, using only the regularity of the potential and not the fractional-Laplacian Euler-Lagrange equations itself. In the second part of this work, we prove the existence of exact self-similar blow-up solutions to the Vlasov-Manev equation, with initial data arbitrarily close to ground states. This construction is based on a suitable variational problem with equimeasurability constraint.
We introduce a new micro/macro decomposition of collisional kinetic equations which naturally incorporates the exact space boundary conditions. The idea is to write the distribution fonction $f$ in all its domain as the sum of a Maxwellian adapted to the boundary (which is not the usual Maxwellian associated with $f$) and a reminder kinetic part. This Maxwellian is defined such that its 'incoming' velocity moments coincide with the 'incoming' velocity moments of the distribution function. Important consequences of this strategy are the following. i) No artificial boundary condition is needed in the micro/macro models and the exact boundary condition on $f$ is naturally transposed to the macro part of the model. ii) It provides a new class of the so-called 'Asymptotic preserving' (AP) numerical schemes: such schemes are consistent with the original kinetic equation for all fixed positive value of the Knudsen number $\eps$, and if $\eps \to 0 $ with fixed numerical parameters then these schemes degenerate into consistent numerical schemes for the various corresponding asymptotic fluid or diffusive models. Here, the strategy provides AP schemes not only inside the physical domain but also in the space boundary layers. We provide a numerical test in the case of a diffusion limit of the one-group transport equation, and show that our AP scheme recovers the boundary layer and a good approximation of the theoretical boundary value, which is usually computed from to the so-called Chandrasekhar function.
We address the following inverse problem in quantum statistical physics: does the quantum free energy (von Neumann entropy + kinetic energy) admit a unique minimizer among the density operators having a given local density $n(x)$? We give a positive answer to that question, in dimension one. This enables to define rigourously the notion of local quantum equilibrium, or quantum Maxwellian, which is at the basis of recently derived quantum hydrodynamic models and quantum drift-diffusion models. We also characterize this unique minimizer, which takes the form of a global thermodynamic equilibrium (canonical ensemble) with a quantum chemical potential.
We study the Schrödinger-Poisson system on the unit sphere $\SS^2$ of $\RR^3$, modeling the quantum transport of charged particles confined on a sphere by an external potential. Our first results concern the Cauchy problem for this system. We prove that this problem is regularly well-posed on every $H^s(\SS ^2)$ with $s>0$, and not uniformly well-posed on $L^2(\SS ^2)$. The proof of well-posedness relies on multilinear Strichartz estimates, the proof of ill-posedness relies on the construction of a counterexample which concentrates exponentially on a closed geodesic. In a second part of the paper, we prove that this model can be obtained as the limit of the three dimensional Schrödinger-Poisson system, singularly perturbed by an external potential that confines the particles in the vicinity of the sphere.