Many functions on normalized FoldedLightCurve failed with `UnitConversionError`: `copy()`, `select_flux()`, `normalize()`, etc. about lightkurve HOT 2 OPEN

class FoldedLightCurve(LightCurve):
        def bin(
        self,
        time_bin_size=None,
        time_bin_start=None,
        time_bin_end=None,
        n_bins=None,
        aggregate_func=None,
        bins=None,
        binsize=None,
    ):
        if (hasattr(self, 'normalize_phase')) & (self.normalize_phase == True):
            with self._delay_required_column_checks():
                normalized_phase = self.time.copy()
                phase = TimeDelta(normalized_phase * self.period)
                self.remove_column("time")
                self.add_column(phase, name="time", index=0)


        *Regular bin operation from LightCurve*


        # Reset at the end of the function
        if (hasattr(self, 'normalize_phase')) & (self.normalize_phase == True):
            
            # Reset the 'time' column of the folded object to be the normalized value
            with self._delay_required_column_checks():
                phase = self.time.copy()
                normalized_phase = phase / self.period
                self.remove_column("time")
                self.add_column(normalized_phase, name="time", index=0)
                self.sort("time")

            # Reset the 'time' column of the folded and binned object to be the normalized value
            with returned_object._delay_required_column_checks():
                phase = returned_object.time.copy()
                normalized_phase = phase / self.period #returned_object.normalized_phase.copy()
                returned_object.remove_column("time")
                returned_object.add_column(normalized_phase, name="time", index=0)
                returned_object.sort("time")