Towards Efficient and Scalable Electronic Structure Methods for the Treatment of Relativistic Effects and Molecular Response