  We study the volatility time series of 1137 most traded stocks in the US
stock markets for the two-year period 2001-02 and analyze their return
intervals $\tau$, which are time intervals between volatilities above a given
threshold $q$. We explore the probability density function of $\tau$,
$P_q(\tau)$, assuming a stretched exponential function, $P_q(\tau) \sim
e^{-\tau^\gamma}$. We find that the exponent $\gamma$ depends on the threshold
in the range between $q=1$ and 6 standard deviations of the volatility. This
finding supports the multiscaling nature of the return interval distribution.
To better understand the multiscaling origin, we study how $\gamma$ depends on
four essential factors, capitalization, risk, number of trades and return. We
show that $\gamma$ depends on the capitalization, risk and return but almost
does not depend on the number of trades. This suggests that $\gamma$ relates to
the portfolio selection but not on the market activity. To further characterize
the multiscaling of individual stocks, we fit the moments of $\tau$, $\mu_m
\equiv <(\tau/<\tau>)^m>^{1/m}$, in the range of $10 < <\tau> \le 100$ by a
power-law, $\mu_m \sim <\tau>^\delta$. The exponent $\delta$ is found also to
depend on the capitalization, risk and return but not on the number of trades,
and its tendency is opposite to that of $\gamma$. Moreover, we show that
$\delta$ decreases with $\gamma$ approximately by a linear relation. The return
intervals demonstrate the temporal structure of volatilities and our findings
suggest that their multiscaling features may be helpful for portfolio
optimization.
