def qD_func(g, a):
    """
    Parameter
    ---------
    g   : float
        ratio between exchange and dipolar interaction (supposed g << 1)
    a   : float
        lattice constant of the cubic unit cell

    Return
    ------
    qD  : float
        dipolar wavevector of the system
    """
    return g**0.5 / a