  In the last decades, the experimental research on Bose-Einstein
interferometry has received much attention due to promising technological
implications. This has thus motivated the development of numerical simulations
aimed at solving the time-dependent Gross-Pitaevskii equation and its reduced
one-dimensional version to better understand the development of
interference-type features and the subsequent soliton dynamics. In this work,
Bohmian mechanics is considered as an additional tool to further explore and
analyze the formation and evolution in real time of the soliton arrays that
follow the merging of two condensates. An alternative explanation is thus
provided in terms of an underlying dynamical velocity field, directly linked to
the local phase variations undergone by the condensate along its evolution.
Although the reduced one-dimensional model is considered here, it still
captures the essence of the phenomenon, rendering a neat picture of the full
evolution without diminishing the generality of the description. To better
appreciate the subtleties of free versus bound dynamics, two cases are
discussed. First, the soliton dynamics exhibited by a coherent superposition of
two freely released condensates is studied, discussing the peculiarities of the
underlying velocity field and the corresponding flux trajectories in terms of
both the peak-to-peak distance between the two initial clouds and the addition
of a phase difference between them. In the latter case, an interesting
correspondence with the well-known Aharonov-Bohm effect is found. Then, the
recurrence dynamics displayed by the more general case of two condensates
released from the two opposite turning points of a harmonic trap is considered
in terms of the distance between such turning points. [...]
