Dark soliton collisions in a toroidal Bose-Einstein condensate

We study the dynamics of two gray solitons in a Bose-Einstein condensate confined by a toroidal trap with a tight confinement in the radial direction. Gross-Pitaevskii simulations show that solitons can be long living objects passing through many collisional processes. We have observed quite different behaviors depending on the soliton velocity. Very slow solitons, obtained by perturbing the stationary solitonic profile, move with a constant angular velocity until they collide elastically and move in the opposite direction without showing any sign of lowering their energy. In this case the density notches are always well separated and the fronts are sharp and straight. Faster solitons present vortices around the notches, which play a central role during the collisions. We have found that in these processes the solitons lose energy, as the outgoing velocity turns out to be larger than the incoming one. To study the dynamics, we model the gray soliton state with a free parameter that is related to the soliton velocity. We further analyze the energy, soliton velocity and turning points in terms of such a free parameter, finding that the main features are in accordance with the infinite one-dimensional system.Comment: 15 pages, 11 figures. Accepted in PR

Two-mode effective interaction in a double-well condensate

We investigate the origin of a disagreement between the two-mode model and the exact Gross-Pitaevskii dynamics applied to double-well systems. In general this model, even in its improved version, predicts a faster dynamics and underestimates the critical population imbalance separating Josephson and self-trapping regimes. We show that the source of this mismatch in the dynamics lies in the value of the on-site interaction energy parameter. Using simplified Thomas-Fermi densities, we find that the on-site energy parameter exhibits a linear dependence on the population imbalance, which is also confirmed by Gross-Pitaevskii simulations. When introducing this dependence in the two-mode equations of motion, we obtain a reduced interaction energy parameter which depends on the dimensionality of the system. The use of this new parameter significantly heals the disagreement in the dynamics and also produces better estimates of the critical imbalance.Comment: 5 pages, 4 figures, accepted in PR

Benchmarking the Variational Reduced Density Matrix Theory in the Doubly Occupied Configuration Interaction Space with Integrable Pairing Models

The variational reduced density matrix theory has been recently applied with great success to models within the truncated doubly occupied configuration interaction space, which corresponds to the seniority zero subspace. Conservation of the seniority quantum number restricts the Hamiltonians to be based on the SU(2) algebra. Among them there is a whole family of exactly solvable Richardson-Gaudin pairing Hamiltonians. We benchmark the variational theory against two different exactly solvable models, the Richardson-Gaudin-Kitaev and the reduced BCS Hamiltonians. We obtain exact numerical results for the so-called PQGT N-representability conditions in both cases for systems that go from 10 to 100 particles. However, when random single-particle energies as appropriate for small superconducting grains are considered, the exactness is lost but still a high accuracy is obtained.Fil: Rubio García, A.. Instituto de Estructura de la Materia; España. Consejo Superior de Investigaciones Científicas; EspañaFil: Alcoba, Diego Ricardo. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Física de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Física de Buenos Aires; ArgentinaFil: Capuzzi, Pablo. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Física de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Física de Buenos Aires; ArgentinaFil: Dukelsky, J.. Consejo Superior de Investigaciones Científicas; España. Instituto de Estructura de la Materia; Españ

Vortex nucleation processes in rotating lattices of Bose-Einstein condensates ruled by the on-site phases

We study the nucleation and dynamics of vortices in rotating lattice potentials where weakly linked condensates are formed with each condensate exhibiting an almost axial symmetry. Due to such a symmetry, the on-site phases acquire a linear dependence on the coordinates as a result of the rotation, which allows us to predict the position of vortices along the low density paths that separate the sites. We first show that, for a system of atoms loaded in a four-site square lattice potential, subject to a constant rotation frequency, the analytical expression that we obtain for the positions of vortices of the stationary arrays accurately reproduces the full three-dimensional Gross-Pitaevskii results. We then study the time-dependent vortex nucleation process when a linear ramp of the rotation frequency is applied to a lattice with sixteen sites. We develop a formula for the number of nucleated vortices which turns to have a linear dependence on the rotation frequency with a smaller slope than that of the standard estimate which is valid in absence of the lattice. From time-dependent Gross-Pitaevskii simulations we further find that the on-site populations remain almost constant during the time evolution instead of spreading outwards, as expected from the action of the centrifugal force. Therefore, the time-dependent phase difference between neighboring sites acquires a running behavior typical of a self-trapping regime. We finally show that, in accordance with our predictions, this fast phase-difference evolution provokes a rapid vortex motion inside the lattice. Our analytical expressions may be useful for describing other vortex processes in systems with the same on-site axial symmetry

Structure of vortices in two-component Bose-Einstein condensates

We develop a three-dimensional analysis of the phase separation of two-species Bose-Einstein condensates in the presence of vorticity within the Thomas-Fermi approximation. We find different segregation features according to whether the more repulsive component is in a vortex or in a vortex-free state. An application of this study is aimed at describing systems formed by two almost immiscible species of rubidium-87 that are commonly used in Bose-Einstein condensation experiments. In particular, in this work we calculate the density profiles of condensates for the same conditions as the states prepared in the experiments performed at JILA [Matthews et al., Phys. Rev. Lett. 83, 2498 (1999)]Comment: 4 pages, 3 figure

Barrier effects on the collective excitations of split Bose-Einstein condensates

We investigate the collective excitations of a single-species Bose gas at T=0 in a harmonic trap where the confinement undergoes some splitting along one spatial direction. We mostly consider onedimensional potentials consisting of two harmonic wells separated a distance 2 z_0, since they essentially contain all the barrier effects that one may visualize in the 3D situation. We find, within a hydrodynamic approximation, that regardless the dimensionality of the system, pairs of levels in the excitation spectrum, corresponding to neighbouring even and odd excitations, merge together as one increases the barrier height up to the current value of the chemical potential. The excitation spectra computed in the hydrodynamical or Thomas-Fermi limit are compared with the results of exactly solving the time-dependent Gross-Pitaevskii equation. We analyze as well the characteristics of the spatial pattern of excitations of threedimensional boson systems according to the amount of splitting of the condensate.Comment: RevTeX, 12 pages, 13 ps figure

Random-phase approximation study of collective excitations in the Bose-Fermi mixed condensate of alkali-metal gases

We perform Random Phase Approximation (RPA) study of collective excitations in the bose-fermi mixed degenerate gas of Alkali-metal atoms at T=0. The calculation is done by diagonalization in a model space composed of particle-hole type excitations from the ground state, the latter being obtained from the coupled Gross-Pitaevskii and Thomas-Fermi equations. We investigate strength distributions for different combinations of bose and fermi multipole ($L$) operators with $L=0,1,2,3$. Transition densities and dynamical structure factors are calculated for collective excitations. Comparison with the sum rule prediction for the collective frequency is given. Time dependent behavior of the system after an external impulse is studied.Comment: 28 pages, 13 figures, submitted to Phys. Rev.

Overlap functions in correlation methods and quasifree nucleon knockout from 16^{16}O

The cross sections of the ($e,e'N$) and ($\gamma,p$) reactions on $^{16}$O are calculated, for the transitions to the $1/2^{-}$ ground state and the first $3/2^{-}$ excited state of the residual nucleus, using single-particle overlap functions obtained on the basis of one-body density matrices within different correlation methods. The electron-induced one-nucleon knockout reaction is treated within a nonrelativistic DWIA framework. The theoretical treatment of the ($\gamma,p$) reaction includes both contributions of the direct knockout mechanism and of meson-exchange currents. The results are sensitive to details of the different overlap functions. The consistent analysis of the reaction cross sections and the comparison with the experimental data make it possible to study the nucleon--nucleon correlation effects.Comment: 26 pages, LaTeX, 5 Postscript figures, submitted to PR