![Page 1: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/1.jpg)
The Lorentz Integral Transform (LIT) The Lorentz Integral Transform (LIT) method and its connections with other method and its connections with other
approachesapproaches
First proposed inV. D. Efros, W. Leidemann and G. Orlandini, Phys. Lett. B338, 130 (1994)
Recent Topical Review:V. D. Efros, W. Leidemann, G. Orlandini and N. Barnea “The Lorentz Integral Transform (LIT) method and its applications to perturbation induced reactions” arXiv:0708.2803
![Page 2: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/2.jpg)
The LIT methodThe LIT method
it is an it is an ab initioab initio method for method for nuclear dynamics nuclear dynamics calculationscalculations and in particular forand in particular for continuumcontinuum dynamicsdynamics
it is general enough to be applied to strong as well as electroit is general enough to be applied to strong as well as electroweak weak reactionsreactions of inclusive as well as exclusive nature of inclusive as well as exclusive nature
the applications so far have been to the applications so far have been to electromagnetic and weakelectromagnetic and weak reactions on light nuclei.reactions on light nuclei.
![Page 3: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/3.jpg)
There are many examples in physics where one uses “integral transform approaches”
T = K F
Integral transform approaches
object of interestaccessible object
There are many classes of problems that are difficult to solve in their original representations. An integral transform "maps" an equation from its original "domain" into another domain. Manipulating and solving the equation in the target domain is sometimes much easier than manipulation and solution in the original domain. The solution is then mapped back to the original domain with the inverse of the integral transform.
![Page 4: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/4.jpg)
T = K F
The instrument “integrates” F with the form of its “window” K and gives T
In experimental physics this is common in that one wantsto extract information on the observable from the datawhich is obtained by means of the instrumentation.
![Page 5: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/5.jpg)
Laplace Kernel
∫†xd3xΦ ∫ e – S() d
In Condensed Matter Physics: = Density OperatorS() = Dynamical Structure FunctionΦis obtained with Monte Carlo Methods
In QCD = quark or gluon creation operatorS() = Hadronic Spectral Function”Borel mass”Φ is obtained by OPE QCD sum rules or Lattice
In Nuclear Physics:= Charge or current density operatorS() = R() “Response” Function (to e.w. excitation)Φ is obtained with Monte Carlo Methods
In theoretical physics:
“Euclidean Correlation Function”
![Page 6: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/6.jpg)
One is able to calculate (or measure) Φ but wants S(),
which is the quantity of direct physical meaning.Problem:The “inversion” of Φ may be an “ill posed problem”
Φ = ∫ dω K(ω,σ) S(ω )
![Page 7: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/7.jpg)
Definition of “well-posed” problems in making mathematical models of physical phenomena ( by Hadamard).: 1. A solution exists 2. The solution is unique 3. The solution depends continuously on the inputs in some reasonable topology .
Continuum problems must often be discretized in order to obtain a numerical solution. While in terms of functional analysis such problems are typically continuous, they may suffer from numerical instability when solved with finite precision i.e. a small error in the initial data can result in much larger errors in the answers i.e. they become “ill-posed” problems
They need to be re-formulated for numerical treatment. Typically this involves including additional assumptions, such as smoothness of solution. This process is known as regularization.
The regularization procedure may be unstable
![Page 8: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/8.jpg)
It is well known that the numerical inversion of the Laplace Transform is a terribly ill-posed problem
![Page 9: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/9.jpg)
a “good” Kernel has to satisfy two requirements
1) one must be able to calculate the integral transform
2) one must be able to invert the transform, minimizing instabilities
![Page 10: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/10.jpg)
What is the perfect Kernel?
![Page 11: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/11.jpg)
What is the perfect Kernel?
the delta-function!
![Page 12: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/12.jpg)
What would be the “perfect” Kernel?
the delta-function!
in fact
Φ S ∫ S d
![Page 13: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/13.jpg)
the LIT method is based on the idea to use one of the so-called “representations of the delta-function”:
it turns out that a very good Kernel is the Lorentzian function
Γ
0
Φ ∫ [ S d
![Page 14: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/14.jpg)
The Lorentz Kernel satisfies the two requirements !
N.1. one can calculate the integral transform
N.2 one is able to invert the transform, minimizing instabilities
![Page 15: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/15.jpg)
Illustration of requirement Illustration of requirement N.1:N.1: one can calculate the integral one can calculate the integral transformtransform
![Page 16: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/16.jpg)
Suppose we want an R() defined as (for example for perturbation induced inclusive reactions)
![Page 17: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/17.jpg)
Closure = 1
Using the Lorentz Kernel one has the following theorem: Φ (ω
0,Γ) =
![Page 18: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/18.jpg)
Closure = 1
Using the Lorentz Kernel one has the following theorem: Φ (ω
0,Γ) =
![Page 19: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/19.jpg)
Summarizing: choosing the 2parameter kernel L(
) the theorem based on
closure states that the integral transform Φ(ω
0,Γ) is given by:
ω0
Γ
where
Φ(ω0,Γ) =
![Page 20: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/20.jpg)
The LIT in practice:The LIT in practice:
is found solving for fixed Γ and many ω0
1.
![Page 21: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/21.jpg)
the overlap is calculated
3.
2.
the transform is inverted
![Page 22: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/22.jpg)
main point of the LIT :main point of the LIT :
Schrödinger-like equation with a source
S =
![Page 23: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/23.jpg)
main point of the LIT :main point of the LIT :Schrödinger-like equation with a source
The solution is unique and has bound state asymptotic behavior
Theorem:
![Page 24: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/24.jpg)
main point of the LIT :main point of the LIT :Schrödinger-like equation with a source
one can apply bound state
methods
The solution is unique and has bound state asymptotic behavior
Theorem:
![Page 25: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/25.jpg)
The LIT methodThe LIT method reduces the reduces the continuumcontinuum problem to a problem to a bound bound
state state problemproblem needs needs onlyonly a “good” method for a “good” method for bound statebound state
calculations (FY, HH, NCSM, ...???)calculations (FY, HH, NCSM, ...???) applies both to applies both to inclusiveinclusive reactions reactions
(straightforward!) and to (straightforward!) and to exclusiveexclusive ones ones has beenhas been benchmarked benchmarked in “directly solvable” in “directly solvable”
systems (A=2,3)systems (A=2,3)
![Page 26: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/26.jpg)
A very good method to solve A very good method to solve bound states:bound states:
similar idea as for the No Core Shell Model method similar idea as for the No Core Shell Model method i.e. use of Effective Interaction i.e. use of Effective Interaction applied to H.H. instead of H.O.applied to H.H. instead of H.O.
avoids the avoids the ΩΩH.O.H.O.
parameter dependenceparameter dependence
fast convergencefast convergence can be applied to A>3 can be applied to A>3
N.Barnea, W.Leidemann, G.O. PRC61(2000)054001
the Effective Interaction in Hyperspherical Harmonics method (EIHH)
![Page 27: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/27.jpg)
For A>3 no other benchmark is possible! NO viable solution of the scattering problem beyond the 3-body break up threshold for A=4 and larger.
The LIT approach is at present the only viable one!
Benchmarks:
with R() calculated with traditional differential equation algorithm in A=2 [ V.D.Efros et al.Phys. Lett. B338, 130 (1994) ]
with R() calculated with Faddeev solutions in the continuum in A=3 [J.Golak et al. [Nucl. Phys. A707 (2002) 365]
![Page 28: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/28.jpg)
Practical calculation of Φ
1. Eigenvalue method:
sum of Lorentzians around ευN
can be expanded on localized functions
![Page 29: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/29.jpg)
Practical calculation of Φ
2. Lanczos method
![Page 30: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/30.jpg)
Illustration of requirement Illustration of requirement N.2: N.2: one can invert the integral one can invert the integral transform minimizing instabilitiestransform minimizing instabilities
![Page 31: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/31.jpg)
Inversion of the LIT: the regularization method
![Page 32: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/32.jpg)
Inversion of the LIT: the regularization method
works quite well with bell-shaped kernels!
![Page 33: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/33.jpg)
Phys Lett. B338 (1994) 130
RL (
)
[MeV
-1]
MeV
test on the Deuteron: R() is the longitudinal (e,e') response function
q =2.3 fm-1
![Page 34: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/34.jpg)
Importance of the regularization:
![Page 35: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/35.jpg)
Remember the practical calculation of Φ
Φ is a sum of Lorentzians around ε
υN
the eigenvalue method implies that
![Page 36: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/36.jpg)
Remember the practical calculation of Φ
Φ is a sum of Lorentzians around ε
υN
the eigenvalue method implies that
![Page 37: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/37.jpg)
Remember the practical calculation of Φ
Φ is a sum of Lorentzians around ε
υN
the eigenvalue method implies that
However, here comes the problem of the Continuum! A regularization is needed!
![Page 38: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/38.jpg)
test on deuteron photodisintegration
h.o. basis:fixed MeV
“true”
Nho
=2400
Nho
=150
![Page 39: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/39.jpg)
test on deuteron photodisintegration
h.o. basis:fixed
“true”
MeV
MeV
MeV
![Page 40: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/40.jpg)
Γ=10 MeV + inversion / regularization
Nho
=150 is enough for accuracies at the % level!!
MeV
ΦN
ho=2400
Nho
=150
![Page 41: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/41.jpg)
ConclusionsConclusions
it allows to calculate reactions to the “far” it allows to calculate reactions to the “far” continuum where the many-body scattering continuum where the many-body scattering problem (all channels!) is not solvableproblem (all channels!) is not solvable
only only bound statebound state technique is needed technique is needed
the LIT represents an accurate viable method on the way from
ab initio NUCLEAR STRUCTURE
ab initio NUCLEAR REACTIONS
also for A>3 (see Bacca's talk)
![Page 42: The Lorentz Integral Transform (LIT) method and its connections](https://reader033.vdocument.in/reader033/viewer/2022051404/586a1ea71a28ab9c7d8b7065/html5/thumbnails/42.jpg)
ConclusionsConclusions
since the LIT is calculated numerically a since the LIT is calculated numerically a regularizationregularization procedure is demanded to procedure is demanded to solve the integral equation (inversion of solve the integral equation (inversion of the LIT)the LIT)
there is no discretization of the continuum: the LIT equation is bound-state like
the bell shaped kernel makes the the bell shaped kernel makes the regularizationregularization procedure “inexpensive”, procedure “inexpensive”, and allows to control instabilities. and allows to control instabilities.