mcfacts.setup.diskstars_hillspheremergers

mcfacts.setup.diskstars_hillspheremergers.hill_sphere_orb_a(orbs_a_sorted, mass_threshold, smbh_mass, disk_radius_outer)

Find location of edge of Hill sphere for stars

Parameters:

  • orbs_a_sorted (numpy.ndarray) – Locations [r_{g,SMBH}] of stars in the disk with float type

  • mass_threshold (numpy.ndarray) – Minimum mass [M_sun] to not merge, as a function of orbs_a_sorted with float type

  • smbh_mass (float) – Mass [M_sun] of supermassive black hole

  • disk_radius_outer (float) – Outer radius [r_{g,SMBH}] of the disk

Returns:

delta_orbs_a – Locations [r_{g,SMBH}] of the Hill spheres for stars in the disk

Return type:

numpy.ndarray

mcfacts.setup.diskstars_hillspheremergers.mass_threshold_merger(disk_star_num, disk_stars_mass_min, smbh_mass, P_m, P_r, disk_stars_orb_a, disk_radius_outer)
disk_star_numint

number of stars in the initial draw

disk_stars_mass_minfloat

minimum mass [M_sun] considered

smbh_massfloat

mass of SMBH, M_sun

P_mfloat

exponent for mass cdf, assuming it is in the form P(> m_min) = (m_min/m)^P_m

P_rfloat

exponent for disk orb_aation cdf, assuming the form P(r) = (r_orb_aation/disk_radius_outer)^P_r

r_orb_aationnumpy array

semi-major axis of stellar orbit around SMBH, R_sun (for now?)

disk_radius_outerfloat

Outer radius [r_{g,SMBH}] of disk

Returns: mass_threshold : numpy.ndarray

minimum mass [M_sun] for stars to not merge

Notes

Eqn is (M_min^(P_m/(P_m - 1/3))) / ((3*M_smbh)^(1/3(P_m - 1/3))) * (N_star * P_r)^(1/(P_m - 1/3)) * (orb_a/R)^(P_r/(P_m - 1/3))