 solitons. As opposed to the bright soliton case explored in I, however, the real and imaginary parts of $p_j$ are not independent because of the constraints (3.2c). Actually, it may be parameterized either by the velocity of the $j$th soliton or by a single angular variable, as will see in section 4. An auxiliary variable $\tau$ introduced in (3.2a) will be used conveniently in performing the proof of (2.10). It can be set to zero after all the calculations have been completed.

[**Remark 3.1.**]{} The tau function $g$ given by (3.1b) is represented by the determinant of an $(N+1)\times (N+1)$ matrix. It can be rewritten by the determinant of an $N\times N$ matrix. To show this, we multiply the $(N+1)$th column of $g$ by $z_{k,t}^*/\rho^2$ and subtract it from the $k$th column for $k=1, 2,..., N$ to obtain $$g=\left|\left(\delta_{jk}-{\kappa+{\rm i}p_k^*\over p_j+p_k^*}\,{p_j\over p_k^*}\,z_jz_k^*\right)_{1\leq j,k\leq N}\right|. \eqno(3.3)$$ Although the tau function from (3.1b) is used in the proof of the dark $N$-soliton solution, an equivalent form (3.3) will be employed in section 4 to analyze the structure of the solution.

[**Remark 3.2.**]{} The complex parameters $p_j$ subjected to the constraints (3.2c) exist only if the condition $\kappa(1+\kappa\rho^2)>0$ is satisfied, as confirmed easily by putting $p_j=a_j+{\rm i}b_j$ with real $a_j$ and $b_j$. We will show in section 3.6 that this condition is closely related to the stability of the plane wave solution of the FL equation.

[*3.2. Notation and basic formulas for determinants*]{}

Before entering into the proof of the dark $N$-soliton solution, we first define the matrices associated with the dark $N$-soliton solution and then provide some basic formulas for determinants
solitons. As opposed to the bright soliton case explored in I, however, the veridical and fanciful parts of $ p_j$ are not independent because of the constraint (3.2c). Actually, it may be parameterized either by the speed of the $ j$th soliton or by a single angular variable, as will see in part 4. An auxiliary variable $ \tau$ precede in (3.2a) will be used conveniently in do the proof of (2.10). It can be set to zero after all the calculation have been complete. 

 [ * * Remark 3.1. * * ] { } The tau function $ g$ give by (3.1b) is constitute by the antigenic determinant of an $ (N+1)\times (N+1)$ matrix. It can be rewritten by the determinant of an $ N\times N$ matrix. To show this, we reproduce the $ (N+1)$th column of $ g$ by $ z_{k, t}^*/\rho^2 $ and subtract it from the $ k$th column for $ k=1, 2, ..., N$ to obtain $ $ g=\left|\left(\delta_{jk}-{\kappa+{\rm i}p_k^*\over p_j+p_k^*}\,{p_j\over p_k^*}\,z_jz_k^*\right)_{1\leq j, k\leq N}\right|. \eqno(3.3)$$ Although the tau function from (3.1b) is practice in the proof of the black $ N$-soliton solution, an equivalent form (3.3) will be employed in section 4 to analyze the structure of the solution. 

 [ * * Remark 3.2. * * ] { } The complex parameters $ p_j$ subjugate to the constraints (3.2c) exist only if the condition $ \kappa(1+\kappa\rho^2)>0 $ is satisfied, as confirmed easily by putting $ p_j = a_j+{\rm i}b_j$ with real $ a_j$ and $ b_j$. We will show in section 3.6 that this condition is closely related to the stability of the plane wave solution of the FL equation. 

 [ * 3.2. Notation and basic formula for determinants * ] { } 

 Before entering into the validation of the blue $ N$-soliton solution, we foremost define the matrices associated with the dark $ N$-soliton solution and then provide some basic formulas for determinants
 sopitons. As opposed to the bright soliton case eeplored in I, hodever, the real and imaginary pqrts if $p_j$ are not independdnt becaude of thw coistraints (3.2c). Actuemly, it may be iaramztxrized either bi the velocidy of the $j$th vouicon or by a single angular variable, ws will sfe in section 4. An stxiljary variable $\tau$ introduced in (3.2a) sill be used convenirntly in performing the prlof lf (2.10). It can be set ho zero aftgd ajo the calculxtions havt yeen complejed.

[**Remark 3.1.**]{} The tau function $g$ gixen bv (3.1b) is reprgssnhgd by the deverminwnt of an $(N+1)\tlkes (N+1)$ katrix. Ot can be rewrltten by the determinant of ai $N\times N$ matrix. To show thiv, ce multiply the $(N+1)$th cilymn ox $g$ ty $z_{y,r}^*/\rhu^2$ ahd shbtrach iv from the $i$th column dor $k=1, 2,..., N$ to obtain $$b=\lqdt|\left(\delta_{jk}-{\iappa+{\ri y}p_k^*\over p_j+p_k^*}\,{p_j\over p_k^*}\,z_jz_k^*\right)_{1\leq j,k\leq N}\rjght|. \eqno(3.3)$$ Although the rau function from (3.1b) id used in the proof of the dark $N$-soliton solution, an equivdlent worn (3.3) dulp be employed in section 4 to analyze the strusfute of the solutijn.

[**Remark 3.2.**]{} Tne cpiplex parametgrs $p_j$ subnected to the conshraints (3.2c) exust only yf tne condition $\kappa(1+\kappa\rho^2)>0$ us satisfied, qs confirmed easilv by putting $p_j=a_k+{\rm i}n_j$ with real $a_j$ and $b_j$. Ce wilm show in sfction 3.6 tgxt this conditiov ix wlosely related to the stafility of the plane dave solutyon of the FL equation.

[*3.2. Notation anf basnc fosmulas for determinants*]{}

Before entering invp the proof pf thv dark $N$-sjlitok solution, we fyrst define thg matricef assuciated winh the dack $N$-soliton folution and djen provide vome basyc firmuoas for aeterminants
solitons. As opposed to the bright soliton in however, the and imaginary parts because the constraints (3.2c). it may be either by the velocity of the soliton or by a single angular variable, as will see in section 4. auxiliary variable $\tau$ introduced in (3.2a) will be used conveniently in performing the of It be to zero after all the calculations have been completed. [**Remark 3.1.**]{} The tau function $g$ given (3.1b) is represented by the determinant of an (N+1)$ matrix. It can rewritten by the determinant of $N\times matrix. To this, multiply $(N+1)$th column of by $z_{k,t}^*/\rho^2$ and subtract it from the $k$th column for $k=1, 2,..., N$ to obtain $$g=\left|\left(\delta_{jk}-{\kappa+{\rm i}p_k^*\over p_k^*}\,z_jz_k^*\right)_{1\leq j,k\leq Although the function (3.1b) used in the the dark $N$-soliton solution, an equivalent be employed in section 4 to analyze the of the [**Remark 3.2.**]{} The complex parameters $p_j$ to the constraints (3.2c) exist only if the $\kappa(1+\kappa\rho^2)>0$ is satisfied, as confirmed easily by putting $p_j=a_j+{\rm i}b_j$ with real $a_j$ and $b_j$. show in section 3.6 this condition is related the of plane wave of the FL equation. [*3.2. Notation and basic formulas for determinants*]{} entering into the proof of the dark $N$-soliton solution, we the associated with the $N$-soliton solution and then some formulas for determinants
 solitons. As opposed to the briGht soliton Case eXplOreD iN I, hoWeveR, the real and imaGInarY parts of $p_j$ are not indepeNdent BeCAuse OF tHe conStraintS (3.2C). ACTUalLy, It May Be PArAmeteRizEd eitheR by the veloCitY oF the $j$th solitON oR by a single AngUlar variable, As wIll see In SecTIon 4. An AuxIliarY variaBLe $\tau$ iNtroduced In (3.2A) Will be USed convENIeNtly In performing the prOOf OF (2.10). It can be set to zEro aftEr ALl THE caLcuLations havE bEen coMPleted.

[**REMaRK 3.1.**]{} tHe tAU function $g$ givEn by (3.1b) is reprESenTed by tHe DetERminanT of an $(n+1)\tIMes (n+1)$ matrix. It caN be rEwritten bY the deTErminanT Of an $N\tiMes N$ maTriX. To Show THiS, wE muLtIPly THe $(n+1)$th COluMn of $g$ by $z_{K,t}^*/\RhO^2$ and sUbtrACT IT froM thE $k$th ColumN for $k=1, 2,..., N$ to obtaiN $$g=\lEft|\lEFt(\dElta_{jK}-{\kappA+{\rm i}P_k^*\Over p_J+p_k^*}\,{p_j\oVer p_k^*}\,Z_jZ_k^*\right)_{1\leq j,k\leq n}\rigHt|. \eqno(3.3)$$ AltHouGh The TaU funcTIon froM (3.1b) iS usEd in the Proof of THe dArK $n$-SOlIton solution, an equiVaLENt Form (3.3) will Be emplOYeD iN Section 4 tO aNalYze tHE StrucTure OF tHe solutiOn.

[**RemaRK 3.2.**]{} THe Complex PaRameteRs $P_j$ sUbjEcted TO the ConstrAints (3.2c) exIst onLY if the conditioN $\Kappa(1+\kappa\rho^2)>0$ IS sATIsFIed, aS coNfirmed easiLy by PUttiNg $p_j=A_J+{\rM i}b_J$ With rEal $a_j$ AnD $B_j$. wE will show in section 3.6 tHaT this cOnditIon is closely rElated to thE STAbility oF the PLaNE wave solution oF the Fl equation.

[*3.2. NOTation anD basiC formulaS for deterMINants*]{}

BefOre EntEriNg iNTO tHe proof of the dARK $N$-soLiTon soluTioN, we firsT deFinE thE maTrIces assocIated witH tHe DaRk $n$-soLiton SOlution aNd TheN pRovIde soME basic FormuLas fOr DeTErmInants
 solitons. As opposed to t he brightsolit oncas eexpl ored  in I, however ,  the  real and imaginary pa rts o f$ p_j$ ar e not  indepe n de n t  be ca us e o ft he  cons tra ints (3 .2c). Actu all y,  it may be p a ra meterizedeit her by the v elo city o fthe $j$th  so liton  or by a sing le angula rv ariabl e , as wi l l  s ee i n section 4. An a u xi l iary variable$\tau$  i n tr o d uce d i n (3.2a) w il l beu sed con v en i e n tly in performing  the proofo f ( 2.10).  I t c a n be s et to  z e roafter all t he c alculatio ns hav e  been c o mpleted .

[** Rem ark  3.1 . ** ]{ } T he tau fu nct i on$g$ give nby  (3.1 b) i s r e pres ent ed b y the  determinantofan $ ( N+1 )\tim es (N +1)$  m atrix . It c an be  r ewritten by the  det erminantofan  $N \t imesN $ matr ix.  To  show t his, we mul ti p l y  t he $(N+1)$th colum no f  $ g$ by $z _{k,t} ^ */ \r h o^2$ and  s ubt ract i t fro m th e  $ k$th col umn fo r  $ k= 1, 2,.. .,  N$ to  o bta in$$g=\ l eft| \left( \delta_{ jk}-{ \ kappa+{\rm i}p _ k^*\over p_j+ p _k ^ * }\ , {p_j \ov er p_k^*}\, z_jz _ k^*\ righ t )_ {1\ l eq j, k\leq  N } \r i ght|. \eqno(3.3)$$Al though  thetau functionfrom (3.1b ) i s used i n th e  p r oof of the dar k $N$ -soliton s o lution,an eq uivalent  form (3. 3 )  will be  em plo yed  in s ec tion 4 to ana l y ze t he  struct ure  of the  so lut ion .
[* *Remark 3 .2.**]{}  T he  c om ple x par a meters $ p_ j$su bje ctedt o theconst rain ts  ( 3 .2c ) exist on l y  ifth econd iti on  $\ka ppa( 1 +\k appa\rh o^2)>0$ i s s a tisf ie d,  as con firmed easily  b y putting$p _j= a_j+{\ r m  i}b_j$with real $a_j$ and $b_ j $. We w ill  show  insection 3 .6that t his condit ion is  clos el y r e l atedt o  t hest ability of t heplane  w avesolutio n of the FL equati o n.
[*3.2. Notat ion  and b as icf or m ula sf ord e terminants*]{}
Before en te r in g into the pro of  of the  dark $ N$-so l iton so lution, w e first d ef inet h e m atrices as sociated  with the dark$ N$ -soli ton  solut io n a nd th en pro v ide  some  basic  f ormula s for  d etermina nts
 solitons._As opposed_to the bright soliton_case explored_in_I, however,_the_real and imaginary_parts of $p_j$_are not independent because_of the constraints_(3.2c)._Actually, it may be parameterized either by the velocity of the $j$th soliton or_by_a single_angular_variable,_as will see in section_4. An auxiliary variable $\tau$_introduced in_(3.2a) will be used conveniently in performing the_proof_of (2.10). It_can be set to zero after all the calculations_have been completed.

[**Remark 3.1.**]{} The tau_function $g$ given_by_(3.1b)_is represented by the_determinant of an $(N+1)\times (N+1)$ matrix._It can be rewritten by the_determinant of an $N\times N$ matrix. To_show this, we multiply the $(N+1)$th_column of $g$ by $z_{k,t}^*/\rho^2$_and subtract_it from the $k$th column_for $k=1, 2,...,_N$ to_obtain $$g=\left|\left(\delta_{jk}-{\kappa+{\rm i}p_k^*\over_p_j+p_k^*}\,{p_j\over p_k^*}\,z_jz_k^*\right)_{1\leq j,k\leq N}\right|. \eqno(3.3)$$ Although_the tau function_from (3.1b) is used in the_proof_of the dark_$N$-soliton_solution,_an equivalent_form (3.3) will_be_employed in_section_4 to analyze the structure of_the_solution.

[**Remark 3.2.**]{} The complex parameters $p_j$ subjected_to the constraints (3.2c)_exist_only if the condition_$\kappa(1+\kappa\rho^2)>0$ is satisfied, as confirmed_easily by putting $p_j=a_j+{\rm i}b_j$ with_real $a_j$_and $b_j$._We will show in section 3.6 that this condition is closely_related to the stability of the_plane wave solution of_the FL_equation.

[*3.2._Notation and basic_formulas_for determinants*]{}

Before_entering into the proof of the dark_$N$-soliton solution,_we first define the matrices associated_with the dark $N$-soliton_solution_and then provide some basic formulas_for determinants