Optical line shapes of molecular aggregates: hierarchical equations of motion method.