It is common knowledge that a dark soliton can be excited in an ultra-cold atomic gas by means of the phase imprinting method. We show that, for a superfluid fermionic mixture, the standard phase imprinting procedure applied to both components fails to create a state with symmetry properties identical to those of the dark soliton solution of the Bogoliubov-de Gennes equations. To produce a dark soliton in the BCS regime, a single component of the Fermi mixture should be phase imprinted only.PACS numbers: 67.85. Lm, 03.75.Ss, 03.75.Lm Solitons, or solitary waves, are solutions of non-linear wave equations that can propagate without change of shapes. Electromagnetic solitons have been intensively studied in non-linear optics [1]. Ultra-cold atomic gases offer a playground for investigation of matter-wave solitons. At low temperature Bose atomic gases form BoseEinstein condensates (BEC) which, in the mean field approximation, can be described by a single-particle nonlinear Gross-Pitaevskii equation (GPE) [2]. Depending on the sign of the s-wave scattering length of atoms, the GPE can possess bright or dark soliton solutions [3,4]. Both kinds of solitons have been created in a laboratory [5][6][7][8]. Signatures of the quantum nature of solitons -beyond the mean-field GPE -have been predicted [9-14], but have not been observed so far in ultra-cold atomic gases.At low temperature a two-species Fermi gas undergoes a transition to a superfluid phase if the particle interactions are attractive. In the weak coupling BardeenCooper-Schrieffer (BCS) regime, the system is described by a set of non-linear Bogoliubov-de Gennes equations [15]. These equations describe the ground state of the atomic gas but they can also describe a dark soliton solution where particle densities are nearly the same as for the ground state case but the BCS pairing function possesses a phase flip at the position of the soliton [16][17][18]. Similar solutions appear also in the theory of conducting polymers where, however, the order parameter is real [19]. On the BEC side of the BCS-BEC crossover regime, the BCS pairing function can be identified with the condensate wave-function of a molecular condensate corresponding to the dark soliton solution of the GPE [18,20].Dark solitons in Bose gases are excited experimentally by means of a phase imprinting method where half of the cloud acquires a phase π after a short interaction with a laser radiation [5,6]. A similar procedure was applied in a superfluid Fermi mixture [21] resulting in a local disturbance of the atomic density which oscillated in an harmonic trap much more slowly than predicted for a dark soliton [22][23][24]. In a recent experiment [25] it has been shown that the state created by means of the phase imprinting method evolves very quickly to a so-called vortex soliton -see also theoretical analysis in Refs. [26,27]. In the present article we show that, in order to create a state of a superfluid Fermi gas with symmetry properties identical to those of a stationary dark soliton solutio...