We prove an L p (Ω) × L q (Ω) × L r (Ω) → L 1 (Ω × (0, ∞)) embedding for triples of elliptic operators in divergence form with complex coefficients and subject to mixed boundary conditions on Ω, and for triples of exponents p, q, r ∈ (1, ∞) mutually related by the identity 1/p+1/q +1/r = 1. Here Ω is allowed to be an arbitrary open subset of R d . Our assumptions involving the exponents and coefficient matrices are rather naturally expressed in terms of a condition known as p-ellipticity. Furthermore, we give applications to paraproducts and square functions associated with the corresponding operator semigroups.