 evolution of the star and [brown dwarf ]{}population within the cluster from their birth to the age of the Pleiades depending on their initial kinematics.

Note that we only present the results for one set of initial conditions but we also performed other simulations using different numbers of objects $N$, which means in particular using various seed numbers to initialize the cluster model. The results are similar to those described in the following sections.

Results and discussion {#results}
======================

Characteristic radii {#radii}
--------------------

Over the timescale of the simulations the half mass radius does not change much through the simulations. It goes up to $\sim3.2-3.4$ pc after 12 crossing times which is consistent with observational results obtained for the Pleiades ($r_{h}=3.6$ pc; Pinfield et al. 1998). This result does not depend on the initial substellar velocity dispersion which is indeed expected considering the fact that [brown dwarf ]{}represent only 3% of the cluster mass.

A nominal core radius $r_{c}$ is calculated in [Nbody2]{} using the definition of the density radius given by Casertano & Hut (1985) slightly modified in order to obtain a convergent result using a smaller sample ($n\simeq N/2$). It is determined by the rms expression (Aarseth 2001) $$r_{c}= \left( \frac{\sum_{i=1}^{n} |\mathbf{r}_{i} -
\mathbf{r}_{d}|^{2} \, \rho_{i}^{2}}{\sum \rho_{i}^{2}} \right)^{1/2},$$ where $\mathbf{r}_{i}$ is the three-dimensional position vector of the $i$th star and $\mathbf{r}_{d}$ denotes the coordinates of the density centre. The density estimator $\rho_{i}=3\,M_{5}/(4\pi r_{6}^{3})$ is defined with respect to the sixth nearest particle $r_{6}$ and takes into account the total mass ($M_{5}$) of the five nearest neighbours. \[See Casertano & Hut (1985) for definitions of density centre and density estimator.\] In our simulations, $r_{c}=1.3$ pc at $t=0$ and $r_{c}\sim0.8$ pc at $t=12\,t_{cr}\simeq 120$ Myr. However, to be compared to the observational core
evolution of the star and [ brown dwarf ] { } population within the cluster from their parturition to the old age of the Pleiades depending on their initial kinematics. 

 Note that we only stage the results for one set of initial circumstance but we also performed early model using unlike numbers of objects $ N$, which means in particular use various semen numbers to initialize the bunch model. The results are similar to those report in the following sections. 

 consequence and discussion { # results } 
 = = = = = = = = = = = = = = = = = = = = = = 

 Characteristic radius { # radii } 
 -------------------- 

 Over the timescale of the simulations the half mass radius does not change a lot through the simulations. It goes up to $ \sim3.2 - 3.4 $ pc after 12 crossing times which is consistent with experimental consequence obtained for the Pleiades ($ r_{h}=3.6 $ pc; Pinfield et al. 1998). This result does not depend on the initial substellar velocity dispersion which is indeed expected considering the fact that [ brown dwarf ] { } typify only 3% of the cluster mass. 

 A nominative effect spoke $ r_{c}$ is calculated in [ Nbody2 ] { } using the definition of the density radius given by Casertano & Hut (1985) slightly modified in order to receive a convergent result using a smaller sample ($ n\simeq N/2 $). It is determined by the rms expression (Aarseth 2001) $ $ r_{c}= \left (\frac{\sum_{i=1}^{n } |\mathbf{r}_{i } - 
 \mathbf{r}_{d}|^{2 } \, \rho_{i}^{2}}{\sum \rho_{i}^{2 } } \right)^{1/2},$$ where $ \mathbf{r}_{i}$ is the three - dimensional place vector of the $ i$th star and $ \mathbf{r}_{d}$ denotes the coordinates of the density centre. The concentration estimator $ \rho_{i}=3\,M_{5}/(4\pi r_{6}^{3})$ is define with esteem to the sixth nearest particle $ r_{6}$ and takes into account the total mass ($ M_{5}$) of the five nearest neighbor. \[See Casertano & Hut (1985) for definitions of density centre and density estimator.\ ] In our pretense, $ r_{c}=1.3 $ pc at $ t=0 $ and $ r_{c}\sim0.8 $ personal computer at $ t=12\,t_{cr}\simeq 120 $ Myr. However, to be compared to the observational core
 evllution of the star and [nrown dwarf ]{}population xithin fhe clusger from their birth to the ege if tht Pleiades dependine on theig initial kintmatics.

Note that xs only ixesenf the cesults for one set of inhtial conditiots bbt we also performed other simulatiogs usinb fifferent numbgrs og objsbtw $N$, which means in particular using narious seed numbrrs to initialize the clusher lodel. The results wre similar to eyose describdd in the following sedtions.

Results and discussion {#rerults}
======================

Eharacterisjnx rwgii {#radii}
--------------------

Ovxr the timescale on the shmulatipns the half msss raeius does not change kuch through the symulationv. Nt goes up to $\sim3.2-3.4$ pc qfrer 12 wrosving rimds shmch is cojsiatent with observatiobal results obtainec sir the Pleiadss ($r_{h}=3.6$ [c; Pinfield et al. 1998). This result does not dtpend on the initial substeloar velocity dispersiln which ys indeed expected considering the fact that [browt dwacf ]{}rekrcfdbt only 3% of the cluster mass.

A nominal core radita $t_{c}$ is calculated in [Nbody2]{} uxijg jhe definition of thz dsnsity radius givej by Cafertabo & Hut (1985) fligntly modified in order to ovtain a convvrgebt result using a dmaller sam'le ($n\someq M/2$). It is determined by tke rms expression (Aarseth 2001) $$f_{c}= \left( \frac{\sum_{i=1}^{v} |\mstvbf{r}_{i} -
\maufbf{r}_{d}|^{2} \, \rho_{i}^{2}}{\sum \rhj_{i}^{2}} \right)^{1/2},$$ xhere $\mathbf{f}_{i}$ ix the ehree-dimendional position vector ov the $i$dh star anf $\mathbf{r}_{d}$ denotes the coordinavxs of the denxidy bentre. Thz denslty estimator $\rro_{i}=3\,M_{5}/(4\pi r_{6}^{3})$ is deyined wich resoect to thv sixth nxarest partisle $r_{6}$ and takad into accouit the toeal nass ($M_{5}$) of tfd five nearest neighbougs. \[See Casertano & Hut (1985) fpr asfinitions of dzusuty centre and cenritr vstmmatow.\] In our simunatiuns, $t_{c}=1.3$ pc at $t=0$ and $r_{g}\sio0.8$ pc at $t=12\,t_{cr}\simeq 120$ Myr. Hmwevsr, to be compared yo the obsgrvationaj core
evolution of the star and [brown dwarf the from their to the age their kinematics. Note that only present the for one set of initial conditions we also performed other simulations using different numbers of objects $N$, which means particular using various seed numbers to initialize the cluster model. The results are to described the sections. Results and discussion {#results} ====================== Characteristic radii {#radii} -------------------- Over the timescale of the simulations half mass radius does not change much through simulations. It goes up $\sim3.2-3.4$ pc after 12 crossing which consistent with results for Pleiades ($r_{h}=3.6$ pc; et al. 1998). This result does not depend on the initial substellar velocity dispersion which is indeed considering the [brown dwarf only of cluster mass. A radius $r_{c}$ is calculated in [Nbody2]{} of the density radius given by Casertano & (1985) slightly in order to obtain a convergent using a smaller sample ($n\simeq N/2$). It is by the rms expression (Aarseth 2001) $$r_{c}= \left( \frac{\sum_{i=1}^{n} |\mathbf{r}_{i} - \mathbf{r}_{d}|^{2} \, \rho_{i}^{2}}{\sum \rho_{i}^{2}} $\mathbf{r}_{i}$ is the three-dimensional vector of the star $\mathbf{r}_{d}$ the of the centre. The density estimator $\rho_{i}=3\,M_{5}/(4\pi r_{6}^{3})$ is defined with respect to sixth nearest particle $r_{6}$ and takes into account the total of five nearest neighbours. Casertano & Hut (1985) definitions density centre and density our $r_{c}=1.3$ and pc $t=12\,t_{cr}\simeq 120$ Myr. However, be compared to the observational
 evolution of the star and [browN dwarf ]{}popuLatioN wiThiN tHe clUsteR from their birtH To thE age of the Pleiades depenDing oN tHEir iNItIal kiNematicS.

noTE ThaT wE oNly PrESeNt the ResUlts for One set of inItiAl Conditions buT We Also perforMed Other simulatIonS using DiFfeREnt nuMbeRs of oBjects $n$, Which mEans in parTiCUlar usINg varioUS SeEd nuMbers to initialize THe CLuster model. The ResultS aRE sIMIlaR to Those descrIbEd in tHE followINg SECTioNS.

Results and diScussion {#resULts}
======================

characTeRisTIc radiI {#radiI}
--------------------

OVEr tHe timescale Of thE simulatiOns the HAlf mass RAdius doEs not cHanGe mUch tHRoUgH thE sIMulATiOns. iT goEs up to $\siM3.2-3.4$ pC aFter 12 cRossING TImes WhiCh is ConsiStent with obseRvaTionAL reSults ObtaiNed fOr The PlEiades ($R_{h}=3.6$ pc; PInField et al. 1998). This reSult Does not dePenD oN thE iNitiaL SubsteLlaR veLocity dIspersiON whIcH IS InDeed expected considErING tHe fact thAt [browN DwArF ]{}RepresenT oNly 3% Of thE CLusteR masS.

a nOminal coRe radiUS $r_{C}$ iS calculAtEd in [NbOdY2]{} usIng The deFInitIon of tHe densitY radiUS given by CasertANo & Hut (1985) slightly MOdIFIeD In orDer To obtain a coNverGEnt rEsulT UsIng A SmallEr samPlE ($N\sIMeq N/2$). It is determined bY tHe rms eXpresSion (Aarseth 2001) $$r_{c}= \Left( \frac{\suM_{I=1}^{N} |\Mathbf{r}_{i} -
\MathBF{r}_{D}|^{2} \, \Rho_{i}^{2}}{\sum \rho_{i}^{2}} \rigHt)^{1/2},$$ wheRe $\mathbf{r}_{i}$ IS the threE-dimeNsional pOsition veCTOr of the $i$Th sTar And $\MatHBF{r}_{D}$ denotes the coORDinaTeS of the dEnsIty centRe. THe dEnsIty EsTimator $\rhO_{i}=3\,M_{5}/(4\pi r_{6}^{3})$ is DeFiNeD wIth RespeCT to the siXtH neArEst PartiCLe $r_{6}$ and Takes Into AcCoUNt tHe total MAsS ($m_{5}$) Of thE fIvE neaResT nEighbOurs. \[sEe CAsertanO & Hut (1985) for deFinITionS oF dEnsity cEntre and densiTy Estimator.\] IN oUr sImulatIONs, $r_{c}=1.3$ pc at $T=0$ and $r_{c}\sim0.8$ pc at $t=12\,t_{cr}\simeq 120$ MYR. HoweveR, to Be comPareD to the obsErvAtionaL coRE
 evolution of the star and  [brown dw arf ] {}p opu la tion  wit hin the cluste r  fro m their birth to the a ge of  t h e Pl e ia des d ependin g  o n the ir  i nit ia l  k inema tic s.

Not e that weonl ypresent ther es ults for o neset of initi alcondit io nsb ut we  al so pe rforme d  other  simulati on s  using differe n t  n umbe rs of objects $N$ ,  w h ich means in p articu la r  u s i ngvar ious seednu mbers to init i al i z e  th e  cluster mode l. The resu l tsare si mi lar to tho se de sc r ibe d in the fo llow ing secti ons.
R esultsa nd disc ussion  {# res ults } 
= == === == = === = == === = === 

Charac te ri sticradi i { # radi i}---- ----- -----------
Ove r th e  ti mesca le of  the  s imula tionsthe h al f mass radius d oesnot chang e m uc h t hr ought he sim ula tio ns. Itgoes up to$\ s i m 3. 2-3.4$ pc after 12  c r o ss ing time s whic h  i sc onsisten twit h ob s e rvati onal re sults ob tained fo rthe Ple ia des ($ r_ {h} =3. 6$ pc ;  Pin fieldet al. 1 998). This result do e s not dependo nt h ei niti alsubstellarvelo c itydisp e rs ion which  is i nd e ed expected considerin gthe fa ct th at [brown dwa rf ]{}repr e s e nt only3% o f  t h e cluster mass .

Anominal co r e radius  $r_{ c}$ is c alculated i n [Nbody 2]{ } u sin g t h e  d efinition oft h e de ns ity rad ius  givenbyCas ert ano  &  Hut (198 5) sligh tl ymo di fie d ino rder toob tai na c onver g ent re sultusin gas mal ler sam p le ( $n\s im eq  N/2 $).  I t isdete r min ed by t he rms ex pre s sion  ( Aa rseth 2 001) $$r_{c}=  \ left( \fra c{ \su m_{i=1 } ^ {n} |\ma thbf{r}_{i} -
\mathbf{r } _{d}|^{ 2}\, \r ho_{ i}^{2}}{\ sum  \rho_ {i} ^ {2}} \ right) ^{1/2 }, $$w h ere $ \ m at hbf {r }_{i}$ ist h e t hree- di mens ional p osition vector oft he$i$th star an d $ \mat h b f{ r}_ { d} $  de no t est h e coordinates o f the dens it y  c entre. The den si ty esti mator $ \rho_ { i}=3\,M _{5}/(4\p i r_{6}^{ 3} )$ i s def ined withrespectto the si x th ne a re st pa rti cle $r _{ 6}$  andtakesi nto  acco unt th etotalmass($ M_{5}$)of the five nearest nei ghbour s. \[ See  Casertan o & Hut  (1985) f or d efinitions  of  de nsity  ce n tre a nd d e ns ity estim ator . \] In our si mul a t io ns, $r_{c}= 1 . 3 $ p c at$t= 0 $ and$r_{ c}\sim0.8$ pc at$ t=12\,t_{cr}\s imeq 1 20$  My r . Ho we ver, to be com par ed t o the ob se rvational c ore
 evolution_of the_star and [brown dwarf_]{}population within_the_cluster from_their_birth to the_age of the_Pleiades depending on their_initial kinematics.

Note that_we_only present the results for one set of initial conditions but we also performed_other_simulations using_different_numbers_of objects $N$, which means_in particular using various seed_numbers to_initialize the cluster model. The results are similar_to_those described in_the following sections.

Results and discussion {#results}
======================

Characteristic radii {#radii}
--------------------

Over the_timescale of the simulations the half_mass radius does_not_change_much through the simulations._It goes up to $\sim3.2-3.4$ pc_after 12 crossing times which is_consistent with observational results obtained for the_Pleiades ($r_{h}=3.6$ pc; Pinfield et al._1998). This result does not_depend on_the initial substellar velocity dispersion_which is indeed_expected considering_the fact that_[brown dwarf ]{}represent only 3% of_the cluster mass.

A_nominal core radius $r_{c}$ is calculated_in_[Nbody2]{} using the_definition_of_the density_radius given by_Casertano_& Hut_(1985)_slightly modified in order to obtain_a_convergent result using a smaller sample ($n\simeq_N/2$). It is determined_by_the rms expression (Aarseth_2001) $$r_{c}= \left( \frac{\sum_{i=1}^{n} |\mathbf{r}_{i}_-
\mathbf{r}_{d}|^{2} \, \rho_{i}^{2}}{\sum \rho_{i}^{2}} \right)^{1/2},$$ where_$\mathbf{r}_{i}$ is_the three-dimensional_position vector of the $i$th star and $\mathbf{r}_{d}$ denotes the coordinates_of the density centre. The density_estimator $\rho_{i}=3\,M_{5}/(4\pi r_{6}^{3})$ is_defined with_respect_to the sixth_nearest_particle $r_{6}$_and takes into account the total mass_($M_{5}$) of_the five nearest neighbours. \[See Casertano_& Hut (1985) for_definitions_of density centre and density estimator.\]_In our simulations, $r_{c}=1.3$ pc at_$t=0$ and $r_{c}\sim0.8$ pc at_$t=12\,t_{cr}\simeq_120$_Myr. However, to be compared_to the observational core