tests:sphtri:tidal
Differences
This shows you the differences between two versions of the page.
| Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
| tests:sphtri:tidal [2014/10/30 09:24] – justin | tests:sphtri:tidal [2022/10/24 12:26] (current) – external edit 127.0.0.1 | ||
|---|---|---|---|
| Line 1: | Line 1: | ||
| ===== Tidally stripped models ===== | ===== Tidally stripped models ===== | ||
| - | For these models, we set up a spherical two component (stars/dark matter) " | + | For these models, we set up a spherical two component (stars/dark matter) " |
| < | < | ||
| Line 17: | Line 17: | ||
| \Phi = \Phi_M + \Phi_L | \Phi = \Phi_M + \Phi_L | ||
| </ | </ | ||
| + | |||
| < | < | ||
| \Phi_M(R,z) = -\frac{GM_{\rm disk}}{\sqrt{R^2 + (a+\sqrt{b^2+z^2})^2}} | \Phi_M(R,z) = -\frac{GM_{\rm disk}}{\sqrt{R^2 + (a+\sqrt{b^2+z^2})^2}} | ||
| </ | </ | ||
| + | |||
| < | < | ||
| \Phi_L(R,z) = \frac{v_0^2}{2}\ln\left(R_c^2 + R^2 + \frac{z^2}{q_0^2}\right) + {\rm const.} | \Phi_L(R,z) = \frac{v_0^2}{2}\ln\left(R_c^2 + R^2 + \frac{z^2}{q_0^2}\right) + {\rm const.} | ||
| Line 28: | Line 30: | ||
| < | < | ||
| M_{\rm disk} = 5 \times 10^{10}\, | M_{\rm disk} = 5 \times 10^{10}\, | ||
| + | </ | ||
| + | < | ||
| v_0 = 220\,{\rm km/s} | v_0 = 220\,{\rm km/s} | ||
| + | </ | ||
| + | < | ||
| r_t = 8\,{\rm kpc} | r_t = 8\,{\rm kpc} | ||
| + | </ | ||
| + | < | ||
| a = 4\,{\rm kpc} | a = 4\,{\rm kpc} | ||
| + | </ | ||
| + | < | ||
| b = 0.5\,{\rm kpc} | b = 0.5\,{\rm kpc} | ||
| + | </ | ||
| + | < | ||
| q_0 = 1.0 | q_0 = 1.0 | ||
| </ | </ | ||
| - | Four orbits were chosen as follows; each aligned at 40 degrees to the x-z plane; each evolved for 5 Gyrs. | + | Four orbits were chosen as follows; each aligned at 40 degrees to the x-z plane; each evolved for 5 Gyrs using a version of the [[http:// |
| ^ Orbit ^ x(kpc) ^ y(kpc) ^ z(kpc) ^ vx(km/s) ^ vy(km/s) ^ vz(km/s) ^ [peri, | ^ Orbit ^ x(kpc) ^ y(kpc) ^ z(kpc) ^ vx(km/s) ^ vy(km/s) ^ vz(km/s) ^ [peri, | ||
| Line 50: | Line 62: | ||
| The following give data files for each. These include .txt files of **all** stars for the final disrupted satellite centred on the origin. This amounts to three data sets if using only projected data (i.e. one along each of the x,y,z projections). In addition, the raw Gadget data files are included. These have not been centred. PyNbody python tools for loading and manipulating these Gadget data can be found [[https:// | The following give data files for each. These include .txt files of **all** stars for the final disrupted satellite centred on the origin. This amounts to three data sets if using only projected data (i.e. one along each of the x,y,z projections). In addition, the raw Gadget data files are included. These have not been centred. PyNbody python tools for loading and manipulating these Gadget data can be found [[https:// | ||
| - | * .txt files :: x(kpc), | + | * .txt files :: x(kpc), |
| - | * .dat files (Gadget) :: < | + | * _err.txt files :: The same as above but assuming 2km/s Gaussian velocity errors. |
| + | * .dat files (Gadget) :: < | ||
| + | * Files are numbered by the sampling e.g. _10000 is 10,000 stars and by the random draw from the full distribution e.g. _10000_0. **As with the spherical mocks, you should model by default including the 2km/s errors and using the _0 then _1, _2, _3 files in order**. (i.e. if you can only afford to run a single tidal mock for each orbit, then you should run the _10000_0_err file projected along the z-axis.) | ||
| ^ Orbit ^ Data files ^ | ^ Orbit ^ Data files ^ | ||
| - | | Orbit1 | {{: | + | | Orbit1 | {{: |
| - | | Orbit2 | {{: | + | | Orbit2 | {{: |
| - | | Orbit3 | {{: | + | | Orbit3 | {{: |
| - | | Orbit4 | {{: | + | | Orbit4 | {{: |
tests/sphtri/tidal.1414661070.txt.gz · Last modified: 2022/10/24 12:26 (external edit)