Consider the following system of $n$ equations:

\begin{equation}f_j^2 = x_j^2\sum_{i=1}^n A_{ij} f_i \tag{$\star$} \end{equation}

where $A_{ij}\geq 0$ are known constants and where $x_j>0$ for all $j$. I am interested in the concavity properties (log-concavity or quasi-concavity) of $f_j$ as a function of the vector $x$. In general, $f_j$ is not concave.

Under what conditions on the $A_{ij}$'s is $f_j$ a log-concave or quasi-concave function of the vector $x$?

A few observations:

- The system has a trivial solution at $f=0$. Taking the square root on both sides of the equation, we can use Theorem 3 in this note to find that there is a unique non-trivial positive solution $f$ for a given $x$.
- The same note shows that we can reach this solution by iterating on the mapping. Therefore, if we start with log-concave (for instance) functions and keep iterating we should converge to a log-concave solution if the mapping preserves log-concavity. Sadly, summation does not preserve log-concavity (or quasi-concavity).
- By iterating on the mapping we get the following nested radicals structure $$f_j=x_j\sqrt{\sum_lA_{lj}x_l\sqrt{\sum_iA_{il}x_i\sqrt{\dots}}}$$
- $f$ is homogeneous of degree 2 in $x$. There are some results about the concavity of homogenous polynomial (see for instance the top of page 127) but again they do not apply directly here. Maybe another version of these results would?
- In simple cases in which we can solve the problem by hand, I find log-concave solutions. For instance, with $A_{ij}=1$ if $i=j$ and 0 otherwise we get $f_j=x_j^2$ which is log-concave. I am hoping that this is a general property of the solution.

Motivation:

- I am interested in the optimization problem $$ \max_{0\leq x_j\leq 1,\,\forall j} \log\left(\sum_{j=1}^n f_j(x)\right)+\log\left(L-\sum_{j=1}^n x_j\right) $$ where $f_j$ satisfies equation $\star$ above and where $L>n$.
- Understanding the concavity of $f_j$ is a step in understanding the concavity of the objective function.
- It is possible to replace the equality in $\star$ by $\leq$ and consider an optimization problem with respect to both $f$ and $x$ with the modified equation $\star$ as a constraint.
- One restriction that would natural in the context of my problem would be to set $A_{ii}=0$ for all $i$.

Please feel free to assume additional assumptions on the problem if needed.