Search

The Homotopy Analysis Method (HAM) is a powerful technique which allows to derive approximate solutions of both ordinary and partial differential equations. We propose to use a variational approach based on the Least Action Principle (LAP) in order to improve the efficiency of the HAM when applied to Lagrangian systems. The extremization of the action is achieved by varying the HAM parameter, therefore controlling the accuracy of the approximation. As case studies we consider the harmonic oscillator, the cubic and the quartic anharmonic oscillators, and the Korteweg-de Vries partial differential equation. We compare our results with those obtained using the standard approach, which is based on the residual error square method. We see that our method accelerates the convergence of the HAM parameter to the exact value in the cases in which the exact solution is known. When the exact solution is not analytically known, we find that our method performs better than the standard HAM for the cases we have analyzed. Moreover, our method shows better performance when the order of the approximation is increased and when the nonlinearity of the equations is stronger., Comment: 26 pages, 7 figures, 15 tables and 3 appendices

We study the effective dynamics of ferromagnetic spin chains in presence of long-range interactions. We consider the Heisenberg Hamiltonian in one dimension for which the spins are coupled through power-law long-range exchange interactions with exponent $\alpha$. We add to the Hamiltonian an anisotropy in the $z$-direction. In the framework of a semiclassical approach, we use the Holstein-Primakoff transformation to derive an effective long-range discrete nonlinear Schr\"odinger equation. We then perform the continuum limit and we obtain a fractional nonlinear Schr\"odinger-like equation. Finally, we study the modulational instability of plane-waves in the continuum limit and we prove that, at variance with the short-range case, plane waves are modulationally unstable for $\alpha < 3$. We also study the dependence of the modulation instability growth rate and critical wave-number on the parameters of the Hamiltonian and on the exponent $\alpha$.

We study analytical solutions of the Fractional Boussinesq Equation (FBE), which is an effective model for the Fermi-Pasta-Ulam (FPU) one-dimensional lattice with long-range couplings. The couplings decay as a power-law with exponent s, with 1 < s < 3, so that the energy density is finite, but s is small enough to observe genuine long-range effects. The analytic solutions are obtained by introducing an ansatz for the dependence of the field on space and time. This allows to reduce the FBE to an ordinary differential equation, which can be explicitly solved. The solutions are initially localized and they delocalize progressively as time evolves. Depending on the value of s the solution is either a pulse (meaning a bump) or an anti-pulse (i.e., a hole) on a constant field for 1 < s < 2 and 2 < s < 3, respectively., Comment: 10 pages, 2 figures. The paper is accepted in JSTAT, the special issue "New Trends in Nonequilibrium Statistical Mechanics: Classical and Quantum Systems (nesmcq18)"

We investigate the existence and propagation of solitons in a long-range extension of the quartic Fermi-Pasta-Ulam (FPU) chain of anharmonic oscillators. The coupling in the linear term decays as a power-law with an exponent greater than 1 and less than 3. We obtain an analytic perturbative expression of traveling envelope solitons by introducing a Non Linear Schrodinger (NLS) equation for the slowly varying amplitude of short wavelength modes. Due to the non analytic properties of the dispersion relation, it is crucial to develop the theory using discrete difference operators. Those properties are also the ultimate reason why kink-solitons may exist but are unstable, at variance with the short-range FPU model. We successfully compare these approximate analytic results with numerical simulations., Comment: 10 pages, 4 figures

We study instabilities and relaxation to equilibrium in a long-range extension of the Fermi-Pasta-Ulam-Tsingou (FPU) oscillator chain by exciting initially the lowest Fourier mode. Localization in mode space is stronger for the long-range FPU model. This allows us to uncover the sporadic nature of instabilities, i.e., by varying initially the excitation amplitude of the lowest mode, which is the control parameter, instabilities occur in narrow amplitude intervals. Only for sufficiently large values of the amplitude, the system enters a permanently unstable regime. These findings also clarify the long-standing problem of the relaxation to equilibrium in the short-range FPU model. Because of the weaker localization in mode space of this latter model, the transfer of energy is retarded and relaxation occurs on a much longer time-scale.

We study the eigenstates of two opposite spin fermions on a one-dimensional lattice with finite range interaction. The eigenstates are projected onto the set of Fock eigenstates of the noninteracting case. We find antiresonances for symmetric eigenstates, which eliminate the interaction between two symmetric Fock states when satisfying a corresponding selection rule., Comment: 4 pages, 3 figures

We investigate the spectrum and eigenstates of a Bose-Hubbard chain containing two bosons with fixed boundary conditions. In the noninteracting case the eigenstates of the system define a two-dimensional normal-mode space. For the interacting case weight functions of the eigenstates are computed by perturbation theory and numerical diagonalization. We identify paths in the two-dimensional normal-mode space which are rims for the weight functions. The decay along and off the rims is algebraic. Intersection of two paths (rims) leads to a local enhancement of the weight functions. We analyze nonperturbative effects due to the degeneracies and the formation of two-boson bound states., Comment: 10 pages, 11 figures. With minor corrections. Accepted in Physica D

We study the spectrum and eigenstates of the quantum discrete Bose-Hubbard Hamiltonian in a finite one-dimensional lattice containing two bosons. The interaction between the bosons leads to an algebraic localization of the modified extended states in the normal mode space of the noninteracting system. Weight functions of the eigenstates in the space of normal modes are computed by using numerical diagonalization and perturbation theory. We find that staggered states do not compactify in the dilute limit for large chains., Comment: 7 pages, 7 figures. Minor changes and additional comments. Acepted in Physical Review B