Lyα is a powerful astrophysical probe. Not only is it ubiquitous at high redshifts, it is also a resonant line, making Lyα photons scatter. This scattering process depends on the physical conditions of the gas through which Lyα propagates, and these conditions are imprinted on observables such as the Lyα spectrum and its surface brightness profile. In this work, we focus on a less-used observable capable of probing any scattering process: polarization. We implement the density matrix formalism of polarization into the Monte Carlo radiative transfer code tlac. This allows us to treat it as a quantum mechanical process where single photons develop and lose polarization from scatterings in arbitrary gas geometries. We explore static and expanding ellipsoids, biconical outflows, and clumpy multiphase media. We find that photons become increasingly polarized as they scatter and diffuse into the wings of the line profiles, making scattered Lyα polarized in general. The degree and orientation of Lyα polarization depends on the kinematics and distribution of the scattering H i gas. We find that it generally probes spatial or velocity space asymmetries and aligns itself tangentially to the emission source. We show that the mentioned observables, when studied separately, can leave similar signatures for different source models. We conclude by revealing how a joint analysis of the Lyα spectra, surface brightness profiles, and polarization can break these degeneracies and help us extract unique physical information on galaxies and their environments from their strongest, most prominent emission line.