We introduce two exponential-type integrators for the “good” Boussinesq equation. They are of orders one and two, respectively, and they require lower spatial regularity of the solution compared to classical exponential integrators. For the first integrator, we prove first-order convergence in H^r for solutions in H^{r+1} with r>1/2. This new integrator even converges (with lower order) in H^r for solutions in H^r, i.e., without any additional smoothness assumptions. For the second integrator, we prove second-order convergence in H^r for solutions in H^{r+3} with r>1/2 and convergence in L^2 for solutions in H^3. Numerical results are reported to illustrate the established error estimates. The experiments clearly demonstrate that the new exponential-type integrators are favorable over classical exponential integrators for initial data with low regularity.

Let $f(x)=\sum_{n=0}^{\infty}\frac{1}{n!}q^{n(n-1)/2}x^n$ ($0<q<1$) be the deformed exponential function. It is known that the zeros of $f(x)$ are real and form a negative decreasing sequence $(x_k)$ ($k\ge 1$). We investigate the complete asymptotic expansion for $x_{k}$ and prove that for any $n\ge1$, as $k\to \infty$, \begin{align*} x_k=-kq^{1-k}\Big(1+\sum_{i=1}^{n}C_i(q)k^{-1-i}+o(k^{-1-n})\Big), \end{align*} where $C_i(q)$ are some $q$ series which can be determined recursively. We show that each $C_{i}(q)\in \mathbb{Q}[A_0,A_1,A_2]$, where $A_{i}=\sum_{m=1}^{\infty}m^i\sigma(m)q^m$ and $\sigma(m)$ denotes the sum of positive divisors of $m$. When writing $C_{i}$ as a polynomial in $A_0, A_1$ and $A_2$, we find explicit formulas for the coefficients of the linear terms by using Bernoulli numbers. Moreover, we also prove that $C_{i}(q)\in \mathbb{Q}[E_2,E_4,E_6]$, where $E_2$, $E_4$ and $E_6$ are the classical Eisenstein series of weight 2, 4 and 6, respectively.

No abstract uploaded!

In the present paper, we prove the existence of solutions $(\lambda_1,\lambda_2,u,v)\in\R^2\times H^1(\R^3,\R^2)$ to systems of coupled Schr\"odinger equations $$ \begin{cases} -\Delta u+\lambda_1u=\mu_1 u^3+\beta uv^2\quad &\hbox{in}\;\R^3\\ -\Delta v+\lambda_2v=\mu_2 v^3+\beta u^2v\quad&\hbox{in}\;\R^3\\ u,v>0&\hbox{in}\;\R^3 \end{cases} $$ satisfying the normalization constraint $ \displaystyle\int_{\R^3}u^2=a^2\quad\hbox{and}\;\int_{\R^3}v^2=b^2, $ which appear in binary mixtures of Bose-Einstein condensates or in nonlinear optics. The parameters $\mu_1,\mu_2,\beta>0$ are prescribed as are the masses $a,b>0$. The system has been considered mostly in the case of fixed frequencies $\lambda_1,\lambda_2$. When the masses are prescribed, the standard approach to this problem is variational with $\lambda_1,\lambda_2$ appearing as Lagrange multipliers. Here we present a new approach based on the fixed point index in cones, bifurcation theory, and the continuation method. We obtain the existence of normalized solutions for any given $a,b>0$ for $\beta$ in a large range. We also have a result about the nonexistence of positive solutions which shows that our existence theorem is almost optimal. Especially, if $\mu_1=\mu_2$ we prove that normalized solutions exist for all $\beta>0$ and all $a,b>0$.

No abstract uploaded!