Double Schubert polynomials for the classical Lie groups

For each infinite series of the classical Lie groups of type B, C or D, we introduce a family of polynomials parametrized by the elements of the corresponding Weyl group of infinite rank. These polynomials represent the Schubert classes in the equivariant cohomology of the corresponding flag variety. They satisfy a stability property, and are a natural extension of the (single) Schubert polynomials of Billey and Haiman, which represent non-equivariant Schubert classes. When indexed by maximal Grassmannian elements of the Weyl group, these polynomials are equal to the factorial analogues of Schur Q- or P-functions defined earlier by Ivanov.