Euclidean to Minkowski Bethe–Salpeter amplitude and observables

We propose a method to reconstruct the Bethe–Salpeter amplitude in Minkowski space given the Euclidean Bethe–Salpeter amplitude – or alternatively the light-front wave function – as input. The method is based on the numerical inversion of the Nakanishi integral representation and computing the corresponding weight function. This inversion procedure is, in general, rather unstable, and we propose several ways to considerably reduce the instabilities. In terms of the Nakanishi weight function, one can easily compute the BS amplitude, the LF wave function and the electromagnetic form factor. The latter ones are very stable in spite of residual instabilities in the weight function. This procedure allows both, to continue the Euclidean BS solution in the Minkowski space and to obtain a BS amplitude from a LF wave function.


Introduction
Among the methods in quantum field theory and quantum mechanical approaches to relativistic few-body systems, the Bethe-Salpeter (BS) equation [1], based on first principles and existing already for more than 60 years, remains rather popular.
The solution of the BS equation in Euclidean space is much more simple than in Minkowski one. Rather often the Euclidean solution is enough, mainly when one is interested in finding the binding energy. However, in many cases ( e.g. for calculating electromagnetic (EM) form factors) one needs the Minkowski space amplitude or, equivalently, the Euclidean one with complex arguments [2,3]. During the recent years, new methods to find the Minkowski solution for the two-body bound state were developed and proved a e-mail: carbonell@ipno.in2p3.fr b e-mail: tobias@ita.br c e-mail: karmanov@sci.lebedev.ru their efficiency, at least for simple kernels like the oneboson exchange (OBE). Some of them [4][5][6][7][8][9] are based on the Nakanishi integral representation [10,11] of the BS amplitude and provide, at first, the Nakanishi weight function, which then allows one to restore the BS amplitude and the light-front (LF) wave function.
Other methods, based on the appropriate treatment of the singularities in the BS equation, provide the Minkowski solution directly, both for bound and scattering states [12].
The straightforward numerical extrapolation of the solution from Euclidean to Minkowski space is very unstable and has not achieved any significant progress of practical interest. However, the above mentioned Nakanishi integral representation provides a more reliable method. It is valid for the Euclidean E as well as for the Minkowski M solutions and both solutions are expressed via one and the same Nakanishi weight function g. This integral representation, given in detail in the next section, can be symbolically written in the form and where K E and K M denote, respectively, the Nakanishi twodimensional integral kernels in Euclidean and Minkowski spaces and g the weight function for the two-body bound state. Solving the integral equation (1) relative to g and substituting the result into Eq. (2), one can in principle determine the Minkowski amplitude M starting with the Euclidean one E . This strategy seems applicable since, in contrast to the direct extrapolation E → M , it uses the analytical properties of the BS amplitude which are implemented in the Nakanishi representation.
A similar integral relation exists expressing the LF wave function ψ in terms of the Nakanishi weight function g [6,7]: One of the most remarkable interests of the Nakanishi representation is that it constitutes a common root for the three different formal objects ( E , M , ψ L F ) allowing one to determine any of them once known any other, and thus providing a link between two theoretical approaches -Bethe-Salpeter and LF dynamics -which are in principle and in practice quite different.
Since the LF wave function can also be found independently either by solving the corresponding equation [13] (see for instance [14]) or by means of quantum field theory (QFT) inspired approaches like the discrete light cone quantization [15] or the basis LF quantization approach [16], one can pose the problem of finding the Nakanishi weight function from Eq. (3). This equation is only one-dimensional. Its solution is more simple, it requires matrices of smaller dimension and therefore it is more stable than the solution of the 2D equation (1).
It is worth noticing that the same method can be applied, far beyond the OBE dynamics, to the complete QFT dynamical content of the lattice calculations. Therein, the full (not restricted to the OBE kernel) Euclidean BS amplitude can be obtained and, via the Nakanishi integral, the corresponding BS amplitude in Minkowski space can be calculated as well as the observables. The Nakanishi representation was used in [17] to calculate the parton distribution amplitude for the pion.
The first results of our research in this field were published in [18]. We have found that the solution of the integral equation (1) relative to g was rather unstable. In order to find numerical solution, this integral equation was discretized using two quite different methods (Gauss quadratures and splines) and it resulted into a linear system. The instability of its the solution dramatically increases with the matrix dimension.
It turned out that Eq. (1) is a Fredholm equation of the first kind, which mathematically is a classical example of an ill-posed problem. On the other hand, this equation has a unique solution. The challenge here is to use an appropriate method allowing us to find this solution. These mathematical methods are developed and well known. Using them and knowing either E or ψ L F we aim to extract g and find from it M and the corresponding observables. This will demonstrate that this procedure is feasible to get the solution.
The appearance of mathematically well defined but numerically ill-posed problems is not a rare exception. It is, for instance, manifested when using the Stieltjes and Lorentz integral transforms [19][20][21] to solve the scattering few-body problems with bound state boundary conditions or when extracting the general parton distributions (GPD) from the experimental data [22].
In this work we will proceed according to the following steps.
In a first step, we will solve numerically equations (1) and (3) in a toy model where E , ψ L F and g are known analytically. The comparison of the numerical solutions for g with the analytical one will tell us the reliability of our method in solving the one-dimensional and two-dimensional Fredholm first kind equations.
In a second step, we will go over a more realistic dynamical case. There, E is obtained by solving the Euclidean BS equation. The LF wave function is given by Eq. (3), where g is calculated by using the methods introduced in [6,7], to solve the Minkowski BS equation, and based on the LF projection and Nakanishi representation. In order to keep trace of its origin we will denote it by g M .
By solving Eqs. (1) and (3), we find g in two independent ways -denoted respectively g E and g L F -to be compared with each other as well as with g M . Our methods should be consistent as far as all these weight functions are reasonably close to each other.
This wave function ψ L F , found by projecting he BS amplitude in the LF, and the ψ L F one obtained from solving the LF equation, are practically indistinguishable from each other. We find also the Minkowski BS amplitude M and using it, we calculate observables, namely, EM form factor and momentum distribution, represented by the LF wave function. We also calculate the form factor independently, expressing it via the LF wave function ψ L F . It turned out that the form factors calculated by these two ways are very close to each other and practically insensitive to instabilities of g's remaining after their suppression by the method we use. That demonstrates that indeed, knowing the Euclidean BS amplitude and using the methods developed in the present paper, one can calculate electroweak observables.
In Sect. 2 we present the formulas for the Nakanishi integral representation, both for the BS amplitude and the LF wave function. A change of variables (mapping) introduced to simplify the integration domain is described in Sect. 3. The numerical method for solving the discretized equation is presented in Sect. 4. The analytically solvable model, on which we will test the numerical solutions, comparing them with the analytical ones, is presented in Sect. 5. In Sect. 6 we study the stability of these solutions. In Sect. 7 we find numerically the Nakanishi weight function using as input the Euclidean BS amplitude and the LF wave function found from the OBE interaction kernel. These results are applied to calculate EM form factors in Sect. 8. Finally, Sect. 9 contains a general discussion and the conclusions.

Nakanishi representation
The BS amplitude in Minkowski space M (k, p), for an Swave two-body bound state with constituents masses m ant total mass M, depends in the rest frame p = (M, 0) on two variables k v and k 0 . We represent the four-vector k as k = (k 0 , k) and denote k v = |k|.
The Nakanishi representation for this amplitude reads [10,11] M (k v , k 0 ) where κ 2 = m 2 − M 2 4 and g is the Nakanishi weight function. The power of the denominator is in fact an arbitrary integer and we have chosen 3 for convenience. This is just the equation which is symbolically written in (2).
In the Euclidean space, after the replacement k 0 = ik 4 , this formula is rewritten as which was represented symbolically in (1). One can also express through g the LF wave function as [6,7]: This equation corresponds to (3). As usually found in the literature, the LF wave function ψ L F is considered as depending on variables k 2 ⊥ , x. They are related to γ, z by γ = k 2 ⊥ , z = 2x − 1 (see e.g. [13]).
The Euclidean BS amplitude E (k v , k 4 ) can easily be found from the corresponding equation with a given kernel (OBE, for instance). As mentioned, it can also be found in lattice calculations [23], which are much difficult numerically but include all the Quantum Field Theory dynamics.
We recall that alternatively to Eq. (3) the LF wave function ψ L F can be found by solving a 2D equation in the LF dynamics [13] or by using other QFT inspired methods like DLCQ [15] or BLFQ [16].
We will therefore assume that one of the functions E or ψ L F is known and can be used as input for solving equations (5) or (6) relative to Nakanishi weight function g. Once obtained in this way, g can be used to calculate any of the remaining quantities in the triplet (ψ, E , M ) and corresponding observables.

Mapping
Though, the time-like momentum variable varies in the range −∞ < k 4 < ∞, we can reduce the problem to the halfinterval 0 < k 4 < ∞ and real arithmetic, by assuming Furthermore, we will make in Eq. (5) the following mapping: Equation (5) takes then the form: where For the normal solutions g(x , −z ) = g(x , z ), which comes from the symmetry of E with respect to k 4 , which separates out the abnormal solutions for the two identical boson case.
After introducing in Eq. (6) for the variables γ, γ the mapping defined in (7), we get where Here z plays the role of a parameter.

Solving Eqs. (8) and (10)
We will look for the solution of Eq. (8) by expanding it on Gegenbauer polynomials in both variables: where is a standard (non-normalized) Gegenbauer polynomial. Whereas, the polynomial G n (x) are orthonormalized: We substitute g(z, x) from (12) in Eq. (8), calculate the integrals numerically and validate the equation (8) transforms into the following linear system: where Equation (14) is a N × N system of linear equations of the type B = A C with the inhomogeneous term given by the Euclidean BS amplitude B ≡ i j and as unknowns the array C ≡ c i j of coefficients of the expansion (12). The solution of this system C = A −1 B, will provide the coefficients c i j and by this the Nakanishi weight function g(x, z) at any point. The solution of Eq. (10) is found similarly. However, since z is a parameter, the solution is found for a fixed z and its decomposition is one-dimensional: We substitute g in Eq. (10), calculate the integral numerically and validate the equation in the N Gaussian quadrature points {x j } in the interval 0 < x < 1. Equation (10) transforms into the inhomogeneous linear system (the parameter z is omitted): where with L(x, x ; z) defined in (11). For a given z, the system (17) of N linear equations is solved, and once determined the coefficients c i , Eq. (16) provides the solution g(x, z) everywhere.
Concerning the number of points N used in the discretization of the integral equations there exists a "plateau of stability", corresponding to an optimal value of N . For a small values of N , the accuracy is not enough but, as we will see below, by increasing N the solution becomes oscillatory and unstable. This is just a manifestation of the above mentioned fact that Eqs. (8) and (10)  However, direct numerical methods do not allow one to find the solution in practice, since they lead to unstable results. To avoid instability, one can, of course, keep N small enough. However, for small N the expansions (12) and (16) give a very crude reproduction of the unknown g. Therefore, to find a solution, we will use a special mathematical method -the Tikhonov regularization method (TRM) [25].
Note that another method (the maximum entropy method) was recently proposed [24] to solve Eq. (5) and successfully applied, at least in the case of monotonic g's.
We will follow here the standard and straightforward way explained above: by discretization of the integral equation we turn it into a matrix equation and solve it by inverting the matrix, however, regularizing the inversion problem. The TRM [25] allows us to find a stable solution for sufficiently large N . Namely, following Ref. [25], we will solve an approximate minimization problem, i.e., we will find C providing the minimum of In the normal form, that corresponds to the replacement of the equation The solution of this equation reads where I is the identity matrix. In our previous work [18] we also used a regularization procedure, but in a more naive form. The equation AC = B was replaced not by A † AC + C = A † B, but by (A+ I )C = B. Then, instead of the solution (19), we got This increases the stability but not so strongly like (19), as we are going to illustrate in the following.

Analytically solvable model
As the first step, we will check the methods and study the appearance of numerical instabilities by solving the simplest one-dimensional equation (6) in a model where the functions ψ L F and g are known analytically. We will compare the numerical solution with the analytical one. In Eq. (6) the values z, m and κ are the parameters which can be chosen arbitrary (keeping the kernel non-singular). We can also put an arbitrary factor at the front of the integral in r.h.-side, which only changes the normalization of ψ L F . Using this freedom, we choose z 2 m 2 + κ 2 (1 − z 2 ) = 1 and replace the factor 1 4 (1 − z 2 ) also by 1. Then Eq. (6) obtains the form As a solution we take Its substitution into (21) provides the l.h.-side: In the mapping variables Eq. (21) is rewritten where the kernel L 1 , the inhomogeneous term ψ L F , and the solution g read We will consider ψ L F , given by Eq. (25) as input, solve numerically, Eq. (24) in the form of the decomposition (16) and compare the solution with the analytical g given in (27).

Studying stability
In Ref. [18] we have found that the inversion of the Nakanishi integral is rather unstable relative to the increase of the number of points N , using the regularization method given by Eq. (20). Therefore, solving Eq. (17), we will first study the onset of instability and its suppression by Tikhonov regularization (19). The solution of the non-regularized equation (17), by using Eq. (20) with = 0 and N = 11, is shown in Fig. 1. It coincides, within the thickness of the lines, with the exact one, given by Eq. (27). However, when increasing N up to N = 14, we can observe (Fig. 2)  Further increase of N results in extremely strong oscillations. The solution for N = 16 and = 0, strongly oscillates and dramatically differs from the exact one, as can be seen in Fig. 3. For N = 32 and = 0, the solution oscillates even  (19) by (20).
Our study shows that using the TRM method there is almost no dependence of the results in a rather wide limits. However, when is taken relatively large (e.g. N = 16, = 10 −4 ), the numerical solution can sensibly differs from the exact one. On the other hand, as mentioned, when using too small one recovers the problem of oscillations, as is seen in Fig. 4  The regularized solution for N = 32, = 10 −10 , found by Eq. (19) and TRM, again coincides with the exact one For higher N the situation is similar. Like in the case N = 32, the solution for N = 64, = 0, found by C = A −1 B strongly oscillates and it strongly differs from the exact one, even stronger than in Fig. 3. The calculated value g(0) ≈ −10 7 instead of g(0) = 1 is again completely wrong. The regularized solution for N = 64, = 10 −10 , found by Eq. (19), also coincides with the exact one within the thickness of lines, like it is in Fig. 1. However, the strong oscillations appear earlier, at ≤ 10 −13 (in contrast to ≤ 10 −16 for N = 32).
We compare now the two ways of regularization given by Eqs. (20) and (19). In Fig. 5 the solution given by Eq. (20) for N = 16, = 10 −10 is shown. It oscillates, while the solution found by using TRM, Eq. (19), with the same N and (not shown), coincides with the analytical one. As mentioned above, the solution for N = 16 and = 10 −18 (Fig. 4) reveals moderate oscillations, which disappear at > 10 −17 . The reason which makes the determinant very small and turns the solution of the linear system in an ill-conditioned problem is the presence of very small eigenvalues of the kernel of the Fredholm integral equation (24). As the dimension of the matrix increases, more small eigenvalues are present, and the eigenstates are oscillating functions, that mixes with the solution obtained within a numerical accuracy. These contributions to the solution are oscillatory, building the pattern seen in Figs. 3 and 5, by increasing dimension of the matrix equation. The amplitude of the eigenvector contribution to the solution increases as the eigenvalue decreases, making the oscillations of the solution divergent. The regularization by cuts the contribution of the small eigenvalues, and to some extend the numerical stability can be found, at the expense of a finite . Both regularization methods from Eq. (19) (Tikhonov method) and (20) improve the numerical solution. The Tikhonov regularization with the reduction to the normal form provides larger eigenvalues as compared to (20), as is evident from the stability analysis that in practice allows much smaller 's before considerable oscillations of the solution appear. See for example Fig. 5, where the method (20) provides huge oscillations. Whereas the TRM allows to be as small as 10 −18 , to be compared with the solution with = 10 −10 . The calculations presented in this section demonstrate that the Nakanishi representation, at least, for the LF wave function, Eq. (6), can be indeed inverted numerically, that is, be solved relative to the Nakanishi weight function g. Though this representation, considered as an equation for g, is an illposed problem, special methods, in particular those based on the Tikhonov regularization (19), allow one to solve it.
In the next section we will show that the Nakanishi representation (5) for the Euclidean BS amplitude can also be inverted. This will be shown not in a toy model, but for the BS solution with a OBE kernel.

OBE interaction
Let us now consider the dynamical case of two spinless particle of unit mass (m = 1), interacting via OBE kernel with exchanged boson mass μ = 0.5, and forming a bound state of total mass M = 1.0.
All the necessary solutions in this model -LF wave function ψ L F , Euclidean BS amplitude E and Nakanishi weight function g -have been computed by solving the corresponding equations. The weight function g was found independently from the LF wave function and BS solutions by solving the equation derived in [6,7,9] with the same OBE kernel. We will use for ψ L F and g the results [9, G. Salmè, private communication] and for E -our own calculations.
We will take profit from the simplicity of the Euclidean solution and extract the Nakanishi weight function g E from E by inverting the Nakanishi representation in Euclidean space (5) via the Tikhonov regularization method and compare it with g found from an equation derived from the Minkowski space BS equation [9, G. Salmè, private communication]. The quality of the solution g E will be checked by computing the LF wave function ψ L F and comparing it with ψ L F found via g provided by [G. Salmè, private communication]. Schematically the above extraction procedure is represented as E → g E → ψ L F . Another possibility to extract g is ψ L F → g L F → E . Namely, starting with the LF wave function, get g L F and calculate with it the Euclidean E and compare this result with the initial E . We will also compare the EM form factors calculated initially via ψ L F and, independently, via M and finally expressed via g. We will see that the observables are almost insensitive to the residual uncertainties in either g E or g L F , which survive after suppressing the instabilities by the Tikhonov regularization (19). This ensures reliable results for the observables, when one uses the Euclidean BS amplitude as input, calculates g by inverting the Nakanishi integral and then with the extracted g goes to the observables. The Nakanishi weight function g E obtained by inverting Eq. (8) with N x = 4, N z = 2 and the TRM with = 10 −10 , is shown in Fig. 6 (dashed line). Solid line represents the direct solution g -denoted below as g F SV -computed in [9, G. Salmè, private communication] by Frederico-Salmè-Viviani by solving a dynamical equation for g and normalized in a different way than E . To compare both solutions for g, we normalize g E so that the two solutions coincide (and equal to −1) at x = 0, z = 0.3.
To check the quality of g E , solution of Eq. (8), we substitute it in r.h.-side of Eq. (8), calculate E and compare it with the input E . The result is shown in Fig. 7. We display by the solid curve the input Euclidean BS amplitude c 2 E as a function of x for a fixed value z = 0.3. The dashed curve corresponds to the BS amplitude E (x, z = 0.3) calculated via g E by Eq. (8). Up to a factor c 2 = 0.88 they are very close to each other.
The corresponding LF wave functions calculated from g E and g F SV by using Eq. (10) are shown in Fig. 8. We multiply ψ L F calculated by Eq. (10) by a normalization factor. Both ψ L F 's well coincide with each other despite the difference between g E and g F SV seen in Fig. 6.
The solution g E for N x = 5, N z = 3, = 10 −10 multiplied by a normalization factor, is shown in Fig. 9. For comparison, we add the solution for N x = 4, N z = 2, = 10 −10 , shown by a dashed line in Fig. 6. In spite of the visible difference between the dashed and solid curves, these two solutions give coinciding LF wave functions and the Euclidean BS amplitudes (after equivalent normalizations), similarly to ones shown in Figs. 7 and 8.

Calculating EM form factor
Once we have the Nakanishi weight function g by Eqs. (4) and (6), we can find both the BS amplitude in Minkowski space and the LF wave function. The EM form factor can be expressed via both of them. In this way, we obtain two expressions in terms of g. We will use both to calculate the form factor and we will compare the results.
The electromagnetic vertex is expressed in terms of the BS amplitude by: where Q 2 = −q 2 , q is the four-momentum transfer, and p = p + q. Substituting M in the form (4) The form factor is expressed via the LF wave function as follows (see e.g. Eq. (6.14) from [13]): where Q 2 ⊥ = Q 2 . Substituting in (30) the LF wave function ψ L F (k ⊥ , x) determined by Eq. (6), one finds [see Eq. (26) from G. Salmè, private communication]: The EM form factor calculated via LF wave function by Eq.
(31), for N x = 4, N z = 2 (dashed curve) and for N x = 5, N z = 3 (solid curve), i.e., via the solutions for g E given in Fig. 8, is shown in Fig. 10. We see that the apparently distinct solutions shown in Fig. 6 does not result in a considerable difference of the corresponding form factors. Notice that the differences represent only a 5% deviation at Q 2 ∼ 10 m 2 The same form factor, calculated via the same Nakanishi weight functions, but in the Minkowski BS framework by Eq. (29) is shown in Fig. 11. The difference between the form factors corresponding to the two solutions for g E shown in Fig. 6 is larger than in Fig. 10, though it is still not so significant.  It is worth to make the two following remarks. ( i) Though the form factors F L F (Q 2 ) and F BSM (Q 2 ) turned out to be very close to each other (as is seen from the present calculations and from G. Salmè, private communication), they should not coincide exactly. Their deviation seen in Fig. 12 is caused not only by numerical uncertainties, but also by physical reasons. This deviation cannot be completely removed by more precise calculations. From point of view of the Fock decomposition of the state vector, the form factor F L F (Q 2 ) is determined only by the two-body contribution in the state vector, whereas F BSM (Q 2 ), determined by the two-body BS amplitude, includes implicitly the effect of higher LF Fock states. ( ii) The contribution of the higher Fock components can be found knowing only the two-body component. Indeed, inverting Eq. (6), we extract g L F from the two-body LF wave function ψ L F . It is the same g as enters in the Minkowski BS amplitude (4) and in the form factor (29), including the higher sector contributions.
As just mentioned, the form factor (29) (calculated with g L F extracted from the two-body LF wave function) includes implicitly not only the two-body component contribution but also the higher Fock components. So, this higher Fock sector contribution is found from the two-body component which was taken as input. Though it seems a little bit paradoxical (the possibility to get information as regards higher Fock states from the two-body one), this is, probably, a manifestation of self-consistency of the Fock decomposition embedded in the Nakanishi integral representation. In the fieldtheoretical framework (in which only the BS amplitude can be defined via the Heisenberg field operators [1]), the number of particles is not conserved and the existence of one Fock component requires the existence of other ones. The set of them, corresponding to different numbers of particles, ensures also the correct transformation properties of the full state vector | p , since the Fock components are transformed by dynamical LF boosts in each other.
The fact that the difference between form factors determined by the two-body LF wave function and the form factor including higher components is small was found also in the Wick-Cutkosky model [27]. The contribution of many-body components with n ≥ 3 reaches 36% in the full normalization integral F(0) only for the coupling constant so huge that the total mass of the bound system tends to zero.
A final remark: the stability of the form factors obtained from g E , though quite unexpected in view of the sizable difference of the Nakanishi weight functions (see e.g. Fig. 8), is in fact appropriate. As we have discussed in the analytical example, the reason for the instability of the solution of the linear system is the contribution from the very small eigenvalues. Therefore the difference between the g E 's comes from the corresponding eigenvectors which induce the observed oscillations. However, when using the Minkowski BS amplitude to compute the form factor, the contributions from the eigenvectors with small eigenvalues are damped in the same way they were enhanced in the inversion. It is very likely that this result follows from the fact that the spectra of the Nakanishi kernels in Minkowski and Euclidean spaces are the same.

Discussion and conclusion
We have demonstrated by explicit calculations, that the Nakanishi representations of the Bethe-Salpeter amplitude in Euclidean space ( E ) and the light-front wave function (ψ L F ) can be numerically inverted. If one of these two quantities is known, one can easily calculate the other one as well as the Bethe-Salpeter amplitude in Minkowski space ( M ) and associated observables, like the electromagnetic form factor.
We have developed an analytically solvable model, in the framework of which we compared the accuracy of the numerical solution with the analytical one. Though the inversion of a Fredholm integral equation of the first kind providing g is an ill-posed problem, it can be solved with satisfactory precision by using appropriate methods.
Our best results are obtained with the Tikhonov regularization procedure. This method is rather efficient and it allowed us to successfully overcome the instabilities of the solution which we found in a previous work [18] when computed g from E . In addition, it turns out that the LF wave function and observables are insensitive to the residual uncertainty of g which remains after stabilizing the solution by the regularization. The uncertainties on the LF wave function are within less than 1%, that is much smaller than for g's.
The Nakanishi weight function g was expanded in terms of Gegenbauer polynomials. A small number of terms in this decomposition is required to control the instabilities. This method is very efficient for describing monotonic behaviors but it is not sufficient to reproduce more involved structure of g like the ones provided by the dynamical model. The latter ones generate considerable uncertainties in the inversion procedure. Increase of the number of terms in the decomposition gives more flexibility, but, at the same time, it enhances the numerical instability of the solution. We believe that it would be useful to study other basis or discretization methods, better reflecting the particular functional form of g.
Our approach can find interesting applications to extract Minkowski amplitudes from an Euclidean theory, like for instance lattice QCD. The Euclidean BS amplitude is being currently computed there with the full dynamical contents of the theory. If one is able to extract from it the corresponding Nakanishi weight function g, one can access to time-like form factors, momentum distributions, GPD's and TMD's which are not accessible in a direct way. This could considerably simplify the study of the Minkowski space structure of hadrons from lattice QCD ab initio calculations.