Adiabatic contraction wrapper potential¶
- class galpy.potential.AdiabaticContractionWrapperPotential(amp=1.0, pot=None, baryonpot=None, method='cautun', f_bar=0.157, rmin=None, rmax=50.0, ro=None, vo=None)[source]¶
AdiabaticContractionWrapperPotential: Wrapper to adiabatically contract a DM halo in response to the growth of a baryonic component. Use for example as:
dm= AdiabaticContractionWrapperPotential(pot=MWPotential2014[2],baryonpot=MWPotential2014[:2])
to contract the dark-matter halo in MWPotential2014 according to the baryon distribution within it. The basic physics of the adiabatic contraction is that a fraction f_bar of the mass in the original potential pot cools adiabatically to form a baryonic component baryonpot; this wrapper computes the resulting dark-matter potential using different approximations in the literature.
- __init__(amp=1.0, pot=None, baryonpot=None, method='cautun', f_bar=0.157, rmin=None, rmax=50.0, ro=None, vo=None)[source]¶
NAME:
__init__
PURPOSE:
initialize a AdiabaticContractionWrapper Potential
INPUT:
amp - amplitude to be applied to the potential (default: 1.)
pot - Potential instance or list thereof representing the density that is adiabatically contracted
baryonpot - Potential instance or list thereof representing the density of baryons whose growth causes the contraction
method= (‘cautun’) Type of adiabatic-contraction formula:
‘cautun’ for that from Cautun et al. 2020 (2020MNRAS.494.4291C),
‘blumenthal’ for that from Blumenthal et al. 1986 (1986ApJ…301…27B 1986ApJ…301…27B)
‘gnedin’ for that from Gnedin et al. 2004 (2004ApJ…616…16G)
f_bar= (0.157) universal baryon fraction; if None, calculated from pot and baryonpot assuming that at rmax the halo contains the universal baryon fraction; leave this at the default value unless you know what you are doing
rmin= (None) minimum radius to consider (default: rmax/2500; don’t set this to zero)
rmax= (50.) maximum radius to consider (can be Quantity)
ro, vo= standard unit-conversion parameters
OUTPUT:
(none)
HISTORY:
2021-03-21 - Started based on Marius Cautun’s code - Bovy (UofT)