Abstract. We define a family of quasi-isometry invariants of groups called higher divergence functions, which measure isoperimetric properties "at infinity." We give sharp upper and lower bounds on the divergence functions for right-angled Artin groups, using different pushing maps on the associated cube complexes. In the process, we define a class of RAAGs we call orthoplex groups, which have the property that their Bestvina-Brady subgroups have hard-to-fill spheres. Our results give sharp bounds on the higher Dehn functions of Bestvina-Brady groups, a complete characterization of the divergence of geodesics in RAAGs, and an upper bound for filling loops at infinity in the mapping class group.