Differences

This shows you the differences between two versions of the page.

--- tests:collision:mock_data [2016/10/12 15:01] – v.henault-brunet
+++ tests:collision:mock_data [2022/10/24 12:27] (current) – external edit 127.0.0.1
@@ Line 1: / Line 1: @@
 ====== Mock data ======
+[[:tests:collision:mock_data:Challenge_1|Challenge 1]]: Single mass clusters, isolated and tidally limited.
-===== Challenge 4: Method comparison: N-body model of M4 =====
+[[:tests:collision:mock_data:Challenge_2|Challenge 2]]: Multi-mass clusters, no stellar evolution
-The mock data is a snapshots of large $N$-body simulation of the globular cluster M4 are [[http://datashare.is.ed.ac.uk/handle/10283/618?show=full|available here]]. Full lists of data for all particles (mass, position, velocity, stellar type, stellar radius, V magnitude, B-V colour) are given at almost 400 instants in the lifetime of the model. The simulation is described in the [[http://arxiv.org/abs/1409.5597]], which should be cited when making use of this data.
+[[:tests:collision:mock_data:Challenge_3|Challenge 3]]: Clusters in tidal fields with stellar evolution
+[[:tests:collision:mock_data:Challenge_4|Challenge 4]]: Comparison of mass modelling methods using N-body model of M4.
-We selected the snapshot at $t = 12023.9$ Myr and assume that radial/LOS velocities are available for all stars brighter than V=15 (N=1143). A distance modulus of 11.35 mag was adopted (i.e. $D=1862.1$ pc). The cluster properties are:
-  - $M = 64255.4\,M_\odot$
-  - $L_{V} = 34304.8\,L_{V,\odot}$
-  - $M/L_V=1.87$
-  - $r_{\rm h} = 3.21\,{\rm pc}$
-We do the comparison in 3 steps:
+The definition of the stellar types (KSTAR) used in the ''NBODY6'' snapshots can be found in the [[:tests:collision:mock_data:appendix|Appendix]].
-  - Surface brightness profile + RVs
-  - Surface brightness profile + RVs + PMs
-  - Surface brightness profile + RVs + PMs + Gaia
-The data for the 3 fitting steps :
-  - Step 1:
-    - Surface brightness: {{:data:M4_sb.dat.gz}}
-    - Discrete LOS velocities: {{:data:M4_rv.dat.gz}}
-    - LOS velocity dispersion profile: {{:data:M4_rv_disp.dat.gz}}
-  - Step 2:
-    - Discrete proper motions: {{:data:M4_pm_mock.dat.gz}}
-    - Proper motion dispersion profiles: {{:data:M4_pm_disp_mock.dat.gz}}
-For methods that allow it, you can also use constraints from a mock observed mass function at a projected radius of 300 arcseconds from the centre: {{:tests:collision:M4_mass_function.tab.gz}}
-A number density profile for stars brighter than the main-sequence turnoff (V < 17 mag; corresponding to a turnoff mass of $0.82 M_{\odot}$) is also available ({{:tests:collision:M4_number_density.dat.gz}}). This can be used as an alternative/complement to the surface brightness profile.
-The mass profile in projection and the $M/L_V$ ratio in projection can be found here:
-  - {{:data:m4_Sigma.dat.gz}}
-===== Challenge 1: single mass clusters =====
-==== Isolated models: ====
-The N-body models can be described as:
-  - Initial conditions: Plummer (1911), N = 32768, all stars the same mass
-  - No primordial binaries, no central black hole
-  - Isolation
-The data has the following format. Note that the first column can be used to recognise binaries (MN=2). The single components of the binaries are not given.
-^ $M\times N$ ^ $X$ ^ $Y$ ^ $Z$ ^ $V_x$ ^ $V_y$ ^ $V_z$ ^
-|  [NBODY] |  [NBODY]  |||    [NBODY]   |||
-  - {{:data:32k_logt4.tab.gz}} $T=10^4$
-  - {{:data:32k_logt5.tab.gz}} $T=10^5$
-  - {{:data:32k_logt6.tab.gz}} $T=10^6$
-  - {{:data:32k_logt7.tab.gz}} $T=10^7$
-==== Tidally limited models: ====
-The N-body models can be described as:
-  - Initial conditions: Plummer (1911), N = 65536, all stars the same mass
-  - No primordial binaries, no central black hole
-  - Circular orbit in a weak tidal field due to a point-mass galaxy with initially r_jacobi/r_h = 100
-The model was ran until complete dissolution (roughly 6e5 N-body times) with Sverre Aarseth's ''NBODY6''. Two-body relaxation drives the evolution: core collapse occurs at roughly T = 1.2e4 and the cluster expands until it fills the Roche-volume roughly half-mass the evolution (T = 3e5). More details about this run can be found [[http://adsabs.harvard.edu/abs/2012MNRAS.422.3415A|Alexander & Gieles (2012)]].
-Below are 3 snapshots at interesting moments of the evolution. The Heggie & Mathieu (1986) N-body units are used: G=M=r_vir=1 (i.e. the mass of individual stars is m=1/65536). The 6 columns contain:
-^ $X$ ^ $Y$ ^ $Z$ ^ $V_x$ ^ $V_y$ ^ $V_z$ ^
-|  [NBODY]  |||    [NBODY]   |||
-  - {{:data:PL_EQ_N64K_RJRH100_T012102.gz}} : In a core minimum just after core collapse [NEW: 19 Aug, 16:15]
-  - {{:PL_EQ_N64K_RJRH100_T013650.gz}} : In a core maximum
-  - {{:PL_EQ_N64K_RJRH100_T323790.gz}} : When ~75% of the stars is lost and the cluster is Roche-filling
-Update 14-Oct-2014: More snapshots for this model (roughly) equally spaced by $5\times10^4$ $N$-body times covering the entire life-cycle:
-  - {{:data:PL_EQ_N64K_RJRH100_T050010.gz}}
-  - {{:data:PL_EQ_N64K_RJRH100_T100020.gz}}
-  - {{:data:PL_EQ_N64K_RJRH100_T150020.gz}}
-  - {{:data:PL_EQ_N64K_RJRH100_T200020.gz}}
-  - {{:data:PL_EQ_N64K_RJRH100_T250020.gz}}
-  - {{:data:PL_EQ_N64K_RJRH100_T300030.gz}}
-  - {{:data:PL_EQ_N64K_RJRH100_T350520.gz}}
-  - {{:data:PL_EQ_N64K_RJRH100_T400020.gz}}
-  - {{:data:PL_EQ_N64K_RJRH100_T450510.gz}}
-  - {{:data:PL_EQ_N64K_RJRH100_T500010.gz}}
-  - {{:data:PL_EQ_N64K_RJRH100_T550500.gz}}
-Update: 29 Okt 2014: New version of the 10 snapshots above:
-  - Removed the individual components of binaries, and added the binary com pos and vel in the end of the file
-  - New first column with MxN = 1 for single stars and MxN = 2 for binaries
-  - New column (8) = 1 if r<rt
-  - New column (9) = 1 if E_Jacobi < E_crit
-Update 29-Nov-2014:
-   - Fixed bug in energy computation
-   - New column (8): phi (= specific potential)
-   - New column (9); E_J = jacobi energy (see e.g. Fukushige & Heggie (2000), below equation 3)
-   - Added top line with: N, rt, E_crit
-First line: N, rt, E_crit
-^ $M\timesN$ | $X$ ^ $Y$ ^ $Z$ ^ $V_x$ ^ $V_y$ ^ $V_z$ ^ $\phi$ ^ $E_J$ ^ r < rt ^ E<Ecrit ^
-| [NBODY] |  [NBODY]  |||    [NBODY]   |||[NBODY]|[NBODY]|0 or 1 | 0 or 1 |
-  - {{:data:PL_EQ_N64K_RJRH100_T050010.new.gz}}
-  - {{:data:PL_EQ_N64K_RJRH100_T100020.new.gz}}
-  - {{:data:PL_EQ_N64K_RJRH100_T150020.new.gz}}
-  - {{:data:PL_EQ_N64K_RJRH100_T200020.new.gz}}
-  - {{:data:PL_EQ_N64K_RJRH100_T250020.new.gz}}
-  - {{:data:PL_EQ_N64K_RJRH100_T300030.new.gz}}
-  - {{:data:PL_EQ_N64K_RJRH100_T350520.new.gz}}
-  - {{:data:PL_EQ_N64K_RJRH100_T400020.new.gz}}
-  - {{:data:PL_EQ_N64K_RJRH100_T450510.new.gz}}
-  - {{:data:PL_EQ_N64K_RJRH100_T500010.new.gz}}
-  - {{:data:PL_EQ_N64K_RJRH100_T550500.new.gz}}
-For the last snapshot a table with specific energy and the z-component of the specific angular momentum vector can be found here:
-  - {{:data:PL_EQ_N64K_RJRH100_T323790.EJZ.gz}} : When ~75% of the stars is lost and the cluster is Roche-filling
-Note: the initial Jacobi radius of this model was $r_{\rm J}= 78.17$, such that the angular frequency of the orbit is $\Omega = 8.354\times 10^{-4}$ and the critical energy $E_{\rm crit} = -7.469\times 10^{-3}$ at T=323790 .
-Illustration of the model evolution, moments of the snapshots are marked with dashed lines:
-{{:data:collisional_singlemass_radii.png?600}}
-Properties of the clusters:
-^ Cluster ^ Mass ^ $r_{\rm c}$ ^ $r_{\rm h}$ ^ $r_{\rm J}$ ^
-|1        | 0.975|$5.25\times 10^{-3}$ |1.143| 77.3|
-|2        | 0.953|$8.61\times 10^{-3}$ |1.334| 76.6|
-|3        | 0.238|$0.199$              |6.871| 48.3|
-===== Challenge 2: Multi-mass clusters (no stellar evolution) =====
-Data: Snapshots of simulations with a mass function.
-   - N <~= $10^5$, initial half-mass radius 2.25 pc, Henon isochrone model
-   - Initial half-mass relaxation time =~ 350 Myr
-   - no primordial binaries
-   - Galaxy = singular isothermal sphere with Vc = 220 kms/s
-   - Orbit: circular orbit at RG = 4 kpc
-The $10^5$ stars were evolved with SSE (Hurley et al. 2000) to an age of 12 Gyr assuming a metallicity of [Fe/H] = -2 **before** the N-body model was run. Then an assumption was made about the retention fraction of neutron stars (NSs) and black holes (BHs). The snapshots contain all the stars within the tidal radius.
-The properties of the snapshots are:
-^ Name ^ f_retention NSs ^ f_rentention BHs ^ age ^ N(NS)_remaining ^ N(BH)_remaining ^
-| cluster1 | 1 | 1 | 0.25 Gyr | 613 | 171 |
-| cluster1 | 1 | 1 | 0.5 Gyr  | 591 | 146 |
-| cluster1 | 1 | 1 | 0.75 Gyr | 559 | 134 |
-| cluster1 | 1 | 1 | 1 Gyr    | 548 | 119 |
-| cluster7 | 0.1 | 0.1 | 2 Gyr | 51 | 0 |
-Format of the snapshots are:
-^ $KSTAR$ ^ $m$^ $X$ ^ $Y$ ^ $Z$ ^ $V_x$ ^ $V_y$ ^ $V_z$ ^
-|   | [Msun] ||[pc]|    [km/s]   |||
-  - {{:data:gc1_1gyr.gz}} :
-  - {{:data:gc1_2gyr.gz}} :
-  - {{:data:gc1_3gyr.gz}} :
-  - {{:data:gc1_4gyr.gz}} :
-  - {{:data:gc1_5gyr.gz}} :
-  - {{:data:gc1_6gyr.gz}} :
-  - {{:data:gc1_7gyr.gz}} :
-  - {{:data:gc1_8gyr.gz}} :
-  - {{:data:gc1_9gyr.gz}} :
-  - {{:data:gc1_10gyr.gz}} :
-  - {{:data:gc1_11gyr.gz}} :
-No retention:
-  - {{:data:gc2_1gyr.gz}} :
-  - {{:data:gc2_2gyr.gz}} :
-  - {{:data:gc2_3gyr.gz}} :
-  - {{:data:gc2_4gyr.gz}} :
-  - {{:data:gc2_5gyr.gz}} :
-  - {{:data:gc2_6gyr.gz}} :
-  - {{:data:gc2_7gyr.gz}} :
-  - {{:data:gc2_8gyr.gz}} :
-  - {{:data:gc2_9gyr.gz}} :
-  - {{:data:gc2_10gyr.gz}} :
-  - {{:data:gc2_11gyr.gz}} :
-  - {{:data:gc2_12gyr.gz}} :
-The definition of the stellar types (KSTAR) used in the ''NBODY6'' snapshots can be found in the [[:tests:collision:appendix|Appendix]].
-===== Challenge 3. Clusters in tidal fields with stellar evolution =====
-(Simulations ran and kindly made available by Holger Baumgardt)\\
-Here we consider 2 clusters which are slightly more realistic:
-  - IC: King (1966) W_0 = 5 model, N = 131072, Kroupa (2001) mass function between 0.1-15 Msun (no black-holes).
-  - No primordial binaries, no central black hole, circular orbit in logarithmic halo with V = 220 km/s.
-  - Z = 0.001
-  - Stellar evolution and mass-loss according to Hurley et al. (2000, 2002)
-  - Two Galactocentric radii: 8.5 kpc and 15 kpc.
-Below are 2 snapshots at an age of roughly 10 Myr, 100 Myr, 1Gyr and 12 Gyr. The columns are the same as in Challenge 2.
-  - {{:data:W05_N131K_RG8.5_FEH-0.0_T10.gz}} UPDATED! Thursday August 22
-  - {{:data:W05_N131K_RG8.5_FEH-0.0_T100.gz}} UPDATED! Thursday August 22
-  - {{:data:W05_N131K_RG8.5_FEH-0.0_T1000.gz}} UPDATED! Thursday August 22
-  - {{:W05_N131K_RG8.5_FEH-0.0_T12000.gz}}
-  - {{:data:W05-N131K_RG15_FEH-0.0.T10.gz}} NEW! Tuesday August 20
-  - {{:data:W05-N131K_RG15_FEH-0.0.T100.gz}} NEW! Tuesday August 20
-  - {{:data:W05-N131K_RG15_FEH-0.0.T1000.gz}} NEW! Tuesday August 20
-  - {{:W05_N131K_RG15_FEH-0.0_T12000.gz}}
-Final list of snapshots used in Sollima et al. in prep, columns are:
-^ $m$    ^ $X$ ^ $Y$ ^ $Z$ ^ $V_x$ ^ $V_y$ ^ $V_z$ ^ kstar ^ logL ^ logTeff ^
-| [Msun] |[pc] |[pc] |[pc] | [km/s |[km/s] |[km/s] |       |      |         |
-  - {{:data:W5rh11.5R8.5-t0.gz}}
-  - {{:data:W5rh11.5R8.5-t10024.05.gz}}
-  - {{:data:W5rh11.5R8.5-t1019.394.gz}}
-  - {{:data:W5rh11.5R8.5-t11009.46.gz}}
-  - {{:data:W5rh11.5R8.5-t2004.809.gz}}
-  - {{:data:W5rh11.5R8.5-t3024.204.gz}}
-  - {{:data:W5rh11.5R8.5-t4009.618.gz}}
-  - {{:data:W5rh11.5R8.5-t5029.013.gz}}
-  - {{:data:W5rh11.5R8.5-t6014.427.gz}}
-  - {{:data:W5rh11.5R8.5-t7033.822.gz}}
-  - {{:data:W5rh11.5R8.5-t8019.236.gz}}
-  - {{:data:W5rh11.5R8.5-t9004.651.gz}}
-  - {{:data:W5rh1R8.5-t0.gz}}
-  - {{:data:W5rh1R8.5-t10020.gz}}
-  - {{:data:W5rh1R8.5-t1020.gz}}
-  - {{:data:W5rh1R8.5-t11010.gz}}
-  - {{:data:W5rh1R8.5-t120.gz}}
-  - {{:data:W5rh1R8.5-t12000.gz}}
-  - {{:data:W5rh1R8.5-t13020.gz}}
-  - {{:data:W5rh1R8.5-t2010.gz}}
-  - {{:data:W5rh1R8.5-t210.gz}}
-  - {{:data:W5rh1R8.5-t30.gz}}
-  - {{:data:W5rh1R8.5-t300.gz}}
-  - {{:data:W5rh1R8.5-t3000.gz}}
-  - {{:data:W5rh1R8.5-t4020.gz}}
-  - {{:data:W5rh1R8.5-t420.gz}}
-  - {{:data:W5rh1R8.5-t5010.gz}}
-  - {{:data:W5rh1R8.5-t510.gz}}
-  - {{:data:W5rh1R8.5-t60.gz}}
-  - {{:data:W5rh1R8.5-t6000.gz}}
-  - {{:data:W5rh1R8.5-t7020.gz}}
-  - {{:data:W5rh1R8.5-t8010.gz}}
-  - {{:data:W5rh1R8.5-t90.gz}}
-  - {{:data:W5rh1R8.5-t9000.gz}}
-The definition of the stellar types (KSTAR) used in the ''NBODY6'' snapshots can be found in the [[:tests:collision:appendix|Appendix]].
-===== Appendix =====
-==== Stellar types in ''NBODY'' ====
-There 23 possible stellar types (KSTAR) in ''NBODY''
-       Low main sequence (M < 0.7).
-       Main sequence.
-       Hertzsprung gap (HG).
-       Red giant.
-       Core Helium burning.
-       First AGB.
-       Second AGB.
-       Helium main sequence.
-       Helium HG.
-       Helium GB.
-       Helium white dwarf.
-       Carbon-Oxygen white dwarf.
-       Oxygen-Neon white dwarf.
-       Neutron star.
-       Black hole.
-       Massless supernova remnant.
-       Circularizing binary (c.m. value).
-       Circularized binary.
-       First Roche stage (inactive).
-       Second Roche stage.