An alternative route towards completely recovering nuclear quantum effects, while still maintaining a trajectory-based representation of nuclear dynamics, is provided by the socalled quantum trajectory methods (QTMs). In past years, such QTMs have been used to describe the quantum nuclear behaviour of different model and molecular systems. In this work, our aim is to focus on the theoretical and numerical aspects of one particular quantum trajectory formalism, and its application to adiabatic and nonadiabatic problems.