def power_law(deg_sep, ALPHA=1.0):
    """
    Applies a power law model.

    WEIGHT = DEG_SEP^(-ALPHA)
    ALPHA: Parameter to control power law decay rate. 
    """
    return (deg_sep + 1) ** (- ALPHA)