A Polarization-Decomposed Method for Simulating Inhomogeneous Birefringence in Laser-Interferometric Gravitational-Wave Detectors

arXiv:2606.08943v1 Announce Type: new Abstract: Birefringence in test mass substrates is an emerging limitation for current and future laser-interferometric gravitational-wave detectors, particularly as detectors move toward higher circulating power, cryogenic operation, and crystalline optical materials. Spatially varying birefringence alters both the polarization state and spatial mode content of the intracavity field, reducing interference contrast and coupling into length and alignment control signals. Accurate modeling of these effects is complicated by the fact that most frequency-domain simulation tools employ scalar modal propagation and lack native support for polarization and two-dimensional substrate maps. In this work, we present a practical and general method for simulating inhomogeneous birefringence without modifying existing simulation frameworks. The approach represents the two polarization components as independent scalar fields and introduces their coupling through an equivalent triple-Mach-Zehnder construction that reproduces the Jones matrix of a birefringent medium. We demonstrate the method using realistic birefringence maps of the KAGRA sapphire input test masses. The technique is compatible with any frequency-domain interferometer model and enables efficient birefringence studies for next-generation gravitational-wave detectors.