Search

In this article we are interested for the numerical computation of spectra of non-self adjoint quadratic operators, in two and three spatial dimensions. Indeed, in the multidimensional case very few results are known on the location of the eignevalues. This leads to solve nonlinear eigenvalue problems. In introduction we begin with a review of theoretical results and numerical results obtained for the one dimensional case. Then we present the numerical methods developed to compute the spectra (finite difference discretization) for the two and three dimensional cases. The numerical results obtained are presented and analyzed. One difficulty here is that we have to compute eigenvalues of strongly non-self-adjoint operators which are unstable. This work is in continuity of a previous work in one spatial dimension.

We consider the time dependent Schr\"odinger equation with a coupling spin-orbit in the semi-classical regime $\hbar\searrow 0$ and large spin number $\spin\rightarrow +\infty$ such that $\hbar^\delta\spin=c$ where $c>0$ and $\delta>0$ are constant. The initial state $\Psi(0)$ is a product of an orbital coherent state in $L^2(\R^d)$ and a spin coherent state in a spin irreducible representation space ${\mathcal H}_{2\spin +1}$. For $\delta <1$, at the leading order in $\hbar$, the time evolution $\Psi(t)$ of $ \Psi(0)$ is well approximated by the product of an orbital and a spin coherent state. Nevertheless for $1/2<\delta<1$ the quantum orbital leaves the classical orbital. For $\delta=1$ we prove that this last claim is no more true when the interaction depends on the orbital variables. For the Dicke model, we prove that the orbital partial trace of the projector on $\Psi(t)$ is a mixed state in $L^2(\R)$ for small $t>0$.

We consider a modulated magnetic field, $B(t) = B_0 +\varepsilon f(\omega t)$, perpendicular to a fixed plane, where $B_0$ is constant, $\varepsilon>0$ and $f$ a periodic function on the torus ${\mathbb T}^n$. Our aim is to study classical and quantum dynamics for the corresponding Landau Hamiltonian. It turns out that the results depend strongly on the chosen gauge. For the Landau gauge the position observable is unbounded for "almost all" non resonant frequencies $\omega$. On the contrary, for the symmetric gauge we obtain that, for "almost all" non resonant frequencies $\omega$, the Landau Hamiltonian is reducible to a two dimensional harmonic oscillator and thus gives rise to bounded dynamics. The proofs use KAM algorithms for the classical dynamics. Quantum applications are given. In particular, the Floquet spectrum is absolutely continuous in the Landau gauge while it is discrete, of finite multiplicity, in symmetric gauge.

This paper is devoted to the construction of approximations of the propagator associated with a semi-classical matrix-valued Schr\"odinger operator with symbol presenting smooth eigenvalues crossings. Inspired by the approach of the theoretical chemists Herman and Kluk who propagated continuous superpositions of Gaussian wave-packets for scalar equations, we consider frozen and thawed Gaussian initial value representations that incorporate classical transport and branching processes along a hopping hypersurface. Based on the Gaussian wave-packet framework, our result relies on an accurate analysis of the solutions of the associated Schr\"odinger equation for data that are vector-valued wave-packets. We prove that these solutions are asymptotic to wavepackets at any order in terms of the semi-classical parameter.

In the phase space $\R^{2d}$, let us denote $\{A,B\}$ the Poisson bracket of two smooth classical observables and $\{A, B\}_\circledast $ their Moyal bracket, defined as the Weyl symbol of $i[ A, B]$, where $ \hat A$ is the Weyl quantization of $A$ and $[ \hat A, \hat B]= \hat A \hat B- \hat B \hat A$ (commutator). In this note we prove that if a smooth Hamiltonian $H$ on the phase space $\R^{2d}$, with derivatives of moderate growth, satisfies $\{A,H\}= \{A, H\}_\circledast$ for any smooth and bounded observable $A$ then $H$ must be a polynomial of degree at most 2. This is related with the Groenewold-van Hove Theorem \cite{Gotay, Groen, vHove} concerning quantization of polynomial observables.

We review some recent results obtained for the time evolution of wave packets for systems of equations of pseudo-differential type, including Schr{\"o}dinger ones, and discuss their application to the approximation of the associated unitary propagator. We start with scalar equations, propagation of coherent states, and applications to the Herman-Kluk approximation. Then we discuss the extension of these results to systems with eigenvalues of constant multiplicity or with smooth crossings.

We analyze the propagation of coherent states through general systems of pseudodifferential form associated with Hamiltonian presenting codimension one eigen-value crossings. In particular, we calculate precisely the non adiabatic effects of the crossing in terms of a transition operator.

We prove an abstract theorem giving a $\langle t\rangle^\epsilon$ bound ($\forall \epsilon>0$) on the growth of the Sobolev norms in linear Schr\"odinger equations of the form $i \dot \psi = H_0 \psi + V(t) \psi $ when the time $t \to \infty$. The abstract theorem is applied to several cases, including the cases where (i) $H_0$ is the Laplace operator on a Zoll manifold and $V(t)$ a pseudodifferential operator of order smaller then 2; (ii) $H_0$ is the (resonant or nonresonant) Harmonic oscillator in $R^d$ and $V(t)$ a pseudodifferential operator of order smaller then $H_0$ depending in a quasiperiodic way on time. The proof is obtained by first conjugating the system to some normal form in which the perturbation is a smoothing operator and then applying the results of \cite{MaRo}.

We prove a reducibility result for a quantum harmonic oscillator in arbitrary dimensions with arbitrary frequencies perturbed by a linear operator which is a polynomial of degree two in $x_j$, $-i \partial_j$ with coefficients which depend quasiperiodically on time., Comment: fixed some misprints

In this paper we consider time dependent Schr{\"o}dinger linear PDEs of the form i$\partial$t$\psi$ = L(t)$\psi$, where L(t) is a continuous family of self-adjoint operators. We give conditions for well-posedness and polynomial growth for the evolution in abstract Sobolev spaces. If L(t) = H + V (t) where V (t) is a perturbation smooth in time and H is a self-adjoint positive operator whose spectrum can be enclosed in spectral clusters whose distance is increasing, we prove that the Sobolev norms of the solution grow at most as t $\epsilon$ when t $\rightarrow$ $\infty$, for any $\epsilon$ \textgreater{} 0. If V (t) is analytic in time we improve the bound to (log t) $\gamma$ , for some $\gamma$ \textgreater{} 0. The proof follows the strategy, due to Howland, Joye and Nenciu, of the adiabatic approximation of the flow. We recover most of known results and obtain new estimates for several models including 1-degree of freedom Schr{\"o}dinger operators on R and Schr{\"o}dinger operators on Zoll manifolds.

In this article we are interested for the numerical study of nonlinear eigenvalue problems. We begin with a review of theoretical results obtained by functional analysis methods, especially for the Schrodinger pencils. Some recall are given for the pseudospectra. Then we present the numerical methods and results obtained for eigenvalues computation with spectral methods and finite difference discretization, in infinite or bounded domains. Comparison with theoretical results is done. The main difficulty here is that we have to compute eigenvalues of strongly non-self-adjoint operators which are very unstable.

Our main goal in this paper is to prove existence (and uniqueness) of the quantum propagator for time dependent quantum Hamiltonians $\hat H(t)$ when this Hamiltonian is perturbed with a quadratic white noise $\dot{\beta}\hat K$. $\beta$ is a continuous function in time $t$, $\dot \beta$ its time derivative and $K$ is a quadratic Hamiltonian. $\hat K$ is the Weyl quantization of $K$. For time dependent quadratic Hamiltonians $H(t)$ we recover, under less restrictive assumptions, the results obtained in \cite{bofu, du}.In our approach we use an exact Hermann Kluk formula \cite{ro2} to deduce a Strichartz estimate for the propagator of $\hat H(t) +\dot \beta K$. This is applied to obtain local and global well posedness for solutions for non linear Schr\"odinger equations with an irregular time dependent linear part.