I am facing a non-linear, discrete optimization problem, which I can formulate in this abstract manner: I have a certain non-analytic non-linear real-valued function $f:S \to \mathbb{R}$ which takes values on a certain set $S \subseteq \mathbb{Z}^9$. The goal is to find $\theta^* \in \mathbb{Z}^9$ corresponding to the minimum of the function. The total number of possible configuration of the system is high but finite ($\sim10^7$).
Recent quantum annealing techniques allow to take advantage of quantum effects to solve computationally hard optimization problems. Quantum annealers are designed to minimize Quadratic Unconstrained Binary Optimization (QUBO) problems, i.e. the cost function is of the form
$$F(x)=\sum_{i<j}J_{ij}x_ix_j+\sum_{i}h_{i}x_i$$
with $x\in\{0,1\}^n$ and $J_{ij},h_{i}\in \mathbb{R}$.
How do I remap my problem into a QUBO problem? I can associate each configuration of my discrete problem into a binary number, but what does the matrix $J$ represent in this context? What about my non-analytic function $f$ which has to be minimized? Has anyone attempted this in a similar way and how?
Here https://www.hindawi.com/journals/complexity/2017/8404231/ I found the best formalization of the issue.