Establishing whether or not the famous first order corner at small quark masses exists in the Columbia plot is one of the major open issues in studies of the phase diagram of QCD. We delve into this problem and present results from our study of the chiral limit in 3-flavor QCD using the Highly Improved Staggered Quark (HISQ/tree) action. We investigate four quark masses, which in the continuum correspond to pion masses in the range 80 MeV to 140 MeV. In our simulations, the temporal lattice size, $N_\tau$, is fixed to be equal to 8, and we explore three different aspect ratios $N_\sigma /N_\tau =3,\ 4$ and $5$. In the pion mass range explored by us, we do not find any direct evidence of the existence of a first order phase transition. We find the quark mass and volume dependence of the chiral observables to be well-described by the finite size scaling functions belonging to $3$-$d$, $O(2)$ universality class. We determine the chiral phase transition temperature at this value of finite lattice spacing to be $T_c = 98_{-6}^{+3}~{\rm MeV}$.