\section{Introduction}
Among the chemically peculiar (CP) stars, metallic-line A-type stars (Am stars for short) attract the attention of the astronomical researchers, because of their extremely strong or weak element lines, high binary frequency and other characteristics. Almost all Am stars are A or F type stars. Am stars show following common characteristic: slow rotation compared with the normal A or F type stars (less than 120 $km s^{-1}$)\citep{1973ApJ...182..809A}; the spectra will show some remarkably characteristics as follows: weaker Ca II, K line  \citep{1940ApJ....92..256T,1948ApJ...107..107R} but enhanced Sr II line; the scandium and calcium elements are under-abundant while the iron-group elements, Y, Ba, and the rare earth elements are over-abundant \citep{1970PASP...82..781C}. The researchers found that almost all Am stars (more than $90\%$) are the components of binaries \citep{1961ApJS....6...37A,1965ApJS...11..429A,2010ASPC..435..257H}. It is generally accepted that the largest rotational velocity of Am stars is about $120 km s^{-1}$ \citep{1973ApJ...182..809A,2000ApJ...544..933A}. There is a lack of very short periods ($P_{orb} < 1.2 days$) Am type binaries  \citep{1996A&A...313..523B}, and the general explanation is synchronism in such systems would force the primary to rotate faster than 120 $km s^{-1}$. Nevertheless, Am peculiarity may have no significant correlation with rotational velocity \citep{2004IAUS..224..209M,2005A&A...442..563M}. In addition, there are
some normal A0-A3 type slow rotators without reasonable explanation \citep{1973ApJ...182..809A,1978ApJS...37..371W,2007A&A...463..671R}. It is not clear yet whether slow rotation should be one necessary condition for the formation of an Am star, or slow rotation is an individual result of the Am phenomenon. It seems like the Am peculiarity can also depend on evolutionary status (or age)
\citep{2000A&A...354..216B,2004IAUS..224..209M,2005A&A...442..563M}, atmospheric parameters \citep{1998A&A...330..651K,2000A&AS..144..203H}or orbital elements in a binary system \citep{1996A&A...313..523B,1997A&A...326..655B,1998A&AS..128..497I,2004IAUS..224..749F} as well. In classical Am stars, the spectral types inferred from the Ca II K line (Sp(K)) are early than that inferred from the metal lines (Sp(m)) with about five or more spectral subclasses. While stars in which the difference of spectral
type ($\Delta=Sp(K)-Sp(m)$) is less ($1 \leq \Delta < 5$) are often referred to as ‘marginal Am’(or Am:) stars.  Many Am stars show anomalous luminosity effect, i.e. the luminosity criteria in certain spectral regions will indicate that the target should be a giant or even supergiant star, whereas in other regions the luminosity will indicate it as a dwarf or even lower luminosity star.

Some theories were presented to expound the origin of Am phenomenon. The radiative diffusion occurred in a strong magnetic field will likely lead to chemical peculiarities of Ap stars \citep{1970ApJ...160..641M}, while when the magnetic field is absent, the diffusion will caused Am phenomenon \citep{1971A&A....13..263W}. The slow rotation will further assists the diffusion to segregate the elements \citep{1983ApJ...269..239M}. The spin braking may be a efficient process for the Am phenomenon in close binaries, but it
is difficult to explain that phenomenon occurred in single stars. Furthermore, the significantly weak but observable magnetic fields were found rarely in Am stars (e.g. Sirius A \citep{2011A&A...532L..13P}, $\beta$ UMa and $\theta$ Leo \citep{2016A&A...586A..97B}). Taking the radiative acceleration and atomic diffusion into consideration, \citet{2000ApJ...529..338R} developed the stellar evolutionary models. These models can produced resemble abundance anomalies resemble to that of Am stars with much larger value. We should note that the standard stellar evolution theory, the accretion processes effects of interstellar and the circum stellar gases should all be taken into account for the investigate of the atmospheric chemical abundances of Am stars in the binary systems. As proposed, the accretion processes may strongly influence the surface abundances, such as the mass transfer from the evolved companion \citep{1965ApJ...142..423F,1989ApJ...345..998P}. \citet{1965ApJ...142..423F} proposed that the peculiar A and B type stars have evolved into an advanced phase, Such stars have evolved through the giant phase and returned to the vicinity of the main sequence. In other words, such a process has occurred in one component of the close binary, mass transfer has happened from the evolved component to the peculiar star.

Near contact binaries (NCBs) were defined as eclipsing binary with following common characteristics: continuous EB-type light variations, facing surface less than 0.1 orbital radius apart, short periods (less than one day) and one or two component at or near their Roche lobes (RLs) \citep{1994MmSAI..65...95S}. NCBs actually were classified into the following subclasses: semi-detached, marginal contact and marginal detached systems \citep{2006MNRAS.367..423Z}. The mass transfer is frequently happened in NCBs. NCBs have been significant targets to study the transition  between the tidal-locked detached stage and W UMa type overcontact stage, occurring in binary systems.

TYC 6408-989-1 (ASAS J235103-1904.5, NSVS 14636842, GSC 06408-00989) is an EB type Am type eclipsing binary with a period of 0.470796 day \citep{1991A&AS...89..429R,1971AJ.....76..338S,2009A&A...498..961R,2014A&A...564A..69S}.
The spectral type of the primary component is A4m \citep{2014A&A...564A..69S}. \citet{1971AJ.....76..338S} found some early type stars including this target near the south galactic pole. The period of TYC 6408-989-1 is shortest among the well known Am type eclipsing binaries \citep{2009A&A...498..961R,2012A&A...537A.117T,2015AJ....150..154T}.
\citet{2014A&A...564A..69S} presented the light curve of this target obtained by the Super WASP(SWASP) survey without photometric solutions. No detailed investigation on this target has been carry out.

The first $BV$ bands light curves and low resolution spectra of TYC 6408-989-1 were obtained in this study. The detailed light curves analysis and spectral fitting were presented. We collected all well known Am type eclipsing binary with absolute parameters from the literature. The relationships between these parameters are summarized. Furthermore, the evolutionary stage of these systems and the reason of the chemical peculiarities are discussed.

\section{OBSERVATION AND DATA REDUCTION}

\subsection{PHOTOMETRIC OBSERVATIONS}
The first $BV$ bands light curves of Am type eclipsing binary TYC 6408-989-1 were observed with the 70 \emph{cm} Sino-Thai telescope at Lijiang Observing Station of Yunnan Observatories (YNO), Chinese Academy of Sciences (CAS) on 2015 December 9, 11, 13, 24. The PHOT package of IRAF was used to process all the observations. The comparison star and the check star were chosen to determine the differential magnitudes. The information of the comparison and check stars along with the target TYC 6408-989-1 are listed in Table \ref{table1}. In the table, the V band magnitude of variable star and the comparison star Ch1 were taken from \citet{2008AJ....136.1067H} and \citet{2014AJ....148...81M}, respectively. The G band magnitudes of other comparison stars and check stars were acquired from the Gaia DR2 data\citep{2018A&A...616A...1G}.
The observed CCD image of TYC 6408-989-1 was shown in Figure \ref{fig1}, in which the variable star was marked with 'V', and the comparison, check stars were marked as C1, C2, C3 and Ch1, Ch2, Ch3, respectively. The mean flux of comparison stars C1,C2,C3 (regarded as the comparison flux) and the mean flux of the check stars Ch1, Ch2, Ch3 (regarded as the check flux) were used to obtain the final differential magnitudes of TYC 6408-989-1, which show the light curves in Figure 2. In the figure, different colors marked the data observed on different nights. The magnitude differences between the comparison stars and the check stars are show in the bottom of the figure.

\begin{figure}
  \begin{center}
   \includegraphics[width=12cm]{ms2020-0214fig1.eps}
  \end{center}
  \caption{Observed CCD image of TYC 6408-989-1. 'Variable star' was marked with 'V'. C1, C2, C3 marked the comparison stars and  Ch1, Ch2, Ch3 represented the check stars. }
\label{fig1}
\end{figure}

\begin{figure}
  \begin{center}
   \includegraphics[width=12cm]{ms2020-0214fig2.eps}
  \end{center}
  \caption{Light curves of TYC 6408-989-1 on $BV$ bands observed with 70 \emph{cm} Sino-Thai telescope.}
\label{fig2}
\end{figure}

\begin{table}
  \caption{Information of TYC 6408-989-1, the comparison stars and the check stars.}
   \centering
\small
\label{table1}
\begin{tabular}{lclll}
\hline
 Targets & Name &${\alpha}_{2000}$  &  ${\delta}_{2000}$ & Mag    \\
\hline
Variable(V)  & TYC 6408-989-1& ${23}^{h}{51}^{m}{03}^{s}.5689$ & ${-19}^{\circ}{04}^{'}{29}^{''}.1076$ & V=12.064 \\
\hline
\multirow{3}{*}{Comparison(C)}
&C1(UCAC2 24500817) & ${23}^{h}{51}^{m}{06}^{s}.2789$   & ${-19}^{\circ}{11}^{'}{40}^{''}.3554$& V=12.390  \\
&C2(Gaia DR2 2390141539218676096)  & ${23}^{h}{50}^{m}{59}^{s}.6221$ & ${-19}^{\circ}{01}^{'}{24}^{''}.8279$ & G=13.8351 \\
&C3(Gaia DR2 2390132433888003200)&  ${23}^{h}{50}^{m}{38}^{s}.6209$ & ${-19}^{\circ}{09}^{'}{17}^{''}.6696$  & G=13.6923  \\
\hline
\multirow{3}{*}{Check(Ch)}
&Ch1(Gaia DR2 2390135560624200320) & ${23}^{h}{50}^{m}{58}^{s}.2628$ & $ {-19}^{\circ}{02}^{'}{42}^{''}.6827$ & G=14.0021  \\
&Ch2(Gaia DR2 2390130509742663808)& ${23}^{h}{51}^{m}{05}^{s}.3924$ & $ {-19}^{\circ}{10}^{'}{47}^{''}.3473$ & G=15.4127  \\
&Ch3(Gaia DR2 2390130475382925568)& ${23}^{h}{51}^{m}{03}^{s}.3342$ & $ {-19}^{\circ}{11}^{'}{11}^{''}.0596$ & G=14.966  \\
\hline
\end{tabular}
\end{table}

\subsection{SPECTRAL OBSERVATIONS}
The spectrum of TYC 6408-989-1 was observed on 2016 December 23 with the Beijing Faint Object Spectrograph and Camera (BFOSC) mounted to the 2.16 \emph{m} telescope of Xinglong station of National Astronomical Observatories of China (NAOC), Chinese Academy of Sciences (CAS). The low-dispersion spectrometer BFOSC and grating G7 were used during the observations. The slit width and line dispersion of grating G7 are 1.8 arcsec and 95 ${\AA} mm^{-1}$, respectively. The observable wavelength range is $4000-6800 {\AA}$. The observations process and spectra extraction were done using IRAF. The fluxes were normalized and the atmospheric absorption lines were corrected. In such a low resolution, the observed spectra only show the spectral lines of the primary component. The observed spectrum was shown in the upper of Figure\ref{fig3} with a black line. The University of Lyon Spectroscopic analysis Software (ULySS) \citep{2009A&A...501.1269K} was employed to acquire the atmospheric parameters through the full spectra fitting with the model spectra generated by an interpolator with the ELODIE library \citep{2001A&A...369.1048P}.
The fit spectrum was shown in the upper of Figure \ref{fig3} with a red line. The obtained atmospheric parameters are as follows: $T_{eff}\!=\!6990\!\pm\!117 K$, $\log g\!=\!4.25\!\pm\!0.26 cm/s^2$, $[Fe/H]\!=\!-0.55\!\pm\!0.03 dex$.
The primary star contributes the most light to the total system, the above atmospheric parameters were commonly applied to show the atmospheric characteristics of the primary star.

The spectroscopic data of TYC 6408-989-1 are available in European Southern Observatory (ESO) archives. The data were observed with the spectrometer Faint Object Spectrograph and Camera (EFOSC2) \citep{1984Msngr..38....9B} mounted at the New Technology Telescope (NTT) at La Silla Paranal observatory. During the observation, a 600 $g mm^{−1}$ grating (i.e. grism 14) was used. Grism 14 has a resolving power of 7.54 ${\AA}$ and a spectral range of $3095-5058$ ${\AA}$. The slit width and the dispersion of grism 14 are 1.0 $arcsec$ and 0.93 ${\AA} pixel^{-1}$. The observations process and spectra extraction were done by us using IRAF. An overview of the ESO spectrum was shown in the bottom of Figure \ref{fig3}. The ESO spectra were used to identify the spectra type through the MKCLASS program \citep{2014AJ....147...80G}, which was designed to classify stellar spectra on the MK Spectral Classification system in a way similar to humans—by direct comparison with the MK classification standards. The standard library $libr18$ was used in our classification progress, which consists of 1.8 ${\AA}$ resolution rectified spectra with a spectral range from $3800-4600 {\AA}$. The spectra in standard library $libr18$ were obtained with a 1200 $g mm^{−1}$ grating on the GM spectrograph of the Dark Sky Observatory (Appalachian State University).

A comparison of the blue-violet spectrum ($3800-4600$ ${\AA}$) of TYC 6408-989-1 (kA3hF1mA5 IV-V) and part of the spectrum of two MK standards Bet Leo (A3 V) and HD 23194 (A5 V) are shown in Figure \ref{fig4}. It is clearly that the Ca II K line of TYC 6408-989-1 is slightly weaker in strength than that of the A3V and A5V MK standards. The spectral type derived from the MKCLASS is kA3hF1mA5, and the index of spectral type difference $\Delta$ is 2. Therefore, we conclude that TYC 6408-989-1 is a marginal Am (Am:) star with a spectral type of kA3hF1mA5 IV-V.

\begin{figure}[htp!]
  \begin{center}
   \includegraphics[width=12cm]{ms2020-0214fig3.eps}
  \end{center}
  \caption{Observed spectrum of TYC 6408-989-1. The black and red line in the upper panel represent the observed and fitted spectra, respectively. The blue line in the bottom panel shows an overview of the ESO spectrum.}
\label{fig3}
\end{figure}

\begin{figure}[htp!]
  \begin{center}
   \includegraphics[width=12cm]{ms2020-0214fig4.eps}
  \end{center}
  \caption{A comparison of the blue-violet spectrum ($3800-4600$ ${\AA}$) of TYC 6408-989-1 (kA3hF0mA5 IV-V) and part of the spectrum of two MK standards Bet Leo (A3 V) and HD 23194 (A5 V). Note the difference in Ca II K line strength between TYC 6408-989-1 and the two MK standards.}
\label{fig4}
\end{figure}

\section{LIGHT CURVE ANALYSES}
The first $BV$ bands light curves of TYC 6408-989-1 were phased with following linear ephemeris:
\begin{eqnarray}
Min \emph{I}& = &HJD\,2457366.05772( .00037)+0.470796^{d}\times E
\end{eqnarray}
In the equation, the period was obtained from \citep{2014A&A...564A..69S}.
The Wilson-Devinney (W-D) program
\citep{1971ApJ...166..605W,1990ApJ...356..613W,2012AJ....144...73W} was used for the light curves analysis. The temperature of the primary component was estimated as $T_1=6990K$ through the spectral fit. The light curves with all modes of WD have been tried. The finally convergent solutions were achieved with mode 3 (over-contact mode). The gravity-darkening coefficients of both components are $g_1=g_2=0.32$ and bolometric albedos are $A_1=A_2=0.5$ \citep{1967ZA.....65...89L,1969AcA....19..245R}. The bandpass limb-darkening coefficients \citep{1993AJ....106.2096V} and the logarithmic bolometric coefficients were applied during the WD modelling.
There are some adjustable parameters in mode 3: the monochromatic luminosity of the primary star, $L_{1B}$, $L_{1V}$; the orbital inclination, i; the mean temperature of the secondary star, $T_2$; and the dimensionless potential of the primary star,$\Omega_1=\Omega_2$. The third light $L_3$ was an adjustable parameter, but no converged result was obtained. The value of the third light is negative and keeps decreasing, which means that there is no third light can be detected through the light curve analysis.

The mass ratio search (q-search) was a common technique to acquire the mass ratio.
The photometric solutions were obtained with some assumed mass ratio values from 0.01 to 1 with the differential correction program. The step of the search is 0.01. The relation between the sums of weighted square deviations ($\sum{(O-C)_i}^2$) and mass ratio ($q$) are shown in Figure \ref{fig5}. The minimal value achieved at $q=0.27$, which means that the solution at $q=0.27$ is the best fit. Then the mass ratio $q$ was adjustable. The final converged photometric solutions are listed in Table \ref{table2}. The uncertainty in the mass ratio in the table was given by the WD program with the standard method. In addition, the q-search curve bottom can be used to confirm the uncertainty of the mass ratio. It can be seen from Figure \ref{fig5}, the q-search bottom is clearly narrow around 0.270(3)(visually, from 0.25 to 0.28) and that the uncertainty is approximately 0.04. The theoretical light curves of TYC 6408-989-1 were shown in Figure \ref{fig6} with black lines. The standard deviations of the residuals on $BV$ bands both are about 0.012 magnitudes. The geometric structure in 3D view is plotted in Figure \ref{fig7}.

\begin{figure}[htp!]
 \begin{center}
  \includegraphics[width=12cm]{ms2020-0214fig5.eps}
  \end{center}
  \caption{$\sum - q$ curve of TYC 6408-989-1.}
  \label{fig5}
\end{figure}

\begin{figure}[htp!]
 \begin{center}
  \includegraphics[width=14cm]{ms2020-0214fig6.eps}
  \end{center}
  \caption{The theoretical light curves of TYC 6408-989-1.}
  \label{fig6}
\end{figure}

\begin{figure}[htp!]
 \begin{center}
     \includegraphics[width=14cm]{ms2020-0214fig7.eps}
  \end{center}
   \caption{Geometrical structure of TYC 6408-989-1 in 3D view.}
   \label{fig7}
\end{figure}

\begin{table}
\caption{Photometric solutions of TYC6408-989-1.}
\normalsize
\begin{center}
\label{table2}
\begin{tabular}{lc|lc}
\hline
Parameters       & Values& Parameters &Values  \\
\hline
Mode                             & Mode 3       &   ${L_{1}/L_{total}}_V$   & 0.93683(16)      \\
$g_{1}=g_{2}$                    & 0.32(fixed)  &   $f$                     & 0.003(29)  \\
$A_{1}=A_{2}$                    & 0.5(fixed)   &   $r_{1}(pole)$           & 0.46460(71)  \\
q ($M_2/M_1$ )                   & 0.270(3)   &   $r_{1}(side)$           & 0.50088(94)  \\
$i(^{\circ})$                    &  69.31(26)   &   $r_{1}(back)$           & 0.5255(11)  \\
$T_{1}(K)   $                    & 6990         &   $r_{2}(pole)$          & 0.2541(32)  \\
$T_{2}(K)$                       & 5089(50)     &   $r_{2}(side)$           &  0.2646(38)  \\
$\Omega_{1}$                     & 2.3985(49)   &   $r_{2}(back)$          &  0.2975(68)   \\
$\Omega_{2}$                     & 2.3985       &   $R_2/R_1$              & 0.5480(55)  \\
${L_{1}/L_{total}}_B$             & 0.959822(89)&   $\sum(O-C)^2$         & 0.0054          \\
\hline
\end{tabular}
\end{center}
\end{table}

\section{DISCUSSION AND CONCLUSION}
The $BV$ light curves of TYC 6408-989-1 were obtained firstly. The light curves show $\beta$ Lyrae characteristics. The atmospheric parameters of the primary star are estimated through the low resolution spectral fit as follows:  $T_{eff}\!=\!6990\!\pm\!117 K$, $\log g\!=\!4.25\!\pm\!0.26 cm/s^2$, $[Fe/H]\!=\!-0.55\!\pm\!0.03 dex$. The light curves analysis were done using WD program. The contact degree factor of TYC 6408-989-1 is 0.003(29), which is extremely low and the error is quite large. That means that TYC 6408-989-1 should be a marginal contact binary, in which both components has nearly filled their critical RLs, just like UX Eri \citep{2007AJ....134.1769Q}, AS Ser \citep{2008AJ....136..337Z}, DD Com \citep{2010AJ....140..215Z} and DI Hya \citep{2017PASP..129c4201L}.
TYC 6408-989-1 is one of the three well known Am type marginal contact binaries, the others are V1073 Cyg \citep{2018RAA....18...20T} and V2782 Ori \citep{2019PASP..131h4203T}. The orbital inclination of this target is $i=69.31(26)^\circ$. The mass ratio was estimated as about $q=0.270(3)$ through the mass ratio search progress. TYC 6408-989-1 was concluded as a marginal Am (Am:) star with a spectral type of kA3hF1mA5 IV-V. The difference between Sp(K) and Sp(M) is two subclasses ( i.e. $\Delta=2$). TYC 6408-989-1 consists of an F1 type primary star with $M_1=1.56(12)M_\odot$ \citep{2000PhT....53j..77C}. The spectra of the secondary star should be K0 type, which can be inferred from the temperature. Based on the photometric solutions, the absolute parameters of TYC 6408-989-1 were calculated as follows:
$M_2=0.42(10)M_\odot$, $R_1=1.581(3)R_\odot$, $R_2=0.874(2)R_\odot$, $L_1=5.360(378)L_\odot$, $L_2=0.459(3)L_\odot$. TYC 6408-989-1 should be a B sub-type contact binaries with temperature difference between the two stars (which is about 1860K) larger than 1000K \citep{1979ApJ...231..502L,2004A&A...426.1001C}, like V2787 Ori\citep{2019PASP..131h4203T}.

Using the Kepler's third law and the photometric solutions with mode 3, the mean densities of the primary and secondary star were calculated as $\rho_1=0.388\rho_\odot$ and $\rho_2=0.636\rho_\odot$, respectively. The corresponding logarithmic values are $\log \rho_1/\rho_\odot=-0.411$ and $\log \rho_2/\rho_\odot=-1.97$. The mean densities of both components (especially the secondary star) are significantly lower than main sequence stars with the same spectral type \citep{2000PhT....53j..77C}, which indicates that the components of TYC 6408-989-1 may have evolved away from the Zero Age Main Sequence (ZAMS) line and the secondary star may have a higher degree of evolution. The luminosity class of the target from IV to V indicates that the primary should be an evolved subgiant.

The components in very short period NCB systems like TYC 6408-989-1 are likely rotating synchronously. For a circular orbit, following equation can be used to obtain the synchronism rotational velocity:
$v = 50.6(R/R_\odot)/(P/d) km s^{-1}$\citep{2007MNRAS.380.1064C}, in which $R$ is the radius of the considered component, $P$ is the orbital period, $v$ is its equatorial rotational velocity. The synchronous rotation velocity $v \sin i$ of the primary  should be about 160 $km s^{-1}$. TYC 6408-989-1 is a very especial Am type eclipsing binary with very short period (less than one day), marginal Am peculiarity and quit large rotational velocity, which make it became a challenge to the cut-off of rotation velocities and periods of Am stars, like V1073 Cyg ($v \sin i = 150 km s^{-1}$)\citep{1996A&A...313..523B}.

According to the thermal relaxation oscillation theory \citep{1976ApJ...205..208L,1976ApJ...205..217F,1977MNRAS.179..359R,1979ApJ...231..502L}, W UMa type systems must undergo oscillations around the state of marginal contact. TYC 6408-989-1 should be a particular target like other marginal contact binaries lying in the TRO stage \citep{2002ApJ...568.1004Q,2009AJ....137.3574Z,2012AJ....144...37Z}
Such targets are very important for investigation of the formation between detached and contact phase of binaries, which are quite rarely in observation. The light curves of these systems show $\beta$ Lyrae-type light-variation characteristics. In the evolutionary process, the systems may be in a semidetached phase with the more massive star filling the RL, and then evolve into a marginal-contact phase with poor thermal contact. TYC 6408-989-1 should be lying on the poor thermal contact stage. It is another promising candidate locating in such rare evolutionary stage, it will evolve to become a shallow contact binary with true thermal contact(e.g.,\citet{2012AJ....144..178L,2013AJ....146...38Q, 2014ApJS..212....4Q} and \citet{2019PASP..131a4202L}).

Some eclipsing Am binaries with known parameters were collected and tabulated in Table \ref{table3}. The temperature-luminosity diagram of Am type EBs are plotted in Figure \ref{fig8}. In the figure, the ZAMS and terminal-age main sequence (TAMS) lines are obtained from \cite{1992A&AS...96..269S} at Z=0.020. From the figure, we can see that both components of Am type EBs are located in the main sequence. According to \citet{1965ApJ...142..423F}, one component in the Am type binary may has evolved through the giant phase, it may have striped most of the hydrogen shell and the stellar hell were exposed leading to higher temperature, then it will returned to main sequence. The statistical result support the above theory well. The accretion processes (such as mass transfer) may be very important for the formation of Am peculiarity in binaries. The mass transfer from the evolved component of Am type eclipsing binaries will recombine the chemical elements in the atmosphere of the companion stars and may result in chemical anomalies composition, such as the marginal contact binary V1073 Cyg \citep{2018RAA....18...20T}, V2787 Ori \citep{2019PASP..131h4203T} and TYC 6408-989-1, which are lying on the thermal relaxation oscillation state \citep{1976ApJ...205..208L,1977MNRAS.179..359R,1979ApJ...231..502L}. Mass transfer between the components may paly an important role in the formation and evolution of Am stars in binaries. More observation evidence are need to support this discussion. The researches on individual Am type binaries including the binaries with very short periods (less than one day) will be very important and indispensable for exploring the formation and evolution of Am star.

\begin{figure}[htp!]
 \begin{center}
  \includegraphics[width=14cm]{ms2020-0214fig8.eps}
  \end{center}
  \caption{The temperature-luminosity diagram of Am type EBs.}
  \label{fig8}
\end{figure}

\begin{landscape}
\begin{table}[]
\centering
\caption{Parameter of some known Am type eclipsing binaries.}
\label{table3}
\begin{tabular}{llllllllllll}
\hline
Star           & Renson & Period     & q$(M_2/M_1)$ & $M_1$       & $M_2$       & $T_1$     & $T_2$     & $L_1$       & $L_2$       & Type        & Ref        \\
-              &        & (days)     &              & $(M_\odot)$ & $(M_\odot)$ & (K)       & (K)       & $(L_\odot)$ & $(L_\odot)$ &             &            \\
\hline
V1073 Cyg      & 56830  & 0.7858506  & 0.303        & 1.810(4)    & 0.549(1)    & 7300      & 6609(18)  & 16.51       & 3.74        & W UMa       & (1)        \\
V2787 Ori      & 11387  & 0.810979   & 0.120(2)     & 1.44        & 0.17(1)     & 6993(82)  & 5386(29)  & 12.93       & 0.71        & W UMa       & (2)        \\
TYC 6408-989-1 & 61280  & 0.470796   & 0.230(3)     & 1.56(12)      & 0.42(10)        & 6990      & 5089(50) )  & 5.36(38)        & 0.459(3)        & W UMa       & This study \\
TX Her         & 44140  & 2.05984    & 0.895        & 1.61        & 1.44        & 8180      & 7536(20   & 9.44        & 6.56        & Algol       & (3)        \\
YZ Cas         & 1140   & 4.4672235  & 0.585(2)     & 2.263(12)   & 1.325(7)    & 9520(120) & 6880(240) & 46.989      & 3.556       & Algol       & (4)(5)     \\
WW Aur         & 12320  & 2.52501941 & 0.9235(27)   & 1.964(7)    & 1.814(7)    & 7960(420) & 6411(410) & 13.458      & 10.544      & Algol       & (6)        \\
GZ CMa         & 15390  & 4.80085    & 0.909        & 2.200(25)   & 2.000(25)   & 8810      & 8694      & 33.496      & 21.627      & Algol       & (7)    \\
V624 Her       & 45460  & 3.894977   & 0.8326       & 2.270(14)   & 1.870(13)   & 8147      & 7943      & 36.308      & 16.982      & Algol       & (8)        \\
DV Boo         & 35946  & 3.782624   & 0.75         & 1.64(2)     & 1.23(2)     & 7370(80)  & 6410(74)  & 10.715      & 3.090       & Algol       & (9)        \\
V885 Cyg       & 50855  & 1.69478781 & 1.114        & 2.005(29)   & 2.234(26)   & 8375(150) & 8150(150) & 23.988      & 45.709      & Algol       & (10)       \\
EI Cep         & 57180  & 8.4393522  & 0.9483       & 1.7716(66)  & 1.6801(62)  & 6750(100) & 6950(100) & 15.596      & 11.350      & Algol       & (11)       \\
EE Peg         & 57410  & 2.628      & 0.6186       & 2.15(2)     & 1.33(1)     & 8709      & 6456      & --          & --          & Algol       & (12)       \\
SW CMa         & --     & 10.091988  & 0.9397       & 2.239(14)   & 2.104(18)   & 8200(150) & 8100(150) & 36.813      & 19.099      & Algol       & (13)       \\
HW CMa         & --     & 21.1178329 & 0.9665(43)   & 1.721(11)   & 1.781(12)   & 7560(150) & 7700(150) & 7.907       & 8.710       & Algol       & (13)       \\
AN And         & 60340  & 3.219566   & 0.53         & 2.48        & 1.32        & 8200      & 6330      & 49.2        & 17.5        & $\beta$ Lyr & (14)       \\
V501 Mon       & --     & 7.02       & 0.8865       & 1.6455(43)  & 1.4588(25)  & 7510(100) & 7000(90)  & 11.940      & 5.534       & Algol       & (15)       \\
\hline
\end{tabular}

\footnotesize{
\centering
(1)\citet{2018RAA....18...20T} (2)\citet{2019PASP..131h4203T}  (3)\citet{2019RAA....19...94Z}   (4)\citet{2014MNRAS.438..590P} (5)\citet{2000astro.ph.11388L}
(6)\citet{2005MNRAS.363..529S} (7)\citet{1985AJ.....90.1324P}  (8)\citet{1984AJ.....89.1057P}\\(9)\citet{2004MNRAS.352..708C} (10)\citet{2004AJ....128.1324L}
(11)\citet{2000AJ....119.1942T} (12)\citet{1984ApJ...281..268L}(13)\citet{2012A&A...537A.117T}
(14)\citet{1978JRASC..72..263T} (15)\citet{2015AJ....150..154T}
}
\end{table}
\end{landscape}

\begin{acknowledgements}
We greatly appreciate and feel great thankfulness for the hard work in the fight against COVID-19 to the Chinese people and all the people worldwide. We would like to thank the editor and the referee very much for the very valuable and useful comments that helped improve this paper. This work is partly supported by the National Natural Science Foundation of China (Nos. 11922306, U1831109), and the CAS Interdisciplinary Innovation Team. We acknowledge the support of the staff of the Xinglong 2.16 \emph{m} telescope. This work was partially supported by the Open Project Program of the Key Laboratory of Optical Astronomy, National Astronomical Observatories, Chinese Academy of Sciences. We acknowledge the kindly help of Pro. Gray Richard O.

CCD photometric observations of the system were obtained with the 70 \emph{cm} telescope at Lijiang station, Yunnan Observatories. The spectral observations were obtained with the 2.16 \emph{m} telescope at Xinglong station of the National Astronomical Observatory. The spectra classification were based on observations made with ESO Telescopes at the La Silla Paranal Observatory under programme 083.D-0472(A).
\end{acknowledgements}

\newpage


\section{#1}\setcounter{equation}{0}}
\newcommand{\subsectiono}[1]{\subsection{#1}\setcounter{equation}{0}}
\renewcommand\O{{\mathcal{O}}}
\newcommand{\be}[1]{ \begin{equation}\label{#1} }
\newcommand{\end{equation}}{\end{equation}}
\newcommand{\bea}[1]{\begin{eqnarray}\label{#1} }
\newcommand{\end{eqnarray}}{\end{eqnarray}}
\newcommand{\partial}{\partial}
\newcommand{\nonumber}{\nonumber}
\newcommand{\Delta}{\Delta}
\newcommand{\tilde{\Delta}}{\tilde{\Delta}}
\newcommand{\tilde{\xi}}{\tilde{\xi}}
\newcommand{\omega}{\omega}
\newcommand{\bar{\omega}}{\bar{\omega}}
\newcommand{\refb}[1]{(\ref{#1})}
\newcommand{\<}{\langle}
\renewcommand{\>}{\rangle}
\renewcommand{\(}{\left(}
\renewcommand{\)}{\right)}
\newcommand{\mbox{SL}(2,\mathbb{R})}{\mbox{SL}(2,\mathbb{R})}
\renewcommand{\a}{\alpha}
\newcommand{\tilde{\a}}{\tilde{\a}}
\renewcommand{\b}{\beta}
\renewcommand{\t}{\tau}
\newcommand{\sigma}{\sigma}
\newcommand{|0\rangle_\a}{|0\rangle_\a}
\newcommand{|0\rangle_c}{|0\rangle_c}

\begin{document}

\title{Rindler on the Worldsheet}

\author{Arjun Bagchi}
\email{abagchi@iitk.ac.in}
\affiliation{Indian Institute of Technology, Kanpur 208016, India.}

\author{Aritra Banerjee}
\email{aritra.banerjee@apctp.org}
\affiliation{Asia Pacific Center for Theoretical Physics, Postech, Pohang 37673, Korea.}

\author{Shankhadeep Chakrabortty.}
\email{s.chakrabortty@iitrpr.ac.in}
\affiliation{Indian Institute of Technology Ropar, Rupnagar, Punjab 140001, India.}

\begin{abstract}
We construct the tensionless limit of bosonic string theory in terms of a family of worldsheets with increasing acceleration and show that the null string emerges in the limit of infinite acceleration when the Rindler horizon is hit. We discover a novel phenomenon we call null string complementarity, which gives two distinct observer dependent pictures of the emergence of open string physics from closed strings in the tensionless limit. The closed string vacuum as observed by the inertial worldsheet turns into a D-instanton in the tensionless limit, while in the complementary picture from the accelerated worldsheet one sees the emergence of a D-25 brane. We finally discuss approaching the Rindler horizon through time evolution at constant acceleration and also show how an open string picture arises very naturally.
\end{abstract}

\maketitle

\noindent {\em{\underline{Introduction}}}.
The study of physics in accelerated frames of reference is an intriguing and fruitful venture. Accelerated observers in Minkowski spacetimes view the Rindler metric and experience a horizon. The entire Minkowski diamond not accessible to them anymore and in this frame of  reference, physics is thermal. Due to the unavailability of information from beyond the Rindler horizon, the density matrix of this accelerated observer is related to that of an inertial Minkowski observer by partially tracing over the inaccessible degrees of freedom. All of this is of course well understood and is very useful for black hole physics since the near-horizon limit of a black hole typically yields a Rindler spacetime. Our discussions in this paper are rather unique. We want to understand aspects of Rindler physics {\em on the worldsheet of a closed string}. Our motivation for doing so is also rather novel, as we elaborate below.

String theory is currently the most promising of avenues for formulating a theory of quantum gravity. One of its primary endearing features is the paucity of tuneable parameters. The free theory has only one, the length of the fundamental string ($\ell_s$). $\ell_s\to0$ reduces the string to a point particle and the theory to the well understood Einstein's theory of general relativity. In this work, we are interested by the other extreme limit, where $\ell_s\to\infty$ \cite{Schild:1976vq}. This bizarre limit corresponds to the ultra-stringy regime of string theory that is very different to Einstein's theory. The quantum version of this theory would be ``very stringy" quantum gravity. This limit also takes the tension of the fundamental string to zero. Our objective in this paper is to {\em formulate the decreasing string tension in terms of accelerated string worldsheets}.

The tensionless string is a null string with a degenerate worldsheet metric. We will show in this paper that the worldsheet analogue of hitting the Rindler horizon leads to the formation of the null string. Rindler observers can approach the horizon in two distinct ways: time evolution at a fixed acceleration or evolution in acceleration at a constant time. The limit from the tensile to the tensionless string can also be formulated in terms of increasingly accelerated worldsheets or time-evolution on a constant acceleration worldsheet. We shall follow both routes with interesting consequences, the most intriguing amongst which is what we call {\em null string complementarity}. Depending on whether the observer sits on an inertial worldsheet and observes an accelerated one, or vice-versa, she sees the emergence of different complementary boundary states as the closed string becomes tensionless. The inertial worldsheet sees the accelerated closed string vacuum evolve into a spacetime point, a D-instanton, while the infinitely accelerated observer sees the inertial vacuum grow into a spacefilling D-25 brane. No one observer has access to both pictures.

\smallskip

\noindent {\em{\underline{Rindler Physics}}}. Accelerated observers in Minkowski spacetimes moving on trajectories $x^2 - t^2 = \kappa^{-2}$ ($\kappa$ is proper acceleration) describe Rindler space with metric
\be{}
ds_{R}^2 = e^{2a\xi}(-d\eta^2+d\xi^2).
\end{equation}
Minkowski and Rindler spacetimes are linked by:
\be{Rmap}
t = \frac{1}{a}e^{a\xi}\sinh a\eta,~~ x= \frac{1}{a}e^{a\xi}\cosh a\eta.
\end{equation}
where $\kappa =a e^{-a\xi}$, $a$ is the redefined acceleration. \refb{Rmap} is for the right Rindler wedge only ($R: |t|< x, x>0$). There is an equivalent left wedge ($L: |t|< x, x<0$), for which Eqs \refb{Rmap} pick up negative signs.
\begin{figure}[h]
\centering
  \includegraphics[width=50mm]{rindlerspace.pdf}
  \caption{Equal time slices in Rindler spacetimes}
  \label{fig1}
\end{figure}

We now consider a massless scalar field theory in Rindler spacetime \cite{Birrell:1982ix}. Since flat and Rindler backgrounds are conformally related, the equations of motion (EOM) for the field are same:
\be{}
\Box_{t,x}\phi =0=\Box_{\eta,\xi}\phi.
\end{equation}
For the Minkowski solution, now defined on a cylinder $(\sigma,\t)$, the wave equation is solved by
\bea{tensilemode}
  \phi (\sigma, \t) &=&\phi_0+\sqrt{2\a'}\a_0\t \\ && +i\sqrt{\frac{\a'}{2}}\sum_n\Big[\frac{\a_n}{n} e^{-in(\t+\sigma)}+\frac{\tilde{\a}_n}{n} e^{-in(\t-\sigma)}\Big] \nonumber
   \end{eqnarray}
where oscillators satisfy $[\a_n,\a_m]=n\delta_{n+m}$ and annihilate the Minkowski vacuum $|0\>_M$. To match with usual scalar field modes, we write \refb{tensilemode}:
\bea{}
\phi (\sigma, \t) &=& \phi_0+\sqrt{2 \alpha'} \alpha_0\tau \\ && + \sqrt{2 \pi \alpha'}\sum_{n>0}[\a_n u_n+\a_{-n} u_n^*+\tilde{\a}_n \tilde{u}_n+\tilde{\a}_{-n} \tilde{u}_n^*] \nonumber
\end{eqnarray}
where $u_n = \frac{i e^{-in(\tau+\sigma)}}{\sqrt{4\pi}n},~~ \tilde{u}_n=\frac{i e^{-in(\tau-\sigma)}}{\sqrt{4\pi}n}$.
Similarly, we write the Rindler mode expansion as
\bea{}
\phi(\xi, \eta) &=& \phi_0+\sqrt{ 2 \alpha'}\beta_0\xi  \\&&+ \sqrt{2 \pi\alpha'}\sum_{n>0}[\b_n U_n+\b_{-n} U_n^*+\tilde{\b}_n \tilde{U}_n+\tilde{\b}_{-n} \tilde{U}_n^*] \nonumber
\end{eqnarray}
where the mode functions are now defined as
\be{}
U_n = \frac{i e^{-in(\xi+\eta)}}{\sqrt{4\pi}n},~~ \tilde{U}_n=\frac{i e^{-in(\xi-\eta)}}{\sqrt{4\pi}n}.
\end{equation}
The oscillators $(\b,\tilde\b)$ now act on a new vacuum $|0\rangle_R$. Importantly, $U_n$ is only defined in the L wedge (hence called $U_n^{(L)}$) and $\tilde{U}_n$ only in the R wedge ($U_n^{(R)}$), unlike flat space. For continuing between wedges one needs to define smearing functions. The combinations analytically continued in both wedges take the form \cite{Unruh:1976db}:
\be{}
U_n^{(R)} - e^{-\frac{\pi n}{a}}U_{-n}^{(L)*}, \quad U_{-n}^{(R)*} - e^{\frac{\pi n}{a}}U_{n}^{(L)}.
\end{equation}
Using these combinations of modes lead us to Bogoliubov transformations between the two sets of oscillators in Rindler and Minkowski space with a form \cite{Unruh:1976db}:
\begin{subequations}\label{R2M}
\bea{}
\b_{n} &=& \frac{e^{\pi n/2a}}{\sqrt{2\sinh\frac{\pi n}{a}}}~\a_{n}-\frac{e^{-\pi n /2a}}{\sqrt{2\sinh\frac{\pi n}{a}}}~\tilde{\a}_{-n},\\
\tilde{\b}_{n}&=& -\frac{e^{-\pi n/2a}}{\sqrt{2\sinh\frac{\pi n}{a}}}~\a_{-n}+\frac{e^{\pi n/2a}}{\sqrt{2\sinh\frac{\pi n}{a}}}~\tilde{\a}_{n}.
\end{eqnarray}
\end{subequations}
Eq. \refb{R2M} forms the backbone of our analysis in this paper.

\smallskip

\noindent {\em{\underline{Intrinsic look at Tensionless strings}}}.
The starting point of our recapitulation is the action \cite{Isberg:1993av}
\be{LST}
S_{\text{ILST}} = \int d^2 \xi \, V^\a V^\b \partial_\a X^\mu \partial_\b X^\nu \eta_{\mu \nu}
\end{equation}
\refb{LST} can be obtained from the Polyakov action for the bosonic string as tension $T \to0$ \cite{Isberg:1993av}, where vector densities $V^\a$ replace the degenerate worldsheet metric. Like in tensile string theory, \refb{LST} enjoys worldsheet diffeomorphism symmetry and needs to be gauge fixed. In $V^\alpha = (1, 0)$ gauge, the residual symmetry is
\bea{bms}
&& [L_n, L_m] = (n-m) L_{n+m} + c_L\delta_{n+m,0} (n^3-n), \cr
&& [L_n, M_m] = (n-m)M_{n+m} + c_M\delta_{n+m,0} (n^3-n), \cr
&& [M_n, M_m]=0.
\end{eqnarray}
This is the BMS$_3$ algebra (here with $c_L=c_M=0$), which also arises as the asymptotic symmetries of 3d flat spacetimes at its null boundary \cite{Barnich:2006av}, and has been used to construct a notion of Minkowskian holography following \cite{Bagchi:2010eg}. For tensionless strings BMS$_3$ replaces the two copies of the Virasoro algebra that dictate the construction of tensile strings \cite{Bagchi:2013bga}.

In $V^\alpha = (1, 0)$ gauge, the EOM of $V^\a$ give constraints while the EOM for $X$ takes a simple form:
\be{xeom}
\ddot{X}^\mu=0; \quad {\text{Constraints:}} \ \dot{X}\cdot X'=0, \, \dot{X}^2=0.
\end{equation}
With closed string boundary conditions $X^\mu(\tau,\sigma)=X^\mu(\tau,\sigma+2\pi)$, EOM can solved by \cite{Bagchi:2015nca}
\be{mode}
X^{\mu}(\sigma,\tau)=x^{\mu}+\sqrt{\frac{c'}{2}}B^{\mu}_0\tau+i\sqrt{\frac{c'}{2}}\sum_{n\neq0}\frac{1}{n} \left(A^{\mu}_n-in\tau B^{\mu}_n \right)e^{-in\sigma}.
\end{equation}
In the above, $c'$ is a length-scale introduced for dimensional consistency. This expansion also leads to the constraints in the form
\be{lmab}
L_n= \frac{1}{2} \sum_{m} A_{- m}\cdot B_{m+n}, \ M_n= \frac{1}{2} \sum_{m} B_{-m}\cdot B_{m+n}.
\end{equation}
$A,B$ are not the usual harmonic oscillator modes:
\bea{AB}
[A_m,A_n]= [B_m,B_n]=0; [A_m,B_n]= 2 m\delta_{m+n}.
\end{eqnarray}
The algebra of constraints leads to the BMS$_3$ algebra as before. We transform $(A,B)$ into a harmonic oscillator basis:
\be{CC}
2C^{\mu}_n = ({A}^{\mu}_n+B^{\mu}_{n}), \quad 2\tilde{C}^{\mu}_n = (-{A}^{\mu}_{-n}+B^{\mu}_{-n}).
\end{equation}
$C, \tilde{C}$ now have canonical commutation relations analogous to tensile $\a$ oscillators. Mode expansion in terms of these $C$ modes reads \cite{Bagchi:2020fpr}
\bea{cexpansion}
X^\mu(\sigma,\t)&=&x^\mu+2\sqrt{\frac{c'}{2}}C^\mu_0\t +i\sqrt{\frac{c'}{2}} \\ && {\hspace{-1cm}} \times \sum_{n\neq 0} \frac{1}{n}\left[(C^\mu_n-\tilde{C}^\mu_{-n})-in\t (C^\mu_n+\tilde{C}^\mu_{-n})\right]e^{-in\sigma} \nonumber
\end{eqnarray}
with zero modes $C^\mu_0=\tilde{C}^\mu_0=\sqrt{\frac{c'}{2}}k^\mu.$

\smallskip

\noindent {\em{\underline{Tensionless strings as a Carrollian limit}}}.
In the discussion above, the string tension was put exactly to zero. Now we describe a limiting procedure on the worldsheet coordinates that takes the tension to zero \cite{Bagchi:2013bga,Bagchi:2015nca}:
\be{URlim}
\sigma \to \sigma, \ \t \to \epsilon \t, \ \a'\to c'/\epsilon, \ \epsilon \to 0.
\end{equation}
This sends the worldsheet speed of light to zero and is called an Ultra-Relativistic (UR) or a Carrollian limit. In this limit, the worldsheet becomes a 2d Carrollian manifold \cite{Duval:2014uoa, Duval:2014lpa}, with a degenerate metric which is the defining feature of a tensionless or a null string. The worldsheet symmetry generators contract
\be{vir2bms}
L_n= \L_n - \bar{{\mathcal{L}}}_{-n}, \ M_n = \epsilon(\L_n + \bar{{\mathcal{L}}}_{-n}),
\end{equation}
(here $\L_n, \bar{{\mathcal{L}}}_n$  generate the tensile Virasoro algebra) and close to form BMS$_3$. Comparing tensile modes (analogue of \refb{tensilemode}) and the tensionless expansions \refb{mode} we get
\be{ABa}
A_n^{\mu} = \frac{1}{\sqrt{\epsilon}} \left( \a_n^\mu - \tilde{\a}_{-n}^\mu \right), \quad B_n^{\mu} = {\sqrt{\epsilon}} \left( \a_n^\mu + \tilde{\a}_{-n}^\mu \right).
\end{equation}
All classical physics of the tensionless string can be reproduced by following this UR limit. The Carrollian limit is usual perceived as a limit on velocities and hence an infinite boost. The tensionless limit is thus an infinite boost on the Riemannian worldsheet of a tensile string that turns it into a degenerate Carrollian worldsheet.

Interestingly, switching to the language of  $C$ oscillators, we find the emergence of a worldsheet Bogoliubov transformation:
\bea{c1}
C^\mu_n =\frac{1}{2} \left( \sqrt{\epsilon} + \frac{1}{\sqrt{\epsilon}} \right) \a^\mu_n+\frac{1}{2} \left( \sqrt{\epsilon} - \frac{1}{\sqrt{\epsilon}} \right) \tilde{\a}^\mu_{-n}, \cr
\tilde{C}^\mu_n =\frac{1}{2} \left( \sqrt{\epsilon} - \frac{1}{\sqrt{\epsilon}} \right) \a^\mu_{-n}+\frac{1}{2} \left( \sqrt{\epsilon} + \frac{1}{\sqrt{\epsilon}} \right) \tilde{\a}^\mu_{n},
\end{eqnarray}
Although the $C$ oscillators have been defined near $\epsilon\to0$, curiously at $\epsilon=1$ they reduce to the $\a$ oscillators. This encourages us to define a flow valid throughout the parameter space $\epsilon \in [0,1]$. For any evolving oscillator $C(\epsilon)$ interpolating between $\a$ for $\epsilon=1$ and \refb{c1} near $\epsilon\to0$, the vacua defined by the flow $|0(\epsilon)\>$ changes continuously with $\epsilon$:
\be{oe}
|0(\epsilon)\>: \, C_n(\epsilon) |0(\epsilon)\> = \tilde{C}_n(\epsilon) |0(\epsilon)\> = 0, \, \forall n>0.
\end{equation}
An evolution in boost can only lead to changes in physics (e.g. change in vacuum structure, changes in spectrum) in the limit of infinite boosts. This $\epsilon$ evolution changes the vacuum continuously and hence cannot be thought of as an evolution in boosts. As we will now see, this evolution in parameter-space is very naturally explained by accelerating string worldsheets. The identification \refb{c1} stays valid near $\epsilon\to0$, while a map for the whole of the parameter space emerges through acceleration.

\smallskip

\noindent {\em{\underline{Reaching the Horizon I: Evolving in acceleration}}}. We now build the string equivalent of a Rindler observer approaching the Rindler horizon by considering a family of worldsheets with increasing values of acceleration. In the limit of large acceleration, the Bogoliubov coefficients \refb{R2M} become:
\bea{binfty}
\b_{n}^{\infty} = \frac{1}{2}\left(\sqrt{\frac{\pi n}{2a}}+\sqrt{\frac{2a}{\pi n}}  \right) \a_{n}+\frac{1}{2}\left(\sqrt{\frac{\pi n}{2a}}-\sqrt{\frac{2a}{\pi n}}\right) \tilde{\a}_{-n}\cr
\tilde{\b}_{n}^{\infty}=\frac{1}{2}\left(\sqrt{\frac{\pi n}{2a}}-\sqrt{\frac{2a}{\pi n}}\right)\a_{-n}+ \frac{1}{2}\left(\sqrt{\frac{2a}{\pi n}}+\sqrt{\frac{\pi n}{2a}}\right)\tilde{\a}_{n}.\cr
\end{eqnarray}
This limit takes us very near the lightcone. Comparing \refb{binfty} with \refb{c1}, we see that we can make the identification
\be{atoe}
C_n = \b_{n}^{\infty},  \ \tilde{C}_n= \tilde{\b}_{n}^{\infty}, \ \epsilon = \frac{\pi n}{2a}
\end{equation}
The limit of zero tension is thus the limit of infinite acceleration
\be{}
\epsilon\to 0 \Rightarrow a \to \infty
\end{equation}
This is because the equivalence has to hold for all $n$. {\footnote{The above expressions of course don't hold for $n=0$ since the zero modes are equal \refb{cexpansion}.}}

The evolution in parameter space alluded to earlier is thus clearly an evolution in terms of accelerated worldsheets defined for all values of acceleration. This picture of a family of accelerated worldsheets ties in nicely with our earlier description of the UR limit as the limit of infinite boost and the limit of infinite acceleration both land up on the horizon of Rindler spacetime, which in the string analogue is equivalent to the null string.

The flow in acceleration is a flow from the tensile to the tensionless string. Hence {\em increasing acceleration amounts to decreasing string tension}, with $a=0$ being the tensile theory and $a\to \infty$ the tensionless null string. This flow is now described for all values of acceleration giving us a complete interpolating solution. The string oscillator construction through the flow is described in terms of the interpolating $\b$ oscillators, defined in \refb{R2M}. At $a=0$, from \refb{R2M}, we see that $\b$ reduces to tensile string $\a$ oscillators. For intermediate values of $a$, i.e. $0<a<\infty$, we have the $\b$ oscillators. Very near the lightcone, as $a\to \infty$, the $\b$ oscillators take the form of the tensionless $C$ oscillators.
\begin{subequations}
\bea{}
a=0 &:& \quad \{ \b_n, \tilde{\b}_n \} \to \{ \a_n, \tilde{\a}_n \} \\
0<a<\infty &:& \quad \{ \b_n (a) , \tilde{\b}_n(a) \} \\
a\to \infty &:& \quad \{ \b_n, \tilde{\b}_n \} \to \{C_n, \tilde{C}_n\}
\end{eqnarray}
\end{subequations}
We will now use these oscillators to understand the evolution of vacuum structure of the string.

\smallskip

\noindent {\em{\underline{Structure of accelerated vacua}}}. The vacuum conditions on these accelerated worldsheets are:
\bea{bcond}
&& \b^\mu_k|0(a)\rangle = (\alpha^\mu_k +\tanh \theta_k\ \tilde{\alpha}^\mu_{-k})|0(a)\rangle=0,\  k>0; \nonumber \\
&& \tilde{\b}^\mu_k|0(a)\rangle = (\tilde{\alpha}^\mu_k+\tanh\theta_k\ \alpha^\mu_{-k})|0(a)\rangle=0.
\end{eqnarray}
where $\tanh{\theta_k}= -\exp{(-\frac{\pi k}{a})}$. Written in terms of the tensile closed string vacuum $|0\>_\a$, the accelerated vacuum is a squeezed state:
\be{vacrltn}
|0(a)\rangle= \prod_{k=1}^\infty \frac{1}{\cosh\theta_k}\exp\left[-\frac{\tanh{\theta_k}}{k}\ \a_k^\dagger\cdot\tilde{\a}_k^\dagger\right]|0\rangle_\alpha.
\end{equation}
As $a\to\infty$, the map to $\epsilon$ \refb{atoe} emerges. We find $\tanh \theta = \frac{\epsilon-1}{\epsilon+1}$, so $\tanh\theta\to -1$ as $\epsilon\to 0$, and the resulting limiting vacuum state $|0\rangle_c$ is
\bea{zc2za}
|0\rangle_c = \lim_{a\to\infty} |0(a)\rangle = \frac{1}{\mathcal{N}} \prod_{n=1}^\infty \exp\left[\frac{1}{n} \a_n^\dagger\cdot\tilde{\a}_n^\dagger\right]|0\rangle_\alpha.
\end{eqnarray}
Here $\mathcal{N}$ is a normalisation constant. Relations between $C$ and $\a$ are invertible and $|0\rangle_\a$ can be expressed in terms of $|0(a)\>$. This would be the string equivalent of a Rindler observer looking at her Minkowski counterpart. As $a\to\infty$, we find:
\be{za2zc}
|0\rangle_\a = \frac{1}{\mathcal{N'}} \prod_{n=1}^\infty \exp\left[- \frac{1}{n} C_n^\dagger\cdot\tilde{C}_n^\dagger\right]|0\rangle_c.
\end{equation}

\begin{figure}[h]
\centering
  \includegraphics[width=90mm]{Rws.jpeg}
  \caption{ Increasing accelerated worldsheets. }
  \label{fig2}
\end{figure}

\smallskip

\noindent {\em{\underline{Boundary states and Null string complementarity}}}. A tensile closed string field $X(\sigma,\t)$ which maps the worldsheet to spacetime is given by a mode expansion analogous to \refb{tensilemode}. D-branes arise as boundary states on the closed string worldsheet CFT. For a boundary located at $\t=0$ on the worldsheet, the possibilities are \cite{Blumenhagen:2009zz}:
\begin{subequations}\label{bst}
\bea{}
\text{N:} \ \partial_\sigma X(\sigma, \t) |B_1\> = 0 \equiv (\a_n + \tilde{\a}_{-n}) |B_1\> = 0. \\
\text{D:} \ \partial_\t X(\sigma, \t) |B_2\> = 0 \equiv (\a_n - \tilde{\a}_{-n}) |B_2\> = 0.
\end{eqnarray}
\end{subequations}
where N and D stand for Neumann and Dirichlet conditions respectively.
This can be solved explicitly to obtain:
\begin{subequations}
\bea{}
\hspace{-0.5cm} \text{N:} \ |B_1\> &=& {\mathcal{N}_1} \prod_{n=1}^\infty \exp\left[-\frac{\a_n^\dagger\cdot\tilde{\a}_n^\dagger}{n} \right]|0\rangle_\alpha, \\
 \text{D:} \ |B_2\> &=& {\mathcal{N}_2} \prod_{n=1}^\infty \exp\left[\frac{\a_n^\dagger\cdot\tilde{\a}_n^\dagger}{n} \right]|0\rangle_\alpha.
\end{eqnarray}
\end{subequations}
We see that in terms of $\a$ oscillators and the $\a$ vacuum \refb{zc2za}, $|0\rangle_c$ is a {\em Dirichlet boundary state} in all directions, while in terms of the $C$ oscillators and the $C$ vacuum \refb{za2zc}, $|0\rangle_\a$ is a {\em Neumann boundary state} (again in all directions). Thus an open string description emerges from the closed string vacuum as tension goes to zero. This can be further clarified by looking at the oscillators. Using \refb{c1} the conditions for the $C$ vacuum \refb{oe} translates into
\bea{}
\left[\left( \sqrt{\epsilon} + \frac{1}{\sqrt{\epsilon}} \right)\frac{\a^\mu_n}{2} + \left( \sqrt{\epsilon} - \frac{1}{\sqrt{\epsilon}} \right) \frac{\tilde{\a}^\mu_{-n}}{2}\right]  |0\>_c = 0 \nonumber
\end{eqnarray}
In the strict limit $\epsilon\to0$, we end up with
\be{}
(\a^\mu_n - \tilde{\a}^\mu_{-n}) |0\>_c=  0.
\end{equation}
In terms of the usual string vacuum $|0\rangle_\a$, this zero tension groundstate $|0\rangle_c$ is thus a D-instanton, which is a Dirichlet boundary state in all spacetime directions \cite{Polchinski:1994fq}. An analogous calculation yields
\be{}
(C^\mu_n + \tilde{C}^\mu_{-n}) |0\rangle_\a = 0
\end{equation}
Thus, from the point of view of the $C$ observer, the tensile string vacuum develops into a D-25 brane.

We now physically describe this process, the one we will call the {\em null string complementarity}. For an observer in tensile string theory with vacuum $|0\>_\a$ looking at strings with decreasing tension, the completely tensionless string appears as a spacetime point, a D-instanton. There is a complementary point of view of accelerated vacua looking at $|0\>_\a$ vacuum. To the observers in this continuous 1-parameter family of vacua, the usual closed string looks more and more distorted and ultimately the tensionless observer looking at the usual string sees a spacefilling D-25 brane \cite{Bagchi:2019cay}. This is a closed to open string transition. The complementary picture of the formation of the D-instanton fits in rather wonderfully with Rindler worldsheets. This is a ``dual" picture of the formation of an open string from a closed string in the tensionless limit, as seen by observer sitting in the $\a$-vacuum.

We pictorially depict the above process in Fig \ref{fig2}. The ``inertial" closed string worldsheet is the cylinder on the extreme left with acceleration $a=0$. For increasing accelerations $a_i$ $(a_1<a_2<a_3)$, the worldsheet can be given by increasingly distorted hyperboloids. Ultimately, at $a\to\infty$, the worldsheet becomes the lightcone. The boundary states in \refb{bst} are defined at $\t=0$, hence to understand their formation we consider the $\t=0$ cross-sections depicted at the bottom of Fig 2. Increasingly accelerated worldsheets result in circles of lower and lower radius, until at $a\to\infty$, we get a point. This spacetime point is what is the D-instanton described above mathematically. The complementary picture is that when viewed from the $|0\rangle_c$, $|0\rangle_\a$ becomes a longer and longer string, gradually filling up all of spacetime to form a D-25 brane when the tension goes to zero \cite{Bagchi:2019cay}.

\smallskip
\noindent {\em{\underline{Reaching the Rindler Horizon II}}}. Finally, we discuss reaching the Rindler horizon at constant acceleration by evolving in time. We are interesting in Rindler time. So this is a limit $\eta\to \infty$. We will equivalently view this as
\be{}
\eta \to \eta, \ \xi \to \epsilon \xi, \, \epsilon\to 0.
\end{equation}
To understand this limit, we rewrite the 2d conformal generators in Rindler spacetime (we put $a=1$):
\be{}
\L_n, \bar{{\mathcal{L}}}_n = \pm\frac{ i^n }{2} e^{n (\xi - \eta)}(\partial_\eta \mp \partial_\xi).
\end{equation}
In the limit $\epsilon\to0$, we get
\bea{}
 L_n &=& \L_n - \bar{{\mathcal{L}}}_{-n} = i^n e^{-n \eta}(\partial_\eta - n \xi \partial_\xi), \cr
 M_n &=& \epsilon(\L_n + \bar{{\mathcal{L}}}_{-n}) = - i^n e^{-n \eta} \partial_\xi
\end{eqnarray}
This closes to form the classical part of the BMS algebra \refb{bms} (i.e. $c_L=c_M=0$), as expected. This is again thus the null string, which we had expected. A detailed analysis about aspects of Rindler physics on constant accelerated worldsheets would be presented elsewhere \cite{Bagchi:upcoming}.

\begin{figure}[t]
\centering
  \includegraphics[width=80mm]{katakata.pdf}
  \caption{Equal time-slices of a Rindler worldsheet}
  \label{fig3}
\end{figure}

We now present a particularly intriguing picture that arises naturally on constant acceleration worldsheets. Notice that in Rindler spacetime, as depicted in Fig \ref{fig1}, constant Rindler time $(\eta)$ slices are straight lines through the origin with increasing slope on the Rindler R-wedge, depicted by $\eta_i$ ($\eta_1<\eta_2<\eta_3$). At $\eta\to \infty$, this hits the lightcone. On the L-wedge however, time runs backwards and the same slices are obtained by continuing R-wedge lines backwards into the third quadrant. For the string worldsheet at constant acceleration, the analogous picture is Fig \ref{fig3}. Increasing $\eta$ planes intersects the constant hyperboloid at increasing angles. The $\eta$-evolution of the closed string is shown in the boxes below. The circular closed string at the initial $\eta=0$ slice gets deformed as $\eta$ evolves. The tension decreases, the string gets longer and longer as given by the ellipses of increasing eccentricity. Ultimately, when $\eta\to \infty$, the lightcone is hit, the cross-section becomes a straight-line (an ellipse with eccentricity $=1$). The BMS algebra appears on the worldsheet as shown above. The string becomes tensionless and transitions into an open string from a closed one.

\smallskip
\noindent {\em{\underline{Conclusions}}}.
In summary, we have shown that Rindler physics on the worldsheet of the string captures many intriguing aspects of the tensionless string, including the emergence of open strings from closed strings and the remarkable null string complimentarity. Our results would be particularly interesting for studies of strings near black holes where the near-horizon Rindler spacetime would induce a Rindler structure on the worldsheet. It has been speculated that strings grow long and hence become tensionless near black hole horizons \cite{Susskind:1993aa}. We have provided the tools for understanding this phenomenon from the worldsheet in this work. Details on this, as well as connections to thermal physics and entanglement on the worldsheet, would be elaborated on in \cite{Bagchi:upcoming}.

\smallskip
\noindent {\em{\underline{Acknowledgements}}}. This project has been a long time in the making, and we are grateful to a number of places for hospitality over the last couple of years especially Tsinghua University, Beijing and ULB Brussels, where the two ABs spent interesting times discussing and debating points that have culminated in this work.

We are grateful to Tarek Anous, Rudranil Basu, Diptarka Das, Daniel Grumiller, Jelle Hartong, Gautam Mandal, Alejandro Rosabal and especially Joan Simon for helpful discussions and comments. We thank Ritankar Chatterjee and Pulastya Parekh for collaboration on related work.

A Bagchi is partially supported by a Swarnajayanti fellowship from the Department of Science and Technology and the Science and Engineering Research Board (SERB) India, and other SERB grants: EMR/2016/008037, ERC/2017/000873, MTR/2017/000740. A Banerjee is supported by the Korea Ministry of Education, Science and Technology, ICT and Future Planning, Gyeongsangbuk-do and Pohang City. SC is partially supported by the ISIRD grant 9-252/2016/IITRPR/708.


\section{Introduction}
Machine learning models are increasingly replacing traditional modeling systems because of better predictive performance and scalability. This has resulted in an explosion of the number of machine learning models being used in decision-making systems for a broad spectrum of activities such as credit lending, candidate recruitment etc. This high rate of adoption also means that machine learning models have a significant impact on people and the society at large. Consequently, it is essential to ensure that these models are well-regulated.

Building machine learning models in highly regulated domains such as finance, healthcare etc. often comes with its own share of additional challenges. For instance, in the finance industry, it is essential to be transparent in the decision-making process. These challenges can also arise in the form of legal requirements such as the Equal Credit Opportunity Act (ECOA)\cite{ECOA} which makes it unlawful for any creditor to discriminate against any applicant on the basis of race, sex, color etc.

Regulatory bodies all around the world such as European Union through its General Data Protection and Regulation (GDPR) \cite{GDPR} have been keen on developing ways to make machine learning models fair, accountable, transparent and explainable. A host of bias-mitigation approaches \cite{calmon2017optimized}\cite{FATB}\cite{hardt2016equality}\cite{KamishimaAAS12} have been proposed in recent times that ensure that models are not discriminatory against a specific population. While most of the initial work involved ensuring fairness in the training data being used to train these models, these methods often suffered a loss in model performance. More sophisticated bias mitigation strategies\cite{Zafar}\cite{zhang2018mitigating} have been proposed recently and have been demonstrated to be effective in ensuring fairness, while also not compromising too much on performance. However, most of these bias mitigation strategies are designed for specific classes of models such as neural networks, which often violate the regulatory requirement of being transparent or for models such as logistic regression and support vector machines, which are typically inferior in performance and scalability. Thus, on the one hand, we have simple, interpretable and high-performing machine learning algorithms such as XGBoost that are widely preferred in the finance industry and on the other hand, we have sophisticated bias mitigation strategies that are designed to work with a completely different family of machine learning algorithms such as neural networks which may not be suitable for usage in regulated domains.

Our goal is to bridge this gap by formulating equivalent bias mitigation strategies for more practical algorithms. Specifically, in this paper, we introduce a bias-mitigation scheme for XGBoost \cite{chen2016xgboost}. XGBoost has the desirable advantages of being flexible, explainable and scalable while also providing state-of-the-art performance in most supervised learning tasks. However, XGBoost differs significantly from other convex-margin based classifiers such as Support Vector Machines (SVM), logistic regression, neural networks etc. in the way in which the model parameters are updated. This is the primary reason why existing bias-mitigation methods are not compatible with XGBoost.

Our main contributions in this paper are
\begin{itemize}
    \item the formulation of a regularization-based in-process bias mitigation technique that squarely fits into XGBoost's greedy tree building algorithm
    \item empirical comparison with other state-of-the-art bias mitigation strategies on common benchmarks typically used in the bias-mitigation literature
\end{itemize}

The paper is structured as follows. In Section 2, we provide a brief overview of the existing state-of-the-art bias-mitigation methods. In Section 3, we define our bias-mitigation framework, derive gradient and hessian values for model building and show how it fits into XGBoost's existing framework and requires very little in terms of modification. In Section 4, we show our experimental results on common benchmark datasets in the fairness literature and compare our framework's performance and fairness against state-of-the-art bias mitigation strategies.

\section{Preliminaries}
As noted in the introduction, bias-mitigation is a growing area of research, thanks to the increasing interest in machine learning model regulation. Typically, these regulations prohibit discrimination in decision-making against a certain population characterized by \textit{sensitive attributes} such as race, gender, age, marital status etc. While one might be tempted to assume that excluding these sensitive attributes in the model building process (which is a common practice in the industry) would result in fair models, \cite{disc-aware} demonstrated that this is not always the case. Often times, other proxy variables that are correlated with these sensitive attributes carry sufficient information to induce bias in the model. Additionally, the datasets used to train these models and the processes that generate these datasets might themselves be inherently biased. Hence, we need a sophisticated bias-mitigation algorithm to overcome this problem.

\subsection{Notation and Metrics}
In this paper, we consider a supervised learning task and we adopt the following notation. The model is trained using $n$ training samples denoted by $(x_i, s_i, y_i)_{i=1}^{n}$. For the $i^{th}$ training sample, $x_i$ is the feature vector , ${s_i}$ is the binary indicator for the sensitive attribute (such as sex or race), $y_i$ is the binary target label and $\hat{y_i}$ is the score produced by the model. Note that in the case of binary classification, $\hat{y_i}$ needs to be transformed using the sigmoid function before being interpreted as the output probabilities for the two classes. In order to quantify the level of bias exhibited by the model, we need a fairness metric. Multiple fairness metrics have been proposed and there isn't a one-size-fits-all metric. In this paper, we will be focusing on the following metric.

\begin{definition}
The disparate impact for a model is the ratio of the positive prediction rates of the minority and the majority groups.
\begin{align}
    DI = \frac{P(\hat{Y} = 1 | S = 0)}{P(\hat{Y} = 1 | S = 1)}
\end{align}
\end{definition}
Disparate impact (DI) is a well-established measure of fairness and is often associated with the 80\% rule that is frequently cited in legislation. We consider this metric due to its prevalence in the regulatory law. Throughout this paper, we will be using DI as the metric for fairness and accuracy as the metric for measuring model performance.
\subsection{Prior Work}Prior work on bias mitigation can be broadly classified into three categories: pre-processing, in-processing and post-processing methods. Pre-processing methods \cite{calmon2017optimized}\cite{kamiran2012data} typically project the data into a feature space with fair representations. In-processing methods \cite{FATB}\cite{KamishimaAAS12}\cite{zhang2018mitigating} involve changing the training procedure in order to make the model predictions fair. Post-processing \cite{hardt2016equality}\cite{kamiran2012decision} methods typically transform the model outputs to ensure fairness. Of these three categories, in-processing methods offer maximum robustness and flexibility. Existing in-processing methods can further be classified into four categories: optimization in a space constrained by a fairness metric \cite{Zafar}, a regularized objective function on an unconstrained space where the regularizer is typically a function of the model output and the sensitive feature \cite{KamishimaAAS12}, an adversarial learning set-up where an adversary attempts to identify the correlation between a sensitive attribute and the predictor model's output while the predictor model's goal is to maximize performance and simultaneously fooling the adversary \cite{zhang2018mitigating} and finally, designing meta-algorithms \cite{celis2019classification}.

Both pre-processing and post-processing methods have been widely applied for black-box models. However they are often inflexible and result in degradation of model performance. In-processing methods on the other hand provide robust bias mitigation with a relatively lower performance degradation. \cite{Zafar} proposed an in-process method where the search space for the model is constrained by a fairness metric, namely the co-variance between the sensitive features and the signed distance of the instance from the model's decision boundary. However, their method is applicable only to the family of convex-margin classifiers and not to algorithms such as XGBoost. Kamishima et al.\cite{KamishimaAAS12} proposed the prejudice remover which is a regularization based bias-mitigation strategy. The idea is to add a regularizer that captures the mutual information between $Y$ and $S$. Once again, their method cannot be extended to algorithms such as XGBoost. \cite{zhang2018mitigating} proposed an adversarial setup as described earlier. However, their method suffers from poor convergence characteristics and it is often difficult to tune the adversarial system.

Literature on fairness in ensemble models, especially in boosted tree models is rather limited. To the best of our knowledge, only \cite{fish2015fair} and \cite{FATB} consider fairness in a boosting setup. \cite{fish2015fair} were the first to perform a case study of fairness for Adaboost. Their approach involved pre-processing and post-processing methods such as random reshuffling which can incur additional performance degradation. \cite{FATB} proposed fair adversarial gradient tree boosting where the predictor from \cite{zhang2018mitigating} was a decision tree model. However, the issues of convergence from an adversarial setup still remain, rendering the method to be often impractical.

\subsection{Gradient Boosted Decision Trees}
Gradient Boosted Decision Trees (GBDT), introduced by \cite{friedman2001greedy}, is a boosting framework consisting of a collection of weak learners which are shallow decision trees. \cite{chen2016xgboost} proposed eXtreme Gradient Boosting (XGBoost) as a scalable end-to-end tree boosting algorithm. XGBoost has enjoyed widespread adoption by data scientists in the machine learning community. XGBoost has gained particular interest in finance owing to the fact that it is flexible, scalable and explainable.

The GBDT setup for a dataset $D = \{(x_i, y_i)\}$ involves $K$ additive functions put together to make a prediction. Formally, a GBDT model consists of $K$ trees each represented as $f_t(x)$ built at the $t^{th}$ boosting round. The prediction function is then defined as
\begin{align}
    \hat{y_i} =\sum_{t=1}^{K} f_t(x_i) \nonumber
\end{align}

The trees are built in a greedy manner by optimizing the following objective function
\begin{align}
    \mathcal{L}^{t} = \sum_{i=1}^{n} l(y_i, \hat{y}^{(t-1)}_i + f_t(x)) + \Omega(f_t)
\end{align}
where $l$ is an appropriate loss function that depends on the task at hand and $\Omega$ is a regularizer for the tree structure. For classfication tasks, a common choice is the cross-entropy loss between $y_i$ and $\hat{y_i}$, while the squared error is used for regression tasks.

The key contribution of \cite{chen2016xgboost} is the reformulation of this optimization problem as follows. The objective function in the above equation can be approximated using the Taylor expansion as
\begin{align}
    \mathcal{L}^{t} = \sum_{i=1}^{n} l(y_i, \hat{y}^{(t-1)}_i) + g_if_t(x_i) + \frac{1}{2}h_if_t(x_i)^2 + \Omega(f_t)
\end{align}
where $g_i = \nabla_{\hat{y}^{(t-1)}_i} l(y_i, \hat{y}^{(t-1)}_i)$ and $h_i = \nabla^2_{\hat{y}^{(t-1)}_i} l(y_i, \hat{y}^{(t-1)}_i)$. This objective function is then transformed from the space of $f_t$ to the space of node weights $w_j$ of the trees, which results in the following.
\begin{equation}
    \mathcal{L}^{t} = \sum_{j=1}^{T}\left[w_j\left(\Sigma_{i\in I_j}g_i\right) + \frac{1}{2}w_j^2\left(\Sigma_{i\in I_j} h_j\right)\right] + \Omega(f_t)
\end{equation}
where $I_j$ is the set of indices of the samples that fall in the leaf $j$. The best split is computed as the split value that optimizes this objective function.

\section{Proposed Framework}
Our proposed approach involves the use of a fairness regularizer that aims to remove correlation between the sensitive attribute and the target value, thereby ensuring model fairness. The extent to which the regularizer affects the model is controlled by a hyperparameter.

\subsection{Fairness Regularizer}
Using the notation that we introduced earlier, we have a set of training samples $D = \{(x_i, s_i, y_i)\}$. For convenience, let us assume, without loss of generality $s_i = 1$ represents instance $i$ belonging to the majority group and $s_i = 0 $ represents instance $i$ belonging to the minority group. Similarly, $y_i = 1$ represents a \textit{favourable outcome} (such as approval of a credit application). Let $\hat{y_i}$ be the raw leaf score produced by the model for the $i^{th}$ instance and  $\sigma(z) = \frac{1}{1+e^{-z}}$ be the classic sigmoid function. We propose the following regularizer.
\begin{align}
\mathcal{R}^t = \sum_{i=1}^{n} s_i\  log\left(\sigma(\hat{y}^{(t)}_i)\right) + (1-s_i) \ log\left(1 - \sigma(\hat{y}^{(t)}_i)\right)
\end{align}
This regularizer is the negative cross-entropy between the sensitive feature and the model prediction and thus, by minimizing this quantity, we decrease the correlation between them. This in turn leads to a decrease in the bias of the model.
\subsection{Gradient and Hessian for the Regularized Objective}
The regularized objective function for a supervised classification task will now be the sum of the  classical cross-entropy loss between the model predictions and the ground truth labels and the negative cross-entropy between the model predictions and the sensitive feature.
\begin{align}
    \bar{\mathcal{L}^{t}} = &\sum_{i=1}^{n} - y_i\  log\left(\sigma(\hat{y}^{(t)}_i)\right) - (1-y_i) \ log\left(1 - \sigma(\hat{y}^{(t)}_i)\right) + \Omega(f_t) \  \nonumber \\
    &+ \mu \sum_{i=1}^{n} s_i\  log\left(\sigma(\hat{y}^{(t)}_i)\right) + (1-s_i) \ log\left(1 - \sigma(\hat{y}^{(t)}_i)\right) \nonumber
\end{align}
The hyperparameter $\mu$ determines the strength of the regularizer: the higher the regularizer strength the higher the fairness score of the model. This gives us fine-grained control over the level of fairness we desire.

We re-trace the steps of \cite{chen2016xgboost} and reformulate this objective from the space of functions $f_t$ to the space of node weights $w_j$, by computing the gradient $\bar{g}_i$ and hessian $\bar{h}_i$ for the new objective function as follows. We drop the superscript $(t-1)$ for convenience of notation.

\begin{align}
    \bar{g}_i =\  &\nabla_{\hat{y}_i} \left(\sum_{i=1}^{n} - y_i\  log\left(\sigma\left(\hat{y}_i\right)\right) - (1-y_i) \ log\left(1 - \sigma\left(\hat{y}_i\right)\right)\right) \nonumber \\
     &\ + \ \mu\ \nabla_{\hat{y}_i} \left(\sum_{i=1}^{n} s_i\  log\left(\sigma\left(\hat{y}_i\right)\right) + (1-s_i) \ log\left(1 - \sigma\left(\hat{y}_i\right)\right)\right) \nonumber \\
    \bar{g}_i =\  & \sigma\left(\hat{y}_i\right) - y_i + \mu \left(s_i - \sigma\left(\hat{y}_i\right) \right) \nonumber
\end{align}
Similarly, we can derive $\bar{h}_i$ and we obtain
\begin{align}
\bar{h}_i =\  &\nabla^2_{\hat{y}_i} \left(\sum_{i=1}^{n} - y_i\  log\left(\sigma\left(\hat{y}_i\right)\right) - (1-y_i) \ log\left(1 - \sigma\left(\hat{y}_i\right)\right)\right) \nonumber \\
     &\ + \ \mu\ \nabla^2_{\hat{y}_i} \left(\sum_{i=1}^{n} s_i\  log\left(\sigma\left(\hat{y}_i\right)\right) + (1-s_i) \ log\left(1 - \sigma\left(\hat{y}_i\right)\right)\right) \nonumber \\
    \bar{h}_i=\  & (1-\mu)\sigma\left(\hat{y}_i\right) \left(1 - \sigma\left(\hat{y}_i\right)\right) \nonumber
\end{align}

The rest of the tree building process remains the same as XGBoost except that we use the new $\bar{g}_i$ and $\bar{h}_i$ instead. It is worth noting that comparing $\bar{g}_i$ and $\bar{h}_i$ with $g_i$ and $h_i$ from the original XGBoost formulation, we get the following relationships for the gradient and hessian.
\begin{align}
    \bar{g}_i\ =\  &g_i + \mu \left(s_i - \sigma\left(\hat{y}_i\right) \right) \\
    \bar{h}_i\ =\  &h_i \left(1 - \mu \right)
\end{align}

This simple relationship between the original $g_i, h_i$ and our proposed $\bar{g}_i, \bar{h}_i$ is what makes our approach appealing since it can be directly implemented using the custom objective feature of XGBoost. In the next section, we describe the experimental setup that we used and compare our approach to the current state-of-the art methods. Additionally, we also provide insights on how the fairness of the model changes as we increase $\mu$.

\section{Experimental setup}

Throughout this section, we take the following approach. The best hyperparameter settings for the XGBoost model such as \texttt{max-depth}, \texttt{num-rounds}, \texttt{learning-rate} have been identified as the ones that maximize model accuracy, without the fairness regularizer in place ($\mu = 0$). The same setting of hyperparameters is then used and models are built with different values of $\mu$ and the corresponding accuracy and fairness metrics are measured and reported.

The fairness metric we report is the disparate impact (DI) defined in Section 2. The datasets we use are the standard benchmark datasets that are used in the bias-mitigation literature, which are described below.

The first dataset is the UCI Adult Income dataset \cite{Dua:2019} where the goal is to train a model that can predict if an individual makes more than \$50K as income, given a set of features such as age, capital gains and capital losses. The dataset also contains the sensitive attribute \texttt{sex} which takes on two values \texttt{\{Male, Female\}}. The dataset comprises of more males than females, thus making \texttt{Male} the majority population. Fairness here would imply that the model predictions do not discriminate against \texttt{Female}.

The second dataset that we consider is the COMPAS recidivism dataset \cite{COMPAS}. The model is trained to predict if an individual is likely to re-offend in the future. The dataset contains 13 features about 7000 individuals. This was one of the hallmark datasets that was used in the first major debate on the fairness of machine learning models. The sensitive attribute being considered here is \texttt{race}. Once again, fairness here would mean that no particular race is discriminated against.

In addition to the above, we also consider the two other datasets that were analyzed in \cite{FATB} namely the Bank and the Default datasets. The Default dataset \cite{taiwancc} comprises of 23 features about 30,000 Taiwanese credit card users with class
labels which state whether an individual will default on
payments. The sensitive attribute that is being considered for this dataset is \texttt{sex}.

The Bank dataset \cite{bankdata} consists of 16 features of about 45000 clients of a Portuguese banking institution. The goal of the task is to predict if the client has subscribed to a term deposit. The sensitive attribute is \texttt{age} after it has been encoded in a binary format indicating if a customer is between 33 and 60 years old or not.

\begin{table}
  \caption{Benchmark dataset statistics}
  \label{tab:freq}
  \begin{tabular}{crc}
    \toprule
    Name & Number of rows&Sensitive Attribute\\
    \midrule
    Adult & $\sim32000$ & sex\\
    COMPAS & $\sim7000$ & race\\
    Default & $\sim30000$ & sex\\
    Bank & $\sim45000$ & age\\
  \bottomrule
\end{tabular}
\end{table}

\begin{table}
  \caption{Drop in accuracy to achieve DI $\geq 80\%$}
  \label{tab:di_stats}
  \begin{tabular}{ccccc}
    \toprule
    Dataset & FairXGB & Grari\cite{FATB} & Kamishima\cite{KamishimaAAS12} & Zhang\cite{zhang2018mitigating}\\
    \midrule
    Adult & 4.4\% & \textbf{1.9}\% & 3.0\% & 2.3\%\\
    COMPAS & \textbf{1.0}\% & 4.6\% & 3.7\% & 1.0\%\\
    Default & \textbf{0.0}\% & 0.7\% & 1.0\% & \textbf{0.0}\%\\
    Bank & \textbf{0.5}\% & 0.6\% & 0.7\% & 0.6\%\\
  \bottomrule
\end{tabular}
\end{table}

We compare our work to three in-process bias mitigation methods - prejudice remover \cite{KamishimaAAS12}, the fair adversarial gradient tree boosting \cite{FATB} and adversarial debiasing \cite{zhang2018mitigating}. For prejudice remover, we use the implementation provided by \cite{KamishimaAAS12} for training the model. For \cite{FATB}, we re-use the numbers reported by them as the benchmark since the hyperparameters are unknown and hence we cannot reproduce their results. Since different model families would provide different accuracies to begin with, we measure the drop in accuracy from the vanilla model in order to obtain a DI of at least 80\%, rather than the absolute value of the accuracy itself. The 80\% DI target is arbitrary and is often a useful rule of thumb that comes from the 80-20 rule\cite{EEOC}. The comparisons in the drop in accuracy incurred by our approach against \cite{FATB}, \cite{KamishimaAAS12} and \cite{zhang2018mitigating} is shown in Table \ref{tab:di_stats}.

Our method outperforms \cite{FATB}, \cite{zhang2018mitigating} and \cite{KamishimaAAS12} on all but the Adult dataset, where the drop in accuracy is more pronounced. When compared against \cite{KamishimaAAS12}, our method and \cite{FATB} incur a smaller dip in accuracy, thus showing their effectiveness. An interesting observation is that for the Default dataset, the vanilla XGBoost model that we trained was already satisfying the critera of $DI \geq 80\%$ before any bias mitigation was applied. The same was observed in the case of the vanilla neural network model for \cite{zhang2018mitigating}. We believe that the higher loss in accuracy for the Adult dataset could be explained by low DI of the vanilla model. Thus, the best hyperparameters for the vanilla model need not necessarily be optimal for all values of $\mu$. This is supported by the fact that the adversarial methods \cite{FATB} and \cite{zhang2018mitigating} incur a lower loss in accuracy for the Adult dataset because they use a more complex multi-layer perceptron adversary. A more nuanced method for tuning hyperparameters is hence required for non-adversarial methods and we defer this to future work.

We also plot the variation of DI and accuracy with respect to the weight of the fairness regularizer $\mu$, to visualise the effect of increasing $\mu$ on the DI and accuracy metrics. For each $\mu$, we pick the classifier that achieves the highest accuracy and report it along with the corresponding disparate impact. This is in contrast to some of the prior analyses where the "best" models are chosen as those with the best DI. We believe that studying the fairness of the best-accuracy model is more practical for modelers. In Fig. \ref{fig:adult_compas}, we show the plot for the Adult and COMPAS datasets. The Adult dataset requires a weight in the range between $\mu=0.6$ and $\mu=0.7$ to reach the acceptable range of DI which is anything greater than 0.8, whereas the COMPAS dataset reaches the acceptable range at a lower value of $\mu$, between $\mu=0.2$ and $\mu=0.5$. Similarly, the plots for the Bank and Default datasets are shown in Fig.\ref{fig:bank_default}.

\begin{figure}
  \includegraphics[width=\columnwidth]{adult_compas.png}
  \caption{Disparate Impact of the maximum-accuracy classifier for different values of $\mu$ for the Adult Income Dataset and the COMPAS dataset}
  \Description{}
    \label{fig:adult_compas}

\end{figure}

\begin{figure}

  \includegraphics[width=\linewidth]{bank_default.png}
  \caption{Disparate Impact of the maximum-accuracy classifier for different values of $\mu$ for the Bank dataset and the Default dataset}
  \Description{DI vs thresh}
  \label{fig:bank_default}
\end{figure}

\section{Conclusion}
In this paper, we have described an extension of XGBoost that can be used to build fair machine learning models. Our choice of the regularizer for fairness makes it easy to be incorporated into XGBoost with minimal changes, while also providing fine-grained control of the level of fairness that needs to be imposed. Furthermore, we have compared our method with the current state-of-the-art bias mitigation strategies on common benchmark datasets. While we have only considered the cross-entropy loss between $\hat{y}_i$ and $s_i$, our framework is applicable to other continuous and differentiable regularizer functions as well. Hence, our proposal helps bridge the gap between fairness researchers and practitioners in the finance community.

As future work it would be interesting to tackle the other challenge that is typically faced in well-regulated domains - privacy. There have been methods such as differential privacy that have been proposed for secure sharing of sensitive features to modelers to build their models. Adopting such a methodology for XGBoost would be a good addition that would go a long way. Another interesting direction to pursue would be the monitoring of the regularized objective in order to gain insights on the fairness-accuracy tradeoffs. XGBoost's inherent support to track an evaluation metric could be reused for this task.

It would also be useful to pursue the handling of polyvalent sensitive attributes (such as race which can take on many values such as Asian, White, Hispanic, African-American etc.).

\bibliographystyle{ACM-Reference-Format}